No Arabic abstract
We introduce a method to obtain one-dimensional collective variables for studying rarely occurring transitions between two metastable states separated by a high free energy barrier. No previous information, not even approximated, on the path followed during the transition is needed. The only requirement is to know the fluctuations of the system while in the two metastable states. With this information in hand we build the collective variable using a modified version of Fishers linear discriminant analysis. The usefulness of this approach is tested on the metadynamics simulation of two representative systems. The first is the freezing of silver iodide into the superionic $alpha$-phase, the second is the study of a classical Diels Alder reaction. The collective variable works very well in these two diverse cases.
We propose a rigorous construction of a 1D path collective variable to sample structural phase transformations in condensed matter. The path collective variable is defined in a space spanned by global collective variables that serve as classifiers derived from local structural units. A reliable identification of local structural environments is achieved by employing a neural network based classification. The 1D path collective variable is subsequently used together with enhanced sampling techniques to explore the complex migration of a phase boundary during a solid-solid phase transformation in molybdenum.
Computing accurate reaction rates is a central challenge in computational chemistry and biology because of the high cost of free energy estimation with unbiased molecular dynamics. In this work, a data-driven machine learning algorithm is devised to learn collective variables with a multitask neural network, where a common upstream part reduces the high dimensionality of atomic configurations to a low dimensional latent space, and separate downstream parts map the latent space to predictions of basin class labels and potential energies. The resulting latent space is shown to be an effective low-dimensional representation, capturing the reaction progress and guiding effective umbrella sampling to obtain accurate free energy landscapes. This approach is successfully applied to model systems including a 5D Muller Brown model, a 5D three-well model, and alanine dipeptide in vacuum. This approach enables automated dimensionality reduction for energy controlled reactions in complex systems, offers a unified framework that can be trained with limited data, and outperforms single-task learning approaches, including autoencoders.
We present extensive new textit{ab intio} path integral Monte Carlo results for the momentum distribution function $n(mathbf{k})$ of the uniform electron gas (UEG) in the warm dense matter (WDM) regime over a broad range of densities and temperatures. This allows us to study the nontrivial exchange--correlation induced increase of low-momentum states around the Fermi temperature, and to investigate its connection to the related lowering of the kinetic energy compared to the ideal Fermi gas. In addition, we investigate the impact of quantum statistics on both $n(mathbf{k})$ and the off-diagonal density matrix in coordinate space, and find that it cannot be neglected even in the strongly coupled electron liquid regime. Our results were derived without any nodal constraints, and thus constitute a benchmark for other methods and approximations.
The computational study of conformational transitions in RNA and proteins with atomistic molecular dynamics often requires suitable enhanced sampling techniques. We here introduce a novel method where concurrent metadynamics are integrated in a Hamiltonian replica-exchange scheme. The ladder of replicas is built with different strength of the bias potential exploiting the tunability of well-tempered metadynamics. Using this method, free-energy barriers of individual collective variables are significantly reduced compared with simple force-field scaling. The introduced methodology is flexible and allows adaptive bias potentials to be self-consistently constructed for a large number of simple collective variables, such as distances and dihedral angles. The method is tested on alanine dipeptide and applied to the difficult problem of conformational sampling in a tetranucleotide.
We present extensive new emph{ab initio} path integral Monte Carlo (PIMC) simulations of normal liquid $^3$He without any nodal constraints. This allows us to study the effects of temperature on different structural properties like the static structure factor $S(mathbf{q})$, the momentum distribution $n(mathbf{q})$, and the static density response function $chi(mathbf{q})$, and to unambiguously quantify the impact of Fermi statistics. In addition, the dynamic structure factor $S(mathbf{q},omega)$ is rigorously reconstructed from imaginary-time PIMC data, and we find the familiar phonon-maxon-roton dispersion that is well known from $^4$He and has been reported previously for two-dimensional $^3$He films [Nature textbf{483}, 576-579 (2012)]. The comparison of our new results for both $S(mathbf{q})$ and $S(mathbf{q},omega)$ to neutron scattering measurements reveals an excellent agreement between theory and experiment.