A Simple Derivation of the Birch–Murnaghan Equations of State (EOSs) and Comparison with EOSs Derived from Other Deﬁnitions of Finite Strain

: Eulerian ﬁnite strain of an elastically isotropic body is deﬁned using the expansion of squared length and the post-compression state as reference. The key to deriving second-, third- and fourth-order Birch–Murnaghan equations-of-state (EOSs) is not requiring a di ﬀ erential to describe the dimensions of a body owing to isotropic, uniform, and ﬁnite change in length and, therefore, volume. Truncation of higher orders of ﬁnite strain to express the Helmholtz free energy is not equal to ignoring higher-order pressure derivatives of the bulk modulus as zero. To better understand the Eulerian scheme, ﬁnite strain is deﬁned by taking the pre-compressed state as the reference and EOSs are derived in both the Lagrangian and Eulerian schemes. In the Lagrangian scheme, pressure increases less signiﬁcantly upon compression than the Eulerian scheme. Di ﬀ erent Eulerian strains are deﬁned by expansion of linear and cubed length and the ﬁrst- and third-power Eulerian EOSs are derived in these schemes. Fitting analysis of pressure-scale-free data using these equations indicates that the Lagrangian scheme is inappropriate to describe P-V-T relations of MgO, whereas three Eulerian EOSs including the Birch–Murnaghan EOS have equivalent signiﬁcance. second-power Lagrangian, ﬁrst-power Eulerian, and third-power Eulerian ﬁnite strains.


Introduction
Density distributions are one of the most fundamental properties to describe planetary interiors. Materials within planetary interiors are under high pressure and, therefore, have higher densities than under ambient conditions. Material density measurements using high-pressure experimental techniques have a limited pressure range and obtained pressure-density data must often be extrapolated to higher pressures. For this procedure, it is necessary to fit available experimental data to a certain formula, which is referred to as the "equation of state" (EOS).
Various EOS formulas have been proposed. Among them, the third-order Birch-Murnaghan EOS is the most frequently used. This EOS is constructed on the basis of finite elastic strain theory in the Eulerian scheme. When strains are infinitesimal, the dimensions of a body decrease linearly with pressure. However, this situation is not the case when strains are finite because matter becomes increasingly incompressible with pressure. A theory to treat this phenomenon is, therefore, required: the finite elastic strain theory. The Eulerian scheme describes compression using the post-compression state as reference. In contrast, the other scheme, namely, the Lagrangian scheme describes compression using the pre-compressed state as reference.

Eulerian Finite Strain
The foundation of the Birch-Murnaghan EOS is Eulerian finite strain, and we therefore begin with its introduction. An important aspect of the Eulerian scheme is that the reference is defined by the post-compression state, whereas the Lagrangian scheme uses the pre-compression state. A second important point is that changes are expanded in squared length before and after compression.
Let us consider a cubic body with edge length X 0 . Its volume is accordingly V 0 = X 3 0 . This cube is uniformly compressed to edge length X and accordingly to a volume of V = X 3 ( Figure 1). The edge length before compression is expressed by the reference length and length change, or displacement u, as: Minerals 2019, 9, x FOR PEER REVIEW 2 of 18 Historically, finite elastic strain in the Euler scheme was treated by Murnaghan [1] who introduced Eulerian finite strain (see Equation (9)) and expressed pressure as a quadratic function of Eulerian finite strain. Birch [2] extended the theory of Murnaghan [1] to derive a prototype of the EOS, which is referred to as the Birch-Murnaghan EOS. They proposed an EOS using tensors in three dimensions, but other simpler arguments are available.
Poirier [3] plainly derived the Eulerian finite strain and Birch-Murnaghan EOS starting from the infinitesimal length squared with tensors in the three-dimensional space. However, if we ignore elastic anisotropy for matter compressed uniformly, there is no need to use a tensor or any special reason to start from infinitesimal length squared. Anderson [4] also derived the Eulerian finite strain and Birch-Murnaghan EOS. His derivation is, however, not easy to follow, either.
This article presents a more simple derivation of Eulerian finite strain and the Birch-Murnaghan EOS so that its essence can be easily understood. For comparison, other EOSs are constructed based on the Lagrangian scheme and other definitions of finite strain, and the usefulness and applicability of the Birch-Murnaghan EOS is discussed. Because we treat isothermal EOS's in this article, all differentiation is carried out at constant temperature. A partial derivative of a quantity with respect to some parameter at constant temperature is, therefore, simply referred to as a derivative of .

