A lowest-order constrained variational (LOCV) method, with modified conditions of healing on the two-body Jastrow wave function, is investigated in calculations for the ground-state energy levels of many-body spin-polarized atomic deuterium. Results are obtained for the , and phases, corresponding to equal occupations of one, two or three nuclear spin states. Estimates for the optimum healing conditions are obtained by comparison of LOCV results with current Monte Carlo benchmarks. The nature of the phases, i.e., quantum gas or quantum liquid, is discussed, the energy of the phase in our calculations always occurring above the gas-liquid interphase for healing conditions within a physically acceptable range.