improved hard sphere radial distribution function in the cris ......benjamin j. cowen & john h....
TRANSCRIPT
A u t h o r s :
Sandia National Laboratories is a multimission
laboratory managed and operated by National
Technology & Engineering Solutions of Sandia,
LLC, a wholly owned subsidiary of Honeywell
International Inc., for the U.S. Department of
Energyโs National Nuclear Security
Administration under contract DE-NA0003525.
Improved hard sphere radial distribution function in the CRIS equation of state model
Ben jamin J. Cowen & John H. Car pen ter
1
Background
โชAccurate equation of state (EOS) models for various materials are needed as input into high-level shock physics codes
โชThe EOS is typically defined in terms of the Helmholtz free energy, which is often split into 3 terms:
๐ด ๐, ๐ = ๐ด0 ๐ + ๐ด๐ ๐, ๐ + ๐ด๐ ๐, ๐
โชThe cold curve is used as an input into the CRIS model, which can be obtained from theory, experiment, or both
โชThe output of the CRIS model includes the first two terms of the Helmholtz free energy equation: ๐ด0 ๐ +๐ด๐ ๐, ๐ , which can then be added to the thermal electronic term for a full EOS
โชThe CRIS model has previously been used to develop the equation of state for metals (Au, Mo, Al, Ta, Pb, Ti, Cu, W), gases (H, D2, N, O, C, CO, CH4, Xe, Ar), and other materials (Basalt, Ice, CaCO3, SiO2)
2
Cold Curve
Thermal Ionic
Thermal Electronic
โขCalculates thermodynamic properties by expanding about a hard sphere reference system
โขFirst we start with the free energy of the hard sphere reference fluid ๐ด0
โขThen we add the first order term: ๐ 0 =โ2๐
3๐0๐๐๐ ๐, ๐ ๐ ๐, ๐ ๐โ1๐๐
โข ๐ ๐, ๐ = scaled cold curve, used as the interaction potential
โข ๐ ๐, ๐ = radial distribution function (RDF)
โข ๐ = packing fraction, ๐0 = hard sphere diameter
โข ๐๐ = cutoff radius determined by solving the normalization condition: 1 =โ2๐
30๐๐ ๐ ๐, ๐ ๐โ1๐๐
โขUp to first order, we have: าง๐ด = ๐ด0 + ๐ ๐ 0
โขWe can then calculate the hard sphere diameter by minimizing: ๐ าง๐ด
๐๐ ๐,๐,๐= 0
โขThe 2nd order corrections are defined as:โข Fluctuation correction:
โข ฮ๐ด1 = โ๐๐๐โ2๐
3๐๐0๐๐ ๐๐
๐๐๐โ เดฅ๐
๐๐ ๐, าง๐ ๐ าง๐ ๐ าง๐
โข ๐ าง๐ =๐(เดฅ๐๐ง)
๐เดฅ๐
โข Soft sphere correction:
โข ฮ๐ด2 = โ๐๐๐2๐๐ 0๐0 ๐โ
๐
2๐ ๐, ๐โ ๐2๐๐
โขThe Helmholtz free energy of the system is defined as: ๐ด = ๐ด0 +๐ ๐ 0 + ฮ๐ด1 + ฮ๐ด2
The CRIS Model3
KERLEY, G. I. (1980). PERTURBATION THEORY AND THE THERMODYNAMIC PROPERTIES OF FLUIDS. II. THE CRIS MODEL. THE JOURNAL OF CHEMICAL PHYSICS, 73(1), 478-486.
Challenges with the CRIS model4
๐ ๐ฅ, ๐ =๐ฅ๐0(1, ๐)
1 + ๐ถ ๐ ๐ฅ โ 1 3 ๐ถ ๐ =1 + ฯ๐=1
3 ๐ต๐๐๐๐
๐
3 1 โ๐๐๐
Z= ๐๐
๐๐๐= 4๐๐0 1, ๐ + 1 = 1 +
3๐
๐๐โ๐+ ฯ๐=1
4 ๐๐ด๐๐
๐๐
๐
ALDER, B. J., & WAINWRIGHT, T. E. (1960). STUDIES IN MOLECULAR DYNAMICS. II. BEHAVIOR OF A SMALL NUMBER OF ELASTIC SPHERES. THE JOURNAL OF CHEMICAL PHYSICS, 33(5), 1439-1451.
KOLAFA, J., LABรK, S., & MALIJEVSKร, A. (2004). ACCURATE EQUATION OF STATE OF THE HARD SPHERE FLUID IN STABLE AND METASTABLE REGIONS. PHYSICAL CHEMISTRY CHEMICAL PHYSICS, 6(9), 2335-2340.
โข ๐ ๐ฅ, ๐ does not represent the hard sphere RDF well over the range of integration
โข ฮA2 integrates x<1, and there is a pole in ๐ ๐ฅ, ๐ in this regime
โข The fit to Z is good for ๐ < 0.5, but is slightly inaccurate for ๐ > 0.5
Improving the RDF5
โข ๐ ๐ฅ, ๐ is fit to a 5th order polynomial for x โ [1:xM ]: ๐ ๐ฅ, ๐ = ๐0 1, ๐ ฯ๐,๐=0๐,๐=5
๐ถ๐๐ ๐ฅ โ 1 ๐ ๐
๐๐
๐
โข Coefficients are fixed due to asymptotic limits:
โข ๐ถ00 = 1 and ๐ถ๐0 = 0 for ๐ > 0 since ๐ ๐ฅ,๐
๐0(1,๐)= 1 when ๐ = 0 and ๐ถ0๐ = 0 for ๐ > 0 since
๐ ๐ฅ,๐
๐0(1,๐)= 1 when ๐ฅ = 1
โข ๐ถ11 = โ9
2๐๐, ๐ถ21 =
3
2๐๐, ๐ถ31 =
1
2๐๐, ๐ถ41 = ๐ถ51 = 0 from the density expansion as ๐ โ 0
โข 20 unconstrained parameters remain
โข Z is refit (first 3 coefficients are the same) so that it is consistent with the ๐(๐ฅ, ๐) approximation: ๐0(1, ๐) =1
4๐
3๐
๐๐โ๐+ ฯ๐=1
1โ4,10,14 ๐๐ด๐๐
๐๐
๐
Effects of Changing the RDF6
VERLET, L., & WEIS, J. J. (1972). EQUILIBRIUM THEORY OF SIMPLE LIQUIDS. PHYSICAL REVIEW A, 5(2), 939.
POTOFF, J. J., & PANAGIOTOPOULOS, A. Z. (2000). SURFACE TENSION OF THE THREE-DIMENSIONAL LENNARD-JONES FLUID FROM HISTOGRAM-REWEIGHTING MONTE CARLO SIMULATIONS. THE JOURNAL OF CHEMICAL PHYSICS, 112(14), 6411-6415.
CAILLOL, J. M. (1998). CRITICAL-POINT OF THE LENNARD-JONES FLUID: A FINITE-SIZE SCALING STUDY. THE JOURNAL OF CHEMICAL PHYSICS, 109(12), 4885-4893.
โข We used the Lennard-Jones (LJ) potential as a test case, since we know the exact potential
โข At lower densities, the relative effect of the RDF is not as important since the ideal gas term is
dominant
โข The new RDF improves the pressure, and is more accurate in the liquid regime than the old RDF
โข The new RDF makes the estimate of the critical point worse
Vapor DomeIsotherms
Summary & Conclusions
โขThe CRIS model has been used for decades at institutions around the world for developing equations of state for various materials
โขHowever, the CRIS model inaccurately models the RDF for the hard spheres
โขWe improved the RDF using a 5th order polynomial function and added
higher order terms to ๐๐
๐๐๐to make it consistent
โขThe improved RDF more accurately models the isotherms of the LJ fluid, compared to molecular dynamics
โขUnfortunately, the critical point of the LJ fluid is more over-estimated with the new RDF
โขSince the critical point is extremely sensitive to small changes in the free energy, future work may involve fitting the RDF to have the EOS better match the vapor dome, while maintaining accurate isotherms
7
Any Questions?
8
Backup Slides
CRIS ฮA2 Integration10
โขฮ๐ด2 = โ๐๐๐2๐๐ 0๐0 ๐โ
๐
2๐ ๐, ๐โ ๐2๐๐
Isotherms of the LJ Fluid11
๐ด = ๐ด0 + ๐ ๐ 0 ๐ด = ๐ด0 +๐ ๐ 0 + ฮ๐ด1 ๐ด = ๐ด0 + ๐ ๐ 0 + ฮ๐ด1 + ฮ๐ด2
VERLET, L., & WEIS, J. J. (1972). EQUILIBRIUM THEORY OF SIMPLE LIQUIDS. PHYSICAL REVIEW A, 5(2), 939.
Vapor Dome of the LJ Fluid12
๐ด = ๐ด0 + ๐ ๐ 0 ๐ด = ๐ด0 +๐ ๐ 0 + ฮ๐ด1 ๐ด = ๐ด0 + ๐ ๐ 0 + ฮ๐ด1 + ฮ๐ด2
POTOFF, J. J., & PANAGIOTOPOULOS, A. Z. (2000). SURFACE TENSION OF THE THREE-DIMENSIONAL LENNARD-JONES FLUID FROM HISTOGRAM-REWEIGHTING MONTE CARLO SIMULATIONS. THE JOURNAL OF CHEMICAL PHYSICS, 112(14), 6411-6415.
CAILLOL, J. M. (1998). CRITICAL-POINT OF THE LENNARD-JONES FLUID: A FINITE-SIZE SCALING STUDY. THE JOURNAL OF CHEMICAL PHYSICS, 109(12), 4885-4893.