Yu Zhoua,
Lan-Ting Shia,
A-Kun Liangb,
Zhao-Yi Zeng*c,
Xiang-Rong Chen*a and
Hua-Yun Gengd
aCollege of Physics, Sichuan University, Chengdu 610064, China. E-mail: xrchen@scu.edu.cn
bDepartamento de Física Aplicada-ICMUV-MALTA Consolider Team, Universitat de València, Burjassot (Valencia) 46100, Spain
cCollege of Physics and Electronic Engineering, Chongqing Normal University, Chongqing 400047, China. E-mail: zhaoyizeng@126.com
dNational Key Laboratory for Shock Wave and Detonation Physics Research, Institute of Fluid Physics, CAEP, Mianyang 621900, China
First published on 13th January 2021
The particle swarm optimization algorithm and density functional theory (DFT) are extensively performed to determine the structures, phase transition, mechanical stability, electronic structures, and thermodynamic properties of lanthanide phosphates (LaP and LaAs) in the pressure range of 0 to 100 GPa. Two novel high-pressure structures of LaP and LaAs are first reported here. It is found that LaX (X = P, As) undergo a phase transition from NaCl-type structure (Fm3m) to CsCl-type structure (P4/mmm) at 19.04 GPa and 17.22 GPa, respectively. With the elevation of the pressure, C2/m-LaP and Imma-LaAs are the most stable structures up to 70.08 GPa and 85.53 GPa, respectively. Finally, the analysis of the elastic constants and hardness confirms that the C2/m-LaP possesses hardness values up to 23.24 GPa due to the strong covalent P–P bonding and ionic La–P bonding, indicating that it is a potential hard material.
High pressure is an effective tool to discover new phenomena and explore new materials. Presently, there are some reports on the structural phase transition and electronic properties of LaX (X = P, As) under high pressure.12,16–19 At ambient conditions, LaX (X = P, As) are cubic structures of NaCl type (B1 phase) with space group of Fm3m. Under lower pressure, they undergo a structure phase transition into a tetragonal distorted structure of CsCl type (PT phase) with a space group of P4/mmm. In theory, Charifi et al.12 showed that LaX (X = P, As) transformed from cubic structure to tetragonal structure in the range of 14.6 to 22.6 GPa for LaP and 11.2 to 19.6 GPa for LaAs, respectively. Vaitheeswaran et al.20 investigated the electronic structure and the relative stability of LaP and LaAs at high pressure. They found that the phase transition from rocksalt to tetragonal occurred around 18 GPa and 11.2 GPa in the case of LaP and LaAs, respectively. In experiments, high-pressure synchrotron X-ray diffraction experiment revealed that the phase transition was found at around 24 GPa for LaP16 and the pressure-volume relationship extracted from the synchrotron X-ray diffraction pattern showed that the phase transition of LaAs occurred at around 20 GPa.19 Similar transformations have also been observed in other rare earth family members, such as PrAs,21 NdAs,21 SmAs,21 and PrSb.22 Recently, Khalid et al.4 found that the hydrostatic pressure can induce the p–d crossing of LaAs, showing a topological phase transition at about 7 GPa. However, there are few studies on lanthanide phosphates LaX (X = P, As) under higher pressures.
The purpose of this work is to study the structural, dynamical, mechanical, thermodynamic properties and electronic structures of lanthanide phosphates LaX (X = P, As) under high pressures. Using CALYPSO method,23,24 the structures of LaX (X = P, As) are investigated under pressures up to 100 GPa. It is found that the Fm3m-LaX (X = P, As) phases undergo structural phase transitions to P4/mmm-LaX (X = P, As) phases at lower pressure (19.04 GPa and 17.22 GPa). Moreover, two novel structural phases C2/m-LaP and Imma-LaAs have been proved to be the most stable structures under higher pressures.
To verify the stable phase at high pressure, the variations of enthalpy differences of LaP and LaAs with pressures are presented in Fig. 2. The relationships between volume and pressure are shown in Fig. 3. Our results show that, when the pressure increases, the cubic structures (Fm3m-LaP/LaAs) will be transformed into tetragonal structures (P4/mmm-LaP/LaAs). The corresponding phase transition pressures of LaP and LaAs are 19.04 GPa and 17.22 GPa, respectively, with the volume collapse are about 10.2% and 9.8%, which are consistent with other available researches.19,20,37 In the inset of Fig. 2a and b, the enthalpy differences are additionally presented to show clearly the phase transition at higher pressures. It is found that the lowest enthalpies for LaP and LaAs are different. For LaP, C2/m has the lowest enthalpy, while for LaAs, the enthalpy of Imma is lower than that of C2/m, so C2/m is the most stable phase of LaP and Imma is the stable phase of LaAs at 100 GPa. The corresponding transformation pressures are 70.08 GPa and 85.53 GPa, and the changes in volume collapse are about 3.1% and 1.3%, respectively. Therefore, the sequence of phase transition is Fm3m → P4/mmm → C2/m for LaP and Fm3m → P4/mmm → Imma for LaAs, respectively.
Fig. 3 Calculated volumes as a function of pressure for (a) LaP and (b) LaAs, where V0 represents the equilibrium volume of Fm3m-LaX (X = P, As). |
The lattice parameters, atomic coordinates and relative formation enthalpies of the three phases at several pressures are listed in Table 1. The formation enthalpies are defined by ΔH = H(LaX) − H(La) − H(X), where X represents P and As atom, respectively. H is the enthalpy per formula unit of the most stable structure under a given pressure. The discussions of Fm3m structures of LaP and LaAs have been mentioned in our previous studies.38,39 In the P4/mmm phase, which contains two atoms per primitive unit-cell, each La atom is surrounded by four P/As atoms. The bond lengths of La–P and La–As are 2.80 and 2.87 Å, respectively. The calculated lattice constants of LaP are a = 3.507 Å and c = 3.004 Å at 20 GPa, well consistent with the experimental value (a = 3.460 Å, c = 3.007 Å)16 and theoretical lattice constants (a = 3.530 Å, c = 3.059 Å calculated by Vaitheeswaran et al.,20 a = 3.711 Å, c = 3.82 Å calculated by Charifi et al.12). The lattice constants of LaAs are a = 3.684 Å and c = 3.108 Å at 20 GPa, consistent well with the available experimental data (a = 3.625 Å, c = 3.06 Å)19 and other theoretical value (a = 3.645 Å, c = 3.171 Å).20 Upon compression, the P4/mmm-LaP/LaAs are predicted to be the monoclinic structure C2/m for LaP and the orthorhombic structure Imma for LaAs. For the C2/m-LaP and Imma-LaAs structures, the obtained equilibrium lattice parameters are: a = 4.063 Å, b = 3.589 Å, c = 8.202 Å and a = 4.063 Å, b = 3.589 Å, c = 8.201 Å, respectively. We noticed that there is no experimental and theoretical information on the structural behavior of LaP and LaAs under high pressures. These results are of great significance for further understanding of the high-pressure characteristics of lanthanide phosphates.
LaP | Structures | Pressures | Lattice parameters (Å, deg.) | Atomic coordinates | ΔH | |
---|---|---|---|---|---|---|
Fm3m | 0 | a = b = c = 6.050 | La | (0.000, 0.000, 0.000) | −2.026 | |
α = β = γ = 90 | La | (0.000, 0.500, 0.500) | ||||
La | (0.500, 0.000, 0.500) | |||||
La | (0.500, 0.500, 0.000) | |||||
P | (0.500, 0.500, 0.500) | |||||
P | (0.500, 0.000, 0.000) | |||||
P | (0.000, 0.500, 0.000) | |||||
P | (0.000, 0.000, 0.500) | |||||
P4/mmm | 50 | a = b = 3.407, c = 2.892 | La | (0.000, 0.000, 0.500) | −0.863 | |
α = β = γ = 90 | P | (0.500, 0.500, 0.000) | ||||
C2/m | 100 | a = 4.063, b = 3.589, c = 8.202 | La | (0.859, 0.500, 0.609) | −0.137 | |
α = β = 90, γ = 137 | La | (0.141, 0.500, 0.391) | ||||
La | (0.358, 0.000, 0.609) | |||||
La | (0.641, 0.000, 0.391) | |||||
P | (0.907, 0.000, 0.157) | |||||
P | (0.093, 0.000, 0.843) | |||||
P | (0.407, 0.500, 0.157) | |||||
P | (0.593, 0.500, 0.843) |
LaAs | Structures | Pressures | Lattice parameters (Å, deg.) | Atomic coordinates | ΔH | |
---|---|---|---|---|---|---|
Fm3m | 0 | a = b = c = 6.151 | La | (0.500, 0.500, 0.500) | −1.846 | |
α = β = γ = 90 | La | (0.500, 0.000, 0.000) | ||||
La | (0.000, 0.500, 0.000) | |||||
La | (0.000, 0.000, 0.500) | |||||
As | (0.000, 0.000, 0.000) | |||||
As | (0.000, 0.500, 0.500) | |||||
As | (0.500, 0.000, 0.500) | |||||
As | (0.500, 0.500, 0.000) | |||||
P4/mmm | 50 | a = b = 3.506, c = 2.928 | La | (0.000, 0.000, 0.000) | −0.686 | |
α = β = γ = 90 | As | (0.500, 0.500, 0.500) | ||||
Imma | 100 | a = 4.063, b = 3.589, c = 8.201 | La | (0.000, 0.750, 0.872) | −0.255 | |
α = β = γ = 90 | La | (0.500, 0.750, 0.628) | ||||
La | (0.500, 0.250, 0.372) | |||||
La | (0.000, 0.250, 0.128) | |||||
As | (0.500, 0.250, 0.883) | |||||
As | (0.000, 0.250, 0.617) | |||||
As | (0.000, 0.750, 0.383) | |||||
As | (0.500, 0.750, 0.117) |
LaP | LaAs | |||||
---|---|---|---|---|---|---|
Fm3m (0 GPa) | P4/mmm (50 GPa) | C2/m (100 GPa) | Fm3m (0 GPa) | P4/mmm (50 GPa) | Imma (100 GPa) | |
C11 | 162.02 | 355.55 | 640.99 | 164.64 | 340.87 | 602.52 |
C22 | 646.96 | 452.63 | ||||
C33 | 381.12 | 444.72 | 392.67 | 595.21 | ||
C44 | 33.65 | 105.32 | 373.69 | 29.46 | 100.68 | 232.98 |
C55 | 162.44 | 165.22 | ||||
C66 | 199.68 | 251.30 | 207.64 | 308.71 | ||
C12 | 27.39 | 225.75 | 112.16 | 22.84 | 230.87 | 359.86 |
C13 | 158.09 | 199.50 | 159.31 | 124.87 | ||
C23 | 286.86 | 162.29 | ||||
C14 | ||||||
C15 | 48.81 | |||||
C25 | −77.01 | |||||
C35 | 35.48 | |||||
C46 | −53.87 |
For Fmm structures:
C11 > 0, C44 > 0, C11 > |C12|, C11 + 2C12 > 0 |
For P4/mmm structures:
C11 > |C12|, 2C132 < C33 (C11 + C12), C44 > 0, C66 > 0 |
For Imma structures:
C11 > 0, C44 > 0, C55 > 0, C66 > 0, C11C22 > C122, C11C22C33 + 2C12C13C23 − C11C232 − C22C132 − C33C122 > 0 |
For C2/m structures:
C11 > 0, C22 > 0, C33 > 0, C44 > 0, C55 > 0, C66 > 0, (C22 + C33) − 2C23 > 0, (C33C55 − C352) > 0, (C44C66 − C462) > 0, C11 + C22 + C33 + 2(C12 + C13 + C23) > 0, C22(C33C55 − C352) + 2C23C25C35 − C232C55 − C252C33 > 0, 2[C15C25(C33C12 − C13C23) + C15C35(C22C13 − C12C23) + C25C35(C11C23 − C12C13)] − [C152(C22C33 − C232) + C252(C11C33 − C132) + C352(C11C22 − C122)] + C55C11C22C33 − C11C232 − C22C132 − C33C122 + 2C12C13C23 > 0 |
It is found that the calculated elastic constants for LaX (X = P, As) (including C2/m-LaP) all satisfy the above conditions at pressure, confirming their mechanical stabilities. The calculated elastic constants for Fm3m-LaX (X = P, As) are basically consistent with the previous theoretical values.9,42,43 For P4/mmm-LaX (X = P, As), the values of C11 and C33 are similar, indicating that they have an almost isotropic linear incompressibility. For the higher-pressure phases C2/m-LaP and Imma-LaAs, their elastic constants (C11, C22 and C33) are larger than those of other phases, suggesting that they are more difficult to compress.
In addition to checking the mechanical stability, the elastic properties can also be used to explain various phenomena, such as inter-atomic potentials, phonon spectra, specific and thermal expansion.44 Table 3 gives the results of bulk modulus B (GPa), shear modulus G (GPa), G/B ratio, Young's modulus E (GPa), Poisson's ratio and Vickers hardness Hv (GPa) of LaX (X = P, As). It can be seen that C2/m-LaP and Imma-LaAs possess the larger values of bulk, shear and Young's modulus, implying the potential hard properties. Our results show that P4/mmm-LaP and P4/mmm-LaAs exhibit ductility, and others are brittle. The type of the bonding can be measured by Poisson ratio (the threshold value of Poisson's ratio is 0.25,45 the covalent bonding (ν < 0.25) corresponds to ionic materials and the ionic bonding (ν ≥ 0.25) corresponds to ionic materials. The Poisson's ratios of P4/mmm-LaX (X = P, As) and Imma-LaAs are larger than 0.25, suggesting that they are ionic crystals and the ionic bond component is dominant. While for the Poisson's ratio of C2/m-LaP, its small value implies that there are a lot of strong covalent bonds.
LaP | B (GPa) | G (GPa) | G/B | E (GPa) | ν | Hv (GPa) |
---|---|---|---|---|---|---|
Fm3m | 72.26 | 44.59 | 0.62 | 110.95 | 0.244 | 7.48 |
P4/mmm | 241.51 | 111.61 | 0.46 | 290.00 | 0.289 | 9.79 |
C2/m | 323.69 | 204.38 | 0.63 | 506.53 | 0.239 | 23.24 |
LaAs | B (GPa) | G (GPa) | G/B | E (GPa) | ν | Hv (GPa) |
---|---|---|---|---|---|---|
Fm3m | 70.11 | 42.24 | 0.60 | 105.53 | 0.249 | 6.88 |
P4/mmm | 241.43 | 107.23 | 0.44 | 279.95 | 0.292 | 8.92 |
Imma | 324.10 | 186.17 | 0.57 | 468.76 | 0.259 | 19.24 |
The hardness can be estimated by Vickers hardness Hv proposed by Chen et al.46 This empirical model is expressed as Hv = 2(k2G)0.585 − 3 (where Hv is the Vickers hardness, G is the shear modulus, and k is the value of G/B), which has been successfully applied to a wide range of material systems.47–51 For C2/m-LaP and Imma-LaAs, the hardness values are 23.24 GPa and 19.24 GPa, respectively. The hardness are much higher than those of other phases. Especially for C2/m-LaP phase whose hardness is similar to those of ReB2 (ref. 52) and WB2,27 indicating that it is a potential hard material.
For engineering materials, elastic anisotropy is an important factor affecting their applications. Fig. 4 shows the direction dependence of the Young's modulus for LaX (X = P, As) and the corresponding slice in the xy, yz, xz planes. The cubic phases of LaP and LaAs exhibit the isotropic Young's modulus. It can be seen that the anisotropy of Young's modulus of LaP-P4/mmm and LaAs-P4/mmm is basically the same. Structural anisotropy can cause the orientation difference of these compounds. The maximum values of Young's modulus are 625 GPa along xy plane for C2/m-LaP and 536 GPa along yz plane for Imma-LaAs. Accordingly, their minimum values are 388 GPa and 312 GPa, respectively. The ratios of Emax/Emin are 1.61 and 1.72, respectively, indicating the comparable anisotropy of Young's modulus.
Fig. 4 Three-dimensional and two-dimensional surface representations of Young's modulus for (a) LaP-Fm3m, (b) LaP-P4/mmm, (c) LaP-C2/m, (d) LaAs-Fm3m, (e) LaAs-P4/mmm, and (f) LaAs-Imma. |
Fig. 5 Phonon dispersions and the total and partial phonon density of states (DOSs) for (a) LaP-P4/mmm, (b) LaAs-P4/mmm, (c) LaP-C2/m and (d) LaAs-Imma. |
Phonon velocities are closely related to the phonon transportation, the phonon velocities of LaX (X = P, As) at different pressures are obtained from the phonon dispersions. The two phonon velocities (vp is the phase velocity and vg is the group velocity) are defined as:53,54 νp = ω(k)/k,νg = ∇kω(k), where k is the distance between two points in the Brillouin zone, ω(k) is the dispersion relation. Considering the anisotropy of P4/mmm-LaP and P4/mmm-LaAs, Fig. 6a–d shows the phonon phase velocities and group velocities along [100] and [001] directions. As we can see, the phase velocities of optical branches are much larger than those of acoustic modes, but in contrast, the group velocities of acoustic branches are about three times as large as those of optical branches. This phenomenon can also be found in the phonon velocities of C2/m-LaP and Imma-LaAs, as shown in the Fig. 6e and f. These results can be understood from the equations of phase velocity and group velocity. The phase velocity is directly proportional to the phonon frequency, and the phase velocities of the optical branches are larger than those of the acoustic branches because of their larger phonon frequencies. While the group velocities are obtained from the slope of the phonon dispersion, and the acoustic branches are steeper near the Γ point, so their group velocities are larger than those of the optical branches.
Fig. 6 The calculated phonon phase velocities and group velocities of LaP-P4/mmm (a) and (b), LaAs-P4/mmm (c) and (d), LaP-C2/m (e) and LaAs-Imma (f). |
Fig. 7 shows the atomic displacements of LaP and LaAs at different pressures. The displacements of La and P/As atoms are indicated by yellow and grey arrows, respectively. For P4/mmm-LaP/LaAs, the vibration direction of La atom is opposite to that of P/As atoms. For Imma-LaAs, the structure is a layered along the c direction. The two adjacent layers are composed of the same atoms, forming a folded layered structure. The vibration directions of different atoms in each layer are opposite. The most intricate one is LaP-C2/m, similar to the combination of two interactive vibrations. We assume that the primitive cell is divided into left and right parts along the b axis, and the vibration direction of different atoms in each part is different.
Fig. 7 Main atomic displacements of (a) P4/mmm-LaP, (b) P4/mmm-LaAs, (c) LaP-C2/m and (d) LaAs-Imma. The displacements of La and P/As atoms are indicated by yellow and grey arrows, respectively. |
Fig. 8 Calculated total and projected DOS for (a) LaP-P4/mmm, (b) LaAs-P4/mmm, (c) LaP-C2/m and (d) LaAs-Imma. The vertical dash lines represent the Fermi level. |
In order to further explore the contribution of each orbit near Fermi level, Fig. 9 presents the projected band structures of LaX with the spectral weight of La-f, La-d orbitals and P/As-p orbitals. The results also demonstrate the metallic characteristics of these high-pressure structures due to the absence of the band gap at the Fermi level. The observation of the projected band structures shows that P4/mmm-LaP and P4/mmm-LaAs have similar electronic distributions near Fermi level, which is in good agreement with the results of density of state (DOS). The valence band is mainly composed of mixed states of P/As-p and La-d states, while the conduction band is mainly composed of La-d states with a significant contribution from La-f and P/As-p states. For the C2/m-LaP and Imma-LaAs, their band structures remain metallic; however, the contributions of La-f orbitals near the Fermi surface are larger than those of P4/mmm-LaP and P4/mmm-LaAs.
To further study the chemical bonds properties of the La–P/As combinations, the electronic localization function (ELF) of LaX (X = P, As) are calculated as depicted in Fig. 10, which can imply the localizing degree of electrons in the certain regions.58 ELF values can provide useful information about the chemical bonding properties by qualitatively defining the bonding nature to iconicity, covalency or metallicity.44 The localizing degree of electrons in the related regions gradually increases from 0 to 1. It can be seen that the electronic localization function of P4/mmm-LaX (X = P, As) gradually changes between La and As/P atoms, which means there exist P–P covalent bonds and La–P ionic bonds. In the Fig. 10b and c, the electronic localization function for C2/m-LaP and Imma-LaAs are similar, while the overlap of electronic localization function for La–P and P–P atoms in C2/m-LaP is much larger in Imma-LaAs. This indicates the La–P ionic bonds and P–P covalent bonds in LaP-C2/m are stronger, that is to say, C2/m-LaP has larger bonding strength, so the hardness is higher than other phases. The results of electronic localization function of each phase of LaX are consistent with the analysis of density of states. The bonding strength and the proportion of the covalent play a key role in the incompressibility and hardness of the compounds.
This journal is © The Royal Society of Chemistry 2021 |