Welcome to visit Communications in Theoretical Physics,
Condensed Matter Theory

Correlation between mechano-electronic features and scattering rates using deformation potential theory

  • Meryem Ziati , 1 ,
  • Sanaa Lahlali 2 ,
  • Hamid Ez-Zahraouy 1
Expand
  • 1Laboratory of Condensed Matter and Interdisciplinary Sciences, `Unité de recherche labellisée CNRST’ Faculty of Sciences, Mohammed V University, Rabat, Morocco
  • 2 Mohammed VI Polytechnic University (UM6P), BenGuerir 43150, Morocco

Received date: 2023-08-03

  Revised date: 2023-10-29

  Accepted date: 2023-11-23

  Online published: 2024-02-07

Copyright

© 2024 Institute of Theoretical Physics CAS, Chinese Physical Society and IOP Publishing

Abstract

This research paper provides valuable insight into the electronic, mechanical and transport properties of the Sr2RuO2F2 compound. The study shows that the Sr2RuO4 compound exhibits a metallic ground state and that the energy gap widens with oxygen substitution with fluorine. The concept of absolute deformation potential and its correlation with band energies and strains is explained using deformation potential theory. The paper also examines the mechanical features of Sr2RuO2F2 using the Voigt–Reuss–Hill approximation method and analyzes its elastic constants, bulk modulus and shear modulus, indicating flexibility and suitability for optoelectronic applications. The role of acoustic phonons in scattering rates and carrier mobility in Sr2RuO2F2 and its potential for phonon-mediated superconductivity is investigated. The intrinsic resistivity of electrons and holes under strain and its potential impact on superconductivity and electrical resistivity are also discussed. The insight provided by this study contributes to the current understanding of Sr2RuO2F2, and its potential applications.

Cite this article

Meryem Ziati , Sanaa Lahlali , Hamid Ez-Zahraouy . Correlation between mechano-electronic features and scattering rates using deformation potential theory[J]. Communications in Theoretical Physics, 2024 , 76(2) : 025702 . DOI: 10.1088/1572-9494/ad0f14

1. Introduction

Sr2RuO4 has drawn significant attention due to its unique electronic properties, including its unconventional superconductivity, anisotropic transport properties and non-trivial topology [13]. The underlying physics of these properties has been the focus of numerous experimental and theoretical studies [4].
Experimental studies have played a pivotal role in elucidating the electronic properties of Sr2RuO4. For example, measurements of the specific heat, resistivity and thermal conductivity have provided important information about the electronic and thermal transport properties of this material [5, 6]. In addition, studies of its magnetic behavior, such as the observation of a small magnetic moment, have led to the development of theoretical models aimed at explaining the origin of magnetism [7]. Theoretical studies based on density functional theory (DFT) have also been critical in advancing our understanding of the electronic properties of Sr2RuO4 [8]. DFT is a powerful tool for predicting the electronic structure of materials and investigating their properties. In the case of Sr2RuO4, DFT calculations have been used to predict its electronic band structure and Fermi surface, which are important for understanding its electronic properties [9]. Furthermore, DFT calculations have been used to study the effects of external perturbations, such as strain and doping, on the electronic properties of Sr2RuO4 [10].
The concept of the deformation potential theory (DPT) has emerged as a powerful tool for understanding the electronic properties of Sr2RuO4 [11]. Deformation potential is a measure of the change in energy of a band due to deformation of the crystal lattice. This theory has been used to explain the unusual anisotropic transport properties of Sr2RuO4, as well as its superconductivity [12]. Specifically, it has been suggested that the superconductivity in Sr2RuO4 arises from an odd-parity spin-triplet pairing mechanism [13], which is intimately connected to the lattice deformation. Recent developments in experimental techniques, such as scanning tunneling microscopy and angle-resolved photoemission spectroscopy, have allowed for direct measurement of the electronic properties of Sr2RuO4, providing new insight into its underlying physics [14]. For example, these techniques have revealed the presence of chiral p-wave superconductivity, which is consistent with the theoretical predictions based on the DPT [15, 16]. The combination of experimental and theoretical studies has led to significant progress in understanding the electronic properties of Sr2RuO4 [17]. The concept of the DPT has emerged as a powerful tool for explaining the unique electronic properties of this material, and its relationship with the observed superconductivity and anisotropic transport properties has been extensively investigated [18]. Further studies of Sr2RuO4 are likely to provide new insight into the physics of unconventional superconductivity and topological materials [19, 20].
This scientific paper focuses on the investigation of the uniaxial absolute deformation potentials (ADPs), which are fundamental physical properties that describe the relationship between uniaxial strain and shifts of individual band edges. These deformation potentials are critical for various applications, such as determining natural band offsets between mismatched material systems and understanding electron–phonon interactions. The paper is structured in a meticulous manner, with section 2 providing a comprehensive overview of the geometric and electronic properties using DFT. Section 3 elucidates the uniaxial deformation potentials and the associated computational methods. Section 4 discusses the Voigt–Reuss–Hill (VRH) approximation method and its use in evaluating the mechanical properties of materials. Section 5 explores the role of acoustic phonons in Sr2RuO2F2, a quaternary alloy semiconductor with potential applications in developing accurate transport property models. Section 6 examines the intrinsic resistivity of Sr2RuO2F2, a conductive material, and how it is affected by factors such as charge carrier type and mobility. The final section summarizes the key findings and provides a forward-looking perspective. Overall, this paper offers valuable insight into the uniaxial ADPs and their implications for various scientific and technological applications.

2. Geometric and electronic properties: DFT

The electronic band gap of a crystalline solid is the decisive factor in classifying it as a metal, semiconductor or insulator. Simulating electronic properties, especially for semiconductors, is crucial to understanding their physical behavior, which regulates their transport characteristics. Nonetheless, substituting techniques and strain engineering are extremely promising processes for enhancing the physical features of a target compound, expanding semiconductor applications and boosting thermoelectric devices.
The foregoing findings are intended to reflect the purpose of the present research. In fact, Ruddlesden–Popper-type superconductor Sr2RuO4 was refined as a tetragonal crystal with the space group I4/mmm (N°. 139). It consists of alternating monolayers stacked along the c-axis, as is evident from its layered crystal structure (see figure 1). As expected, based on static calculations, the relaxed state parameters attained using the generalized gradient approximation (a = b =3.84 Å; c = 12.78 Å) are congruent with earlier theoretical reports as well as experimental findings [2123]. The electronic properties of Sr2RuO4 along the high-symmetry direction were investigated through the equilibrium lattice parameters via the Yukawa screened hybrid functional and spin–orbit (SO) coupling. Formally, the SO effect is still prevalent and yields adjustments to the total energy and its derivatives. As inner-shell electrons are drawn closer to the nucleus, their kinetic energy rises, and relativistic effects become crucial. In fact, the strength of the SO coupling increases quickly with the atomic number Z. For light elements, it is frequently possible to ignore or approximate using the scalar relativistic parts of the Dirac equation [24]. Nevertheless, even with only light elements, the SO interaction can effectively describe some physical properties [25]. In the second row, for heavier elements and transition metals, to some extent, the SO effect becomes important for structural and dynamical features. Meanwhile, Yukawa screening has been shown to be more successful in predicting band gaps of materials due to its screened Coulomb exchange interaction [26]. This is the inherent reason that SO coupling and hybrid functional were applied in the present scientific report. As illustrated in table 1, our outcomes notably exhibit a zero gap, and the combined contributions of the YS-PBE0 approach and SO coupling suggest a metallic ground state of Sr2RuO4. These results are fully consistent with other theoretical discoveries [27].
Figure 1. Schematic view of Sr2RuO2F2 structure.
Table 1. Calculated equilibrium parameters: lattice constants a (Å), c (Å) and atomic volume V (Å) 3.
Sr2RuO4 Sr2RuO2F2
V (Å)3 189.45 193.85
a (Å) 3.84 3.98
b/a 1 1
c (Å) 12.78 13.23
Eg YS-PBE0 0 0
YS-PBE0 + SOC 1.586 1.518
Extrinsic defects (i.e. doping, substituting, etc) lead to a local disturbance of the crystal lattice periodicity. A defect is described by its electronic configuration (i.e. by the wave functions ψi associated with bound electrons i and also by the energy levels Ei of these electrons) and atomic configuration, including the distortions of the surrounding crystal lattice. The study of point defects is usually paramount. The presence of imperfections is typically regarded as an efficient means to enhance the physical features of a target compound, especially the electronic and transport properties. The electronic behavior can be deduced, at least in principle, from the electronic structure (ψi, Ei), and this was the aim of our research. From the relevant literature, we highlighted the studied compound. The fractional atomic coordinates of the two non-equivalent oxygen atoms were defined as follows: O (1): 4c (0, 0.5, 0) and O (2): 4e (0, 0, 0.1607) [28]. The strontium and ruthenium reside in Sr: 4e (0, 0, 0.3527) and Ru: 2a (0, 0, 0), respectively. Our findings suggest that the least total energy was produced when substituting oxygen atoms with fluorine atoms at O (1) sites (i.e. creating point defects), thus being energetically the most advantageous. Volume–energy optimization was performed to ascertain the structural stability of unconstrained Sr2RuO2F2 using Murnaghan’s equation of state [29]. Under free conditions, the tetragonal symmetry of the primitive unit cell persists despite the embedding of fluorine atoms. Somehow, the proposed substitution pushes the bands and electronic states from the Fermi level to lower energy levels, creating a forbidden region, where the density vanishes, as depicted in figure 2.
Figure 2. Calculated density of states and electronic band gaps of Sr2RuO2F2 with the YS-PBE0 approach.
The electronic relocation was systematic, due to the electronic amendment and enhancement of interaction intensities between atoms after doping. This led to the widening of the energy gap in both the majority and minority spins. The upper part of the valence band and the bottommost of the conduction band are inclined at different symmetry points, attributed to an indirect band gap of the material (see figure 3).
Figure 3. Calculated band structure of Sr2RuO2F2 with and without spin–orbit coupling (SOC).
For instance, several physical quantities, notably the band locations, were strongly affected by the choice of functional. Indeed, the band gap for Sr2RuO2F2 at the YS-PBE0 level of theory (1.586 eV) showed a slight mismatch when SOC was applied (1.518 eV), as quantified in table 1. Since ruthenium (Ru) is a heavier element than strontium (Sr) and oxygen (O), it was expected that Ru would exhibit stronger SOC. However, the incorporation of SOC terms had no discernible effect on the primary electronic features, such as the electronic band gap. This proves that ruthenium is not ‘heavy enough’ to induce qualitative changes. Nonetheless, we have taken full SOC effects into account when calculating the band structures and electronic band gaps. Under free conditions, the SO separation energy was exceptionally small, as illustrated in figure 4. We attribute this smaller value to the minor effect of the nuclear charge (Zeff) of the Ru atom and the weak variation of the charge gradient in an elementary crystal system such as Sr2RuO2F2. Nonetheless, the SOC effect is greatest in the conduction band due to the orbital composition of these bands, which is typically contributed by heavy atoms (such as ruthenium). In the following, strain engineering was quantitatively applied along the c-axis, which significantly affects the electronic features in Sr2RuO2F2. For convenience, we replicated a compressed and expanded region of the same material, which is clearly visible in figure 3.
Figure 4. Change in spin splitting at conduction and valence bands under strain.
The semiconducting behavior persists, and so a linear relationship was observed between the external strain and the corresponding bandgap values (see figure 2). Stretching provides effective tunability in the band gap and returns it to 1.75 eV at ϵz = + 0.08. The bandgap energy folds back to 1.15 eV at ϵz =0.08 upon spreading the compressive strain. The outcomes obtained from the YS−PBE0 method were fairly homogeneous and proportional to those obtained from YS−PBE0 + SOC. The effective electron mass ratio and the effective masses of holes can be readily calculated from the band structure using the formula ${m}^{* }={\overline{h}}^{2}{\left({\partial }^{2}\,E\left(K\right)/\partial {K}^{2}\right)}^{-1}$ [29]. In the equilibrium state (i.e. unstrained Sr2RuO2F2), the topmost valence bands exhibit heavy masses, which are detrimental to their mobility. The mass of the upper valence band is as low as ≃9.10−31Kg and falls continuously in the compressed region, as shown in figure 5. In contrast, the pressure pushes the electron effective masses upwards to a value of ≃ 1.75.10−31Kg and a huge increase was witnessed. Neither experimental nor theoretical findings are currently available for comparison. Nonetheless this discovery yields a theoretical prediction, which ought to fill a gap and contribute to new research.
Figure 5. Effective masses at room temperature of strained and unconstrained Sr2RuO2F2.

3. Uniaxial ADPs: DPT

In an attempt to explain the interaction between electrons and acoustic phonons and to circumvent this problem, Bardeen and Shockley [30] introduced the DPT [31]. The ADP concept reflects the electron–phonon interaction and outlines the mismatch between electron energy levels in the same stretched and compressed solid. The former concept was coined to describe the strain-induced shift of one-electron energy levels in the structure [32]. Indeed, conduction and valence band displacements caused by the lattice deformations due to the sonic waves significantly affect carrier mobility. Carrier mobilities of their lathe were sketched in terms of scattering by the deformation potential. Research demonstrated that the ADPs are simple derivatives of the energies with respect to the strains by making straightforward assumptions [33] and can be defined as,
$\begin{eqnarray}{D}_{i}=\frac{{\rm{\Delta }}{E}_{i}}{{\rm{\Delta }}{\rm{l}}{\rm{n}}V},\end{eqnarray}$
where ${D}_{{i}}$ is the deformation potential of a given state i, ${\rm{\Delta }}{E}_{{i}}$ is the shift of its energy under strain and usually Ei depicts the eigenvalue of the Kohn–Sham Hamiltonian. The relative change of volume is related to the normal strains in x, y and z directions and is given by [33],
$\begin{eqnarray}{\rm{\Delta }}\,{\rm{l}}{\rm{n}}V=\frac{{\rm{\Delta }}V}{V}={\varepsilon }_{xx}+{\varepsilon }_{yy}+\varepsilon .\end{eqnarray}$
For completeness, ADP is often derived from strained superlattice computations on the Sr2RuO2F2 3D semiconductor. Theoretical ascertainment of the ADPs has been fairly controversial and involves a substantial number of inspections, renditions and assumptions [34]. The uniaxial ADP highlights the long-range electrostatic effects and can be described by [33],
$\begin{eqnarray}{D}_{zz}^{i}=\frac{{\rm{\Delta }}{E}_{i}}{{\varepsilon }_{z}}.\end{eqnarray}$
Quantitatively, in this research, ϵz represents the external uniaxial strain applied along the c-axis. The conduction and valence band ADPs (${{\rm{D}}}^{{\rm{VBM}}}\,{\rm{and}}\,{{\rm{D}}}^{{\rm{CBM}}}$) allow for the direct calculation of the bandgap deformation potential [33]:
$\begin{eqnarray}{D}_{zz}^{\mathrm{Gap}}=\,\frac{{\rm{\Delta }}{E}_{\mathrm{Gap}}}{{\rm{\Delta }}\,{\rm{l}}{\rm{n}}V}={D}^{\mathrm{CBM}}-{D}^{\mathrm{VBM}}.\end{eqnarray}$
The fundamental derivatives of the band energies with respect to strains at the ground states are ADPs. The latter can only be used to predict the absolute energies of the electronic states for different strain conditions. For convenience, we will later demonstrate the correlation between band energies and strains by illustrating the linear regime of this dependence in the intrinsic case of Sr2RuO2F2. To begin, we investigate the ADPs of uniaxially strained Sr2RuO2F2 using the YS-PBE0 approach and SO coupling. The energies of the topmost valence band at the Λ, N and P points and the bottommost of the conduction band at the H, N and Γ positions in the Brillouin zone, as a function of strain, are listed in table 2.
Table 2. ADPs (in eV) of the VBM at the $\bigwedge $, N and P, and of the CBM at the H, N and Γ points.
Approach Uniaxial absolute deformation potentials
${\boldsymbol{AD}}{{\boldsymbol{P}}}_{\wedge }^{{\bf{VBM}}}$ ${\boldsymbol{AD}}{{\boldsymbol{P}}}_{{\boldsymbol{N}}}^{{\bf{VBM}}}$ ${\boldsymbol{AD}}{{\boldsymbol{P}}}_{{\boldsymbol{P}}}^{{\rm{VBM}}}$ ${\boldsymbol{AD}}{{\boldsymbol{P}}}_{{\boldsymbol{H}}}^{{\bf{CBM}}}$ ${\boldsymbol{AD}}{{\boldsymbol{P}}}_{{\boldsymbol{N}}}^{{\bf{CBM}}}$ ${\boldsymbol{AD}}{{\boldsymbol{P}}}_{{\boldsymbol{\Gamma }}}^{{\bf{CBM}}}$
YS-PBE0 1.027 −0.863 1.141 4.250 3.278 4.242
YS-PBE0+SOC 0.861 −0.432 0.857 4.298 3.201 4.213
The ADPs were obtained by linearly extrapolating the band energy edges as a function of axial strains, as shown in figure 6. The findings indicate that the correlation can be assumed to be linear for strains between ϵz = −0.08 and ϵz = + 0.08.
Figure 6. Energies of the VBM and CBM of Sr2RuO2F2 as a function of the uniaxial strain in the z direction. Solid and dotted lines are linear fits using values calculated for strains of +8% and −8%.
The band edge positions of Sr2RuO2F2 are highly responsive to strain and evolve towards different values in a consistent manner. The valence band energies (VBM) at the N point steadily decrease in the stretched region (i.e. when the lattice parameter grows). Despite this behavior, the CBM energies at the N point progressively increase, and this trend persists even when external strain is applied (see figure 6).
Note that the change in valence and conduction states under deformation can be due to the following physical trend.
i

(i) The kinetic energy effects push the VBM energies upwards under strain at the P and $\bigwedge $ points. It seems that the kinetic energy of our system is commensurate with the reciprocal lattice vector $\vec{k}\,$ for which the related plane waves ψk have the Bravais lattice periodicity. For example, if we change our system to a different material or change particles within it, then the system will vibrate at different frequencies until a new set of particle positions is found in which the same sort of periodic motion occurs [34, 35].

ii

(ii) When discussing the effect of bond length on the energy level of a state, it is important to consider whether the state has bonding (π) or antibonding (π*) nature. Bonding states have a lower energy level as the bond length increases, while antibonding states have a higher energy level. This is because bonding states have electrons that occupy the same region of space, resulting in a lower energy level due to the mutual attraction of electrons. On the other hand, antibonding states have electrons that occupy different regions of space, resulting in a higher energy level due to the mutual repulsion of electrons [34, 35].

Regarding the quaternary alloy Sr2RuO2F2, the state at the N point of the valence band (VBM) exhibits metal-like characteristics and has a bonding interaction with the F and O orbitals. As the alloy is expanded, the energy of this state experiences a significant decrease using both the YS-PBE0 approach and SO coupling. This is due to the fact that as the distance between atoms increases, the bond length (i.e. stretching) also increases, leading to a decrease in energy. The bonding nature of the state contributes to an increase in energy when the bond length shortens, but as the bond length increases, the energy also decreases. (see figure 6). The kinetic energy of electrons within a substance is intimately related to the motion of these subatomic particles. As the distance between atoms within a compound expands (i.e. bond length under stretch), the energy of these electrons correspondingly decreases due to the greater spatial separation between atoms. This reduction in energy is a direct consequence of the decline in the momentum of the electrons as the atoms are farther apart, and the opportunities decline in the momentum of the electrons as the atoms are farther apart, and the opportunities for interaction with the nuclei are accordingly diminished. The correlation between the potential energy of electrons in a material and the electrostatic interactions between these electrons and nuclei is a complex and fascinating phenomenon. As the distance between nuclei, otherwise known as bond length, increases, the potential energy of the electrons correspondingly rises as a result of the weakening of the electrostatic interactions. This highlights the intricate interplay between the various components of Sr2RuO2F2 quaternary alloy and their effect on the overall energy state. The kinetic energy of the electrons decreases at a greater rate than the increase in potential energy, resulting in a net diminution of the aggregate energy of the system. This salient phenomenon, intrinsic to the bonding state at the valence band maximum (VBM), is the very explanation for the diminutive ADP observed in table 2. In general terms, the effect of antibonding interactions on the energy edges of a material can be significant, particularly in the conduction band (CBM). As mentioned earlier, the antibonding states have a lower electron density between the nuclei, which leads to higher repulsion and increased energy. This elevation of energy can propel the state into the conduction band, thereby rendering it amenable for electrical conduction [36]. In Sr2RuO2F2 semiconductor, the antibonding interactions arise from the overlap of the electronic wave functions of the Ru-4d orbitals in the RuO2 planes with the F-2p orbitals in bonds and the electrons of the F layers. The antibonding interactions are formed between the electrons of the Ru–O and the Ru–F bonds. The Born–Oppenheimer approximation is a fundamental principle in quantum mechanics that explains the relationship between the positions of nuclei and electrons in a substance, and the resulting energy levels [37]. According to this principle, when the bond length lengthens (e.g. from ϵz = −0.08 to ϵz = + 0.08), the energy of the antibonding states at N, H and Γ points also rises, as shown in figure 6. This is because the electrons in these antibonding states, which have higher energy than the bonding states, are farther apart, resulting in higher potential energy. As illustrated in figure 6, the conduction band (CBM) electrons are more sensitive to external strains/deformations, and their energy levels change more readily in response to these changes. The energy edge is significantly affected by strain; it increased steadily under pressure/stretch and abruptly jumped to ≈3.77 eV, ≈1.87 eV and ≈1.85 eV at N, H, and Γ points respectively, when ϵz = + 0.08 (see figure 6). This results in the ADP being larger in the conduction band than in the valence band, as listed in table 2. The previously observed behavior persists even upon applying SO coupling.
The bandgap deformation potential, a fundamental characteristic in the assessment of the mechanical attributes of semiconductors, serves to quantitatively evaluate the alteration in the band gap of a substance in response to an imposed strain or deformation. To calculate the bandgap deformation potential of the Sr2RuO2F2 quaternary alloy, a linear fit is performed on the electronic band gaps at various strains, as illustrated in figure 7. This process yields a correlation between the electronic band gap and strain, thereby enabling the computation of the bandgap deformation potential, listed in table 3. The upsurge in the gaps (${{E}}_{{\rm{gap}}}^{{N}-{\rm{\Gamma }}},$ ${{E}}_{{\rm{gap}}}^{\wedge -{H}}\,{\rm{and}}\,{{E}}_{{\rm{gap}}}^{{N}-{N}}$) under pressure/stretch is a consequence of the alteration in the electronic structure of Sr2RuO2F2 brought about by the modification in the interatomic spacing. Specifically, when the Sr2RuO2F2 quaternary alloy is stretched by an amount ϵz, the bond lengths grow, which leads to a shrinkage in the overlap between the valence and conduction bands. This reduction in overlap results in an enlargement of the band gap, perceived in figure 7, which is reflected by an escalation in the gap deformation potential (${{\rm{DP}}}_{{\rm{gap}}}$). When comparing the results obtained using the YS—PBE0 and YS—PBE0 + SO functionals, it is important to note that the two methods predict similar gap deformation potentials for the majority of cases. The significance of the gap deformation potential as a metric of the alteration in the band gap, consequent to mechanical strain, cannot be overstated. In order to comprehend the mechanical characteristics of semiconductors, it is imperative to establish a correlation between the strain and the band gap. Thus, the gap deformation potential serves as a crucial indicator in this regard.
Figure 7. Band gap deformation potentials. Band energies result from a linear extrapolation of the values calculated at − 8% and + 8% strain.
Table 3. ADPs (in eV)${\,{DP}}_{{\rm{gap}}}^{{N}-{\rm{\Gamma }}},$ ${{DP}}_{{\rm{gap}}}^{\wedge -{H}}$ [indirect gap], ${{DP}}_{{\rm{gap}}}^{{N}-{N}}$ [direct gap].
Approach Band gap deformation potentials
${\boldsymbol{D}}{{\boldsymbol{P}}}_{{\bf{gap}}}^{{\boldsymbol{N}}-{\boldsymbol{\Gamma }}}$ ${\boldsymbol{D}}{{\boldsymbol{P}}}_{{\bf{gap}}}^{{\boldsymbol{N}}-{\boldsymbol{N}}}$ ${\boldsymbol{D}}{{\boldsymbol{P}}}_{{\bf{gap}}}^{\wedge -{\boldsymbol{H}}}$
YS-PBE0 4.791 3.871 3.991
YS-PBE0+SOC 4.739 3.726 4.049

4. Mechanical features: VRH approximation

The VRH approximation method is a widely employed technique within the field of materials science for ascertaining the mechanical properties of various materials, such as their stiffness and strength [38]. This method leverages a weighted average of the results obtained through the harnessing of the Voigt, Reuss and Hill approximations, which are based on distinct assumptions regarding the microstructure of the material in question. The implementation of this method allows for a more precise estimation of elastic properties since it takes into account the inherent heterogeneity of materials at a microstructural level, resulting in a comprehensive evaluation of the material’s mechanical features. The Sr2RuO2F2 quaternary alloy, with a tetragonal crystal structure, is the subject of our extensive research regarding its mechanical properties. The stiffness matrix will be defined by the nine independent elastic constants specific to that material and crystal structure [39, 40]:
$\begin{eqnarray}{C}_{{\rm{t}}{\rm{e}}{\rm{t}}{\rm{r}}{\rm{a}}}=\,\left(\begin{array}{c}{C}_{11}\\ .\\ .\\ 0\\ 0\\ 0\end{array}\,\begin{array}{c}{C}_{12}\\ {C}_{22}\\ .\\ 0\\ 0\\ 0\end{array}\,\begin{array}{c}{C}_{13}\\ {C}_{23}\\ {C}_{33}\\ 0\\ 0\\ 0\end{array}\,\begin{array}{c}0\\ 0\\ 0\\ {C}_{44}\\ 0\\ 0\end{array}\,\begin{array}{c}0\\ 0\\ 0\\ 0\\ {C}_{55}\\ 0\end{array}\,\begin{array}{c}0\\ 0\\ 0\\ 0\\ 0\\ {C}_{66}\end{array}\right).\end{eqnarray}$
Pursuant to the Born–Huang stability criteria, which are both necessary and sufficient for a tetragonal system, a structure that is mechanically sound must fulfill the following requisites [39, 40]:
$\begin{eqnarray}\begin{array}{l}\left({C}_{11}-{C}_{12}\right)\gt 0,\left({C}_{11}+{C}_{33}-2{C}_{13}\right)\gt 0;\\ {C}_{11}\gt 0,{C}_{22}\gt 0,{C}_{44}\gt 0,{C}_{66}\gt 0;\\ \left(2\,{C}_{11}+{C}_{33}+2\,{C}_{12}+4{C}_{13}\right)\gt 0.\end{array}\end{eqnarray}$
The elastic constants describe the material’s ability to withstand deformation when a load is applied. The assessed elastic constants for the strained Sr2RuO2F2 perovskite reveal that it is mechanically stable under strain. This can be confirmed by appraising the elastic constants against the Born criteria, as previously mentioned. The findings in table 4 reveal that all the elastic constants of Sr2RuO2F2 satisfy these criteria. From the elastic constants, we can also derive various moduli that give insight into the mechanical properties of the material. The bulk modulus (B) sets forth the material’s resistance to changes in volume under pressure and can be deduced from B = (C11 + 2C12)/3 [40].
Table 4. Calculated elastic constants for Sr2RuO2F2 with YS-PBEO and YS-PBEO + SOC.
Strain YS-PBE0 YS-PBE0 + SOC
C11 C12 C33 C66 C11 C12 C33 C66
−0.08 722.4 303.8 1486.5 21.1 721.2 1481.4 300.7 21.3
−0.06 471.1 190.2 1330.3 15.9 648.4 1326.8 189.1 16.1
−0.04 662.9 355.8 764.7 13.9 663.0 764.9 353.7 13.9
−0.02 607.8 328.7 594.9 11.7 616.1 584.9 359.2 11.3
0.00 304.6 59.7 546.4 4.5 559.6 542.9 308.8 4.9
0.02 562.3 332.5 575.0 8.1 562.7 531.6 325.7 9.1
0.04 565.4 352.5 512.1 2.6 572.2 521.9 323.1 4.5
0.06 560.3 338.9 398.8 4.3 573.6 399.4 323.1 4.5
0.08 484.6 263.4 291.8 3.7 618.8 281.3 361.2 3.7
The calculated values of strained Sr2RuO2F2 span from BV = 288.34 GPa (ϵz = + 0.08) to BV = 450.93 GPa (ϵz = −0.08) through YS-PBE0. Young’s modulus (E) portrays the material’s resistance to longitudinal stress and can be determined by E = 9BG/(3B + G) [40]. The large E suggests that the material is very stiff and resistant to deformation when a force is applied. Values ranging from EV = 217.29 GPa (ϵz = + 0.08) to EV = 490.02 GPa (ϵz = −0.08) are acquired for strained Sr2RuO2F2 through YS-PBE0. The Poisson’s ratio (ν) describes the material’s tendency to expand or contract in the transverse direction when under longitudinal stress, and it can be computed with the relation ν = (3B –2G)/ (2(3B+G)) [40]. It extends from ${\nu }_{V}$ = 0.32 (ϵz = −0.08) to ${{\rm{\nu }}}_{{\rm{V}}}$ = 0.37 (ϵz = + 0.08) through YS-PBE0. The shear modulus (G) states the material’s resistance to shearing forces and can be extrapolated from G = C44 [40]. G ranges from 79.05 GPa (ϵz = + 0.08) to 185.77 GPa (ϵz = −0.08) through YS-PBE0, which is related to Young’s modulus and Poisson’s ratio by G = E/2(1 + ν) [40]. The low bulk modulus and shear modulus yielded under stretch by an amount ϵz = + 0.08 indicate that strained Sr2RuO2F2 is flexible, making it suitable for optoelectronic applications such as solar cells. The Voigt, Reuss and Hill models were all applied to estimate the elastic properties of the Sr2RuO2F2 quaternary alloy, but they each make different assumptions about the behavior of the researched compound. These assumptions are based on the material’s structure and the way that it deforms under strain. For instance, the Voigt model assumes that the material is isotropic (i.e. has the same properties in all directions) and elastic (i.e. returns to its original shape after being deformed), while the Reuss model assumes that the material is isotropic and plastic (i.e. permanently deforms under strain). Since all three models make different assumptions, they yield slightly different results. However, these results are still considered to be close to each other, as listed in tables 5 and 6.
Table 5. Calculated Bulk modulus and shear modulus for Sr2RuO2F2 with YS-PBE0 and YS-PBE0 + SOC.
Strain YS-PBE0 YS-PBE0 + SOC
Bulk modulus Shear modulus Bulk modulus Shear modulus
BV BR BH GV GR GH BV BR BH GV GR GH
−0.08 450.93 428.31 439.62 185.77 57.89 121.89 450.39 428.31 439.62 185.77 57.89 121.83
−0.06 418.20 380.57 399.38 163.17 48.24 105.71 418.20 380.57 399.38 163.17 48.24 105.71
−0.04 387.33 386.57 386.57 117.10 42.05 79.58 387.33 386.57 386.95 117.10 42.05 79.58
−0.02 358.59 347.23 352.91 90.89 35.39 63.14 358.59 347.23 352.91 90.89 35.39 63.14
0.00 166.55 156.33 161.44 84.01 17.81 50.91 333.11 329.21 331.16 48.89 18.93 51.91
0.02 310.51 305.28 307.89 95.47 28.94 62.20 310.51 305.28 307.89 95.47 28.94 62.20
0.04 291.12 277.96 284.54 95.29 14.76 55.02 387.33 389.57 386.95 117.10 42.05 79.58
0.06 273.28 244.43 258.86 88.77 17.25 53.01 418.20 380.57 399.38 163.17 48.24 105.71
0.08 288.34 218.71 253.52 79.05 14.49 46.77 450.93 428.31 439.62 185.77 57.89 121.83
Table 6. Calculated Young's modulus and Poisson ratio for Sr2RuO2F2 with YS-PBE0 and YS-PBE0 + SOC.
Strain YS-PBE0 YS-PBE0 + SOC
Young's modulus Poisson ratio Young’s modulus Poisson ratio
EV ER EH ${{\boldsymbol{\nu }}}_{{\boldsymbol{V}}}$ ${{\boldsymbol{\nu }}}_{{\boldsymbol{R}}}$ ${{\boldsymbol{\nu }}}_{{\boldsymbol{H}}}$ EV ER EH ${{\boldsymbol{\nu }}}_{{\boldsymbol{V}}}$ ${{\boldsymbol{\nu }}}_{{\boldsymbol{R}}}$ ${{\boldsymbol{\nu }}}_{{\boldsymbol{H}}}$
−0.08 490.02 166.17 334.58 0.32 0.44 0.37 490.02 166.17 334.17 0.32 0.44 0.37
−0.06 433.18 138.86 291.41 0.33 0.44 0.38 433.18 138.86 291.41 0.33 0.44 0.38
−0.04 319.15 121.75 223.42 0.36 0.45 0.40 319.15 121.75 223.42 0.36 0.45 0.40
−0.02 251.44 102.69 178.77 0.38 0.45 0.42 251.44 102.69 178.77 0.38 0.45 0.42
0.00 215.76 215.76 138.21 0.28 0.45 0.36 234.73 55.73 148 0.38 0.47 0.43
0.02 259.78 84.15 174.83 0.36 0.45 0.41 259.78 84.15 174.83 0.36 0.45 0.41
0.04 257.76 43.50 155.08 0.35 0.47 0.41 319.15 121.74 223.42 0.36 0.45 0.40
0.06 240.29 50.57 148.87 0.35 0.47 0.40 433.18 138.86 291.41 0.33 0.44 0.38
0.08 217.29 42.54 132.18 0.37 0.47 0.41 490.02 166.14 334.58 0.32 0.44 0.37
Ductile and brittle materials possess distinct features and are employed in various applications [41]. Ductile materials possess the capability to undergo plastic deformation prior to fracture, whereas brittle materials tend to exhibit sudden breakage with minimal or no plastic deformation. The mechanical properties of the strained Sr2RuO2F2 quaternary alloys were evaluated employing the Pugh criteria, as outlined in table 7. This methodology entails the calculation of the Pugh ratio, obtained by dividing the bulk modulus by the shear modulus (BH/GH) [40]. Materials exhibiting a high BH/GH ratio are considered to possess ductile characteristics, whereas those with a low BH/GH ratio are deemed brittle. Our study reveals that the strained compounds exhibited ductile behavior, as evidenced by their BH/GH ratios, which were typically greater than 1.75. It is worth noting that these materials (i.e. between ϵz = −0.08 and ϵz = + 0.08 including the ground state) are relatively inexpensive and possess high strength and toughness. Our research has led to the prediction of ductile materials that exhibit exceptional strength and toughness features, making them suitable for various high-demand applications such as those found in the construction, transportation and energy sectors [42, 43]. The Cauchy relationship, also known as the Cauchy–Born rule, is a mathematical relationship that describes the deformation of a material in response to an applied strain. The relationship is represented by the equation Cc = C12C44, where Cc is the compliance tensor, C12 is the elastic compliance tensor and C44 is the shear compliance tensor [39]. C12 represents the Sr2RuO2F2 ability to deform in response to a compressive or tensile strain, while C44 represents its ability to deform in response to a shear strain. The difference between these two values, C12C44, is a measure of the material’s ductility or brittleness. Our study on the mechanical response of strained Sr2RuO2F2 quaternary alloys via the Cauchy relationship has revealed positive values of C12C44 (see table 7), a clear indication of ductility.
Table 7. Cauchy and Pugh ratio for Sr2RuO2F2 with YS-PBE0 and YS-PBE0 + SOC.
Strain YS-PBE0 YS-PBE0 + SOC
Cauchy Pugh Cauchy Pugh
− 0.08 242.20 3.52 239.10 3.61
− 0.06 132.30 3.64 131.50 3.77
− 0.04 302.30 4.86 300.10 5.02
− 0.02 134.40 5.43 309.70 5.58
0.00 51.90 6.02 265.70 6.37
+ 0.02 293.50 4.68 289.20 4.94
+ 0.04 304.90 5.46 296.7 5.17
+ 0.06 301.30 5.00 285.10 4.88
+ 0.08 231.20 4.35 328.9 6.89
These results indicate that these alloys are capable of undergoing significant deformation prior to failure under applied strain. These findings align with the conclusions derived from the Pugh criteria, further supporting the assertion that Sr2RuO2F2 quaternary alloys possess highly desirable ductile characteristics. As such, these materials are well-suited to a broad range of engineering applications. It is worth noting that a material exhibiting a negative value of C12C44 is considered to be brittle since it will fracture prior to undergoing substantial deformation under applied pressure/stretch [40, 41].

5. Scattering rates and carrier mobility of Sr2RuO2F2

Acoustic phonons are a fundamental aspect of the physics of crystalline materials, and they correspond to the quantum mechanical excitations of the lattice vibrations [44]. In the case of the Sr2RuO2F2 compound, the crystal lattice is composed of Ru, F and O atoms that are arranged in a layered structure. These phonon excitations can propagate through the lattice in a variety of directions, including along the c-axis, which is perpendicular to the layers (i.e. perpendicular direction (0 0 1)), as well as within the layers themselves (i.e. parallel direction (1 1 0)). The study of acoustic phonons in this material can provide valuable insight into its properties and behavior, and further our understanding of the underlying physics of crystalline materials. As the phonon wave propagates through Sr2RuO2F2, the local lattice can undergo compression and expansion. This is due to the fact that the phonon wave is a collective motion of the atoms in the lattice, where the atoms oscillate about their equilibrium positions in phase with each other. The pressure and stretch of the local lattice leads to a perturbation of the potential energy of the bands in the crystal, and this perturbation is related to the deformation potential (e.g. in this case, the change in the potential energy of the electrons can be calculated by multiplying the deformation potential constants with the corresponding strain). The perturbation of the potential energy causes the scattering of carriers, such as electrons and holes. The scattering rate depends on the overlap of the initial and final states of the carriers, which is determined by the phonon wave’s momentum and polarization [45]. The electronic properties of the quaternary alloy Sr2RuO2F2, such as electrical resistivity and mobility, can be affected by the scattering of carriers. This has been shown in subsequent sections of our research. Complementary, the phonon–electron interaction in Sr2RuO2F2 can also cause the phonon-mediated superconductivity [46]. Under these circumstances, the phonons mediate the attraction between the electrons, and can lead to the formation of Cooper pairs, which are the fundamental building blocks of superconductivity [47]. The scattering rate is consistently the rate at which carriers scatter (or interact) with each other. The relaxation time is the amount of time it takes for a system to return to equilibrium after being perturbed. The two are related in that the relaxation time is inversely proportional to the scattering rate. It is arguable that if the relaxation time is known, the scattering rate can be calculated by taking the reciprocal of the relaxation time [48]. The relaxation time, represented by τ, is a crucial parameter that sheds light on the energy dissipation dynamics of carriers in semiconductor materials. The inverse proportionality between relaxation time and carrier mobility highlights the intrinsic relationship between the ease of carrier movement and the rate of energy dissipation. The order of magnitude of relaxation time in the nanosecond range for electrons and holes in semiconductors underscores its significance in determining the electrical and transport features of these materials.
The results of this study provide strong evidence for the feasibility of using the Sr2RuO2F2 semiconductor in the development of a highly accurate model for transport properties. To confirm this potential, we conducted a thorough evaluation of the electronic relaxation time in both the (0 0 1) and (1 1 0) directions, by means of mathematically rigorous formulas to ensure the validity and precision of our findings [49]:
$\begin{eqnarray}\tau =\,\frac{{\left(8\pi \right)}^{1/2}\,{ \hbar }^{4}\,{C}_{ij}}{3\,{({m}^{* }{k}_{{\beta }}T)}^{\displaystyle \frac{3}{2}}\,{E}^{2}}.\end{eqnarray}$
The reduced Planck constant is represented by ħ, e is the elementary charge and Cij is the elastic matrix constant obtained by VRH approximation. For the (0 0 1) direction (⊥), C33 is selected when calculating carrier mobility and electronic relaxation time. Forthe  (1 1 0) direction (∥), C66 is employed. kβ is the Boltzmann constant; T is the temperature (set at 300 K) and E is the deformation potential constant, ascertained by the band edge shift with respect to the strain as previously performed. Relaxation times of carriers are illustrated at 300 K in figure 8 [(a) perpendicular direction and (c) parallel direction]. The relaxation time of holes (i.e. positively charged carriers) and electrons (i.e. negatively charged carriers), which pertains to the duration it takes for these carriers to relinquish their energy and attain a state of equilibrium within a material, was observed to be significantly longer in the perpendicular direction (${\tau }_{\perp }$) as opposed to the parallel direction (${\tau }_{\parallel }$) of Sr2RuO2F2 at room temperature (ca. 300 K). This phenomenon can be attributed to the intricate interplay between the material’s electronic structure and the interaction of the holes and electrons with the crystal lattice. Being a layered compound, Sr2RuO2F2 exhibits an electronic confinement of the holes and electrons primarily to the RuF layers. The longer relaxation time in the perpendicular direction implies that the carriers experience a greater difficulty in scattering and losing their energy as they traverse through the layers, as opposed to within the layers. Intuitively, these effects can be attributed to the presence of strong electronic correlations within the material, which enhance the stability of the electronic states and reduce the carrier’s ability to scatter and lose energy. In reality, the unique crystal structure of Sr2RuO2F2 could also play a role, as the presence of specific types of bonding or lattice defects may impact the carrier’s mobility and scattering.
Figure 8. Relaxation time: [(a) perpendicular direction and (c) parallel direction] using YS-PBE0.
The modulation of the relaxation time in Sr2RuO2F2 under external strain is the result of its intricate crystal structure and complex inter-atomic interactions (i.e. covalent bonding, electrostatic forces and electron transfer) between Ru-d, F-p, Sr-p and O-p orbitals.
The Sr2RuO2F2 quaternary alloy is defined by the presence of ruthenium ions occupying a Ru+4 oxidation state, which participate in forming strong covalent bonds with the O−2 oxidation state oxygen ions. These covalent bonds lead to the establishment of a robust crystal lattice. Fluorine, as a highly electronegative element, also plays a crucial role in determining inter-atomic interactions through its participation in electrostatic bonds. The Sr atoms provide stability to the lattice through the formation of ionic bonds with the oxygen ions. Strain induces changes in lattice dimensions, leading to alterations in the magnitude of these inter-atomic interactions and phonon scattering processes.
Upon tensile strain, the lattice expands, resulting in a reduction of scattering and a consequent acceleration of the relaxation rate, shortening the relaxation time. The relaxation periods of electrons and holes in the taut zone of (001) and (110) Sr2RuO2F2 shorten as strain magnitudes alter. In the case of (001) Sr2RuO2F2, the electron relaxation time fell from 70 to 41 fs and hole relaxation time underwent a reduction from 8 to 2.5 fs as ϵz progressed from ϵz = + 0.00 to ϵz = + 0.08. Analogously, in the instance of (110) Sr2RuO2F2, the electron relaxation time was curtailed from 1.3 to 0.5 fs, while the hole relaxation time went from 7.10−2 to 2.10−2 fs as ϵz regressed from its unrestrained state to + 0.08. Conversely, compressive strain leads to lattice contraction, which enhances inter-atomic interactions and scattering processes, slows down the relaxation rate, and extends the relaxation time in both perpendicular and parallel directions. The fluctuations seen in figures 8 and 9 are correlated with alterations in the rate of phonon scattering within the crystal lattice, which is related to the type of inter-atomic interactions and the scattering process. The SOC manifests itself in a manner that is congruent with that of the hybrid functional approach (YS-PBE0) due to the fact that both methods incorporate the interactions between electrons and account for the relativistic impacts on the electron configuration. This assertion is substantiated by the graphical representation of the data (see figure 9). When scrutinizing the relationship between temperature and relaxation time in Sr2RuO2F2 (as depicted in figures 8 and 9), it becomes evident that this material is progressively exhibiting a higher degree of thermal excitation and a diminished resistance to alteration. This phenomenon is a result of the interplay between thermal energy and fluctuations within the Sr2RuO2F2 quaternary alloy, which excites its electrons and renders them more susceptible to alterations in their surroundings. The discernible decrease in relaxation time in both the perpendicular and parallel directions serves to underscore the material’s intensifying thermal activity, thus leading to a more expeditious reversion to a state of equilibrium.
Figure 9. Relaxation time: [(a) perpendicular direction and (c) parallel direction] using YS-PBE0 + SOC.
A comprehensive study of electron scattering in semiconductors has illuminated the critical relationship between the mobility of carriers, including electrons and holes, and their effective mass and relaxation time, which are intricately tied to the interactions between the carriers and lattice vibrations. To advance our knowledge of transport properties, we propose an analytical framework that delves into the mobility of electrons and holes in both perpendicular (⊥) and parallel (∥) orientations. This framework is intended to isolate the intrinsic properties of the carriers [50], providing a more nuanced and sophisticated understanding of the underlying transport properties. Our research constitutes a valuable addition to the growing body of knowledge in the field of semiconductor physics and further advances our understanding of this complex and multifaceted phenomenon. The alteration in the mobility of charge carriers in Sr2RuO2F2 observed through DFT and DPT calculations can be attributed to a complex interplay of multiple mechanisms in condensed matter physics [51]. Electron and hole mobility was significantly affected by strain in (⊥) orientation; it waned steadily under stretch and abruptly jumped to 2.60 $\times $ 102 cm2 V−1 s−1, 39.77 cm2 V−1 s−1 in the pressure zone at ϵz = + 0.08, respectively, as illustrated in figure 10. Analogously, uniaxial stretching efficiently stirred the electron transport properties in the (∥) direction.
Figure 10. Electron and hole mobility: [(a) perpendicular direction and (b) parallel direction] using YS-PBE0.
The outcomes suggest that the tensile strain progressively curtailed the (110) relaxation times of electrons and holes and reduced their mobilities to 5.10 cm2 V−1 s−1, 5 $\times $ 10−2 cm2 V−1 s−1 at ϵz = + 0.08, respectively. Conversely, a drastic linear upward trend in carrier mobility was noticed when the strain was lengthened from ϵz = −0.02 to ϵz = −0.08, as displayed in figure 10. Indeed, it is of the utmost significance to note that the data acquired through the YS-PBE0 approach in figure 10 exhibit a quantifiable correspondence with the results depicted in figure 11 where SOC was taken into consideration. It is evident that comparable behavior was manifested in both instances. A key factor contributing to the decline in mobility under tensile strain is the modification of the band structure of Sr2RuO2F2, as previously elucidated (refer to figure 3). Tensile strain caused alterations in the density of states and effective masses of the charge carriers, leading to changes in the scattering and transport properties and a decline in mobility. In addition, tensile strain introduced lattice defects and distortions into the crystal structure of Sr2RuO2F2, which serve as scattering centers for the charge carriers and reduce their mobility. The electronic correlation effects in Sr2RuO2F2 can also be modified by tensile strain, leading to further alterations in the scattering and transport properties of the charge carriers. On the other hand, pressure intensifies the electron–phonon interaction. This is due to the fact that the increase in pressure can cause modifications to the electron–phonon coupling, thereby influencing the charge carrier mobility. A robust coupling can result in heightened mobility by decreasing the scattering rate of the carriers. This phenomenon highlights the delicate interplay between compression, electron–phonon coupling and charge carrier mobility in solid-state systems. Furthermore, the susceptibility of the band gap in Sr2RuO2F2 is contingent upon the intricacies of its crystal lattice, which undergoes alteration upon application of compression. As previously elucidated, the pressure serves to reduce the band gap, enabling electrons and holes to traverse the material with greater ease, thus augmenting their intrinsic mobility.
Figure 11. Electron and hole mobility: [(c) perpendicular direction and (d) parallel direction] using YS-PBE0 + SOC.

6. Intrinsic resistivity of electrons and holes under strain

Intrinsic resistivity refers to a property of a material that describes its inherent resistance to the flow of electric current. It is an intrinsic characteristic of the material itself and does not depend on the material’s geometry or shape. Intrinsic resistivity is usually known to be dependent on various factors including the type and mobility of charge carriers present in the material. In this context, it is intriguing to note that the resistivity of Sr2RuO2F2 exhibits the opposite behavior as the mobilities of electrons and holes, as observed in figure 12. In conventional metals (e.g. to Sr2RuO4) [52], electrical conductivity is primarily governed by the high mobility of free electrons, while in semiconductors (e.g. to Sr2RuO2F2), electrical conductivity is facilitated by both electrons and holes [53], with the latter characterized by lower mobility. As a result, the resistivity of semiconductors is influenced by the mobility of both charge carriers. With regard to Sr2RuO2F2, there is an assumption that the perturbations inflicted upon the Fermi surface (i.e. the surface in momentum space that separates the occupied and unoccupied states) as a result of strain could potentially be undermining superconductivity. Superconductivity, a fascinating phenomenon characterized by the elimination of electrical resistance and the attainment of perfect conductivity below a critical temperature, could be hindered if the strain impinges upon the Fermi surface in a manner that hinders the essential prerequisites for superconductivity to transpire. A consequence of this would be the augmentation of electrical resistivity in (⊥) and (∥) orientations (see the stretch region), as illustrated in figure 12.
Figure 12. Intrinsic resistivity: [(a), (c) perpendicular direction and (b), (d) parallel direction] using YS-PBE0 and SOC.

7. Conclusion

The research presented here focuses on various properties of the Ruddlesden Popper-type superconductor Sr2RuO4 and its quaternary alloy Sr2RuO2F2. The electronic properties of Sr2RuO4 were studied using hybrid functional and SOC, which showed that the compound has a metallic ground state. Point defects, specifically oxygen substitution with fluorine, were investigated for their effects on the electronic structure, which resulted in a widening of the energy gap in both majority and minority spins. The ADPs of uniaxially strained Sr2RuO2F2 were also analyzed, and the results showed that the band edge positions of the material are highly responsive to strain. The VRH approximation method was used to determine the mechanical properties of the material, and the nine independent elastic constants were found to be close to each other. Finally, the study investigated the acoustic phonons and phonon–electron interaction in Sr2RuO2F2, which can cause phonon-mediated superconductivity. The electronic relaxation time of the material was found to be significantly longer in the perpendicular direction than in the parallel direction at room temperature, and the resistivity of the material was found to be influenced by the mobility of both charge carriers. The perturbations that may be inflicted upon the Fermi surface due to strain could potentially undermine superconductivity and lead to an increase in electrical resistivity in both orientations. Understanding these properties is crucial for the design and optimization of electronic devices, and further research in this area is needed.

The authors are grateful to Prof. P Blaha and Prof. K Schwarz at Vienna Technical University for providing the WIEN2k package, which was instrumental in conducting the calculations for this study. The authors also thank the WIEN2k group for valuable discussions and insight that contributed to the development of this research.

1
Moon S J 2006 Electronic structures of layered perovskite Sr2MO4 (M = Ru, Rh, And Ir) Phys. Rev. B 74 113104

DOI

2
Mazin I I Papaconstantopoulos D A Singh D J 2000 Tight-binding Hamiltonians for Sr-filled ruthenates: application to the gap anisotropy and Hall coefficient in Sr2ruo4 Phys. Rev. B 61 5223

DOI

3
Lutchyn R M Nagornykh P Yakovenko V M 2009 Frequency and temperature dependence of the anomalous ac Hall conductivity in a chiral px + ipy superconductor with impurities Phys. Rev. B 80 104508

DOI

4
Kashiwaya S 2019 Time-reversal invariant superconductivity of Sr2RuO4 revealed by Josephson effects Phys. Rev. B 100 094530

DOI

5
Tanatar M A 2001 Thermal conductivity of superconducting Sr2RuO4 in oriented magnetic fields Phys. Rev. B 63 064505

DOI

6
Hassinger E 2017 Vertical line nodes in the superconducting gap structure of Sr2RuO4 Phys. Rev. X 7 011032

DOI

7
Cao G 1997 Magnetic and transport properties of single-crystal Ca2RuO4: relationship to superconducting Sr2RuO4 Phys. Rev. B 56 R2916

DOI

8
Yazdani-Hamid M E Mehrabi V Exchange-correlation effects on the structural, electronic, and optical properties of Sr2RuO4 using DFTcalculations SSRN

DOI

9
Mackenzie A P 1996 Calculation of thermodynamic and transport properties of Sr2RuO4 at low temperatures using known Fermi surface parameters Physica C 263 P. 510 515

DOI

10
Sunko V 2019 Direct observation of a uniaxial stress-driven Lifshitz transition in Sr2RuO4 npj Quant. Mater. 4 46

DOI

11
Conyers E Vogt E 1956 Transport and deformation-potential theory for many-valley semiconductors with anisotropic scattering Phys. Rev. 101 944

DOI

12
Deng X Haule K Kotliar G 2016 Transport properties of metallic ruthenates: a DFT+ DMFT investigation Phys. Rev. Lett. 116 256401

DOI

13
Maeno Y 2011 Evaluation of spin-triplet superconductivity in Sr2RuO4 J. Phys. Soc. Jpn. 81 011009

DOI

14
Damascelli A 2001 Fermi surface of Sr2RuO4 from angle resolved photoemission J. Electron. Spectrosc. Relat. Phenom. 114 641 646

DOI

15
Huang W 2015 Nontopological nature of the edge current in a chiral p-wave superconductor Phys. Rev. B 91 094507

DOI

16
Pustogow A 2019 Constraints on the superconducting order parameter in Sr2RuO4 from oxygen-17 nuclear magnetic resonance Nature 574 72 75

DOI

17
Kallin C 2012 Chiral P-wave order in Sr2RuO4 Rep. Prog. Phys. 75 042501

DOI

18
Liu Y-C 2017 Theory of the evolution of superconductivity in Sr2RuO4 under anisotropic strain npj Quant. Mater. 2 12

DOI

19
Steppke A 2017 Strong peak in Tcof Sr2RuO4 under uniaxial pressure. Science 355 Eaaf9398

DOI

20
Mackenzie A P 2020 A personal perspective on the unconventional superconductivity of Sr2RuO4 J. Supercond. Novel Magn. 33 177 182

DOI

21
Ziati M Bekkioui N Ez-Zahraouy H 2022 Correlation between carrier mobility and effective mass in Sr2RuO4–XFx (X = 2) under uniaxial strain using the Yukawa screened PBE0 hybrid functional J. Phys. Chem. Solids 161 110409

DOI

22
Jin R 2000 Experimental observation of the selection rule in Josephson coupling between In and Sr2RuO4 Europhys. Lett. 51 341

DOI

23
Strand H 2019 Magnetic susceptibility in the Hund’s metal Sr2RuO4 from the dynamical two-particle vertex Aps March Meeting Abstracts B46 005

24
Gross E K U Dreizler R M 1986 LDA Density Approximations In Quantum Chemistry And Solid State Physics New York Plenum 353 379

25
Kane C L Mele E J 2005 Quantum spin Hall effect in graphene Phys. Rev. Lett. 95 226801

DOI

26
Henderson T M Paier J Scuseria, Gustavo E 2011 Accurate treatment of solids with the HSE screened hybrid Phys. Status Solidi (B) 248 767 774

DOI

27
Wang Q H 2013 Theory of superconductivity in a three-orbital model of Sr2RuO4 Europhys. Lett. 104 17013

DOI

28
Ziati M Hamid. E-Z 2021 Effect of uniaxial strain on structural, electronic and optical properties of Sr2RuO4−XFx: a DFT study Chem. Phys. 544 111132

DOI

29
Tyuterev V G Vast N 2006 Murnaghan’s equation of state for the electronic ground state energy Comput. Mater. Sci. 38 350 353

DOI

30
Giustino F 2017 Electron-phonon interactions from first principles Rev. Mod. Phys. 89 015003

DOI

31
Bardeen J Shockley W J P R 1950 Deformation potentials and mobilities in non-polar crystals Phys. Rev. 80 72

DOI

32
Sekiba D 2007 Strain-induced change in electronic structure of Cu (100) Phys. Rev. B 75 115404

DOI

33
Wiktor J Pasquarello A 2016 Absolute deformation potentials of two-dimensional materials Phys. Rev. B 94 245411

DOI

34
Faghaninia A Ager Iii Joel W Lo, Cynthia S 2015 Ab initio electronic transport model with explicit solution to the linearized Boltzmann transport equation Phys. Rev. B 91 235123

DOI

35
Ves S 1990 Cubic ZnS under pressure: optical-absorption edge, phase transition, and calculated equation of state Phys. Rev. B 42 9113

DOI

36
Brod M K Anand S Snyder G J 2022 The importance of avoided crossings in understanding high valley degeneracy in half-Heusler thermoelectric semiconductors Adv. Electron. Mater. 8 2101367

DOI

37
Harkins W D 1917 The evolution of the elements and the stability of complex atoms. I. A new periodic system which shows a relation between the abundance of the elements and the structure of the nuclei of atoms JACS 39 856 879

DOI

38
Kanoun M B Goumri-Said S Reshak A H 2009 Theoretical study of mechanical, electronic, chemical bonding and optical properties of Ti2SnC, Zr2SnC, Hf2SnC and Nb2SnC. Comput. Mater. Sci. 47 491 500

DOI

39
Ziati M Ez-Zahraouy H 2021 Mechanical properties and thermodynamic parameters of Sr2RuO4 and Sr2RuO2F2compounds under pressure and temperature effects: Voigt–Reuss–Hill approximations and Debye model J. Phys. Opt. Sci. 3 2 8

DOI

40
Du J 2021 Cerium-based lead-free chalcogenide perovskites for photovoltaics Phys. Rev. B 104 235206

DOI

41
Kamran S Chen K Chen L 2009 Ab initio examination of ductility features of fcc metals Phys. Rev. B 79 024106

DOI

42
Lee J Mahendra S Alvarez Pedro J J 2010 Nanomaterials in the construction industry: a review of their applications and environmental health and safety considerations ACS Nano 4 3580 3590

DOI

43
Zheng X F 2014 A review of thermoelectrics research–recent developments and potentials for sustainable and renewable energy applications Renew. Sustain. Energy Rev. 32 486 503

DOI

44
Grundmann M 2010 Physics Of Semiconductors Berlin Springer 775 776

45
Waldecker L 2017 Momentum-resolved view of electron-phonon coupling in multilayer WSe2 Phys. Rev. Lett. 119 036803

DOI

46
Sentef Michael A 2016 Theory of light-enhanced phonon-mediated superconductivity Phys. Rev. B 93 144506

DOI

47
Wang F Lee D-H 2011 The electron-pairing mechanism of iron-based superconductors Science 332 200 204

DOI

48
Song Q 2017 Ab initio study of electron mean free paths and thermoelectric properties of lead telluride Mater. Today Phys. 2 69 77

DOI

49
Liu J 2019 Carrier mobility and relaxation time in BiCuSeO Phys. Lett. A 383 125990

DOI

50
Li H 2023 Recent progress on synthesis, intrinsic properties and optoelectronic applications of perovskite single crystals Adv. Funct. Mater. 33 24 2214339

DOI

51
Guo L Liang P 2020 Condensed matter physics in time crystals New J. Phys. 22 075003

DOI

52
Chester G V Thellung A 1959 On the electrical conductivity of metals Proc. Phys. Soc. 73 745

DOI

53
Sirringhaus H Sakanoue T Chang J F 2012 Charge-transport physics of high-mobility molecular semiconductors Phys. Status Solidi B 249 9 1655 1676

DOI

Outlines

/