Search

This paper series aims to establish a complete correspondence between fine-grained (FG) and coarse-grained (CG) dynamics by way of excess entropy scaling (introduced in Paper I). While Paper II successfully captured translational motions in CG systems using a hard sphere mapping, the absence of rotational motions in single-site CG models introduces differences between FG and CG dynamics. In this third paper, our objective is to faithfully recover atomistic diffusion coefficients from CG dynamics by incorporating rotational dynamics. By extracting FG rotational diffusion, we unravel, for the first time reported to our knowledge, a universality in excess entropy scaling between the rotational and translational diffusion. Once the missing rotational dynamics are integrated into the CG translational dynamics, an effective translation-rotation coupling becomes essential. We propose two different approaches for estimating this coupling parameter: the rough hard sphere theory with acentric factor (temperature-independent) or the rough Lennard-Jones model with CG attractions (temperature-dependent). Altogether, we demonstrate that FG diffusion coefficients can be recovered from CG diffusion coefficients by (1) incorporating "entropy-free" rotational diffusion with translation-rotation coupling and (2) recapturing the missing entropy. Our findings shed light on the fundamental relationship between FG and CG dynamics in molecular fluids. [ABSTRACT FROM AUTHOR]

Thermal energy storage and utilization has been widely concerned due to the intermittency, renewability, and economy of renewable energy. In this paper, the potential energy function of binary Na2CO3–K2CO3 salt was first constructed using the Deep Potential GENerator (DPGEN) enhanced sampling method. Deep potential molecular dynamics simulations were performed to calculate the thermal properties and structural evolution of binary carbonates. The results show that as the temperature increases from 1073 to 1273 K, the viscosity and thermal conductivity decrease from 5.011 mPa s and 0.502 W/(m K) to 2.526 mPa s and 0.481 W/(m K), respectively. The decrease in viscosity is related to the distance and interaction between the molten salt ions. In addition, the diffusion coefficients, energy barriers, ionic radius, angular distribution function, and coordination number of molten salt were calculated and analyzed. The CO32− exhibits a stable planar triangular structure. The ionic radius of Na+ is smaller than that of K+, which makes Na+ suffer less spatial hindrance during motion and has a higher diffusion coefficient. The energy barriers that Na+ needs to overcome to escape the Coulomb force is greater than that of K+ ions, so molten salt containing Na+ may possess greater heat storage potential. We believe that the potential function constructed with DPGEN enhanced sampling strategy can provide more convincing results for predicting the thermal properties of molten salts. This paper aims to provide a technical route to develop the novel complex molten salt phase change material for thermal energy storage. [ABSTRACT FROM AUTHOR]

This work follows a companion article, which will be referred to as Paper I [Campeggio et al., J. Chem. Phys. 158, 244104 (2023)] in which a quantum-stochastic Liouville equation for the description of the quantum–classical dynamics of a molecule in a dissipative bath has been formulated in curvilinear internal coordinates. In such an approach, the coordinates of the system are separated into three subsets: the quantum coordinates, the classical relevant nuclear degrees of freedom, and the classical irrelevant (bath) coordinates. The equation has been derived in natural internal coordinates, which are bond lengths, bond angles, and dihedral angles. The resulting equation needs to be parameterized. In particular, one needs to compute the potential energy surfaces, the friction tensor, and the rate constants for the nonradiative jumps among the quantum states. While standard methods exist for the calculation of energy and dissipative properties, an efficient evaluation of the transition rates needs to be developed. In this paper, an approximated treatment is introduced, which leads to a simple explicit formula with a single adjustable parameter. Such an approximated expression is compared with the exact calculation of transition rates obtained via molecular dynamics simulations. To make such a comparison possible, a simple sandbox system has been used, with two quantum states and a single internal coordinate (together with its conjugate momentum). Results show that the adjustable parameter, which is an effective decoherence time, can be parameterized from the effective relaxation times of the autocorrelation functions of the conjugated momenta of the relevant nuclear coordinates. [ABSTRACT FROM AUTHOR]

In a recent paper [Aschi et al., ChemPhysChem 6, 53 (2005)], we characterized, by means of theoretical-computational procedures, the electronic excitation of water along the typical liquid state isochore (55.32 mol/l) for a large range of temperature. In that paper we were able to accurately reproduce the experimental absorption maximum at room temperature and to provide a detailed description of the temperature dependence of the excitation spectrum along the isochore. In a recent experimental work by Marin et al. [J. Chem. Phys. 125, 104314 (2006)], water electronic excitation energy was carefully analyzed in a broad range of density and temperature, finding a remarkable agreement of the temperature behavior of the experimental data with our theoretical results. Here, by means of the same theoretical-computational procedures (molecular dynamics simulations and the perturbed matrix method), we investigate water electronic absorption exactly in the same density-temperature range used in the experimental work, hence, now considering also the absorption density dependence. Our results point out that, (1) for all the densities and temperatures investigated, our calculated absorption spectra are in very good agreement with the experimental ones and (2) the gradual maxima redshift observed increasing the temperature or decreasing the density has to be ascribed to a real shift of the lowest X→A electronic transition, supporting the conclusions of Marin et al. [ABSTRACT FROM AUTHOR]

Biased and accelerated molecular simulations (BAMS) are widely used tools to observe relevant molecular phenomena occurring on time scales inaccessible to standard molecular dynamics, but evaluation of the physical time scales involved in the processes is not directly possible from them. For this reason, the problem of recovering dynamics from such kinds of simulations is the object of very active research due to the relevant theoretical and practical implications of dynamics on the properties of both natural and synthetic molecular systems. In a recent paper [A. Rapallo et al., J. Comput. Chem. 42, 586–599 (2021)], it has been shown how the coupling of BAMS (which destroys the dynamics but allows to calculate average properties) with Extended Diffusion Theory (EDT) (which requires input appropriate equilibrium averages calculated over the BAMS trajectories) allows to effectively use the Smoluchowski equation to calculate the orientational time correlation function of the head–tail unit vector defined over a peptide in water solution. Orientational relaxation of this vector is the result of the coupling of internal molecular motions with overall molecular rotation, and it was very well described by correlation functions expressed in terms of weighted sums of suitable time-exponentially decaying functions, in agreement with a Brownian diffusive regime. However, situations occur where exponentially decaying functions are no longer appropriate to capture the actual dynamical behavior, which exhibits persistent long time correlations, compatible with the so called subdiffusive regimes. In this paper, a generalization of EDT will be given, exploiting a fractional Smoluchowski equation (FEDT) to capture the non-exponential character observed in the relaxation of intramolecular distances and molecular radius of gyration, whose dynamics depend on internal molecular motions only. The calculation methods, proper to EDT, are adapted to implement the generalization of the theory, and the resulting algorithm confirms FEDT as a tool of practical value in recovering dynamics from BAMS, to be used in general situations, involving both regular and anomalous diffusion regimes. [ABSTRACT FROM AUTHOR]

