Abstract
The path integral method is formulated by means of the extended Thouless parametrization. The ab initio Hamiltonian and the Zeeman Hamiltonian are used in our model systems. Numerical calculation of physical properties of molecules is performed in terms of the Monte Carlo method. The dependence of the magnetization of a hydrogen molecule on the division for the partition function, the temperature and the external magnetic field is discussed.