Electronic, magnetic, optical and thermoelectric properties of Ca₂Cr_1−xNi_xOsO₆ double perovskites

Abstract

With the help of density functional theory calculations, we explored the recently synthesized double perovskite material Ca₂CrOsO₆ and found it to be a ferrimagnetic insulator with a band gap of ∼0.6 eV. Its effective magnetic moment is found to be ∼0.23 μ_B per unit cell. The proposed behavior arises from the cooperative effect of spin–orbit coupling and Coulomb correlation of Cr-3d and Os-5d electrons along with the crystal field. Within the ferrimagnetic configuration, doping with 50% Ni in the Cr-sites resulted in a half-metallic state with a total moment of nearly zero, a characteristic of spintronic materials. Meanwhile, the optical study reveals that both ε₁^xx and ε₁^zz decrease first and then increase rapidly with increasing photon energy up to 1.055 eV. We also found optical anisotropy up to ∼14 eV, where the material becomes almost optically isotropic. This material has a plateau like region in the σ_xx and σ_zz parts of the optical conductivity due to a strong 3d–5d interband transition between Cr and Os. In addition, we performed thermoelectric calculations whose results predict that the material might not be good as a thermoelectric device due to its small power factor.

---

Introduction

In recent years, research on double perovskites (DPs) with chemical formula A₂BB′O₆ (where the B and B′ cations are transition metals and A is an alkaline or rare earth metal) has gained significant interest due to their novel properties arising from the combination of crystal field, spin orbit coupling (SOC) and Coulomb correlation (U). Transition metal (TM) DPs are found to show remarkable properties such as structural stability, high charge mobility, finite band gap, superconductivity, half-metallicity (HM), piezoelectricity, thermoelectricity, etc., which can be exploited in modern technological devices.^1–5 Reports on the realization of room-temperature colossal magnetoresistance and a HM state in Sr₂FeMoO₆ and Sr₂FeReO₆,^6,7 multiferroicity in Bi₂NiMnO₆,⁸ and magnetodielectricity in La₂NiMnO₆ (ref. 9) have motivated researchers to rigorously study these types of materials. In this regard, DPs may help in overcoming today’s major challenges through the development of efficient energy converters and storage devices. TM doped DPs exhibit HM behaviour,^10–13 in which the material exhibits a conducting state in one of the spin channels while the other spin channel is insulating or semiconducting. HM behaviour has been exhibited by diverse groups of materials such as Heusler alloys, shandites, the Ruddlesden–Popper series, etc.,^14–21 and is a key property for potential applications in spintronics.^22,23 The majority of experimentally determined HMs are ferromagnets (FMs) such as La_1−xSr_xMnO₃, and some are ferrimagnets (FiMs) with finite values of magnetic moments.^6,7,17,18

Recently, osmium perovskites have been observed to show novel properties. For instance, NaOsO₃ was found to be magnetically driven by a metal–insulator (MI) transition,²⁴ while a ferroelectric type structural transition was observed in LiOsO₃.²⁵ The CsOs₂O₇, RbOs₂O₇ and KOs₂O₇ pyrochlores exhibit unusual superconductivity.²⁶ In DPs such as Sr₂CuOsO₆ and Sr₂NiOsO₆, a MI state was identified,²⁷ while HM antiferromagnetism (HMAFM) was predicted in Sr₂CrOsO₆ and Sr₂CrRuO₆,¹ and for Ba₂NiOsO₆, FM with a Dirac–Mott insulating state was observed near 100 K for the bulk material and anomalous quantum Hall behavior was observed on the surface.^28,29 In Sr₂ScOsO₆, a transition to AFM was observed at 92 K (ref. 30) and a transition to FiM was observed at 725 K.³¹ Similarly, Ca₂MOsO₆ (where M = Mn, Fe, Co, and Ni) displays a FiM state with Curie temperature (T_C) of 305, 320, 145, and 175 K, respectively.^32–35 From density functional theory (DFT) calculations, PrSrMgIrO₆ was predicted to be an HMAFM with nearly zero effective magnetic moment per unit cell due to the combined effect of Coulomb correlation, U, and SOC.¹¹ A few DPs such as Ca₂MgOsO₆, and Sr₂MgOsO₆ have finite gaps at the Fermi level (E_F) due to the correlation effect and are reported to be AFM Mott insulators,³⁶ while Sr₂NiOsO₆, Sr₂FeOsO₆, Sr₂NiRuO₆, and SrLaBB′O₆ (B = Ni, Fe; B′ = Os, Ru)^37,38 show Mott insulating states under the combined effect of U and SOC.

Electron doping has also been considered in FiM DPs such as La_xCa_2−xCrWO₆, which demonstrates that doping generally increases the T_C.³⁹ DPs are known also as potential candidates for optoelectronic and photovoltaic device applications. The optically active region for these perovskites lies within 3 to 15 eV which is in the range of UV and IR spectra.^40–42 Frequencies suitable for optoelectronics can be absorbed in this region. Bi₂FeMnO₆ is one such example with a band gap of ∼0.8 eV that has been predicted as a suitable material for optical devices.⁴³

