Next Article in Journal
The Marine Bacterial Genus Euzebya Is Distributed Worldwide in Terrestrial Environments: A Review
Previous Article in Journal
Notch Fatigue Life Research Based on Critical Distance Theory
Previous Article in Special Issue
Magnetocaloric Properties and Critical Behaviour of the Sm2Ni17 Compound
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Principal Hugoniots of Promethium, Terbium, Thulium, Lutetium, and Actinium in a Wide Pressure Range

Los Alamos National Laboratory, P.O. Box 1663, Los Alamos, NM 87545, USA
*
Author to whom correspondence should be addressed.
Appl. Sci. 2023, 13(17), 9643; https://doi.org/10.3390/app13179643
Submission received: 11 April 2023 / Revised: 16 August 2023 / Accepted: 21 August 2023 / Published: 25 August 2023

Abstract

:
We present the analytic forms of the principal Hugoniots of actinium (Ac) and the lanthanide promethium (Pm), which have both never been measured or calculated before, as well as those of terbium (Tb), thulium (Tm), and lutetium (Lu), the three least studied of the remaining lanthanides. They are based on our new analytic model of principal Hugoniot. A comparison of the five Hugoniots to our own independent theoretical calculations demonstrates very good agreement in every case, but each of the Hugoniots of Tb, Tm, and Ac from the TEFIS database, which ours are also compared to, appear to violate Johnson’s theoretical constraint 4 < η max < 7 for the maximum compression ratio η max , which corresponds to the Hugoniot turnaround point. Possible reason for this behavior of the TEFIS Hugoniots is briefly discussed.

1. Introduction

Rare earth 4 f elements, or lanthanides, play a key role in modern day technological applications, such as computer memories and permanent magnets. Similarly, 5 f element actinides, being the backbone of nuclear fission technologies for the production of energy, also find applications in many non-power strategic fields, from space exploration to medical diagnostic and treatments. Specifically, the first member of the 5 f -series, actinium (Ac), is a potential therapeutic agent for cancer and infectious diseases [1,2]. Additionally, rare-earth superhydrides have recently been discovered to have near room temperature T c for superconductors above megabar pressures (1 Mbar = 100 GPa) [3,4,5]. As such, understanding and predicting the material behavior of lanthanides is crucial. This can be done using equation of state (EOS) information; however, both 4 f -electron lanthanides and 5 f -electron actinides present a unique challenge in performing first-principles based ab initio calculations due to their complex electronic and lattice structures. Therefore, experimental EOS measurements are indispensable in providing benchmark data to tightly constrain the ab initio calculations. Research efforts on Ac, as well as some lanthanides, have been hampered by the lack of supply and the high costs of current production methods, especially for Ac, and the radioactivity of promethium (Pm). In fact, Pm is the only radioactive element among the lanthanides, and the second radioactive element after technetium (Tc) that has stable neighboring elements. Thus, the EOS data on both Pm and Ac are scarce and virtually non-existent. Lutetium (Lu), in contrast, has been the subject of intense studies, mainly because it is the last element in the lanthanide series and has a completely filled 4 f shell. However, even with the availability of several EOS studies [6], there is only one published study of the shock compression of Lu, by Al’tshuler et al. [7]. Similarly, there is only one published study of the shock compression of both terbium (Tb) and thulium (Tm), by Carter et al. [8].
The shock Hugoniot, both the principal and second-shock one, is a continuous curve which describes the locus of all possible thermodynamic states a substance can exist in behind a shock propagating through it. It is usually given as a projection onto a two-dimensional plane and is described in terms of two coordinates which can be chosen among particle velocity behind the shock front ( U p ) , shock velocity ( U s ) , pressure, internal energy ( E ) , temperature ( T ) , density ( ρ ) , and degree of compression ( η ) called “compression” for simplicity. Up to P of a few hundred GPa, the Hugoniot is known from experiments, and in the vast majority of cases is described by a linear relation U s = C + B U p , between the two main variables, U s and U p , in terms of which all other relevant quantities can be expressed via the Rankine–Hugoniot (RH) relations (which represent the conservation of mass, momentum, and energy, respectively):
η ρ ρ 0 = U s U s U p , P = P 0 + ρ 0 U s U p , E E 0 = 1 2 ( P + P 0 ) 1 ρ 0 1 ρ ,
where the subscript 0 indicates the initial (unshocked) state. For the principal Hugoniot, P 0 = 0 .
Over decades of research, shock Hugoniot measurements have been one of the most important sources of EOS data. The use of planar shock waves to determine the EOS of condensed materials to very high pressure ( P ) began in 1955 with the classic papers of Walsh and Christian (1955) [9] and Bancroft et al. (1956) [10]. Walsh and Christian described the use of in-contact explosives to determine dynamic pressure-volume relations for metals and compared these to the then available static compression data. Bancroft et al. described the first polymorphic phase change discovered in a solid via shock waves, namely, iron. Two years later, Al’tshuler et al. (1958) [11] reported the first data for iron to P of several Mbar, essentially in excess of P in the center of the Earth. Since that time, the EOSs of virtually hundreds or even thousands of condensed materials have been studied, including elements, compounds, alloys, rocks and minerals, polymers, fluids, and porous media. These studies have employed both conventional and nuclear explosive sources, as well as impactors launched with a range of guns to speeds of order of 10 km/s. Recently, with the establishment of the NIF Gigabar platform at the Lawrence Livermore National Laboratory [12], pressures of order of 1 Gbar ( 10 5 GPa, or 100 TPa), deep into the atomic pressure regime, have been routinely achieved in their experiments [13] (in this regime, P exceeds P a , , the so-called atomic pressure required to significantly distort core electron orbitals and estimated as P a = E H / r Bohr 3 = 294 Mbar, where E H is the Hartree energy and r Bohr is the Bohr radius).
The lanthanides series includes lanthanum (La) and the 14 following elements in which the 4 f -electron shell gets progressively filled. This shell is half-filled for europium (Eu) and fully filled for ytterbium (Yb), which is responsible for both of them being divalent, in contrast to the remaining lanthanides, which are all trivalent, similar to their transition metal counterparts, scandium (Sc) and yttrium (Y). Because the 4 f -electrons are localized deep in the lanthanide atoms (specifically, in their Xe-cores), they do not participate in metallic bonding [14], which is dominated by the 5 p , 5 d , and 6 s valence electrons. The typical electron structure of a lanthanide atom is [Xe] 4 f n 5 d 0 6 s 2 or [Xe] 4 f n 1 5 d 1 6 s 2 .
The dominance of the s p d -valence electrons in the bonding is responsible for the lanthanides’ phase diagrams to vary slowly across the series so that their relations to one another is clearly seen [15]. The well established fact about the lanthanides is that, with increasing P , virtually all of them undergo the same sequence of phase transformations from one hexagonal polytype to the next; it is characterized by the varying hexagonal layer stacking of the corresponding crystal structure: hexagonal close-packed (hcp, AB) → α -Sm (9R, ABCBCACAB) → double-hcp (dhcp, ABAC) → face-centered cubic (fcc, ABC) → distorted fcc (dfcc) →… Evans et al. [16] showed that the dfcc phase of praseodymium (Pr) is h R 24 , another hexagonal polytype with 24 atoms per unit cell, i.e., with a stacking of 24 ABC-layers. h R 24 , along with two alternative crystal structures, namely, orthorhombic o I 16 , and o S 8 , were identified as dfcc for other lanthanides. Recently, a possible extension of the above phase transition sequence to include post-dfcc phases was suggested in Ref. [17].
The actinides form a similar series, from Ac to lawrencium (Lr), in which the 5 f -electron shell gets progressively filled. In this respect, a comparison of the phase diagrams of the two series is useful and informative, but in contrast to the lanthanides, the structural systematics for the actinides, in terms of a sequence of phase transformations on increasing P, is much less obvious for the light actinides, but becomes more obvious for the heavier actinies, starting with americium (Am) [18]: dhcp → fcc → α -U →… Note that, again, in contrast to the lanthanides, α -U is not one of the hexagonal polytypes.
The principal Hugoniots of lanthanides have been experimentally investigated in several studies, e.g., refs. [7,19], but only a few substances at a time. The systematic study of the Hugoniots for the lanthanide series as a whole was carried out in Refs. [8,20]. Specifically, 13 lanthanides out of 15 were studied in [8] (in addition to Sc and Y), and 9 in [20] (in addition to Sc, Y, and hafnium (Hf)). Both sets of data are generally consistent with each other: the numerical values of the parameters of the least-square linear fits to the corresponding U s - U p data agree to within ∼15%. Here, U p and U s are particle and shock velocities, respectively. The unique feature of both sets of data is that for each substance, except cerium (Ce), the adequate description of the data in the whole U p region of up to ∼3.6 km/s requires a two-segment fit; that is, the lower and upper portions of the data sets are described by two distinct linear segments such that the combined two-segment graph is continuous, but not smooth, as the slope of the second segment (∼1.2–1.6) is ∼1.5 times higher than that of the first one (∼0.8–1.0). In each case, the change of slope happens at a compression of ∼1.6; the corresponding transition P increases from ∼20 GPa for the light lanthande members to ∼40 GPa for the heavy ones [8]. Eu and Yb are two exceptions, for which the transition P is ∼10 GPa. Additionally, the Eu data show the existence of a two-wave structure in the shock front, indicating a significant volume change at the transition, whereas in most of other cases, the data indicate that the volume change associated with this transition is either zero or very small [8]. Such U s - U p behavior is typical of a solid–solid phase transformation or melting. In any event, the U s - U p slope being ∼50% higher indicates that the compressibility of the emerging phase is changing more rapidly (the value of the P-derivative of the bulk modulus is higher). Direct evidence of phase transitions and melting along the principle Hugoniot in several lanthanide members using in-situ laser shock diffraction at the Dynamic Compression Sector was recently presented in [21].
In this work, we derive the analytic forms of the principal Hugoniots of Ac and Pm, which have both never been measured or calculated before, as well as those of Tb, Tm, and Lu, the three least studied of the remaining lanthanides. Indeed, the existing Hugoniot data on each of these three substances can only be found in a single literature source, namely, Ref. [8] for both Tb and Tm and Ref. [7] for Lu.

