اشترك بالحزمة الذهبية واحصل على وصول غير محدود شمرا أكاديميا
تسجيل مستخدم جديدAtom-in-jellium equations of state for cryogenic liquids
336
0
0.0
(
0
)
اسأل ChatGPT حول البحث
ﻻ يوجد ملخص باللغة العربية
Equations of state (EOS) calculated from a computationally efficient atom-in-jellium treatment of the electronic structure have recently been shown to be consistent with more rigorous path integral Monte Carlo (PIMC) and quantum molecular dynamics (QMD) simulations of metals in the warm dense matter regime. Here we apply the atom-in-jellium model to predict wide-ranging EOS for the cryogenic liquid elements nitrogen, oxygen, and fluorine. The principal Hugoniots for these substances were surprisingly consistent with available shock data and Thomas-Fermi (TF) EOS for very high pressures, and exhibited systematic variations from TF associated with shell ionization effects, in good agreement with PIMC, though deviating from QMD and experiment in the molecular regime. The new EOS are accurate much higher in pressure than previous widely-used models for nitrogen and oxygen in particular, and should allow much more accurate predictions for oxides and nitrides in the liquid, vapor, and plasma regime, where these have previously been constructed as mixtures containing the older EOS.
قيم البحث
اقرأ أيضاً
Recent path-integral Monte Carlo and quantum molecular dynamics simulations have shown that computationally efficient average-atom models can predict thermodynamic states in warm dense matter to within a few percent. One such atom-in-jellium model ha
s typically been used to predict the electron-thermal behavior only, although it was previously developed to predict the entire equation of state (EOS). We report completely atom-in-jellium EOS calculations for Be, Al, Si, Fe, and Mo, as elements representative of a range of atomic number and low-pressure electronic structure. Comparing the more recent method of pseudo-atom molecular dynamics, atom-in-jellium results were similar: sometimes less accurate, sometimes more. All these techniques exhibited pronounced effects of electronic shell structure in the shock Hugoniot which are not captured by Thomas-Fermi based EOS. These results demonstrate the value of a hierarchical approach to EOS construction, using average-atom techniques with shell structure to populate a wide-range EOS surface efficiently, complemented by more rigorous 3D multi-atom calculations to validate and adjust the EOS.
Ru and Rh are interesting cases for comparing equations of state (EOS), because most general purpose EOS are semi-empirical, relying heavily on shock data, and none has been reported for Ru. EOS were calculated for both elements using all-electron at
om-in-jellium theory, and cold compression curves were calculated for the common crystal types using the multi-ion pseudopotential approach. Previous EOS constructed for these elements used Thomas-Fermi (TF) theory for the electronic behavior at high temperatures, which neglects electronic shell structure; the atom-in-jellium EOS exhibited pronounced features from the excitation of successive electron shells. Otherwise, the EOS matched surprisingly well, especially considering the lack of experimental data for Ru. The TF-based EOS for Ru may however be inaccurate in the multi-terapascal range needed for some high energy density experiments. The multi-ion calculations predicted that the hexagonal close-packed phase of Ru remains stable to at least 2.5 TPa and possibly 10 TPa, and that its c/a should gradually increase to the ideal value. A method was devised to estimate the variation in Debye temperature from the cold curve, and thus estimate the ion-thermal EOS without requiring relatively expensive dynamical force calculations, in a form convenient for adjusting EOS or phase boundaries. The Debye temperature estimated in this way was similar to the result from atom-in-jellium calculations. We also predict the high-pressure melt loci of both elements.
Although usually considered as a technique for predicting electron states in dense plasmas, atom-in-jellium calculations can be used to predict the mean displacement of the ion from its equilibrium position in colder matter, as a function of compress
ion and temperature. The Lindemann criterion of a critical displacement for melting can then be employed to predict the melt locus, normalizing for instance to the observed melt temperature or to more direct simulations such as molecular dynamics (MD). This approach reproduces the high pressure melting behavior of Al as calculated using the Lindemann model and thermal vibrations in the solid. Applied to Fe, we find that it reproduces the limited-range melt locus of a multiphase equation of state (EOS) and the results of ab initio MD simulations, and agrees less well with a Lindemann construction using an older EOS. The resulting melt locus lies significantly above the older melt locus for pressures above 1.5,TPa, but is closer to recent ab initio MD results and extrapolations of an analytic fit to them. This study confirms the importance of core freezing in massive exoplanets, predicting that a slightly smaller range of exoplanets than previously assessed would be likely to exhibit dynamo generation of magnetic fields by convection in the liquid portion of the core.
Atom-in-jellium calculations of the electron states, and perturbative calculations of the Einstein frequency, were used to construct equations of state (EOS) from around $10^{-5}$ to $10^7$g/cm$^3$ and $10^{-4}$ to $10^{6}$eV for elements relevant to
white dwarf (WD) stars. This is the widest range reported for self-consistent electronic shell structure calculations. Elements of the same ratio of atomic weight to atomic number were predicted to asymptote to the same $T=0$ isotherm, suggesting that, contrary to recent studies of the crystallization of WDs, the amount of gravitational energy that could be released by separation of oxygen and carbon is small. A generalized Lindemann criterion based on the amplitude of the ion-thermal oscillations calculated using atom-in-jellium theory, previously used to extrapolate melt curves for metals, was found to reproduce previous thermodynamic studies of the melt curve of the one component plasma with a choice of vibration amplitude consistent with low pressure results. For elements for which low pressure melting satisfies the same amplitude criterion, such as Al, this melt model thus gives a likely estimate of the melt curve over the full range of normal electronic matter; for the other elements, it provides a useful constraint on the melt locus.
Atom-in-jellium calculations of the Einstein frequency in condensed matter and of the equation of state were used to predict the variation of shear modulus from zero pressure to ~$10^7$ g/cm$^3$, for several elements relevant to white dwarf (WD) star
s and other self-gravitating systems. This is by far the widest range reported electronic structure calculation of shear modulus, spanning from ambient through the one-component plasma to extreme relativistic conditions. The predictions were based on a relationship between Debye temperature and shear modulus which we assess to be accurate at the o(10%) level, and is the first known use of atom-in-jellium theory to calculate a shear modulus. We assessed the overall accuracy of the method by comparing with experimental measurements and more detailed electronic structure calculations at lower pressures.
سجل دخول لتتمكن من نشر تعليقات
التعليقات
جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
