New Ternary Superconducting Compound LaRu2As2: Physical Properties from DFT Calculations
NNew Ternary Superconducting Compound LaRu As : Physical Properties from DFT Calculations M. A. Hadi , M. S. Ali , S. H. Naqib and A. K. M. A. Islam Department of Physics, University of Rajshahi, Rajshahi-6205, Bangladesh Department of Physics, Pabna University of Science and Technology, Pabna-6600, Bangladesh International Islamic University Chittagong, Chittagong-4203, Bangladesh
Abstract
In this paper, we have presented the density functional theory (DFT) based calculations performed within the first-principles pseudopotential method to investigate the physical properties of the newly discovered superconductor LaRu As for the first time. The optimized structural parameters are in good agreement with the experimental results. The calculated independent elastic constants ensure the mechanical stability of the compound. The calculated Cauchy pressure, Pugh’s ratio as well as Poisson’s ratio indicate that LaRu As should behave as a ductile material. Due to low Debye temperature, LaRu As may be used as a thermal barrier coating (TBC) material. The new compound should exhibit metallic nature as its valence bands overlap considerably with the conduction bands. LaRu As is expected to be a soft material and easily machineable because of its low hardness value of 6.8 GPa. The multi-band nature is observed in calculated Fermi surface. A highly anisotropic combination of ionic, covalent, and metallic interactions is expected in accordance with charge density calculations. Keywords : New superconductor LaRu As ; DFT calculations; Mechanical properties; Electronic features PACS:
1. Introduction
A great variety of compounds with RT X stoichiometry take on the body-centered tetragonal ThCr Si type structure (space group I
4 / mmm , No. 139) [1].
In this stoichiometry, R is symbolized for a rare earth, alkaline earth or alkali element, T refers a transition metal and X represents p-metal atoms, namely B, P, Si, As or Ge. This large family includes more than 700 members of 122-type intermetallics (so-called 122 phases) and exhibits an outstanding collection of physical and chemical properties [2]. The structure of RT X stoichiometry offers an innate multilayer system in which the planes of R ions are taken apart from T metallic layers by X atomic sheets, comprising of R–X–T–X–R stacking of individual basal planes along the c-axis. The ThCr Si type structure exhibits a strong uniaxial anisotropy with keeping c-axis as an anisotropy axis. In other words, the layer of edge sharing TX tetrahedron that is parallel to the ab plane and separated by the basal plane of R atoms characterizes the structure of RT X . In this structure, the R atomic sheets separate two adjacent [T X ] blocks. In sequence, T ions inside [T X ] blocks build up a square lattice pressed between two X sheets shifted in order that each T is enclosed by a distorted X tetrahedron TX . A wide range of variation occurs in the interlayer distances between X atoms when either R or T is altered for a particular p-metal atom [2,3]. This group of materials have attracted a lot of attention in recent years because of the discovery of high-temperature superconductivity in iron arsenide AFe As (A = K, Ca, Sr, Ba, etc.) [4,5]. As Ru lies in the same group together with Fe in the periodic table, the superconductivity in the level of iron arsenide is anticipated for Ru-related compounds, and several studies have been devoted to these materials [6-10]. Motivated by this, Guo et al. [11] have recently synthesized the polycrystalline samples of ThCr Si -type LaRu As and reported the bulk superconductivity in LaRu As . Their temperature dependent resistivity measurement ensures a superconducting transition of LaRu As with an onset T c of 7.8 K. This measurement also estimates the upper critical field H c2 (0) at zero temperature and obtains a value of 1.6 T. Further, DC magnetic susceptibility measurement detects a bulk superconducting Meissner transition at 7.0 K, and the isothermal magnetization measurement signifies that LaRu As is a type-II superconductor. Surprisingly, LaRu As has the highest transition temperature T c among the iron-free transition metal pnictides with the ThCr Si -type crystal structure. In this paper, the study of LaRu As has been carried out via the ground state electronic structure calculations using the plane-wave pseudopotential approach within the density functional formalism. The aim of the present paper is to investigate the ground state properties including structural, elastic, * Corresponding author:[email protected] lectronic and bonding characters of the newly discovered 122 phase compound by calculating its total energy at optimized volume. To the best of our knowledge, this is the first ab initio calculations completed on this new superconductor LaRu As . The rest of this paper is divided into three sections. In Section 2, a brief description of the method of computation used in this study has been presented. The results obtained for the structural, elastic, electronic, and bonding properties of LaRu As are analyzed in Section 3. Finally, Section 4 consists of the main conclusions of the present work.
2. Method of calculations
The approach used to investigate the physical properties of LaRu As as implemented in the Cambridge Serial Total Energy Package (CASTEP) code [12] is as follows: a set of Kohn-Sham equations are solved by means of the plane-wave pseudopotential method [13] within the strong electron ensemble Density Functional Theory (DFT) [14]. Exchange and correlation energies are described by a nonlocal correction for LDA in the form of GGA [15]. The wavefunctions are allowed to expand in a plane wave basis set followed by the periodic boundary conditions and Bloch’s Theorem [16]. The electron-ion potential is treated by means of first-principles pseudopotentials within Vanderbilt-type ultrasoft formulation [17]. BFGS energy minimization technique [18] best for crystalline materials is used to find, self-consistently, the electronic wavefunctions and consequent charge density. In addition, the density mixing [19,20] scheme is employed. For k -points sampling integration over the first Brillouin, the Monkhorst–Pack scheme [21] is used. The parameters are chosen for present calculations are as follows: energy cutoff 500 eV (340 eV for elastic properties calculations), k -point mesh 16
7, the difference in total energy per atom within 5 × 10 -6 eV, maximum ionic Hellmann-Feynman force within 0.01 eV/Å, maximum stress within 0.02 GPa, and maximum ionic displacement within 5 × 10 -4 Å.
3. Physical properties
In this section the investigated physical properties of the 122 phase superconductor LaRu As are presented and analyzed; among them are the structural and elastic properties, Debye temperature, Band structure and DOS, and then Mulliken population and Vickers hardness and the electron charge density as well as Fermi surface of LaRu As under zero pressure. The body centered tetragonal LaRu As structure shown in Fig. 1 crystallizes with space group I mmm (No. 139). Every atomic species in this crystal occupy only one crystallographic position and only one internal parameter z X settles the relative position of As inside the unit cell. There are two formula units in each unit cell. The 2a, 4d and 4e atomic positions are inhabited by the La, Ru, and As atoms, respectively. The negatively charged [Ru As ] blocks and positively charged La layers take positions by turns along the c -direction and form the crystal structure. The optimized lattice constants and z X parameter of LaRu As are given in Table 1 along with those obtained for some other 122 phase compounds. The computed lattice parameters accord well with the measured values [11]. Moreover, the calculated structural parameters are slightly larger ( a ~ 0.23%, c ~ 2.94% and V ~ 3.41%) than the experimental values, which is a general trend inherent to GGA calculations. Table 1 . Structural properties of LaRu As in comparison with those for the related ThCr Si type compounds. Compounds a (Å) c (Å) V (Å ) c / a z X References LaRu As As As Fig. 1.
Crystal structure of new 122 phase superconductor LaRu As . The mechanical behavior of solids can be described successfully by the elastic constants that are related to materials’ response to an applied stress. The finite strain theory assigned in CASTEP code is well established for calculating the elastic constants successfully for a range of materials including metallic systems [24]. In this method, a given uniform deformation ε j is applied and then the resulting stress σ i is calculated. By choosing an appropriate functional deformation, elastic constants are then evaluated by solving the linear equation, σ i = C ij ε j . The well-established Voight-Reuss-Hill (VRH) approximations validated in many metallic and insulating materials [25-27] are applied to calculate the polycrystalline bulk elastic properties namely, bulk modulus B and shear modulus G from calculated C ij . Further, the equations, Y = (9 GB )/(3 B + G ) and v = (3 B – 2 G ) /(6 B + 2 G ) are employed to determine the Young’s modulus Y and Poisson’s ratio v , respectively. The values obtained for single crystal elastic constants ( C , C , C , C , C and C ) and bulk elastic properties ( B , G , Y , B / G and v ) of tetragonal crystal system LaRu As are listed in Table 2 along with the results found in literature for La-based isostructural compounds [28]. The mechanical stability conditions for tetragonal crystals [29] are as follows: C > 0, C > 0, C > 0, C > 0, C – C > 0, C + C – 2 C > 0, 2( C + C ) + C + 4 C > 0. Satisfying these criteria the newly synthesized 122 phase appears as mechanically stable compound. The elastic constants C and C explain the linear compression resistance along the a- and c-directions, respectively. It is evident that the obtained elastic constants C and C are very large in comparison to other elastic constants, implying that the LaRu As crystal should be so much incompressible under uniaxial stress along both a- and c-directions. Again, the elastic constant C is much larger than C , meaning that the incompressibility along the a-direction should be much higher than that along c-direction. In fact, the bonds aligned to the a-axis contribute a dominating effect on C making it much larger than C . Since C + C > C , the bonding in the (001) plane should exhibit more elastic rigidity than that along the c-axis and the elastic tensile modulus should be higher on the (001) plane than that along the c-axis. The elastic constant C indicates the ability to resist the shear distortion in (100) plane, whereas the elastic constant C correlates to the resistance to shear in the <110> direction [30]. Since C > C , the new compound should be more competent to resist the shear distortion in the <110> direction than in the (100) plane. The shear anisotropy factor A , defined as A = 2 C /( C – C ) is 1.24, signifying that the shear elastic properties of the (001) plane in LaRu As depend on the shear directions significantly. The concept of failure mode i.e., ductile or brittle nature, of materials is important in study of mechanical behaviors of solids. A material can be classified as either brittle or ductile for most practical situations. To quantify the failure state of solids, Cauchy pressure, and Pugh and Poisson’s ratios are used as powerful tools. A brittle material changes its volume easily and a ductile material can be easily distorted under the action of external loads . The Cauchy pressure [31] defined as ( C – C ) identifies a solid as either brittle or ductile one. When Cauchy pressure is positive (negative), the aterial is prone to be ductile (brittle). In accordance with the Pugh’s ratio [32], the high (low) B / G value signifies the ductile (brittle) nature of the materials. The threshold value of B / G is found to be 1.75, which separates the ductile materials from brittle ones. Frantsevich et al. [33] separated the ductile materials from brittle materials in terms of Poisson’s ratio. This rule proposes v ~ 0.26 as the border line that separates the brittle and ductile materials. A material having the Poisson’s ratio greater than 0.26 will be ductile, otherwise the material will be brittle. The Cauchy pressure ( C – C ) of LaRu As is positive, its Pugh’s ratio B / G is 1.79 which is greater than1.75 and its Poisson’s ratio v is 0.265 which is greater than 0.26. As a result, the compound LaRu As should behave as a ductile material. It should be mentioned that there are some papers [34–37] in which B / G = 2.0 ( G / B = 0.5) and v ~ 0.33 are used as critical values for separating the brittle materials from ductile ones. On the other hand, many authors [38–44] used the current values presented in this paper to judge brittle/ductile nature of the materials. In fact, Vaitheeswaran et al. [39] prefer B / G = 1.75 as threshold value for defining the brittle/ductile nature of the materials and showed that Pugh’s critical value corresponds to v = 0.26 as the border line between the ductile and brittle materials. These criteria are frequently used to identify the brittle material from ductile ones. Table 2.
Single crystal elastic constants C ij , polycrystalline bulk modulus B , shear modulus G , and Young’s modulus Y , in GPa, Pugh’s ratio G / B , and Poisson’s ratio v of LaRu As . Compounds C C C C C C B G Y B / G v A Ref. LaRu As
244 86 82 113 62 98 113 63 159 1.79 0.265 1.24 [This] LaNi Ge
152 87 40 126 87 90 62 74 159 0.84 0.073 2.77 [28] LaNi P
202 116 15 102 116 97 62 93 186 0.67 0.002 2.26 [28]
For describing the elastic behavior of solids, bulk and shear moduli are two important parameters. Fracture and plastic deformation are two essential issues related to solid materials. Bulk modulus B measures the resistance to fracture and shear modulus G evaluates the resistance to plastic deformation of polycrystalline materials [45]. From dislocation theory [46] it is known that a weak relationship is observed between bulk modulus and hardness. On the contrary, a good relationship exists between hardness and shear modulus [47]. Certainly, the hardness is less susceptible to the bulk modulus in comparison to the shear modulus. Therefore, as can be seen from Table 2, the hardness of LaRu As is expected to be small and LaRu As should be soft and easily mechinable. Many important issues related to physical properties of solids can be addressed with Young’s modulus. The resistance in opposition to longitudinal tension can be estimated with Young’s modulus Y . The Young’s modulus also manipulates the thermal shock resistance of solids because of having the relation between critical thermal shock coefficient R and Y . In fact, the critical thermal shock coefficient R is found to be changed inversely proportional to the Young’s modulus Y [48]. The larger the R value, the better the thermal shock resistance. The thermal shock resistance is a fundamental parameter for thermal barrier coating (TBC) materials selection. In comparison to a promising thermal barrier coating material Y Al O [49], the small Y value of LaRu As signifies that it should have a reasonable resistance to thermal shock. Debye temperature D is a distinctive temperature of crystals at which the highest-frequency mode (and hence all modes) of vibration is excited. It controls many physical properties of materials and is used to make a division between high- and low-temperature regions for a solid. To distinguish the classical and quantum behavior of phonons, the Debye temperature also defines a border line between them. When the temperature T of a solid is elevated over D , every mode of vibrations is expected to be associated with an energy equal to k B T . At T < D , all high-frequency modes seem to be sleep. The vibrational excitations at low temperature are due to acoustic vibrations. The Debye temperature, D , is determined by using one of the standard methods depending on the elastic moduli [50]. This method guides us to calculate the Debye temperature from the average sound velocity (cid:1874) (cid:3040) by the following equation: (cid:2016) (cid:2888) = ℎ/(cid:1863) (cid:2886) [(3(cid:1866)/4(cid:2024))(cid:1840) (cid:2885) (cid:2025)/(cid:1839)] (cid:2869)/(cid:2871) (cid:1874) (cid:3040) , where h stands for the Planck’s constant, k B represents the Boltzmann’s constant, n is the number of atoms per formula unit, M is symbolized for molar mass, A is the Avogadro’s number and denotes the mass density of solid. The average sound velocity within a material is obtained by (cid:1874) (cid:3040) = [1/3(1/(cid:1874) (cid:3039)(cid:2871) + 2/(cid:1874) (cid:3047)(cid:2871) )] (cid:2879)(cid:2869)/(cid:2871) . Here v l and v t indicate the longitudinal and transverse sound velocities in crystalline solids and are expressed as (cid:1874) (cid:3039) =[(3(cid:1828) + 4(cid:1833))/3(cid:2025)] (cid:2869)/(cid:2870) and (cid:1874) (cid:3047) = [(cid:1833)/(cid:2025)] (cid:2869)/(cid:2870) . The Debye temperature in cooperation with longitudinal, transverse and average sound velocities calculated within the present formalism is listed in Table 3. The values of v l , v t , v m , and D found in literature [28] are calculated by using elastic moduli B and G within the Voigt approximation. Here we converted these values into aggregate values using Hill approximations [27]. As a rule of thumb, a lower Debye temperature involves in a lower phonon thermal conductivity. The examined phase LaRu As has the lower Debye temperature in comparison with other two isostructural compounds as well as a candidate material for thermal barrier coating Y Al O [49] and hence has a lower thermal conductivity. Accordingly, LaRu As should have opportunity to use as a thermal barrier coating (TBC) material. Table 3 . Calculated density ( in gm/cm ), longitudinal, transverse and average sound velocities ( v l , v t , and v m in km/s) and Debye temperature ( D in K) of LaRu As . Compounds v l v t v m D Refs. LaRu As Ge P Fig. 2.
Electronic structures of LaRu As ; (a) electronic band structure (left panel), (b) total and partial densities of states (right panel). Understanding of electronic features (band structure, DOS, etc) is important to explain many physical phenomena, namely optical spectra and transport properties of solids.
The full picture of energy bands and band gaps of a solid is known as electronic band structure or simply band structure. Actually, in solid-state and condensed matter physics, the band structure defines certain ranges of energy that are allowed for electron s within a solid, and the ranges of energy that are not allowed for any electrons. The number of bands indicates the available number of atomic orbitals in the unit cell. The calculated energy band structure with high symmetry points in the first Brillouin zone using equilibrium lattice parameters is shown in Fig. 2a. The dashed line in the band structure corresponds to the energy of the highest filled state. This energy level is known as Fermi level, E F . Depending on the position of E F , D e n s it y o f s t a t e s ( s t a t e s p e r un it ce ll p e r e V ) Energy (eV) -20 -16 -12 -8 -4 0 40246810 As 4sAs 4pRu 5sRu 4p Ru 4dLa 6s
La 5p
La 5dTotal E F High symmetry points E n e r gy ( e V ) -20-18-16-14-12-10-8-6-4-2024 Z A M Z R X E F everal important distinctions regarding the expected electrical conductivity of a material can be made. The Fermi level of new 122 phase superconductor takes position just above the valence bands maximum near the point. The new compound should behave metallic in nature since its valence bands cross the Fermi level and overlap noticeably with the conduction bands. In addition, no band gap is found in the Fermi level. The near-Fermi bands show a complicated ‘mixed’ character, combining the quasi-flat bands with a series of high-dispersive bands intersecting the Fermi level. The detailed features of band structure can be explained with the calculated total and partial densities of states. The calculated total and partial densities of states are shown in Fig. 2(b), where the vertical broken line denotes the Fermi level, E F . At the left of the Fermi level a deep valley is observed, which is known as pseudogap EP. The delocalized electrons occupy the levels above E p and the materials containing delocalized electrons turn into metalized. It is seen that the Fermi level of LaRu As lies in a region of relatively low DOS. Consequently, the new compound should be stable in perspective of electronic features as stated in the free electron model. [51]. The new 122 phase superconductor LaRu As should exhibit metallic conductivity due to having finite value of TDOS at E F . The atom resolved partial DOS gives more insights on the electronic structure. Mainly Ru 4d states contribute to the DOS at the Fermi level along with little contribution from La 5d and As 2p states. At E F the calculated TDOS has a value of 5.4 states/eV-unit cell. The large TDOS is a sign of high metallicity the compound. The TDOS at the Fermi level is quite high compared to some other recently discovered ternary superconducting compounds [52, 53]. This may be responsible for the higher superconducting transition temperature in LaRu As , compared to, for example, in ternary silicides [54]. The lowest lying valence bands extended from -17.9 to -16.3 eV arises entirely from La 5p states. These valence bands are observed to be separated by a wide prohibited energy gap of ~3.3 eV form the next valence bands lying in the energy range from -13.0 eV to -10.6 eV. These valence bands are again separated by a forbidden energy gap of width 3.4 eV from the highest valence bands. The highest valence bands contain three visible peaks and broaden in the energy range from -6.6 eV to Fermi level. The left peak originates mainly from As 4p and Ru 4d states. The middle and highest peak are due to equal contribution from Ru 4d and As 4p states. The first peak near to Fermi level forms due to contribution coming from Ru 4d states. In the energy range between -6.6 eV and the Fermi level, the covalent bonding between comprising elements is expected. This is due to the reason that the states are truly degenerate regarding angular momentum and lattice site. Because of the difference in the value of electronegativity among the constituent elements, few ionic interactions can be predicted within LaRu As as well. The allocation of electrons in several fractional manners among the various parts of bonds can be explained with Mulliken population analysis. The overlap population provides a relation between covalency/ionicity of bonding and bond strength. Sanchez-Portal et al. [55] developed a method for population analysis in CASTEP code using a projection of the plane wave states onto a localized basis. The Mulliken scheme [56] is then applied for population analysis of the resulting projected states. Analyzing the electronic structure calculations with Linear Combination of Atomic Orbital (LCAO) basis sets, this technique is used extensively. Using minimal basis Mulliken scheme [56,57] provides two fundamental quantities regarding atomic bond population are the effective charge and the bond order values between a pair of bonding atoms as follows: (cid:1843) (cid:3080)∗ = (cid:3533) (cid:3533) (cid:3533) (cid:1829) (cid:3036)(cid:3080)∗(cid:3041) (cid:1829) (cid:3037)(cid:3081)(cid:3041) (cid:1845) (cid:3036)(cid:3080),(cid:3037)(cid:3081)(cid:3037),(cid:3081)(cid:3041) (cid:3042)(cid:3030)(cid:3030)(cid:3036),(cid:3080) (cid:2025) (cid:3080)(cid:3081) = (cid:3533) (cid:3533) (cid:1829) (cid:3036)(cid:3080)∗(cid:3041) (cid:1829) (cid:3037)(cid:3081)(cid:3041) (cid:1845) (cid:3036)(cid:3080),(cid:3037)(cid:3081)(cid:3036),(cid:3037)(cid:3041) (cid:3042)(cid:3030)(cid:3030) where i , j denote the orbital quantum numbers, n indicates the band index, (cid:1829) (cid:3036)(cid:3080)∗(cid:3041) as well as (cid:1829) (cid:3037)(cid:3081)(cid:3041) are the eigenvector coefficients of the wave function and (cid:1845) (cid:3036)(cid:3080),(cid:3037)(cid:3081) is the overlap matrix between atoms and . or understanding the bonding character in crystal systems the effective valence charge and Mulliken atomic population are two important parameters. The effective valence is the difference between the formal ionic charge and Mulliken charge on the anion species within a crystal. It measures the level of covalency/ionicity of chemical bonds. A chemical bond becomes an ideal ionic bond if there is present zero effective valences. In contrast, the positive effective valence indicates the presence of covalent bond and the high values signifies the high level of covalency in chemical bonds. Table 4 lists the effective valence, which reveals the high degree of covalency in chemical bonds inside the new superconductor LaRu As . Table 5 tabulates the obtained bond overlap populations between nearest neighbors in this crystal. When the values of overlap population tend to be zero a trivial interaction occurs between the electronic populations of the two atoms. Actually, a weak bond occurs when small Mulliken population exists and such bond makes a negligible contribution to materials hardness. With a low overlap population a bond exhibits high degree of ionicity, whereas with a high value the chemical bond possesses the high level of covalency. Positive and negative bond overlap populations lead to bonding and antibonding states, respectively. Table 4.
Population analysis of LaRu2As2. Species Mulliken Atomic populations Effective valence Charge (e) s P d Total Charge (e) As 1.57 3.31 0.00 4.88 0.12 -- Ru 2.52 6.89 7.36 16.76 – Table. 5.
Calculated bond and total Vickers hardness v H , H v (in GPa) of LaRu As along with bond number n , bond length d (Å), bond volume b v (Å ) and bond and metallic populations P , P . Bond n d P P b v v H H v Ru–Ru 2 2.96437 0.98 0. 11832 13.23 8.62 6.8 Ru–La 8 3.43822 1.47 0. 11832 20.64 6.44
Hardness assesses the ability of a material to refuse the plastic deformation. The amount of force per unit area acted opposite to plastic deformation takes part in estimating the hardness of a material. F. M. Gao [58] formulated an equation to determine the hardness of non-metallic materials using Mulliken population within first-principles approach. There is no direct relation between delocalized metallic bonding and hardness of materials [59]. For this reason the above method cannot predict the hardness of materials with partial metallic bonding. Gou et al. [60] developed an equation considering a correction for metallic bonding to calculate the bond hardness of mixed metallic crystals as follows: (cid:1834) (cid:3049)(cid:3091) = 740 (cid:4672)(cid:1842) (cid:3091) − (cid:1842) (cid:3091) ′ (cid:4673) ((cid:1874) (cid:3029)(cid:3091) ) (cid:2879)(cid:2873)/(cid:2871) where Mulliken overlap population of the -type bond is presented by (cid:1842) (cid:3091) , (cid:1842) (cid:3091) ′ is a symbol for metallic population and is calculated by using the unit cell volume V and the number of free electrons in a cell (cid:1866) (cid:3033)(cid:3045)(cid:3032)(cid:3032) = ∫ (cid:1840)((cid:1831))(cid:1856)(cid:1831) (cid:3006) (cid:3255) (cid:3006) (cid:3265) as (cid:1842) (cid:3091) ′ = (cid:1866) (cid:3033)(cid:3045)(cid:3032)(cid:3032) /(cid:1848) , E P represents the energy at pseudogap, and (cid:1874) (cid:3029)(cid:3091) stands for the volume of a bond of -type, which can be calculated from the bond length (cid:1856) (cid:3091) of type and the number of bonds (cid:1840) (cid:3029)(cid:3100) of type per unit volume through (cid:1874) (cid:3029)(cid:3091) = ((cid:1856) (cid:3091) ) (cid:2871) / ∑ [((cid:1856) (cid:3091) ) (cid:2871) (cid:1840) (cid:3029)(cid:3092) ] (cid:3092) . For the complex multiband crystals the hardness can be calculated as a geometric average of all bond harnesses as follows [61,62]: nnvV HH /1 ])([= where n represents the number of - type bonds, which compose the real multiband crystals. The evaluated Vickers hardness in consideration of positive and reasonable populations between nearest neighbors for the new compound LaRu As is given in Table 5. The hardness value of 6.8 GPa for aRu As is very less than that of diamond (96 GPa), which is the hardest material known so far. Therefore, LaRu As is a soft material and it should be easily machineable. To visualize the nature of chemical bonding in LaRu As , the electron charge density distribution is calculated and the contour of electron charge density is presented in Fig. 3. In the adjacent scale, the blue and red colors signify the low and high electron densities, respectively. Atom with large electronegativity (electronic charge) exerts a pull on electron density towards itself [63]. Because of large difference in electronegativity and radius of atoms, the electronic charges around Ru (2.20) and As (2.18) are greater than that around La (1.10). The higher electronegativity of Ru and As exhibits strong accumulation of electronic charge while relatively low charge density indicates weak charge accumulation of La. Typical ionic picture based on standard oxidation states of atoms can help to explain the bonding nature of compounds. In LaRu As , the oxidation states of atoms are La , Ru and As . The sum of the oxidation numbers must be zero for a neutral compound such as LaRu As . In LaRu As structure, the [Ru As ] blocks are separated by La atomic sheets. The charge states of these sheets/blocks are as follows: [La] and [(Ru ) (As )2] . It can be said that a charge transfer will take place from [La] sheets to [Ru As ] blocks. The existence of ionic bonds between Ru-As within [Ru As ] blocks can be guessed from negative and positive charge balance at the positions of relevant atoms. From atom resolved density of states it is observed that the hybridization between Ru 4d and As 4p states lead to formation of covalent Ru-As bonds. Moreover, metallic type Ru-Ru bonds are expected to exist in [Ru As ] blocks due to overlapping of Ru 4d states near Fermi level [Fig. 2b]. Surprisingly, no As-As bonds exist between two adjacent [Ru As ] blocks in LaRu As . But such bonds are found to be present in (Ca/Ba)Fe As iron pnictides [64]. The possible reason for this contravention may be explained as: the low strength of Fe-As bonds in iron pnictides in comparison to Ru-As bonds in LaRu As favors to form As-As bonds in (Ca/Ba)Fe As . The above discussion reflects that the chemical bonding in LaRu As can be described as a highly anisotropic combination of ionic, covalent and metallic interactions. Fig. 3.
Charge density (left) and Fermi surface (right) of LaRu As . The Fermi surface topology has been calculated, which is presented in Fig. 3 (right). As we observe, the near-Fermi surface band pictures have a complicated ‘mixed’ character: simultaneously with quasi-flat bands and high-dispersive bands intersects the Fermi level. These features concede a multi-sheet Fermi surface consisting of two quasi-two-dimensional electron-like sheets in the corners of the Brillouin zone. The nearest sheet is also parallel to A-M direction and the distant sheet is half-
La La La La La La La La La La La La Ru Ru Ru Ru Ru Ru Ru Ru Ru Ru Ru Ru As As As As As As As As As As As As As As As As e/ Å c - a x i s old and intersected at the centre of the sheet by the Γ-M line. Between these two sheets a concave sheet appears and the Γ-M line serves as the axis of this tree-dimensional sheet. The central electron-like sheet is very complicated and has four wings shaped of open half-tube cutting along own axis. These four wings form a red-cross shape keeping each arm along Z-R direction. The existence of the multiband nature suggests that the new compound LaRu As could be a class of multiple-gap superconductor. Accurate value of electron-phonon coupling constant ( ) is important to evaluate the superconducting transition temperature T c . The method assigned to QUANTUM ESPRESSO program [65] can accurately compute the value directly. We stress that a double-delta function integration over a dense net of electron and phonon vectors ( k and q ) is required to compute the electron-phonon coupling constant with significant computational resources. Moreover, this method is more complicated for a crystal system containing many atoms in its unit cell like LaRu As . All of these prohibited us to carry on with this method. We have an alternative method in which the equation (cid:2870) (cid:1863) (cid:3003)(cid:2870) (cid:1840)((cid:1831) (cid:3007) )] provides the electron-phonon coupling constant [66]. But in the present situation we have no experimental values for electronic specific-heat coefficient and theoretically calculated value yields a lower value of [67]. So, we have only one indirect method using the McMillan’s equation: (cid:2019) = 1.04 + (cid:2020) ∗ ln (cid:3436) (cid:2016) (cid:2888) (cid:3030) (cid:3440)(1 − 0.62(cid:2020) ∗ )ln (cid:3436) (cid:2016) (cid:2888) (cid:2913) (cid:3440) − 1.04 The repulsive Coulomb potential * is assigned a value in the range of 0.10-0.15 and the choice of * is fairly arbitrary. Using known T c = 7.8 K [11] with our calculated Debye temperature D = 335 K in McMillan’s equation, we have = 0.64 and 0.76 for μ * = 0.10 and 0.15, respectively, which implies that LaRu As should be a moderately coupled BCS superconductor.
4. Conclusions
DFT methodologies were used within the first-principles pseudopotential method to calculate some physical properties of newly discovered superconductor LaRu As . The evaluated lattice constants show sound agreement with the experimental results. The obtained single crystal elastic constants obey the mechanical stability conditions for the new tetragonal system. The LaRu As crystal shows high uniaxial elastic anisotropy. The newly synthesized compound is expected to resist the shear distortion in the <110> direction than in the (100) plane. Cauchy pressure, Poisson’s ratio as well as Pugh’s ratio indicate that LaRu As should be a ductile material. The relatively small young’s modulus indicates that LaRu As should attain reasonable resistance to thermal shock. The low Debye temperature also suggests that this compound should be favorable for use as a thermal barrier coating material. The calculated electronic structures reveal that the intra-atomic bonding in new 122 phase compound LaRu As may be explained as a mixture of ionic, covalent and metallic interactions. The hardness value of 6.8 GPa implies that this material should be soft and easily machineable. The obtained Fermi surface exhibits the multi-band nature, suggesting that the new ternary compound LaRu As could be a kind of multiple-gap superconductor. The calculated indicates the strong electron-phonon coupling in LaRu As . Very recently, it came to our notice that Rahaman and Rahman has done some ab initio calculations on LaRu As and LaRu P compounds [68]. Some of the basic features of those calculations agree quite well with the ones presented in this paper. To conclude, we expect that this study will stimulate further experimental and theoretical research activities on the newly synthesized LaRu As compound. eferences [1] P. Villars, L.D. Calvert (Eds.), Pearson’s Handbook of Crystallographic Data for Intermetallic Phases, 2nd Edition, ASM International, Materials Park, OH, 1991. [2] I. B. Shameem Banu, M. Rajagopalan, Mohammed Yousuf, P. Shenbagaraman, J. Alloys Comp. 288 (1999) 88–95. [3] S. Baran, M. Bałanda, Ł. Gondek, A. Hoser, K. Nenkov, B. Penc, A. Szytuła, J. Alloys Comp. 507 (2010) 16–20. [4] M . Rotter, M . Tegel, and D . Johrendt , Phys. Rev. Lett. 101, (2008) 107006. [5] K. Sasmal, B. L. B. Lorenz, A. M. Guloy, F. Chen, Y.-Y. Xue, and C.-W. Chu, Phys. Rev. Lett. 101 (2008) 107007. [4] M . Rotter, M . Tegel, and D . Johrendt , Phys. Rev. Lett. 101, (2008) 107006. [5] K. Sasmal, B. L. B. Lorenz, A. M. Guloy, F. Chen, Y.-Y. Xue, and C.-W. Chu, Phys. Rev. Lett. 101 (2008) 107007. [6] W. Jeitschko, R. Glaum, L.Boonk, J. Solid State Chem. 69 (1987) 93–100. [7] G. P. Meisner, Physica B&C 108 (1981) 763–764. [8] I. Shirotani, T, Uchiumi, K, Ohno et al., Phys Rev B 56 (1997) 7866–7869. [9] H. Barz, H. C. Ku, G. P. Meisner et al., Proc. Natl. Acad. Sci. USA 77 (1980) 3132–3134. [10] G. P. Meisner, H. C. Ku, H. Barz, Mater. Res. Bull. 18 (1983) 983–991. [11] Q. Guo, B.-J. Pan, J. Yu, B.-B. Ruan, D.-Y. Chen, X.-C. Wang, Q.-G. Mu, G.-F. Chen, H.-A Ren, Sci. Bull. 61 (2016) 921–924. [12] M. Segall, M. Probert, C. Pickard, P. Hasnip, S. Clark, K. Refson, J. R. Yates, M. Payne, Z. Kristallogr. 220 (2005) 567–570. [13] M. C. Payne, M. P. Teter, D. C. Allan, T. A. Arias, J. D. Joannopoulos, Rev. Mod. Phys. 64 (1992) 1045–1097. [14] N. Marzari, D. Vanderbilt,M. C. Payne, Phys. Rev. Lett. 79 (1997) 1337–1340. [15] J.P. Perdew, S. Burke, M. Ernzerhof, Phys. Rev. Lett. 77 (1996) 3865. [16] F. Bloch, Z. Phys. 52 (1928) 555–600. [17] D. Vanderbilt, Phys. Rev. B 41 (1990) 7892–7895. [18] T. H. Fischer and J. Almlof, J. Phys. Chem. 96 (1992) 9768. [19] G. Kresse, J. Furthmuller, Phys. Rev. B 54 (1996) 11169–11186. [20] D. Alfe, Comp. Phys. Commun. 118 (1999) 31–33. [21]
H. J. Monkhorst, J. D. Pack, Phys. Rev. B 13 (1976) 5188. [22] I.R. Shein, A.L. Ivanovskii, Solid State Communications 149 (2009) 1860–1865. [23] R. Nath, Y. Singh, and D. C. Johnston, Phys. Rev. B 79 (2009) 174513. [24] F. D. Murnaghan, Finite Deformation of an Elastic Solid, Wiley, New York, 1951. [25] W. Voigt, Lehrbuch der Kristallphysik, Taubner, Leipzig, 1928. [26] A. Reuss, Berechnung der Fließgrenze von Mischkristallen auf Grund der Plastizitätsbedingung für Einkristalle, Z. Angew, Math. Mech . (1929) 49–58. [27] R. Hill, Proc. Phys. Soc. A
65 (1952) 349–354. [28] I. R. Shein, A. L. Ivanovskii, Intermetallics 26 (2012) 1-7. [29] M. Born, K. Huang and M. Lax, Am. J. Phys. 23 (1995) 474. [30] H. Koc, H. Ozisik, E. Deligöz, A. M. Mamedov and E. Ozbay, J. Mol. Model. 20 (2014) 2180. [31] D.G. Pettifor, Mat. Sci. Technol. 8 (1992) 345–349. [32] S.F. Pugh, Phil. Mag. 45 (1954) 823–843. [33] I. N. Frantsevich, F. F. Voronov, and S.A. Bokuta, Elastic constants and elastic moduli of metals and insulators handbook, Naukova Dumka, Kiev, 1983, pp. 60–180. [34] I.R. Shein, A.L. Ivanovskii, Physica B 410 (2013) 42–48 [35] I.R. Shein and A.L. Ivanovskii, Phys. Stat. Solidi B 248 (2011) 228–232. [36] I.R. Shein and A.L. Ivanovskii, JETP Letters, 91 (2010) 410–414. [37] M. Romero and R. Escamilla, Comp. Mat. Sci. 81 (2014) 184–190. [38] P. Ravindran, Lars Fast, P. A. Korzhavyi, B. Johansson, J. Wills, and O. Eriksson, J. Appl. Phys. 84 (1998) 4892. 39] G. Vaitheeswaran, V. Kanchana, A. Svane and A. Delin, J. Phys.: Cond. Mat. 19 (2007) 326214. [40] S. Aydin and A. Özcan, Turk. J. Phys. 34 (2010) 1 – 12. [41] B. Dong, K. Liu, X.-L. Zhou, J. Tan, X.-C. Mao and J. Chang, Phys. Stat. Solidi B, (2015); DOI 10.1002/pssb.201552545 [42] C. Li, C.Wang, F. Zhang, D. Ma, B. Wang, G.Wu, J. Phys. Chem. Solids 78 (2015) 28–34. [43] Z. Wu, E. Zhao, J. Phys. and Chem. Solids 69, (2008), 2723. [44] N. Benayad, D. Rached, R. Khenata, F. Litimein, A. H. Reshak, M. Rabah and H. Baltache, Mod. Phys. Lett. B 25 (2011) 747–761. [45] D. Varshney and S. Shriya, Phys Chem Minerals 40 (2013) 521–530. [46] J. H. Westbrook and H. Conrad, Edited Conference Proceedings, The Science of Hardness Testing and Its Research Applications, ASM, Metals Park, Ohio, 1973. [47] H.-Y. Chung, M. B. Weinberger, J. M. Yang, S. H. Tolbert, and R. B. Kaner, Appl. Phys. Lett. 92 (2008) 261904. [48] X. Wang, H. Xiang, X. Sum, J. Liu, F. Hou and Y. Zhou, J. Mater. Sci. Tech. 31 (2015) 369–374. [49] Y. Zhou, H. Xiang, X. Lu, Z. Feng, Z. Li, J. Advan. Ceram. 4 (2015) 83–93. [50] O. L. Anderson, J. Phys. Chem. Solids 24 (1963) 909. [51] S. Nagel and J. Tauc, Phys. Rev. Lett.,
M. A. Alam, M. A. Hadi, M. T. Nasir, M. Roknuzzaman, F. Parvin, M. A. K. Zilani, A. K. M. A. Islam, S. H. Naqib, J. Supercond. Nov. Magn. (2016); DOI 10.1007/s10948-016-3619-7 [54]
D. Hirai, R. Kawakami, O. Magdysyuk, R. E. Dinnebier, A. Yaresko, H. Takagi, J. Phys. Soc. Japan , 103703 (2014). [55] D. Sanchez-Portal, E. Artacho, and J. M. Soler, Solid State Commun., 95, 685 (1995). [56] R. S. Mulliken, J. Chem. Phys. , 1833 (1955). [57] W. Y. Ching and P. Rulis, Electronic Structure Methods for Complex Materials -The Orthogonalized Linear Combination of Atomic Orbitals, (Oxford Univ. Press, 2012). [58] F. M. Gao, Phys. Rev. B , 132104 (2006). [59] J. H. Westbrook and H. Conrad, The Science of Hardness Testing and Its Research Applications (ASM, Ohio, 1973). [60] H. Gou, L. Hou, J. Zhang and F. Gao, Appl. Phys. Lett. , 241901(2008). [61] A. Szymański and J. M. Szymański, Hardness Estimation of Minerals Rocks and Ceramic Materials, th edition, (Elsevier, Amsterdam, 1989). [62] V. M. Glazov and V. N. Vigdorovid, Hardness of Metalsedition, (Elsevier, Amsterdam, 1989). [62] V. M. Glazov and V. N. Vigdorovid, Hardness of Metals