Predicting Melt Curves of Energetic Materials Using Molecular Models

Tow, Garrett M.; Larentzos, James P.; Sellers, Michael S.; Lísal, Martin; Brennan, John K.

doi:10.1002/prep.202100363

Predicting Melt Curves of Energetic Materials Using Molecular Models

Journal Article · Mon Apr 04 00:00:00 EDT 2022 · Propellants, Explosives, Pyrotechnics

DOI:https://doi.org/10.1002/prep.202100363· OSTI ID:1981430

Tow, Garrett M. ^[1]; Larentzos, James P. ^[1]; Sellers, Michael S. ^[2]; Lísal, Martin ^[3]; ^[1]

US Dept. of Homeland Security (DHS), Aberdeen Proving Ground, MD (United States). Army Research Laboratory
Booz Allen Hamilton Inc., McLean, VA (United States)
Czech Academy of Sciences (CAS), Prague (Czech Republic); Jan Evangelista Purkyně University (Czech Republic)

Here, in this work, the solid–liquid coexistence curves of classical fully flexible atomistic models of α-RDX and β-HMX were calculated using thermodynamically rigorous methodologies that identify where the free energy difference between the phases is zero. The free energy difference between each phase at a given state point was computed using the pseudosupercritical path (PSCP) method, and Gibbs–Helmholtz integration was used to evaluate the solid–liquid free energy difference as a function of temperature. This procedure was repeated for several pressures to determine points along the coexistence curve, which were then fit to the Simon–Glatzel functional form. While effective, this method is computationally expensive. An alternative approach is to compute the melting point at a single pressure via the PSCP method, and then use the Gibbs–Duhem integration technique to trace out the coexistence curve in a more computationally economical manner. Both approaches were used to determine the coexistence curve of α-RDX. The Gibbs–Duhem integration method was shown to generate a melt curve that is in good agreement with the PSCP-derived melt curve, while only costing ~10 % of the computational resources used for the PSCP method. For α-RDX, the predicted melting temperature increases significantly more for a given increase in pressure when compared to available experimental data.

View Accepted Manuscript (DOE)

Research Organization:: National Energy Technology Laboratory (NETL), Pittsburgh, PA, Morgantown, WV (United States)

Sponsoring Organization:: USDOE; US Army Research Laboratory (USARL); US Army Research Office (ARO)

OSTI ID:: 1981430

Journal Information:: Propellants, Explosives, Pyrotechnics, Journal Name: Propellants, Explosives, Pyrotechnics Journal Issue: 8 Vol. 47; ISSN 0721-3115

Publisher:: WileyCopyright Statement

Country of Publication:: United States

Language:: English

References (55)

