atomistic calculations of dynamic compression of materialso semiclassical quantum nuclear effects in...
TRANSCRIPT
![Page 1: Atomistic Calculations of Dynamic Compression of Materialso Semiclassical quantum nuclear effects in classical molecular dynamics: An initial approach . How does a shock wave in a](https://reader033.vdocument.in/reader033/viewer/2022042319/5f087bc07e708231d4223c27/html5/thumbnails/1.jpg)
Atomistic Calculations of Dynamic Compression of Materials
July 2, 2013
Tingting Qi, Evan Reed
Department of Materials Science and Engineering Stanford University
![Page 2: Atomistic Calculations of Dynamic Compression of Materialso Semiclassical quantum nuclear effects in classical molecular dynamics: An initial approach . How does a shock wave in a](https://reader033.vdocument.in/reader033/viewer/2022042319/5f087bc07e708231d4223c27/html5/thumbnails/2.jpg)
OUTLINE
o What’s a shock wave and how is it used?
o Tractable atomistic shock simulations (Multi-scale Shock Technique)
o When classical MD isn’t good enough o Semiclassical quantum nuclear effects in classical
molecular dynamics: An initial approach
![Page 3: Atomistic Calculations of Dynamic Compression of Materialso Semiclassical quantum nuclear effects in classical molecular dynamics: An initial approach . How does a shock wave in a](https://reader033.vdocument.in/reader033/viewer/2022042319/5f087bc07e708231d4223c27/html5/thumbnails/3.jpg)
How does a shock wave in a material look?
Typical laser-driven shocks exhibit pressures up to a few Mbar (few million atm) and temperatures up to several eV (few 10,000 K)."
possible ablation
drive pulse
1 mm
Jerry Forbes, Shock-Wave Compression of Condensed Matter: A Primer, Springer, 2012. "
HOW DOES A SHOCK WAVE IN A MATERIAL LOOK?
![Page 4: Atomistic Calculations of Dynamic Compression of Materialso Semiclassical quantum nuclear effects in classical molecular dynamics: An initial approach . How does a shock wave in a](https://reader033.vdocument.in/reader033/viewer/2022042319/5f087bc07e708231d4223c27/html5/thumbnails/4.jpg)
Shock compression of materials enables study of p,T states higher than statically achievable"
Proposed SiO2 phase diagram"
(100 Gpa = 1 million atmospheres)"
SHOCK COMPRESSION ENABLES STUDY OF HIGHER P, T STATES THAN STATIC METHODS
![Page 5: Atomistic Calculations of Dynamic Compression of Materialso Semiclassical quantum nuclear effects in classical molecular dynamics: An initial approach . How does a shock wave in a](https://reader033.vdocument.in/reader033/viewer/2022042319/5f087bc07e708231d4223c27/html5/thumbnails/5.jpg)
Shock waves are routinely generated in a controlled lab setting using projectiles from guns, lasers (including table-top scale), explosives, and
other means. !
Pressures >> 10 Mbar achievable, much greater than ~1-3 Mbar static pressures in diamond anvil cells. (Earth core is ~3.8 Mbar). !
Shock compression of materials enables study of p,T states higher than statically achievable"
SHOCK COMPRESSION ENABLES STUDY OF HIGHER P, T STATES THAN STATIC METHODS
![Page 6: Atomistic Calculations of Dynamic Compression of Materialso Semiclassical quantum nuclear effects in classical molecular dynamics: An initial approach . How does a shock wave in a](https://reader033.vdocument.in/reader033/viewer/2022042319/5f087bc07e708231d4223c27/html5/thumbnails/6.jpg)
1. The shock propagation speed (vs) exceeds the sound speed in the pre-shock material (c0)"
2. The shock propagation speed is less than the sound propagation speed behind the shock (c1+u1)"
Mechanical shock wave stability criteria"
Shock compression of materials enables study of p,T states higher than statically achievable"
CONDITIONS FOR SHOCK STABILITY
![Page 7: Atomistic Calculations of Dynamic Compression of Materialso Semiclassical quantum nuclear effects in classical molecular dynamics: An initial approach . How does a shock wave in a](https://reader033.vdocument.in/reader033/viewer/2022042319/5f087bc07e708231d4223c27/html5/thumbnails/7.jpg)
• Conservation of mass, momentum, and energy for a material give the 1D Euler equations."
• Find steady state solutions by considering the variables to be functions of "
• Integration and some manipulation yields local equations for thermodynamic variables in a shock wave"
particle velocity"
uniaxial stress"
energy per unit mass"
(Rayleigh line)"
€
u − u0 = vs − u0( ) 1− ρ0ρ
$
% &
'
( )
€
˜ e − ˜ e 0 = p01ρ0
−1ρ
$
% &
'
( ) +
u0 − vs( )2
21− ρ0
ρ
$
% &
'
( )
2
€
p − p0 = u0 − vs( )2ρ0 1−ρ0ρ
$
% &
'
( )
x − vst €
dρdt
+ ρ∂u∂x
= 0
€
dudt
+ ˜ v ∂p∂x
= 0
€
d ˜ e dt
+ p d ˜ v dt
= 0
Shock compression of materials enables study of p,T states higher than statically achievable"
A DESCRIPTION OF THE THERMODYNAMIC STATES BEHIND THE SHOCK FRONT
![Page 8: Atomistic Calculations of Dynamic Compression of Materialso Semiclassical quantum nuclear effects in classical molecular dynamics: An initial approach . How does a shock wave in a](https://reader033.vdocument.in/reader033/viewer/2022042319/5f087bc07e708231d4223c27/html5/thumbnails/8.jpg)
uniaxial stress" (Rayleigh line)"
€
p − p0 = u0 − vs( )2ρ0 1−ρ0ρ
$
% &
'
( )
x − vst €
dρdt
+ ρ∂u∂x
= 0
€
dudt
+ ˜ v ∂p∂x
= 0
€
d ˜ e dt
+ p d ˜ v dt
= 0
Shock compression of materials enables study of p,T states higher than statically achievable"
SHOCK STABILITY: THE GEOMETRICAL INTERPRETATION
Jerry Forbes, Shock-Wave Compression of Condensed Matter: A Primer, Springer, 2012. "
![Page 9: Atomistic Calculations of Dynamic Compression of Materialso Semiclassical quantum nuclear effects in classical molecular dynamics: An initial approach . How does a shock wave in a](https://reader033.vdocument.in/reader033/viewer/2022042319/5f087bc07e708231d4223c27/html5/thumbnails/9.jpg)
Shock compression is not isentropic"
Entropy change during (small) shock"compression:"
€
u − u0 = vs − u0( ) 1− ρ0ρ
$
% &
'
( )
x − vst
Isentropic compression can occur in the limit of an infinite number of small amplitude shocks.!
Shock compression of materials enables study of p,T states higher than statically achievable"
SHOCK COMPRESSION IS NOT ISENTROPIC
Jerry Forbes, Shock-Wave Compression of Condensed Matter: A Primer, Springer, 2012. "
![Page 10: Atomistic Calculations of Dynamic Compression of Materialso Semiclassical quantum nuclear effects in classical molecular dynamics: An initial approach . How does a shock wave in a](https://reader033.vdocument.in/reader033/viewer/2022042319/5f087bc07e708231d4223c27/html5/thumbnails/10.jpg)
The high strain rates, temperatures and stress in shocked materials introduce extreme atomic scale distortions and electronic excitation: !
Must use an atomic description!
€
Atomic equation of motion!
Shock compression of materials enables study of p,T states higher than statically achievable"
THE EFFECTS OF SHOCK COMPRESSION CAN BE STUDIED USING MOLECULAR DYNAMICS
![Page 11: Atomistic Calculations of Dynamic Compression of Materialso Semiclassical quantum nuclear effects in classical molecular dynamics: An initial approach . How does a shock wave in a](https://reader033.vdocument.in/reader033/viewer/2022042319/5f087bc07e708231d4223c27/html5/thumbnails/11.jpg)
A variety of atom-based energy models exist
Slower More accurate, more transferrable
Fewer assumptions
Faster Less accurate, less transferrable
More assumptions
Quantum chemical methods (CI, MP2,
MCSCF, etc.)
Density functional theory
Tight-binding Analytical models
101 102 103 104 106 109
Quantum energy models like DFT are (hopefully) generally more
accurate than empirical models, but at the price of increased computational expense
![Page 12: Atomistic Calculations of Dynamic Compression of Materialso Semiclassical quantum nuclear effects in classical molecular dynamics: An initial approach . How does a shock wave in a](https://reader033.vdocument.in/reader033/viewer/2022042319/5f087bc07e708231d4223c27/html5/thumbnails/12.jpg)
For analytical potentials: • Typical simulations:
• Thousands to millions of atoms and 10-100 ps timescales • Big simulations (need parallel computer):
• Billions of atoms, 10 nanosecond timescales For DFT:
• Typical simulations: • Fewer than 100 atoms and 1-10 ps
• Big simulations (need parallel computer): • 100-1000 atoms and 100 ps
CURRENT PRACTICAL LIMITATIONS OF MD
![Page 13: Atomistic Calculations of Dynamic Compression of Materialso Semiclassical quantum nuclear effects in classical molecular dynamics: An initial approach . How does a shock wave in a](https://reader033.vdocument.in/reader033/viewer/2022042319/5f087bc07e708231d4223c27/html5/thumbnails/13.jpg)
Multi-scale method: "Approximate hydrodynamics, few atoms,
more expensive (accurate) potentials including quantum methods."
Computational work ~t"
Direct shock MD method:"Exact hydrodynamics, but requires
cheap (poor) atomic potentials."Computational work ~t2 at best!
Thermodynamic constraints of a steady shock in a continuum are applied to a MD simulation as if the shock were passing over it."
Reed, Fried, Joannopoulos, Phys. Rev. Lett. 90, 235503 (2003)."
A MULTI-SCALE METHOD FOR SIMULATION OF MICROSCOPIC DYNAMICS IN SHOCK WAVES
![Page 14: Atomistic Calculations of Dynamic Compression of Materialso Semiclassical quantum nuclear effects in classical molecular dynamics: An initial approach . How does a shock wave in a](https://reader033.vdocument.in/reader033/viewer/2022042319/5f087bc07e708231d4223c27/html5/thumbnails/14.jpg)
"Conservation of mass, momentum, and energy for a continuum material give the 1D Navier-Stokes equations."
uniaxial stress"
energy per unit mass"
(Rayleigh line)"
€
˜ e − ˜ e 0 = p01ρ0
−1ρ
$
% &
'
( ) +
u0 − vs( )2
21− ρ0
ρ
$
% &
'
( )
2
€
p −µ˜ ˙ v ˜ v − p0 = u0 − vs( )2
ρ0 1− ρ0
ρ
$
% &
'
( )
€
dρdt
+ ρ∂u∂x
= 0
€
dudt
+ ˜ v ∂∂x
p −µ∂u∂x
$
% &
'
( ) = 0
€
d ˜ e dt
+ p −µ∂u∂x
$
% &
'
( )
d ˜ v dt
= 0
These relations hold everywhere in the wave (i.e., they are more than jump conditions).!
"Steady solutions are:"
Equations of motion for the atoms and volume of a computational cell are chosen to satisfy the steady solution conditions.!
There are 1 or 2 empirical parameters (computational cell mass and NS viscosity).!Reed, Fried, Joannopoulos, “Tractable dynamical studies of single and double shock compression,”
Phys. Rev. Lett. 90, 235503 (2003)."
A DESCRIPTION OF THERMODYNAMIC STATES BEHIND THE SHOCK FRONT
![Page 15: Atomistic Calculations of Dynamic Compression of Materialso Semiclassical quantum nuclear effects in classical molecular dynamics: An initial approach . How does a shock wave in a](https://reader033.vdocument.in/reader033/viewer/2022042319/5f087bc07e708231d4223c27/html5/thumbnails/15.jpg)
Molecular dynamics MSST equations of motion"
Computational cell volume:"
Atom positions (scaled):"
thermal stress"
cold stress" Hugoniot external potential"
damping (empirical viscosity)"
interatomic forces"
fictitious force due to
computational cell motion"
energy from volume damping"
shock speed is specified"
MSST EQUATIONS OF MOTION
![Page 16: Atomistic Calculations of Dynamic Compression of Materialso Semiclassical quantum nuclear effects in classical molecular dynamics: An initial approach . How does a shock wave in a](https://reader033.vdocument.in/reader033/viewer/2022042319/5f087bc07e708231d4223c27/html5/thumbnails/16.jpg)
"Conservation of mass, momentum, and energy for a continuum material give the 1D Navier-Stokes equations."
uniaxial stress"
energy per unit mass"
(Rayleigh line)"
€
˜ e − ˜ e 0 = p01ρ0
−1ρ
$
% &
'
( ) +
u0 − vs( )2
21− ρ0
ρ
$
% &
'
( )
2
€
p −µ˜ ˙ v ˜ v − p0 = u0 − vs( )2
ρ0 1− ρ0
ρ
$
% &
'
( )
€
dρdt
+ ρ∂u∂x
= 0
€
dudt
+ ˜ v ∂∂x
p −µ∂u∂x
$
% &
'
( ) = 0
€
d ˜ e dt
+ p −µ∂u∂x
$
% &
'
( )
d ˜ v dt
= 0
These relations hold everywhere in the wave (i.e., they are more than jump conditions).!
"Steady solutions are:"
Equations of motion for the atoms and volume of a computational cell are chosen to satisfy the steady solution conditions.!
There are 1 or 2 empirical parameters (computational cell mass and NS viscosity).!Reed, Fried, Joannopoulos, “Tractable dynamical studies of single and double shock compression,”
Phys. Rev. Lett. 90, 235503 (2003)."
A DESCRIPTION OF THERMODYNAMIC STATES BEHIND THE SHOCK FRONT
![Page 17: Atomistic Calculations of Dynamic Compression of Materialso Semiclassical quantum nuclear effects in classical molecular dynamics: An initial approach . How does a shock wave in a](https://reader033.vdocument.in/reader033/viewer/2022042319/5f087bc07e708231d4223c27/html5/thumbnails/17.jpg)
Excellent agreement is demonstrated between
MSST and direct approach for shocks in
Lennard-Jones, amorphous C (Tersoff potential) and other
materials.!
Reed, Maiti, Fried, Phys. Rev. E 81, 016607 (2010)."
MSST AGREES WITH DIRECT METHOD IN AMORPHOUS LENNARD-JONES
![Page 18: Atomistic Calculations of Dynamic Compression of Materialso Semiclassical quantum nuclear effects in classical molecular dynamics: An initial approach . How does a shock wave in a](https://reader033.vdocument.in/reader033/viewer/2022042319/5f087bc07e708231d4223c27/html5/thumbnails/18.jpg)
Snapshots of shocked nitromethane"
5 ps:!Hot unreacted liquid.!
8 ps:!Hydrogen transfer occurs,
consistent with density functional theory
simulations at similar thermodynamic
conditions: M. R. Manaa, et al., J. Chem. Phys. 120,
10146 (2004).!!
96 ps:!A dynamic
composition exists including water and numerous transient metastable species. !
Reed, Manaa, Fried, Glaesemann, Joannopoulos,” Nature Physics 4, 72 (2008).!
FIRST MICROSCOPIC PICTURE OF DETONATION: NITROMETHANE
![Page 19: Atomistic Calculations of Dynamic Compression of Materialso Semiclassical quantum nuclear effects in classical molecular dynamics: An initial approach . How does a shock wave in a](https://reader033.vdocument.in/reader033/viewer/2022042319/5f087bc07e708231d4223c27/html5/thumbnails/19.jpg)
We’ve been treating them as point particles, e.g. E=E({rj}). Is this assumption valid? The de Broglie wavelength is a rough measure of the length scale over which a quantum particle can be localized. When the distance between particles is less than the de Broglie wavelength, quantum effects are important.
ARE THE NUCLEI PARTICLES OR WAVES
![Page 20: Atomistic Calculations of Dynamic Compression of Materialso Semiclassical quantum nuclear effects in classical molecular dynamics: An initial approach . How does a shock wave in a](https://reader033.vdocument.in/reader033/viewer/2022042319/5f087bc07e708231d4223c27/html5/thumbnails/20.jpg)
• Typical distances between atoms range from 1 to 3 Angstroms (H2 bond length is ~0.74 Angstroms) • When the quantum wavelength exceeds the atomic separation, quantum effects play a role. • Electrons can be treated as classical point particles at sufficiently high temperatures and low densities
atom mass temperature de Broglie wavelegnth (Angstroms) H 1 300 1.44 Quantum! He 4 300 0.72 Li 7 300 0.55 Si 28 300 0.27 U 238 300 0.09
Liquid states: H 1 20 5.59 Quantum! He 4 4 6.25 Quantum!
Electron 0.000535294 300 62.37 Quantum! Electron 0.000535294 3000 19.72 Quantum! Electron 0.000535294 100000 3.42 Quantum!
SOME DE BROGLIE WAVELENGTH ESTIMATES
![Page 21: Atomistic Calculations of Dynamic Compression of Materialso Semiclassical quantum nuclear effects in classical molecular dynamics: An initial approach . How does a shock wave in a](https://reader033.vdocument.in/reader033/viewer/2022042319/5f087bc07e708231d4223c27/html5/thumbnails/21.jpg)
• Classical particles obey Boltzmann statistics while quantum particles obey Fermi or Bose statistics
• Quantum particles have zero-point energy: fluctuations even at T=0
• Classical particles evolve according to Newton’s equation while quantum particles obey a wave equation
• Quantum particles can tunnel
CLASSICAL VERSUS QUANTUM PARTICLES: SOME CONSEQUENCES
![Page 22: Atomistic Calculations of Dynamic Compression of Materialso Semiclassical quantum nuclear effects in classical molecular dynamics: An initial approach . How does a shock wave in a](https://reader033.vdocument.in/reader033/viewer/2022042319/5f087bc07e708231d4223c27/html5/thumbnails/22.jpg)
Post classical MD corrections of shock temperatures for quantum nuclear effects
show temperature shifts of hundreds of K for organic
molecular liquids."
THERE IS A NEED FOR INCORPORATION OF QUANTUM NUCLEAR EFFECTS IN SHOCK
CALCULATIONS
Methane
Water
N. Goldman, E. Reed, et. al. J. Chem. Phys., 131: 204103, 2009
![Page 23: Atomistic Calculations of Dynamic Compression of Materialso Semiclassical quantum nuclear effects in classical molecular dynamics: An initial approach . How does a shock wave in a](https://reader033.vdocument.in/reader033/viewer/2022042319/5f087bc07e708231d4223c27/html5/thumbnails/23.jpg)
QUANTUM NUCLEAR EFFECTS SIGNIFICANTLY IMPACT HEAT CAPACITY
• Heat capacity may differ from classical values by over 50%"• This alters Hugoniot temperatures from classical values"
"
0 500 1000 1500 2000 2500Temperature (K)
40
60
80
100
120
Hea
t Cap
acity
(J m
ol-1
K-1
)Classical
QM
Expt (Gurvich et al.)
"Heat capacity of liquid methane"
![Page 24: Atomistic Calculations of Dynamic Compression of Materialso Semiclassical quantum nuclear effects in classical molecular dynamics: An initial approach . How does a shock wave in a](https://reader033.vdocument.in/reader033/viewer/2022042319/5f087bc07e708231d4223c27/html5/thumbnails/24.jpg)
QUANTUM NUCLEAR EFFECTS CAN IMPACT MOL. AND PHASE STABILITY
• Bose-Einstein vibration energy spectrum can shift molecular formation energies by order tens of percent "
• Zero point energy for HN3 formation is 10 kcal/mol"• Reactions with many bonds breaking or forming can have large quantum
nuclear energy components"
"
0 50 100 1500
0.1
0.2
0.3E av
g (eV
) for
har
mon
ic o
scill
ator
0 50 100 1500
0.1
0.2
0.3
Vibrational frequency (THz)
H-C-H bend
C-C stretch
QM QM
Classical
Classical
T1 = 200 K T2 = 1000 K
C-H stretch
H-H stretch
C-C stretch
H-C-H bend
C-H stretch
H-H stretch
![Page 25: Atomistic Calculations of Dynamic Compression of Materialso Semiclassical quantum nuclear effects in classical molecular dynamics: An initial approach . How does a shock wave in a](https://reader033.vdocument.in/reader033/viewer/2022042319/5f087bc07e708231d4223c27/html5/thumbnails/25.jpg)
THE QUANTUM NUCLEAR EFFECT WISH LIST
• We would like to incorporate quanutum nuclear effects into shock simulations to describe"• Quantum heat capacity (to describe shock temperatures)"• Bose-Einstein spectrum (to describe chemistry and the
chemical equilibrium)"• Wave-particle effects (e.g. tunneling for transport and kinetics)"
• Wish list:"• Calculations will be self-consistent and on the fly"
• Not a post-classical MD temperature correction method"• Same computational cost (or nearly so) as classical MD"
• Path integral approaches are expensive for lower temperatures"
• Easy to use (like classical MD) and publically available "
"A new method that accomplishes most of these goals:"
QBMSST: Quantum Bath Multi-Scale Shock Technique!
![Page 26: Atomistic Calculations of Dynamic Compression of Materialso Semiclassical quantum nuclear effects in classical molecular dynamics: An initial approach . How does a shock wave in a](https://reader033.vdocument.in/reader033/viewer/2022042319/5f087bc07e708231d4223c27/html5/thumbnails/26.jpg)
In equation of motion, introduce both a random force R(t) and a dissipative force The random force spectrum is:
Θ~(ω,T ) = kBT
€
Θ~(ω,T ) =
12ω + ω[exp(ω / kBT )−1]
−1.
< R(t)R(t +τ )>= 2mγ Θ(ω,T )−∞
+∞
∫ exp[−iωτ ]dω2π.
H. Dammak, et. al. Phys. Rev. Lett. 103: 190601, (2009). M. Ceriotti, et. al. Phys. Rev. Lett. 103, 030603 (2009).
For quantum, Bose-Einstein ensemble averages
For classical (canonical) ensemble averages
Quantum thermal bath:!!
Langevin-type thermostats have been proposed to couple atoms to a Bose-Einstein thermal bath:!
mr(t) = f (r(t))+ R(t)−mγ r(t),
LANGEVIN THERMOSTAT APPROACH TO SEMICLASSICAL QUANTUM THERMAL BATH
H. Dammak, et. al. Phys. Rev. Lett. 103: 190601, (2009).
![Page 27: Atomistic Calculations of Dynamic Compression of Materialso Semiclassical quantum nuclear effects in classical molecular dynamics: An initial approach . How does a shock wave in a](https://reader033.vdocument.in/reader033/viewer/2022042319/5f087bc07e708231d4223c27/html5/thumbnails/27.jpg)
0 5 10 15 20Frequency (THz)
Phon
on d
ensi
ty o
f sta
tes (
arb.
uni
t) 10K QB10K classical
Liquid methane vibrational velcotiy
autocorrelation spectrum computed using classical MD
(red) and the quantum bath of Dammak et al.
(red)."
VIBRATION SPECTRUM OF LIQUID METHANE
0 50 100 1500
0.1
0.2
0.3
E avg (e
V) f
or h
arm
onic
osc
illat
or
0 50 100 1500
0.1
0.2
0.3
Vibrational frequency (THz)
H-C-H bend
C-C stretch
QM QM
Classical
Classical
T1 = 200 K T2 = 1000 K
C-H stretch
H-H stretch
C-C stretch
H-C-H bend
C-H stretch
H-H stretch
![Page 28: Atomistic Calculations of Dynamic Compression of Materialso Semiclassical quantum nuclear effects in classical molecular dynamics: An initial approach . How does a shock wave in a](https://reader033.vdocument.in/reader033/viewer/2022042319/5f087bc07e708231d4223c27/html5/thumbnails/28.jpg)
Heat capacity is corrected within the quantum thermal bath!
Gurvich et al. “Thermodynamic Properties of Individual Substances”, 4th ed. 1989, Vols. 1 and 2
0 500 1000 1500 2000 2500Temperature (K)
40
60
80
100
120
Hea
t Cap
acity
(J m
ol-1
K-1
)
Classical
QM
Expt (Gurvich et al.)
Zero-Point Energy 0 300 600 900 1200
Temperature (K)0
300
600
900
1200
Cla
ssic
al te
mpe
ratu
re 2
KE/
(3N
k B) (K
)QBClassical
QB (l.c. compressed 10%)
METHANE EXPERIMENTAL HEAT CAPACITY IS REPRODUCED BY QUANTUM BATH
Liquid methane heat capacity"
EAM Aluminum exhibits pressure-dependent zero-point
energy"
Langevin quantum bath describes the temperature and
pressure dependence of quantum heat
capacity.!
![Page 29: Atomistic Calculations of Dynamic Compression of Materialso Semiclassical quantum nuclear effects in classical molecular dynamics: An initial approach . How does a shock wave in a](https://reader033.vdocument.in/reader033/viewer/2022042319/5f087bc07e708231d4223c27/html5/thumbnails/29.jpg)
Volume EOM"
(QB)MSST EQUATIONS OF MOTION
Atom EOM"
Constraint quantity"
Quantum temperature
EOM"
γ and η are new parameters to be
chosen.!
T. Qi, E. J. Reed, J. Phys. Chem. A 116, 10451 (2012).
![Page 30: Atomistic Calculations of Dynamic Compression of Materialso Semiclassical quantum nuclear effects in classical molecular dynamics: An initial approach . How does a shock wave in a](https://reader033.vdocument.in/reader033/viewer/2022042319/5f087bc07e708231d4223c27/html5/thumbnails/30.jpg)
We combine quantum thermal bath with MSST (QB-MSST) for Methane!EFFICIENT GENERATION OF COLORED NOISE
White noise
Nf is the number of discrete points used to represent the Bose-
Einstein spectrum
QBMSST incorporates quantum nuclear effects with performance within factor of 2 of classical MSST.
Spectral filter
Random force
J. Barrat, et. al. J. Stat. Phys. 144: 679, 2011
![Page 31: Atomistic Calculations of Dynamic Compression of Materialso Semiclassical quantum nuclear effects in classical molecular dynamics: An initial approach . How does a shock wave in a](https://reader033.vdocument.in/reader033/viewer/2022042319/5f087bc07e708231d4223c27/html5/thumbnails/31.jpg)
We combine quantum thermal bath with MSST (QB-MSST) for Methane!
0 5 10 15 20 25 30Time (ps)
0
300
600
900
1200Te
mpe
ratu
re (K
)
QB-MSST
MSST
ΔT, QB-MSST
ΔT, MSST
~
~
6 KM/S SHOCK IN LIQUID METHANE
ReaxFF potential of Mattson et al, 216-molecules, 6 km/s shock speed, initial temperature: 111K, initial density: 0.432 g/cc.
Deviation from the constraint energy (ΔT) fluctuatates around zero for QBMSST.
![Page 32: Atomistic Calculations of Dynamic Compression of Materialso Semiclassical quantum nuclear effects in classical molecular dynamics: An initial approach . How does a shock wave in a](https://reader033.vdocument.in/reader033/viewer/2022042319/5f087bc07e708231d4223c27/html5/thumbnails/32.jpg)
Experimental data from: W. Nellis, et al., J. Chem. Phys. 75: 3055, 1981; H. Radousky, et al., J. Chem. Phys. 93: 8235, 1990
0.5 0.6 0.7 0.8 0.9 1 1.1 1.2Density (g/cc)
0
1000
2000
3000
4000
5000
Tem
pera
ture
(K)
QB-MSST (ReaxFF)MSST (ReaxFF)Radousky (expt)Goldman (DFT)Goldman (DFT, uncorrected)Li (DFT)
0.5 0.6 0.7 0.8 0.9 1 1.1 1.2Density (g/cc)
10
20
30
40
Pres
sure
(GPa
)
QB-MSST (ReaxFF)MSST (ReaxFF)Nellis (expt)Goldman (DFT, uncorrected)Li (DFT)
QBMSST LIQUID METHANE HUGONIOT COMPARISON
QBMSST Hugoniot differs from MSST by up to 5 GPa and 1000K at a given densitiy
QBMSST and MSST simulation durations of 300 ps
T. Qi, E. J. Reed, J. Phys. Chem. A 116, 10451 (2012).
![Page 33: Atomistic Calculations of Dynamic Compression of Materialso Semiclassical quantum nuclear effects in classical molecular dynamics: An initial approach . How does a shock wave in a](https://reader033.vdocument.in/reader033/viewer/2022042319/5f087bc07e708231d4223c27/html5/thumbnails/33.jpg)
QB-MSST enables self-consistent predictions for reaction phase boundaries!
100150200
08
16
0
50
100
0
30
60
2000 2500 3000 3500Temperature (K)
0
8
16
Num
ber o
f mol
ecul
es p
er c
ompu
tatio
nal c
ell CH4
CH3
H
H2
C2H6
MSSTQB-MSST
0
10
20
0
10
20
0
10
20
0
10
20
0 20 40 60 80 1000
10
20
0
10
20
0
10
20
0 20 40 60 80 1000
10
20
Num
bers
of h
ydro
carb
on p
er c
ompu
tatio
nal c
ell
9 km/s
10 km/s
11 km/s
12.2 km/s
Molar mass (g mol-1)
9 km/s
10 km/s
MSST
11 km/s
12.2 km/s
QB-MSST
CH3
C2H6
CH3
C2H6
T = 1560 K
T = 2050 K
T = 2630 K
T = 3300 K
T = 2120 K
T = 2540 K
T = 3020 K
T = 3710 K
p = 20.2 GPa
p = 25.9 GPa
p = 32.2 GPa
p = 40.5 GPa
p = 19.7 GPa
p = 25.3 GPa
p = 31.3 GPa
p = 39.6 GPa
SELF-CONSISTENT NATURE OF QBMSST CAN SHIFT CHEMICAL EQUILIBRIUM
The onset of methane dissociation occurs 15 GPa (~35%) lower pressure and 800 K lower T for QBMSST.
QBMSST and MSST simulation durations of 300 ps
T. Qi, E. J. Reed, J. Phys. Chem. A 116, 10451 (2012).
![Page 34: Atomistic Calculations of Dynamic Compression of Materialso Semiclassical quantum nuclear effects in classical molecular dynamics: An initial approach . How does a shock wave in a](https://reader033.vdocument.in/reader033/viewer/2022042319/5f087bc07e708231d4223c27/html5/thumbnails/34.jpg)
• Quantum nuclear effects in liquid methane:"• Shift Hugoniot temperature by over 30% (over 1000 K)"• Shift onset of chemistry to ~35% lower pressure"
• We have developed a semi-classical approach to self-consistent quantum nuclear effects in classical MD shock calculations"• Runs only ~20% slower than classical MD"
SUMMARY
T. Qi, E. J. Reed, J. Phys. Chem. A 116, 10451 (2012).
![Page 35: Atomistic Calculations of Dynamic Compression of Materialso Semiclassical quantum nuclear effects in classical molecular dynamics: An initial approach . How does a shock wave in a](https://reader033.vdocument.in/reader033/viewer/2022042319/5f087bc07e708231d4223c27/html5/thumbnails/35.jpg)
Acknowledgements
• Tingting Qi (now at LLNL)
Computational resources for this work:
• National Energy Research Science Computing Center (NERSC)
• National Nanotechnology Infrastructure Network (NNIN)
Conclusions ACKNOWLEDGEMENTS