Eulerian Finite Strain
The foundation of the Birch-Murnaghan EOS is Eulerian finite strain, and we therefore begin with its introduction. An important aspect of the Eulerian scheme is that the reference is defined by the post-compression state, whereas the Lagrangian scheme uses the pre-compression state. A second important point is that changes are expanded in squared length before and after compression.
Let us consider a cubic body with edge length . Its volume is accordingly = . This cube is uniformly compressed to edge length and accordingly to a volume of = ( Figure 1). The edge length before compression is expressed by the reference length and length change, or displacement , as: Because of compression, < 0.  The initial edge length of the cube X 0 is reduced to the final edge length X = X 0 + u. The volume of the cube accordingly decreases from V 0 = X 3 0 to V = X 3 .
Because of compression, u < 0. As mentioned, the change in squared edge length by compression is expanded as: (2) Displacement should be proportional to the reference edge length because of uniform compression and is expressed using a proportional constant at a given compression, c, as: The quantity c is referred to as the strain in the linear elasticity. Using Equation (3), Equation (2) becomes: Finite strain in the Eulerian scheme, ε E2 , is defined as: In this article, finite strain is referred to as the second-power Eulerian strain because it is defined by expanding the second power of length. The change of squared length (Equation (2)) is expressed using the Eulerian finite strain (Equation (5)) as: The ratio of the edge length before compression to the reference edge length is expressed by Eulerian finite strain as: The ratio of the pre-compression volume to the reference volume is also expressed by Eulerian finite strain as: Eulerian finite strain is therefore expressed by the volume ratio as: The ε is negative when a body is compressed. To make strain positive by compression, we introduce f E2 as: The quantity f E2 instead of ε E2 is usually referred to as Eulerian finite strain.

The Second-Order Birch-Murnaghan EOS
In isothermal EOSs, pressure, P, is expressed as a function of volume, V. From thermodynamics, pressure is the volume derivative of Helmholtz energy, F, as: The Helmholtz free energy of matter should increase with compression and may be expressed by a series of the Eulerian finite strain: Because the absolute value of F is arbitrary, the coefficient of the first term in Equation (12) can be a 0 = 0. Because pressure should be zero in an uncompressed state, f E2 = 0, we have: The coefficient of the second term in Equation (12) is therefore a 1 = 0.
Here, we truncate Equation (12) to the f E2 2 term for the first approximation and therefore have: By substituting Equation (14) into Equation (11), we have From the definition of finite strain (Equation (10)), the volume derivative of Eulerian finite strain is: As shown in Appendix A, the coefficient a 2 is given by: where K T 0 is the isothermal bulk modulus at standard temperature. By substituting Equations (10), (16), and (17) into Equation (15), we have the second-order Birch-Murnaghan EOS: This equation contains a subtraction formula between the 7/3 and 5/3 powers of V 0 /V. The difference of these two powers, 7/3 − 5/3 = 2/3 is owing to the definition of Eulerian finite strain (Equation (10)), and the power of 5/3 in the second term is because of the volume derivative of Eulerian finite strain (Equation (16)).

