Time Recovery for a Complex Process Using Accelerated Dynamics

Fractional Extended Diffusion Theory to capture anomalous relaxation from biased/accelerated molecular simulations

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.

Diversity of Behavior after Collisions of Sn and Si Nanoparticles Found Using a New Density Functional Tight-Binding Method

We present a new approach to studying nanoparticle collisions using density functional based tight binding (DFTB). A novel DFTB parametrization has been developed to study the collision process of Sn and Si clusters (NPs) using molecular dynamics (MD). While bulk structures were used as training sets, we show that our model is able to accurately reproduce the cohesive energy of the nanoparticles using density functional theory (DFT) as a reference. A surprising variety of phenomena are revealed for the Si/Sn nanoparticle collisions, depending on the size and velocity of the collision: from core-shell structure formation to bounce-off phenomena.

Characterization of amorphous LixSi structures from ReaxFF via accelerated exploration of local minima

Motivated by the abundant experimental work in the area of Li-ion batteries, in the present work we characterize via computer simulations the structure of Si-Li amorphous alloys in a wide range of compositions. Using a reactive force field we propose a novel accelerated exploration of local minima to obtain amorphous structures close to equilibrium. The features of this system analyzed for different alloy compositions are the partial radial distribution functions g(r), the first and second nearest neighbour coordination numbers and the short-order structure. The complex structure of the second peak of the Si-Li g(r) is elucidated using a cluster-connection analysis.