An effective variational non-Born-Oppenheimer method is applied to calculate the ground state of the HD molecule in a strong magnetic field. The Hamiltonian used in the calculations is obtained by subtracting the operator representing the kinetic energy of the center-of-mass motion from the total laboratory-frame Hamiltonian. Orbital basis sets are used for the deuteron, the proton, and the electrons. Based on the calculated expectation values, it is determined that, with increasing field strength, the bond length decreases and the alignment of the molecule with the field increases. (C) 2015 Elsevier B.V. All rights reserved.