The Third-Order Birch-Murnaghan EOS
The concept of the third-order Birch-Murnaghan EOS is almost identical to that of the second-order equation. The difference is that the Helmholtz free energy expressed by the Eulerian finite strain (Equation (12)) is truncated not up to the second term but to the third term as: By substituting Equation (19) into Equation (11) and differentiating it with respect to volume, we have: where ξ 1 = 3a 3 /2a 2 . As presented in the Appendix A, the parameter ξ 1 is given by: where K T 0 is the pressure derivative of the isothermal bulk modulus at standard temperature. By substituting the definition of Eulerian finite strain (Equation (10)), its volume derivative (Equation (16)), the parameter a 2 (Equation (17)), and parameter m 1 (Equation (21)) into Equation (20), we have the third-order Birch-Murnaghan EOS as: The second term in the first square bracket appears because of the truncation of the Helmholtz free energy to the higher-order (third) term. The form of the curly bracket is owing to the Eulerian finite strain (Equation (10)).
The third-order equation (Equation (22)) becomes identical to the second-order equation (Equation (18)) when On the other hand, if K T 0 is neglected as K T 0 = 0, the third-order equation should differ from the second-order equation.

The Fourth-Order Birch-Murnaghan EOS
For the fourth-order EOS, Equation (12) is truncated up to the fourth term as: The pressure is then expressed as: where ξ 2 = 2a 4 /a 2 . As given in the Appendix A, the parameter ξ 2 is: where K T 0 is the second pressure derivative of the isothermal bulk modulus at standard temperature. Similar to the lower-order EOSs, we have the fourth-order Birch-Murnaghan EOS as: This fourth-order equation (Equation (27)) becomes identical to the third-order equation when Again, if K T 0 is neglected as K T 0 = 0, the fourth-order equation should differ from the third-order equation.

Equations of States from Other Finite Strain Definitions
Through deriving the Birch-Murnaghan EOS, we may have some questions in mind. One is why the Eulerian scheme is more frequently used to describe the compression of bodies instead of the Lagrangian scheme. Another question is why is a change in squared length considered; why not other powers of length? Murnaghan's [1] argument leaves one to reflect that the use of squared length was based on Pythagorean theorem but there is no physical reason. To answer these questions, we derive EOSs from finite strain using the following different definitions:

1.
Lagrangian scheme, expansion of squared length, referred to as the second-power Lagrangian EOS; 2.
Eulerian scheme, linear expansion of length, referred to as the first-power Eulerian EOS; 3.
Eulerian scheme, expansion of cubed length, referred to as the third-power Eulerian EOS.

The Second-Power Lagrangian EOS
Let us come back to Figure 1. In the Lagrangian scheme, the reference is taken as the pre-compressed state. Hence, the length after compression, X, is expressed using the reference length, X 0 , and displacement: The change in squared edge length is: The displacement is expressed by the reference length as: Note that Equation (31) differs from Equation (3) unless X X 0 , namely when strain is finite. Using Equation (31), Equation (30) becomes: We now define the Lagrangian finite strain, ε L2 , as: Like the Eulerian scheme, Lagrangian finite strain is expressed using the volume ratio: For convenience, the quantity f L2 below is hereafter referred to as the second-power Lagrangian finite strain: EOSs are constructed based on this strain, which is hereafter referred to as the second-power Lagrangian EOS. The Birch-Murnaghan EOS is hereafter referred to as the second-power Eulerian EOS.
Like the second-order second-power Eulerian EOS, the Helmholtz free energy is expressed with the second-power Lagrangian finite strain squared as: The expression for pressure in the Lagrangian scheme is identical to the Eulerian scheme as: The volume derivative of the second-power Lagrangian finite strain is given by: The formula of parameter a 2 in the Lagrangian scheme is identical to that in the Eulerian scheme (Equation (17)). As a result, the second-order second-power Lagrangian EOS is given by: For the third-order EOS, the Helmholtz free energy is expanded to the third term of f L2 as: By differentiating Equation (33) with respect to volume, pressure is expressed as: The expression of m 1 for the second-power Lagrangian finite strain differs from that for the second-power Eulerian strain and is given by: From Equations (17), (28), (34), and (35), we have the third-order second-power Lagrangian EOS as: The third-order equation (Equation (43)) becomes identical to the second-order equation (Equation (39)) when: In contrast to the Eulerian scheme, the neglect of K T 0 as K T 0 = 0 is identical to the neglect of a term higher than third order in the Helmholtz free energy (Equation (36)).