Recently, Morrow et al.⁴⁴ synthesized a new DP material Ca₂CrOsO₆ whose T_C was measured to be much higher than room temperature, i.e., 490 K. This study further reports that the atomic sites occupied by the two transition metals, i.e., Cr and Os cations, are in the ratios of 76% and 24%, respectively, which is a signature of antisite disorder in the system. Here, we are particularly motivated to explore the electronic and related properties of Ca₂CrOsO₆ due to (i) its high T_C above room temperature, (ii) the FiM insulating state, and (iii) the small effective magnetic moment. Upon electron doping, the material is expected to close the band gap in one spin-channel, giving rise to HM state; compensate the total moments nearly to zero; and the increase in T_C.³⁹ These features are expected to be important in new DP spintronic devices that work at room temperature.

In this work, we have carried out DFT calculations and found that Ca₂CrOsO₆ is a FiM insulator with a total magnetic moment, μ_tot = 0.23 μ_B per unit cell. When Ni is partially substituted in the Cr-sites, it contributes an additional five electrons to the system. This significantly changes the electronic as well as magnetic behaviour, giving rise to a nearly compensated HMAFM state. The optical study on Ca₂CrOsO₆ shows a strong optical interband transition between the Cr-3d and Os-5d bands in a plateau like region. We further consider the transport properties and calculate the Seebeck coefficient, thermal conductivity and power factor of the parent material Ca₂CrOsO₆.

Crystal structures and computational details

Ca₂CrOsO₆ crystallizes in the monoclinic structure with space group P2₁/n (14), and contains CrO₆ and OsO₆ octahedra [see Fig. 1]. The experimental lattice parameters chosen for the calculations are a = 5.3513 Å, b = 5.4561 Å, and c = 7.6204 Å, with β = 90.092°. The inter-octahedral Cr–O–Os bond angles are 153.3°, 152.6° and 153.8°, with average Cr–O and Os–O bond lengths of 1.972 Å and 1.954 Å respectively.

Fig. 1 Crystal structure of double perovskite Ca₂CrOsO₆.

We performed DFT calculations to explore the electronic, magnetic and related properties of 3D and 2D materials as reported recently.^45–48 The electronic and magnetic properties were calculated using the full-potential linearized augmented plane wave (FP-LAPW) method as implemented in the WIEN2k code.⁴⁹ The standard generalized gradient approximation (PBE-GGA)⁵⁰ was used, incorporating U of 4 eV for Cr (Ni) and 1.5 eV for Os.^35–37 The SOC was included via the second variational step.⁵¹ The magnetic anisotropy energy (MAE) was calculated by self consistent calculations in full relativistic mode considering five different orientations of the magnetization. The atomic sphere radii R_MT were 2.14, 1.94, 1.99, 2.01 and 1.64 Bohr for Ca, Cr, Ni, Os and O respectively, and a set of 500 k-points within the full Brillouin zone was used which corresponds to an 8 × 8 × 6 k-mesh. The convergence was set to 1 mRy per a.u. which gives reliable results for materials that contain transition metals.⁵²

In order to consider different magnetic configurations, the symmetry of the crystal has been lowered to P1 which corresponds to a total of 20 inequivalent atoms in which three types of oxygen atoms are present and form octahedra with Cr and Os. The tilt and rotation of the CrO₆ and OsO₆ octahedra result in Cr–O–Os bond angles within 152.6° to 153.8°. This variation plays a significant role in the magnetic properties.⁴⁴

For the optical property calculations, the full potential local orbital (FPLO) code⁵³ (version 18.00) has been used, from which the velocity matrix elements for the optical parameters were extracted. In order to minimize the major fluctuations in the lower energy range, a 16 × 16 × 16 k-mesh was used. The electrical transport properties have been calculated using the Boltzmann theory⁵⁴ and the relaxation time approximation as implemented in the Boltztrap code,⁵⁵ interfaced with WIEN2k.⁴⁹ The electrical conductivity and power factor are calculated with respect to time relaxation, τ; the Seebeck coefficient is independent of τ.

Results and discussion

Electronic and magnetic properties

We begin our discussion with the magnetic ground state. For this, we have considered four types of magnetic configurations, one FM state (FM1-↑↑↑↑), two AFM states (AF1-↑↓↑↓, AF2-↑↓↓↑), and one FiM state (FiM-↑↑↓↓). This was done after noting that the nonmagnetic (NM) state of the crystal was found to be the most unstable state, with the highest energy. The FiM configuration is found to be the magnetic ground state with the lowest energy among the other magnetic configurations. The relativistic calculations considered for spin quantisation align along the [100], [010], [001], [110] and [111] directions. For the easy axis along [001], the calculated MAE is 12.75 meV per unit cell.