2. Principal Hugoniot in a Wide Pressure Range

We will construct the analytic models of the principal Hugoniots of the five substances, namely, Pm, Tb, Tm, Lu, and Ac, using the analytic framework established in our previous publication [22]. In this framework, a wide P range is divided into three regimes and the Hugoniot is constructed in each of these regimes and then interpolated smoothly between them. These regimes are: (i) the low-P regime in which the Hugoniot is described by U s = C + B U p , where the values of C and B come from the experiment; (ii) the intermediate-P regime (discussed in more detail below) where the Hugoniot is described by the Thomas–Fermi–Kalitkin (TFK) model [23,24] U s = c + b U p + a U p 2 , with the values of c , b, and a determined virtually for all Zs (Z being the atomic number) [23,24,25]; and (iii) the high-P regime in which the Hugoniot is described by the Debye-Hückel-Johnson (DHJ) model [26,27,28]. The only assumption made was that the principal Hugoniot is goverened by some function U s = U s ( U p ) (which is linear at low P and quadratic at intermediate P ; its high-P form was established in [22]) and that this function is continuous and smooth (the first derivative d U s / d U p is continuous) at all U p . Then, it follows from the RH relations [22] that η , P, and E are all continuous and smooth as well. No other assumption, and no additional free parameter except the six mentioned above, C ,   B ,   c ,   b ,   a , and Z , are required for the construction of the analytic model of the principal Hugoniot. In order to match the next regime, the linear form of the low-P regime is modified into U s = C + B U p + A U p 2 , where A is an additional 7th parameter which introduces a very small non-linearity and is obtained using the formula A = a ( B b ) 2 / ( 4 ( c C ) ) [22]. Since a = O ( 10 5 s / km ) and c C = O ( 1 km / s ) (see, e.g., Table 1 below), the value of A depends on B b . Typically, this difference is of the order of 0.1; hence, A 10 2 s/km. However, if this difference is small, so is A , but U p * in Equation (1), which determines the low-P–med-P transition point, is large, thus pushing the upper boundary of the low-P regime closer to the turnaround point. This is, e.g., the case of thulium considered in what follows for which U p * 160.5 km/s is 2 / 3 of U p max 233.4 km/s at the turnaround point (see Table 1); typically, U p * 0.1 U p max .
The names of the three P regimes of the principal Hugoniot may sound confusing and may not correspond to those adapted in high-pressure research in general and phase diagram and EOS studies in particular. Specifically, it is generally adopted that high-P corresponds to pressures in excess of ∼100 GPa. In our case, the low-P–med-P transition point corresponds to shock velocities of ∼10 km/s (see, e.g., Table 1); with ambient density of ∼1–10 g/cm 3 , corresponding to a pressure of ∼100–1000 GPa. Thus, our med-P regime is analogous to the more familiar high-P range of EOS studies. Our med-P–high-P transition point corresponds to P 1–10 Gbar ( 10 5 10 6 GPa).

