Recent Progress in Materials is an international peer-reviewed Open Access journal published quarterly online by LIDSEN Publishing Inc. This periodical is devoted to publishing high-quality papers that describe the most significant and cutting-edge research in all areas of Materials. Its aim is to provide timely, authoritative introductions to current thinking, developments and research in carefully selected topics. Also, it aims to enhance the international exchange of scientific activities in materials science and technology.
Recent Progress in Materials publishes original high quality experimental and theoretical papers and reviews on basic and applied research in the field of materials science and engineering, with focus on synthesis, processing, constitution, and properties of all classes of materials. Particular emphasis is placed on microstructural design, phase relations, computational thermodynamics, and kinetics at the nano to macro scale. Contributions may also focus on progress in advanced characterization techniques.          

Main research areas include (but are not limited to):
Characterization & Evaluation of Materials
Metallic materials 
Inorganic nonmetallic materials 
Composite materials
Polymer Materials
Sustainable Materials and Technologies
Special types of Materials
Macro-, micro- and nano structure of materials
Environmental interactions, process modeling
Novel applications of materials

Archiving: full-text archived in CLOCKSS.

Rapid publication: manuscripts are undertaken in 13 days from acceptance to publication (median values for papers published in this journal in 2021, 1-2 days of FREE language polishing time is also included in this period).

Current Issue: 2022  Archive: 2021 2020 2019
Open Access Research Article

Optical, Charge Transport, Thermal, Magnetic, Plasmonic, and Quantum Mechanical Properties of Iridium

William E Vargas 1,*, Federico Muñoz-Rojas 1, Esteban Avendano 1, Victoria Quirós-Cordero 2, Marcela Hernández-Jiménez 1

  1. Centro de Investigación en Ciencia e Ingeniería de Materiales, Universidad de Costa Rica, San José, Costa Rica

  2. School of Materials Science and Engineering, College of Engineering, Georgia Institute of Technology, Atlanta, Georgia, United States

Correspondence: William E. Vargas

Academic Editor: Vardan Apinyan

Special Issue: Quantum Mechanics in Solid State Systems

Received: July 30, 2022 | Accepted: September 23, 2022 | Published: October 12, 2022

Recent Progress in Materials 2022, Volume 4, Issue 4, doi:10.21926/rpm.2204019

Recommended citation: Vargas WE, Muñoz-Rojas F, Avendano E, Quirós-Cordero V, Hernández-Jiménez M. Optical, Charge Transport, Thermal, Magnetic, Plasmonic, and Quantum Mechanical Properties of Iridium. Recent Progress in Materials 2022; 4(4): 019; doi:10.21926/rpm.2204019.

© 2022 by the authors. This is an open access article distributed under the conditions of the Creative Commons by Attribution License, which permits unrestricted use, distribution, and reproduction in any medium or format, provided the original work is correctly cited.


Spectrophotometry has been widely used to retrieve the dielectric function of a bulk iridium sample using an extended version of the Drude–Lorentz model. The parameters of the model are optimized using a spectral-projected-gradient-method-assisted acceptance-probability-controlled simulated annealing approach. Furthermore, optimized values of Drude parameters corresponding to the optical response of electrons and holes (scattering frequency of electrons, the ratio between scattering frequencies of holes and electrons, the ratio between effective masses of electrons and holes, the ratio between the number densities of holes and electrons, and electron volume plasma frequency) are used to evaluate charge transport and magnetic properties. These include static and dynamic conductivities, intrinsic mean free paths, the effective mass of charge carriers and their number densities, Fermi velocities and energies, densities of states at Fermi energies, mobilities, specific heats, Hall’s coefficient, thermal conductivities, charge carrier coupling constant, paramagnetic and diamagnetic susceptibilities, and the number of Bohr magnetons. In addition, optimized resonance energy values of the Lorentz contribution to the dielectric function were compared with the background information provided by density-functional-theory calculations for iridium. A decomposition of the energy loss function was used as the starting point to calculate the effective numbers of bound electrons involved in interband transitions, as well as the densities of states at the final energies of the sets of transitions considered. The Drude–Lorentz model involves charge carrier parameters for both electrons and holes, as well as the resonance energies correlating with the energies associated with quantum transitions. To a large extent, several physical quantities calculated from optimized parameters exhibit values close to those obtained from measurements or by applying other models, including quantum mechanics formulations.

Graphical abstract

Click to view original image


Charge transport; reflectivity; density of states; magnetic susceptibility; Drude–Lorentz model; electron localization; iridium; density functional theory; simulated annealing

1. Introduction

The transition metal iridium (Ir) has been increasingly used in the industry due to its physical properties, including high mechanical strength, low oxidation at ambient conditions, corrosion resistant, large work function, high melting point, and high electrical conductivity. It has been widely applied in thermocouples and as an encapsulating material to stabilize and provide mechanical protection to radioisotopes used as nuclear fuels [1], in microelectronics for designing and implementing dynamic- and ferroelectric- random access memory electrodes [2], as gate electrodes in metal-oxide-semiconductor field-effect transistors [3], and as ohmic contact in diamond-based electronic devices [4], as coatings in spatial instrumentation to protect against ultraviolet radiation [5], in gas sensors [6] and catalytic processes [7,8]. In addition, Ir has been implicated in the design and characterization of metallic photonic-band-gap crystals [9]. Iridium oxide has been used in optically variable electrochromic materials used in smart windows [10], as an electrode catalytic material in water electrolysis [11], and in several other industrial and technological applications [12].

The optical, charge-carrier transport, and magnetic properties, as well as band structure calculations of Ir have not been well studied. Reports on the optical properties of Ir are scarce, with the first study reported in 1910, which focused on metal mirrors containing the measured spectral reflectivity of Ir for near- and mid-infrared wavelengths [13]. In 1966, a study on the optical constants and emissivity of heated semicircular filaments of refractory metals, including Ir, for near-infrared and visible photon energies was published [14]. A year later, the optical constants of Ir films were obtained from the inversion of reflection measurements in a vacuum ultraviolet spectral range [15]. Moreover, the reflection and transmission of ultraviolet radiation by Ir thin films became a popular research topic [16]. Similarly, optical constants retrieved from reflectivity measurements for near-infrared, visible, and near-ultraviolet energies were reported in 1972 [17]. Weaver et al. used an extended spectral range [18,19,20] to obtain the dielectric function of bulk Ir from absorptivity and reflectivity measurements using a wide spectrum ranging from near-infrared to vacuum ultraviolet energies (0.2 to 40 eV). The optical constants of electron-beam evaporated Ir films were obtained from extreme ultraviolet to X-rays spectral ranges [21]. Refractive index and extinction coefficients of dc magnetron sputtered/atomic-layer-deposited Ir thin films were obtained using ellipsometry in the spectral range from the middle/near-infrared to vacuum/near ultraviolet [22,23]. Recently, the Ir optical constants in the visible, ultraviolet, and X-ray energy ranges, obtained from inversion of reflected electron energy loss spectra, have been published [24] (2 to 200 eV). Furthermore, the literature on band structure calculations for Ir is scarce. Characterizations of cross sections of the Fermi surface were initially reported [25,26], with band structure diagrams obtained using a Relativistic Augmented Plane Wave method [27]. Noffke and Fritsche [28] reported the total density of states (DOS) and the contributions from s, p, d, and f states, with similar results reported in more recent publications [29,30,31].

The primary purpose of the present study is to model the dielectric function (DF) of Ir for a polycrystalline bulk sample with known dielectric function obtained from absorptivity measurements at low temperature (for energies between 0.2 and 4.4 eV), and from reflectivity spectra at room temperature (RT) for energies between 3 and 40 eV [18,19,20]. The modeling was performed using an extended version of the Drude–Lorentz model, incorporating the contribution of electrons and holes in the Drude term, whose parameters are calculated from an optimization procedure based on a Spectral Projected Gradient Method (SPGM)-assisted Acceptance Probability Controlled Simulated Annealing (APCSA) approach. The same approach has been applied to Rh [32]. Using the optimized values of Drude parameters and the relative number density of conduction electrons, a set of derived physical quantities were evaluated for conduction electrons and holes, including effective masses, relaxation times, mobilities, Fermi energies, DOS at the Fermi energies, Fermi velocities, mean free paths, thermal conductivities, specific heats, magneton numbers, paramagnetic and diamagnetic susceptibilities, Hall’s coefficient, and hole and electron conductivities. We use our density-functional-theory (DFT) calculations, performed with the inclusion of the spin-orbit coupling, to interpret the Fermi energies and corresponding DOS values as background evaluations of the spectral dependence of the DOS and the band structure diagram for Ir. The resonance frequencies involved in the Lorentz contribution to the dielectric function were initially estimated from reported photometric studies. Because most spectral measurements were recorded at RT, we considered that the dielectric function, optimized parameters, and those that are derived corresponded to their values at RT.

Section 2 summarizes the formulation of the model, describing the parameters to be optimized and those obtained as derived quantities. Section 3 reports the characterization of Ir from the point of view of band structure calculations. Section 4 presents the modeling of the dielectric function of the sample previously mentioned whose dielectric function was obtained from the literature. Section 5 analyzes the approach from the band structure calculations to estimate the density of states and Fermi energies of itinerant electrons and holes charge carriers to compare with the evaluations based on the SPGM-APCSA optimization. Decompositions of bulk and surface energy loss functions (ELFs) are considered in Section 6. Evaluations of effective numbers of charge-carriers and electrons participating in each interband transition are reported in Section 7. Section 8 describes the decomposition of the Lorentzian ELF performed to estimate the density of states at the final energies corresponding to each set of transitions. A final summary and conclusions are given in the last Section.

2. Dielectric Function Model and Optimization Parameters

An extended version of the Drude–Lorentz (DL) model is obtained when we consider the contribution of conduction electrons and holes in the Drude term of the DF (ε = ε1 + iε2) and incorporate the polarization and absorption due to bound electrons [33,34]:

\[ \varepsilon=\varepsilon_{h f}-\frac{1}{\Omega}\left(\frac{\Omega_{p e}^{2}}{\Omega+i \gamma_{o e}}+\frac{\Omega_{p h}^{2}}{\Omega+i \gamma_{o h}}\right)+\frac{1}{z_{e}} \sum_{j=1}^{K} \frac{f_{j} \Omega_{p e}^{2}}{\Omega_{j}^{2}-\Omega^{2}-i \Omega \gamma_{j}}, \tag{1} \]

where εhf is the high-frequency dielectric constant associated with the high energy interband transitions spectrally located beyond the energy range considered; Ω = hω/2π = hν, where ν (ω) is the frequency (angular frequency) of the incident light; Ωpe = hνpeph = hνph), with ωpe = 2πνpeph = 2πνph) as the volume plasma frequency of the collective oscillation of electrons (holes) whose scattering frequency is γoeoh), with h as Planck's constant; and ωpe = (nee2/εome)1/2 with εo as the free space permittivity, e as the electron’s charge, ne as the number density of conduction electrons, and me as their effective mass. The scattering frequencies, γoe and γoh for electrons and holes, respectively, are related to corresponding relaxation times: τoe = h/2πγoe and τoh = h/2πγoh. The relative number density of conduction electrons, i.e., the number of conduction electrons per atom, is ze = ne/N with N as the number density of metal ions. N can be calculated from a, the lattice constant, and the number of atoms per conventional unit cell. Iridium has a face-centered cubic (fcc) structure with a = 0.3839 nm at RT [35], with four atoms per conventional unit cell [N = 4/Vcc with Vcc = a3 as the volume of the conventional unit cell]. It gives N = 7.07 × 1022 atoms/cm3 for Ir.