Protein folding in confined media has attracted wide attention over the past 15 years due to its importance to both in vivo and in vitro applications. It is generally believed that protein stability increases by decreasing the size of the confining medium, if the medium’s walls are repulsive, and that the maximum folding temperature in confinement is in a pore whose size D0 is only slightly larger than the smallest dimension of a protein’s folded state. Until recently, the stability of proteins in pores with a size very close to that of the folded state has not received the attention it deserves. In a previous paper [L. Javidpour and M. Sahimi, J. Chem. Phys. 135, 125101 (2011)], we showed that, contrary to the current theoretical predictions, the maximum folding temperature occurs in larger pores for smaller α-helices. Moreover, in very tight pores, the free energy surface becomes rough, giving rise to a new barrier for protein folding close to the unfolded state. In contrast to unbounded domains, in small nanopores proteins with an α-helical native state that contain the β structures are entropically stabilized implying that folding rates decrease notably and that the free energy surface becomes rougher. In view of the potential significance of such results to interpretation of many sets of experimental data that could not be explained by the current theories, particularly the reported anomalously low rates of folding and the importance of entropic effects on proteins’ misfolded states in highly confined environments, we address the following question in the present paper: To what extent the geometry of a confined medium affects the stability and folding rates of proteins? Using millisecond-long molecular dynamics simulations, we study the problem in three types of confining media, namely, cylindrical and slit pores and spherical cavities. Most importantly, we find that the prediction of the previous theories that the dependence of the maximum folding temperature Tf on the size D of a confined medium occurs in larger media for larger proteins is correct only in spherical geometry, whereas the opposite is true in the two other geometries that we study. Also studied is the effect of the strength of the interaction between the confined media’s walls and the proteins. If the walls are only weakly or moderately attractive, a complex behavior emerges that depends on the size of the confining medium. [ABSTRACT FROM AUTHOR]

The first paper of this series [J. Chem. Phys. 158, 034103 (2023)] demonstrated that excess entropy scaling holds for both fine-grained and corresponding coarse-grained (CG) systems. Despite its universality, a more exact determination of the scaling relationship was not possible due to the semi-empirical nature. In this second paper, an analytical excess entropy scaling relation is derived for bottom-up CG systems. At the single-site CG resolution, effective hard sphere systems are constructed that yield near-identical dynamical properties as the target CG systems by taking advantage of how hard sphere dynamics and excess entropy can be analytically expressed in terms of the liquid packing fraction. Inspired by classical equilibrium perturbation theories and recent advances in constructing hard sphere models for predicting activated dynamics of supercooled liquids, we propose a new approach for understanding the diffusion of molecular liquids in the normal regime using hard sphere reference fluids. The proposed "fluctuation matching" is designed to have the same amplitude of long wavelength density fluctuations (dimensionless compressibility) as the CG system. Utilizing the Enskog theory to derive an expression for hard sphere diffusion coefficients, a bridge between the CG dynamics and excess entropy is then established. The CG diffusion coefficient can be roughly estimated using various equations of the state, and an accurate prediction of accelerated CG dynamics at different temperatures is also possible in advance of running any CG simulation. By introducing another layer of coarsening, these findings provide a more rigorous method to assess excess entropy scaling and understand the accelerated CG dynamics of molecular fluids. [ABSTRACT FROM AUTHOR]

In 2012, the molecular structure of a new, broad class of ice-binding proteins, classified as "domain of unknown function" (DUF) 3494, was described for the first time. These proteins have a common tertiary structure and are characterized by a very wide spectrum of antifreeze activity (from weakly active to hyperactive). The ice-binding surface (IBS) region of these molecules differs significantly in its structure from the IBS of previously known antifreeze proteins (AFPs), showing a complete lack of regularity and high hydrophilicity. The presence of a regular, repeating structural motif in the IBS region of hitherto known AFP molecules, combined with the hydrophobic nature of this surface, promotes the formation of an ice-like ordering of the solvation water layer and, as a result, facilitates the process of transformation of this water layer into ice. It is, therefore, surprising that the newly discovered DUF3494 class of proteins clearly breaks out of this characteristic. In this paper, using molecular dynamics simulations, we analyze the solvation water structure of the IBS region of both DUF3494 family molecules and AFPs. As we show, although the IBS of DUF3494 molecules does not form an ice-like water structure in the solvation layer, this is compensated by the formation of the equivalent of "anchored clathrate water," in the form of a relatively large number of water molecules bound to the surface of the protein molecule and providing potential binding sites for it to the ice surface. [ABSTRACT FROM AUTHOR]

A central paradigm of polymer physics states that chains in melts behave like random walks as intra- and interchain interactions effectively cancel each other out. Likewise, θ-chains, i.e., chains at the transition from a swollen coil to a globular phase, are also thought to behave like ideal chains, as attractive forces are counterbalanced by repulsive entropic contributions. While the simple mapping to an equivalent Kuhn chain works rather well in most scenarios with corrections to scaling, random walks do not accurately capture the topology and knots, particularly for flexible chains. In this paper, we demonstrate with Monte Carlo and molecular dynamics simulations that chains in polymer melts and θ-chains not only agree on a structural level for a range of stiffnesses but also topologically. They exhibit similar knotting probabilities and knot sizes, both of which are not captured by ideal chain representations. This discrepancy comes from the suppression of small knots in real chains, which is strongest for very flexible chains because excluded volume effects are still active locally and become weaker with increasing semiflexibility. Our findings suggest that corrections to ideal behavior are indeed similar for the two scenarios of real chains and that the structure and topology of a chain in a melt can be approximately reproduced by a corresponding θ-chain. [ABSTRACT FROM AUTHOR]

This paper investigates the dynamic behavior of EthylChlorophyllide A (EChlideA) molecules in a methanol solution confined within a 4 nm silica nanoslit, using molecular dynamics simulations over a duration of 1 ms. Three systems, containing 1, 2, and 4 solutes, were studied at 298 K. The results demonstrate that EChlideA molecules predominantly adsorb onto the silica surfaces, driven by specific interactions between chlorin ring's methyl group and the hydroxyl groups of the silica. This adsorption leads to stable binding, particularly in less crowded environments, as indicated by the potential of mean force analysis. Higher molecular concentrations, such as those with four EChlideA molecules, introduce variation in binding strength due to molecular aggregation and complex interactions. The orientation analysis reveals that the chlorin ring tends to align parallel to the surface, requiring rotational adjustments during surface diffusion. In addition, solvent coordination around the Mg ion remains consistent under bulk conditions, although with some variation in higher concentrations. This study also highlights a decrease in linear diffusion and an increase in rotational relaxation times for EChlideA molecules within the confined nanoslit, reflecting the influence of molecular concentration and arrangement on their dynamics. These findings provide valuable insights into the role of surface interactions, molecular orientation, and solvent coordination in confined environments, offering implications for the design of nanoscale systems. [ABSTRACT FROM AUTHOR]

Liquid–liquid phase transitions hold a unique and profound significance within condensed matter physics. These transitions, while conceptually intriguing, often pose formidable computational challenges. However, recent advances in neural network (NN) potentials offer a promising avenue to effectively address these challenges. In this paper, we delve into the structural transitions of liquid CdTe, CdS, and their alloy systems using molecular dynamics simulations, harnessing the power of an NN potential named LaspNN. Our investigations encompass both pressure and temperature effects. Through our simulations, we uncover three primary liquid structures around melting points that emerge as pressure increases: tetrahedral, rock salt, and close-packed structures, which greatly resemble those of solid states. In the high-temperature regime, we observe the formation of Te chains and S dimers, providing a deeper understanding of the liquid's atomic arrangements. When examining CdSxTe1−x alloys, our findings indicate that a small substitution of S by Te atoms for S-rich alloys (x > 0.5) exhibits a structural transition much different from CdS, while a large substitution of Te by S atoms for Te-rich alloys (x < 0.5) barely exhibits a structural transition similar to CdTe. We construct a schematic diagram for liquid alloys that considers both temperature and pressure, providing a comprehensive overview of the alloy system's behavior. The local aggregation of Te atoms demonstrates a linear relationship with alloy composition x, whereas that of S atoms exhibits a nonlinear one, shedding light on the composition-dependent structural changes. [ABSTRACT FROM AUTHOR]