Finite Strains and EOSs from Linear Length
The change of length itself, namely the first power of length, is considered here in the Eulerian scheme and given from Equations (1) and (3) as: The first-power Eulerian strain, f E1 , is therefore defined as: Similar to the change in squared length, the first-power Eulerian strain is expressed by the volume change V 0 /V as: The volume derivative of the first-power Eulerian finite strain is given by: The parameter a 2 for the first-power Eulerian strain is identical to those of the second-power strain in the Eulerian and Lagrangian schemes. The second-order first-power Eulerian EOS is, therefore, given by: The parameter ξ 1 is given by: The third-order first-power Eulerian EOS is therefore given by:

Finite Strains and EOSs from Cubed Length
We now consider the change in length cubed, namely the third-power of length in the Eulerian scheme, which is given from Equations (1) and (3) as: The third-power Eulerian strain, f E3 , is therefore defined as: Similar to those with the change in squared length, the first-power Eulerian strain is expressed by volume change V 0 /V: The volume derivative of the third-power Eulerian finite strain is given by: In the same way as shown previously, the second-order third-power Eulerian EOS is given by: The parameter ξ 1 is given by: The third-order third-power Eulerian EOS is therefore given by: 4. Discussion Figure 2 shows pressures obtained by the second-, third-and fourth-order Birch-Murnaghan EOSs of NaCl in the B1-structure, Au, and MgO, which are materials frequently used as pressure standards in high-pressure experiments. The isothermal bulk moduli of these three materials and their first and second pressure derivatives at ambient temperature and zero pressure used for construction of Figure 2 are summarized in Table 1. Note that these EOS parameters were all obtained by the latest studies of sound velocity measurements [5][6][7]. Table 1. Bulk moduli and their pressure derivatives of frequently used pressure standard materials.

Material
K T 0 (GPa) The second-and third-order Birch-Murnaghan EOSs give significantly different pressures for NaCl and Au under high compression owing to large deviations of K T 0 from 4 (5.14 and 6.56). At V/V 0 = 0.75, the second-and third-order EOSs yield 12.1 and 14.3 GPa for NaCl and 85 and 112 GPa for Au, respectively, with differences of 18% and 41% in these cases. Birch [8] justified the validity of the Birch-Murnaghan EOS by identifying the second-and third-order EOSs within experimental uncertainty (e.g., K T 0 4). However, modern experimental results have demonstrated that K T 0 K S 0 and significantly differs from 4 in many materials (see Table 11 in Bass [9]), which therefore compromises the validity of the second-order Birch-Murnaghan EOS. Figure 2a also shows that the fourth-order equation gives slightly lower pressures than the third-order equation in the case of NaCl. At V/V 0 = 0.75, the fourth-order equation gives a pressure of 13.7 GPa, which is about 4% lower than that from the third-order equation. For the third-and fourth-order equations to give identical pressures, the second pressure derivative of the isothermal bulk modulus must not be K T 0 = −0.392 GPa −1 but rather K T 0 = −0.267 GPa −1 . Figure 2b,c also show the fourth-order Birch-Murnaghan EOSs of Au and MgO whose second pressure derivatives of the isothermal bulk modulus at zero pressure are set to K T 0 = 0. The pressure at the bottom of the Earth's mantle is 136 GPa. To obtain such pressure by the third-order Birch-Murnaghan EOSs of Au and MgO, the compression must be V/V 0 = 0.736 and V/V 0 = 0.678, respectively. Under these compressions, the fourth-order Birch-Murnaghan EOSs with K T 0 = 0 give 160 and 152 GPa. Thus, approximations that ignore higher-order terms of finite strain and equate higher-order derivatives of the isothermal bulk modulus to zero are not identical in the Birch-Murnaghan EOSs.
The differences in pressure given by the two kinds of EOSs above correspond to 470 and 330 km in the case of Au and MgO, respectively. The depth of the mantle phase transition at the base of the Earth's mantle is particularly complicated to argue. Although Hirose et al. [10] determined the post-spinel transition boundary in MgSiO 3 using Au and MgO pressure scales to interpret the phase boundary pressure as the D" layer, their argument contains errors on the order of a few hundred km simply by neglecting K T 0 as K T 0 = 0.

Eulerian Versus Lagrangian Schemes
An essential difference between the Eulerian and Lagrangian finite strains is that the volume ratio is inversed with respect to one another (Equations (10) and (35)). Because the reference state is before compression in the Lagrangian scheme, the post-compression volume is put in the numerator of the strain formula (Equation (35)). As shown in Figure 3, the Lagrangian finite strain increases nearly linearly upon compression to V/V 0 = 0.2 and then approaches a finite value (0.5) with further compression ( V/V 0 → 0 ). As a result, the Helmholtz free energy and its volume derivative, namely pressure, only moderately increase with compression. However, it is important to note that pressure approaches infinity with compression to zero volume because the volume derivative of the finite strain and one term of the EOS are proportional to the −1/3 power of V/V 0 (Equations (38), (39), and (43)).  In contrast to the Lagrangian scheme, the Eulerian reference state is after compression and the post-compression volume is put in the denominator of the strain formula (Equation (10)). As shown in Figure 3, all Eulerian finite strains rapidly increase and then diverge to infinity as ⁄ → 0. As a result, the Helmholtz free energy and pressure also rapidly increase to infinity with ⁄ → 0. This behavior in the Eulerian scheme is in better agreement than that in the Lagrangian scheme. It is also important to note that the rate of pressure increase is very similar in both schemes under low compression. In other words, both schemes give the same results if the strain is infinitesimal. Figure 2a-c show pressures given by the second-and third-order (second-power) Lagrangian EOSs as well as the Eulerian EOS's. For comparison, we also show pressures given by the simplest EOS, which is obtained by integrating the definition of the isothermal bulk modulus as: with green curves in Figure 2a-c. As previously discussed, finite strains given in the Eulerian scheme increase much more rapidly with compression than those of the Lagrangian scheme. The second-order second-power Lagrangian In contrast to the Lagrangian scheme, the Eulerian reference state is after compression and the post-compression volume is put in the denominator of the strain formula (Equation (10)). As shown in Figure 3, all Eulerian finite strains rapidly increase and then diverge to infinity as V/V 0 → 0 . As a result, the Helmholtz free energy and pressure also rapidly increase to infinity with V/V 0 → 0 . This behavior in the Eulerian scheme is in better agreement than that in the Lagrangian scheme. It is also important to note that the rate of pressure increase is very similar in both schemes under low compression. In other words, both schemes give the same results if the strain is infinitesimal. Figure 2a-c show pressures given by the second-and third-order (second-power) Lagrangian EOSs as well as the Eulerian EOS's. For comparison, we also show pressures given by the simplest EOS, which is obtained by integrating the definition of the isothermal bulk modulus as: with green curves in Figure 2a-c. As previously discussed, finite strains given in the Eulerian scheme increase much more rapidly with compression than those of the Lagrangian scheme. The second-order second-power Lagrangian EOS gives almost identical pressures as the integration of a constant bulk modulus (Equation (37)). This means that the second-order Lagrangian EOS gives pressures without considering an increase in bulk modulus with compression. This is reasonable because the third-order Lagrangian EOS becomes identical to the second-order Lagrangian EOS when K T 0 is neglected as K T 0 = 0. In contrast, the Eulerian EOSs give substantially higher pressures than when integrating a constant bulk modulus (Equation (59)). As can be seen from the fact that the second-and third-order (second-power) Eulerian EOSs become identical in the case of K T 0 = 4, the Eulerian EOSs implicitly contain the effects of an increased bulk modulus with compression.
The EOSs in both schemes converge with increasing order owing to better approximations of the Helmholtz free energy by the finite strains. This is the case for MgO (Figure 2c), however not for Au (Figure 2b). The small difference between the third-and fourth-orders of EOSs of NaCl (Figure 2a) does not imply that expansion to higher orders effectively helps conversion of the two schemes. These observations imply an essential problem with the Lagrangian scheme. Because the finite strains based on higher powers of length increase more rapidly with compression, the EOSs based on higher powers of length are expected to give higher pressures. This is the case for materials with high K" such as Au. As is seen in Figure 4b, the first-, second-and third-power EOSs of Au gives pressures of 73, 76, and 78 GPa at V/V 0 = 0.8, respectively. The differences of the first-and third-power EOSs from the Birch-Murnaghan (second-power) EOS are both about −4% and +3%. The discrepancy between different powers of EOSs increases with compression, giving pressures of 170, 184, and 194 GPa, respectively, with differences of about −8% and +5%. The finite strains more rapidly increase under compression with increasing power of V 0 /V from 1/3 to 2/3 and then to 1 (Equations (10), (47), and (54)), as shown in Figure 3. Because these three EOSs are obtained in the same way except for the power of length to define the finite strain, we cannot say that the Birch-Murnaghan EOS provides substantially more accurate pressures than the others.

Equations of State Obtained from Expansions of Different Powers of Length
As discussed previously, Birch [8] justified the Birch-Murnaghan EOSs by the identity of the second-and third-order EOSs when K T 0 4, which is approximately the case in many kinds of materials. As Equations (51) and (58) show, the second-and third-order Eulerian EOSs based on finite strains defined from the first-and third-powers of length become identical when K T 0 = 3 and K T 0 = 5, respectively. Again, as summarized by Bass [9], the K T 0 of the majority of solids are larger than 4. In the case of Au, K T 0 = 6.56 and the third-power EOS (Equation (58)) should therefore be more appropriate for Au than the Birch-Murnaghan EOS (Equation (22)) if we follow Birch's [8] justification. powers of EOSs increases with compression, giving pressures of 170, 184, and 194 GPa, respectively, with differences of about −8% and +5%. The finite strains more rapidly increase under compression with increasing power of ⁄ from 1/3 to 2/3 and then to 1 (Equations (10), (47), and (54)), as shown in Figure 3. Because these three EOSs are obtained in the same way except for the power of length to define the finite strain, we cannot say that the Birch-Murnaghan EOS provides substantially more accurate pressures than the others.   As discussed previously, Birch [8] justified the Birch-Murnaghan EOSs by the identity of the second-and third-order EOSs when ′ ≅ 4, which is approximately the case in many kinds of materials. As Equations (51) and (58) show, the second-and third-order Eulerian EOSs based on finite strains defined from the first-and third-powers of length become identical when = 3 and = 5, respectively. Again, as summarized by Bass [9], the of the majority of solids are larger than 4.

