|NONE||No temperature bath. This is the default.|
|SCALING||Velocity scaling thermostat .|
|BERENDSEN||Berendsen thermostat .|
|HOOVER||Isokinetic Hoover thermostat .|
|NOSE||Nosé-Hoover type thermostat [237,238,239].|
|T=Real||Temperature [K] of BOMD simulation. Default is 300 K.|
|TAU=Real||Berendsen parameter in ps. Default is 0.5 ps.|
|INT=Integer||Interval for velocity scaling. Default is 1.|
|FREQ=Real||Nosé coupling parameter. Default is 1500.0 cm.|
|NHC=Integer||Number of Nosé-Hoover chain thermostats. Default is 3.|
The option SCALING of the BATH keyword activates a simple velocity scaling . The interval for the velocity scaling can be defined with the option INT. By default, the scaling is performed in each step. Figure 11a shows the instantaneous (red) and average (blue) temperature curves for this thermostat in a BOMD run of Cu. The thermostat temperature was set to 300 K.
As this figure shows, velocity scaling is very efficient to adjust both instantaneous and average temperatures to the thermostat value. The freezing of the instantaneous temperature indicates that the phase space distribution is not canonical. Indeed, analysis of the phase space distribution (see 4.7.3) shows a pattern similar to a microcanonical simulation rather than a canonical one . The option BERENDSEN of the BATH keyword selects the BERENDSEN thermostat . This thermostat couples the physical system, namely the molecule, to an external heat bath that is fixed at the desired temperature. The coupling between the molecule and the heat bath can be modified with the TAU option of the BATH keyword. Because of this coupling, the instantaneous temperature can fluctuate around the heat bath temperature. Figure 11b shows the instantaneous (red) and average (blue) temperature curves of the Berendsen thermostat in the Cu simulation. In contrast to the velocity scaling, the fluctuations in the instantaneous temperature remain over the full trajectory length. However, the phase space distribution of the Berendsen thermostat is similar to the one from velocity scaling . Thus, the Berendsen thermostat cannot be recommended for canonical BOMD simulations of molecules. It is, however, well suited for simulated annealing due to its exponential convergence behavior. The option HOOVER of the BATH keyword specifies a thermostat for the so-called constant kinetic temperature dynamics . This method generates correct configurational properties in the canonical ensemble. The instantaneous (red) and average (blue) temperature curves for this thermostat are depicted in Figure 11c, again for the Cu simulation. The repetitive instantaneous temperature fluctuations during the full simulation are characteristic of this thermostat. As a result, it causes only a few different points in phase space to be visited during the full simulation .
The option NOSE of the keyword BATH invokes a Nosé-Hoover (chain) thermostat [237,238,239]. The chain length is specified by the NHC option. By default, three chain links (thermostats) are used. With NHC=1 a pure Nosé-Hoover thermostat is specified. The instantaneous (red) and average (blue) temperature curves for the default Nosé-Hoover chain thermostat with three links are depicted in Figure 11d, again for the Cu simulation. This figure shows that the instantaneous temperature exhibits stochastic fluctuations during the full simulation. Despite these fluctuations, the average temperature of 300 K is reached rather quickly and is stable along the trajectory. The coupling between the Nosé-Hoover thermostat and the real system (molecule), as well as between the thermostats in the chain, is controlled by the FREQ option of the BATH keyword. The larger the FREQ value is, the weaker the coupling is. Because the performance of the Nosé-Hoover (chain) thermostat is sensitive to this coupling frequency, test calculations are recommended. This is particularly important for sensitive properties such as temperature-dependent molecular energies and entropies. The phase space distribution for a canonical N BOMD simulation at 300 K employing the Nosé-Hoover (NHC=1) and Nosé-Hoover chain thermostat with the default value of NHC=3 are depicted in Figure 12. The plotted reduced momenta and coordinates were obtained from a SIMULATION PHASESPACE analysis of the trajectory file. As this figure shows, the Nosé-Hoover chain thermostat yields an almost perfect canonical phase space distribution for this system (see reference  for a detailed discussion). For this reason the Nosé-Hoover chain thermostat is the method of choice for the calculation of molecular properties in canonical ensembles, e.g. the calculation of the temperature dependence of molecular heat capacities.