A nonadiabatic theory of diamagnetic susceptibility of molecules is presented in which the electrons and nuclei are considered to be a united system of charged particles whose motion is simultaneously perturbed by a magnetic field. It is found that on separating out the translational motion of the molecule as a whole, there is certain freedom in choosing the phase of the wave function. Its optimum choice corresponds to the gauge of the vector potential with which two contributions opposite in sign to the magnetic susceptibility-the first order diamagnetism and the second order paramagnetism-have minimum magnitudes. Expressions for nonadiabatic calculations of the diamagnetic susceptibility of atoms and molecules are derived. The diamagnetic contributions to the energy of the hydrogen, helium, and lithium atoms, the hydrogen molecule, the pi(-)mu(-)pi(+)mu(+) and p(-)K(-)p(+)K(+) mesomolecules, and the positronium molecule e(-)e(-)e(+)e(+) are calculated. The nonadiabatic contribution of the nuclear motion to the diamagnetic susceptibility amounts to 0.01-0.1% for ordinary atoms and molecules, is increased by several hundred times on passing to mesomolecules, and reaches 50% for the positronium molecule. (C) 2002 MAIK "Nauka/Interperiodica".