Optimized nested Markov chain Monte Carlo sampling: application to the liquid nitrogen Hugoniot using density functional theory

Journal Article · Thu Jan 01 04:00:00 EST 2009

DOI:https://doi.org/10.1063/1.3295191· OSTI ID:991274

Shaw, Milton Sam ^[1]; Coe, Joshua D ^[1]; Sewell, Thomas D ^[2]

---

1. Los Alamos National Laboratory
2. UNIV OF MISSOURI-COLUMBIA

An optimized version of the Nested Markov Chain Monte Carlo sampling method is applied to the calculation of the Hugoniot for liquid nitrogen. The 'full' system of interest is calculated using density functional theory (DFT) with a 6-31 G* basis set for the configurational energies. The 'reference' system is given by a model potential fit to the anisotropic pair interaction of two nitrogen molecules from DFT calculations. The EOS is sampled in the isobaric-isothermal (NPT) ensemble with a trial move constructed from many Monte Carlo steps in the reference system. The trial move is then accepted with a probability chosen to give the full system distribution. The P's and T's of the reference and full systems are chosen separately to optimize the computational time required to produce the full system EOS. The method is numerically very efficient and predicts a Hugoniot in excellent agreement with experimental data.

• Research Organization: Los Alamos National Laboratory (LANL)
• Sponsoring Organization: DOE
• DOE Contract Number: AC52-06NA25396
• OSTI ID: 991274
• Report Number(s): LA-UR-09-04904; LA-UR-09-4904
• Country of Publication: United States
• Language: English

Related Subjects

71 CLASSICAL AND QUANTUM MECHANICS
GENERAL PHYSICS
CHAINS
DISTRIBUTION
FUNCTIONALS
NITROGEN
PROBABILITY
SAMPLING