Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userDefining velocities for accurate kinetic statistics in the GJF thermostat
192
0
0.0
(
0
)
Added by
Niels Gronbech-Jensen
Publication date
2019
fields
Physics
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
We expand on two previous developments in the modeling of discrete-time Langevin systems. One is the well-documented Gr{o}nbech-Jensen Farago (GJF) thermostat, which has been demonstrated to give robust and accurate configurational sampling of the phase space. Another is the recent discovery that also kinetics can be accurately sampled for the GJF method. Through a complete investigation of all possible finite difference approximations to the velocity, we arrive at two main conclusions:~1) It is not possible to define a so-called on-site velocity such that kinetic temperature will be correct and independent of the time step, and~2) there exists a set of infinitely many possibilities for defining a two-point (leap-frog) velocity that measures kinetic energy correctly for linear systems in addition to the correct configurational statistics obtained from the GJF algorithm. We give explicit expressions for the possible definitions, and we incorporate these into convenient and practical algorithmic forms of the normal Verlet-type algorithms along with a set of suggested criteria for selecting a useful definition of velocity.
rate research
Read More
We make a brief historical review to the moment model reduction to the kinetic equations, particularly the Grads moment method for Boltzmann equation. The focus is on the hyperbolicity of the reduced model, which is essential to the existence of its classical solution as a Cauchy problem. The theory of the framework we developed in last years is then introduced, which may preserve the hyperbolic nature of the kinetic equations with high universality. Some lastest progress on the comparison between models with/without hyperbolicity is presented to validate the hyperbolic moment models for rarefied gases.
Entropy production during the process of thermal phase-separation of multiphase flows is investigated by means of a discrete Boltzmann kinetic model. The entropy production rate is found to increase during the spinodal decomposition stage and to decrease during the domain growth stage, attaining its maximum at the crossover between the two. Such behaviour provides a natural criterion to identify and discriminate between the two regimes. Furthermore, the effects of heat conductivity, viscosity and surface tension on the entropy production rate are investigated by systematically probing the interplay between non-equilibrium energy and momentum fluxes. It is found that the entropy production rate due to energy fluxes is an increasing function of the Prandtl number, while the momentum fluxes exhibit an opposite trend. On the other hand, both contributions show an increasing trend with surface tension. The present analysis inscribes within the general framework of non-equilibrium thermodynamics and consequently it is expected to be relevant to a broad class of soft-flowing systems far from mechanical and thermal equilibrium.
In light of the recently developed complete GJ set of single random variable stochastic, discrete-time St{o}rmer-Verlet algorithms for statistically accurate simulations of Langevin equations, we investigate two outstanding questions: 1) Are there any algorithmic or statistical benefits from including multiple random variables per time-step, and 2) are there objective reasons for using one or more methods from the available set of statistically correct algorithms? To address the first question, we assume a general form for the discrete-time equations with two random variables and then follow the systematic, brute-force GJ methodology by enforcing correct thermodynamics in linear systems. It is concluded that correct configurational Boltzmann sampling of a particle in a harmonic potential implies correct configurational free-particle diffusion, and that these requirements only can be accomplished if the two random variables per time step are identical. We consequently submit that the GJ set represents all possible stochastic St{o}rmer-Verlet methods that can reproduce time-step-independent statistics of linear systems. The second question is thus addressed within the GJ set. Based in part on numerical simulations of complex molecular systems, and in part on analytic scaling of time, we analyze the apparent difference in stability between different methods. We attribute this difference to the inherent time scaling in each method, and suggest that this scaling may lead to inconsistencies in the interpretation of dynamical and statistical simulation results. We therefore suggest that the method with the least inherent time-scaling, the GJ-I/GJF-2GJ method, be preferred for statistical applications where spurious rescaling of time is undesirable.
To study epitaxial thin-film growth, a new model is introduced and extensive kinetic Monte Carlo simulations performed for a wide range of fluxes and temperatures. Varying the deposition conditions, a rich growth diagram is found. The model also reproduces several known regimes and in the limit of low particle mobility a new regime is defined. Finally, a relation is postulated between the temperatures of the kinetic and thermal roughening transitions.
We expand on the previously published Gr{o}nbech-Jensen Farago (GJF) thermostat, which is a thermodynamically sound variation on the St{o}rmer-Verlet algorithm for simulating discrete-time Langevin equations. The GJF method has been demonstrated to give robust and accurate configurational sampling of the phase space, and its applications to, e.g., Molecular Dynamics is well established. A new definition of the discrete-time velocity variable is proposed based on analytical calculations of the kinetic response of a harmonic oscillator subjected to friction and noise. The new companion velocity to the GJF method is demonstrated to yield correct and time-step-independent kinetic responses for, e.g., kinetic energy, its fluctuations, and Green-Kubo diffusion based on velocity autocorrelations. This observation allows for a new and convenient Leap-Frog algorithm, which efficiently and precisely represents statistical measures of both kinetic and configurational properties at any time step within the stability limit for the harmonic oscillator. We outline the simplicity of the algorithm and demonstrate its attractive time-step-independent features for nonlinear and complex systems through applications to a one-dimensional nonlinear oscillator and three-dimensional Molecular Dynamics.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
