The role of the dipole moment orientations in the crystallization tendency of the van der Waals liquids - molecular dynamics simulations

Role of anisotropy in understanding the molecular grounds for density scaling in dynamics of glass-forming liquids

Molecular Dynamics (MD) simulations of glass-forming liquids play a pivotal role in uncovering the molecular nature of the liquid vitrification process. In particular, much focus was given to elucidating the interplay between the character of intermolecular potential and molecular dynamics behaviour. This has been tried to achieve by simulating the spherical particles interacting via isotropic potential. However, when simulation and experimental data are analysed in the same way by using the density scaling approaches, serious inconsistency is revealed between them. Similar scaling exponent values are determined by analysing the relaxation times and pVT data obtained from computer simulations. In contrast, these values differ significantly when the same analysis is carried out in the case of experimental data. As discussed thoroughly herein, the coherence between results of simulation and experiment can be achieved if anisotropy of intermolecular interactions is introduced to MD simulations. In practice, it has been realized in two different ways: (1) by using the anisotropic potential of the Gay-Berne type or (2) by replacing the spherical particles with quasi-real polyatomic anisotropic molecules interacting through isotropic Lenard-Jones potential. In particular, the last strategy has the potential to be used to explore the relationship between molecular architecture and molecular dynamics behaviour. Finally, we hope that the results presented in this review will also encourage others to explore how 'anisotropy' affects remaining aspects related to liquid-glass transition, like heterogeneity, glass transition temperature, glass forming ability, etc.

Peptide Stereochemistry Effects from pKa-Shift to Gold Nanoparticle Templating in a Supramolecular Hydrogel

The divergent supramolecular behavior of a series of tripeptide stereoisomers was elucidated through spectroscopic, microscopic, crystallographic, and computational techniques. Only two epimers were able to effectively self-organize into amphipathic structures, leading to supramolecular hydrogels or crystals, respectively. Despite the similarity between the two peptides' turn conformations, stereoconfiguration led to different abilities to engage in intramolecular hydrogen bonding. Self-assembly further shifted the pKa value of the C-terminal side chain. As a result, across the pH range 4-6, only one epimer predominated sufficiently as a zwitterion to reach the critical molar fraction, allowing gelation. By contrast, the differing pKa values and higher dipole moment of the other epimer favored crystallization. The four stereoisomers were further tested for gold nanoparticle (AuNP) formation, with the supramolecular hydrogel being the key to control and stabilize AuNPs, yielding a nanocomposite that catalyzed the photodegradation of a dye. Importantly, the AuNP formation occurred without the use of reductants other than the peptide, and the redox chemistry was investigated by LC-MS, NMR, and infrared scattering-type near field optical microscopy (IR s-SNOM). This study provides important insights for the rational design of simple peptides as minimalistic and green building blocks for functional nanocomposites.

The origin of the density scaling exponent for polyatomic molecules and the estimation of its value from the liquid structure

In this article, we unravel the problem of interpreting the density scaling exponent for the polyatomic molecules representing the real van der Waals liquids. Our studies show that the density scaling exponent is a weighted average of the exponents of the repulsive terms of all interatomic interactions that occur between molecules, where the potential energy of a given interaction represents its weight. It implies that potential energy is a key quantity required to calculate the density scaling exponent value for real molecules. Finally, we use the well-known method for potential energy estimation and show that the density scaling exponent could be successfully predicted from the liquid structure for fair representatives of the real systems.

Computational Evidence for the Crucial Role of Dipole Cross-Correlations in Polar Glass-Forming Liquids

In this Letter, we analyze the dipole-dipole correlations obtained from the molecular dynamics simulations for strongly and weakly polar model liquids. As a result, we find that the cross-correlations' contribution to the system's total dipole moment correlation function, which is directly measured in the dielectric spectroscopy experiment, is negligible for weakly polar liquids. In contrast, the cross-correlations' term dominates over the self-correlations' term for the examined strongly polar liquid. Consequently, our studies strongly support the interpretation of the dielectric spectra nature of glass-forming liquids recently proposed by Pabst et al.

The role of the diffusion in the predictions of the classical nucleation theory for quasi-real systems differ in dipole moment value

In this paper, we examine the crystallization tendency for two quasi-real systems, which differ exclusively in the dipole moment's value. The main advantage of the studied system is the fact that despite that their structures are entirely identical, they exhibit different physical properties. Hence, the results obtained for one of the proposed model systems cannot be scaled to reproduce the results for another corresponding system, as it can be done for simple model systems, where structural differences are modeled by the different parameters of the intermolecular interactions. Our results show that both examined systems exhibit similar stability behavior below the melting temperature. This finding is contrary to the predictions of the classical nucleation theory, which suggests a significantly higher crystallization tendency for a more polar system. Our studies indicate that the noted discrepancies are caused by the kinetic aspect of the classical nucleation theory, which overestimates the role of diffusion in the nucleation process.

Virial-potential-energy correlation and its relation to density scaling for quasireal model systems

In this paper, we examine the virial- and the potential-energy correlation for quasireal model systems. This correlation constitutes the framework of the theory of the isomorph in the liquid phase diagram commonly examined using simple liquids. Interestingly, our results show that for the systems characterized by structural anisotropy and flexible bonds, the instantaneous values of total virial and total potential energy are entirely uncorrelated. It is due to the presence of the intramolecular interactions because the contributions to the virial and potential energy resulting from the intermolecular interactions still exhibit strong linear dependence. Interestingly, in contrast to the results reported for simple liquids, the slope of the mentioned linear dependence is different than the values of the density scaling exponent. However, our findings show that for quasireal materials, the slope of dependence between the virial and potential energy (resulting from the intermolecular interactions) strongly depends on the interval of intermolecular distances that are taken into account. Consequently, the value of the slope of the discussed relationship, which enables satisfactory density scaling, can be obtained. Interestingly, this conclusion is supported by the results obtained for analogous systems without intermolecular attraction, for which the value the slope of the virial-potential-energy correlation is independent of considered intermolecular distances, directly corresponds to the exponent of the intermolecular repulsion, and finally leads to accurate density scaling.

Coupling between structural relaxation and diffusion in glass-forming liquids under pressure variation

We theoretically investigate structural relaxation and activated diffusion of glass-forming liquids at different pressures using both Elastically Collective Nonlinear Langevin Equation (ECNLE) theory and molecular dynamics (MD) simulations.