The formulation of DF given by Equation (1) has the following advantages with respect to the traditional form of applying DL models: (1) the energy dependence of the Drude term, denoted as D(ε), is self-consistent with the physical behavior required when the angular frequency tends to zero. It allows to obtain the static conductivity and corresponding electrical resistivity once the optimization is performed with a non-null convergent behavior for Im [D(ε)]; (2) the contribution to itinerant holes is accounted for in D(ε), allowing to use the same approach to describe the dielectric function of semiconductors; (3) all parameters involved in the formulation of DF have known physical interpretation consistent with the role played in Equation (1); and (4) the number of parameters depends on the specific material and the spectral range being considered. The number of parameters is not selected under the criterion that it is the minimum requirement to access a minimum of the objective function described later. This number is selected from previous information available in the literature about conduction properties (electrons and/or holes), spectrophotometric measurements, energy loss spectra, and band structure calculations. These four features allow the present formulation to be applied even in the case of large spectral ranges and materials with significant spectral structure in the DF, as is the case of Ir.

Within the context of the optimization method described below, the following three substitutions are done in Equation (1): γoh = ηγoe, Ωph2 = βχΩpe2, and εhf = 1/fo with β = me/mh, χ = nh/ne = zh/ze, and fo with values between zero and unity. The counterparts of ne and ze are nh and zh for holes, respectively, with ne different from nh because Ir is a non-compensated metal [36]. A frequency-dependent relaxation time was assumed to improve the fitting of infrared optical properties of certain metals using the Drude model [37,38,39]. We assumed the same dependence on angular frequency for both electrons and holes: 1/τe = 1/τoe + Λeω2 and 1/τh = 1/τoh + Λhω2, respectively. The corresponding scattering frequencies were γe = γoe [1 + (Ω/Γe)2] and γh = γoh [1 + (Ω/Γh)2], where Γe and Γh are two additional parameters to be optimized: Γk = (hγok/2πΛk)1/2 with the subscript letter k = e for electrons, and k = h for holes. These values for γe and γh substitute those of γoe and γoh in Equation (1), respectively. The seven parameters to be optimized in DF’s Drude term [εD] are Ωpe, γoe, β, χ, η, Γe, and Γh, which determine the contribution of conduction electrons and holes to the dielectric function of the metal. Additionally, the background polarization due to interband transitions beyond the spectral range considered is accounted for as the inverse of the optimized value of fo. In σo = εoωε2(ω), when ω tends to zero, there are no interband transitions contribution in this limit.

The contribution of interband excitations to the dielectric function is approximated by the classical Lorentz summation term in Equation (1), where these excitations are modeled by resonance oscillators, each one associated with a jth population of electrons bound to the metal ions whose resonance frequency (energy) is ωjj), i.e., εL = εL,1 + εL,2 +…+ εL,K with εL,j(Ω) = fjΩpe2/{ze [(Ωj22)-iΩγj]}. The index j goes from 1 to K, with K equal to the number of oscillators considered. In quantum mechanics, the resonance frequencies correspond to electronic transitions between the occupied and vacant states close to critical points in the band structure. The resonance energies (Ωj), Lorentzian widths (γj), and oscillator strengths (fj) appear in the Lorentz term. The light absorption by bound electrons is determined by the imaginary component of the Lorentz term in the DF, which is proportional to fjγj. The Np = 9 + 3K parameters in the DL model of the dielectric function are initially optimized using an Acceptance-Probability-Controlled Simulated Annealing (APCSA) approach [40,41]. We implemented our APCSA code with certain improvements as described previously [33,34]. In addition, we minimize the following merit function using the APCSA method

\[ F=\frac{1}{2 \kappa-N_{p}-1} \sum_{i=1}^{\kappa}\left[\left(\frac{\varepsilon_{_1}\left(\omega_{i}\right)-\bar{\varepsilon}_{1}\left(\omega_{i}\right)}{\bar{\varepsilon}_{1}\left(\omega_{i}\right)}\right)^{2}+\left(\frac{\varepsilon_{2}\left(\omega_{i}\right)-\bar{\varepsilon}_{2}\left(\omega_{i}\right)}{\bar{\varepsilon}_{2}\left(\omega_{i}\right)}\right)^{2}\right], \tag{2} \]

where κ is the number of spectral points where the dielectric function’s values are known, i.e., $\bar{\varepsilon}\left(\omega_{i}\right)=\bar{\varepsilon}_{1}\left(\omega_{i}\right)+i \bar{\varepsilon}_{2}\left(\omega_{i}\right)$. The number of degrees of freedom, i.e., the number of terms that can change in the minimization of F is 2κ minus the number of quantities that are being optimized, including F itself. This definition of the merit function allows us to simultaneously retrieve both ε1(ω) and ε2(ω) with minimal error throughout the spectral range considered. The optimizations were conducted until the solidification condition was satisfied: the lowest F-value for each temperature cycle was compared with the lowest values of the three preceding consecutive temperatures [40]. When the corresponding three relative differences are lower than the specified numerical tolerance Δ, the simulated annealing process is finished. We used Δ = 10-4. In the APCSA approach, the temperature is a control parameter with the same units the merit function has [42], which is updated from the dispersion in the merit function or accepted energy values once each cycle is finished. After this stage of optimization, the optimized values were used to repeat the process but with the application of an SPGM after each temperature cycle. We call this an SPGM-assisted APCSA approach applied to improve the solution. Details of the SPGM have been published elsewhere [43,44]. The Appendix includes additional information on the components of the gradient of F obtained from the explicit expressions of ε1 and ε2.

Once the optimization is performed, a set of derived quantities were evaluated, as shown in Table 1, demonstrating the fundamental nature of the dielectric function to all extent. We assumed that the holes contribute to the paramagnetic and diamagnetic susceptibilities in the same manner as by conduction electrons [45]. In the two-band model, Hall’s coefficient listed in Table 1 was written as RH = -(1 - χq2)/[ezeN(1 + χq)2] with q = β/η [46].

Table 1 List of physical quantities evaluated once the optimization of the DL parameters is finished.

Both electrons and holes contribute to the intrinsic conductivity in fractions given by Pe = 1/[1 + χq] and Ph = 1 - Pe, respectively. In addition to the spectral dependence of ε2(ω), the absorption of light was considered through the dynamic conductivity [σ(ω) = εoωε2(ω)] and the energy loss function as described in Section 6.

To compare with previous reports of Drude parameters, the added Drude dielectric susceptibilities corresponding to the second and third terms in the right side of Equation (1) [χDrude = -Ωpe2/(Ω(Ω + iγe) - Ωph2/(Ω(Ω + iγh)] can be expressed approximately as χeff = -Ωo2/[Ω(Ω + iγeff)]. The effective values of γeff are obtained from γeff = -Ω·Im(χeff)/Re(χeff). As expected, this effective scattering frequency shows a dependence on energy that can be modeled by γeff = γo(1 + (Ω/Γo)2), with regression coefficients close to unity. The fitting of γeff allows us to calculate both γo (τo = h/2πγo) and Γo. This approach to obtaining the γeff- and Γo-values holds when the Γe and Γh, as well as γoe and γoh, have similar orders of magnitude, i.e., when their ratios are about between unity and ten. Otherwise, the frequency-dependence of the effective scattering frequency must be fitted from some approximation of γeff = [γe2 + γh2) + χγh2 + γe2)]/[(Ω2 + γh2) + χ(Ω2 + γe2)] and depending on the orders of magnitudes of γoeoh and Γeh. The value of effective volume plasma energy Ωo was obtained for each spectral point from the relation Ωo = [(-Re(χeff) + Im(χeff))Ω2(1 + α2)/(1 + α)]1/2, with α = γeff/Ω, showing very low dispersion. Its average value was evaluated from the set of Ωo values. The γo, Γo, and Ωo parameters were considered as effective Drude quantities.

3. Band Structure Calculations Based on Density Functional Theory (DFT)

We have initially performed DFT calculations, with no inclusion of spin-orbit coupling (SOC), using a norm-conservative pseudopotential [48] and a dense k-point mesh for the band diagram and DOS spectra by solving the Kohn–Sham equation [49], with a plane wave basis set to calculate the band structure diagram (BSD), charge distribution and projected density of states (DOS) of Ir through the Quantum Espresso computational suite [50] (Figure 1). An exchange–correlation functional within the local density approximation was used as devised by Perdew, Becke, and Ernzerhof (PBE) [51,52], using the Troullier–Martins procedure to generate norm-conserving pseudopotentials [53]. The norm-conservative pseudopotential was used with nine valence electrons per atom of Ir in the primitive unit cell. Unit cell relaxation provides us with a lattice constant of 0.389 nm, close to the experimental value of 0.384 nm [35,54]. Thus, we showed the pertinence of the pseudopotential and exchange-correlation functional choice for our calculations.

Click to view original image

Figure 1 (a) Electronic band structure diagram (BSD) of Ir evaluated by DFT with no spin-orbit coupling. Projected DOS onto orthogonalized atomic orbitals resolved in momentum-space for (b) hybridized sp-orbitals, and (c) d-orbitals. (d) Total and projected DOS summed over all k-points. Values in the horizontal axis are given in states/eV atom. The vertical scale is given in eV, with the horizontal 0 energy indicating the standard Fermi level EF (SFL).

The overestimation of the lattice constant is characterized by DFT-based methods using the PBE exchange-correlation functional [55]. The wavefunction cutoff was set to 1360 eV. The self-consistent calculation was achieved with an accuracy of 1.4 × 10-5 eV. We used the Methfessel–Paxton smearing function [56] to account for the metallicity of Ir, with a width of 0.272 eV and 10 × 10 × 10 mesh k-points to sample the Brillouin zone [57]. We used 500 k-points to obtain the BSD. Atomic positions were relaxed until the converged successive values of the pseudopotentials corresponded to relative differences lower than 10-4. Figure 1(d) shows the DOS of Ir, as evaluated from DFT with no SOC. Similar spectral features have been reported by Noffke and Fritsche, who calculated the DOS at the Fermi level, ρ(EF), in about 1.04 states/eV atom (Figure 1 in [28]). Our evaluations provide around 0.84 states/eV atom with no spin-orbit coupling and 1.06 states/eV atom when the spin-orbit is considered (Figure 2). The values 0.91 and 1.01 states/eV atom have been reported for ρ(EF) [29,31]. The DOS is predominantly determined by the contribution of d-states (Figure 1(c)), with the Fermi energy just below the largest peak of the DOS at energies above the standard Fermi level (SFL). The width for the d-band was close to 10.6 eV, i.e., the d-band extended from about 8.9 eV below the Fermi level to 1.7 eV above. The SFL, with respect to the bottom of the hybridized sp band, was around EF = 10 eV. A heavy d-hole-like sheet extending from the W- to the X- symmetry point, at energies between 1.29 and 1.41 eV above the SFL, and light d-hole-like pockets centered around the X-point, at 0.52 and 0.92 eV above the Fermi level, are displayed in the BSD of Figure 1(a) and Figure 1(c). A similar finding has been reported for Pd [58]. Another band of empty d-states just above the Fermi level was displayed in the L-W direction with energies between 0.42 and 1.29 eV. The presence of these d-hole states correlated with peaks in the projected density of d-states, centered at 0.16 and 1.38 eV just above the SFL, whose maximum values are 0.91 and 0.93 states/eV atom when j = 5/2, respectively, 0.51 and 0.34 states/eV atom when j = 3/2, respectively.

Click to view original image

