Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userHamiltonian adaptive resolution simulation for molecular liquids
145
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Adaptive resolution schemes allow the simulation of a molecular fluid treating simultaneously different subregions of the system at different levels of resolution. In this work we present a new scheme formulated in terms of a global Hamiltonian. Within this approach equilibrium states corresponding to well defined statistical ensembles can be generated making use of all standard Molecular Dynamics or Monte Carlo methods. Models at different resolutions can thus be coupled, and thermodynamic equilibrium can be modulated keeping each region at desired pressure or density without disrupting the Hamiltonian framework.
rate research
Read More
Path integral-based simulation methodologies play a crucial role for the investigation of nuclear quantum effects by means of computer simulations. However, these techniques are significantly more demanding than corresponding classical simulations. To reduce this numerical effort, we recently proposed a method, based on a rigorous Hamiltonian formulation, which restricts the quantum modeling to a small but relevant spatial region within a larger reservoir where particles are treated classically. In this work, we extend this idea and show how it can be implemented along with state-of-the-art path integral simulation techniques, such as ring polymer and centroid molecular dynamics, which allow the approximate calculation of both quantum statistical and quantum dynamical properties. To this end, we derive a new integration algorithm which also makes use of multiple time-stepping. The scheme is validated via adaptive classical--path-integral simulations of liquid water. Potential applications of the proposed multiresolution method are diverse and include efficient quantum simulations of interfaces as well as complex biomolecular systems such as membranes and proteins.
In stationary nonequilibrium states coupling between hydrodynamic modes causes thermal fluctuations to become long ranged inducing nonequilibrium Casimir pressures. Here we consider nonequilibrium Casimir pressures induced in liquids by a velocity gradient. Specifically, we have obtained explicit expressions for the magnitude of the shear-induced pressure enhancements in a liquid layer between two horizontal plates that complete and correct results previously presented in the literature. In contrast to nonequilibrium Casimir pressures induced by a temperature or concentration gradient, we find that in shear nonequilibrium contributions from short-range fluctuations are no longer negligible. In addition, it is noted that currently available computer simulations of model fluids in shear observe effects from molecular correlations at nanoscales that have a different physical origin and do not probe shear-induced pressures resulting from coupling of long-wavelength hydrodynamic modes. Even more importantly, we find that in actual experimental conditions, shear-induced pressure enhancements are caused by viscous heating and not by thermal velocity fluctuations. Hence, isothermal computer simulations are irrelevant for the interpretation of experimental shear-induced pressure enhancements.
In this thesis we have used Quantum Monte Carlo techniques to study two systems that can be regarded as the archetype for neutral strongly interacting systems: 4He, and its fermionic counterpart 3He.More specifically, we have used the Path Integral Ground State and the Path Integral Monte Carlo methods to study a system of two dimensional 3He (2d-3He) and a system of 4He adsorbed on Graphene-Fluoride (GF) and Graphane (GH) at both zero and finite temperature. The purpose of the study of 4He on GF (GH) was the research of new physical phenomena, whereas in the case of 2d-3He it was the application of novel methodologies for the ab-initio study of static and dynamic properties of Fermi systems. In the case of 2d-3He we have computed the spin susceptibility as function of density which turned out to be in very good agreement with experimental data; we have also obtained the first ab-initio evaluation of the zero-sound mode and the dynamic structure factor of 2d-3He that is in remarkably good agreement with experiments. In the case of 4He adsorbed on GF (GH), we determined the zero temperature equilibrium density of the first monolayer of 4He showing also that the commensurate sqrt(3) x sqrt(3) R30 phase is unstable on both substrates; at equilibrium density we found that 4He on GF (GH) is a modulated superfluid with an anisotropic phono-rotonic spectrum; at high coverages we found an incommensurate triangular solid and, on both GF and GH, a commensurate phase at filling factor x= 2/7 that is locally stable or at least metastable. Remarkably, in this commensurate solid phase and for both GF and GH, our computations show preliminary evidence of the presence of a superfluid fraction.
We discuss the use of a Langevin equation with a colored (correlated) noise to perform constant-temperature molecular dynamics simulations. Since the equations of motion are linear in nature, it is easy to predict the response of a Hamiltonian system to such a thermostat and to tune at will the relaxation time of modes of different frequency. This allows one to optimize the time needed to thermalize the system and generate independent configurations. We show how this frequency-dependent response can be exploited to control the temperature of Car-Parrinello-like dynamics, keeping at low temperature the electronic degrees of freedom, without affecting the adiabatic separation from the vibrations of the ions.
It was recently demonstrated that a simple Monte Carlo (MC) algorithm involving the swap of particle pairs dramatically accelerates the equilibrium sampling of simulated supercooled liquids. We propose two numerical schemes integrating the efficiency of particle swaps into equilibrium molecular dynamics (MD) simulations. We first develop a hybrid MD/MC scheme combining molecular dynamics with the original swap Monte Carlo. We implement this hybrid method in LAMMPS, a software package employed by a large community of users. Secondly, we define a continuous time version of the swap algorithm where both the positions and diameters of the particles evolve via Hamiltons equations of motion. For both algorithms, we discuss in detail various technical issues as well as the optimisation of simulation parameters. We compare the numerical efficiency of all available swap algorithms and discuss their relative merits.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