PACKMOL: A package for building initial configurations for molecular dynamics simulations Martínez, L.; Andrade, R.; Birgin, E. G. Journal of Computational Chemistry, Vol. 30, Issue 13 https://doi.org/10.1002/jcc.21224	journal	October 2009
Tandem Molecular Dynamics and Continuum Studies of Shock‐Induced Pore Collapse in TATB Zhao, Puhan; Lee, Sangyup; Sewell, Tommy Propellants, Explosives, Pyrotechnics, Vol. 45, Issue 2 https://doi.org/10.1002/prep.201900382	journal	November 2019
Bemerkungen zur Schmelzdruckkurve Simon, Franz; Glatzel, Gunther Zeitschrift f�r anorganische und allgemeine Chemie, Vol. 178, Issue 1 https://doi.org/10.1002/zaac.19291780123	journal	January 1929
A theoretical investigation into the cooperativity effect on the TNT melting point under external electric field Ren, Fu-de; Shi, Wen-jing; Cao, Duan-lin Journal of Molecular Modeling, Vol. 27, Issue 1 https://doi.org/10.1007/s00894-020-04624-w	journal	January 2021
Thermal Expansion of HMX Deschamps, Jeffrey R.; Frisch, Mark; Parrish, Damon Journal of Chemical Crystallography, Vol. 41, Issue 7 https://doi.org/10.1007/s10870-011-0026-6	journal	February 2011
Melting curve for argon calculated from pure theory Solca, Jan; Dyson, Anthony J.; Steinebrunner, Gerold Chemical Physics, Vol. 224, Issue 2-3 https://doi.org/10.1016/S0301-0104(97)00317-0	journal	December 1997
A generalized crystal-cutting method for modeling arbitrarily oriented crystals in 3D periodic simulation cells with applications to crystal–crystal interfaces Kroonblawd, Matthew P.; Mathew, Nithin; Jiang, Shan Computer Physics Communications, Vol. 207 https://doi.org/10.1016/j.cpc.2016.07.007	journal	October 2016
LAMMPS - a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales Thompson, Aidan P.; Aktulga, H. Metin; Berger, Richard Computer Physics Communications, Vol. 271 https://doi.org/10.1016/j.cpc.2021.108171	journal	February 2022
Sensitivity of pore collapse heating to the melting temperature and shear viscosity of HMX Kroonblawd, Matthew P.; Austin, Ryan A. Mechanics of Materials, Vol. 152 https://doi.org/10.1016/j.mechmat.2020.103644	journal	January 2021
Thermodynamics of HMX Polymorphs and HMX/RDX Mixtures Myint, Philip C.; Nichols, Albert L. Industrial & Engineering Chemistry Research, Vol. 56, Issue 1 https://doi.org/10.1021/acs.iecr.6b03697	journal	December 2016
Computing the Liquidus of Binary Monatomic Salt Mixtures with Direct Simulation and Alchemical Free Energy Methods DeFever, Ryan S.; Maginn, Edward J. The Journal of Physical Chemistry A, Vol. 125, Issue 38 https://doi.org/10.1021/acs.jpca.1c06107	journal	September 2021
Melting Point Prediction of Energetic Materials via Continuous Heating Simulation on Solid-to-Liquid Phase Transition Liu, Yingzhe; Lai, Weipeng; Yu, Tao ACS Omega, Vol. 4, Issue 2 https://doi.org/10.1021/acsomega.8b03597	journal	February 2019
Molecular Dynamics Simulations of the Melting of 1,3,3-Trinitroazetidine Agrawal, Paras M.; Rice, Betsy M.; Zheng, Lianqing The Journal of Physical Chemistry B, Vol. 110, Issue 11 https://doi.org/10.1021/jp056690g	journal	March 2006
Molecular Dynamics Simulations of Hexahydro-1,3,5-trinitro-1,3,5-s-triazine (RDX) Using a Combined Sorescu−Rice−Thompson AMBER Force Field Agrawal, Paras M.; Rice, Betsy M.; Zheng, Lianqing The Journal of Physical Chemistry B, Vol. 110, Issue 51 https://doi.org/10.1021/jp065241t	journal	November 2006
Molecular Dynamics Simulations of the Melting Mechanisms of Perfect and Imperfect Crystals of Dimethylnitramine Zheng, Lianqing; Rice, Betsy M.; Thompson, Donald L. The Journal of Physical Chemistry B, Vol. 111, Issue 11 https://doi.org/10.1021/jp0667184	journal	March 2007
Melting Point Determination from Solid−Liquid Coexistence Initiated by Surface Melting Siavosh-Haghighi, Ali; Thompson, Donald L. The Journal of Physical Chemistry C, Vol. 111, Issue 22 https://doi.org/10.1021/jp070242m	journal	June 2007
Phase Diagram of Hexahydro-1,3,5-trinitro-1,3,5-triazine Crystals at High Pressures and Temperatures Dreger, Zbigniew A.; Gupta, Yogendra M. The Journal of Physical Chemistry A, Vol. 114, Issue 31 https://doi.org/10.1021/jp105226s	journal	August 2010
Simulations of High-Pressure Phases in RDX Munday, Lynn B.; Chung, Peter W.; Rice, Betsy M. The Journal of Physical Chemistry B, Vol. 115, Issue 15 https://doi.org/10.1021/jp112042a	journal	April 2011
Quantum Chemistry Based Force Field for Simulations of HMX Smith, Grant D.; Bharadwaj, Rishikesh K. The Journal of Physical Chemistry B, Vol. 103, Issue 18 https://doi.org/10.1021/jp984599p	journal	May 1999
Free-energy calculations using classical molecular simulation: application to the determination of the melting point and chemical potential of a flexible RDX model Sellers, Michael S.; Lísal, Martin; Brennan, John K. Physical Chemistry Chemical Physics, Vol. 18, Issue 11 https://doi.org/10.1039/C5CP06164D	journal	January 2016
Comparative study of melting points of 3,4-bis(3-nitrofurazan-4-yl)furoxan (DNTF)/1,3,3-trinitroazetidine (TNAZ) eutectic compositions using molecular dynamic simulations Liu, Ning; Zeman, Svatopluk; Shu, Yuan-jie RSC Advances, Vol. 6, Issue 64 https://doi.org/10.1039/C6RA12041E	journal	January 2016
Thermal decomposition and phase transitions in solid nitramines Hall, P. G. Transactions of the Faraday Society, Vol. 67 https://doi.org/10.1039/tf9716700556	journal	January 1971
Molecular dynamics studies of melting and liquid properties of ammonium dinitramide Velardez, Gustavo F.; Alavi, Saman; Thompson, Donald L. The Journal of Chemical Physics, Vol. 119, Issue 13 https://doi.org/10.1063/1.1605380	journal	October 2003
Molecular dynamics study of the melting of nitromethane Agrawal, Paras M.; Rice, Betsy M.; Thompson, Donald L. The Journal of Chemical Physics, Vol. 119, Issue 18 https://doi.org/10.1063/1.1612915	journal	November 2003
Molecular dynamics studies of melting and solid-state transitions of ammonium nitrate Velardez, Gustavo F.; Alavi, Saman; Thompson, Donald L. The Journal of Chemical Physics, Vol. 120, Issue 19 https://doi.org/10.1063/1.1705573	journal	May 2004
Toward a robust and general molecular simulation method for computing solid-liquid coexistence Eike, David M.; Brennecke, Joan F.; Maginn, Edward J. The Journal of Chemical Physics, Vol. 122, Issue 1 https://doi.org/10.1063/1.1823371	journal	January 2005
Molecular dynamics simulations of melting and the glass transition of nitromethane Zheng, Lianqing; Luo, Sheng-Nian; Thompson, Donald L. The Journal of Chemical Physics, Vol. 124, Issue 15 https://doi.org/10.1063/1.2174002	journal	April 2006
The melting point of ice Ih for common water models calculated from direct coexistence of the solid-liquid interface García Fernández, Ramón; Abascal, José L. F.; Vega, Carlos The Journal of Chemical Physics, Vol. 124, Issue 14 https://doi.org/10.1063/1.2183308	journal	April 2006
Atomistic simulation of solid-liquid coexistence for molecular systems: Application to triazole and benzene Eike, David M.; Maginn, Edward J. The Journal of Chemical Physics, Vol. 124, Issue 16 https://doi.org/10.1063/1.2188400	journal	April 2006
Molecular dynamics simulations of melting of perfect crystalline hexahydro-1,3,5-trinitro-1,3,5-s-triazine Zheng, Lianqing; Thompson, Donald L. The Journal of Chemical Physics, Vol. 125, Issue 8 https://doi.org/10.1063/1.2238860	journal	August 2006
Revisiting the Frenkel-Ladd method to compute the free energy of solids: The Einstein molecule approach Vega, Carlos; Noya, Eva G. The Journal of Chemical Physics, Vol. 127, Issue 15 https://doi.org/10.1063/1.2790426	journal	October 2007
Shock-induced melting of (100)-oriented nitromethane: Structural relaxation Siavosh-Haghighi, Ali; Dawes, Richard; Sewell, Thomas D. The Journal of Chemical Physics, Vol. 131, Issue 6 https://doi.org/10.1063/1.3202441	journal	January 2009
A molecular dynamics simulation study of crystalline 1,3,5-triamino-2,4,6-trinitrobenzene as a function of pressure and temperature Bedrov, Dmitry; Borodin, Oleg; Smith, Grant D. The Journal of Chemical Physics, Vol. 131, Issue 22 https://doi.org/10.1063/1.3264972	journal	December 2009
Shock-induced melting of (100)-oriented nitromethane: Energy partitioning and vibrational mode heating Dawes, Richard; Siavosh-Haghighi, Ali; Sewell, Thomas D. The Journal of Chemical Physics, Vol. 131, Issue 22 https://doi.org/10.1063/1.3271349	journal	December 2009
A comparison of methods for melting point calculation using molecular dynamics simulations Zhang, Yong; Maginn, Edward J. The Journal of Chemical Physics, Vol. 136, Issue 14 https://doi.org/10.1063/1.3702587	journal	April 2012
New Monte Carlo method to compute the free energy of arbitrary solids. Application to the fcc and hcp phases of hard spheres Frenkel, Daan; Ladd, Anthony J. C. The Journal of Chemical Physics, Vol. 81, Issue 7 https://doi.org/10.1063/1.448024	journal	October 1984
Error estimates on averages of correlated data Flyvbjerg, H.; Petersen, H. G. The Journal of Chemical Physics, Vol. 91, Issue 1 https://doi.org/10.1063/1.457480	journal	July 1989
Constant pressure molecular dynamics algorithms Martyna, Glenn J.; Tobias, Douglas J.; Klein, Michael L. The Journal of Chemical Physics, Vol. 101, Issue 5 https://doi.org/10.1063/1.467468	journal	September 1994
Melting curves for neon calculated from pure theory Solca, Jan; Dyson, Anthony J.; Steinebrunner, Gerold The Journal of Chemical Physics, Vol. 108, Issue 10 https://doi.org/10.1063/1.475808	journal	March 1998
On fluid-solid direct coexistence simulations: The pseudo-hard sphere model Espinosa, Jorge R.; Sanz, Eduardo; Valeriani, Chantal The Journal of Chemical Physics, Vol. 139, Issue 14 https://doi.org/10.1063/1.4823499	journal	October 2013
Anisotropy in surface-initiated melting of the triclinic molecular crystal 1,3,5-triamino-2,4,6-trinitrobenzene: A molecular dynamics study Mathew, N.; Sewell, Thomas D.; Thompson, Donald L. The Journal of Chemical Physics, Vol. 143, Issue 9 https://doi.org/10.1063/1.4929806	journal	September 2015
Heat Capacities and Entropies of Organic Compounds in the Condensed Phase. Volume III Domalski, Eugene S.; Hearing, Elizabeth D. Journal of Physical and Chemical Reference Data, Vol. 25, Issue 1 https://doi.org/10.1063/1.555985	journal	January 1996
Melting points of alkali chlorides evaluated for a polarizable and non-polarizable model DeFever, Ryan S.; Wang, Haimeng; Zhang, Yong The Journal of Chemical Physics, Vol. 153, Issue 1 https://doi.org/10.1063/5.0012253	journal	July 2020
Molecular dynamics-guided material model for the simulation of shock-induced pore collapse in β-octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (β-HMX) Das, Pratik; Zhao, Puhan; Perera, Dilki Journal of Applied Physics, Vol. 130, Issue 8 https://doi.org/10.1063/5.0056560	journal	August 2021
Exponential-six potential scaling for the calculation of free energies in molecular simulations Sellers, Michael S.; Lísal, Martin; Brennan, John K. Molecular Physics, Vol. 113, Issue 1 https://doi.org/10.1080/00268976.2014.942405	journal	August 2014
Gibbs-Duhem integration: a new method for direct evaluation of phase coexistence by molecular simulation Kofke, David A. Molecular Physics, Vol. 78, Issue 6 https://doi.org/10.1080/00268979300100881	journal	April 1993
Predicted melt curve and liquid-state transport properties of TATB from molecular dynamics simulations Mathew, Nithin; Kroonblawd, Matthew P.; Sewell, Tommy Molecular Simulation, Vol. 44, Issue 8 https://doi.org/10.1080/08927022.2017.1418084	journal	January 2018
Direct Evaluation of Vapour-Liquid Equilibria by Molecular Dynamics using Gibbs-Duhem Integration Lísal, Martin; Vacek, Václav Molecular Simulation, Vol. 17, Issue 1 https://doi.org/10.1080/08927029608024092	journal	March 1996
A Liouville-operator derived measure-preserving integrator for molecular dynamics simulations in the isothermal–isobaric ensemble Tuckerman, Mark E.; Alejandre, José; López-Rendón, Roberto Journal of Physics A: Mathematical and General, Vol. 39, Issue 19 https://doi.org/10.1088/0305-4470/39/19/S18	journal	April 2006
Molecular-dynamics study of lattice-defect-nucleated melting in silicon Phillpot, S. R.; Lutsko, J. F.; Wolf, D. Physical Review B, Vol. 40, Issue 5 https://doi.org/10.1103/PhysRevB.40.2831	journal	August 1989
Molecular-dynamics study of lattice-defect-nucleated melting in metals using an embedded-atom-method potential Lutsko, J. F.; Wolf, D.; Phillpot, S. R. Physical Review B, Vol. 40, Issue 5 https://doi.org/10.1103/PhysRevB.40.2841	journal	August 1989
Melting line of aluminum from simulations of coexisting phases Morris, J. R.; Wang, C. Z.; Ho, K. M. Physical Review B, Vol. 49, Issue 5 https://doi.org/10.1103/PhysRevB.49.3109	journal	February 1994
Rapid estimation of elastic constants by molecular dynamics simulation under constant stress Shinoda, Wataru; Shiga, Motoyuki; Mikami, Masuhiro Physical Review B, Vol. 69, Issue 13 https://doi.org/10.1103/PhysRevB.69.134103	journal	April 2004
The crystal structure of cyclotrimethylenetrinitramine Choi, C. S.; Prince, E. Acta Crystallographica Section B Structural Crystallography and Crystal Chemistry, Vol. 28, Issue 9 https://doi.org/10.1107/S0567740872007046	journal	September 1972
Computer Simulation Using Particles Hockney, R.; Eastwood, J. CRC Press https://doi.org/10.1201/9781439822050	book	January 1988

Similar Records

High-pressure phase diagram of beryllium from ab initio free-energy calculations

Journal Article · Fri Jul 09 00:00:00 EDT 2021 · Physical Review. B · OSTI ID:1850371

A self-contained algorithm for determination of solid-liquid equilibria in an alloy system

Journal Article · Fri Apr 20 00:00:00 EDT 2018 · Computational Materials Science · OSTI ID:1459538

Predicted Melt Curve and Liquid Shear Viscosity of RDX up to 30 GPa

Journal Article · Thu Mar 10 23:00:00 EST 2022 · Propellants, Explosives, Pyrotechnics · OSTI ID:1994044

Related Subjects

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEFENSE
HMX
RDX
melting point
molecular dynamics
phase coexistence

Predicting Melt Curves of Energetic Materials Using Molecular Models

Citation Formats

References (55)

Similar Records

Related Subjects