Figure 2 (a) Electronic BSD of Ir evaluated by DFT with the inclusion of spin-orbit coupling. Projected DOS onto orthogonalized atomic orbitals resolved in momentum-space for (b) hybridized sp-orbitals and (c) d-orbitals. (d) Total and projected DOS summed over all k-points. The values in the horizontal axis are given in states/eV atom. The vertical scale is given in eV, with the horizontal 0 energy indicating the standard Fermi level. The red arrows in the left figure are indicated for use in the discussion.

Itinerant s- or sp-electrons contribute to the electrical conduction in Ir undergoing transitions up to empty d-states of neighboring atoms, just above EF. The three red arrows in Figure 2(a) indicate the plausible transitions of these electrons. This issue is discussed in detail in Section 5. The inclusion of SOC in DFT calculations leads to a more structured DOS, as seen in Figure 2(d). Regarding the BSD, the major change is the splitting of the d-band displayed in the Γ-L direction, which crosses the SFL. At the Γ-symmetry point, the lower d-band displayed in the energy range showed a splitting due to the SOC, which projected toward the Γ-Δ and Γ-L-W directions.

Once the SPGM–APCSA optimization is achieved, the density of states and Fermi energies for itinerant conduction holes and electrons as derived parameters were used and compared these quantum mechanics parameters with corresponding ones were obtained from the BSD. Thus, we will be able to distinguish the transitions involved in the transport of itinerant charge carriers.

4. Dielectric Function of Bulk Iridium

The DF of Ir was obtained by Weaver et al. for a large specular polycrystalline sample normally illuminated [19] (read lines in Figure 3(a)). The purity of the sample was unspecified. The spectral absorptivity measurements (A) at low temperature (4.2 K) were used to obtain the reflectivity spectrum (R = 1 - A) in the energy range of 0.2 to 4.4 eV, and it was directly measured at RT in the energy range of 3.0 to 40 eV. A unique reflectivity spectrum was considered in the spectral range of 0.20 to 40 eV, with good agreement in the region of the overlap of the two reflectivity spectra previously measured. The Kramers–Kronig analysis was used to obtain the DF from inversion of the reflectivity spectrum [20]. The analysis by Weaver et al. includes a report of the spectral dependence of the dynamic conductivity, which is proportional to $\bar\varepsilon _2$, the imaginary part of the retrieved DF (blue line in Figure 3(b)). The following resonance energies, reported in eV, are: Ω2 = 1.05, Ω3 = 1.9, Ω4 = 3.1, Ω5 = 4.12 [19]. In the same study, the imaginary component of the DF shows two additional higher resonance energies at Ω6 = 10.8 and Ω7 = 18.8 eV. Finally, the experimental energy loss function, -Im [$1/ \bar\varepsilon(\omega) $], depicts a resonance energy at Ω8 = 33 eV. These authors comment about a low energy resonance in Ω1 = 0.40 eV reported in [17].

Click to view original image

Figure 3 (a) RT real (ε1) and imaginary (ε2) components of the DF of Ir as fitted from an SPGM-assisted APCSA optimization with eight resonance energies in the Lorentz term (blue lines). Red lines are traced from the experimental values of the DF [20]. Values of ε1 are negative in the bluish spectral ranges. (b) Spectral dependence of the dynamic conductivity, with contributions from charge–carriers (Drude) and interband transitions (Lorentz).

The APCSA optimization was performed by assuming small variations in the energy resonances around the set of indicated Ωj-values, with κ = 200 spectral points and K = 8 resonance energies. The specified number of spectral points was generated by cubic spline interpolation from the 124 data reported previously [20]. The other parameters were optimized by searching in ranges with mobile boundaries, requiring 106 temperature cycles to satisfy the solidification condition. At the end of the first stage of the APCSA optimization, the value of the merit function was F = 0.027. With the application of the SPGM-assisted APCSA strategy, the final merit function reached the value F = 0.010.

4.1 Optimized Transport Parameters

The modeling of the DF was displayed in Figure 3(a), with optimized values of the parameters specified in Table 2 and derived quantities in Table 3. Figure 3(b) shows the spectral dependence of the dynamic conductivity, with contributions from the Lorentz (interband transitions) and Drude (charge-carriers transport) terms to the imaginary component of the DF, whose spectral details resemble those reported in [19] for the total conductivity. The optimized fo-value close to unity implies that the contribution to the dielectric function by interband transitions beyond 40 eV was small (εhf = 1.06), i.e., there existed a small background polarization associated with interband transitions at energies beyond the spectral range considered.

Table 2 SPGM-assisted APCSA optimized DL parameters used to calculate the RT dielectric function of a polycrystalline annealed bulk sample of Ir [20]. Parameters are given in eV, except the dimensionless ze, β, η, χ, fo, and the oscillator strength fj-values. Values in bold correspond to Drude parameters, the rest are the Lorentz parameters, and F is the optimized value of the merit function.

Table 3 Set of derived quantities for the Ir bulk sample whose room temperature DF is depicted in Figure 3(a) with the optimized DL parameters reported in Table 2. Quantities in bold are the effective Drude parameters.

The β-value was close to 8/3, meaning that the average effective mass of conduction electrons was larger than that of the holes (me = 1.64 m and mh = 0.61 m). Similar values of the scattering frequencies of conduction electrons and holes were obtained according to the value close to unity for the η-parameter. Scattering frequencies, γoe, and γoh, were determined by electron-electron and hole-hole interactions, respectively, as well as interactions of charge carriers with phonons. According to the optimized values of Γe and Γh, the frequency dependence of the relaxation times was negligible (Λe = 4.70 and Λh = 4.80 as). This is expected when considering the annealing process applied to the sample used by Weaver et al. [19]. The number density of electrons was larger than that of holes (χ = 0.76), with ze = 0.33. A similar value was reported for Pd, a transition metal characterized by ze = 0.36 [59]. Similar values were reported for other transition metals: 0.38 for Ti, 0.40 for Zr, and 0.44 for Hf [60].

After subtracting the residual contribution, the experimental resistivity values of Ir (ρexp) were reported between 0.005 µΩ cm at low temperatures and 5.0 µΩ cm at room temperature [61]. The value obtained as a quantity derived from the optimization process (ρo = 11.6 µΩ cm) agreed with the order of magnitude for the corresponding static conductivity as extrapolated from the spectral dependence of the dynamic conductivity reported by Weaver et al. [19]. Observation ρo > ρexp suggests the presence of impurities in the sample prepared by Weaver et al. The Hall’s coefficient (RH = 12.9 × 10-11 m3/As) is larger than the values reported for samples at low- and room-temperatures, namely, 3.49 × 10-11 m3/As at 81 K and 3.18 × 10-11 m3/As at 300 K, respectively [62]. These two parameters, ρo and RH, are sensitive to the purity of the sample. The effect of impurities is used in evaluating the parameters through the intrinsic relaxation times of charge carriers. The thermal conductivity appears to be affected by impurities. Its derived value, κth = 67.3 W/mK at T = 300 K, was evaluated using the Lorenz number reported for Ir: L = 2.60 × 10-8 WΩ/K2 [63]. The reported value in the literature is around 150 W/mK [64].

The large values of resistivity and Hall’s coefficient, as compared with the reported data, indicate a significant content of impurities in the sample used by Weaver et al. In addition, the light diffusively reflected by the surface of the sample corresponds to the absorption by the medium when normalizing reflectivity measurements. Although the number density of holes (zh = 0.25) is less, they primarily contribute to the conduction in Ir, as shown by σe and σh values (Pe = 33% and Ph = 67%). The relative effective mass of holes was determined by Hörnfeldt et al. using de Haas-van Alphen oscillation measurements [65], with values between 0.23 and 1.28 depending on the crystallographic plane and specific hole pocket considered. Our calculated average value (mh/m = 0.61) was within this range. Grodski and Dixon have reported values for the relative effective mass between 1.34 and 3.02 for conduction electrons with states in the Γ-point [66]. In addition, the average value calculated from the optimization process (me/m = 1.64) is within this range.

In addition to the effect of impurities and surface roughness, there is another aspect that significantly influences the values of the transport parameters. For example, thermal conductivity value, which is proportional to the electrical conductivity, is sensitive to the limiting behavior of the imaginary component of the DF when the frequency tends to zero. At low energies, the dielectric function was obtained from absorptivity measurements recorded at low temperatures. From the visible to the ultraviolet, measurements were recorded at RT. The net effect of these three factors was reflected in a thermal conductivity derived from the optimization whose value is around 50% lower than the measurements reported in the literature.

4.2 Magnetic Susceptibility and Enhancement by Charge Carrier Exchange Correlation

The reported magnetic susceptibility of Ir is χexp = 0.133 × 10-6 cm3/g at room temperature [67]. A higher value is expected at low temperatures, which is consistent with our calculation derived from the SPGM-assisted APCSA optimization: χm = 0.53 × 10-6 cm3/g. The total molar susceptibility is obtained from the total mass susceptibility as follows: χM = Mχm = 1.02 × 10-4 cm3/mol, with M = 192.217 g/mol as the atomic weight of Ir, and a similar relation for the diamagnetic susceptibilities. The fact that χm > 0.133 × 10-6 cm3/g (the experimental value) is an indication of the presence of traces of ferromagnetic elements in the sample studied by Weaver et al. This also affects the enhancement of the magnetic susceptibility due to electron–phonon interactions. According to Curie’s law for paramagnetic materials and the relation between diamagnetic and paramagnetic susceptibilities, the number of effective Bohr magnetons is given by p = (3 kBoNav)½MT)½B with Nav as Avogadro’s number. The average magnetic moment per atom is µavg = B. The evaluation for Ir results in p = 0.10 at room temperature. If the exchange–correlation effect is included, χM,ec = χM/[1-ρ(EF)IF] where IF is the exchange–correlation integral and ρ(EF)IF is the Stoner parameter [68,69]. For Ir, IF has been reported by Sigalas et al. as 0.287 eV, and with ρ(EF) = 1.01 states/eV atom [31], the evaluations results in χM,ec = 1.44 × 10-4 cm3/mol and p = 0.12.

The electrical conduction is attributed to itinerant electrons and holes charge carriers, implying that these charge carriers partially behave like free charge carriers or like localized ones. The presence of this partial localization or its spatial confinement decreases the magnetic susceptibility, as shown by Cahaya, who extended the known relations of Pauli and Landau susceptibilities by incorporating the following two facts in his derivation [70]: (a) by considering an infinite quantum well, he assumed null values of the charge carrier wave functions at the boundaries of the medium instead of periodic boundary conditions, and (b) an artificial local confinement was introduced through a spatially periodic magnetic field whose spatial period λq was tuned through the wavenumber q = 2π/λq. This field introduced local confinement of charge carriers in those planes perpendicular to the direction of the magnetic field. He obtained the same expression for Pauli and Landau magnetic susceptibilities, times corresponding factors given by

\[ \mathrm{Z}_{P, k}=\frac{1}{2}+\frac{1-x_{k}^{2}}{4 x_{k}} ln \left|\frac{x_{k}+1}{x_{k}-1}\right|, \tag{3a} \]

\[ \mathrm{Z}_{L, k}=\frac{3}{8 x_{k}^{2}}\left[1+x_{k}^{2}-\frac{\left(1-x_{k}^{2}\right)^{2}}{2 x_{k}} ln \left|\frac{x_{k}+1}{x_{k}-1}\right|\right]=\frac{3}{4 x_{k}^{2}}\left[1-\left(1-x_{k}^{2}\right) Z_{P, k}\right], \tag{3b} \]