3. Methods

3.1. New Analytic Model for Principal Hugoniot

Here is a brief summary of the new analytic model, as per Ref. [22]. The new model is based on the following representation of the shock velocity U s as a function of the particle velocity U p over the three U p intervals, which is continuous and smooth from one interval to the next:
U s = C + B U p + A U p 2 , 0 U p U p * = 2 ( c C ) B b U s = c + b U p + a U p 2 , U p * U p U p * * = n U p * U s = d | f | + 4 3 U p + d U p 1 + | f | U p , U p U p * *
where | f | stands for the absolute value of f .
The corresponding expressions for P along the principal Hugoniot are:
P = 4 ρ 0 C 2 η ( η 1 ) { η B ( η 1 ) + [ η B ( η 1 ) ] 2 4 A C ( η 1 ) 2 } 2 , 0 U p U p * P = 4 ρ 0 c 2 η ( η 1 ) { η b ( η 1 ) + [ η b ( η 1 ) ] 2 4 a c ( η 1 ) 2 } 2 , U p * U p U p max = c a P = 4 ρ 0 c 2 η ( η 1 ) { η b ( η 1 ) [ η b ( η 1 ) ] 2 4 a c ( η 1 ) 2 } 2 , U p max U p U p * * P = ρ 0 4 | f | 2 η η 1 12 d η 1 η 4 + 1 1 2 , U p U p * *
which define pressure along the Hugoniot as a continuous and smooth function of compression ratio η ρ / ρ 0 , where ρ is density and the subscript 0 indicates the initial (unshocked) state. In the last equation of (3), the “−” sign corresponds to f > 0 and the “+” sign to f < 0 . Whereas d is always positive, it appears that f is positive for about 2/3 of the periodic table and negative for the remaining 1/3. In particular, f is negative for the entire lanthanide series, as well as for the early actinides, including Ac. Therefore, in our study, we use this equation with the “+” sign. The numerical values of d and f are defined by the formulas [22]:
f = a [ ( 4 3 b ) c a n 3 c ( 1 + n 2 ) ] 3 a [ L + c c a n ( 1 + n 2 ) ] c ( 4 3 b ) n 2 , d = L a [ ( 4 3 b ) c a n 3 c ( 1 + n 2 ) ] 3 a [ L + c c a n ( 1 + n 2 ) ] c ( 4 3 b ) n 2 2 ,
with n = U p * * / U p max , U p * * being the (only positive) solution of Equation (5) below. The functions of Z , d, and f are shown, respectively, in Figure 1 and Figure 2, and are not presented in [22].
In Equations (2) and (3), 0 U p U p * corresponds to the low-P regime, and U p * U p U p * * and U p U p * * to the med-P and high-P regimes, respectively. The value of U p * * is found from Equation [22]:
c + b 4 3 U p * * + a U p * * 2 2 = L b 4 3 + 2 a U p * * ,
where [22]:
L = 875 Z 2.4 A 875 Z 2.4 2 Z + 0.006 Z 2 km s 2 ,
where A is the atomic mass.

3.2. The Choice of the Five Sets of the Hugoniot Parameters

Whereas the parameters c , b, and a for the med-P regime are available for all substances, and, therefore, d and f for the high-P regime can be calculated for all substances from Equation (4) with n = U p * * / U p max , where U p max = c / a and U p * * satisfies Equation (5), the parameters C and B for the low-P regime must be specified (and then A = a ( B b ) 2 / ( 4 ( c C ) ) [22]). Here, we use the original version of the analytic Hugoniot model [22], which describes the entire low-P regime in terms of a single ( C , B , A ) parameter set, regardless of whether the actual low-P portion of the Hugoniot is multi-segment (two or more linear segments forming a continuous broken line, just as for the lanthanides), or multi-linear (two or more quasi-parallel segments, which is the case of, e.g., zirconium [29]), or a combination of both. In our case of the four lanthanides, it is logical to describe the low-P regime in terms of the corresponding parameters of the upper segment of the two-segment portion in order to match the upper segment directly with the med-P regime. These values come from [8] for Tb and Tm, and from [7] for Lu. Since no ( C , B ) data exist on Pm, we obtain both values from smooth interpolation of the corresponding values for other members across the entire lanthanide series based on the results of [8,20]. For Ac, for which no ( C , B ) data exist either, assuming that it is a single-phase material (in fact, Ac is fcc, and there is no evidence for a phase transition to any other solid phase [18]), two parameters can be calculated as [26,27] C = B s ( ρ 0 ) / ρ 0 and B = ( B s ( ρ 0 ) + 1 ) / 4 , with bulk modulus B s ( ρ 0 ) and its P-derivative B s ( ρ 0 ) coming from the corresponding EOS. The best fit to the ab initio data on the EOS of Ac of Ref. [30] of the third-order Birch–Murnaghan form results in B s ( ρ 0 ) = 23.7 ± 3.7 GPa and B s ( ρ 0 ) = 3.9 ± 0.3 ; therefore, with ρ 0 = 10.05 g/cm 3 , C = 23.7 / 10.05 = 1.535 km/s, and B = ( 3.9 + 1 ) / 4 = 1.225 .  Table 1 contains all the five complete sets of the Hugoniot parameters.