In the parent material Ca₂CrOsO₆, Cr has a charge of +3 with a 3d³ configuration, involving the occupancy of the three t_2g states in the spin up channel and leaving the e_g states empty. As a result, in the density of states (DOS), the three occupied t_2g states lie in the valence region while the two empty e_g states move to the conduction region. In the spin down channel, all the Cr-d bands are unoccupied and thus lie far into the conduction region. In the case of Os, its charge is +5 with a 5d³ configuration. Therefore, only three electrons from Os will occupy the t_2g states in the spin down channel, leaving all the other states to shift to the conduction region.

The electronic DOS and band structures of Ca₂CrOsO₆ and Ca₂Cr_0.5Ni_0.5OsO₆ within GGA+U+SOC are shown in Fig. 2 and 3, respectively. Starting with the parent material Ca₂CrOsO₆, a band gap of ∼0.59 eV has been calculated (Fig. 2 and 3 (left)), signifying that Ca₂CrOsO₆ is an insulator. The results thus obtained are in good agreement with experimental and theoretical reports.^44,56 The major contribution to the total DOS around E_F (Fig. 2(a) (left)) is mainly from the Cr-3d and Os-5d orbitals hybridizing strongly with the O-2p orbitals (Fig. 2(d) (left)).

Fig. 2 The total and partial density of states (DOS) of Ca₂CrOsO₆ (left) and Ca₂Cr_0.5Ni_0.5OsO₆ (right). (a) Total DOS, (b) partial DOS of Cr-3d and Ni-3d, (c) PDOS of Os-5d, (d) PDOS of O-2p states for spin up (↑) and spin down (↓) channels within the GGA+U+SOC functional. The vertical dotted line indicates E_F = 0.

Fig. 3 The band structures of Ca₂CrOsO₆ (left) and Ca₂Cr_0.5Ni_0.5OsO₆ (right) within GGA+U+SOC. The horizontal dotted line indicates E_F = 0.

On the other hand, for the doped material Ca₂Cr_0.5Ni_0.5OsO₆, one Ni atom is substituted in a Cr-site in Ca₂CrOsO₆. This will add five extra electrons to the system as Ni with a charge of +2 has a 3d⁸ configuration. Out of five additional electrons, two occupy the e_g states of the spin up channel, while the remaining three are bound to go to the spin down channel to occupy the t_2g states. These three electrons of Ni in the spin down channel generate a repulsive effect in the Os-t_2g band, which causes the Os band to shift from the valence region towards the conduction region by crossing E_F in the spin down channel (see Fig. 2(c) (right)). This gives rise to a metallic state in the spin down channel, while the material remains insulating in the spin up channel. Thus, with an insulating state in the spin up channel and a metallic state in the spin down channel, the Ca₂Cr_0.5Ni_0.5OsO₆ material is predicted to be a HM. It is interesting to note that the HM state remains robust for U_Os as large as 4 eV. The DOS and band structures for the parent and the doped material within GGA and GGA+U are available in the ESI section for comparison (see Fig. S1–S4†).

We also considered the magnetic behaviour of Ca₂Cr_1−xNi_xOsO₆. As tabulated in Table 1, for x = 0, the individual magnetic moments of the Cr and Os atoms are 2.52 μ_B and −1.59 μ_B respectively, with an effective magnetic moment of 0.23 μ_B per unit cell (2 f.u.). Similarly, for x = 0.5, the individual magnetic moments of the Cr, Ni and Os atoms are 2.52 μ_B, 1.64 μ_B and −1.42 μ_B respectively, with an effective magnetic moment of 0.21 μ_B per unit cell. Due to the partial charge transfer from Cr and Os to oxygen, the oxygen atoms are spin polarized in parallel with the Os ions and gain a sizable magnetic moment of −0.07 μ_B, consistent with the isosurface plot shown in Fig. 4 (left). Here, the polarisation is mainly found in the 2p orbitals and this hybridized moment in oxygen increases the resultant magnetic moment.

Table 1 Calculated spin magnetic moments (in μ_B) of B (Cr/Ni), B′ (Os) and the 3 types of oxygen atoms, net moments (per unit cell) and band gaps, E_g, (in eV). The calculated orbital moments at the B and B′ sites are shown within parentheses for Ca₂CrOsO₆ and Ca₂Cr_1−xNi_xOsO₆

Site	Ca2CrOsO6				Ca2Cr0.5Ni0.5OsO6
	GGA	GGA+U	GGA+U+SOC	Expt.	GGA	GGA+U	GGA+U+SOC
B	2.20	2.52	2.52(0.05)	2.5	2.19/1.32	2.52/1.65	2.52(0.05)/1.64(0.14)
B′	−1.45	−1.66	−1.59(0.14)	1.26	−1.28	−1.48	−1.42(0.16)
O1	−0.10	−0.11	−0.10		−0.08	−0.09	−0.08
O2	−0.09	−0.10	−0.10		−0.03	−0.05	−0.05
O3	−0.09	−0.10	−0.10		−0.05	−0.07	−0.07
Net	0	0	0.23	0.2	0.0	0	0.21
Eg	0.53	0.65	0.59		Metallic	Metallic	Metallic

---

Fig. 4 Isosurfaces of spin magnetization density at ±0.21 e Å⁻³ with red (blue) for spin up (down): (left) parent material Ca₂CrOsO₆, (right) material with 50% Cr replacement Ca₂Cr_0.5Ni_0.5OsO₆.

To understand the charge transfer effects among Cr, Ni, Os and O atoms, spin density isosurfaces are drawn as shown in Fig. 4 (right). As observed for the parent compound Ca₂CrOsO₆ (see Fig. 4 (left)), the spin density contributions for Cr and Os are due to 3d and 5d states. The observed characteristics are due to the combination of d_xy, d_xz and d_yz (i.e., t_2g) states that are indicative of d³ configurations. In addition, the oxygen-2p states become spin-polarized due to the Cr-3d and Os-5d states. When one of the Cr atoms is replaced by Ni, the nature of the spin density changes significantly (see Fig. 4 (right)). Actually, Ni with a charge of +2 has a d⁸ configuration. Thus, with five electrons fully occupying one spin channel, the remaining three electrons move to another spin channel to fill the t_2g states. As a result, the remaining visible states are the e_g states which are clearly seen in the isosurfaces of the Ni (see Fig. 4 (right)) atoms.

From Table 1, the spin magnetic moment in Os is reduced due to the strong Os–O hybridization, and SOC induces an orbital moment. Thus, the SOC is found to be responsible for the small reduction in spin magnetic moment. The alignment of the orbital magnetic moment of Cr is the same as that of the total spin magnetic moment, showing that the 3d orbital is either semi occupied or more than half filled. Whereas, in the Os-5d state, the spin orientations are antiparallel as compared to the orbital magnetic moments due to the less than half filled states, which is in accordance with Hund’s third rule.⁵⁷ Hence, the orbital magnetic moments on Cr and Os are parallel which increases the net magnetic moment in Ca₂CrOsO₆.

Optical properties

In this work, we have investigated the optical properties of the double perovskite Ca₂CrOsO₆, considering the polarization directions E_xx and E_zz, as E_xx = E_yy due to the tetragonal symmetry. The optical current j_a generated by an electric field E_b is given by: j_a = σ_abE_b. For Ca₂CrOsO₆, we consider only E_x and E_z polarization as E_x = E_y due to the tetragonal symmetry, leaving only two independent x and z components of the optical element. The complex dielectric function ε(ω) = ε₁(ω) + iε₂(ω) describes the behavior of materials under incident light, where ε₁(ω) and ε₂(ω) are the real and imaginary parts of the dielectric function, respectively. The imaginary part ε₂(ω) is obtained from the electronic structure through the joint density of states and the momentum matrix elements between the occupied and unoccupied states:where p is the momentum operator, kn is the eigenfunction with eigenvalue E_kn and f(kn) is the Fermi distribution function. The real part of the dielectric function ε₁ is given by the Kramers–Kronig dispersion relations.⁵⁸

The other optical parameters like optical conductivity σ(ω) and electron loss function L(ω) are directly related to the dielectric function as follows.⁵⁹

The real part of the complex dielectric function ε₁(ω) expresses the electronic polarizability of materials due to incident photons. Fig. 5 (left) shows the variation of the dielectric function (real and imaginary) as a function of photon energy. The calculated static dielectric constants for x and z polarization (ε₁^xx(0) and ε₁^zz(0)) are 15.94 eV and 12.43 eV, respectively (Table 2). Both ε₁^xx and ε₁^zz decrease first and then increase rapidly with increasing photon energy up to 1.055 eV. With a small increase in photon energy above 0.0 eV, ε₁^zz decreases more as compared to ε₁^xx and then increases again. At a certain energy, both ε₁^xx and ε₁^zz show a maximum peak where ε₁^xx surpasses ε₁^zz in magnitude. The plot also indicates the presence of optical anisotropy up to ∼14 eV, and after that the material becomes almost optically isotropic. Both ε₁^xx and ε₁^zz become zero at ∼7.6 eV which may be due to the occurrence of plasmonic type oscillations. We know that the imaginary part ε₂(ω) of the complex dielectric function is related to optical absorption. The dependence of ε₂(ω) on photon energy is presented in Fig. 5 (left).

Fig. 5 (Left) The real and imaginary parts of the dielectric function ε(ω) for the interband part. The intraband part of ε(ω) diverges at ω = 0; (middle) total, intraband and interband loss function L(ω); (right) total, interband and individual spin contributions to optical conductivity σ(ω). Here, for Ca₂CrOsO₆, σ_xx(ω) and σ_zz(ω) contribute.

Table 2 Calculated optical parameters at the main peak value for the Ca₂CrOsO₆ compound

Optical parameters	ε1^xx(0)	ε1^zz(0)	L(ω)	σ1^xx(ω)	σ1^zz(ω)
eV	15.94	12.43	14.11	7.79	7.32

---

From the figure, we notice a rapid increase in ε₂(ω) beyond the threshold energy which corresponds to the calculated energy band gap and observed a maximum peak at 1.0–1.5 eV. After that, the value of ε₂(ω) decreases with severe fluctuations. A significant plateau like region is obtained in the energy range from 3.0 eV up to 10.0 eV. This can also be verified by the optical conduction plot which shows active absorption of photon energy due to direct electron transitions. The variation of the electron energy loss function L(ω) versus photon energy is also plotted in Fig. 5 (middle). L(ω) gives a description of the energy loss due to the scattering of a fast electron travelling in the material. The electron energy loss is very high in the low energy region in the case of intraband transitions. However, for interband transitions, L(ω) slowly increases up to 12.0 eV and then abruptly increases beyond 13.0 to 14.0 eV. The low value of electron energy loss is due to the fact that L(ω) is inversely related to ε₂(ω). The large peak in the energy loss spectra at around 13.0 eV to 14.0 eV corresponds to a plasmon frequency ω_p at 13.0–14.0 eV. We also calculated the optical conductivity σ(ω) as shown in Fig. 5 (right). The only independent components of σ(ω) in this material are σ_xx(ω), σ_zz(ω), σ_yx(ω) and σ_zx(ω). σ_yx(ω) and σ_zx(ω) are zero, and only the σ_xx(ω) and σ_zz(ω) components contribute. The optical conductivity of Ca₂CrOsO₆ has a plateau like region at 3–12 eV rather than a sharp peak like feature. This indicates a strong optical interband transition between the Cr-3d and Os-5d bands in this region.

Thermoelectric properties

The insulating behavior of Ca₂CrOsO₆ prompted us to investigate the thermoelectric properties. The performance of thermoelectric materials is given by a dimensionless parameter called the figure of merit, i.e., ZT = S²σT/κ, where S, σ and κ represent the Seebeck coefficient, electrical conductivity, and thermal conductivity, respectively.^60–62 Conventional materials such as Bi₂Te₃ and PbTe have shown a promising figure of merit, however, these alloys suffer from toxic compositions, low thermal stability, and expensive constituents. In this regard, oxide materials could be a promising alternative owing to their environmental friendliness, non-toxicity, good thermal stability, and low-cost compositions.^63,64 Recently, good thermoelectric performances have been reported for double perovskite oxides.^65,66 Motivated by these findings, we investigated the thermoelectric properties of Ca₂CrOsO₆ utilizing the Boltztrap code which is used for computing the electrical transport coefficients S and σ, and thereby S²σ, within the rigid band approximation (RBA) and constant relaxation time approach (CRTA). The CRTA assumes the Seebeck coefficient to be independent of relaxation time, τ, and the electrical conductivity and power factor (PF) are expressed with respect to τ. These two approaches have been successful in the theoretical investigation of new thermoelectric materials by many groups.^67–70

Fig. 6 shows the Seebeck coefficient, electrical conductivity and power factor with respect to relaxation time, as a function of temperature for Ca₂CrOsO₆ for the spin up and spin down channels, and the total contribution. The parameters are calculated for both the spin channels according to the two current model.⁷¹ Within this model, the total electrical conductivity and Seebeck coefficients are given by σ = σ(↑) + σ(↓) and S = [S(↑)σ(↑) + S(↓)σ(↓)]/[σ(↑) + σ(↓)] where (↑) and (↓) represent the coefficients for spin up and spin down respectively. Utilizing this model, the total Seebeck coefficient obtained shows an increasing trend in the temperature range 300–900 K. The positive values are in the range 43–83 μV K⁻¹ and are indicative of p-type charge carriers. However, the Seebeck coefficient values are not significant as compared to benchmark values (150 μV K⁻¹).⁷² For evaluating the electrical conductivity and power factor, one needs to incorporate the relaxation time which is not trivial to calculate. The relaxation time is very sensitive as it may vary for different systems at different doping levels and temperatures.⁷³ Thus, the actual value of the relaxation time is quite essential for determining the transport properties. However, in previous theoretical studies, a relaxation time of 5 × 10⁻¹⁵ s has been used for evaluating the transport properties.⁷⁴ Utilizing the same value, we obtained good values of total electrical conductivity in the range 429–404 S cm⁻¹ in the temperature range 300–900 K. These values can be attributed to the greater density of states at the VBM owing to the degenerate bands. The degenerate bands increase the density of states which results in a larger number of charge carriers. The greater density of states at the VBM close to the Fermi level as compared to the conduction band minimum suggests that hole doping would be beneficial for improving the transport properties.⁷⁵

Fig. 6 (a) Seebeck coefficient, (b) electrical conductivity, and (c) power factor of Ca₂CrOsO₆ as a function of temperature for spin up (black), spin down (red) and total contribution (green). The electrical conductivity and power factor are plotted with respect to the relaxation time.