where xk = q/2kF,k, with kF,k as the Fermi momentum of the kth charge carrier type. These two parameters tend to unity when q tends to zero, i.e., the known expressions of paramagnetic and diamagnetic susceptibilities are recovered as the local confinement disappears. Otherwise, the values of ZP and ZL are lower than unity, i.e., they are diminution factors, ZP,k being the Lindhard function [71]. In the absence of an external magnetic field, other effects contribute to a spatially periodic or quasi-periodic internal magnetic field: thermal fluctuations of atomic nuclei around their equilibrium positions [72] and small variations in the quasi-stationary electrostatic screening potential due to thermal fluctuations in the position of the center of mass of the electronic gas of charge carriers [73].

To apply Cahaya’s formalism linked to the initially calculated magnetic Pauli and Landau susceptibility components for both electrons and holes, we wrote the corresponding mean free paths in terms of the spatial period of the magnetic field: Lk = NP/L,kλq,P/Z,k and xk becomes xk = πNP/L,k/LkkF,k with kF,k = (3π2nk)1/3 [kF,e = 8.89 nm, kF,h = 8.11 nm]. Table 4 shows the results. The fraction of each contribution to the total susceptibility (FP,k and FL,k with k = e and h) is calculated from the results of the optimization, after applying the enhancement due to exchange correlations.

Table 4 Diminution factors for both conduction electrons and holes contributing to the Pauli and Landau magnetic susceptibilities. The lattice constant is a = 0.389 nm.

The four contributions giving the measured susceptibility are evaluated by assuming that the same fractions can be associated with the reported experimental values. For each component, the equation │FP/L,k·χexp - ZP/L,k·χk│ = ε, with ε close to or lower than 10-10, has been solved for χ1 = χe,ec, χ2 = χh,ec, χ3 = χe,ec', and χ4 = χh,ec'.

4.3 Charge-Carrier Heat Capacity and the Coupling Constant

The reported value of the charge–carrier heat capacity coefficient is γexp = 3.14 mJ/K2 mol [74], which includes the contributions from electron-phonon coupling (EPC), and exchange correlation (EC) between charge carriers. The value obtained from the optimization process was significantly lower, γhcc = 1.72 mJ/K2 mol, about half of the reported one, implying that the contributions from the EPC and EC are large in Ir (see Table 4 in [64]). The total heat capacity coefficient γ, is given by the relation γ = π2 kB2ρ(EF)(1 + λ)/3 [75] with contributions from charge carriers and electron-phonon interactions. From γ = γexp = γhcc(1 + λ) with λ as the coupling strength of the charge carriers both to the lattice phonons and to other charge carriers [76], λ = (γexphcc)/γhcc = 0.83 is obtained.

These evaluations do not account for the effect of EC between charge carriers at the Fermi surface neither EPC in the value of γhcc obtained from the optimization. McMillan estimated λ for Ir as 0.34 [75]. The electronic specific heat is given by C = γT, where T is the absolute temperature, and γ is the specific heat coefficient, which can also be represented as γ = (πkB)2nm/(hkF/2π),2 where kF = (3π2n)1/3 is the Fermi momentum, n is the number density of charge carriers, and m is the effective mass of each one. If the previous equation for γ is applied to an ideal electron gas, with no interactions between charge carriers or between them and the phonons, γ $\to$ γo,k = (πkB)2nkm/(hkF,k/2π)2 with k = e or h. In the context of the present optimization approach, γk = (πkB)2nkmk/(hkF,k/2π)2. Consequently, γko,k = mk/m. This relation is valid in the limit of T $\to$ 0, where the Fermi momentum remains unchanged by many-body effects according to the Luttinger theorem [77,78]. The optimized values of relative effective masses of electrons and holes become enhancement (λk > 0) or diminution (λk < 0) factors [mk/m ≡ 1 + λk] that must be applied to the first estimations of the corresponding specific heat coefficients. The corrected values are γhcc,e $\to$ (1 + λe)·γhcc,e = 2.10 mJ/K2 mol and γhcc,h $\to$ (1 + λh)·γhcc,h = 0.27 mJ/K2 mol. The recalculated total heat capacity coefficient is γhcc = 2.37 mJ/K2 mol. The evaluation of Papaconstantopoulos, based on tight-binding calculations, results in 2.20 mJ/K2 mol [79]. The corrected value of the electron-phonon coupling parameter is λ = 0.32, which is close to that reported by McMillan [75].

4.4 Optimized Interband Energy Resonances

The optimized values of resonance energies associated with interband transitions are close to those initially used as restricted values in the first APCSA optimization. Once the second stage of optimization was performed, Ω1 changes from 0.40 to 0.42 eV with a 5% variation with respect to the initial value [Ω1: 0.40 $\to$ 0.42 (5%)]. For other resonances, the following variations existed: Ω2: 1.05 $\to$ 1.21 (15%); Ω3: 1.90 $\to$ 2.07 (9%); Ω4: 3.10 $\to$ 3.22 (4%); Ω5: 4.12 $\to$ 4.38 (6%); Ω6: 10.8 $\to$ 11.1 (3%); Ω7: 18.8 $\to$ 19.4 (3%); and Ω8: 33.0 $\to$ 31.2 (5%). These optimized resonance energies must be considered effective because each one could correspond to convolution of several sets of transitions between bands, with similar energy gaps between initial and final bands. We will again consider these resonance energy values when calculating the densities of states of final states involved in these transitions, according to the formulation given in Section 8.

5. Densities of States and Fermi Energies of Conduction Holes and Electrons

The values of the density of states of itinerant conduction electrons and the corresponding Fermi energy, derived from the SPGM-assisted APCSA optimization, are ρ(Ee) = 0.27 states/eV atom and Ee = 1.84 eV, respectively. These values can be approached from the band structure diagram in the following manner: the contribution to electrical conduction in metal transitions is associated with electrons transitioning from s- to d-bands [80]. Figure 4(a) depicts a zoomed image of the density of states displayed in Figure 2(d) and shows in detail the density of states of s-electrons below the SFL. The composition of this s-DOS is highlighted with number labels for the nine peaks displayed. Each jth peak can be modeled as a Gaussian function of amplitude Hj and deviation σj. The composition of the fourth peak has not been resolved. It will be considered as a single peak in the following analysis. The area of each peak was Aj = (2π)1/2σjHj. Each σj-value was obtained from that point where the height of the peak decreased to the value Hje-1/2. The areas obtained are used to calculate the weighted average for the density of states and the Fermi energy. Those transitions corresponding to the largest peaks closer to the SFL contribute more significantly to these averages.

Click to view original image

Figure 4 (a) DOS below the standard Fermi level for s-states whose composition has been highlighted with number labels for nine peaks displayed. (b) DOS above the standard Fermi level for d-states (blue and red colors for j = 5/2 and j = 3/2, respectively). The four peaks initially considered in the analysis are numbered.

Therefore, we considered the third, fourth, and seventh peaks with the following values: σ3 = 0.29, σ4 = 0.42, and σ7 = 0.30 eV; H3 = 0.21, H4 = 0.31, and H7 = 0.15 eV-1. The areas are A3 = 0.15, A4 = 0.33, and A7 = 0.11. The weighted average value of the DOS of these itinerant conduction electrons was <ρ(Ee)> = (A3H3 + A4H4 + A7H7)/(A3 + A4 + A7) = 0.25 states/eV atom. The average Fermi energy is <Ee> = Nσ(A3σ3 + A4σ4 + A7σ7)/(A3 + A4 + A7) = 1.82 eV with Nσ = 5 in this case. Given the position on the energy of each peak, Ej, the weighted average of the Fermi energy was calculated from the energy range around each peak, which was approximated from Ej - Nσσj/2 to Ej + Nσσj/2. This explains the factor Nσ used to evaluate the average Fermi energy. The drawback of this method was the use of Nσ as a fitting parameter. When its value was close to 3, it would be interpreted as an indication of the presence of a Gaussian peak. When Nσ takes a value close to 10, it implies that the tails of the peak decrease more slowly than those of Gaussian profiles. In this case, a Lorentzian behavior is expected for energies lower than Ej - σj and larger than Ej + σj. Another method to interpret the information provided by the band structure diagram is to consider only the dominant peak. In this case, A4 = 0.33 is equal to the average number of electrons per atom participating in the transitions as itinerant charge carriers. With <ρ(Ee)>ΔE4 = A4, one has that <ρ(Ee)> = A4/Nσσ4 = 0.26 states/eV atom with Nσ = 3. The corresponding Fermi energy can be obtained from <Ee> = 3A4/[2<ρ(Ee)>] = 1.90 eV. These weighted averages are close to those obtained from the APCSA optimization. The correlation between these three peaks considered and the band structure diagram shown in Figure 4(a) can be used to discriminate the transitions involved. The third peak correlated to the transitions at the Σ-point, between the two parallel bands: from 7.77 eV below to that standard Fermi level (see the right red arrow in Figure 2(a)). The dominant fourth peak was contributed by transitions from the same parallel bands; however, it also included transitions at the X-symmetry point, from the band 6.97 eV below the Fermi level up to empty states of the d-hole pockets (see the red arrow in the middle of Figure 2(a)). The seventh peak was determined by transitions close to the K-symmetry point, from 3.90 eV below the Fermi level up to the parallel band just above the SFL (see the left red arrow in Figure 2(a)).

For conduction holes, the derived values of the density of states and the corresponding Fermi energy are ρ(Eh) = 0.09 states/eV atom and Eh = 4.09 eV, respectively. The contribution of these charge carriers to conduction is attributed to intra-band transitions in the d-band, just above the Fermi level and involving the d-hole pockets and/or the d-hole sheet mentioned in Section 3. Figure 4(b) shows the d-DOS for two values of the quantum number j associated with the spin-orbit coupling (j = 3/2 and 5/2). In addition, we numbered those peaks as well resolved: the first two peaks correspond to the d-band with j = 5/2, and the other two peaks are displayed for the d-band with j = 3/2. After analyzing the contributions of each peak to the average density of states and the average number of holes participating in the intra-band transitions, the only alternative that approached the values obtained from the APCSA optimization corresponded to the fourth peak characterized by H4 = 0.34 states/eV atom and σ4 = 0.27 eV. The corresponding area is A4 = 0.23 andthe average density of states is <ρ(Eh)> = A4/Nσσ4 = 0.085 states/eV atom with Nσ = 10. The large value of Nσ could be attributed to a Lorentzian characteristic of this peak. The average Fermi energy is <Eh> = 3A4/[2<ρ(Eh)>] = 4.06 eV. The agreement with those values obtained from the SPGM-assisted APCSA optimization was satisfactory. The dominant intra-band d-hole transitions correlated with the fourth peak were those involving the d-hole-like sheet shown between X- and W-symmetry points, at about 1.30 eV above the SFL, toward the d-band displayed between the W- and L-symmetry points.

6. Bulk and Surface Energy Loss Functions