4. Results

The five principal Hugoniots of the substances discussed in this work are shown in Figure 3, Figure 4, Figure 5, Figure 6 and Figure 7, in which the five Hugoniots are compared to both independent theoretical calculations and the available experimental data.
In addition to the available experimental data (for Tb, Tm, and Lu), for the purpose of comparison to the new analytic model, we have carried out theoretical Hugoniot calculations using: (i) the relativistic Green’s function quantum average atom code Tartarus [31,32] for Pm and Lu; and (ii) the Thomas–Fermi model with corrections [33] for Tb, Tm, and Ac, for which Lambert’s orbital-free molecular dynamics (OFMD) code used, e.g., in a study on the transport properties of lithium hydride [34], was modified appropriately. These simulations are similar to those of Ref. [35] for platinum. The ( V , P ) points that map out the Hugoniot (presented in the η -P coordinates upon conversion of the ( V , P ) points into the ( η , P ) ones) are found from the RH relation for internal energy: E E 0 = P ( V 0 V ) / 2 . We also used three Hugoniots from the TEFIS database [36,37] for Tb, Tm, and Ac for the comparison to the new model, as well as to the Thomas–Fermi model with corrections. We note that, although Lu is an example of the application of Tartarus to a real material considered in detail in Ref. [32], the Lu Hugoniot shown in Figure 4 is not presented in [32].
Similar to P = P ( η ) , E along the Hugoniot can be calculated using the corresponding RH relation (1) for E ,
E = E 0 + P 2 ρ 0 1 1 η ,
and Equations (3) for P = P ( η ) in the corresponding U p intervals. Figure 8 shows E = E ( η ) for Ac (we assumed E 0 = 0 , for simplicity).
It is clearly seen that the forms of P = P ( η ) and E = E ( η ) are very similar to each other. The form of T = T ( η ) is also expected to be similar to both P = P ( η ) and E = E ( η ) . However, the accurate calculation of T = T ( η ) requires the knowledge of the specific heat along the Hugoniot; hence, this calculation goes well beyond the scope of this work. We can, however, estimate the value of T , based on that of E , and compare it to the corresponding value from the TEFIS tables. For an ideal gas (which the system well above the turnaround point represents), E 3 / 2 k B T , at the upper limit of E = E ( η ) k B T 2 / 3 E . The uppermost TEFIS point in Figure 7 corresponds to η 5.04 . At this η , the TEFIS table of T = T ( η ) gives T = 10 4 eV (a table of E = E ( η ) does not exist), and our value of E is E ( η = 5.04 ) 3.8 × 10 4 eV/mol; therefore, the above relation gives T 2.5 × 10 4 eV. Thus, our value of T is about 2.5 times larger than TEFIS’ one, which directly corresponds with the fact that our P is larger than TEFIS’ one by roughly the same amount, as can be clearly seen in Figure 7.

5. Discussion

As Figure 3, Figure 4, Figure 5, Figure 6 and Figure 7 clearly demonstrate, in each of the five cases, the agreement between the new analytic model and independent data is very good except for the three Hugoniots from the TEFIS database. As a matter of fact, the three TEFIS Hugoniots appear to violate Johnson’s constraint 4 < η max < 7 , which is based on rigorous theoretical grounds. (We have determined that this is the case for the vast majority of TEFIS Hugoniots for other substances.)
It is worthwhile to dwell on Johnson’s constraint in some more detail. In Ref. [38], Johnson derives the formula for the maximum compression on the princial Hugoniot:
η max = 4 ( 1 + 7 C ) 1 + 4 C ,
where C 0.011 A Z 4.2 / [ ρ 0 ( 1 + Z ) 4 ] . As can be clearly seen in Figure 4 of [22], this η max and the one given by the new analytic model are in very good agreement with each other, which lends further support to both formulations. Johnson’s constraint now follows directly from (8):
η max = 4 + 28 C 1 + 4 C < 7 + 28 C 1 + 4 C = 7 ,
and:
η max = 4 + 28 C 1 + 4 C > 4 + 16 C 1 + 4 C = 4 .
We believe, the reason for this behavior of TEFIS Hugoniots is the way these Hugoniots are constructed; specifically, the scheme used in Refs. [39,40] for the interpolation between the portion of the Hugoniot below and across the turnaround point described by the TFK model, and that above the turnaround point, which incorporates electron shell effects, such shell effects are clearly seen in Figure 4, Figure 5, and Figure 7 for Tb, Tm, and Ac, respectively. In other words, the violation of Johnson’s theoretical constraint by the TEFIS Hugoniots may be the artifact of the interpolation scheme used for their construction. Otherwise, for P below and above the turnaround point, 1 < η < 4 and 4.5 < η < 6 , agreement between the TEFIS Hugoniots and both the new analytic model and the Thomas–Fermi model with corrections is generally good.
Let us note that in a more realistic case of a shock compression of a substance beyond the corresponding turnaround point, the Hugoniot must be modeled by taking into account the well-known effects of the contribution of both the equilibrium radiation of hot plasma and relativistic effects [41,42,43,44]. Indeed, at a turnaround point, U p 200 –250 km/s (see Table 1), which constitutes < 0.1 % of the speed of light, and at the med-P–high-P transition point, U p 1700 ± 200 km/s, about 0.5% of the speed of light. Hence, as the system enters the high-P regime, relativistic effects are expected to start manifesting themselves and eventually to dominate the evolution of the system at even higher U p . Our model can, in principle, be modified to incorporate these effects. For instance, taking into account the radiation-dominated (or the so-called strong shock) regime can be done by replacing equations for U s ( U p ) and P ( η ) describing the high-P regime with their counterparts stemming from the physics of a photon gas. This goes beyond the scope of the present work, but will be undertaken in one of our subsequent studies.
The new model discussed in this work does not incorporate potential electronic shell effects in the high-P regime. If present, they manifest themselves in terms of some “irregularities”, that is, changes of the sign of d P / d η of the continuous line P = P ( η ) over some (small) regions of P ; such “irregularities” are seen, e.g., in the Tartarus Hugoniots of both Pm and Lu around 10 5 10 6.5 GPa in Figure 3 and Figure 6, respectively, since Tartarus takes electronic shell effects into account explicitly. Similar irregularities are present in the TEFIS Hugoniots, but the corresponding changes of the sign of d P / d η are very abrupt. This is because these Hugoniots are constructed the way that causes the violation of Johnson’s theoretical constraint, as discussed above. In contrast to Tartarus, the orbital-free procedure of the OFMD code treats all electrons on an equal footing, albeit approximately, with no distinction between bound and ionized electrons. This is why any shell-ionization related effects are absent in the corresponding Hugoniots, as seen in Figure 4, Figure 5, and Figure 7. Such effects cannot be predicted by the new model, but, if firmly established, they can be added to the model by considering additional region(s) of P described by the corresponding U s = U s ( U p ) functional forms. We plan to undertake such an addition of electronic shell effects to the new model in one of our subsequent studies on this subject.

