Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userThermodiffusion in model nanofluids by molecular dynamics simulations
167
0
0.0
(
0
)
Authors
Guillaume Galliero
Ask ChatGPT about the research
No Arabic abstract
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.
rate research
Read More
One popular approach to incorporating experimental data into molecular simulations is to restrain the ensemble average of observables to their experimental values. Here I derive equations for the equilibrium distributions generated by restrained ensemble simulations and the corresponding expected values of observables. My results suggest a method to restrain simulations so that they generate distributions that are minimally perturbed from the unbiased distributions while reproducing the experimental values of the observables within their measurement uncertainties.
We present a molecular dynamics simulation method for the computation of the solubility of organic crystals in solution. The solubility is calculated based on the equilibrium free energy difference between the solvated solute and its crystallized state at the crystal surface kink site. In order to efficiently sample the growth and dissolution process, we have carried out well-tempered Metadynamics simulations with a collective variable that captures the slow degrees of freedom, namely the solute diffusion to and adsorption at the kink site together with the desolvation of the kink site. Simulations were performed at different solution concentrations using constant chemical potential molecular dynamics and the solubility was identified at the concentration at which the free energy values between the grown and dissolved kink states were equal. The effectiveness of this method is demonstrated by its success in reproducing the experimental trends of solubility of urea and naphthalene in a variety of solvents.
In the replica-exchange molecular dynamics method, where constant-temperature molecular dynamics simulations are performed in each replica, one usually rescales the momentum of each particle after replica exchange. This rescaling method had previously been worked out only for the Gaussian constraint method. In this letter, we present momentum rescaling formulae for four other commonly used constant-temperature algorithms, namely, Langevin dynamics, Andersen algorithm, Nos{e}-Hoover thermostat, and Nos{e}-Poincar{e} thermostat. The effectiveness of these rescaling methods is tested with a small biomolecular system, and it is shown that proper momentum rescaling is necessary to obtain correct results in the canonical ensemble.
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.
We introduce a scheme for deriving an optimally-parametrised Langevin dynamics of few collective variables from data generated in molecular dynamics simulations. The drift and the position-dependent diffusion profiles governing the Langevin dynamics are expressed as explicit averages over the input trajectories. The proposed strategy is applicable to cases when the input trajectories are generated by subjecting the system to a external time-dependent force (as opposed to canonically-equilibrated trajectories). Secondly, it provides an explicit control on the statistical uncertainty of the drift and diffusion profiles. These features lend to the possibility of designing the external force driving the system so to maximize the accuracy of the drift and diffusions profile throughout the phase space of interest. Quantitative criteria are also provided to assess a posteriori the satisfiability of the requisites for applying the method, namely the Markovian character of the stochastic dynamics of the collective variables.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