The bulk electron energy loss function can be written in the form $L^{(B)}(\Omega)=-{Im}[1 / \varepsilon(\Omega)]=-A_{o}^{(B)}(\Omega) {Im}\left[1 / \varepsilon_{D}(\Omega)\right]-B_{o}^{(B)}(\Omega) \operatorname{Im}\left[1 / \varepsilon_{L}(\Omega)\right]$, with εD and εL as the Drude and Lorentz contributions, respectively [81,82]. In this way, $L_{D}^{(B)}(\Omega)=-A_{o}^{(B)}(\Omega) \operatorname{Im}\left[1 / \varepsilon_{D}(\Omega)\right]$ and $L_{L}^{(B)}(\Omega)=-B_{o}^{(B)} \operatorname{Im}\left[1 / \varepsilon_{L}(\Omega)\right]$ are the bulk ELFs corresponding to the charge-carriers and bound electrons, respectively. The real values of $A^{(B)}_o$ and $B^{(B)}_o$ coefficients are obtained at each Ω-point by solving the previous equality for $Im[1 / \varepsilon({\mit\Omega})]$ along with that for the real component $\left\{\operatorname{Re}[1 / \varepsilon(\Omega)]=A_{o}^{(B)}(\Omega) \operatorname{Re}\left[1 / \varepsilon_{D}(\Omega)\right]+B_{o}^{(B)}(\Omega) \operatorname{Re}\left[1 / \varepsilon_{L}(\Omega)\right]\right\}$. This allows a system of two equations with two unknowns at each spectral point. A similar approach is followed to perform the decomposition of the surface ELF:$L^{(S)}(\Omega)=-\operatorname{Im}\{1 /[\varepsilon(\Omega)+1]\}$ in terms of $L_{D}^{(S)}(\Omega)=-A_{o}^{(S)}(\Omega) \operatorname{Im}\left\{1 /\left[\varepsilon_{D}(\Omega)+1\right]\right\}$ and $L_{L}^{(S)}(\Omega)=-B_{o}^{(S)}(\Omega) \operatorname{Im}\left\{1 /\left[\varepsilon_{L}(\Omega)+1\right]\right\}$ [83]. Figure 5(a) displays the results for the bulk ELF values, and Figure 5(b) corresponds to the surface ELF spectra. Based on the total bulk and surface ELF spectra, Weaver et al. assigned the low-energy peaks to excitations of volume and surface plasmons spectrally located at 7.8 and 7.2 eV, respectively [19]. From our deconvolution, the corresponding resonance energies of these plasmons were calculated to be 7.6 and 7.0 eV (see bold labels 1 in Figure 5). From the SPGM-assisted APCSA, the effective volume plasmon energy corresponding to its resonance frequency was calculated as Ωo = 7.7 eV. Other two couples of volume- and surface-plasmons are displayed: the first at 14.0 and 13.3 eV (see bold labels 2) and the second at 36.4 and 26.2 eV (see bold labels 3). In addition, Weaver et al. interpreted that the largest peak of the total bulk ELF was due to the superposition of a plasmon excitation and absorption due to interband transitions. Recently measured ELF of Ir showed absorption peaks at 8.9, 15.4, and 31.0 eV, with the other two peaks at 54.2 and 67.0 eV [24]. The decomposition performed showed the spectral superposition of both plasmon excitations and electronic transitions. It is difficult to distinguish the true characteristic of the peaks just from the total ELF spectra, as stated by Weaver et al. Additional low-energy volume and surface plasmons are displayed at energies of 1.3 and 0.8 eV, respectively. These oscillations are probably attributed to the hybridization of collective excitations of bound electrons participating in interband transitions that accounted for the first three oscillators. In addition, a surface plasmon was obtained at 3.9 eV. Its corresponding volume counterpart was not resolved due to the superposition with the tail of the volume excitation 1 at energies lower than 7.6 eV.

Click to view original image

Figure 5 Bulk (a) and surface (b) energy loss function (ELF) for the Ir sample whose modeled DF is displayed in Figure 3(a). The subscripts D and L are used to specify the ELF for the charge-carriers (Drude) and bound electrons (Lorentz). Each spectrum is dominated by the presence of three peaks, as indicated by the labels of the figures.

Bound electrons related to the fourth and fifth oscillators could be involved in the collective oscillations and correlated with these low-energy peaks of the bulk and surface ELFs. The low-energy structure in the ELF of Ir was slightly observed from the analysis reported in [84], with the predominant peaks at 7.9, 14.8, and 32.2 eV.

7. Effective Numbers of Electrons and Densities of States

Considering the bulk ELF, the effective total number of bulk electrons and holes per atom participating as charge carriers and bound electrons per atom involved in interband transitions in the spectral range considered can be evaluated from the partial sum rule [85]

\[ Z_{e f f}=\frac{2}{\pi \Omega_{p e}^{2}} \int\limits_{0}^{\Omega_{\max }} \Omega \cdot \operatorname{Im}\left[-\frac{1}{\varepsilon(\Omega)}\right] d \Omega=\frac{2}{\pi \Omega_{p e}^{2}} \int\limits_{0}^{\Omega_{\max }} \Omega \cdot L^{(B)}(\Omega) d \Omega, \tag{4} \]

with Ωmax = 40 eV in our case. The integration is performed with the Simpson rule and Richardson extrapolation [86]. The optimized parametric form of the DF is used to extrapolate $L^{(B)}(\Omega)$from Ωmin = 0.20 eV to zero frequency. The evaluation was carried out as Zeff = 25.48. Similar sum rules hold for the charge-carriers and bound electrons per atom. They can be evaluated using the decomposition of the ELF to obtain the effective number of charge-carriers contributing to conduction [$Z_{e f f}^{(c c)}$], as well as the effective number of bound electrons participating in interband transitions [$Z_{e f f}^{(i t)}$]:

\[ Z_{e f f}^{(c c)}=\frac{2}{\pi \Omega_{p e}^{2}} \int\limits_{0}^{\Omega_{\max }} \Omega \cdot L_{D}^{(B)}(\Omega) d \Omega, \tag{5a} \]

\[ Z_{e f f}^{(i t)}=\frac{2}{\pi \Omega_{p e}^{2}} \int\limits_{0}^{\Omega_{\max }} \Omega \cdot L_{L}^{(B)}(\Omega) d \Omega. \tag{5b} \]

Furthermore, the ELF associated with the bulk Lorentz contribution to the DF can be decomposed to display the relevant absorption peaks of each oscillator and the presence of collective oscillations of bound electrons. The starting point is the following equation

\[ \frac{1}{\varepsilon_{_L}(\Omega)}=\frac{A_{j}^{(B)}(\Omega)}{\varepsilon_{_{L \ , \ j}}(\Omega)}+\frac{B_{j}^{(B)}(\Omega)}{\sum_{i=j+1}^{K} \varepsilon_{_{L \ , \ i}}(\Omega)}, \tag{6} \]

with j = 1,2, …, K-1, where the real coefficients $A_{j}^{(B)}$ and $B_{j}^{(B)}$ are obtained by solving Equation (6) for their real and imaginary components. From them, the ELFs of each one of the first j-1 oscillators are given by

\[ L_{_{L \ , \ j}}^{(B)}(\Omega)=-\left[\prod_{i=0}^{j-1} B_{i}^{(B)}(\Omega)\right] A_{j}^{(B)}(\Omega) \cdot \operatorname{Im}\left[\frac{1}{\varepsilon_{_{L \ , \ j}}(\Omega)}\right], \tag{7a} \]

and for the last Kth oscillator

\[ L_{L, K}^{(B)}(\Omega)=-\left[\prod_{i=0}^{K-2} B_{i}^{(B)}(\Omega)\right] B_{K-1}^{(B)}(\Omega) \cdot \operatorname{Im}\left[\frac{1}{\varepsilon_{_{L \ , \ j}}(\Omega)}\right], \tag{7b} \]

In this way, $L_{L}^{(B)}(\Omega)=\sum_{j=1}^{K} L_{L, j}^{(B)}(\Omega)$. The effective number of electrons per atom participating in each interband transition can be evaluated from the decomposition of the Lorentz ELF:

\[ Z_{e f f}^{(j)}=\frac{2}{\pi \Omega_{p e}^{2}} \int\limits_{0}^{\Omega_{\max }} \Omega \cdot L_{L, j}^{(B)}(\Omega) d \Omega . \tag{8} \]

Other parameter that one can evaluate from the ELF is the dimensionless square of the dipole-matrix-element (DME) for all possible inelastic scattering processes within the spectral range considered [85,87,88]:

\[ M_{t o t}^{2}=\frac{2 R_{\infty}}{\pi \Omega_{p e}^{2}} \int\limits_{0}^{\Omega_{\max }} \operatorname{Im}\left[-\frac{1}{\varepsilon(\Omega)}\right] d \Omega=\frac{2 R_{\infty}}{\pi \Omega_{p e}^{2}} \int\limits_{0}^{\Omega_{\max }} L^{(B)}(\Omega) d \Omega, \tag{9} \]

where R is the Rydberg constant (R = 13.6057 eV). This DME parameter has been defined as a dimensionless one by dividing the summation of individual DMEs between eaoEo where ao is Bohr’s radius and Eo is the modulus of the amplitude of the incident field (see Equation 2.11 in [88]). From the decompositions performed, the square of DMEs corresponding to the charge-carriers, bound electrons, and individual oscillators can also be evaluated:

\[ M_{D}^{2}=\frac{2 R_{\infty}}{\pi \Omega_{p e}^{2}} \int\limits_{0}^{\Omega_{\max }} L_{D}^{(B)}(\Omega) d \Omega, \tag{10a} \]

\[ M_{L}^{2}=\frac{2 R_{\infty}}{\pi \Omega_{p e}^{2}} \int\limits_{0}^{\Omega_{\max }} L_{L}^{(B)}(\Omega) d \Omega, \tag{10b} \]

\[ M_{L,j}^{2}=\frac{2 R_{\infty}}{\pi \Omega_{p e}^{2}} \int\limits_{0}^{\Omega_{\max }} L_{L,j}^{(B)}(\Omega) d \Omega, \tag{10c} \]

The evaluations reported in Table 5 show that $Z_{e f f}=Z_{e f f}^{(c c)}+Z_{e f f}^{(i t)}$, $Z_{e f f}^{(c c)}$≅ ze + zh = 0.586, $M_{tot}^{2}=M_{D}^{2}+M_{L}^{2}$, and

\[ Z_{e f f}^{(i t)}=\sum_{j=1}^{K} Z_{ {eff }}^{(j)}, \quad M_{L}^{2}=\sum_{j=1}^{K} M_{L, j}^{2}. \tag{10d} \]

A comparison of the $Z_{e f f}^{(j)}$ values with the corresponding fj and γj parameters show that a larger absorption width and the oscillator strength, which is proportional to the square of the dipole-matrix-element associated with the specific probability of transition, is associated with a larger effective number of electrons per atom participating in the jth interband transitions.

Table 5 Set of derived quantities corresponding to effective numbers of charge-carriers (bound electrons) per atom doing intraband (interband) transitions in the spectral range considered, 0.2 to 40 eV, as well as the square of the dimensionless dipole-matrix-elements

8. Densities of States at the Final States Involved in Interband Transitions

For a given jth transition between bands, the absorption rate of energy can be obtained from Fermi’s golden rule times the absorbed energy Ωj [89]. The golden rule states that the transition rate is proportional to the density of states of the final states involved in the transition and the square matrix element of the Hamiltonian associated with the perturbation of the system due to the interaction between the propagating electric field and the induced dipoles. We have approached this square Hamiltonian matrix element in terms of (eaoEo)2Mj2. In classical electrodynamics, the absorption rate of energy corresponding to the jth oscillator, at the resonance frequency ωj, and within the framework devised to decompose the Lorentz ELF in terms of those corresponding to each oscillator, is proportional to the corresponding jth ELF [81]. By considering the formulation of these two approaches for the rate of absorption as equal, the following equation is obtained

\[ \frac{(2 \pi)^{2}}{\mathrm{~h}}\left[M_{j}^{2}\left(e a_{o} E_{o}\right)^{2}\right] \rho\left(E_{j}\right) \Delta E_{j}=\varepsilon_{o} E_{o}^{2}\left(\frac{2 \pi \Omega_{j}}{h}\right) \operatorname{Im}\left[-\frac{C_{j}^{(B)}\left(\Omega_{j}\right)}{\varepsilon_{L, j}\left(\Omega_{j}\right)}\right], \tag{11} \]

