Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userSimulating Met-Enkephalin With Population Annealing Molecular Dynamics
333
0
0.0
(
0
)
Added by
Henrik Christiansen
Publication date
2020
fields
Physics
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
Met-enkephalin, one of the smallest opiate peptides and an important neurotransmitter, is a widely used benchmarking problem in the field of molecular simulation. Through its range of possible low-temperature conformations separated by free-energy barriers it was previously found to be hard to thermalize using straight canonical molecular dynamics simulations. Here, we demonstrate how one can use the recently proposed population annealing molecular dynamics scheme to overcome these difficulties. We show how the use of multi-histogram reweighting allows one to accurately estimate the density of states of the system and hence derive estimates such as the potential energy as quasi continuous functions of temperature. We further investigate the free-energy surface as a function of end-to-end distance and radius-of-gyration and observe two distinct basins of attraction.
rate research
Read More
Physics students now have access to interactive molecular dynamics simulations that can model and animate the motions of hundreds of particles, such as noble gas atoms, that attract each other weakly at short distances but repel strongly when pressed together. Using these simulations, students can develop an understanding of forces and motions at the molecular scale, nonideal fluids, phases of matter, thermal equilibrium, nonequilibrium states, the Boltzmann distribution, the arrow of time, and much more. This article summarizes the basic features and capabilities of such a simulation, presents a variety of student exercises using it at the introductory and intermediate levels, and describes some enhancements that can further extend its uses. A working simulation code, in HTML5 and JavaScript for running within any modern Web browser, is provided as an online supplement.
Driving an inertial many-body system out of equilibrium generates complex dynamics due to memory effects and the intricate relationships between the external driving force, internal forces, and transport effects. Understanding the underlying physics is challenging and often requires carrying out case-by-case analysis. To systematically study the interplay between all types of forces that contribute to the dynamics, a method to generate prescribed flow patterns could be of great help. We develop a custom flow method to numerically construct the external force field required to obtain the desired time evolution of an inertial many-body system, as prescribed by its one-body current and density profiles. We validate the custom flow method in a Newtonian system of purely repulsive particles by creating a slow motion dynamics of an out-of-equilibrium process and by prescribing the full time evolution between two distinct equilibrium states. The method can also be used with thermostat algorithms to control the temperature.
Population annealing is a recent addition to the arsenal of the practitioner in computer simulations in statistical physics and beyond that is found to deal well with systems with complex free-energy landscapes. Above all else, it promises to deliver unrivaled parallel scaling qualities, being suitable for parallel machines of the biggest calibre. Here we study population annealing using as the main example the two-dimensional Ising model which allows for particularly clean comparisons due to the available exact results and the wealth of published simulational studies employing other approaches. We analyze in depth the accuracy and precision of the method, highlighting its relation to older techniques such as simulated annealing and thermodynamic integration. We introduce intrinsic approaches for the analysis of statistical and systematic errors, and provide a detailed picture of the dependence of such errors on the simulation parameters. The results are benchmarked against canonical and parallel tempering simulations.
The dynamical arrest of gels is the consequence of a well defined structural phase transition, leading to the formation of a spanning cluster of bonded particles. The dynamical glass transition, instead, is not accompanied by any clear structural signature. Nevertheless, both transitions are characterized by the emergence of dynamical heterogeneities. Reviewing recent results from numerical simulations, we discuss the behavior of dynamical heterogeneities in different systems and show that a clear connection with the structure exists in the case of gels. The emerging picture may be also relevant for the more elusive case of glasses. We show, as an example, that the relaxation process of a simple glass-forming model can be related to a reverse percolation transition and discuss further perspective in this direction.
In this work we investigate the transient solidification of a Lennard-Jones liquid using non-equilibrium molecular dynamics simulations and continuum heat transfer theory. The simulations are performed in slab-shaped boxes, where a cold thermostat placed at the centre of the box drives the solidification of the liquid. Two well-defined solid fronts propagate outwards from the centre towards the ends of the box until solidification is completed. A continuum phase change model that accounts for the difference between the solid and the liquid densities is formulated to describe the evolution of the temperature and the position of the solidification front. Simulation results for a small and a large nanoscale system, of sizes $30.27$,nm and $60.54$,nm, are compared with the predictions of the theoretical model. Following a transient period of $sim$20-40 ps and a displacement of the solidification front of 1-2.5 nm we find that the simulations and the continuum theory show good agreement. We use this fact to combine the simulation and theoretical approaches to design a simple procedure to calculate the latent heat of the material. We also perform simulations of the homogeneous freezing process, i.e. in the absence of a temperature gradient and at constant temperature, by quenching the liquid at supercooled temperatures. We demonstrate that the solidification rate of homogenous freezing is much faster than the one obtained under a thermal gradient for systems of the same size subject to the same thermostat temperature. Our study and conclusions should be of general interest to a wide range of atomistic solids.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
