Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userMass-zero constrained molecular dynamics for electrode charges in simulations of electrochemical systems
71
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Classical molecular dynamics simulations have recently become a standard tool for the study of electrochemical systems. State-of-the-art approaches represent the electrodes as perfect conductors, modelling their responses to the charge distribution of electrolytes via the so-called fluctuating charge model. These fluctuating charges are additional degrees of freedom that, in a Born-Oppenheimer spirit, adapt instantaneously to changes in the environment to keep each electrode at a constant potential. Here we show that this model can be treated in the framework of constrained molecular dynamics, leading to a symplectic and time-reversible algorithm for the evolution of all the degrees of freedom of the system. The computational cost and the accuracy of the new method are similar to current alternative implementations of the model. The advantage lies in the accuracy and long term stability guaranteed by the formal properties of the algorithm and in the possibility to systematically introduce additional kinematic conditions of arbitrary number and form. We illustrate the performance of the constrained dynamics approach by enforcing the electroneutrality of the electrodes in a simple capacitor consisting of two graphite electrodes separated by a slab of liquid water.
rate research
Read More
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.
Efficient implementations of the classical molecular dynamics (MD) method for Lennard-Jones particle systems are considered. Not only general algorithms but also techniques that are efficient for some specific CPU architectures are also explained. A simple spatial-decomposition-based strategy is adopted for parallelization. By utilizing the developed code, benchmark simulations are performed on a HITACHI SR16000/J2 system consisting of IBM POWER6 processors which are 4.7 GHz at the National Institute for Fusion Science (NIFS) and an SGI Altix ICE 8400EX system consisting of Intel Xeon processors which are 2.93 GHz at the Institute for Solid State Physics (ISSP), the University of Tokyo. The parallelization efficiency of the largest run, consisting of 4.1 billion particles with 8192 MPI processes, is about 73% relative to that of the smallest run with 128 MPI processes at NIFS, and it is about 66% relative to that of the smallest run with 4 MPI processes at ISSP. The factors causing the parallel overhead are investigated. It is found that fluctuations of the execution time of each process degrade the parallel efficiency. These fluctuations may be due to the interference of the operating system, which is known as OS Jitter.
The quantum many-body problem in condensed phases is often simplified using a quasiparticle description, such as effective mass theory for electron motion in a periodic solid. These approaches are often the basis for understanding many fundamental condensed phase processes, including the molecular mechanisms underlying solar energy harvesting and photocatalysis. Despite the importance of these effective particles, there is still a need for computational methods that can explore their behavior on chemically relevant length and time scales. This is especially true when the interactions between the particles and their environment are important. We introduce an approach for studying quasiparticles in condensed phases by combining effective mass theory with the path integral treatment of quantum particles. This framework incorporates the generally anisotropic electronic band structure of materials into path integral simulation schemes to enable modeling of quasiparticles in quantum confinement, for example. We demonstrate the utility of effective mass path integral simulations by modeling an exciton in solid potassium chloride and electron trapping by a sulfur vacancy in monolayer molybdenum disulfide.
In this work, a new algorithm is proposed to compute single particle (infinite dilution) thermodiffusion using Non-Equilibrium Molecular Dynamics simulations through the estimation of the thermophoretic force that applies on a solute particle. This scheme is shown to provide consistent results for simple Lennard-Jones fluids and for model nanofluids (spherical non-metallic nanoparticles + Lennard-Jones fluid) where it appears that thermodiffusion amplitude, as well as thermal conductivity, decrease with nanoparticles concentration. Then, in nanofluids in the liquid state, by changing the nature of the nanoparticle (size, mass and internal stiffness) and of the solvent (quality and viscosity) various trends are exhibited. In all cases the single particle thermodiffusion is positive, i.e. the nanoparticle tends to migrate toward the cold area. The single particle thermal diffusion 2 coefficient is shown to be independent of the size of the nanoparticle (diameter of 0.8 to 4 nm), whereas it increases with the quality of the solvent and is inversely proportional to the viscosity of the fluid. In addition, this coefficient is shown to be independent of the mass of the nanoparticle and to increase with the stiffness of the nanoparticle internal bonds. Besides, for these configurations, the mass diffusion coefficient behavior appears to be consistent with a Stokes-Einstein like law.
Recent experiments have observed that the chemical and photophysical properties of molecules can be modified inside an optical Fabry-Perot microcavity under collective vibrational strong coupling (VSC) conditions, and such modification is currently not well understood by theory. In an effort to understand the origin of such cavity induced phenomena, some recent studies have focused on the effect of the cavity environment on the nonlinear optical response of the molecular subsystem. Here, we use a recently proposed protocol for classical cavity molecular dynamics (CavMD) simulations to numerically investigate the linear and nonlinear response of liquid carbon dioxide under such VSC conditions following an optical pulse excitation. We find that applying a strong pulse of excitation to the lower hybrid light-matter state, i.e., the lower polariton (LP), can lead to an overall molecular nonlinear absorption which is enhanced by up to two orders of magnitude relative to the excitation outside the cavity. This polariton-enhanced multiphoton absorption also causes an ultrashort LP lifetime (0.2 ps) under strong illumination. Unlike usual polariton relaxation processes -- whereby polaritonic energy transfers directly to the manifold of singly excited vibrational dark states -- under the present mechanism, the LP transfers energy directly to the manifold of higher vibrationally excited dark states; these highly excited dark states subsequently relax to the manifold of singly excited states with a lifetime of tens of ps. Because the present mechanism is generic in nature, we expect these numerical predictions to be experimentally observed in different molecular systems and in cavities with different volumes.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