with $\Delta E_{j}=E_{j}-E_{b, j}$ where $E_{j}$ is the energy of the final state, and $E_{b, j}$ is the energy in the initial state of the transition, with $\Delta E_{j} \simeq \Omega_{j}$. Moreover

\[ C_{j}^{(B)}\left(\Omega_{j}\right)=\left[\prod_{i=0}^{j-1} B_{i}^{(B)}\left(\Omega_{j}\right)\right] A_{j}^{(B)}\left(\Omega_{j}\right), \tag{12} \]

as indicated above. In this way, we obtained the following expression for the density of states at the final states involved in the transitions:

\[ \rho\left(\Omega_{j}\right)=\frac{\varepsilon_{o}}{2 \pi\left(e a_{o}\right)^{2} M_{j}^{2}} \operatorname{Im}\left[-\frac{C_{j}^{(B)}\left(\Omega_{j}\right)}{\varepsilon_{_{L \ , \ j}}\left(\Omega_{j}\right)}\right]\left(2 V_{p c}\right), \tag{13} \]

where $\rho\left(E_{j}\right) \simeq \rho\left(\Omega_{j}\right)$. The volume of the primitive unit cell (Vpc) was introduced in Equation (13) to obtain the density of states in units of states/energy atom, and a factor of 2 was added to account for the two accessible spin states of each electron. Table 6 shows the evaluations of the densities of states.

Table 6 Densities of states (in states/eV atom) at final states associated with each effective set of interband transitions

We finished our analysis with a comparison of the estimated ρ(Ωj)-values with those obtained from the DFT calculations reported in Section 3. Figure 6 and Figure 7 show the BSD displayed through an energy range larger than that used in Figure 1 and Figure 2. The right sides of these figures show the DOS using a logarithmic scale, which is appropriate to observe the orders of magnitudes of the densities of states when final states correspond to s, p, d, or f states, first for DTF calculations with no spin-orbit coupling and subsequently accounting for it, respectively. Through what follows, we compared the values obtained for ρ(Ωj), as reported in Table 6, with corresponding values ρ(Ej) estimated from the BSD when the energy of the final states (Ej) was projected toward the DOS in Figure 6 or Figure 7. We showed to aim that the ρ(Ωj)-values are consistent in order of magnitude with those of ρ(Ej) when considering plausible transitions displayed in the BSD. Although our analysis is not intended to be exhaustive, it provides a significant degree of confidence in the method devised to obtain the ρ(Ωj)-values. We first considered the low-energy infrared transitions. According to the BSD displayed in Figure 2(a), the first resonance energy Ω1 = 0.42 eV could arise from transitions between parallel bands displayed when approaching the L-symmetry point from the left, in the Λ-L direction (see Figure 8(a)). The light absorption due to these transitions was weak according to the low value of f1γ1. From the APCSA optimization, the fjγj values were between 0.17 and 66.3, with f1γ1 = 0.38. Electrons with energies between around 0.55 eV below and 0.53 eV above the SFL could participate in these transitions, showing that the largest density of states for the transitions with final states above the Fermi level (see the small peak of the s-DOS green line just above the SFL in Figure 6), i.e., ΔE1 was about 0.54 eV. This largest peak value was close to ρ(E1) = 0.9 × 10-3 states/eV atom, which is of the same order of magnitude as ρ(Ω1) = 0.001 states/eV atom if one considered that the electrons are transitioning into p-states. The bound electrons involved in these transitions could be sp- or d-electrons, according to Figure 2(b) and Figure 2(c), depending on two factors: the number of d-electrons participating in these transitions is significantly larger than that of sp-electrons; however, d-states are more localized than sp and p states.

Click to view original image

Figure 6 The BSD of Ir (left figure) and its total DOS (black line in right figure) and projected densities of states for s (violet line, j = 1/2), p (green and pink lines, j = 3/2 and j = 1/2, respectively), d (clear brown and yellowish lines, j = 5/2 and j = 3/2, respectively), and s with no spin-orbit coupling (blue line) states, in units of states/eV atom. The vertical scale on the left side is the energy in eV. The 0 energy (red dashed line) indicates the SFL. Except for the blue line of the DOS, the calculations included spin-orbit coupling. The numbered red arrows are used throughout the discussion in the text to indicate plausible interband transitions.

Click to view original image

Figure 7 The BSD of Ir (left figure) and its total DOS (black line in right figure) and projected densities of states for s (green line), p (violet), d (yellowish), and f (red) states, in units of states/eV atom. The vertical scale on the left side is the energy in eV. The 0 energy (red dashed line) indicates the SFL. The spin-orbit coupling has not been accounted for. The red arrow with label 8 is used through the discussion in the text to indicate a plausible interband transition.

Click to view original image

Figure 8 Amplified views of the BSD to show in detail the bands involved in the plausible transitions of bound electrons. The vertical scales are given in eV.

Regarding the second set of inter-band transitions, whose energy difference was close to Ω2 = 1.21 eV with f2γ2 = 0.73, these could be located when going from the Δ- to the X-symmetry point, close to the middle point where two almost superposed bands cross the SFL having a band just below them and another just above (see Figure 8(b)). At the point where the two superimposed bands crosse the SFL Eb,2 = -1.35 eV, E2 = 0 eV, and ρ(E2) = 0.028 states/eV atom if we assumed that these excitations involved sp-electrons transitioning to s-states with no spin-orbit coupling (see the blue line in Figure 7), according to the low value of f2γ2 = 0.73. At the left side of the point where the two superposed bands cross the SFL, E2 = 0.78 eV, Eb,2 = 0 eV, and ρ(E2) = 0.055 states/eV atom. These two estimations of ρ(E2) have the same order of magnitude as the value obtained from the APCSA approach, ρ(Ω2) = 0.010 states/eV atom.

The density of states of the third set of transitions was evaluated in ρ(Ω3) = 0.032 states/eV atom, with f3γ3 = 0.17 and Ω3 = 2.07 eV. These could be attributed to transitions in the K-symmetry point, also involving the two parallel bands displayed from K to Σ (see Figure 8(c)). The energy below the SFL, with available empty states just above, can be estimated from Eb,3 = -1.60 eV, and the energy just above the SFL in the K point was close to E3 = 0.48 eV above the SFL, i.e., ΔE3 = 2.08 eV, and the value of ρ(E3) was between 0.028 and 0.041 states/eV atom, also assuming that sp-electrons were excited to s-states with no spin-orbit coupling, as observed, ρ(E3) $\simeq $ ρ(Ω3). Other transitions that could be contributing to this set are displayed in the L-symmetry point, between the bands 0.92 eV below and 1.04 above the SFL (see Figure 8(e)). If the transitions of these sp-electrons are into p- or s-states with no spin-orbit coupling, the density of states would be the order of 0.001 or about 0.095 states/eV atom, respectively.

For the transition sets associated with Ω4 = 3.22 eV, the density of final states was calculated as ρ(Ω4) = 0.005 states/eV atom. According to Figure 8(d), this set of transitions could occur in the W-point between the band closer to the SFL, at Eb,4 = 1.74 eV below, and the upper band at E4 = 1.30 eV above (ΔE4 = 3.04 eV). The projection of the E4-energy toward the DOS spectrum indicates that the density of states would be between 0.0004 and 0.032 states/eV atom, depending on the characteristic of the final states: p-states with j = 1/2 and j = 3/2 or s-states with no coupling between the spin and angular momentum of the electrons, respectively. The intermediate value of f4γ4 = 8.38 and the bands displayed in Figure 2(b) and (c) at the W-symmetry point suggested that both sp- and d-electrons participated in these transitions.

The fifth set of transitions corresponded to an energy difference of about Ω5 = 4.38 eV, with a density of states ρ(Ω5) = 0.046 states/eV atom for the final states. Transitions in the L-point between bands separated by the energy difference ΔE5 = E5-Eb,5 = 4.48 eV, with E5 = 1.04 eV above the SFL and Eb,5 = 3.44 eV below it (see Figure 8(e)). The projection of the E5-value toward the DOS figure showed that the density of states could be around ρ(E5) = 0.088 states/eV atom involving the transitions of sp-electrons to s-states with no spin-orbit coupling, which would be consistent with the low value of f5γ5 = 0.24. Other transitions that contribute to this set could be occurring up to the band 0.46 eV above the SFL. In that case, ρ(E5) = 0.041 states/eV atom.

We finally considered the transitions at the vacuum ultraviolet, which are characterized by large values of fjγj. The effective sixth set of transitions was characterized by a resonance energy Ω6 = 11.10 eV with a density of states ρ(Ω6) = 0.015 states/eV atom. According to the BSD depicted in Figure 8(f), plausible transitions at the L-point whose energy difference was ΔE6 = E6 - Eb,6 = 11.12 eV, with E6 = 6.82 eV above the Fermi level and Eb,6 = 4.30 eV below, and a corresponding density of states close to ρ(E6) = 0.013 states/eV which corresponded to d-electrons transitioning to d-states (with j = 3/2). This is consistent with the BSDs of Figure 2(b) and Figure 2(c) and with a large value of f6γ6 = 10.08.

For the seventh set of transitions, Ω7 = 19.43 eV with ρ(Ω7) = 0.006 states/eV atom. The BSD showed plausible transitions at the L-direction from the lower band at Eb,7 = 7.73 eV below the Fermi level up to the upper band at E7 = 13.24 eV above, i.e., ΔE7 = 20.97 eV (see Figure 8(g)). At this E7-spectral position, the DOS in Figure 6 indicated that final d-states (with j = 5/2) were involved in these transitions with a density of states close to ρ(E7) = 0.004 states/eV atom and with initial d-states.

For the eighth set of transitions, we used the information displayed in Figure 7. The resonance frequency was estimated as Ω8 = 31.17 eV with the density of states ρ(Ω8) = 0.015 states/eV atom. According to Figure 8(h), the BSD showed plausible transitions at the Γ-point from the lower band at Eb,8 = 10.71 eV below the Fermi level up to the upper band at E8 = 19.11 eV above the Fermi level, i.e., ΔE8 = 29.82 eV. At this E8-spectral position, the DOS figure indicated that final f-states were involved in these transitions of d-electrons, with a density of states close to ρ(E8) = 0.070 states/eV.

When the transitions involved only p-electrons, the density of states approached the 0.1 states/eV atom value. Close to this Γ-point, in the Δ-direction, there was a section where one could see parallel bands contributing to this effective set of transitions, with decreasing densities of states up to about 0.010 states/eV atom. According to the value obtained for ρ(Ω8), these transitions contributed more significantly and involved s- and d-electrons.

We finished the analysis correlating the spectral composition of the DF’s imaginary component with the set of interband transitions displayed in previous figures. Figure 9(a) shows the spectrum of ε2 with an indication of spectral positions of the resonance energies and the symmetry points or parallel bands involved in the quantum transitions of bound electrons. Figure 9(b) displays the decomposition of the Lorentz contribution to the DF as provided by each set of transitions.

Click to view original image

Figure 9 (a) Imaginary component of Ir’s DF indicating that the spectral positions of the resonances energies correlated with the transitions indicated in the inset: RP (LP) means right (left) side of the Λ (Δ) point involving parallel (quasi-parallel) bands, and SP means symmetry point. The bands involved in the transitions have been indicated in the two previous figures. Dots are the original data taken from [20], and the solid line displays the evaluation from the set of optimized parameters specified in Table 2. (b) The contribution of each set of interband transitions to the DF of Ir.

9. Summary and Conclusions

A framework to model the dielectric function of transition metals has been devised, with Ir as the material considered in this work. The model is an extended version of the Drude–Lorentz model, which incorporates the contribution of holes as charge-carriers. The complete set of model parameters was optimized using a simulated annealing approach assisted by a projected gradient method. The physical meaning of the model parameters was emphasized, going beyond a simple fitting procedure. Several derived quantities were obtained once the optimization procedure was finished, and additional optical, charge transport, magnetic, and quantum mechanical properties of the material were reported. Background knowledge was provided by Density Functional Theory calculations, including spin-orbit coupling. Within this background, a consistent interpretation of Fermi energies and density of states at Fermi energies of holes and electrons was provided. In addition to the parametric spectral description of the dielectric function, other optical parameters are obtained from the optimization process or as derived quantities: plasma frequency of s-conduction electrons, effective plasma frequency of both charge carriers (holes and electrons), intrinsic relaxation times of holes and electrons, high-frequency dielectric constant, the number density of electrons and holes, and their effective masses. The optical response of bound electrons was required to optimize the values of oscillator strengths, resonance energies, and absorption widths. Furthermore, several charge transport properties were obtained as derived quantities both for electrons and holes: Fermi velocities, intrinsic mean free paths, mobilities, static conductivity, dynamic conductivity contributions, Hall’s coefficient, intrinsic resistivity, and heat capacities. The magnetic properties evaluated included Pauli’s paramagnetic and Landau’s diamagnetic susceptibilities for electrons and holes. In addition, a large set of quantum mechanical parameters were evaluated: charge-carriers Fermi energies and corresponding densities of states at the Fermi energy. A first decomposition of the energy loss function allowed us to calculate the effective numbers of itinerant charge carriers participating in inter- or intra-band transitions through the spectral range considered, as well as the total effective numbers of bound electrons doing interband transitions. From a second decomposition, focused on the Lorentz contribution to the ELF, a model was devised to obtain the energy loss functions associated with each set of transitions. From these oscillator ELFs, the characteristic of the participating electrons (s, p, d, or f) was established, and the density of states at the final states involved in these transitions was calculated. Comparisons of these densities of states with the data provided by the DFT calculations showed good agreement in terms of their orders of magnitude.


The authors thank the support given by the Universidad de Costa Rica to carry out this research. They also thank Professor Gunnar A. Niklasson, at Uppsala University, Sweden, for reading the manuscript and for his valuable comments.

Author Contributions

Implementation of the SPGM + APCSA programs (WE Vargas), application of the SPGM + APCSA method for modeling the DF of Ir (WE Vargas, M Hernández-Jiménez, and V Quirós-Cordero), DFT calculations (F Muñoz-Rojas and E Avendano), and manuscript preparation (WE Vargas).

Competing Interests

The authors have declared that no competing interests exist.

Appendix: Optimization by a Spectral Projected Gradient Method

The SPGM requires initial input values for the parameters being optimized. They are obtained from previous applications at each temperature cycle of the APCSA method to minimize the merit function F specified by Equation (2), which can be written as $F=F_{o} \sum_{i=1}^{\kappa}\left(F_{1, i}^{2}+F_{2, i}^{2}\right)$ with

\[ F_{1, i}\left(\Omega_{i}\right)=\left[\frac{\varepsilon_{_1}\left(\Omega_{i}\right)}{\bar{\varepsilon}_{1}\left(\Omega_{i}\right)}\right]-1 \text { and } F_{2, i}\left(\Omega_{i}\right)=\frac{\varepsilon_{2}\left(\Omega_{i}\right)}{\bar{\varepsilon}_{2}\left(\Omega_{i}\right)}-1 \]

where the constant factor is given by $F_{o}=\left(2 \kappa-N_{p}-1\right)^{-1}$. The explicit forms of ε1 and ε2 are obtained from Equation (1), namely,

\[ \varepsilon_{_1}(\Omega)=\frac{1}{f_{o}}-\frac{\Omega_{p e}^{2}}{\Omega^{2}+\gamma_{e}^{2}}-\frac{\chi \beta \Omega_{p e}^{2}}{\Omega^{2}+\gamma_{h}^{2}}+\varepsilon_{1, L}, \quad \varepsilon_{2}(\Omega)=\frac{\gamma_{e} \Omega_{p e}^{2}}{\Omega\left(\Omega^{2}+\gamma_{e}^{2}\right)}+\frac{\chi \beta \gamma_{h} \Omega_{p e}^{2}}{\Omega\left(\Omega^{2}+\gamma_{h}^{2}\right)}+\varepsilon_{2, L} \tag{A1} \]


\[ \varepsilon_{1, L}=\frac{1}{z_{e}} \sum_{j=1}^{K} \frac{f_{j} \Omega_{p e}^{2}\left(\Omega_{j}^{2}-\Omega^{2}\right)}{\left(\Omega_{j}^{2}-\Omega^{2}\right)^{2}+\left(\Omega \gamma_{j}\right)^{2}}, \quad \varepsilon_{2, L}=\frac{1}{z_{e}} \sum_{j=1}^{K} \frac{f_{j} \Omega_{p e}^{2} \Omega \gamma_{j}}{\left(\Omega_{j}^{2}-\Omega^{2}\right)^{2}+\left(\Omega \gamma_{j}\right)^{2}}. \tag{A2} \]

The chain rule was applied to evaluate the components of the gradient, which requires the previous evaluation of $d F_{1, i} / d \varepsilon_{1} \text { and } d F_{2, i} / d \varepsilon_{2}$. If p is one of the Np parameters under optimization,

\[ \frac{d F}{d p}=2 F_{o} \sum_{i=1}^{\kappa} F_{1, i}\left(\Omega_{i}\right)\left(\frac{d F_{1, i}}{d \varepsilon_{_1}}\right)\left(\frac{d \varepsilon_{_1}}{d p}\right)_{\Omega_{i}}+2 F_{o} \sum_{i=1}^{\kappa} F_{2, i}\left(\Omega_{i}\right)\left(\frac{d F_{2, i}}{d \varepsilon_{2}}\right)\left(\frac{d \varepsilon_{_2}}{d p}\right)_{\Omega_{i}}, \tag{A3} \]

with $d F_{1, i} / d \varepsilon_{_1}=1 / \bar{\varepsilon}_{_1}\left(\Omega_{i}\right) \text { and } d F_{2, i} / d \varepsilon_{_2}=1 / \bar{\varepsilon}_{_2}\left(\Omega_{i}\right)$. The non-null derivates of ε1 and ε2 respect to each optimization parameter are specified in Table A1.

Table A1 Derivatives with respect to parameters related to the Drude contribution to the dielectric function, and with respect to ze parameter involved in the Lorentz term.

For the optimization parameters involved in the Lorentz term, we used the following expressions for the corresponding derivates related to the jth oscillator’s contribution:

