pwtools.mttk.particle_nhc_masses

pwtools.mttk.particle_nhc_masses(freq, temp, nf=None, nnos=4)[source]

Fictitious masses Q_p for the particle NHC to thermostat the atoms.

Abinit: qmass for ionmov 13.

Parameters:
  • freq (scalar) – Typical phonon frequency [Hz].

  • temp (scalar) – Temperature [K].

  • nf (scalar) – Number of degrees of freedom. Usually ndim*natoms for no constranis with ndim==3 (x,y,z).

  • nnos (scalar) – Length of the NH chain.

Returns:

  • [q1, qi …] (list of length nnos)

  • q1 (mass for the 1st thermostat in the chain)

  • qi (masses for the 2nd, 3rd, … thermostat)