The calculated total power factor with respect to the relaxation time is not significant, owing to the low Seebeck coefficient. The values range from 0.08 μW cm⁻¹ K⁻² s⁻¹ at room temperature to 3 μW cm⁻¹ K⁻² s⁻¹. Utilizing the relaxation time as before, i.e., τ = 5 × 10⁻¹⁵ s, the maximum obtained power factor is 15 μW cm⁻¹ K⁻² at 900 K. The proposed value, though not significant, can be further improved by utilizing suitable dopants and synthesis approaches, meriting further experimental investigations.

Conclusions

By means of density functional calculations, we found that the double perovskite material Ca₂CrOsO₆ is a ferrimagnetic insulator. Upon substantial doping with Ni (50%) in the Cr-sites, the composition Ca₂Cr_0.5Ni_0.5OsO₆ shows nearly compensated half-metal behaviour which is desirable for spintronic device applications. The optical study of the parent material Ca₂CrOsO₆ reveals that both ε₁^xx and ε₁^zz decrease first and then increase rapidly with increasing photon energy up to 1.055 eV. At a certain energy, both ε₁^xx and ε₁^zz show a maximum peak where ε₁^xx surpasses ε₁^zz in magnitude. We also found that optical anisotropy is present up to ∼14 eV. The promising optical response of Ca₂CrOsO₆ due to Cr-3d and Os-5d interband transitions makes it an important material for photovoltaic and optoelectronic applications. In addition, we investigated the thermoelectric properties, however, the material does not show promise as a thermoelectric device.

Conflicts of interest

There are no conflicts to declare.

Acknowledgements

SRB thanks the Nepal Academy of Science and Technology, Kathmandu for the PhD fellowship. Computations were performed at IFW Dresden, Germany and CDP, Nepal. We thank Ulrike Nitzsche for technical assistance. SRB and MPG thank Manuel Richter and Jeroen van den Brink for fruitful discussions and suggestions to complete this work. MPG acknowledges the Higher Education Reform Project of Tribhuvan University for the start-up grant, and the Alexander von Humboldt Foundation, Germany for the equipment grants.

Notes and references