Table A2 Derivatives with respect to parameters related to the Lorentz contribution to the dielectric function


  1. Hunt LB. A history of iridium. Platinum Metals Rev. 1987; 31: 32-41.
  2. Shirvaliloo S, Kangarloo H. Production of iridium metal thin films for application as electrodes in DRAMs and FRAMs. J Electr Electron Eng. 2015; 3: 35-38. [CrossRef]
  3. Youm M, Sim HS, Jeon H, Kim S, Kim YT. Metal oxide semiconductor field effect transistor characteristics with iridium gate electrode on atomic layer deposited ZrO2 high-k dielectrics. Jpn J Appl Phys. 2003; 42: 5010-5013. [CrossRef]
  4. Wang YF, Chang X, Li S, Zhao D, Shao G, Zhu T, et al. Ohmic contact between iridium film and hydrogen-terminated single crystal diamond. Sci Rep. 2017; 7: 12157. [CrossRef]
  5. Hemphill R, Hurwitz M, Pelizzo MG. Osmium atomic-oxygen protection by an iridium overcoat for increased extreme-ultraviolet grating efficiency. Appl Opt. 2003; 42: 5149-5157. [CrossRef]
  6. Ozaki Y. Temperature and humidity dependence of SnO2-based CO gas sensors modified with iridium and ruthenium. Electrochem Solid-State Lett. 1999; 3: 297-299. [CrossRef]
  7. Blaser HU. Application of iridium catalysts in the fine chemicals industry. In: Iridium complexes in organic synthesis. Weinheim: Wiley‐VCH Verlag GmbH & Co. KGaA; 2008. pp. 1-14. [CrossRef]
  8. Oro A, Claver C. Iridium catalysts for organic reactions. Cham: Springer Nature Switzerland AG; 2021. [CrossRef]
  9. John S, Wang R. Metallic photonic-band-gap filament architectures for optimized incandescent lighting. Phys Rev A. 2008; 78: 043809. [CrossRef]
  10. Wen RT, Niklasson GA, Granqvist CG. Strongly improved electrochemical cycling durability by adding iridium to electrochromic nickel oxide films. ACS Appl Mater Interfaces. 2015; 7: 9319-9322. [CrossRef]
  11. Sachse R, Pflüger M, Velasco-Vélez JJ, Sahre M, Radnik J, Bernicke M, et al. Assessing optical and electrical properties of highly active IrOx catalysts for the electrochemical oxygen evolution reaction via spectroscopic ellipsometry. ACS Catal. 2020; 10: 14210-14223. [CrossRef]
  12. El Sayed AM, Shaban M. Morphological, surface and optical properties of spin-coated IrOx films; influence of spin speed, annealing and (Cr, La) codoping. Ceram Int. 2019; 45: 8460-8470. [CrossRef]
  13. Coblentz WW. The reflecting power of various metals. Bull Bur Stand. 1910; 7: 197-225. [CrossRef]
  14. Barnes BT. Optical constants of incandescent refractory metals. J Opt Soc Am. 1966; 56: 1546-1550. [CrossRef]
  15. Hass G, Jacobus GF, Hunter WR. Optical properties of evaporated iridium in the vacuum ultraviolet from 500 Å to 2000 Å. J Opt Soc Am. 1967; 57: 758-762. [CrossRef]
  16. Samson JAR, Padur JP, Sharma A. Reflectance and relative transmittance of laser-deposited iridium in the vacuum ultraviolet. J Opt Soc Am. 1967; 57: 966-967. [CrossRef]
  17. Kirillova M, Nomerovannaya LV, Noskov MM. Optical properties of iridium at 77 and 295 K. Fiz Metal Metalloved. 1972; 34: 291-296.
  18. Weaver JH. Optical properties of Rh, Pd, Ir, and Pt. Phys Rev B. 1975; 11: 1416-1425. [CrossRef]
  19. Weaver JH, Olson CG, Lynch DW. Optical investigation of the electronic structure of bulk Rh and Ir. Phys Rev B. 1977; 15: 4115-4118. [CrossRef]
  20. Weaver JH. Optical properties of metals. In: Handbook of chemistry and physics. 73rd ed. Boca Raton: CRC Press; 1993. pp. 111-126.
  21. Windt DL, Cash WC, Scott M, Arendt P, Newnam B, Fisher RF, et al. Optical constants of thin films of Ti, Zr, Nb, Mo, Ru, Rh, Pd, Ag, Hf, Ta, W, Re, Ir, Os, Pt, and Au from 24 Å to 1216 Å. Appl Opt. 1988; 27: 246-278. [CrossRef]
  22. Yan L, Woollam JA. Optical constants and roughness study of dc magnetron sputtered iridium films. J Appl Phys. 2000; 92: 4386-4392. [CrossRef]
  23. Lehmuskero A, Kuittinen M, Vahimaa P. Refractive index and extinction coefficient dependence of thin Al and Ir films on deposition technique and thickness. Opt Express. 2007; 15: 10744-10752. [CrossRef]
  24. Yang LH, Menyhárd M, Sulyok A, Tokési K, Ding ZJ. Optical properties and excitation energies of iridium derived from reflection electron energy loss spectroscopy spectra. Appl Surf Sci. 2018; 456: 999-1003. [CrossRef]
  25. Krog-Andersen O. Fermi surfaces and effective masses in fcc transitions metals. Solid State Commun. 1968; 6: 285-290. [CrossRef]
  26. Krogh-Andersen O. Electronic structure of the fcc transition metals Ir, Rh, Pt, and Pd. Phys Rev B. 1970; 2: 883-906. [CrossRef]
  27. Arbman GO, Hornfelt S. Crystal potentials and the band electrons of iridium. J Phys F Metal Phys. 1972; 2: 1033-1045. [CrossRef]
  28. Noffke J, Fritsche L. Band structure calculation and photoemission analysis of iridium. J Phys F Metal Phys. 1982; 12: 921-933. [CrossRef]
  29. Choi WS, Seo SSA, Kim KW, Noh TW, Kim MY, Shin S. Dielectric constants of Ir, Ru, Pt, and IrO2: Contributions from bound charges. Phys Rev B. 2006; 74: 205117. [CrossRef]
  30. Ahmed S, Zafar M, Shakil M, and Choudhary MA. Ab initio study of structural, electronic, and thermal properties of Ir1-xRhx alloys. Condens Matter Phys. 2015; 18: 23601. [CrossRef]
  31. Sigalas MM, Papaconstantopoulos DA. Calculations of the total energy, electron-phonon interaction, and Stoner parameter for metals. Phys Rev B. 1993; 50: 7255-7261. [CrossRef]
  32. Abarca-Quesada DA, Vargas-Quirós CE, Vargas WE. Physical properties of rhodium retrieved from modeling its dielectric function by a simulated annealing approach. OSA Continuum. 2021; 4: 3233-3253. [CrossRef]
  33. Vargas WE. Optical, magnetic, and charge-carriers transport properties of a transition metal: Bulk palladium. Appl Opt. 2017; 56: 6496-6503. [CrossRef]
  34. Vargas WE. Dielectric functions of Pd and Zr transition metals: An application of Drude-Lorentz models with simulated annealing optimization. Appl Opt. 2017; 56: 1266-1276. [CrossRef]
  35. Arblaster JW. Crystallographic properties of iridium. Platinum Metals Rev. 2010; 54: 93-102. [CrossRef]
  36. Volkenshtein NV, Novoselov VA, Startsev VE. The Fermi surface of iridium. Sov Phys JETP. 1970; 31: 862-863.
  37. Theye ML. Investigation of the optical properties of Au by means of thin semitransparent films. Phys Rev B. 1970; 2: 3060-3078. [CrossRef]
  38. Nagel SR, Schnatterly SE. Frequency dependence of the Drude relaxation time in metal films. Phys Rev B. 1974; 9: 1299-1303. [CrossRef]
  39. Eisenhammer T. Photon frequency dependent electron relaxation time in noble metals: Effect of voids. Thin Solid Films. 1995; 270: 55-59. [CrossRef]
  40. Rakic AD, Elazar JM, Djurisic AB. Acceptance-probability-controlled simulated annealing: A method for modeling the optical constants of solids. Phys Rev E. 1995; 52: 6862-6867. [CrossRef]
  41. Rakic AD, Djurisic AB, Elazar JM, Majewski ML. Optical properties of metallic films for vertical-cavity optoelectronic devices. Appl Opt. 1998; 37: 5271-5283. [CrossRef]
  42. Kirkpatrick S, Gelatt CD, Vecchi MP. Optimization by simulated annealing. Science. 1983; 220: 671-680. [CrossRef]
  43. Raydan M. The Barzilai and Borwein gradient method for the large-scale unconstrained minimization problem. SIAM J Optim. 1997; 7: 26-33. [CrossRef]
  44. Ramírez-Porras A, Vargas WE. Transmission of visible light through oxidized copper films: Feasibility of using a spectral projected gradient method. Appl Opt. 2004; 43: 1508-1514. [CrossRef]
  45. Morrish AH. The physical principles of magnetism. New York: IEEE Press; 2001. [CrossRef]
  46. Ziman JM. Electrons and phonons. London: Oxford University Press; 1960.
  47. Batyev EG. Pauli paramagnetism and Landau diamagnetism. Phys Usp. 2009; 52: 1245-1246. [CrossRef]
  48. Hamman DR, Schluter M, Chiang C. Norm-conserving pseudopotentials. Phys Rev Lett. 1979; 43: 1494-1497. [CrossRef]
  49. Kohn W, Sham LJ. Self-consistent equations including exchange and correlation effects. Phys Rev. 1964; 140: A1133-A1138. [CrossRef]
  50. Giannozzi P, Baroni S, Bonini N, Calandra M, Car R, Cavazzoni C, et al. QUANTUM ESPRESSO: A modular and open-source software project for quantum simulations of materials. J Phys Condens Matter. 2009; 21: 395502. [CrossRef]
  51. Perdew JP, Burke K, Ernzerhof M. Generalized gradient approximation made simple. Phys Rev Lett. 1996; 77: 3865-3868. [CrossRef]
  52. Zhang Y, Yang W. Comment on Generalized gradient approximation made simple. Phys Rev Lett. 1998; 80: 890. [CrossRef]
  53. Troullier N, Martins JL. Efficient pseudopotential for plane-wave calculations. Phys Rev B. 1991; 43: 1993-2006. [CrossRef]
  54. Singh HP. Determination of thermal expansion of germanium, rhodium, and iridium by X-rays. Acta Crystallogr Sect A. 1968; 24: 469-471. [CrossRef]
  55. Fischer M, Evers FO, Formalik F, Olejniczak A. Benchmarking DFT-GGA calculations for the structure optimisation of neutral-framework zeotypes. Theor Chem Acc. 2016; 135: 257. [CrossRef]
  56. Methfessel M, Paxton AT. High-precision sampling for Brillouin-zone integration in metals. Phys Rev B. 1989; 40: 3616-3621. [CrossRef]
  57. Monkhorst HJ, Pack JD. Special points for Brillouin-zone integrations. Phys Rev B. 1976; 13: 5188-5192. [CrossRef]
  58. Vargas WE, Clark N, Muñoz-Rojas F, Azofeifa DE, Niklasson GA. Optical, charge transport, and magnetic properties of Pd retrieved from photometric measurements: Approaching the quantum mechanics background. Phys Scr. 2019; 94: 055101. [CrossRef]
  59. Vuillemin JJ, Priestley MG. De Haas-van Alphen effect and Fermi surface in palladium. Phys Rev Lett. 1965; 14: 307-309. [CrossRef]
  60. Reale C. Determination of charge-transport parameters for the group IV metals. Re Bras Fís. 1973; 3: 431-439.
  61. Meaden GT. Electrical resistance of metals. New York: Plenum Press; 1965. [CrossRef]
  62. Hurd CM. The hall effect in metals and alloys. New York: Plenum Press; 1972. [CrossRef]
  63. Powell RW, Tye RP, Woodman MJ. Thermal conductivities and electrical resistivities of the platinum metals. Platinum Met Rev. 1962; 6: 138-143.
  64. Aisaka T, Shimizu M. Electrical resistance, thermal conductivity and thermoelectric power of transition metals at high temperatures. J Phys Soc Japan. 1970; 28: 646-654. [CrossRef]
  65. Hörnfeldt S, Hammarström Å, Carrander K, Björck G. The de Haas-van Alphen effect and the Fermi surface in iridium. J Phys Chem Solids. 1971; 32: 753-760. [CrossRef]
  66. Grodski JJ, Dixon AE. de Hass-van Alphen effect in iridium. Phys Rev B. 1972; 6: 1198-1209. [CrossRef]
  67. Albert HJ, Rubin LR. Magnetic properties of the platinum metals and their alloys. In: Platinum group metals and compounds. Washington: American Chemical Society; 1971. pp. 1-16. [CrossRef]
  68. Vosko SH, Perdew JP. Theory of the spin susceptibility of an inhomogenous electron gas via the density functional formalism. Can J Phys. 1975; 53: 1885-1897. [CrossRef]
  69. Allen PB. The electron-phonon coupling constant. In: Handbook of Superconductivity. New York: Academic Press; 1999. pp. 478-483.
  70. Cahaya AB. Paramagnetic and diamagnetic susceptibility of infinite quantum well. J Mater Sci Geophys Instrum Theor Phys. 2020; 3: 61-67. [CrossRef]
  71. Lindhard J. On the properties of a gas of charge particles. Dan Mat Fys Medd. 1954; 28: 8.
  72. Dorofeev IA, Vinogradov EA. Structure of electromagnetic fields near the surface of an ionic crystal. J Surf Investig. 2012; 6: 796-804. [CrossRef]
  73. Ashcroft NW, Mermin ND. Solid state physics. New York: Saunders College Publishing; 1976.
  74. Furukawa GT, Reilly ML, Gallagher S. Critical analysis of heat-capacity data and evaluation of thermodynamics properties of ruthenium, rhodium, palladium, iridium, and platinum from 0 to 300 K. A survey of the literature data on osmium. J Phys Chem Ref Data. 1974; 3: 163-209. [CrossRef]
  75. McMillan WL. Transition temperature of strong-coupled superconductors. Phys Rev. 1968; 167: 331-344. [CrossRef]
  76. Stewart GR. Measurement of low-temperature specific heat. Rev Sci Instrum. 1983; 54: 1-11. [CrossRef]
  77. Luttinger JM, Ward JC. Ground-state energy of a many-fermion system. II. Phys Rev. 1960; 118: 1417-1427. [CrossRef]
  78. Kas JJ, Blanton TD, Rehr JJ. Exchange-correlation contributions to thermodynamic properties of the homogenous electron gas from a cumulant Green’s function approach. Phys Rev B. 2019; 100: 195144. [CrossRef]
  79. Papaconstantopoulos DA. Handbook of the band structure of elemental solids. From Z = 1 to Z = 112. New York: Springer; 2015. [CrossRef]
  80. Mott NF, Fowler RH. The electrical conductivity of transition metals. Proc R Soc London Ser A. 1936; 153: 699-717. [CrossRef]
  81. Dressel M, Grüner G. Electrodynamics of solids. Cambridge: Cambridge University Press; 2003.
  82. Sun Y, Xu H, Da B, Mao SF, Ding ZJ. Calculation of energy-loss function for 26 materials. Chinese J Chem Phys. 2016; 29: 663-670. [CrossRef]
  83. Madelung O. Introduction to solid-state theory. Berlin: Springer; 1996.
  84. Tanuma S, Powell CJ, Penn DR. Use of sum rules on the energy-loss function for the evaluation of experimental optical data. J Electron Spectros Relat Phenomena. 1993; 62: 95-109. [CrossRef]
  85. Tanuma S, Powell CJ, Penn DR. Inelastic mean free paths of low-energy electrons in solids. Acta Phys Pol A. 1992; 81: 169-186. [CrossRef]
  86. Burden RL, Faires DJ, Burden AM. Numerical analysis. USA: Cengage Learning; 2015.
  87. Fano U. Atomic theory of electromagnetic interactions in dense materials. Phys Rev. 1956; 103: 1202-1218. [CrossRef]
  88. Inokuti M. Inelastic collisions of fast charged particles with atoms and molecules—The Bethe theory revisited. Rev Mod Phys. 1971; 43: 297-347. [CrossRef]
  89. Bransden BH, Joachain CJ. Introduction to quantum mechanics. USA: Prentice Hall; 2000.
Download PDF Download Citation
0 0