Search

We present a method based on dynamical nonequilibrium molecular dynamics (D-NEMD) that allows one to produce rigorous ensemble averages for the transient regimes. We illustrate the method by describing the formation of convective cells within a two-dimensional fluid system of soft disks in which a gravity field and a thermal gradient are present. We analyze two different physical settings, with the thermal gradient orthogonal or parallel to the gravity field. In both settings, we follow the formation of the convective flows from the initial time, when the perturbation is turned on, to the steady state. In the first setting (orthogonal fields) we investigate several different cases, varying the initial stationary ensemble and the perturbing field. We find that the final steady-state convective cell is independent of the specific sequence of perturbation fields, which only affects the transient behavior. In all cases, we find that the convective roll is formed through a sequence of damped oscillations of the local fields (density, temperature, and velocity), superimposed to an overall relaxation toward the local steady-state values. Then, we show how D-NEMD can be applied to the Rayleigh–Bénard (RB) setting (parallel fields). In these conditions, the convective flow only establishes above a threshold, without a preferred verse of rotation. We analyze only the response to the ignition of the gravity field in a stationary system under the action of a vertical thermal gradient. Also in this case we characterize the transient response by following the evolution of the density, temperature, and velocity fields until the steady-state RB convective cell is formed. The observed transients are similar to those observed in the case of orthogonal fields. However, the final steady states are quite different. Finally, we briefly discuss the conditions for the general applicability of the D-NEMD method. [ABSTRACT FROM AUTHOR]

The solvation shell of aqueous formaldehyde has been studied by ab initio molecular dynamics. Two different DFT approaches using BLYP and PBE functionals were explored. The results show only a slightly different mobility in the solvation shells and allow characterization of the hydrogen bonded structure with a H2C = O··HOH hydrogen bond lifetime of ca. 3 ps. Formaldehyde hydrolysis was studied by following the reverse process, methanediol decomposition, by Blue Moon constrained MD showing that four water molecules are directly involved in the reaction and assisted by the whole hydration shell. The total energy of the aqueous methanediol to formaldehyde inter-conversion process is calculated with a barrier height of ca. 95 kJ mol-1 while the corresponding free energy barrier is only ΔG‡ = 46 kJ mol-1 at 300 K. [ABSTRACT FROM AUTHOR]