1. K. W. Lee and W. E. Pickett, Phys. Rev. B: Condens. Matter Mater. Phys., 2008, 77, 115101.
2. Y. P. Liu, H. R. Fuh and Y. K. Wang, Comput. Mater. Sci., 2014, 92, 63.
3. A. H. Reshak, Phys. Chem. Chem. Phys., 2016, 6, 92887.
4. D. P. Rai, A. Shankar, M. P. Ghimire and R. K. Thapa, Comput. Mater. Sci., 2015, 101, 313.
5. M. T. Anderson, K. B. Greenwood, G. A. Taylor and K. R. Poeppelmeier, Prog. Solid State Chem., 1993, 22, 197.
6. Y. Tokura, Rep. Prog. Phys., 2006, 69, 797.
7. K.-I. Kobayashi, T. Kimura, H. Sawada, K. Terakura and Y. Tokura, Nature, 1998, 395, 677.
8. Y. Shimakawa, M. Azuma and N. Ichikawa, Materials, 2011, 4, 153.
9. N. S. Rogado, J. Li, A. W. Sleight and M. A. Subrama- nian, Adv. Mater., 2005, 17, 2225.
10. W. E. Pickett, Phys. Rev. B: Condens. Matter Mater. Phys., 1998, 57, 10613.
11. M. P. Ghimire, L. H. Wu and X. Hu, Phys. Rev. B, 2016, 93, 134421.
12. X. Hu, Adv. Mater., 2012, 24, 294.
13. Y.-M. Nie and X. Hu, Phys. Rev. Lett., 2008, 100, 117203.
14. M. Sargolzaei, M. Richter, K. Koepernik, I. Opahle and H. Eschrig, Phys. Rev. B: Condens. Matter Mater. Phys., 2006, 74, 224410.
15. M. P. Ghimire, J. I. Facio, J.-S. You, L. Ye, J. G. Checkelsky, S. Fang, E. Kaxiras, M. Richter and J. van den Brink, Phys. Rev. Res., 2019, 1, 032044.
16. W. E. Pickett and J. S. Moodera, Phys. Today, 2001, 54, 39.
17. G. M. Muller, J. Walowski, M. Djordjevic, G.-X. Miao, A. Gupta, A. V. Ramos, K. Gehrke, V. Moshnyaga, K. Samwer, J. Schmalhorst, A. Thomas, A. Hütten, G. Reiss, J. S. Moodera and M. Münzenberg, Nat. Mater., 2009, 8, 56.
18. G. Shibata, M. Kitamura and M. Minohara, et al., npj Quantum Mater., 2018, 3, 3.
19. S. N. Ruddlesden and P. Popper, Acta Crystallogr., 1957, 10, 538.
20. M. P. Ghimire, R. K. Thapa, D. P. Rai, Sandeep, T. P. Sinha and X. Hu, J. Appl. Phys., 2015, 117, 063903.
21. H. van Leuken and R. A. de Groot, Phys. Rev. Lett., 1995, 74, 1171.
22. C. Felser, G. H. Fecher and B. Balke, Angew. Chem., Int. Ed., 2007, 46, 668.
23. M. I. Katsnelson, Y. Yu. Irkhin, L. Chioncel, A. I. Licht-enstein and R. A. de Groot, Rev. Mod. Phys., 2008, 80, 315.
24. S. Calder, V. O. Garlea, D. F. McMorrow, M. D. Lumsden, M. B. Stone, J. C. Lang, J. W. Kim, J. A. Schlueter, Y. G. Shi, K. Yamaura, Y. S. Sun, Y. Tsujimoto and A. D. Christianson, Phys. Rev. Lett., 2012, 108, 257209.
25. Y. G. Shi, Y. F. Guo, X. Wang, A. J. Princep, D. Khalyavin, P. Manuel, Y. Michiue, A. Sato, T. Tsuda, S. Yu, M. Arai, Y. Shirako, M. Akaogi, N. L. Wang, K. Yamaura and A. T. Boothroyd, Nat. Mater., 2013, 12, 1023.
26. Z. Hiroi, J. Yamaura and K. Hattori, J. Phys. Soc. Jpn., 2012, 81, 011012.
27. A. C. Tian, H. C. Wibowo, Z. Loye and M. H. Whangbo, Inorg. Chem., 2011, 50, 4142.
28. H. L. Feng, S. Calder, M. P. Ghimire, Y. Yahua, Y. Shirako, Y. Tsujimoto, Y. Matsushita, Z. Hu, C.-Y. Kuo, L. H. Tjeng, T.-W. Pi, Y.-L. Soo, J. He, M. Tanaka, Y. Katsuya, M. Richter and K. Yamaura, Phys. Rev. B, 2016, 94, 235158.
29. H.-S. Lu and G.-Y. Guo, Phys. Rev. B, 2019, 100, 054443.
30. A. E. Taylor, R. Morrow, D. J. Singh, S. Calder, M. D. Lumsden, P. M. Woodward and A. D. Christianson, Phys. Rev. B: Condens. Matter Mater. Phys., 2015, 91, 100406.
31. Y. Krockenberger, K. Mogare, M. Reehuis, M. Tovar, M. Jansen, G. Vaitheeswaran, V. Kanchana, F. Bultmark, A. Delin, F. Wilhelm, A. Rogalev, A. Winkler and L. Alff, Phys. Rev. B: Condens. Matter Mater. Phys., 2007, 75, 020404.
32. H. L. Feng, M. Arai, Y. Matsushita, Y. Tsujimoto, Y. Guo, C. I. Sathish, X. Wang, Y. H. Yuan, M. Tanaka and K. Yamaura, J. Am. Chem. Soc., 2014, 136, 3326.
33. M. P. Ghimire and X. Hu, Mater. Res. Express, 2016, 3, 106107.
34. R. Morrow, K. Samanta, T. Saha-Dasgupta, J. Xiong, J. W. Freeland, D. Haskel and P. M. Woodward, Chem. Mater., 2016, 28, 3666.
35. H. L. Feng, M. P. Ghimire, Z. Hu, S.-C. Liao, S. Agrestini, J. Chen, Y. Yuan, Y. Matsushita, Y. Tsujimoto, Y. Katsuya, M. Tanaka, H.-J. Lin, C.-T. Chen, S.-C. Weng, M. Valvidares, K. Chen, F. Baudelet, A. Tanaka, M. Greenblatt, L. H. Tjeng and K. Yamaura, Phys. Rev. Mater., 2019, 3, 124404.
36. Y. Yuan, H. L. Feng, M. P. Ghimire, Y. Matsushita, Y. Tsujimoto, J. He, M. Tanaka, Y. Katsuya and K. Yamaura, Inorg. Chem., 2015, 54, 3422.
37. D. K. Yadav, S. R. Bhandari, B. P. Belbase, G. C. Kaphle, D. P. Rai and M. P. Ghimire, Comput. Mater. Sci., 2019, 170, 109168.
38. W. Song, E. Zhao, J. Meng and Z. Wu, J. Chem. Phys., 2009, 130, 114707.
39. S. Geprags, P. Majewski and R. Gross, J. Appl. Phys., 2006, 99, 08J102.
40. M. Nabi and D. C. Gupta, RSC Adv., 2019, 9, 15852.
41. K. A. Parrey, S. A. Khandy, I. Islam, A. Laref, D. C. Gupta, A. Niazi, A. Aziz, S. G. Ansari, R. Khenata and S. Rubab, J. Electron. Mater., 2018, 47, 7.
42. M. Roknuzzaman, C. Zhang, K. Ostrikov, A. Du, H. Wang, L. Wang and T. Tesfamichael, Sci. Rep., 2019, 9, 718.
43. T. Ahmed, A. Chen, D. A. Yarotski, S. A. Trugman, Q. Jia and J.-X. Zhu, APL Mater., 2017, 5, 035601.
44. R. Morrow, J. R. Soliz, A. J. Hauser, J. C. Gallagher, M. A. Susner, M. D. Sumption, A. A. Aczel, J. Yan, F. Yang and P. M. Woodward, J. Solid State Chem., 2016, 238, 46.
45. S. R. Bhandari, R. K. Thapa and M. P. Ghimire, J. Nep. Phys. Soc., 2015, 3, 89.
46. M. Rafique, S. Yong and H. Tan, Phys. E., 2017, 88, 115.
47. M. Rafique, M. A. Unar, I. Ahmed, A. R. Chachar and Y. Shuai, J. Phys. Chem. Solids, 2018, 118, 114.
48. M. Rafique, S. Yong and H. Tan, J. Mater. Chem. C, 2017, 5, 8112.
49. P. Blaha, K. Schwarz, G. K. H. Madsen, D. Kvasnicka and J. Luitz, An Augmented Plane Wave Plus Local Orbitals Program for Calculating Crystal Properties, Technische Universität Wien, Vienna, Austria, 2001, ISBN 3-9501031-1-2.
50. J. P. Perdew, K. Burke and M. Ernzerhof, Phys. Rev. Lett., 1996, 77, 3865.
51. J. Kunes, P. Novak, M. Divis and P. M. Oppeneer, Phys. Rev. B: Condens. Matter Mater. Phys., 2001, 63, 205111.
52. R. Muhammad, M. A. Uqaili, Y. Shuai, M. A. Mahar and I. Ahmed, Appl. Surf. Sci., 2018, 458, 145.
53. K. Koepernik and H. Eschrig, Phys. Rev. B: Condens. Matter Mater. Phys., 1999, 59, 1743 , https://www.fplo.de/.
54. P. B. Allen, Boltzmann theory and resistivity of metals in Quantum Theory of Real Materials, ed. J. R. Chelikowsky and S. G. Louie, Kluwer, Boston, 1996.
55. G. K. H. Madsen and D. J. Singh, Comput. Phys. Commun., 2006, 67, 175.
56. N. Zu, R. Li and R. Ai, J. Magn. Magn. Mater., 2018, 467, 145–149.
57. H. T. Jeng and G. Y. Guo, Phys. Rev. B: Condens. Matter Mater. Phys., 2003, 67, 094438.
58. J. S. Toll, Phys. Rev., 1956, 104, 1760.
59. D. Hoat, J. R. Silva and A. M. Blas, J. Solid State Chem., 2019, 85, 270.
60. H. Sevincli, C. Sevik, T. Cagin and G. Cuniberti, Sci. Rep., 2013, 3, 1228.
61. S. R. Boona, Nature, 2017, 549, 169.
62. G. Tan, S. Hao, J. Zhao, C. Wolverton and M. G. Kanatzidis, J. Am. Chem. Soc., 2017, 139, 6467.
63. P. Roy, V. Waghmare and T. Maiti, RSC Adv., 2016, 6, 54636.
64. M. Saxena, K. Tanwar and T. Maiti, Scr. Mater., 2017, 130, 205.
65. P. Roy, I. Bose and T. Maiti, Integr. Ferroelectr., 2016, 34, 174.
66. T. Sugahara, M. Ohtaki and T. Souma, J. Ceram. Soc. Jpn., 2008, 116, 1278.
67. G. K. H. Madsen, J. Am. Chem. Soc., 2006, 128, 12140.
68. J. Yang, H. Li, T. Wu, W. Zhang, L. Chen and J. Yang, Adv. Funct. Mater., 2008, 18, 2880.
69. M.-S. Lee, F. P. Poudeu and S. D. Mahanti, Phys. Rev. B: Condens. Matter Mater. Phys., 2011, 83, 085204.
70. M. Zeeshan, H. K. Singh, J. van den Brink and H. C. Kandpal, Phys. Rev. Mater., 2017, 1, 075407.
71. S. Sharma and S. K. Panday, Phys. Lett. A, 2015, 379, 2357.
72. T. M. Tritt, Annu. Rev. Mater. Res., 2011, 41, 433.
73. M. Zeeshan, J. van den Brink and H. C. Kandpal, Phys. Rev. Mater., 2017, 1, 074401.
74. K. A. Parrey, S. A. Khandy, I. Islam, A. Laref, D. C. Gupta, A. Niazi, A. Aziz, S. G. Ansari, R. Khenata and S. Rubab, J. Electron. Mater., 2018, 47, 3615.
75. M. Zeeshan, T. Nautiyal, J. van den Brink and H. C. Kandpal, Phys. Rev. Mater., 2018, 2, 065407.

---

Footnote

† Electronic supplementary information (ESI) available. See DOI: 10.1039/c9ra10775d

---