6. Conclusions

Here is a brief summary of the findings of this work. We have presented the principal Hugoniots of actinium and the lanthanide promethium, which have both never been measured or calculated before, as well as those of terbium, thulium, and lutetium, the three least studied of the remaining lanthanides. We used the analytic framework established in our previous publication [22]. All five sets of the relevant parameters are summarized in Table 1. The five principal Hugoniots are compared to the available experimental data (for Tb and Tm only) and independent theoretical calculations in Figure 3, Figure 4, Figure 5, Figure 6 and Figure 7. As can be clearly seen in these figures, in each of the five cases, the agreement between the new analytic model and independent data is very good except for the three Hugoniots from the TEFIS database, which all violate Johnson’s theoretical constraint on the maximum compression. Their behavior is analyzed in some detail in our paper. We anticipate that new experimental measurements could shed more light on a potential systematics of the lanthanide Hugoniots in terms of the values of the parameters C and B , the maximum compression η max , etc. We expect our results to serve as the initial guidance for such experiments.
Our new analytic model of the principal Hugoniot [22] can be used for the validation of the P-V-T EOS by comparing the Hugoniot produced by the EOS to that given by the model. Additionally, the new model itself can be used as a basis for EOS construction. Indeed, if the Grüneisen parameter along the Hugoniot is available, e.g., using the approach discussed in this work, then it can be used in the Mie-Grüneisen-type EOS P P H = γ H ρ ( E E H ) [45], where the subscript “H” implies that the corresponding variable is in shock-compression conditions. This EOS can then be brought in direct correspondence to the more familiar Mie-Grüneisen (M-G) EOS, P P c = γ ρ ( E E c ) , where the subscript “c” implies the cold ( T = 0 ) conditions, since there exists a direct algebraic connection between γ H (of this work) and γ of the M-G EOS [45].
The analytic model developed in our previous study [22] and applied to five substances in this work can be used to calculate of the Hugoniots of other substances, specifically, other lanthanides and/or actinides. In this respect, the analytic knowledge of the regimes of the Hugoniot past the turnaround point is very important. In a very recent publication [46], a team of astronomers report a detailed study of a pair of shock waves produced by a collision of two clusters of galaxies that occurred roughly a billion years ago. The shocks that are associated with cluster mergers are known as radio relics and they can be used to probe the properties of the intergalactic space within the cluster, known as the intracluster medium, as well as intracluster dynamics. The study focused on a particular cluster called Abell 3667, which at least 550 galaxies are associated with and which is still coming together. It was concluded that the shock waves are propagating through it at velocities of ∼1500 km/s, which are 5–6 times larger than the velocities corresponding to turnaround points [22], U s max = 2 c + b c / a b c / a 1.2 U p max 250–300 km/s (see Table 1). Moreover, at such shock velocities, many elements, especially the low-Z ones, will be under the conditions that are beyond the validity of Kalitkin’s parabolic representation (5), i.e., in the high-P regime considered in this work. Hence, to predict the properties of the intracluster medium and to describe intracluster dynamics, an analytic model of the principal Hugoniot in the high-P regime is a must. Once the analytic formulas describing the high-P regime are available (the last lines of both systems of Equations (2) and (3)), the basic mechanical and thermodynamic properties of a material under intergalactic shock, such as the bulk modulus, the Grüneisen parameter, energy, temperature, etc., along the principal Hugoniot can be derived from the RH relations and Equation (16) using the new model and some additional assumptions on the Grüneisen γ .

Author Contributions

Methodology, L.B., D.L.P., S.D.R., S.K.S., C.E.S. and R.S.B.; investigation, L.B., D.L.P., S.K.S. and C.E.S.; writing—original draft preparation, L.B.; writing—review and editing, D.L.P., S.D.R., S.K.S., C.E.S. and R.S.B.; project administration, S.D.R. and R.S.B.; funding acquisition, S.D.R. and R.S.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research was carried out under the auspices of the US DOE/NNSA.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All the data discussed in this work are available from the corresponding author upon request.

Acknowledgments