Understanding how external electric fields (EFs) impact the properties of aqueous molecules is crucial for various applications in chemistry, biology, and engineering. In this paper, we present a study utilizing molecular dynamics simulation to explore how direct-current (DC) and alternative-current (AC) EFs affect hydrophobic (n-triacontane) and hydrophilic (PEG-10) oligomer chains. Through a machine learning approach, we extract a 2-dimensional free energy (FE) landscape of these molecules, revealing that electric fields modulate the FE landscape to favor stretched configurations and enhance the alignment of the chain with the electric field. Our observations indicate that DC EFs have a more prominent impact on modulation compared to AC EFs and that EFs have a stronger effect on hydrophobic chains than on hydrophilic oligomers. We analyze the orientation of water dipole moments and hydrogen bonds, finding that EFs align water molecules and induce more directional hydrogen bond networks, forming 1D water structures. This favors the stretched configuration and alignment of the studied oligomers simultaneously, as it minimizes the disruption of 1D structures. This research deepens our understanding of the mechanisms by which electric fields modulate molecular properties and could guide the broader application of EFs to control other aqueous molecules, such as proteins or biomolecules. [ABSTRACT FROM AUTHOR]

Stable crystalline structures of confined water can be different from bulk ice. In Paper I [T. Yagasaki et al., J. Chem. Phys. 151, 064702 (2019)] of this study, it was shown, using molecular dynamics (MD) simulations, that a zeolite-like ice structure forms in nanobrushes consisting of (6,6) carbon nanotubes (CNTs) when the CNTs are located in a triangle arrangement. The melting temperature of the zeolite-like ice structure is much higher than the melting temperature of ice Ih when the distance between the surfaces of CNTs is ∼0.94 nm, which is the best spacing for the bilayer structure of water. In this paper, we perform MD simulations of nanobrushes of CNTs that are different from (6,6) CNTs in radius. Several new porous ice structures form spontaneously in the MD simulations. A stable porous ice forms when the radius of its cavities matches the radius of the CNTs well. All cylindrical porous ice structures found in this study can be decomposed into a small number of structural blocks. We provide a new protocol to classify cylindrical porous ice crystals on the basis of this decomposition. [ABSTRACT FROM AUTHOR]

Molecular dynamics (MD) is the method of choice for understanding the structure, function, and interactions of molecules. However, MD simulations are limited by the strong metastability of many molecules, which traps them in a single conformation basin for an extended amount of time. Enhanced sampling techniques, such as metadynamics and replica exchange, have been developed to overcome this limitation and accelerate the exploration of complex free energy landscapes. In this paper, we propose Vendi Sampling, a replica-based algorithm for increasing the efficiency and efficacy of the exploration of molecular conformation spaces. In Vendi sampling, replicas are simulated in parallel and coupled via a global statistical measure, the Vendi Score, to enhance diversity. Vendi sampling allows for the recovery of unbiased sampling statistics and dramatically improves sampling efficiency. We demonstrate the effectiveness of Vendi sampling in improving molecular dynamics simulations by showing significant improvements in coverage and mixing between metastable states and convergence of free energy estimates for four common benchmarks, including Alanine Dipeptide and Chignolin. [ABSTRACT FROM AUTHOR]

Interpreting high-dimensional data from molecular dynamics simulations is a persistent challenge. In this paper, we show that for a small peptide, deca-alanine, metastable states can be identified through a neural net based on structural information alone. While processing molecular dynamics data, dimensionality reduction is a necessary step that projects high-dimensional data onto a low-dimensional representation that, ideally, captures the conformational changes in the underlying data. Conventional methods make use of the temporal information contained in trajectories generated through integrating the equations of motion, which forgoes more efficient sampling schemes. We demonstrate that EncoderMap, an autoencoder architecture with an additional distance metric, can find a suitable low-dimensional representation to identify long-lived molecular conformations using exclusively structural information. For deca-alanine, which exhibits several helix-forming pathways, we show that this approach allows us to combine simulations with different biasing forces and yields representations comparable in quality to other established methods. Our results contribute to computational strategies for the rapid automatic exploration of the configuration space of peptides and proteins. [ABSTRACT FROM AUTHOR]

Ninety years ago, Wigner derived the leading order expansion term in ℏ2 for the tunneling rate through a symmetric barrier. His derivation included two contributions: one came from the parabolic barrier, but a second term involved the fourth-order derivative of the potential at the barrier top. He left us with a challenge, which is answered in this paper, to derive the same but for an asymmetric barrier. A crucial element of the derivation is obtaining the ℏ2 expansion term for the projection operator, which appears in the flux-side expression for the rate. It is also reassuring that an analytical calculation of semiclassical transition state theory (TST) reproduces the anharmonic corrections to the leading order of ℏ2. The efficacy of the resulting expression is demonstrated for an Eckart barrier, leading to the conclusion that especially when considering heavy atom tunneling, one should use the expansion derived in this paper, rather than the parabolic barrier approximation. The rate expression derived here reveals how the classical TST limit is approached as a function of ℏ and, thus, provides critical insights to understand the validity of popular approximate theories, such as the classical Wigner, centroid molecular dynamics, and ring polymer molecular dynamics methods. [ABSTRACT FROM AUTHOR]

Vibrational spectroscopy is a useful technique for probing chemical environments. The development of models that can reproduce the spectra of nitriles and azides is valuable because these probes are uniquely suited for investigating complex systems. Empirical vibrational spectroscopic maps are commonly employed to obtain the instantaneous vibrational frequencies during molecular dynamics simulations but often fail to adequately describe the behavior of these probes, especially in its transferability to a diverse range of environments. In this paper, we demonstrate several reasons for the difficulty in constructing a general-purpose vibrational map for methyl thiocyanate (MeSCN), a model for cyanylated biological probes. In particular, we found that electrostatics alone are not a sufficient metric to categorize the environments of different solvents, and the dominant features in intermolecular interactions in the energy landscape vary from solvent to solvent. Consequently, common vibrational mapping schemes do not cover all essential interaction terms adequately, especially in the treatment of van der Waals interactions. Quantum vibrational perturbation (QVP) theory, along with a combined quantum mechanical and molecular mechanical potential for solute–solvent interactions, is an alternative and efficient modeling technique, which is compared in this paper, to yield spectroscopic results in good agreement with experimental FTIR. QVP has been used to analyze the computational data, revealing the shortcomings of the vibrational maps for MeSCN in different solvents. The results indicate that insights from QVP analysis can be used to enhance the transferability of vibrational maps in future studies. [ABSTRACT FROM AUTHOR]

