166 research outputs found
Local density dependent potential for compressible mesoparticles
We focus on finding a coarse grained description able to reproduce the
thermodynamic behavior of a molecular system by using mesoparticles
representing several molecules. Interactions between mesoparticles are modelled
by an interparticle potential, and an additional internal equation of state is
used to account for the thermic contribution of coarse grained internal degrees
of freedom. Moreover, as strong non-equilibrium situations over a wide range of
pressure and density are targeted, the internal compressibility of these
mesoparticles has to be considered. This is done by introducing a dependence of
the potential on the local environment of the mesoparticles, either by defining
a spherical local density or by means of a Voronoi tessellation. As an example,
a local density dependent potential is fitted to reproduce the Hugoniot curve
of a model of nitromethane, where each mesoparticle represents one thousand
molecules
Microscopic calculations of Hugoniot curves of neat TATB and of its detonation products
We compute the Hugoniot curves of both neat TATB and its detonation products
mixture using atomistic simulation tools. To compute the Hugoniot states, we
adapted our "Sampling Constraints in Average" (SCA) method (Maillet et al.,
Applied Math. Research eXpress 2008, 2009) to Monte-Carlo simulations. For neat
TATB, we show that the potential proposed by Rai (Rai et al., J. Chem. Phys.
129, 2008) is not accurate enough to predict the Hugoniot curve and requires
some optimization of its parameters. Concerning detonation products,
thermodynamic properties at chemical equilibrium are computed using a specific
RxMC method (Bourasseau et al., Phys. Chem. Chem. Phys. 13, 2011) taking into
account the presence of carbon clusters in the fluid mixture. We show that this
explicit description of the solid phase immersed in the fluid phase modifies
the chemical equilibrium
Permutation-invariant distance between atomic configurations
We present a permutation-invariant distance between atomic configurations,
defined through a functional representation of atomic positions. This distance
enables to directly compare different atomic environments with an arbitrary
number of particles, without going through a space of reduced dimensionality
(i.e. fingerprints) as an intermediate step. Moreover, this distance is
naturally invariant through permutations of atoms, avoiding the time consuming
associated minimization required by other common criteria (like the Root Mean
Square Distance). Finally, the invariance through global rotations is accounted
for by a minimization procedure in the space of rotations solved by Monte Carlo
simulated annealing. A formal framework is also introduced, showing that the
distance we propose verifies the property of a metric on the space of atomic
configurations. Two examples of applications are proposed. The first one
consists in evaluating faithfulness of some fingerprints (or descriptors), i.e.
their capacity to represent the structural information of a configuration. The
second application concerns structural analysis, where our distance proves to
be efficient in discriminating different local structures and even classifying
their degree of similarity
New Potential Model for Molecular Dynamic Simulation of liquid HF. II -Parameter Optimization for Repulsion-Dispersion potential
International audienceIn order to build a complete potential model to perform classical molecular dynamic simulations of liquid HF, a new optimization method is proposed to obtain transferable parameters for repulsion-dispersion potential on the basis of ab initio reference data. This process is decomposed into two steps. The first step, using the force-matching method, consists in exploring the parameter space and selecting a first potential used as a start point for the second step. This last step consists in optimizing the parameters of the selected potential in order to reproduce reference thermodynamic and structural data. The obtained potential correctly reproduces the radial distribution functions and the pressures of HF liquid over a large range of thermodynamic states
The Sine-Gordon Solitons as a N-Body Problem
We consider the N-soliton solutions in the sine-Gordon model as a N-body
problem. This leads to a relativistic generalization of the Calogero model
first introduced by Ruijsenaars. We show that the fundamental Poisson bracket
of the Lax matrix is quadratic, and the -matrix is a dynamical one. This is
in contrast to the Calogero model where the fundamental Poisson bracket of the
Lax matrix is linear.Comment: 10 pages LATEX SPhT-93-072; LPTHE-93-4
Molecular Simulations of Hugoniots of detonation products mixtures at chemical equilibrium: Microscopic calculation of the Chapman-Jouguet State
International audienceIn this work, we used simultaneously the Reaction Ensemble Monte Carlo (ReMC) method and the Adaptive Erpenbeck Equation Of State (AE-EOS) method to directly calculate the thermodynamical and chemical equilibrium of mixtures of detonation products on the Hugoniot curve. The ReMC method (W. R. Smith and B. Triska, J. Chem. Phys. 100, pp 3019-3027 (1994)) allows to reach the chemical equilibrium of a reacting mixture, and the AE-EOS method (J. J. Erpenbeck, Phys. Rev. A, 46, p 6406 (1992)) constrains the system to satisfy the Hugoniot relation. Once the Hugoniot curve of the detonation products mixture is established, the CJ state of the explosive can be determined. Performing a NPT simulation at P(CJ) , T(CJ) , we then calculate the direct thermodynamic properties and the following derivative properties of the system using a fluctuation method: calorific capacities, sound velocity and Gruneisen coefficient. As the composition fluctuates, and the number of particles is not necessarily constant in this ensemble, a fluctuation formula has been developed to take into account the fluctuations of mole number and composition. This type of calculation has been applied to several usual energetic materials: nitromethane, tetranitromethane, hexanitroethane, PETN and RDX
Molecular Simulations of Shock to Detonation Transition in Nitromethane
An extension of the model described in a previous work of Maillet, Soulard
and Stoltz based on a Dissipative Particule Dynamics is presented and applied
to liquid nitromethane. Large scale non-equilibrium simulations of reacting
nitromethane under sustained shock conditions allow a better understanding of
the shock-to-detonation transition in homogeneous explosives. Moreover, the
propagation of the reactive wave appears discontinuous since ignition points in
the shocked material can be activated by the compressive waves emitted from the
onset of chemical reactions
Neighbors Map: an Efficient Atomic Descriptor for Structural Analysis
Accurate structural analysis is essential to gain physical knowledge and
understanding of atomic-scale processes in materials from atomistic
simulations. However, traditional analysis methods often reach their limits
when applied to crystalline systems with thermal fluctuations, defect-induced
distortions, partial vitrification, etc. In order to enhance the means of
structural analysis, we present a novel descriptor for encoding atomic
environments into 2D images, based on a pixelated representation of graph-like
architecture with weighted edge connections of neighboring atoms. This
descriptor is well adapted for Convolutional Neural Networks and enables
accurate structural analysis at a low computational cost. In this paper, we
showcase a series of applications, including the classification of crystalline
structures in distorted systems, tracking phase transformations up to the
melting temperature, and analyzing liquid-to-amorphous transitions in pure
metals and alloys. This work provides the foundation for robust and efficient
structural analysis in materials science, opening up new possibilities for
studying complex structural processes, which can not be described with
traditional approaches
- …