`ensemble`

(overview)¶

This keyword is used to set up an integration method (an integrator). There are different categories of methods accessible via this keyword, which are described on the following pages:

integrators for path integral molecular dynamics simulations

momentum mirror integrator for simulating compressive shock wave

## Units and suggested parameters¶

The units of temperature and pressure for this keyword are K and GPa, respectively.

The temperature coupling constant `<T_coup>`

means \(\tau_T/\Delta t\), where \(\tau_T\) is the relaxation time of the thermostat and \(\Delta t\) is the time step for integration.
We require \(\tau_T/\Delta t \geq 1\) and a good choice is \(\tau_T/\Delta t \approx 100\).

When \(\tau_T/\Delta t > 100000\), the Berendsen thermostat in `npt_ber`

will be completely ignored, leaving only the Berendsen barostat. In this case, the NPT ensemble reduces to the NPH ensemble that can be useful for melting point calculations using the two-phase method. We have not yet achieved a similar NPH ensemble using the stochastic cell rescaling method.

The pressure coupling constant `<p_coup>`

means \(\tau_p/\Delta t\), where \(\tau_p\) is the relaxation time of the barostat and \(\Delta t\) is the time step for integration.
We require \(\tau_p/\Delta t \geq 1\) and a good choice is \(\tau_p/\Delta t \approx 1000\).

The elastic constants are in units of GPa.

## Caveats¶

One should use one and only one instance of this keyword for each run keyword.