This paper deals with a hybrid finite element method for wave scattering problems in infinite domains. Scattering of waves involving complex geometries, in conjunction with infinite domains is modeled by introducing a mathematical boundary within which a finite element representation is employed. On the mathematical boundary, the finite element representation is matched with a known analytical solution in the infinite domain in terms of fields and their derivatives. The derivative continuity is implemented by using a slope constraint. Drilling degrees of freedom at each node of the finite element model are introduced to make the numerical model more sensitive to the transverse component of the elastodynamic field. To verify the effects of drilling degrees freedom and slope constraints individually, reflection of normally incident P and SV waves on a traction free half space is considered. For P-wave incidence, the results indicate that the use of a slope constraint is more effective because it suppresses artificial reflection at the mathematical boundary. For the SV-wave case, the use of drilling degrees of freedom is effective in reducing numerical error at the irregular frequencies.