The OFMD simulations were performed on the LANL clusters Badger and Snow.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Morgenstern, A.; Apostolidis, C.; Kratochwil, C.; Sathekge, M.; Krolicki, L.; Bruchertseifer, F. An overview of targeted alpha therapy with 225actinium and 213bismuth. Curr. Radiopharm. 2018, 11, 200. [Google Scholar] [CrossRef]
  2. Actinium Presents New Data Demonstrating Effective Lymphodepletion with Lutetium-177 for CAR-T at the 2019 Society of Nuclear M; Bloomberg: New York, NY, USA, 2019; Available online: https://www.bloomberg.com/press-releases/2019-06-25/actinium-presents-new-data-demonstrating-effective-lymphodepletion-with-lutetium-177-for-car-t-at-the-2019-society-of-nuclear-m (accessed on 10 April 2023).
  3. Drozdov, A.P.; Kong, P.P.; Minkov, V.S.; Besedin, S.P.; Kuzovnikov, M.A.; Mozaffari, S.; Balicas, L.; Balakirev, F.F.; Graf, D.E.; Prakapenka, V.B.; et al. Superconductivity at 250 K in lanthanum hydride under high pressures. Nature 2019, 569, 528. [Google Scholar] [CrossRef] [PubMed]
  4. Sun, D.; Minkov, V.S.; Mozaffari, S.; Sun, Y.; Ma, Y.; Chariton, S.; Prakapenka, V.B.; Eremets, M.I.; Balakirev, L.B.; Balakirev, F.F. High-temperature superconductivity on the verge of a structural instability in lanthanum superhydride. Nat. Commun. 2021, 12, 6863. [Google Scholar] [CrossRef]
  5. Li, Z.; He, X.; Zhang, C.; Wang, X.; Zhang, S.; Jia, Y.; Feng, S.; Lu, K.; Zhao, J.; Zhang, J. Superconductivity above 200 K discovered in superhydrides of calcium. Nat. Commun. 2022, 13, 2863. [Google Scholar] [CrossRef] [PubMed]
  6. Samudrala, G.L.; Vohra, Y.K. Structural properties of lanthanides at ultra high pressure. In Handbook on the Physics and Chemistry of Rare Earths; Elsevier B.V.: Amsterdam, The Netherlands, 2013; Volume 43, p. 275. [Google Scholar]
  7. Al’tshuler, L.V.; Bakanova, A.A.; Dudoladov, I.P. Effect of electron structure on the compressibility of metals at high pressure. Sov. Phys. JETP 1968, 26, 1115. [Google Scholar]
  8. Carter, W.J.; Fritz, J.N.; Marsh, S.P.; McQueen, R.G. Hugoniot equation of state of the lanthanides. J. Phys. Chem. Solids 1975, 36, 741. [Google Scholar] [CrossRef]
  9. Walsh, J.M.; Christian, R.R. Equation of state of metals from shock wave measurements. Phys. Rev. 1955, 97, 1544. [Google Scholar] [CrossRef]
  10. Bancroft, D.; Peterson, E.L.; Minshall, S. Polymorphism of iron at high pressure. J. Appl. Phys. 1956, 27, 291. [Google Scholar] [CrossRef]
  11. Al’tshuler, L.V.; Krupnikov, K.; Ledenev, B.N.; Zhuchikhin, V.I.; Brazhnik, M.I. Dynamic compressibility and equation of state of iron under high pressure. Soviet Phys. JETP 1958, 34, 606. [Google Scholar]
  12. Fratanduono, D.E.; Swift, D.; Jenei, A.; Whitley, H. Accessing extreme equation of state conditions on the national ignition facility, LLNL Preprint LLNL-JRNL-739795, 11 October 2017. Available online: https://www.osti.gov/servlets/purl/1455406 (accessed on 10 April 2023).
  13. Kritcher, A.L.; Swift, D.C.; Döppner, T.; Bachmann, B.; Benedict, L.X.; Collins, G.W.; DuBois, J.L.; Elsner, F.; Fontaine, G.; Gaffney, J.A. A measurement of the equation of state of carbon envelopes of white dwarfs. Nature 2020, 584, 51. [Google Scholar] [CrossRef]
  14. Min, B.I.; Jansen, H.J.F.; Oguchi, T.; Freeman, A.J. Local density total energy description of ground and excited state properties of the rare earth metals. J. Magn. Magnet. Mater. 1986, 61, 139. [Google Scholar] [CrossRef]
  15. Johansson, B.; Rosengren, A. Generalized phase diagram for the rare-earth elements: Calculations and correlations of bulk properties. Phys. Rev. B 1975, 11, 2836. [Google Scholar] [CrossRef]
  16. Evans, S.R.; Loa, I.; Lundegaard, L.F.; McMahon, M.I. Phase transitions in praseodymium up to 23 GPa: An X-ray powder diffraction study. Phys. Rev. B 2009, 80, 134105. [Google Scholar] [CrossRef]
  17. McMahon, M.I.; Finnegan, S.E.; Pace, E.J.; Storm, C.V.; Stevenson, M.G.; Macleod, S.G.; Plekhanov, E.; Bonini, N.; Weber, C. New structural systematics in the lanthanide elements at high pressure. Bull. Mater. Sci. 2022, 45, 186. [Google Scholar] [CrossRef]
  18. Young, D.A. Phase Diagrams of the Elements; University of California Press: Berkeley, CA, USA, 1991; 232p. [Google Scholar]
  19. Al’tshuler, L.V.; Bakanova, A.A.; Dudoladov, I.P.; Dynin, E.A.; Trunin, R.F.; Checkin, B.S. Shock adiabats for metals. New data, statistical analysis and general regularities. J. Appl. Mech. Tech. Phys. 1981, 22, 145. [Google Scholar]
  20. Gust, W.H.; Royce, E.B. New electronic interactions in rare-earth metals at high pressure. Phys. Rev. B 1973, 8, 3595. [Google Scholar] [CrossRef]
  21. Sakun, D.; McCoy, C.; Melton, C.; Knudson, M.; Root, S.; Seagle, C. Structural Evolution of f-Electron Systems under Dynamic Compression. In Proceedings of the APS March Meeting 2022, Chicago, IL, USA, 14–18 March 2022; Available online: https://meetings.aps.org/Meeting/MAR22/Content/4178 (accessed on 10 April 2023).
  22. Burakovsky, L.; Preston, D.L.; Ramsey, S.D.; Baty, R.S. Analytic model of principal Hugoniot at all pressures. J. Appl. Phys. 2022, 132, 215109. [Google Scholar] [CrossRef]
  23. Kalitkin, N.N.; Kuzmina, L.V. Quantum-statistical Hugoniots of porous substances. Mat. Model. 1998, 10, 111. [Google Scholar]
  24. Kalitkin, N.N.; Kuzmina, L.V. Wide-Range Characteristic Thermodynamic Curves. In High-Pressure Shock Compression of Solids VII: Shock Waves and Extreme States of Matter; Fortov, V.E., Al’tshuler, L.V., Trunin, R.F., Funtikov, A.I., Eds.; Springer: New York, NY, USA, 2004; p. 116. [Google Scholar]
  25. Kalitkin, N.N.; Kuzmina, L.V. Shock Hugoniots of 83 Substances. Chem. Phys. Rep. 2000, 18, 1913. [Google Scholar]
  26. Johnson, J.D. General Features of Hugoniots, LANL Preprint LA-13137-MS, UC-705. April 1996. Available online: https://www.osti.gov/servlets/purl/230530 (accessed on 10 April 2023).
  27. Johnson, J.D. General Features of Hugoniots—II, LANL Preprint LA-13217-MS, UC-910. January 1997. Available online: https://www.osti.gov/servlets/purl/437718 (accessed on 10 April 2023).
  28. Johnson, J.D. The features of the principal Hugoniot. AIP Conf. Proc. 1998, 429, 27. [Google Scholar]
  29. Greef, C.W. Phase changes and the equation of state of Zr. Model. Simul. Mater. Sci. Eng. 2005, 13, 1015. [Google Scholar] [CrossRef]
  30. Young, D.A.; Cynn, H.; Söderlind, P.; Landa, A. Zero-Kelvin compression isotherms of the elements 1 ≤ Z ≤92 to 100 GPa. J. Phys. Chem. Ref. Data 2016, 45, 043101. [Google Scholar] [CrossRef]
  31. Gill, N.M.; Starrett, C.E. Tartarus: A relativistic Green’s function quantum average atom code. High Energy Density Phys. 2017, 24, 3338. [Google Scholar] [CrossRef]
  32. Starrett, C.E.; Gill, N.M.; Sjostrom, T.; Greeff, C.W. Wide ranging equation of state with Tartarus: A hybrid Green’s function/orbital based average atom code. Comput. Phys. Commun. 2019, 235, 50. [Google Scholar] [CrossRef]
  33. Nikiforov, A.F.; Novikov, V.G.; Uvarov, V.B. Quantum-Statistical Models of Hot Dense Matter. Methods for Computation Opacity and Equation of State; Birkhäuser: Basel, Switzerland, 2005. (In Russian) [Google Scholar]
  34. Burakovsky, L.; Ticknor, C.; Kress, J.D.; Collins, L.A.; Lambert, F. Transport properties of lithium hydride at extreme conditions from orbital-free molecular dynamics. Phys. Rev. E 2013, 87, 023104. [Google Scholar] [CrossRef] [PubMed]
  35. Elkin, V.M.; Mikhaylov, V.N.; Ovechkin, A.A.; Smirnov, N.A. A wide-range multiphase equation of state for platinum. J. Phys. Cond. Mat. 2020, 32, 435403. [Google Scholar] [CrossRef]
  36. Ivanchenko, E.S.; Kalitkin, N.N.; Kuz’mina, L.V. Main Hugoniot Adiabats in the Tefis Database of Thermophysical Properties of Substances (TEFIS). Math. Models Comput. Simul. 2009, 1, 383. [Google Scholar] [CrossRef]
  37. Belov, A.A.; Golovanov, R.V.; Kalitkin, N.N.; Kozlitin, I.A.; Koriakin, P.V.; Kuzmina, L.V. TEFIS Database. Thermodynamic Properties of Substances, Keldysh IAM Preprint 219, Moscow. 2018. Available online: https://keldysh.ru/papers/2018/prep2018_219.pdf (accessed on 10 April 2023).
  38. Johnson, J.D. Bound and estimate for the maximum compression of single shocks. Phys. Rev. E 1999, 59, 3727. [Google Scholar] [CrossRef]
  39. Ivanchenko, E.S.; Kalitkin, N.N.; Kuzmina, L.V. Quantum-statistical calculations of thermodynamics and Hugoniots (codes QSM). Math. Model. 2008, 20, 30. [Google Scholar]
  40. Ivanchenko, E.S.; Kalitkin, N.N. Shock wave data interpretation in the TEFIS database. Math. Models Comput. Simul. 2010, 2, 526. [Google Scholar] [CrossRef]
  41. Kadatskiy, M.A.; Khishchenko, K.V. Comparison of Hugoniots calculated for aluminum in the framework of three quantum-statistical models. J. Phys. Conf. Ser. 2015, 653, 012079. [Google Scholar] [CrossRef]
  42. Kadatskiy, M.A.; Khishchenko, K.V. Shock compressibility of iron calculated in the framework of quantum-statistical models with different ionic parts. J. Phys. Conf. Ser. 2016, 774, 012005. [Google Scholar] [CrossRef]
  43. Kadatskiy, M.A.; Khishchenko, K.V. Theoretical investigation of the shock compressibility of copper in the average-atom approximation. Phys. Plasmas 2018, 25, 112701. [Google Scholar] [CrossRef]
  44. Orlov, N.Y.; Kadatskiy, M.A.; Denisov, O.B.; Khishchenko, K.V. Application of quantum-statistical methods to studies of thermodynamic and radiative processes in hot dense plasmas. Matter Radiat. Extrem. 2019, 4, 054403. [Google Scholar] [CrossRef]
  45. Ramsey, S.D.; Schmidt, E.M.; Boyd, Z.M.; Lilieholm, J.F.; Baty, R.S. Converging shock flows for a Mie-Grüneisen equation of state. Phys. Fluids 2018, 30, 046101. [Google Scholar] [CrossRef]
  46. Fde, G.; Rudnick, L.; Finoguenov, A.; Wittor, D.; Akamatsu, H.; Brüggen, M.; Chibueze, J.O.; Clarke, T.E.; Cotton, W. MeerKAT view of the diffuse radio sources in Abell 3667 and their interactions with the thermal plasma. Astron. Astrophys. 2022, 659, A146. [Google Scholar]