Examination of Equations of State Using Pressure Scale-Free Experimental Data
In this section, we examine the validity of the second-power Lagrangian and first-and third-power Eulerian third-order EOS's in comparison with the third-order Birch-Murnaghan (second-power Eulerian) EOS. Specifically, we attempt to obtain parameters of the Mie-Grüneisen-Debye thermal EOS of MgO on the basis of these isothermal EOSs using pressure scale-free experimental data following the method by Tange et al. [11]. "Pressure scale-free" data were data that can be used to build an equation of state, but no other pressure scales such as pressure values obtained using equations of state of other pressure standard materials were used to obtain those data. Use of pressure scale-free data allows us building an equation of state with avoiding any circular argument. The pressure scale-free experimental data sets used here are zero-pressure and high-temperature thermal expansivity, zero-pressure and high-temperature adiabatic bulk moduli, 300-K and compressed adiabatic bulk moduli, and shock compression. Data sources are summarized in Table 2 of Tange et al. [11]. The following parameters at 300 K and 0 GPa are fixed: lattice volume, adiabatic bulk modulus, thermal expansivity, and isobaric heat capacity. The parameters obtained by these fittings are given in Table 3 of Tange et al. [11]. Note that parameters a and b are adjustable to express the volume dependence of the Grüneisen parameter γ, according to the following formula: Table 2 lists the obtained parameters of the thermal EOSs and Figure 5 shows the reproduction of Hugoniot curves using the obtained parameters. V 0 : lattice volume at 300 K and 0 GPa; K T 0 : isothermal bulk modulus; K T 0 : pressure derivative of the bulk modulus; θ 0 : Debye temperature; γ 0 : Grüneisen parameter; a, b: parameters to express volume dependence of the Grüneisen parameter given in Equation (59). Figure 5 indicates that all Eulerian EOSs provide Hugoniot curves in agreement with the experimental data by Marsh [12] and Duffy and Ahrens [13], whereas the Lagrangian EOS does not. The reason for the failure in construction of an EOS in the Lagrangian scheme is that the parameters a and b in Equation (60) must not be less than zero; namely, thermal pressure or thermal expansivity should not increase with increasing pressure. Because the Lagrangian EOS gives much lower pressures at ambient temperature than the Eulerian EOSs, the thermal pressure must be abnormally high in the Lagrangian scheme to reproduce the shock experiment data. As a conclusion, the Eulerian scheme is more appropriate than the Lagrangian scheme to construct an EOS of real materials. On the other hand, linear and cubed length expansions in the Eulerian scheme provide equivalently appropriate EOSs compared with the EOS derived from the squared length expansion to define finite strain. : pressure derivative of the bulk modulus; : Debye temperature; : Grüneisen parameter; a, b: parameters to express volume dependence of the Grüneisen parameter given in Equation (59).

Figure 5.
Hugoniot curves of MgO reproduced by the various EOSs using pressure scale-free data following Tange et al. [11]. The dashed light-red, thin-solid red, dotted dark-red, and solid blue curves are obtained from the first-, second-and fourth-Eulerian and second-power Lagrangian EOSs, respectively. The light-gray square and dark-gray circle are experimental data obtained by Marsh [12] and Duffy and Ahrens [13], respectively.

Figure 5.
Hugoniot curves of MgO reproduced by the various EOSs using pressure scale-free data following Tange et al. [11]. The dashed light-red, thin-solid red, dotted dark-red, and solid blue curves are obtained from the first-, second-and fourth-Eulerian and second-power Lagrangian EOSs, respectively. The light-gray square and dark-gray circle are experimental data obtained by Marsh [12] and Duffy and Ahrens [13], respectively. modulus at zero pressure. The other assumes that the change in Helmholtz free energy by compression is proportional to the finite strain squared. Let us first present the expression from the isothermal bulk modulus, which is defined as: The partial derivative of pressure with respect to volume at zero pressure is therefore: Second, let us present the expression from the Helmholtz free energy by the Eulerian finite strain. By differentiating Equation (15) with respect to volume at constant temperature, we have: At zero pressure, the strain is zero: f = 0. The partial derivative of pressure with respect to volume at constant temperature and zero pressure is therefore: For Equations (16), (38), (48), and (57), the volume derivative of finite strain at zero pressure is identical for all Eulerian and Lagrangian finite strains as: The essence of this derivation is identical to that of a 2 . We present the second volume derivative of pressure, ∂ 2 P/∂V 2 T , in this proof. We first express the volume second derivative of pressure using the isothermal bulk modulus and its pressure derivative. From the definition of the isothermal bulk modulus, we have At zero pressure, Equation (A7) is: