_{1}

The equation of state and the shear modulus data of sintered Nd <SUB>2</SUB>Fe <SUB>14</SUB>B were investigated up to 140 GPa by the Gruneisen’s model, the volume superposition principle and the Hugoniot’s relations. Then, the results were compared to the prior experiments with a standard deviation of 0.125% from 18 GPa to 78 GP; and then, the loading pressure was extended to higher. Meanwhile, the softening feature has not been observed both in adiabat and shear modulus throughout the interested range.

Creating a simple, usable and reliable computing model of a complex mixture, such as NdFeB, is quite stubborn theoretically under high pressure and high temperature, but is really significant in aerospace engineering and defense industry [

For the sake of simplifying our calculation, the anisotropic Nd_{2}Fe_{14}B was chosen, because of its principle application in modern society with its high magnetic energy product, remanence ratios and coercive force, and we ignored the existence of B in the model.

Based on the shock experimental inertia of Nd_{2}Fe_{14}B within 18 GPa - 78 GPa [

where P, V, E, g, P_{x} and E_{x} are, respectively, total pressure, volume, energy, Gruneisen’s coefficient, cold pressure and energy. Among of them, P_{x} and E_{x} can be specifically represented by the Born-Meyer’s potential [

where Q, q, r_{0} and a_{v} are, respectively, the parameters of cold energy, the initial density and the volume dilatation, d = r/r_{0k} and r_{0k} = r_{0} (1 + 300 a_{n}).

Instead of computing this intermetallic compound, each component, Nd, Fe and B will be allowed to be treated separately, according to their proportion by weight, 26.68%, 72.32% and 0.99% due to the chemical stability of the compound under high pressure and high temperature, as mentioned above. The average formulae are listed below,

where a_{i} denotes weight percentage of each component.

The interaction of particles in Nd_{2}Fe_{14}B will be ignored if the compound is actually treated as an ideal mixture. The compressibility of Nd_{2}Fe_{14}B, therefore, can be evaluated by the Hugoniot of each component itself [

where c_{0} and l_{0} are, respectively, the wave velocity and the characteristic parameter of Nd_{2}Fe_{14}B. At the initial state, they are,

where m_{i} stands for the mass of each independent component in mixture. Sequentially, we take c_{0} = 3.686 and λ_{0} = 1.059, at the same pressures in literature [

It must be emphasized that c_{0} and λ_{0} in Equations (8)-(11) should be demarcated to the absolute zero from the room temperature in calculation by the linear mixed rule [

where g_{0} is the initial Grunasen’s coefficient of Nd_{2}Fe_{14}B, which is taken as 1.492, and ρ_{0k} as 7.451 g/cm^{3} [

The effective shear modulus, G, can be deduced from the SCG model [

where c_{l} and c_{b} are, respectively, the Euler’s longitudinal and bulk wave velocity, and given by,

where n, the Poisson’s ratio, is taken as 0.24, and g is supposed to depend on volume V only. Then we have,

where r is density at final state, and a = 0.822 nm is the lattice constant of Nd_{2}Fe_{14}B.

For convenience reasons in the following calculation and discussion, part of initial parameters of Nd and Fe are listed in

Compare with prior experiments, the calculated results have been listed in

By using Equations (1)-(16), a new D-u relation, D = 3.476 + 1.203 u, of Nd_{2}Fe_{14}B can be obtained, as shown in _{2}Fe_{14}B up to 140 GPa.

Considering the accordance between calculations and experiments below 78 GPa in

Elements | C_{0}/km/S | l_{0} | g_{0} | a_{n}/10^{−}^{5}/K | a_{i}/% | r_{0}/g/cm^{3} | ||
---|---|---|---|---|---|---|---|---|

Nd | 2.2^{ } | 1.015^{ } | 0.57 | 2.994 | 26.68 | 7.003 | 2.219 | 1.010 |

Fe | 3.955^{ } | 1.30^{ } | 1.78 | 3.51 | 72.32 | 7.856 | 3.989 | 1.291 |

^{*}p/GPa | p/GPa | ^{*}D/km/s | D/km/s | ^{*}u/km/s | u/km/s | ^{*}r/g/cm^{3} | r/g/cm^{3} |
---|---|---|---|---|---|---|---|

77.84 | 78 | 5.695 | 5.64 | 1.861 | 1.80 | 11.034 | 10.944 |

50.99 | 53.32 | 5.040 | 5.18 | 1.377 | 1.40 | 10.222 | 10.211 |

41.10 | 39.58 | 4.929 | 4.68 | 1.135 | 1.00 | 9.653 | 9.476 |

27.75 | 26.4 | 4.609 | 4.44 | 0.820 | 0.80 | 9.036 | 9.089 |

19.23 | 18.23 | 4.302 | 4.20 | 0.608 | 0.60 | 8.653 | 8.693 |

^{*}Cited from literature [

the segments are also reasonable from 78 GPa to 140 GPa based on the prior assumption of the stability chemically as loading pressure on Nd_{2}Fe_{14}B.

On the basis of above discussion, the longitudinal, bulk wave velocities and the shear modulus of Nd_{2}Fe_{14}B can be estimated by Equations (14)-(18) as shown together in _{2}Fe_{14}B were unchanged with such rising pressure, which meant no shock softening meanwhile, or no melting occurred from the lowest to the highest pressure region once again.

A simple and effective algorithm for evaluating Nd_{2}Fe_{14}B has been developed via comparing experiments below 78 GPa and predicted towards 140 GPa together with the shear modulus, even no corresponded experimental data accompanied. We believe the deviation in calculation can be further reduced by carefully adjusting parameters in the equations listed above; and the processing method can be consequentially extended to other inert mixtures under high pressure and high temperature. But two rules should be followed:

1) Each component of solid mixture should be inert during shock compression, and

2) Each initial parameter of components should be known in advance.