It is significant to investigate the calcium carbonate (CaCO3) precipitation mechanism during the carbon capture process; nevertheless, CaCO3 precipitation is not clearly understood yet. Understanding the carbonation mechanism at the atomic level can contribute to the mineralization capture and utilization of carbon dioxide, as well as the development of new cementitious materials with high-performance. There are many factors, such as temperature and CO2 concentration, that can influence the carbonation reaction. In order to achieve better carbonation efficiency, the reaction conditions of carbonation should be fully verified. Therefore, based on molecular dynamics simulations, this paper investigates the atomic-scale mechanism of carbonation. We investigate the effect of carbonation factors, including temperature and concentration, on the kinetics of carbonation (polymerization rate and activation energy), the early nucleation of calcium carbonate, etc. Then, we analyze the local stresses of atoms to reveal the driving force of early stage carbonate nucleation and the reasons for the evolution of polymerization rate and activation energy. Results show that the higher the calcium concentration or temperature, the higher the polymerization rate of calcium carbonate. In addition, the activation energies of the carbonation reaction increase with the decrease in calcium concentrations. [ABSTRACT FROM AUTHOR]

In liquid water, there is a natural tendency to form aggregates that consist of water molecules linked by hydrogen bonds. Such spontaneously formed aggregates are surrounded by a "sea" of disordered water molecules, with both forms remaining in equilibrium. The process of creating water aggregates also takes place in the solvation water of proteins, but in this case, the interactions of water molecules with the protein surface shift the equilibrium of the process. In this paper, we analyze the structural properties of the solvation water in antifreeze proteins (AFPs). The results of molecular dynamics analysis with the use of various parameters related to the structure of solvation water on the protein surface are presented. We found that in the vicinity of the active region responsible for the binding of AFPs to ice, the equilibrium is clearly shifted toward the formation of "ice-like aggregates," and the solvation water has a more ordered ice-like structure. We have demonstrated that a reduction in the tendency to create "ice-like aggregates" results in a significant reduction in the antifreeze activity of the protein. We conclude that shifting the equilibrium in favor of the formation of "ice-like aggregates" in the solvation water in the active region is a prerequisite for the biological functionality of AFPs, at least for AFPs having a well-defined ice binding area. In addition, our results fully confirm the validity of the "anchored clathrate water" concept, formulated by Garnham et al. [Proc. Natl. Acad. Sci. U. S. A. 108, 7363 (2011)]. [ABSTRACT FROM AUTHOR]

In this paper, we study the permeation of polyatomic gas molecules through 2D graphene membranes. Using equilibrium molecular dynamics simulations, we investigate the permeation of pure gas compounds (CH4, CO2, O2, N2, and H2) through nanoporous graphene membranes with varying pore sizes and geometries. Our simulations consider the recrossing mechanism, often neglected in previous studies, which has a significant effect on permeation for intermediate pore size to molecular diameter ratios. We find that the permeation process can be decoupled into two steps: the crossing process of gas molecules through the pore plane and the escaping process from the pore region to a neighboring adsorption site, which prevents recrossing. To account for these mechanisms, we use a permeance model expressed as the product of the permeance for the crossing process and the probability of molecule escape. This phenomenological model is extended to account for small polyatomic gas molecules and to describe permeation regimes ranging from molecular sieving to effusion. The proposed model captures the temperature dependence and provides insights into the key parameters of the gas/membrane interaction controlling the permeance of the system. This work lays the foundation for predicting gas permeance and exploring membrane separation factors in 2D materials such as graphene. [ABSTRACT FROM AUTHOR]

This paper presents the molecular dynamics simulation of 1 mol kg−1 LiPF6 in a binary solvent of ethylene carbonate (EC) and dimethylcarbonate, which is a representative electrolyte solution for lithium-ion batteries. The simulation successfully reproduced the diffusion coefficient, ionic conductivity, and shear viscosity as functions of EC content at 300 K, which had been experimentally determined in our previous study. The Yukawa potential was adopted to model intercharge interactions to reduce computational costs, which consequently allowed us to precisely calculate the conductivity and viscosity by directly integrating time-correlation functions without explicitly modeling the molecular polarization. Breaking down microscopic current correlation functions into components revealed that, whereas the cation–anion attractive interaction dominantly impedes the conduction when the EC content is low, it is the cation–cation and anion–anion repulsive interactions that reduce the conductivity at a high EC content. An analysis of the pressure correlations revealed that all components positively contribute to the viscosity in the binary solvent without the electrolyte. On the other hand, negative terms are observed in five out of six cross correlations in the presence of the electrolyte, implying that these correlations negatively contribute to the shear stress and entropy production, both of which are net positive. [ABSTRACT FROM AUTHOR]

As biomolecules vibrate and rotate in the terahertz band, the biological effects of terahertz electromagnetic fields have drawn considerable attention from the physiological and medical communities. Ion channels are the basis of biological electrical signals, so studying the effect of terahertz electromagnetic fields on ion channels is significant. In this paper, the effect of a terahertz electromagnetic field with three different frequencies, 6, 15, and 25 THz, on the Kv1.2 potassium ion channel was investigated by molecular dynamics simulations. The results show that an electromagnetic field with a 15 THz frequency can significantly enhance the permeability of the Kv1.2 potassium ion channel, which is 1.7 times higher than without an applied electric field. By analyzing the behavior of water molecules, it is found that the electromagnetic field with the 15 THz frequency shortens the duration of frozen and relaxation processes when potassium ions pass through the channel, increases the proportion of the direct knock-on mode, and, thus, enhances the permeability of the Kv1.2 potassium ion channel. [ABSTRACT FROM AUTHOR]

The collective coupling of an ensemble of molecules to a light field is commonly described by the Tavis–Cummings model. This model includes numerous eigenstates that are optically decoupled from the optically bright polariton states. Accessing these dark states requires breaking the symmetry in the corresponding Hamiltonian. In this paper, we investigate the influence of non-unitary processes on the dark state dynamics in the molecular Tavis–Cummings model. The system is modeled with a Lindblad equation that includes pure dephasing, as it would be caused by weak interactions with an environment, and photon decay. Our simulations show that the rate of pure dephasing, as well as the number of two-level systems, has a significant influence on the dark state population. [ABSTRACT FROM AUTHOR]

The MACE architecture represents the state of the art in the field of machine learning force fields for a variety of in-domain, extrapolation, and low-data regime tasks. In this paper, we further evaluate MACE by fitting models for published benchmark datasets. We show that MACE generally outperforms alternatives for a wide range of systems, from amorphous carbon, universal materials modeling, and general small molecule organic chemistry to large molecules and liquid water. We demonstrate the capabilities of the model on tasks ranging from constrained geometry optimization to molecular dynamics simulations and find excellent performance across all tested domains. We show that MACE is very data efficient and can reproduce experimental molecular vibrational spectra when trained on as few as 50 randomly selected reference configurations. We further demonstrate that the strictly local atom-centered model is sufficient for such tasks even in the case of large molecules and weakly interacting molecular assemblies. [ABSTRACT FROM AUTHOR]

Inelastic collisions of OH with an inert liquid perfluoropolyether (PFPE) surface have been studied experimentally. A pulsed molecular beam of OH with a kinetic energy distribution peaking at 35 kJ mol−1 was directed at a continually refreshed PFPE surface. OH molecules were detected state-selectively with spatial and temporal resolution by pulsed, planar laser-induced fluorescence. The scattered speed distributions were confirmed to be strongly superthermal, regardless of the incidence angle (0° or 45°). Angular scattering distributions were measured for the first time; their reliability was confirmed through extensive Monte Carlo simulations of experimental averaging effects, described in Paper II [A. G. Knight et al., J. Chem. Phys. 158, 244705 (2023)]. The distributions depend markedly on the incidence angle and are correlated with scattered OH speed, consistent with predominantly impulsive scattering. For 45° incidence, the angular distributions are distinctly asymmetric to the specular side but peak at sub-specular angles. This, along with the breadth of the distributions, is incompatible with scattering from a surface that is flat on a molecular scale. New molecular dynamics simulations corroborate the roughness of the PFPE surface. A subtle but unexpected systematic dependence of the angular distribution on the OH rotational state was found, which may be dynamical in origin. The OH angular distributions are similar to those for kinematically similar Ne scattering from PFPE and hence not strongly perturbed by OH being a linear rotor. The results here are broadly compatible with prior predictions from independent quasiclassical trajectory simulations of OH scattering from a model-fluorinated self-assembled monolayer surface. [ABSTRACT FROM AUTHOR]

Nonequilibrium molecular dynamics simulations have been performed for a system of 512 four-site models representing n-butane molecules. The results are compared to simulations previously reported for 125 molecules to check a recently reported number dependence in such simulations. A small system size dependence was observed at low shear rates, but the rheological behavior is not affected. Considerably greater certainty is obtained using larger systems. [ABSTRACT FROM AUTHOR]

A large number of force fields have been proposed for describing the behavior of liquid water within classical atomistic simulations, particularly molecular dynamics. In the past two decades, models that incorporate molecular polarizability and even charge transfer have become more prevalent, in attempts to develop more accurate descriptions. These are frequently parameterized to reproduce the measured thermodynamics, phase behavior, and structure of water. On the other hand, the dynamics of water is rarely considered in the construction of these models, despite its importance in their ultimate applications. In this paper, we explore the structure and dynamics of polarizable and charge-transfer water models, with a focus on timescales that directly or indirectly relate to hydrogen bond (H-bond) making and breaking. Moreover, we use the recently developed fluctuation theory for dynamics to determine the temperature dependence of these properties to shed light on the driving forces. This approach provides key insight into the timescale activation energies through a rigorous decomposition into contributions from the different interactions, including polarization and charge transfer. The results show that charge transfer effects have a negligible effect on the activation energies. Furthermore, the same tension between electrostatic and van der Waals interactions that is found in fixed-charge water models also governs the behavior of polarizable models. The models are found to involve significant energy–entropy compensation, pointing to the importance of developing water models that accurately describe the temperature dependence of water structure and dynamics. [ABSTRACT FROM AUTHOR]

The present paper is the second part of a series of papers aimed at assessing the accuracy of alchemical computational approaches based on nonequilibrium techniques for solvation free energy of organic molecules in the context of molecular dynamics simulations. In Paper I [Procacci, J. Chem. Phys. 151, 144113 (2019)], we dealt with bidirectional estimates of solvation free energies using nonequilibrium approaches. Here, we assess accuracy and precision of unidirectional estimates with the focus on the Gaussian and Jarzynski estimators. We present a very simple methodology to increase the statistics in the work distribution, hence boosting the accuracy and precision of the Jarzynski unidirectional estimates at no extra cost, exploiting the observed decorrelation between the random variables represented by the Lennard-Jones solute-solvent recoupling or decoupling work and by the electrostatic work due to the charging/discharging of the solute in the solvent. [ABSTRACT FROM AUTHOR]

In this paper, we present a combined experimental and theoretical study that explored the initial sticking of water on cooled surfaces. Specifically, these ultra-high vacuum gas–surface scattering experiments utilized supersonic molecular beam techniques in conjunction with a cryogenically cooled highly oriented pyrolytic graphite crystal, giving control over incident kinematic conditions. The D2O translational energy spanning 300–750 meV, the relative D2O flux, and the incident angle could all be varied independently. Three different experimental measurements were made. One involved measuring the total amount of D2O scattering as a function of surface temperature to determine the onset of sticking under non-equilibrium gas–surface collision conditions. Another measurement used He specular scattering to assess structural and coverage information for the interface during D2O adsorption. Finally, we used time-of-flight (TOF) measurements of the scattered D2O to determine how energy is exchanged with the graphite surface at surface temperatures above and near the conditions needed for gaseous condensation. For comparison and elaboration of the roles that internal degrees of freedom play in this process, we also did similar TOF measurements using another mass 20 incident particle, atomic neon. Enriching this study are precise molecular dynamics simulations that elaborate on gas–surface energy transfer and the roles of molecular degrees of freedom in gas–surface collisional energy exchange processes. This study furthers our fundamental understanding of energy exchange and the onset of sticking and ultimately gaseous condensation for gas–surface encounters occurring under high-velocity flows. [ABSTRACT FROM AUTHOR]

This paper reports on the effects of shear rate and interface modeling parameters on the hydrodynamic slip length (LS) for water–graphite interfaces calculated using non-equilibrium molecular dynamics. Five distinct non-bonded solid–liquid interaction parameters were considered to assess their impact on LS. The interfacial force field derivations included sophisticated electronic structure calculation-informed and empirically determined parameters. All interface models exhibited a similar and bimodal LS response when varying the applied shear rate. LS in the low shear rate regime (LSR) is in good agreement with previous calculations obtained through equilibrium molecular dynamics. As the shear rate increases, LS sharply increases and asymptotes to a constant value in the high shear regime (HSR). It is noteworthy that LS in both the LSR and HSR can be characterized by the density depletion length, whereas solid–liquid adhesion metrics failed to do so. For all interface models, LHSR calculations were, on average, ∼28% greater than LLSR, and this slip jump was confirmed using the SPC/E and TIP4P/2005 water models. To address the LS transition from the LSR to the HSR, the viscosity of water and the interfacial friction coefficient were investigated. It was observed that in the LSR, the viscosity and friction coefficient decreased at a similar rate, while in the LSR-to-HSR transition, the friction coefficient decreased at a faster rate than the shear viscosity until they reached a new equilibrium, hence explaining the LS-bimodal behavior. This study provides valuable insights into the interplay between interface modeling parameters, shear rate, and rheological properties in understanding hydrodynamic slip behavior. [ABSTRACT FROM AUTHOR]

Several enhanced sampling techniques rely on the definition of collective variables to effectively explore free energy landscapes. The existing variables that describe the progression along a reactive pathway offer an elegant solution but face a number of limitations. In this paper, we address these challenges by introducing a new path-like collective variable called the "deep-locally non-linear-embedding," which is inspired by principles of the locally linear embedding technique and is trained on a reactive trajectory. The variable mimics the ideal reaction coordinate by automatically generating a non-linear combination of features through a differentiable generalized autoencoder that combines a neural network with a continuous k-nearest neighbor selection. Among the key advantages of this method is its capability to automatically choose the metric for searching neighbors and to learn the path from state A to state B without the need to handpick landmarks a priori. We demonstrate the effectiveness of DeepLNE by showing that the progression along the path variable closely approximates the ideal reaction coordinate in toy models, such as the Müller-Brown potential and alanine dipeptide. Then, we use it in the molecular dynamics simulations of an RNA tetraloop, where we highlight its capability to accelerate transitions and estimate the free energy of folding. [ABSTRACT FROM AUTHOR]

Protein–calixarenes binding plays an increasingly central role in many applications, spanning from molecular recognition to drug delivery strategies and protein inhibition. These ligands obey a specific bio-supramolecular chemistry, which can be revealed by computational approaches, such as molecular dynamics simulations. In this paper, we rely on all-atom, explicit-solvent molecular dynamics simulations to capture the electrostatically driven association of a phosphonated calix-[4]-arene with cytochome-C, which critically relies on surface-exposed paired lysines. Beyond two binding sites identified in direct agreement with the x-ray structure, the association has a larger structural impact on the protein dynamics. Then, our simulations allow a direct comparison to analogous calixarenes, namely, sulfonato, similarly reported as "molecular glue." Our work can contribute to a robust in silico predictive tool to assess binding sites for any given protein of interest for crystallization, with the specificity of a macromolecular cage whose endo/exo orientation plays a role in the binding. [ABSTRACT FROM AUTHOR]

In order to prepare highly heat-resistant packaging insulation materials, in this paper, bismaleimide/epoxy resin (BMI/EP55) composites with different contents of BMI were prepared by melt blending BMI into amino tetrafunctional and phenolic epoxy resin (at a ratio of 5:5). The microstructures and thermal and electrical properties of the composites were tested. The electrostatic potential distribution, energy level distribution, and molecular orbitals of BMI were calculated using Gaussian. The results showed that the carbonyl group in BMI is highly electronegative, implying that the carbonyl group has a strong electron trapping ability. The thermal decomposition temperature of the composites gradually increased with the increase of BMI content, and the 20% BMI/EP55 composites had the highest heat-resistance index, along with a glass transition temperature (Tg) of >250 °C. At different test temperatures, with increase in the BMI content, the conductivity of epoxy resin composites showed a tendency to first decrease and then increase, the breakdown field strength showed a tendency to first increase and then decrease, and the dielectric constant was gradually decreased. Two trap centers were present simultaneously in the composites, where the shallow trap energy level is the deepest in 20% BMI/EP composites and the deep trap energy level is the deepest in 10% BMI/EP55 composites. Correspondingly, the 10% BMI/EP55 composite had a slower charge decay rate, while the 20% BMI/EP55 had a faster charge decay rate. In summary, the BMI/EP55 composites with high heat resistance and insulating properties were prepared in this study, which provided ideas for preparing high-temperature packaging insulating materials. [ABSTRACT FROM AUTHOR]

The large discrepancy among the nucleation kinetics extracted from experimental measurements and computer simulations and the prediction of the classical nucleation theory (CNT) has stimulated intense arguments about its origin in the past decades, which is crucially relevant to the validity of the CNT. In this paper, we investigate the atomistic mechanism of the nucleation in liquid Al in contact with amorphous substrates with atomic-level smooth/rough surfaces, using molecular dynamics (MD) simulations. This study reveals that the slightly distorted local fcc/hcp structures in amorphous substrates with smooth surfaces can promote heterogeneous nucleation through a structural templating mechanism, and on the other hand, homogeneous nucleation will occur at a larger undercooling through a fluctuation mechanism if the surface is rough. Thus, some impurities, previously thought to be impotent, could be activated in the homogeneous nucleation experiments. We further find that the initial growth of the nucleus on smooth surfaces of amorphous substrates is one order of magnitude faster than that in homogeneous nucleation. Both these factors could significantly contribute to the discrepancy in the nucleation kinetics. This study is also supported by a recent study of the synthesis of high-entropy alloy nanoparticles assisted with the liquid metal Ga [Cao et al., Nature 619, 73 (2023)]. In this study, we established that the boundary existed between homogeneous and heterogeneous nucleation, i.e., the structural templating is a general mechanism for heterogeneous nucleation, and in its absence, homogeneous nucleation will occur through the fluctuation mechanism. This study provides an in-depth understanding of the nucleation theory and experiments. [ABSTRACT FROM AUTHOR]

Introducing nanorods into a polymeric matrix can enhance the physical and mechanical properties of the resulting material. In this paper, we focus on understanding the dispersion and orientation patterns of nanorods in an unentangled polymer melt, particularly as a function of nanorod concentration, using molecular dynamics simulations. The system is comprised of flexible polymer chains and multi-thread nanorods that are equilibrated in the NPT ensemble. All interactions are purely repulsive except for those between polymers and rods. Results with attractive vs repulsive polymer–rod interactions are compared and contrasted. The concentration of rods has a direct impact on the phase behavior of the system. At lower concentrations, rods phase separate into nematic clusters, whereas at higher concentrations more isotropic and less structured rod configurations are observed. A detailed examination of the conformation of the polymer chains near the rod surface shows extension of the chains along the director of the rods (especially within clusters). The dispersion and orientation of the nanorods are a result of the competition between depletion entropic forces responsible for the formation of rod clusters, the enthalpic effects that improve mixing of rods and polymer, and entropic losses of polymers interpenetrating rod clusters. [ABSTRACT FROM AUTHOR]

Fewest-switches surface hopping (FSSH) has emerged as one of the leading methods for modeling the quantum dynamics of molecular systems. While its original formulation was limited to adiabatic populations, the growing interest in the application of FSSH to coherent phenomena prompts the question of how one should construct a complete density matrix based on FSSH trajectories. A straightforward solution is to define adiabatic coherences based on wavefunction coefficients. In this paper, we demonstrate that inconsistencies introduced in the density matrix through such treatment may lead to a violation of positivity. We furthermore show that a recently proposed coherent generalization of FSSH results in density matrices that satisfy positivity while yielding improved accuracy throughout much (but not all) of parameter space. [ABSTRACT FROM AUTHOR]

Using an approach based on the diffusion analog of the Cattaneo–Vernotte differential model, we find the exact analytical solution to the corresponding time-dependent linear hyperbolic initial boundary value problem, describing irreversible diffusion-controlled reactions under Smoluchowski's boundary condition on a spherical sink. By means of this solution, we extend exact analytical calculations for the time-dependent classical Smoluchowski rate coefficient to the case that includes the so-called inertial effects, occurring in the host media with finite relaxation times. We also present a brief survey of Smoluchowski's theory and its various subsequent refinements, including works devoted to the description of the short-time behavior of Brownian particles. In this paper, we managed to show that a known Rice's formula, commonly recognized earlier as an exact reaction rate coefficient for the case of hyperbolic diffusion, turned out to be only its approximation being a uniform upper bound of the exact value. Here, the obtained formula seems to be of great significance for bridging a known gap between an analytically estimated rate coefficient on the one hand and molecular dynamics simulations together with experimentally observed results for the short times regime on the other hand. A particular emphasis has been placed on the rigorous mathematical treatment and important properties of the relevant initial boundary value problems in parabolic and hyperbolic diffusion theories. [ABSTRACT FROM AUTHOR]

In this paper, we investigate prenucleation and heterogeneous nucleation in the liquid Pb/solid Al system as an example of systems with large lattice misfit using molecular dynamics simulation. Solid Pb and Al have a large positive lattice misfit (f) of 18.2% along the densely packed [110] direction. This study reveals that prenucleation occurs at 600 K (an undercooling of 15 K), and a 2-dimensional (2D) ordered structure forms at the interface with a coincidence site lattice (CSL) between the first Pb and first Al layers. The CSL accommodates the major part of the f, and only a small residual lattice misfit (fr) of 1.9% remains. The formation of the CSL transforms the original substrate into a considerably potent nucleant, where the first Pb layer becomes the new surface layer of the substrate. At an undercooling of about 22 K, nucleation proceeds by merging 2D ordered structure through structural templating: the second Pb layer is epitaxial to the CSL Pb layer, the third Pb layer largely accommodates the fr, and the fourth Pb layer is a nearly perfect crystalline plane. Further analysis indicates that the interface with the CSL has a lower interfacial energy than with a cube-to-cube orientation relationship. For the first time, we established that the CSL was an effective mechanism to accommodate the f for systems with a large positive misfits. Heterogeneous nucleation is governed not by a single mechanism (misfit dislocations in Turnbull's model), but instead by various mechanisms depending on f. This study sheds new light on the atomistic mechanism of heterogeneous nucleation. [ABSTRACT FROM AUTHOR]

This paper presents a state-averaged complete active space self-consistent field (SA-CASSCF) in the atomic multipole optimized energetics for biomolecular application (AMOEBA) polarizable water model, which enables rigorous simulation of non-adiabatic molecular dynamics with nonequilibrium solvation effects. The molecular orbital and configuration interaction coefficients of the solute wavefunction, and the induced dipoles on solvent atoms, are solved by minimizing the state averaged energy variationally. In particular, by formulating AMOEBA water models and the polarizable continuum model (PCM) in a unified way, the algorithms developed for computing SA-CASSCF/PCM energies, analytical gradients, and non-adiabatic couplings in our previous work can be generalized to SA-CASSCF/AMOEBA by properly substituting a specific list of variables. Implementation of this method will be discussed with the emphasis on how the calculations of different terms are partitioned between the quantum chemistry and molecular mechanics codes. We will present and discuss results that demonstrate the accuracy and performance of the implementation. Next, we will discuss results that compare three solvent models that work with SA-CASSCF, i.e., PCM, fixed-charge force fields, and the newly implemented AMOEBA. Finally, the new SA-CASSCF/AMOEBA method has been interfaced with the ab initio multiple spawning method to carry out non-adiabatic molecular dynamics simulations. This method is demonstrated by simulating the photodynamics of the model retinal protonated Schiff base molecule in water. [ABSTRACT FROM AUTHOR]

This is the final paper in a series that introduces geodesic molecular dynamics at constant potential energy. This dynamics is entitled NVU dynamics in analogy to standard energy-conserving Newtonian NVE dynamics. In the first two papers [T. S. Ingebrigtsen, S. Toxvaerd, O. J. Heilmann, T. B. Schro\der, and J. C. Dyre, J. Chem. Phys. 135, 104101 (2011); T. S. Ingebrigtsen, S. Toxvaerd, T. B. Schro\der, and J. C. Dyre, ibid. 135, 104102 (2011)], a numerical algorithm for simulating geodesic motion of atomic systems was developed and tested against standard algorithms. The conclusion was that the NVU algorithm has the same desirable properties as the Verlet algorithm for Newtonian NVE dynamics, i.e., it is time-reversible and symplectic. Additionally, it was concluded that NVU dynamics becomes equivalent to NVE dynamics in the thermodynamic limit. In this paper, the NVU algorithm for atomic systems is extended to be able to simulate the geodesic motion of molecules at constant potential energy. We derive an algorithm for simulating rigid bonds and test this algorithm on three different systems: an asymmetric dumbbell model, Lewis-Wahnström o-terphenyl (OTP) and rigid SPC/E water. The rigid bonds introduce additional constraints beyond that of constant potential energy for atomic systems. The rigid-bond NVU algorithm conserves potential energy, bond lengths, and step length for indefinitely long runs. The quantities probed in simulations give results identical to those of Nosé-Hoover NVT dynamics. Since Nosé-Hoover NVT dynamics is known to give results equivalent to those of NVE dynamics, the latter results show that NVU dynamics becomes equivalent to NVE dynamics in the thermodynamic limit also for molecular systems. [ABSTRACT FROM AUTHOR]

Polymer crystallization is a long-standing interesting problem both in fundamental polymer physics and in polymer manufacturing. Fundamentally, the connectivity of the macromolecules provides a unique feature for the study of nucleation and growth of crystals in contrast to the crystallization of smaller molecules. In addition, understanding the crystallization in polymers is industrially important due to the necessity of its control to achieve mechanically durable plastic materials. Molecular dynamics simulations offer a suitable way of studying this phenomenon due to their capability to probe the small time and length scales that are characteristic of nucleation. In this paper, we use a long alkane chain model to study nucleation and the growth of polyethylene crystals both within bulk and nano-droplets whose diameters range from about 6 to 11 nm. It is found that the droplets approach being more spherical during the nucleation whereas they deviate from this shape during the growth regime. Strikingly, a mean first passage time analysis indicates that the nucleation rate per unit volume decreases as the droplet size is increased. Finally, visual inspection of the simulation snapshots suggests that nuclei majorly emerged from the surface of the droplets. [ABSTRACT FROM AUTHOR]

The remarkable success of x-ray free-electron lasers and their ability to image biological macromolecules while outrunning secondary radiation damage due to photoelectrons, by using femtosecond pulses, raise the question of whether this can be done using pulsed high-energy electron beams. In this paper, we use excited state molecular dynamics simulations, with tabulated potentials, for rare gas solids to investigate the effect of radiation damage due to inelastic scattering (by plasmons, excitons, and heat) on the pair distribution function. We use electron energy loss spectra to characterize the electronic excitations responsible for radiation damage. [ABSTRACT FROM AUTHOR]

The mechanism of vibrational relaxation and dissociation in the O2–O system at elevated temperatures is investigated by means of molecular dynamics. The most recent O3 potential energy surfaces (PESs), obtained from the first principles quantum mechanical calculations [Varga et al., J. Chem. Phys. 147, 154312 (2017)], are used to derive a complete set of state-specific rate coefficients of vibrational energy transfer and dissociation. Unlike most of the previous efforts that utilize only the lowest and supposedly most reactive 11A′ O3 PES [A. Varandas and A. Pais, Mol. Phys. 65, 843 (1988)], this paper demonstrates the necessity to account for a complete ensemble of all excited O3 PESs that correlate with O2(X) and O(3P) when high-temperature kinetics is of interest. At the same time, it is found that the Varandas 11A′ O3 PES adequately describes vibrational energy transfer and dissociating dynamics when compared to the most recent 11A′ O3 PES by Varga et al. [J. Chem Phys. 147, 154312 (2017)]. The differences between this new dataset and previous rate coefficients are quantified by the master equation model. [ABSTRACT FROM AUTHOR]

Hydration and, in particular, the coordination number of a metal ion is of paramount importance as it defines many of its (bio)physicochemical properties. It is not only essential for understanding its behavior in aqueous solutions but also determines the metal ion reference state and its binding energy to (bio)molecules. In this paper, for divalent metal cations Ca 2 + , Cd 2 + , Cu 2 + , Fe 2 + , Hg 2 + , Mg 2 + , Ni 2 + , Pb 2 + , and Zn 2 + , we compare two approaches for predicting hydration numbers: (1) a mixed explicit/continuum DFT-D3//COSMO-RS solvation model and (2) density functional theory based ab initio molecular dynamics. The former approach is employed to calculate the Gibbs free energy change for the sequential hydration reactions, starting from [M( H 2 O)2]2+ aqua complexes to [M( H 2 O) 9]2+, allowing explicit water molecules to bind in the first or second coordination sphere and determining the most stable [M( H 2 O) n]2+ structure. In the latter approach, the hydration number is obtained by integrating the ion–water radial distribution function. With a couple of exceptions, the metal ion hydration numbers predicted by the two approaches are in mutual agreement, as well as in agreement with the experimental data. [ABSTRACT FROM AUTHOR]

A Markov state model is a powerful tool that can be used to track the evolution of populations of configurations in an atomistic representation of a protein. For a coarse-grained linear chain model with discontinuous interactions, the transition rates among states that appear in the Markov model when the monomer dynamics is diffusive can be determined by computing the relative entropy of states and their mean first passage times, quantities that are unchanged by the specification of the energies of the relevant states. In this paper, we verify the folding dynamics described by a diffusive linear chain model of the crambin protein in three distinct solvent systems, each differing in complexity: a hard-sphere solvent, a solvent undergoing multi-particle collision dynamics, and an implicit solvent model. The predicted transition rates among configurations agree quantitatively with those observed in explicit molecular dynamics simulations for all three solvent models. These results suggest that the local monomer–monomer interactions provide sufficient friction for the monomer dynamics to be diffusive on timescales relevant to changes in conformation. Factors such as structural ordering and dynamic hydrodynamic effects appear to have minimal influence on transition rates within the studied solvent densities. [ABSTRACT FROM AUTHOR]

The network structure within polymers significantly influences their mechanical properties, including their strength, toughness, and fatigue resistance. All-atom molecular dynamics (AAMD) simulations offer a method to investigate the energy dissipation mechanism within polymers during deformation and fracture; Such an approach is, however, computationally inefficient when used to analyze polymers with complex network structures, such as the common chemically double-networked hydrogels. Alternatively, coarse-grained molecular dynamics (CGMD) models, which reduce the computational degrees of freedom by concentrating a set of adjacent atoms into a coarse-grained bead, can be employed. In CGMD simulations, a coarse-grained force field (CGFF) is a critical factor affecting the simulation accuracy. In this paper, we proposed a data-based method for predicting the CGFF parameters to improve the simulation efficiency of complex cross-linked network in polymers. Here, we utilized a typical chemically double-networked hydrogel as an example. An artificial neural network was selected, and it was trained with the tensile stress–strain data from the CGMD simulations using different CGFF parameters. The CGMD simulations using the predicted CGFF parameters show good agreement with the AAMD simulations and are almost fifty times faster. The data-inspired CGMD model presented here broadens the applicability of molecular dynamics simulations to cross-linked polymers and has the potential to provide insights that will aid the design of polymers with desirable mechanical properties. [ABSTRACT FROM AUTHOR]

The fewest switches surface hopping (FSSH) method proposed by Tully in 1990 [Tully, J. Chem. Phys. 93, 1061 (1990)]—along with its many later variations—forms the basis for most practical simulations of molecular dynamics with electronic transitions in realistic systems. Despite its popularity, a rigorous formal derivation of the algorithm has yet to be achieved. In this paper, we derive the energy-conserving momentum jumps employed by FSSH from the perspective of quantum trajectory surface hopping (QTSH) [Martens, J. Phys. Chem. A 123, 1110 (2019)]. In the limit of localized nonadiabatic transitions, simple mathematical and physical arguments allow the FSSH algorithm to be derived from first principles. For general processes, the quantum forces characterizing the QTSH method provide accurate results for nonadiabatic dynamics with rigorous energy conservation, at the ensemble level, within the consistency of the underlying stochastic surface hopping without resorting to the artificial momentum rescaling of FSSH. [ABSTRACT FROM AUTHOR]

Polymer crystallization has long been a fascinating problem and is still attracting many researchers. Most of the previous simulations are concentrated on clarifying the universal aspects of polymer crystallization using model linear polymers such as polyethylene. We are recently focusing on a nearly untouched but very interesting problem of chiral selecting crystallization in helical polymers. We previously proposed a stepwise approach using two kinds of helical polymers, simple "bare" helical polymers made of backbone atoms only such as polyoxymethylene (POM) and "general" helical polymers containing complicated side groups such as isotactic polypropylene. We have already reported on the crystallization in oligomeric POM-like helix but have observed only weak chiral selectivity during crystallization. In the present paper, we investigate the crystallization of sufficiently long POM-like polymer both from the isotropic melt and from the highly stretched melt. We find in both cases that the polymer shows a clear chiral selecting crystallization. Specifically, the observation of a single crystal growing from the isotropic melt is very illuminating. It shows that the crystal thickness and the crystal chirality are closely correlated; thicker crystals show definite chirality while thinner ones are mostly mixtures of the R- and the L-handed stems. The single crystal is found to have a marked lenticular shape, where the thinner growth front, since being made of the mixture, shows no chiral selectivity. The final chiral crystal is found to be completed through helix reversal processes within thicker regions. [ABSTRACT FROM AUTHOR]

Finding a low dimensional representation of data from long-timescale trajectories of biomolecular processes, such as protein folding or ligand–receptor binding, is of fundamental importance, and kinetic models, such as Markov modeling, have proven useful in describing the kinetics of these systems. Recently, an unsupervised machine learning technique called VAMPNet was introduced to learn the low dimensional representation and the linear dynamical model in an end-to-end manner. VAMPNet is based on the variational approach for Markov processes and relies on neural networks to learn the coarse-grained dynamics. In this paper, we combine VAMPNet and graph neural networks to generate an end-to-end framework to efficiently learn high-level dynamics and metastable states from the long-timescale molecular dynamics trajectories. This method bears the advantages of graph representation learning and uses graph message passing operations to generate an embedding for each datapoint, which is used in the VAMPNet to generate a coarse-grained dynamical model. This type of molecular representation results in a higher resolution and a more interpretable Markov model than the standard VAMPNet, enabling a more detailed kinetic study of the biomolecular processes. Our GraphVAMPNet approach is also enhanced with an attention mechanism to find the important residues for classification into different metastable states. [ABSTRACT FROM AUTHOR]

Hybrid methods that combine molecular dynamics methods capable of analyzing dynamics with Monte Carlo (MC) methods that can efficiently treat thermodynamically stable states are valuable for understanding complex chemical processes in which an equilibrium state is reached through many elementary processes. The hybrid MC (HMC) method is one such promising method; however, it often fails to sample configurations properly from the canonical multimodal distribution due to the rugged potential energy surfaces. In this paper, we extend the HMC method to overcome this difficulty. The new method, which is termed potential scaling HMC (PS-HMC), makes use of an artificially modulated trajectory to propose a new configuration. The trajectory is generated from Hamilton's equations, but the potential energy surface is scaled to be gradually flattened and then recovered to the original surface, which facilitates barrier-crossing processes. We apply the PS-HMC method to three kinds of molecular processes: the thermal motion of argon particles, butane isomerization, and an atom transfer chemical reaction. These applications demonstrate that the PS-HMC method is capable of correctly constructing the canonical ensemble with a multimodal distribution. The sampling efficiency and accepted trajectories are examined to clarify the features of the PS-HMC method. Despite the potential scaling, many reactive atom transfer trajectories (elementary processes) pass through the vicinity of the minimum energy path. Furthermore, we demonstrate that the method can properly imitate the relaxation process owing to the inherent configurational continuity. By comparing the PS-HMC method with other relevant methods, we can conclude that the new method is a unique approach for studying both the dynamic and thermodynamic aspects of chemical processes. [ABSTRACT FROM AUTHOR]