Non-Cartesian coordinates for instantaneous normal mode theory of atomic liquids

Thermodynamic properties of a molecular dipolar liquid using the two-phase thermodynamic approach

A modified version of the original two-phase thermodynamic approach has been extended to evaluate the thermodynamic properties of molecular systems. Its basic assumption states that the density of states can be decomposed into solid-like and gas-like components. The solid part has been approximated by that of a set of harmonic oscillators, whereas a subset composed of rough hard spheres has been considered for the gas part. In this new approach, molecules have been modelled as rotating hard spheres that experience elastic collisions. The technique has been tested on a system made up of dipolar diatomic molecules, and it leads to very good results for total entropy, potential energy reference and heat capacity. Translation and rotation solid components of the overall spectra have been compared to the real part of the instantaneous normal mode vibrational densities of states. Similarities between them reinforce the validity of the two-phase thermodynamic approach.

Single trajectory transport coefficients and the energy landscape by molecular dynamics simulations

The Green-Kubo (GK) method is widely used to calculate the transport coefficients of model liquids by Molecular Dynamics (MD) simulation. A reformulation of GK was proposed by Heyes et al. [J. Chem. Phys. 150, 174504 (2019)], which expressed the shear viscosity in terms of a probability distribution function (PDF) of "single trajectory (ST) viscosities," called "viscuits." This approach is extended here to the bulk viscosity, thermal conductivity, and diffusion coefficient. The PDFs of the four STs expressed in terms of their standard deviations (calculated separately for the positive and negative sides) are shown by MD to be statistically the same for the Lennard-Jones fluid. This PDF can be represented well by a sum of exponentials and is independent of system size and state point in the equilibrium fluid regime. The PDF is not well reproduced by a stochastic model. The PDF is statistically the same as that derived from the potential energy, u, and other thermodynamic quantities, indicating that the transport coefficients are determined quantitatively by and follow closely the time evolution of the underlying energy landscape. The PDFs of out-of-equilibrium supercooled high density states are quite different from those of the equilibrium states.

Free energy and entropy of a dipolar liquid by computer simulations

Thermodynamic properties for a system composed of dipolar molecules are computed. Free energy is evaluated by means of the thermodynamic integration technique, and it is also estimated by using a perturbation theory approach, in which every molecule is modeled as a hard sphere within a square well, with an electric dipole at its center. The hard sphere diameter, the range and depth of the well, and the dipole moment have been calculated from properties easily obtained in molecular dynamics simulations. Connection between entropy and dynamical properties is explored in the liquid and supercooled states by using instantaneous normal mode calculations. A model is proposed in order to analyze translation and rotation contributions to entropy separately. Both contributions decrease upon cooling, and a logarithmic correlation between excess entropy associated with translation and the corresponding proportion of imaginary frequency modes is encountered. Rosenfeld scaling law between reduced diffusion and excess entropy is tested, and the origin of its failure at low temperatures is investigated.

Effect of confinement on the mode dynamics of dipole clusters

Dynamical properties of colloidal clusters composed of paramagnetic beads are presented. The clusters were trapped either in a parabolic trough or in a hard-wall confinement. In order to access the dynamics of the ensembles, the instantaneous normal mode (INM) approach is utilized, which uses cluster configurations as an input. The peaks in the mode spectra weaken when the system size is increased and when the coupling strength is lowered. The short-time diffusive properties of the clusters are deduced using the INM technique. It is found that angular diffusion is always larger than radial diffusion regardless of the shape of the external trap. Further, short-time diffusion seems to be almost independent of the coupling strength in the solid regime, but decreases with increasing packing fraction and size of the ensembles. In general, it is found that diffusion is larger for parabolically confined than for hard-wall trapped clusters.