Figure 1. The value of d in the last lines of Equations (2) and (3) as a function of Z .
Figure 1. The value of d in the last lines of Equations (2) and (3) as a function of Z .
Applsci 13 09643 g001
Figure 2. The value of f in the last lines of Equations (2) and (3) as a function of Z .
Figure 2. The value of f in the last lines of Equations (2) and (3) as a function of Z .
Applsci 13 09643 g002
Figure 3. The principal Hugoniot of Pm: new analytic model vs. theoretical calculations using the relativistic Green’s function quantum average atom code Tartarus [31,32].
Figure 3. The principal Hugoniot of Pm: new analytic model vs. theoretical calculations using the relativistic Green’s function quantum average atom code Tartarus [31,32].
Applsci 13 09643 g003
Figure 4. The principal Hugoniot of Tb: our new analytic model vs. the experimental data of Ref. [8] (Carter et al.), the Tb Hugoniot from TEFIS database (TEFIS), and our theoretical calculations using the Thomas–Fermi model with corrections (Thomas–Fermi Corr.).
Figure 4. The principal Hugoniot of Tb: our new analytic model vs. the experimental data of Ref. [8] (Carter et al.), the Tb Hugoniot from TEFIS database (TEFIS), and our theoretical calculations using the Thomas–Fermi model with corrections (Thomas–Fermi Corr.).
Applsci 13 09643 g004
Figure 5. The principal Hugoniot of Tm: our new analytic model vs. the experimental data of Ref. [8] (Carter et al.), the Tm Hugoniot from TEFIS database (TEFIS), and our theoretical calculations using the Thomas–Fermi model with corrections (Thomas–Fermi Corr.).
Figure 5. The principal Hugoniot of Tm: our new analytic model vs. the experimental data of Ref. [8] (Carter et al.), the Tm Hugoniot from TEFIS database (TEFIS), and our theoretical calculations using the Thomas–Fermi model with corrections (Thomas–Fermi Corr.).
Applsci 13 09643 g005
Figure 6. The principal Hugoniot of Lu: our new analytic model vs. the experimental data of Ref. [7] (Al’tshuler et al.) and our theoretical calculations using Tartarus.
Figure 6. The principal Hugoniot of Lu: our new analytic model vs. the experimental data of Ref. [7] (Al’tshuler et al.) and our theoretical calculations using Tartarus.
Applsci 13 09643 g006
Figure 7. The principal Hugoniot of Ac: our new analytic model vs. the Ac Hugoniot from TEFIS database (TEFIS) and our theoretical calculations using the Thomas–Fermi model with corrections (Thomas–Fermi Corr.).
Figure 7. The principal Hugoniot of Ac: our new analytic model vs. the Ac Hugoniot from TEFIS database (TEFIS) and our theoretical calculations using the Thomas–Fermi model with corrections (Thomas–Fermi Corr.).
Applsci 13 09643 g007
Figure 8. Internal energy along the principal Hugoniot for Ac: the new analytic model (curve) vs. our own theoretical calculations using the Thomas–Fermi model with corrections (symbols).
Figure 8. Internal energy along the principal Hugoniot for Ac: the new analytic model (curve) vs. our own theoretical calculations using the Thomas–Fermi model with corrections (symbols).
Applsci 13 09643 g008
Table 1. Numerical values of the parameters for the analytic model of the principal Hugoniot for Pm, Tb, Tm, Lu, and Ac.
Table 1. Numerical values of the parameters for the analytic model of the principal Hugoniot for Pm, Tb, Tm, Lu, and Ac.
Z ρ 0 CB A × 10 2 cb a × 10 5 Ld f × 10 3 U p * U p max U p * *
Pm617.251.681.34−0.4018593.488431.167987.2095811680316.8604−12.014521.0258219.9681532.46
Tb658.231.721.29−0.1980523.533541.168046.9927812638336.9251−17.092929.7399224.7921574.89
Tm699.321.901.19−6.83647 ×   10 5 3.638141.168346.6795913603820.3201−12.2217160.493233.3811649.44
Lu719.840.981.65−2.158973.658391.168336.5620614089019.4395−11.746311.1213236.1161679.54
Ac8910.051.5351.225−0.1094462.496981.158306.171841850756.05256−5.7186828.8450201.1411913.82
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Burakovsky, L.; Preston, D.L.; Ramsey, S.D.; Sjue, S.K.; Starrett, C.E.; Baty, R.S. Principal Hugoniots of Promethium, Terbium, Thulium, Lutetium, and Actinium in a Wide Pressure Range. Appl. Sci. 2023, 13, 9643. https://doi.org/10.3390/app13179643

AMA Style

Burakovsky L, Preston DL, Ramsey SD, Sjue SK, Starrett CE, Baty RS. Principal Hugoniots of Promethium, Terbium, Thulium, Lutetium, and Actinium in a Wide Pressure Range. Applied Sciences. 2023; 13(17):9643. https://doi.org/10.3390/app13179643

Chicago/Turabian Style

Burakovsky, Leonid, Dean L. Preston, Scott D. Ramsey, Sky K. Sjue, Charles E. Starrett, and Roy S. Baty. 2023. "Principal Hugoniots of Promethium, Terbium, Thulium, Lutetium, and Actinium in a Wide Pressure Range" Applied Sciences 13, no. 17: 9643. https://doi.org/10.3390/app13179643

APA Style

Burakovsky, L., Preston, D. L., Ramsey, S. D., Sjue, S. K., Starrett, C. E., & Baty, R. S. (2023). Principal Hugoniots of Promethium, Terbium, Thulium, Lutetium, and Actinium in a Wide Pressure Range. Applied Sciences, 13(17), 9643. https://doi.org/10.3390/app13179643

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop