Next Article in Journal
Longitudinal Near-Field Coupling between Acoustic Resonators Grafted onto a Waveguide
Previous Article in Journal
FOX, Current State and Possibilities
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

First-Principles Study of the Nonlinear Elasticity of Rare-Earth Hexaborides REB6 (RE = La, Ce)

1
School of Mathematical Sciences and Physics, Jinggangshan University, Ji’an 343009, China
2
Research Center of Laser Fusion, China Academy of Engineering Physics, Mianyang 621900, China
3
School of Materials Science and Engineering, Shanghai Jiaotong University, Shanghai 200240, China
4
School of Materials Science and Engineering, Jiangxi University of Science and Technology, Ganzhou 341000, China
*
Author to whom correspondence should be addressed.
Crystals 2017, 7(11), 320; https://doi.org/10.3390/cryst7110320
Submission received: 30 March 2017 / Revised: 16 October 2017 / Accepted: 20 October 2017 / Published: 25 October 2017

Abstract

:
The complete set of independent second- and third-order elastic constants of rare-earth hexaborides LaB 6 and CeB 6 are determined by the combination method of first-principles calculations and homogeneous deformation theory. The ground-state lattice parameters, second-order elastic constants, and bulk modulus are in reasonable agreement with the available experimental data. The third-order elastic constant of longitudinal mode C 111 has a larger absolute value than other shear modes, showing the contribution to lattice vibrations from longitudinal modes to be greater. The pressure derivatives of the second-order elastic constants related to the third-order elastic constants are calculated to be positive for the two hexaborides, which are consistent with those of their polycrystalline bulk modulus and shear modulus. Furthermore, the effect of pressure on the structural stability, mechanical property, and elastic anisotropy of the two hexaborides are investigated, showing a reduction in mechanical stability and an increase in ductility and anisotropy with increasing pressure.

1. Introduction

Rare-earth hexaborides ( R E B 6 ) have been widely used in various field-electron emitter devices and high-energy optical systems due to their attractive properties of high melting point, high mechanical strength, low work function, low volatility at high temperatures, conductivity, chemical resistance, brightness, small optical size, long service life, and the monoenergetic character of their electrons [1]. These attractive properties are closely related to their crystal structure, which is a simple cubic CsCl (space group Pm 3 ¯ m) structure with boron octahedra at the cube corners and a rare-earth atom occupying the body-center position. In this structure, each rare-earth atom is surrounded by eight boron octahedra, and the boron octahedra are linked together into a three-dimensional network. It is generally accepted that the structural peculiarity should be taken as one of the basic starting points to explain some of the anomalies in the R E B 6 systems. Currently, research focuses on the developments and applications of the R E B 6 nanostructures. Understanding the mechanical behavior of these structures is very important in the development and application stages. It is well-known that nonlinear effects become significant in nanostructural materials. Nonlinear elastic properties are important for describing nonlinear effects in mechanical behavior. Thus, it is very necessary to study nonlinear effects in the elasticity of the hexaborides.
In general, the second-order elastic constants (SOECs) C i j describe the linear elastic stress–strain response for single crystals. Third-order elastic constants (TOECs) C i j k are important quantities to characterize the nonlinear elasticity of the crystals. Both SOECs and TOECs are important parameters to model the mechanical response of crystals under high pressure. Tanaka et al. [2] have made measurements of the transit times of pulses of longitudinal and transverse ultrasonic waves propagating in single crystal LaB 6 at room temperature, and have determined its SOECs from the resultant velocities. Baranovskiy et al. [3] measured the sound velocities along the principal crystallographic axes of LaB 6 at 78 K by employing the phase-frequency method, and evaluated its SOECs and bulk modulus from the resultant sound velocities. Nakamura et al. [4,5] performed sound velocity measurements with an ultrasonic apparatus based on the phase comparison method to investigate the temperature dependence of the SOCEs in LaB 6 and CeB 6 . Goto et al. [6] also studied the temperature dependence of the SOECs in CeB 6 . Lüthi et al. [7] took the ultrasonic and Brillouin scattering measurements to redetermine the SOECs of CeB 6 at room temperature. Moreover, Gürel and Eryiǧit [8] used a first-principles calculations method to study structural, elastic, lattice-dynamical, and thermodynamical properties of LaB 6 and CeB 6 . Tang et al. [9] performed a first-principles study of structural, elastic, and electronic properties of CeB 6 under pressure. Besides, high-pressure phase transitions in LaB 6 and CeB 6 have been investigated experimentally. By using Raman and X-ray diffraction measurements on LaB 6 , Teredesai et al. [10] proposed that a phase transition from cubic to orthorhombic crystal structure occurs around 10 GPa, while Godwal et al. [11] stated that there is no structural phase transitions up to at least 25 GPa. By using X-ray diffraction measurements on CeB 6 , Leger et al. [12] revealed no structural phase transition at ambient temperature up to 20 GPa. Foroozani et al. [13] failed to detect any change in crystal structure up to 85 GPa, but they did not exclude the possibility of structural changes in the pressure range of 85−122 GPa. Nevertheless, the TOECs and related elastic properties of LaB 6 and CeB 6 have not yet been investigated either experimentally or theoretically to the best knowledge of the authors.
In recent years, first-principles calculations based on density functional theory (DFT) have been employed to successfully determine the SOECs, TOECs, and higher-order elastic constants of single crystals by utilizing a series of homogeneous deformation strains applied to a crystalline system to obtain the energy–strain relations [14,15,16]. To study the nonlinear effects in the elasticity of LaB 6 and CeB 6 , we shall use the same method to determine the SOECs and TOECs of both hexaborides. Subsequently, the pressure derivatives of the effective SOECs have been estimated from the obtained values of SOECs and TOECs. The polycrystalline bulk, shear, and Young’s moduli, Poisson’s ratio, and their pressure dependence in LaB 6 and CeB 6 have also been studied along with the elastic anisotropy of the two hexaborides. The paper is organized as follows. In Section 2, a brief description of computational methodology is given. In Section 3, the results we have obtained are presented with available experimental and theoretical values for comparison. Finally, the conclusions are drawn in Section 4.

2. Computational Methods

2.1. First-Principles Total-Energy Calculations

First-principles calculations have been performed by means of the Vienna Ab initio simulation package (VASP) code based on density functional theory (DFT) [17,18,19]. The projector augmented wave (PAW) method was used for describing the ion–electron interaction [20,21]. The generalized gradient approximation (GGA) of Perdew–Burke–Ernzerhof (PBE) was used for evaluating the exchange-correlation energy [22,23]. The standard PAW potentials were used for La, Ce, and B elements. A cutoff energy of 600 eV was chosen for the plane wave basis. A threshold of 10 6 eV per atom on total energy was set for the convergence of electronic self-consistency. A 15 × 15 × 15 Monkhorst–Pack grid of k-point was adopted for sampling the Brillouin zone [24]. Before calculating the elastic constants, the structures of LaB 6 and CeB 6 were fully relaxed with respect to the volume, shape, and internal atomic position until the atomic forces were less than 0.01 eV/Å. To accurately calculate the elastic constants, the linear tetrahedron method was used for the final self-consistent calculations of total-energies.

2.2. SOECs and TOECs of Single Crystal

In this paper, the method of the finite-strain continuum elasticity theory is employed to calculate the SOECs and TOECs. Here we discuss the theory briefly, the details of which have been given in References [25,26,27,28,29]. Let a i be the initial Cartesian coordinates of a material point in the unstrained state. A finite homogeneous deformation carries the material point to the final position with the coordinates x i in the strained state. After introducing the Jacobian deformation gradient
J i j = x i a j ( i , j = 1 , 2 , 3 ) ,
the Lagrangian strain tensor in the stressed state may be defined as
η i j = 1 2 r = 1 3 ( J r i J r j δ i j ) .
The elastic strain energy ( Δ E ) can be expanded in a Taylor series in terms of the strain tensor as
Δ E = V 2 ! i j C i j k l η i j η k l + V 3 ! i j k l m n C i j k l m n η i j η k l η m n + O ( η 4 ) ,
where V is the volume of the unstrained lattice. After applying the Voigt notation (11→1, 22→2, 33→3, 23→4, 13→5, and 12→6) to denote the strain tensor, the strain energy can be rewritten as
Δ E = V 2 ! i j C i j η i η j + V 3 ! i j k C i j k η i η j η k + O ( η 4 ) .
For cubic systems, there are three independent SOECs ( C 11 , C 12 , C 44 ) and six dependent TOECs ( C 111 , C 112 , C 123 , C 144 , C 155 , C 456 ). To obtain the complete set of SOECs and TOECs of LaB 6 and CeB 6 , we used six Lagrangian strain tensors in terms of a single strain parameter ξ . Table 1 gives the relationship between the coefficients A 2 and A 3 and SOECs and TOECs for the six selected strain tensors. For each strain tensor, the strain parameter ξ varied from −0.08 to 0.08 in steps of 0.01. Inserting these strain tensors into Equation (4), the strain energy density Φ can be written as an expansion in the strain parameter ξ as
Φ = Δ E V = 1 2 A 2 ξ 2 + 1 6 A 3 ξ 3 + O ( ξ 4 ) ,
where A 2 and A 3 are the combinations of SOECs and TOECs, respectively. For every deformed configuration, the atomic positions were optimized, and the total-energy was calculated by using first-principles method based on DFT. The strain energy is defined as the total-energy difference between the deformed and the perfect crystals. In this way, the dependencies of the strain energy Δ E on the strain parameter ξ were obtained for each homogeneous deformation. By comparing with the expressions from the finite-strain elasticity theory given in Table 1, the elastic constants were extracted from a polynomial fit to the strain energy versus strain parameter curves under the various strains.

2.3. Pressure Derivatives of the Effective SOECs

When an external hydrostatic pressure is applied to a crystal, the effective SOECs are very useful to describe the nonlinear elastic properties of the crystal. Usually, the effective SOECs under hydrostatic pressure P ( C i j ( P ) ) can be expanded by a Taylor expansion as [15]
C i j ( P ) C i j + d C i j ( P ) d P P = C i j + C i j P ,
where C i j is the first-order pressure derivative, which can be determined from SOECs and TOECs, and can be expressed for cubic systems as [15,25]
C 11 = C 111 + 2 C 112 + 2 C 11 + 2 C 12 C 11 + 2 C 12 , C 12 = 2 C 112 + C 123 C 11 C 12 C 11 + 2 C 12 , C 44 = 2 C 155 + C 144 + C 11 + 2 C 12 + C 44 C 11 + 2 C 12 .

2.4. Pressure Derivatives of Polycrystalline Elastic Moduli

On the basis of the effective elastic constants, the bulk modulus B and shear modulus G for LaB 6 and CeB 6 under different pressure were obtained using the Voigt, Reuss, and Hill approximations [30,31,32]. For the specific case of cubic structures, the Voigt’s and Reuss’s bulk moduli can be expressed as
B V = B R = C 11 + 2 C 12 3 ,
and the Voigt’s and Reuss’s shear moduli are defined as
G V = C 11 C 12 + 3 C 44 5 , G R = 5 ( C 11 C 12 ) C 44 3 ( C 11 C 12 ) + 4 C 44 .
Hill proposed that the effective elastic moduli are the arithmetic averages of the Voigt and Reuss moduli, and thus obtained by
B H = B V + B R 2 , G H = G V + G R 2 ,
where the subscripts “V” and “R” correspond to the Voigt and Reuss bounds, and the subscript “H” represents the Hill averaging method. The Young’s modulus (E) and Poisson’s ratio ( ν ) are given by
E X = 9 B X G X 3 B X + G X , ν X = 3 B X 2 G X 2 ( 3 B X + G X ) ,
where X = V , R , H . The pressure derivative of the bulk and shear moduli can be given by
B H = B V = B R = C 11 + 2 C 12 3 , G V = C 11 C 12 + 3 C 44 5 , G R = 5 [ 3 ( C 11 C 12 ) 2 C 44 + 4 ( C 11 C 12 ) C 44 2 ] [ 3 ( C 11 C 12 ) + 4 C 44 ] 2 , G H = G V + G R 2 .
The pressure derivative of the Young’s modulus E X and the Poisson’s ratio ν X can be given by
E X = 9 ( 3 B X 2 G X + B X G X 2 ) ( 3 B X + G X ) 2 , ν X = 9 ( B X G X B X G X ) 2 ( 3 B X + G X ) 2 .
Under the external hydrostatic pressure P, the polycrystalline elastic moduli Y X ( P ) ( Y = B , G , E , ν ) can be given by
Y X ( P ) = Y X + d Y X ( P ) d P P = Y X + Y X P .

3. Results and Discussion

3.1. Second-Order and Third-Order Elastic Constants of R E B 6 ( R E = La, Ce)

The calculated results of the lattice parameters and the SOECs of LaB 6 and CeB 6 at ground-state are listed in Table 2 along with the published experimental values and other calculated results [2,3,4,5,6,7,8,9,13,33,34]. Comparison of the lattice constants shows that the calculated quantities of the present work are in excellent agreement with the previous experimental and theoretical values. The maximum relative error between our calculated result and the experimental data is 0.07 % ( 0.65 % ) for LaB 6 (CeB 6 ), and that between our and other calculated results correspond to 0.63 % ( 0.96 % ) for LaB 6 (CeB 6 ). The room-temperature lattice constant of LaB 6 was measured as 4.156 Å from the X-ray power diffraction patterns [2], which is larger than that (4.1407 Å) of CeB 6 [34]. It is the same with the results obtained in this study. These can be explained well by the larger atomic radius of La (2.74 Å) compared with Ce (2.70 Å) [35].
The strain energy versus strain parameter curves under the various strains are fitted with a suitable polynomial to obtain the coefficients A 2 and A 3 in the Equation (5) for determining the SOECs and TOECs, as illustrated in Figure 1. The discrete points and the solid lines represent the results obtained from the first principles calculations and fitted polynomials, respectively. Obviously, these dependent curves of the strain energy on the strain have the characteristics of the asymmetry, which is the expected behavior under finite-strain elastic deformation. The strain energy with negative strains is always larger than that with positive strains, and thus the TOECs are typically negative. The fitted curves match well with the first-principles calculations results. For the values of SOECs, it is possible that some constants may be determined from a few polynomial fits (e.g., C 44 from coefficients in f 4 ( η ) , f 5 ( η ) , and f 6 ( η ) ), together with obtaining slightly different results (e.g., for LaB 6 , C 44 = 90.2, 89.1, and 91.0 GPa from f 4 ( η ) , f 5 ( η ) , and f 6 ( η ) ). In such cases, the average value of all results is given in Table 2. Measurements of the SOECs for LaB 6 and CeB 6 exhibit large discrepancies among themselves, partly because of different techniques to obtain them. Especially, the SOEC C 12 of CeB 6 was measured to be large negative in reference [6] while positive in references [5,7]. For LaB 6 , the calculated SOECs of the present study are in reasonable accordance with the previous experimental and theoretical results [2,3,4,8]. For CeB 6 , a negative value of C 12 reported in reference [6] was also not found in the present study. The calculated results of the SOECs are also in reasonable agreement with other experiments [5,7] and previous calculations [8,9]. For the two hexaborides, the SOEC C 11 measures the resistance to linear compression along the uniaxial axes, and the others are mainly related to the non-axial sound propagation. The values of C 11 are significantly larger than the others in both systems, implying that they are difficult to compress under uniaxial stress. In addition, the ground-state SOECs of LaB 6 and CeB 6 can satisfy the three Born stability criteria for the cubic system [36]: C 11 + 2 C 12 > 0 , C 11 C 12 > 0 and C 44 > 0 , and thus their cubic structures are both mechanically stable.
The TOECs allow the determination of anharmonic properties of crystals, such as thermal expansion, interactions of acoustic and thermal phonons, and temperature and pressure dependence of elastic constants. The evaluation of the TOECs is of general interest. The calculated TOECs of LaB 6 and CeB 6 are given in Table 3. All the TOECs of the two hexaborides are negative except for the values of C 123 (=254.3 GPa (304.8 GPa) for LaB 6 (CeB 6 )). The absolute values of their longitudinal mode C 111 (= 2647.1 GPa ( 2568.4 GPa) for LaB 6 (CeB 6 )) were found to be much greater than the corresponding shear modes C 112 , C 123 , C 144 , C 155 , and C 456 , implying the contribution to lattice vibrations from their longitudinal modes is much greater. The highest absolute values of C 111 indicate a pronounced anisotropy in both materials. Unfortunately, no measurements or calculations on the TOECs of LaB 6 and CeB 6 are available for comparison.

3.2. Pressure Derivatives of the Effective Second-Order Elastic Constants of R E B 6 ( R E = La, Ce)

The TOECs above were further used to evaluate the first-order pressure derivatives of LaB 6 and CeB 6 . These values are summarized in Table 4. A linear increase with pressure was observed for C 11 , C 12 , C 44 with pressure derivative C i j of 4.586, 1.895, 0.451 and 4.507, 1.921, 0.085 for LaB 6 and CeB 6 . The pressure-induced variation in the longitudinal mode C 11 was the largest, followed by the shear mode C 12 , and the smallest for the pure shear mode C 44 . For the two hexaborides, the effect of the pressure on the C 12 was obviously smaller than that on the C 11 , but markedly greater than that on the C 44 . Under hydrostatic pressure, the Born stability criteria for the cubic system are given by [9]
K 1 = C ˜ 11 P + 2 C ˜ 12 P > 0 , K 2 = C ˜ 11 P C ˜ 12 P > 0 , K 3 = C ˜ 44 P > 0 ,
with
C ˜ i i P = C i i ( P ) P ( i = 1 , 4 ) , C ˜ 12 P = C 12 ( P ) + P .
In terms of the effective SOECs, these criteria can be expressed as
K 1 = ( C 11 + 2 C 12 ) + ( C 11 + 2 C 12 + 1 ) P > 0 , K 2 = ( C 11 C 12 ) + ( C 11 C 12 2 ) P > 0 , K 3 = C 44 + ( C 44 1 ) P > 0 .
From the pressure derivatives of the SOECs, the first-order pressure derivatives of K 1 and K 2 were calculated as 9.376 (9.349) and 0.691 (0.586) for LaB 6 (CeB 6 ), while that of the corresponding K 3 was calculated to be 0.549 ( 0.915 ). Thus, the K 1 and K 2 values of both materials linearly increase while the K 3 values linearly decrease with increasing pressure. When the pressure applied to LaB 6 (CeB 6 ) is beyond 164.1 GPa (81.7 GPa), the K 3 has a negative value, as shown in Figure 2. This indicates that the cubic structure of LaB 6 (CeB 6 ) can remain mechanically stable up to 164.1 GPa (81.7 GPa). Previous high-pressure phase transitions showed that no structural phase transitions occur up to 25 GPa (85 GPa) for LaB 6 (CeB 6 ) [11,12,13]. Therefore, the results of the present work are basically consistent with those of previous experiments.

3.3. Pressure Derivatives of the Polycrystalline elastic moduli of R E B 6 ( R E = La, Ce)

The polycrystalline bulk (B), shear (G), and Young’s (E) moduli and Poisson’s ratio ( ν ) of LaB 6 and CeB 6 were calculated based on their single crystal SOECs, which are collected in Table 5. The calculated bulk modulus (174.2 GPa) of LaB 6 agrees well with the reported experimental values of 163, 184, 188, 164 ± 2 , 173 ± 7 , and 172 GPa [2,3,4,11,37], and the previous theoretical results of 185 and 180 GPa [3,8]. For CeB 6 , the calculated bulk modulus (174.7 GPa) of the present study are also in good accordance with experimental measurements of 191, 168, 182, and 166 GPa [5,7,12] and other theoretical results of 173 and 166.8 GPa [8,9]. The calculated G value of 120.9 GPa is excellently consistent with that of 121.2 GPa reported in Ref. [9]. The pressure derivatives of the bulk modulus and shear modulus of the two hexaborides are calculated using the pressure derivatives of the single crystal SOECs, and then those of the Young’s modulus and the Poisson’s ratio are also estimated. All the obtained results are presented in Table 6. The calculated B value of 2.792 (2.783) for LaB 6 (CeB 6 ) can compare well with the experimental value of 4.2 ± 1.5 [11] (3.15 [12]). One can see that the B and G are positive for both materials, which agree well with the corresponding C i j s. A similar behavior is found for the E and ν being positive, and is in good accordance with the B and G . The bulk modulus is a measure of the resistance of a material to volume changes. The shear modulus is a measure of the resistance of a material to shear deformation. The bulk modulus and shear modulus play an important role in determining the strength of solids [38]. Usually, a superhard material has a high bulk modulus, high shear modulus, and high shear strength. Thus, we can predict that LaB 6 and CeB 6 have high hardness. The positive values of B and G mean that the elastic moduli B and G can increase gradually with the pressure, showing that the two hexaborides with a simple cubic structure became more difficult to compress and shear as the pressure increased. The Young’s modulus is a measure of the stiffness of a material. The larger the value of E, the stiffer the material. The positive values of E demonstrate that LaB 6 and CeB 6 become more and more stiff as the pressure increases.
Pugh [39] has introduced the ratio between the bulk modulus and the shear modulus ( B / G ) to assess the ductile/brittle behaviors of a material. A high (low) B / G value is correlated with the ductility (brittleness) of the material. The critical value of the brittle-to-ductile transition was observed to be ∼1.75. The B / G ratio was calculated as 1.327 (1.445) for LaB 6 (CeB 6 ) at ground-state, indicating the brittleness of the hexaboride. The pressure derivative ( B / G ) can be obtained from those of the bulk modulus and the shear modulus, which is expressed as ( B / G ) = ( B G B G ) / G 2 . On this basis, the ( B / G ) value was evaluated to be 0.014 (0.018) for LaB 6 (CeB 6 ), showing that the B / G ratio can increase with increasing pressure. The pressure at which a brittle-to-ductile transition happens is predicted as 30.1 GPa (16.5 GPa) for LaB 6 (CeB 6 ). Frantsevich et al. [40] distinguished the ductility/brittleness of the materials in terms of Poisson’s ratio. Generally, a brittle material has a lower Poisson’s ratio than 0.26. The ground-state value of the Poisson’s ratio is consistent with the B / G ratio for LaB 6 (CeB 6 ). The very small and positive values of ν imply a very slow increase of the ν of both materials with the pressure. Pettifor [41] introduced the Cauchy pressure to describe the covalent character of atomic bonding related to the ductile/brittle characteristics of a material. He suggested that larger positive Cauchy pressure corresponds to a ductile material with more metallic bonds, whereas larger negative values indicate a brittle behavior with a more covalent character of bonds. For a cubic system, the Cauchy pressure P c is defined as P c = C 12 C 44 . The P c of LaB 6 (CeB 6 ) was calculated as 65.9 GPa ( 58.6 GPa) in terms of the ground-state SOECs given in Table 2. We can find that the bonding of both materials is covalent with B / G < 1.75 , leading to a brittle behavior. The pressure derivative P c can be obtained from those of the effective SOECs, which is given by P c = C 12 C 44 . The P c value of LaB 6 (CeB 6 ) was estimated to be 1.444 (1.836) based on the values of C 12 and C 44 given in Table 4, showing the increase of the Cauchy pressure with the pressure. Overall, the brittleness of LaB 6 (CeB 6 ) can reduce as the pressure increases, which is consistent with the previous theoretical result [9].

