This thesis presents a systematic study of the orbital evolution, gravitational wave radiation, and merger remnant of the black hole–neutron star binary merger in full general relativity for the first time. Numerical-relativity simulations are performed using an adaptive mesh refinement code, SimulAtor for Compact objects in Relativistic Astrophysics (SACRA), which adopts a wide variety of ze…