3.4. Pressure Derivatives of the Elastic Anisotropy of R E B 6 ( R E = La, Ce)

The elastic anisotropy of a material has an important implication in engineering science due to its high association with the possibility of inducing microcracks in the material [42]. The elastic anisotropy factor for a cubic crystal introduced firstly by Zener [43] is expressed as
A Z = 2 C 44 ( C 11 C 12 ) .
A single crystal with A Z = 1 is isotropic, while values smaller or greater than unity describe the degree of elastic anisotropy. The pressure derivative A Z can be obtained from those of the effective SOECs, which is given by
A Z = [ C 44 ( C 11 C 12 ) C 44 ( C 11 C 12 ) ] ( C 11 C 12 ) 2 .
Subsequently, Chung and Buessem [44] empirically improved the anisotropy factor of Zener by
A C = ( G V G R ) ( G V + G R ) .
A single crystal with A C = 0 is isotropic, otherwise it is anisotropic. The pressure derivative A C can be obtained from those of the elastic modulus G V and G R , which is given by
A C = 2 ( G V G R G V G R ) ( G V + G R ) 2 .
However, it is noteworthy that A Z and A C do not account for the bulk part of the elastic stiffness tensor, and they are only effective for cubic crystals. To consider the contributions of the shear modulus and bulk modulus, Shivakumar et al. [45] proposed a universal anisotropy index that is applicable to various crystal systems, which was expressed as A S = B V / B R + 5 G V / G R 6 . A nonzero value of A S is a measure of the anisotropy. Because the B V is equal to the B R for cubic system, the A S can be simplified as
A S = 5 ( G V G R 1 ) .
The pressure derivative A S can be calculated in terms of those of the elastic modulus G V and G R , which is given by
A S = 5 ( G V G R G V G R ) G R 2 .
In terms of the SOECs, the elastic moduli, and the corresponding pressure derivatives above, we predicted the values and the pressure derivatives of various anisotropy factors for LaB 6 and CeB 6 , as presented in Table 7. The A Z values of less than one and the nonzero A C and A S values show the elastic anisotropy of the two hexaborides at ground-state. Meanwhile, the negative values of A Z and the positive ones of A C and A S values indicate that the corresponding anisotropy factors can be far away from unity with increasing pressure. These results show that the anisotropy of both materials can be enhanced by increasing the pressure.

4. Conclusions

The SOECs and TOECs of rare-earth hexaborides LaB 6 and CeB 6 have been determined by using first-principles calculations and homogeneous deformation methods. The calculated lattice parameters, SOECs, and bulk moduli are in reasonable accordance with the available experimental and theoretical values. All the TOECs of both hexaborides have negative values except for C 123 . The highest absolute values of the longitudinal modes C 111 show that their contribution to lattice vibrations is the greatest, and so also their anisotropy in C 111 . From the calculated elastic constants, the pressure derivatives of the effective SOECs have been investigated along with the pressure effects on the structural stability, mechanical properties, and elastic anisotropy. As the pressure increases, the mechanical stability reduces, and the ductility and anisotropy increase for both hexaborides. The full set of TOECs of both materials have been determined for the first time, and comparison could not be made as the same is not available in literature. The measurements of the TOECs are essential for the research of LaB 6 and CeB 6 in the future.

Acknowledgments

The work is supported by the Natural Science Foundation of China (51661013, 51461020, 11464020) and the Ph.D. Start-up Fund of Natural Science Foundation of Jinggangshan University (JZB15007).

Author Contributions

X.Z. and Y.W. conceived and designed the calculations; Y.Y. and Q.G. performed the calculations; P.O. checked the data; X.Z., S.Z. and Y.W. analyzed the data; X.Z. wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ji, X.H.; Zhang, Q.Y.; Xu, J.Q.; Zhao, Y.M. Rare-earth hexaborides nanostructures Recent advances in materials, characterization and investigations of physical properties. Prog. Solid State Chem. 2011, 39, 51–69. [Google Scholar] [CrossRef]
  2. Tanaka, T.; Yoshimoto, J.; Ishli, M.; Bannai, E.; Kawai, S. Elastic constants of LaB6 at room temperaure. Solid State Commun. 1977, 22, 203–205. [Google Scholar] [CrossRef]
  3. Baranovskiy, A.; Grechnev, G.; Fil, V.; Ignatova, T.; Logosha, A.; Panfilov, A.; Svechkarev, I.; Shitsevalova, N.; Filippov, V.; Eriksson, O. Electronic structure, bulk and magnetic properties of MB6 and MB12 borides. J. Alloys Compd. 2007, 442, 228–230. [Google Scholar] [CrossRef]
  4. Nakamura, S.; Goto, T.; Kasaya, M.; Kunii, S. Electron-strain interaction in valence fluctuation compound SmB6. J. Phys. Soc. Jpn. 1991, 60, 4311–4318. [Google Scholar] [CrossRef]
  5. Nakamura, S.; Goto, T.; Kunii, S.; Iwashita, K.; Tamaki, A. Quadrupole-strain interaction in rare-earth hexaborides. J. Phys. Soc. Jpn. 1994, 63, 623–636. [Google Scholar] [CrossRef]
  6. Goto, T.; Tamaki, A.; Kunii, S.; Nakajima, T.; Fujimura, T.; Kasuya, T.; Komatsubarra, T.; Woods, S.B. Elastic properteis in CeB6. J. Magn. Magn. Mater. 1983, 31–34, 419–420. [Google Scholar] [CrossRef]
  7. Lüthi, B.; Blumenröder, S.; Hillebrands, B.; Zirngiebl, E.; Güntherodt, G.; Winzer, K. Elastic and magnetoelastic effects in CeB6. J. Magn. Magn. Mater. 1985, 58, 321–322. [Google Scholar] [CrossRef]
  8. Gürel, T.; Eryiğithys, R. Abinitio lattice dynamics and thermodynamics of rare-earth hexaborides LaB6 and CeB6. Phys. Rev. B 2010, 82, 104302–104313. [Google Scholar] [CrossRef]
  9. Tang, M.; Liu, L.; Cheng, Y.; Ji, G.J. First-principles study of structural, elastic, and electronic properties of CeB6 under pressure. Front. Phys. 2015, 10, 107104–107109. [Google Scholar] [CrossRef]
  10. Teredesai, P.; Muthu, D.V.S.; Chandrabhas, N.; Meenakshi, S.; Vijayakumar, V.; Modak, P.; Rao, R.S.; Godwal, B.K.; Sikka, S.K.; Sood, A.K. High pressure phase transition in metallic LaB6: Raman and X-ary diffraction studies. Solid State Commun. 2004, 129, 791–796. [Google Scholar]
  11. Godwal, B.K.; Petruska, E.A.; Speziale, S.; Yan, J.; Clark, S.M.; Kruger, M.B.; Jeanloz, R. High-pressure Raman and X-ray diffraction studies on LaB6. Phys. Rev. B 2009, 80, 172104–172107. [Google Scholar] [CrossRef]
  12. Leger, J.M.; Rossat-Mignod, J.; Kunii, S.; Kasuya, T. High pressure compression of CeB6. Solid State Commun. 1985, 54, 995–997. [Google Scholar] [CrossRef]
  13. Foroozani, N.; Lim, J.; Fabbris, G.; Rosa, P.F.S.; Fisk, Z.; Schilling, J.S. Suppression of dense Kondo state in CeB6 under pressure. Physica B 2015, 457, 12–16. [Google Scholar] [CrossRef]
  14. Zhao, J.J.; Winey, J.M.; Gupta, Y.M. First-principles calculations of second- and third-order elastic constants for single crystals of arbitrary symmetry. Phys. Rev. B 2007, 75, 94105–94111. [Google Scholar] [CrossRef]
  15. opuszyński, M.; Majewski, J.A. Ab initio calculations of third-order elastic constants and related properties for selected semiconductors. Phys. Rev. B 2007, 76, 45202–45209. [Google Scholar] [CrossRef]
  16. Wang, H.; Li, M. Ab initio calculations of second-, third-, and fourth-order elastic constants for single crystals. Phys. Rev. B 2009, 79, 224102–224111. [Google Scholar] [CrossRef]
  17. Kresse, G.; Hafner, J. Ab initio molecular dynamics for open-shell transition metals. Phys. Rev. B 1993, 48, 13115–13118. [Google Scholar] [CrossRef]
  18. Kresse, G.; Furthmller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comput. Mater. Sci. 1996, 6, 15–50. [Google Scholar] [CrossRef]
  19. Kresse, G.; Furthmller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B 1996, 54, 11169–11186. [Google Scholar] [CrossRef]
  20. Blöchl, P.E. Projector augmented-wave method. Phys. Rev. B 1994, 50, 17953–17979. [Google Scholar] [CrossRef]
  21. Kresse, G.; Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. Rev. B 1999, 59, 1758–1775. [Google Scholar] [CrossRef]
  22. Perdew, J.P.; Burke, K.; Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 1996, 77, 3865–3868. [Google Scholar] [CrossRef] [PubMed]
  23. Perdew, J.P.; Burke, K.; Ernzerhof, M. Generalized gradient approximation made simple [Phys. Rev. Lett. 77, 3865 (1996)]. Phys. Rev. Lett. 1997, 78, 1396. [Google Scholar] [CrossRef]
  24. Monkhorst, H.J.; Pack, J.D. Special points for Brillouin-zone integrations. Phys. Rev. B 1976, 13, 5188–5192. [Google Scholar] [CrossRef]
  25. Birch, F. Finite elastic strain of cubic crystals. Phys. Rev. 1947, 71, 809–824. [Google Scholar] [CrossRef]
  26. Thurston, R.; Brugger, K. Third-order elastic constants and the velocity of small amplitude elastic waves in homogeneously stressed media. Phys. Rev. 1964, 135, 3. [Google Scholar] [CrossRef]
  27. Brugger, K. Thermodynamic definition of higher order elastic coefficients. Phys. Rev. 1964, 133, 1611–1612. [Google Scholar] [CrossRef]
  28. Hiki, Y. High-order elastic constants of solids. Annu. Rev. Mater. Sci. 1981, 11, 51–73. [Google Scholar] [CrossRef]
  29. Rao, R.R.; Padmaja, A. Effective second-order elastic constants of a strained crystal using the finite strain elasticity theory. J. Appl. Phys. 1987, 62, 440–443. [Google Scholar] [CrossRef]
  30. Voigt, W. Lehrbuch der Kristallphysik; Taubner: Leipzig, Germany, 1928. [Google Scholar]
  31. Reuss, A. Calculation of the flow limits of mixed crystals on the basis of the plasticity of monocrystals. Z. Angew. Math. Mech. 1929, 9, 49–58. [Google Scholar] [CrossRef]
  32. Hill, R. The elastic behaviour of a crystalline aggregate. Proc. Phys. Soc. A 1952, 65, 349–354. [Google Scholar] [CrossRef]
  33. Chen, C.H.; Aizawa, T.; Iyi, N.; Sato, A.; Otani, S. Structural refinement and thermal expansion of hexaborides. J. Alloys Compd. 2004, 366, L6–L8. [Google Scholar] [CrossRef]
  34. Tanaka, K.; Ōnuki, Y. Observation of 4f electron transfer from Ce to B6 in the Kondo crystal CeB6 and its mechanism by multitemperature X-ray diffraction. Acta Cryst. B 2002, 58, 423–436. [Google Scholar] [CrossRef]
  35. Clementi, E.; Raimondi, D.L.; Reinhardt, W.P. Atomic screening constants from SCF functions. II. Atoms with 37 to 86 Electrons. J. Chem. Phys. 1967, 47, 1300–1307. [Google Scholar] [CrossRef]
  36. Mouhat, F.; Coudert, F.X. Necessary and sufficient elastic stability conditions in various crystal systems. Phys. Rev. B 2014, 90, 224104–224107. [Google Scholar] [CrossRef]
  37. Lundström, T.; Lönnberg, B.; Törmä, B.; Etourneau, J.; Tarascon, J.M. An investigation of the compressibility of LaB6 and EuB6 using a high pressure X-ray power diffraction technique. Phys. Scr. 1982, 26, 414–416. [Google Scholar] [CrossRef]
  38. Teter, D.M. Computational alchemy: The search for new superhard materials. MRS Bull. 1998, 23, 22–27. [Google Scholar] [CrossRef]
  39. Pugh, S.F. Relations between the elastic moduli and the plastic properties of polycrystalline pure metals. Philos. Mag. 1954, 45, 823–843. [Google Scholar] [CrossRef]
  40. Frantsevich, I.N.; Voronov, F.F.; Bokuta, S.A. Elastic Constants and Elastic Moduli of Metals and Insulators Handbook; Frantsevich, I.N., Ed.; Naukova Dumka: Kiev, Ukraine, 1983; pp. 60–180. [Google Scholar]
  41. Pettifor, D.G. Theoretical predictions of structure and related properties of intermetallics. Mater. Sci. Technol. 1992, 8, 345–349. [Google Scholar] [CrossRef]
  42. Tvergaard, V.; Hutchinson, J.W. Microcracking in ceramics induced by thermal expansion or elastic anisotropy. J. Am. Ceram. Soc. 1988, 71, 157–166. [Google Scholar] [CrossRef]
  43. Zener, C. Elasticity and Anelasticity of Metals; University of Chicago: Chicago, IL, USA, 1984. [Google Scholar]
  44. Chung, D.H.; Buessem, W.R. The elastic anisotropy of crystals. J. Appl. Phys. 1967, 38, 2010–2012. [Google Scholar] [CrossRef]
  45. Ranganathan, S.I.; Ostoja-Starzewski, M. Universal elastic anisotropy index. Phys. Rev. Lett. 2008, 101, 55504–55507. [Google Scholar] [CrossRef] [PubMed]
Figure 1. The strain–energy relations for (a) LaB 6 and (b) CeB 6 . The discrete points and the solid lines represent the density functional theory (DFT) results and the results of third-order polynomial fitting, respectively.
Figure 1. The strain–energy relations for (a) LaB 6 and (b) CeB 6 . The discrete points and the solid lines represent the density functional theory (DFT) results and the results of third-order polynomial fitting, respectively.
Crystals 07 00320 g001
Figure 2. The pressure dependence of the K 3 values of LaB 6 and CeB 6 .
Figure 2. The pressure dependence of the K 3 values of LaB 6 and CeB 6 .
Crystals 07 00320 g002
Table 1. The coefficients A 2 and A 3 in Equation (5) of the corresponding selected strain tensors as the linear combinations of the second- and third-order elastic constants for cubic crystal [14].
Table 1. The coefficients A 2 and A 3 in Equation (5) of the corresponding selected strain tensors as the linear combinations of the second- and third-order elastic constants for cubic crystal [14].
Strain A 2 A 3
η 1 = ( ξ , 0 , 0 , 0 , 0 , 0 ) C 11 C 111
η 2 = ( ξ , ξ , 0 , 0 , 0 , 0 ) 2 C 11 + 2 C 12 2 C 111 + 6 C 112
η 3 = ( ξ , ξ , ξ , 0 , 0 , 0 ) 3 C 11 + 6 C 12 3 C 111 + 18 C 112 + 6 C 123
η 4 = ( ξ , 0 , 0 , 2 ξ , 0 , 0 ) C 11 + 4 C 44 C 111 + 12 C 144
η 5 = ( ξ , 0 , 0 , 0 , 2 ξ , 0 ) C 11 + 4 C 44 C 111 + 12 C 155
η 6 = ( 0 , 0 , 0 , 2 ξ , 2 ξ , 2 ξ ) 12 C 44 48 C 456
Table 2. Lattice parameter a (in Å) and elastic constants C i j (in GPa) of LaB 6 and CeB 6 .
Table 2. Lattice parameter a (in Å) and elastic constants C i j (in GPa) of LaB 6 and CeB 6 .
CrystalMethoda C 11 C 12 C 44
LaB 6 Present4.154474.124.390.1
Exp. [2]4.156453.318.290.1
Exp. [4,33]4.15694784384
Exp. [3]4.15654634589
The. [8]4.12774663788
CeB 6 Present4.114488.217.974.8
Exp. [5,13]4.1324731681
Exp. [7,34]4.14075081979
Exp. [7] 4725378
Exp. [6] 406 93 78
The. [9]4.1214831075
The. [8]4.1544523498
Table 3. Third-order elastic constants C i j k (in GPa) of LaB 6 and CeB 6 .
Table 3. Third-order elastic constants C i j k (in GPa) of LaB 6 and CeB 6 .
Crystal C 111 C 112 C 123 C 144 C 155 C 456
LaB 6 2647.1 373.2 254.3 167.6 340.4 403.8
CeB 6 2568.4 402.8 304.8 69.6 286.9 408.4
Table 4. Pressure derivatives C i j of the second-order elastic constants of LaB 6 and CeB 6 .
Table 4. Pressure derivatives C i j of the second-order elastic constants of LaB 6 and CeB 6 .
Crystal C 11 C 12 C 44
LaB 6 4.5861.8950.451
CeB 6 4.5071.9210.085
Table 5. Polycrystalline bulk modulus B, shear modulus G, and Young’s modulus E (in GPa), and Poisson’s ratio ν of LaB 6 and CeB 6 .
Table 5. Polycrystalline bulk modulus B, shear modulus G, and Young’s modulus E (in GPa), and Poisson’s ratio ν of LaB 6 and CeB 6 .
CrystalMethodBGE ν B / G
LaB 6 Voigt174.2144.0338.80.176
Reuss174.2118.5289.80.223
Hill174.2131.3314.80.1991.327
CeB 6 Voigt174.7138.9329.40.186
Reuss174.7102.8257.80.254
Hill174.7120.9294.60.2191.445
Table 6. Pressure derivatives of polycrystalline elastic modulus B (bulk modulus), G (shear modulus), and E (Young’s modulus), and Poisson’s ratio ν of LaB 6 and CeB 6 .
Table 6. Pressure derivatives of polycrystalline elastic modulus B (bulk modulus), G (shear modulus), and E (Young’s modulus), and Poisson’s ratio ν of LaB 6 and CeB 6 .
CrystalMethod B G E ν ( × 10 3 ) ( B / G )
LaB 6 Voigt2.7920.8092.6642.645
Reuss2.7920.6182.0892.445
Hill2.7920.7132.3792.5500.014
CeB 6 Voigt2.7830.5682.1652.942
Reuss2.7830.1961.0842.886
Hill2.7830.3821.6362.9180.018
Table 7. Anisotropic factors A of LaB 6 and CeB 6 and their pressure derivatives A ( × 10 4 ).
Table 7. Anisotropic factors A of LaB 6 and CeB 6 and their pressure derivatives A ( × 10 4 ).
Crystal A Z A Z A C A C A S A S
LaB 6 0.401 1.956 0.0972.0021.07624.560
CeB 6 0.318 6.900 0.14910.7001.756147.830

Share and Cite

MDPI and ACS Style

Zeng, X.; Ye, Y.; Zou, S.; Gou, Q.; Wen, Y.; Ou, P. First-Principles Study of the Nonlinear Elasticity of Rare-Earth Hexaborides REB6 (RE = La, Ce). Crystals 2017, 7, 320. https://doi.org/10.3390/cryst7110320

AMA Style

Zeng X, Ye Y, Zou S, Gou Q, Wen Y, Ou P. First-Principles Study of the Nonlinear Elasticity of Rare-Earth Hexaborides REB6 (RE = La, Ce). Crystals. 2017; 7(11):320. https://doi.org/10.3390/cryst7110320

Chicago/Turabian Style

Zeng, Xianshi, Yuanxiu Ye, Shenlin Zou, Qingdong Gou, Yufeng Wen, and Ping Ou. 2017. "First-Principles Study of the Nonlinear Elasticity of Rare-Earth Hexaborides REB6 (RE = La, Ce)" Crystals 7, no. 11: 320. https://doi.org/10.3390/cryst7110320

APA Style

Zeng, X., Ye, Y., Zou, S., Gou, Q., Wen, Y., & Ou, P. (2017). First-Principles Study of the Nonlinear Elasticity of Rare-Earth Hexaborides REB6 (RE = La, Ce). Crystals, 7(11), 320. https://doi.org/10.3390/cryst7110320

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