Epitaxial VNi2 layers are deposited onto MgO(001) and their resistivity ρ measured as a function of layer thickness d = 10.5–138 nm to quantify the resistivity size effect. The layers exhibit resistivity minima at both stoichiometric V:Ni = 1:2 composition and 700 °C growth temperature, which is attributed to electron scattering at V–Ni antisite defects and is described using the ordering parameter within the framework by Rossiter. A cube-on-cube epitaxy of the fcc parent structure on MgO(001) leads to two possible layer orientations for orthorhombic VNi2(010) and VNi2(103), resulting in considerable atomic disorder at domain boundaries, consistent with relatively small x-ray coherence lengths of 8 and 14 nm in-plane and along the growth direction of a 33.5 nm thick layer. In situ ρ vs d measurements yield a bulk resistivity of ρo = 46 ± 2 μ Ω cm and a benchmark quantity of ρoλ = (138 ± 5) × 10−16 Ω m2, where λ is the bulk electron mean free path. Air exposure causes a minor resistivity increase due to 2 ± 1 nm thick surface oxide that perturbs the surface potential. Resistivities at 77 K are Δρ = 16 ± 3 μ Ω cm below those at room temperature. This Δρ is thickness independent and is close to the previously predicted 13.9 μ Ω cm bulk resistivity for VNi2 along [100]. However, the measured bulk resistivity is well above this prediction, which is attributed to electron scattering at domain boundaries/atomic disorder. Consequently, the theoretically predicted superior directional conduction cannot be experimentally confirmed in this study. The overall results indicate that VNi2 is only a promising compound for narrow interconnects if a synthesis scheme can be developed that results in a strong atomic order, a negligible domain boundary density, and a [100] crystalline orientation along the transport direction.

The electrical resistivity ρ of metallic thin films or wires increases as their thickness or width decreases below the electron mean free path.1–3 This increase is referred to as the resistivity size effect and impedes the continued downscaling of devices in integrated circuits because of the resulting increase in power consumption and resistance–capacitance (RC) delays in the metal wires that interconnect transistors.4,5 Electron scattering at surfaces6–9 and grain boundaries10,11 are the main cause for the resistivity size effect, as described by the semiclassical Fuchs and Sondheimer (FS)12,13 and Mayadas and Shatzkes (MS)14 models, respectively. These models predict that both surface scattering and grain boundary scattering result in a resistivity contribution that is proportional to the product of the bulk resistivity ρo times the bulk electron mean free path λ.15 Consequently, the search for candidate materials that conduct well at small dimensions and are promising to replace the current copper interconnects focuses on metals with a small ρoλ product.16–19 However, despite considerable research efforts, replacing Cu (ρoλ = 6.7 × 10−16 Ω m2)20 with a reliable alternative conductor remains a major challenge.

Binary and ternary compounds have been proposed as a possible alternative conductor to replace elemental metals for highly scaled interconnects based on their promising ρoλ descriptor,21 interdiffusion reliability,22 and back-end-of-line-compatible growth temperature.23 Reports on intermetallic compounds including CuAl2,24 NiAl,23 Cu2Mg,25 and CuTi26 indicate a promising resistivity scaling that is comparable to elemental metals. They also exhibit good time-dependent dielectric breakdown reliability and electromigration properties.27,28 Layered MAX phase materials have even lower (below elemental metals) resistivity scaling for transport in the in-plane direction, which has been experimentally confirmed for Ti4SiC329 and V2AlC30 and is attributed to their predicted anisotropic Fermi surfaces.21,31 However, their bulk resistivities are likely too high to be competitive for interconnect applications.29 Recently, first-principles computational methods have revealed additional promising compounds with electron velocities that are concentrated near a single transport axis, providing the possibility for wire conduction with minimal electron scattering at all four (top/bottom/left/right) surfaces of an interconnect line.21,32 One of these promising compounds is VNi2, which has a calculated cohesive energy that is 1.5 times higher than that of Cu, and a low predicted effective ρoλ = 4.5 × 10−16 Ω m2 along the [001] crystalline direction,21 suggesting that it would outperform both the stability and the conductivity of Cu for narrow interconnect lines. Relatively little is known about the VNi2 compound, as previous studies have mainly focused on its thermodynamic and structural properties as well as its synthesis during high-temperature solidification.33,34 Thus, the experimental confirmation of its superior conduction at small dimensions and even a measurement of its bulk resistivity is still lacking. Therefore, the goal of the present study is to measure the resistivity of VNi2 thin films and validate its low resistivity size effect.

In this article, we report on the resistivity scaling of the ordered intermetallic VNi2 compound. Deposition conditions during co-sputtering are tuned to optimize composition, crystalline quality, phase purity, surface smoothness, and atomic ordering. We use in situ, ex situ, and cryogenic transport measurements and determine ρo = 46 ± 2 μ Ω cm and λ = 30 ± 5 nm at 293 K. The measured room temperature resistivity is three times higher than the previously predicted value along [100],21 which is ascribed to structural multiplicity. More specifically, two possible layer orientations VNi2(103) and VNi2(010) share the same atomic positions but different ordering of V and Ni atoms in a slightly distorted face-centered cubic (fcc) parent structure. X-ray diffraction (XRD) cannot distinguish between the two structurally similar orientations. However, the deviation from long-range atomic order causes extra electron scattering and a correspondingly high ρo. Thus, VNi2 is a promising interconnect conductor only if future studies can achieve a single domain or even single-crystal VNi2 growth and demonstrate the high conductivity along the [100] direction.

VNi2 thin films were deposited onto MgO(001) in a three-chamber ultrahigh vacuum DC magnetron sputtering system with a base pressure of <10−9 Torr.35 Single-side polished 10 × 10 mm2 MgO(001) substrates were cleaned in consecutive ultrasonic baths of trichloroethylene, acetone, iso-propyl alcohol, and de-ionized water,36,37 blown-dried with dry N2, mounted on a Mo sample holder, inserted into the deposition system through a load lock, and degassed in vacuum at 1000 °C for 1 h. Experiments using other substrates including A plane, C plane, and R plane Al2O3 yielded VNi2 layers that exhibit multiple crystalline orientations. That is, they are polycrystalline and, therefore, less suited for the intended transport studies and, in turn, are not reported in this paper. 5 cm-diameter 99.9% V and 99.99% Ni targets facing the substrate at −45 and +45 tilts were sputter-cleaned for 5 min with a shutter covering the substrate prior to depositions in 3 mTorr 99.999% pure Ar, employing constant powers of 50 and 57 W applied to V and Ni sources to achieve a stoichiometric 1:2 deposition flux with a deposition rate of 0.11 nm/s. The substrates were continuously rotated to maximize compositional uniformity and were kept at 700 °C, which was the optimal temperature to facilitate epitaxial layer growth, good crystalline quality, low surface roughness, and negligible secondary phase formation. The deposition time was adjusted to obtain a series of VNi2 films with thickness d = 10.5–138 nm, as measured by x-ray reflectivity (XRR). The thickness of the thickest d = 138 nm layer was determined from the deposition rate measured from thinner samples because no XRR oscillations can be resolved for approximately d > 100 nm. After cooling down to room temperature, the samples were transferred in vacuum to the analysis chamber for in situ electrical transport measurements using a linear four-point probe with a 1 mm inter-probe spacing operated with 1–100 mA. Subsequently, the samples were removed from the vacuum and submerged in liquid nitrogen to perform similar four-point probe measurements at 77 K, followed by ex situ room temperature resistivity measurements after warming up by continuously blowing dry N2 gas to the sample surface to minimize water condensation.

X-ray diffraction (XRD) and XRR measurements were taken in a PANalytical X’pert PRO MPD system with a Cu Kα source using 45 kV accelerating voltage and 40 mA current. Symmetric θ–2θ scans were obtained using a parallel beam geometry with an x-ray mirror and a PIXcel solid-state line detector with 0.165 mm active length acting as a point detector. ω rocking curves were obtained with the same optics. Azimuthal φ-scans were acquired at constant χ and 2θ angles using a point source in combination with an x-ray lens yielding a quasi-parallel beam with an equatorial and axial divergence of 0.3°.

Photoelectron spectroscopy (XPS) analyses were performed using Al Kα radiation (1486.6 eV) in a PHI 5000 VersaProbe system with a hemispherical analyzer and an eight-channel detector. The sample surfaces were sputter-etched using a 3 keV 2 mA Ar+ beam incident at 45° relative to the surface normal over a 1 mm2 area, yielding a 7 nm/min sputter rate. Ion beam sputtering was performed in 15 s steps, corresponding to 1.75 nm etch depth per step, followed by the acquisition of high-resolution spectra of V 2p, Ni 2p, and O 1s peaks using a pass energy of 50 eV and a step size of 0.2 eV. The layer compositions were determined from the measured area under the curves after Shirley background subtraction and using the sensitivity factors within the PHI MultiPak software package.

Initial experiments focus on determining the required power to the deposition sources and the suitable substrate temperature to achieve a stoichiometric composition and to maximize atomic ordering, respectively. Figure 1(a) shows the measured resistivity ρ as a function of V concentration from a series of 33.5 nm thick VNi2 layers deposited at Ts = 700 °C with a varying power of 45–60 W to the Ni deposition source to obtain a range of compositions near the desired V:Ni = 1:2 ratio. The composition in the x axis is obtained from XPS depth profiles. An example is shown in Fig. S1. The depth profiles indicate a constant composition throughout the thickness of the layer and a negligible oxygen impurity concentration below the surface contamination, which is attributed to water adsorption and/or surface oxidation during transfer in air from the deposition to the XPS analysis chamber. The plotted resistivity in Fig. 1(a) is a strong function of the composition, exhibiting pronounced minima of 54.4 and 40.2 μ Ω cm at 293 and 77 K, respectively, at a composition of 34.3 at. % V and 65.7 at. % Ni. This measured composition is close to the expected 33.3:66.7 ratio for stoichiometric VNi2, suggesting a resistivity minimum for the stoichiometric VNi2 compound.

We attribute the composition dependence of the resistivity to the level of atomic ordering within the VNi2 compound. The resistivity contribution from electron alloy scattering is described with Rossiter's formula, which relates the residual resistivity to the order parameter S,38,
ρ ( S ) = ρ r 1 S 2 1 A S 2 .
(1)
Here, ρ(S) is the bulk residual resistivity (at zero temperature), ρr is the resistivity of the completely random compound (S = 0) at T = 0, and A is a parameter between 0 and 1 depending upon the position of the Fermi surface in relation to the superlattice Brillouin-zone boundaries.39 The order parameter,40 
S = r α x A y β = r β x B y α ,
(2)
depends on the fractions rα and rβ of α and β sites, respectively, being occupied by the correct (V and Ni) atoms, xA and xB are the atomic fractions of V and Ni within the layer, and yα = 0.333 and yβ = 0.667 denote the fraction of α and β sites, respectively. We note that rα = xA for a completely random alloy, resulting in S = 0. In contrast, if all Ni atoms occupy α sites, then rα = xA/yα for xA < yα and rα = 1 for xA > yα, yielding a maximum S = 1 only at the stoichiometric xA = yα. We expect some disorder even at a stoichiometric composition and, therefore, define Smax < 1 as the order parameter at xA = yα, while S decreases for deviations above and below the stoichiometric composition. Correspondingly, a linear approximation for rα yields
r α = x A ( 1 S max ) + S max min ( x A y α , 1 ) ,
(3)
and plugging rα back into Eq. (2) provides an expression for the composition dependence of S,
S = S max y β [ min ( x A y α , 1 ) x A ] .
(4)

The only unknown parameter in Eq. (4) is the order parameter Smax at stoichiometric composition. Plugging Eq. (4) into Eq. (1) yields the blue curve in Fig. 1(a) through the 77 K datapoints and values for the three fitting parameters Smax = 0.95, ρr = 106 μ Ω cm, and A = 0.80. The latter value is relatively large in comparison with the common range for A = 0.3–0.5,41 the value for ρr suggests a resistivity of 106 μ Ω cm for a VNi2 alloy with random lattice occupation, and Smax = 0.95 corresponds to rα = 0.97 for xA = yα, suggesting that approximately 3% of V atoms occupy Ni-sites at a stoichiometric composition. We note that these fitting parameters are correlated, resulting in considerable uncertainty in their values as, for example, a smaller A and a larger ρr yield a curve that matches the measured data nearly as well. In addition, the fitting procedure neglects the resistivity contribution from electron–phonon scattering at 77 K, which is comparable in magnitude to the experimental uncertainty and is estimated to be 2–3 μ Ω cm assuming that, as known for Ni,42 electron–phonon scattering in perfectly ordered VNi2 at 77 K is only 15% of the value at 293 K. The red datapoints in Fig. 1(a) from the 293 K measurements are approximately 15 μ Ω cm above the 77 K values for all compositions, indicating additive resistivity contributions from electron scattering at phonons and lattice defects. That is, Matthiessen's rule applies, indicating that the compositional variation has a negligible effect on the electronic structure and/or crystalline phase. This is also consistent with our XRD analyses (see Fig. S2), indicating a single VNi2 phase for all layers. We note that the temperature coefficient of resistivity is positive and nearly identical for all samples within the investigated composition range. This indicates that all layers exhibit metallic conduction with negligible carrier localization within the measured temperature range T ≥ 77 K. Correspondingly, coulomb interactions43 and quantum (interference) effects44 are negligible for our investigation as the characteristic temperature45 for carrier localization46,47 is well below the investigated temperature range. The temperature effect is attributed purely to electron–phonon scattering and the composition/disorder effect is purely attributed to scattering at potential perturbations.

Figure 1(b) shows a corresponding plot of the measured ρ as a function of Ts = 20–850 °C for a constant composition of 34.3 at. % V. The resistivity exhibits a minimum for the layer deposited at 700 °C but increases at both lower and higher growth temperatures, which is attributed to kinetic barriers for atomic ordering and a thermodynamic driving force for disorder, respectively. The lower atomic ordering also results in a slightly reduced crystalline quality, as quantified by XRD rocking curves shown in Fig. S3, indicating a 2%–9% higher full width at half maximum (FWHM) of 1.32° and 1.41° for Ts = 400 and 850 °C, respectively, than 1.29° for the layer with the lowest resistivity grown at Ts = 700 °C. Reducing Ts further to room temperature dramatically degrades the crystallinity so that no VNi2 XRD peak can be detected (Fig. S3), resulting in a high resistivity of 250 μ Ω cm that is primarily attributed to the limited crystallinity. This value is 2.4 times larger than the fitted value of ρr = 106 μ Ω cm for the alloy with complete atomic disorder (S = 0) but good crystallinity. In contrast, the measured 85.5 and 65.5 μ Ω cm for Ts = 400 and 850 °C are well below ρr and can, therefore, be primarily attributed to a reduced S < 1. As mentioned above, the reduction in ordering at low temperature (Ts = 400 °C) is attributed to kinetic barriers. At higher temperatures (Ts = 850 °C), the entropy from atomic disorder represents a thermodynamic driving force that causes the atomic long-range order parameter S to decrease exponentially as the temperature is increased from approximately 50% to 100% of the critical transition temperature TC for complete randomization.40,41,48 We estimate the transition temperature for VNi2 to be 1872 K, based on TC= −4W1/kB, where kB is the Boltzmann constant and W1 is the formation energy of a single A–B bond (with respect to the average of A–A and B–B bonds),40 estimated to be 40 meV based on the reported VNi2 formation energy of 0.242 eV per atom.32 Thus, we expect entropy-driven disorder to develop above 50% of TC, which corresponds to 936 K (663 °C). This is in good agreement with our experimental observation of an increasing ρ with Ts > 700 °C. We note that 850 °C is well below the melting point of both V49 and Ni,50 so evaporation is not considered in our work.

Figure 2 shows representative XRD and XRR results from a stoichiometric 33.5 nm thick VNi2 layer deposited on MgO(001) at Ts = 700 °C. The θ–2θ pattern in Fig. 2(a) exhibits a strong doublet feature at 42.92° and 43.04° from MgO 002 substrate reflections of the Cu Kα1 and Kα2 x-rays. The peak at 2θ =50.59° is ascribed to a combination of the 020 and 103 reflections of VNi2. It is the only layer peak that can be detected over the entire measured 2θ = 10°–80° range, implying the absence of detectable misoriented grains or secondary phases. Thus, the pattern suggests a single out-of-plane crystalline orientation, which, however, can be VNi2 020 and/or 103. More specifically, the reported interplanar spacings range from 0.1774 to 0.1788 nm for 020 planes and from 0.1786 to 0.1818 nm for 103 planes such that d020 and d103 vary by only 0.1%–1.7%, corresponding to an expected 2θ difference of 0.12°–0.90°.21,33,34 This 2θ range has a comparable magnitude as the measured 0.7° full width at half maximum (FWHM) of the peak in the measured θ–2θ pattern. Thus, this peak cannot resolve the difference between VNi2 020 and 103. In fact, the peak width corresponds to an out-of-plane x-ray coherence length51 of 14 nm, which is considerably smaller than the layer thickness of 33.5 nm, suggesting atomic disorder that corresponds to a domain size along the growth direction of ∼14 nm. Similarly, the ω rocking curve of the VNi2 020/103 reflection shown in the inset in Fig. 2(a) has a FWHM of 1.3°. This may be attributed to a 1.3° mosaicity (crystalline tilt variability) but also corresponds to an in-plane coherence length of 8.0 nm, suggesting in-plane atomic disorder consistent with 8.0 nm wide VNi2(020) and VNi2(103) domains. Nevertheless, despite these complications due to the structural multiplicity, which is further discussed below, the XRD peak demonstrates the formation of the ordered VNi2 phase. Second, the peak position provides a measured value for the interplanar spacing of 0.1804 nm. This is 1.3% and 0.1% larger than the reported d020 and d103,21,33 indicating in-plane compressive stress, which is attributed to differential thermal contraction during cooling from the Ts = 700 °C growth temperature. Third, the rocking curve indicates good out-of-plane crystalline alignment and suggests epitaxial growth, comparable to previously reported rocking curves for other intermetallic compounds deposited on MgO(001) including CuAl224 and CuTi.26 

The epitaxy is confirmed by XRD azimuthal φ-scans, as shown in Fig. 2(b). The plot shows at the bottom the asymmetric MgO 113 reflections obtained by tilting the sample by χ = 25.23° and fixing 2θ = 74.74°. This yields four peaks at φ = 45°, 135°, 225°, and 315°, as expected from the fourfold symmetric single-crystal MgO(001) substrate. The top pattern is acquired with χ = 54.8° and 2θ = 43.8° to detect VNi2 013 and/or 110 reflections. The VNi2 (013) and (110) planes have nearly identical interplanar spacings and are also expected at the nearly same 54.8° tilt for VNi2 [010] or [301] out-of-plane directions, respectively. The plot shows fourfold symmetric peaks at the same φ values as for the MgO 113 peaks, indicating a single in-plane orientation of the VNi2(010)/(103) layer. The peak width in φ is <1.5° for all four peaks, indicating a strong in-plane crystalline alignment. The combined results in Figs. 1(a) and 1(b) demonstrate the epitaxial growth of VNi2 on MgO(001), where either VNi2(010)||MgO(001) and VNi2[100]||MgO[110] or VNi2(103)||MgO(001) and VNi2[010]||MgO[100], as also illustrated in Fig. 3 and discussed below. Similar XRD results were obtained for all samples in this study of VNi2 layers, confirming their epitaxy.

Figure 2(c) shows the XRR results from the same layer. The measured intensity is plotted as a solid red line in a logarithmic scale as a function of scattering angle 2θ = 0.1°–3.0°. The plot includes the result from data fitting as a dotted purple line, which is offset by a factor of 5 for clarity purposes and describes the measured data well. The fitting procedure yields a thickness d = 33.5 ± 0.7 nm, which is consistent with the expected 35 nm from deposition rate calibrations. It also provides values for a surface oxide thickness of 2.8 nm and a root-mean-square (RMS) surface roughness of 0.9 nm. Similar XRR analyses were conducted for a series of stoichiometric VNi2 layers deposited at Ts = 700 °C, as shown in Fig. S4. They provide values of d = 10.5 ± 0.3, 21.0 ± 0.5, 33.5 ± 0.7, and 77.2 ± 1.5 nm for VNi2 layers, which are covered by 0.4, 1.3, 1.9, and 2.3 nm thick surface oxides with a 0.3, 0.4, 0.9, and 0.1 nm roughness. That is, all layers have surface oxides that are thinner than 3 nm, consistent with XPS depth profile measurements, and have smooth surfaces with a roughness of below 1 nm.

Figure 3 illustrates the epitaxy of VNi2 on MgO(001). The schematic in Fig. 3(a) shows the conventional unit cell of VNi2, which adopts an orthorhombic crystal structure with a = 2.52 Å, b = 3.57 Å, and c = 7.59 Å. This system exhibits almost perfect relationships between lattice constants with b 2 a and c 3 a such that each atom is 12-fold coordinated with a nearest neighbor distance of 2.52 Å. Thus, when disregarding atom types, the structure forms a slightly distorted face-centered cubic (fcc) structure, where the [301]VNi2, [010]VNi2, and [ 3 ¯ 01 ] VNi 2 directions of VNi2 are parallel to the [100]fcc, [010]fcc, and [001]fcc lattice vectors of the fcc parent structure. Figures 3(b) and 3(c) illustrate how these particular relations enable two epitaxial VNi2 layer orientations on MgO(001) with exactly the same atomic positions but a different ordering of V and Ni atoms. More specifically, Fig. 3(b) shows the case for a VNi2(010) layer. The cubes drawn with black lines indicate the conventional unit cell of the fcc structure, which exhibits a cube-on-cube epitaxy with the MgO(001) substrate such that [010]fcc and, therefore, also [010]VNi2 point in the growth direction. The V and Ni atoms (orange and gray) occupy positions above the oxygen atoms (red) of the MgO substrate. Similarly, the atoms of the VNi2(103) layer shown in Fig. 3(c) also sit on top of the oxygen atoms of the substrate, yielding again a cube-on-cube epitaxy of the parent fcc structure on MgO(001). In this case, the out-of-plane growth direction corresponds to the [100]fcc and, in turn, the [301]VNi2 crystalline direction. In summary, both VNi2(010) and VNi2(103) form an epitaxial layer on MgO(001) with a cube-on-cube epitaxy of the parent fcc structure. However, the cubes of the fcc structure, which are drawn with black lines in Figs. 3(b) and 3(c), do not represent repeating unit cells if considering the atom type (V or Ni) such that the atomic positions in the two figures are identical but not the color, which represents the atom type. Correspondingly, the structure of our VNi2/MgO(001) can be described by a single-crystal layer consisting of both VNi2(010) and VNi2(103) domains or alternatively as an epitaxial fcc V–Ni(001) alloy, which exhibits regions that exhibit atomic ordering to form locally a VNi2 compound. We note that the XRD signature for each of these descriptions is a single peak that indicates d-spacing corresponding to the lattice constant of the fcc parent structure, as discussed above. Similar structural ambiguity has previously been reported for pseudo-hexagonal AB2 layers where the ordered phase forms a MoPt2-type crystal structure.33,52 The epitaxial system illustrated in Fig. 3 exhibits a relatively large 15% lattice mismatch. Nevertheless, MgO(001) is the only investigated substrate (see also Sec. II) that leads to epitaxial VNi2 growth. We attribute the epitaxy to the preferential cube-on-cube epitaxy system, which has previously been applied to the growth of epitaxial fcc metals including Cu53,54 and Ni42 on MgO(001), despite similarly large mismatches of 14% and 16%, respectively.

Figure 4 shows the measured resistivity ρ as a function of thickness d = 10.5–138 nm of epitaxial VNi2(010)/(103) layers. The plot includes ρ measured in situ at room temperature (293 K) as red circles, in liquid N2 at 77 K as blue squares, and ex situ after air exposure at 293 K as gray open circles. The in situ ρ increases with decreasing thickness from 50.6 μ Ω cm at d = 138 nm to 103.8 μ Ω cm at d = 10.5 nm. The corresponding increase at 77 K is from 36.4 to 87.4 μ Ω cm. We attribute these increases to electron scattering at the layer surfaces. The resistivity at 77 K is smaller than that at room temperature. This is due to reduced electron–phonon scattering at cryogenic temperatures and is, to first order, independent of other scattering mechanisms such that the red and blue data points in Fig. 4 are separated by an approximately constant Δρ = 12.9–19.6 μ Ω cm, consistent with the 15 μ Ω cm from Fig. 1(a). Δρ represents the difference in the resistivity contribution from electron–phonon scattering at 293 and 77 K. This value is comparable to the previously reported first-principle prediction of 13.9 μ Ω cm for the VNi2 room temperature bulk resistivity along the [100] direction but smaller than the predictions of 25 and 70 μ Ω cm along [010] and [001].21 The plotted datapoints for the room temperature resistivity after air exposure are shifted by 2 nm to the left to account for the reduction in the effective conductive thickness during air exposure, as quantified from the measured 0.4–2.3 nm surface oxide by assuming stoichiometric VNi2 oxide with an estimated density of 4.5 g/cm3. These ex situ values are 1.8%–3.1% above the in situ resistivity, which may be attributed to an increase in diffuse electron surface scattering due to localized surface states9 or charge transfer at oxygen-exposed surfaces7 similar to what has previously been reported for Nb(001),55 Cu(001),56,57 Co(0001),8 Ni(001),42,58 and Mo.59 However, the relative resistivity increase during air exposure of our VNi2 layers is quite small, which we attribute to relatively large resistivity contributions from electron scattering at phonons and atomic disorder, i.e., domain boundaries. We note that we have deposited another VNi2 layer with a nominal thickness of 6.7 nm. However, it is not included in Fig. 4 because it shows an infinitely high resistance caused by a discontinuous microstructure, i.e., incomplete substrate coverage, similar to what has previously been reported for Ir(001) layers with d < 5 nm.18 

The curves in Fig. 4 are obtained by data fitting using the FS model.12,13 For this purpose, we assume surface scattering to be completely diffuse (p = 0), while the bulk electron mean free path λ and the bulk resistivity ρo are free fitting parameters. Setting the specularity parameter p to 0 is used because the parameters p and λ in the FS model are strongly correlated and typically cannot be independently determined from measured ρ vs d data.55,60,61 This analysis approach with p = 0 has previously been employed to quantify the resistivity scaling for various metals such as Ru,16,62 Co,8 W,63 and Ni42 and yields an effective λeff that represents a lower bound of possible λ values. Data fitting of the measured room temperature resistivity yields the red curve and ρo = 46 ± 2 μ Ω cm and λeff = 30 ± 5 nm. This bulk resistivity is approximately three times larger than the previous theoretical prediction of 13.9 μ Ω cm along the [100] direction.21 We attribute this disagreement to a higher measured resistivity caused by electron scattering at domain boundaries (i.e., atomic disorder), as illustrated in Fig. 3, but also note that the predicted resistivity along perpendicular [010] and [001] is 1.8 and 5 times higher than along [100], respectively.21 Applying the same analysis to the ex situ data results in a slightly larger λeff = 32 nm, which is consistent with the measured resistivity increase during air exposure. However, this value is within the error range of the in situ measured λeff, indicating that surface oxidation has a negligible effect on electron scattering in our layers, which may be attributed to the large competing contribution from electron scattering at domain boundaries.

Analysis of the 77 K data yields ρo = 32 ± 2 μ Ω cm and λ = 50 ± 6 nm. The cryogenic bulk resistivity is only 31% below the room temperature value, confirming that a temperature-independent electron scattering mechanism (domain walls/atomic disorder) dominates the overall resistivity. We determine the product ρoλ = (160 ± 10) × 10−16 Ω m2 at 77 K and (138 ± 5) × 10−16 Ω m2 at 293 K. This product is used as a benchmark to determine promising conductors for narrow lines. It is quite large, approximately an order-of-magnitude higher than for most competing metals, indicating that VNi2 exhibits a large resistivity size effect. We note that the ρoλ product is, to first order, independent of the magnitude of bulk scattering since additional scattering causes an increase in ρo but a corresponding decrease in λ, keeping the ρoλ product constant. Therefore, ρoλ is expected to be temperature-independent, consistent with the similar values measured at 77 and 293 K. In addition, the ρoλ is expected to be unaffected by defect scattering at antisite defects and domain boundaries. Correspondingly, a measured ρoλ should represent the intrinsic value for VNi2, independent of the level of atomic ordering. Based on this argument, our large measured ρoλ would imply that VNi2 is not promising as a conductor for narrow interconnects. However, we expect that the level of atomic ordering increases with layer thickness since columnar growth results in an increasing domain size with increasing thickness. As a result, the measured resistivity size effect presented in Fig. 4 is likely exacerbated by electron scattering at atomic disorder and domain walls, which become more pronounced for small d. That is, we expect the resistivity scaling for intrinsic VNi2 to be considerably smaller than the measured ρoλ suggests. Thus, our negative conclusion about the promise of VNi2 as interconnect material is a consequence of insufficient atomic ordering, while a different processing scheme that yields ordered and directional VNi2 may be able to confirm the promising theoretical predictions for this compound.

Sputter deposition of orthorhombic VNi2 on MgO(001) results in epitaxial layers, which exhibit domains with a VNi2(010) or VNi2(103) orientation that are separated by anti-phase boundaries but share the same fcc parent structure that forms a cube-on-cube epitaxy with the substrate. Measurements of the resistivity vs composition indicate a minimum for stoichiometric (V:Ni = 1:2) layers, while deviations from this composition result in a decrease in the order parameter, which leads to electron scattering at antisite defects and domain boundaries and a correspondingly higher resistivity. The highest atomic ordering is achieved at a growth temperature of 700 °C, while lower and higher temperatures cause increased atomic disorder attributed to kinetic barriers for ordering and an entropic driving force for disorder, respectively. The measured ρ increases with decreasing layer thickness d. Data fitting with the Fuchs–Sondheimer model yield a room temperature bulk resistivity ρo = 46 ± 2 μ Ω cm and ρoλ = (138 ± 5) × 10−16 Ω m2, while the corresponding values at 77 K are ρo = 32 ± 2 μ Ω cm and ρoλ = (160 ± 10) × 10−16 Ω m2. The measured resistivity is considerably higher than what was previously predicted. This is attributed to electron scattering at domain boundaries/atomic disorder. In addition, the resistivity scaling as quantified by the ρoλ product is higher than for many competing metals, suggesting that VNi2 is not promising as a conductor for narrow interconnects. However, the high experimental ρoλ may be attributed to a thickness-dependent atomic ordering. In that case, VNi2 is a promising interconnect material if a synthesis scheme can be developed that results in strong atomic order, a negligible domain boundary density, and a [100] crystalline orientation along the transport direction.

See the supplementary material for additional figures presenting x-ray photoelectron spectroscopy, x-ray diffraction, and x-ray reflectivity data.

The authors acknowledge the funding from the Semiconductor Research Corporation (SRC) and DARPA under Task Nos. 3137.19 and 3137.21 and from the NY State Empire State Development’s Division of Science, Technology and Innovation (NYSTAR) through Focus Center-NY–RPI Contract No. C180117.

The authors have no conflicts to disclose.

Minghua Zhang: Conceptualization (equal); Data curation (lead); Formal analysis (lead); Investigation (lead); Visualization (lead); Writing – original draft (lead); Writing – review & editing (equal). Daniel Gall: Conceptualization (equal); Funding acquisition (lead); Project administration (lead); Supervision (lead); Writing – review & editing (equal).

The data that support the findings of this study are available within the article and its supplementary material.

1.
S. M.
Rossnagel
and
T. S.
Kuan
,
J. Vac. Sci. Technol. B
22
,
240
(
2004
).
3.
S.
Salahuddin
,
K.
Ni
, and
S.
Datta
,
Nat. Electron.
1
,
442
(
2018
).
4.
T. S.
Kuan
,
C. K.
Inoki
,
G. S.
Oehrlein
,
K.
Rose
,
Y. P.
Zhao
,
G. C.
Wang
,
S. M.
Rossnagel
, and
C.
Cabral
,
MRS Proc.
612
,
711
(
2000
).
5.
G.
Steinlesberger
,
M.
Engelhardt
,
G.
Schindler
,
W.
Steinhögl
,
A.
Von Glasow
,
K.
Mosig
, and
E.
Bertagnolli
,
Microelectron. Eng.
64
,
409
(
2002
).
6.
T.
Zhou
and
D.
Gall
,
Phys. Rev. B
97
,
165406
(
2018
).
7.
A.
Jog
,
E.
Milosevic
,
P.
Zheng
, and
D.
Gall
,
Appl. Phys. Lett.
120
,
041601
(
2022
).
8.
E.
Milosevic
,
S.
Kerdsongpanya
,
M. E.
McGahay
,
A.
Zangiabadi
,
K.
Barmak
, and
D.
Gall
,
J. Appl. Phys.
125
,
245105
(
2019
).
9.
E.
Milosevic
and
D.
Gall
,
IEEE Trans. Electron Devices
66
,
2692
(
2019
).
10.
T.
Zhou
,
A.
Jog
, and
D.
Gall
,
Appl. Phys. Lett.
120
,
241603
(
2022
).
11.
M.
César
,
D.
Gall
, and
H.
Guo
,
Phys. Rev. Appl.
5
,
054018
(
2016
).
12.
K.
Fuchs
,
Mathematical Proceedings of the Cambridge Philosophical Society
(Cambridge University Press,
1938
), pp.
100
108
.
13.
E. H.
Sondheimer
,
Adv. Phys.
1
,
1
(
1952
).
14.
A. F.
Mayadas
and
M.
Shatzkes
,
Phys. Rev. B
1
,
1382
(
1970
).
15.
D.
Gall
,
J. Appl. Phys.
127
,
050901
(
2020
).
16.
E.
Milosevic
,
S.
Kerdsongpanya
,
A.
Zangiabadi
,
K.
Barmak
,
K. R.
Coffey
, and
D.
Gall
,
J. Appl. Phys.
124
,
165105
(
2018
).
17.
A.
Jog
and
D.
Gall
,
IEEE Trans. Electron Devices
69
,
3854
(
2022
).
18.
A.
Jog
and
D.
Gall
,
J. Appl. Phys.
130
,
115103
(
2021
).
19.
A.
Jog
,
T.
Zhou
, and
D.
Gall
,
IEEE Trans. Electron Devices
68
,
257
(
2021
).
20.
J.
Chawla
and
D.
Gall
,
Appl. Phys. Lett.
94
,
252101
(
2009
).
21.
S.
Kumar
,
C.
Multunas
,
B.
Defay
,
D.
Gall
, and
R.
Sundararaman
,
Phys. Rev. Mater.
6
,
085002
(
2022
).
22.
L.
Chen
,
S.
Kumar
,
M.
Yahagi
,
D.
Ando
,
Y.
Sutou
,
D.
Gall
,
R.
Sundararaman
, and
J.
Koike
,
J. Appl. Phys.
129
,
035301
(
2021
).
23.
J.-P.
Soulié
,
Z.
Tőkei
,
J.
Swerts
, and
C.
Adelmann
, in
2022 IEEE International Interconnect Technology Conference (IITC)
(IEEE,
2022
), pp.
73
75
.
24.
M.
Zhang
and
D.
Gall
,
IEEE Trans. Electron Devices
69
,
5110
(
2022
).
25.
L.
Chen
,
Q.
Chen
,
D.
Ando
,
Y.
Sutou
,
M.
Kubo
, and
J.
Koike
,
Appl. Surf. Sci.
537
,
148035
(
2021
).
26.
M.
Zhang
,
S.
Kumar
,
R.
Sundararaman
, and
D.
Gall
,
J. Appl. Phys.
133
,
045102
(
2023
).
27.
T.
Kuge
,
M.
Yahagi
, and
J.
Koike
, in
2021 IEEE International Interconnect Technology Conference (IITC)
(IEEE,
2021
), pp.
1
3
.
28.
T.
Kuge
,
M.
Yahagi
, and
J.
Koike
, in
2022 IEEE International Interconnect Technology Conference (IITC)
(IEEE,
2022
), pp.
22
24
.
29.
M.
Zhang
,
S.
Kumar
,
R.
Sundararaman
, and
D.
Gall
,
J. Appl. Phys.
130
,
034302
(
2021
).
30.
J. E.
Yoo
,
J. Y.
Sung
,
J. H.
Hwang
,
I.
Maeng
,
S.-J.
Oh
,
I.
Lee
,
J. H.
Shim
,
S. D.
Kim
,
D.-S.
Yoon
, S. Y. Jang, Y. J. Kang, and S. W. Lee,
ACS Appl. Mater. Interfaces
13
,
61809
(
2021
).
31.
K.
Sankaran
,
K.
Moors
,
Z.
Tőkei
,
C.
Adelmann
, and
G.
Pourtois
,
Phys. Rev. Mater.
5
,
056002
(
2021
).
32.
S.
Kumar
,
C.
Multunas
, and
R.
Sundararaman
,
Phys. Rev. Mater.
6
,
125201
(
2022
).
33.
B. C.
Giessen
and
N. J.
Grant
,
J. Less-Common Met.
8
,
114
(
1965
).
34.
C.
Wu
,
X.
Su
,
J.
Wang
,
Y.
Liu
,
H.
Tu
, and
Y.
Zhang
,
J. Electron. Mater.
43
,
4111
(
2014
).
35.
B.
Wang
and
D.
Gall
,
Thin Solid Films
688
,
137165
(
2019
).
36.
P.
Fang
,
C. P.
Mulligan
,
R.
Jia
,
J.
Shi
,
S. V.
Khare
, and
D.
Gall
,
Acta Mater.
226
,
117643
(
2022
).
37.
P.
Fang
,
B.
Wang
, and
D.
Gall
,
Surf. Coat. Technol.
420
,
127333
(
2021
).
38.
P. L.
Rossiter
,
The Electrical Resistivity of Metals and Alloys
(
Cambridge University Press
,
Cambridge
,
1991
).
39.
H.
Lang
,
H.
Uzawa
,
T.
Mohri
, and
W.
Pfeiler
,
Intermetallics
9
,
9
(
2001
).
40.
D. R.
Gaskell
and
D. E.
Laughlin
,
Introduction to the Thermodynamics of Materials
(
CRC Press
,
Boca Raton
,
2017
).
41.
T.
Ikeda
,
H.
Numakura
, and
M.
Koiwa
,
Acta Mater.
46
,
6605
(
1998
).
42.
E.
Milosevic
,
P.
Zheng
, and
D.
Gall
,
IEEE Trans. Electron Devices
66
,
4326
(
2019
).
43.
I. S.
Beloborodov
,
A. V.
Lopatin
,
V. M.
Vinokur
, and
K. B.
Efetov
,
Rev. Mod. Phys.
79
,
469
(
2007
).
44.
I. S.
Beloborodov
,
A. V.
Lopatin
, and
V. M.
Vinokur
,
Phys. Rev. B
70
,
205120
(
2004
).
45.
I. S.
Beloborodov
,
A. V.
Lopatin
, and
V. M.
Vinokur
,
Phys. Rev. B
72
,
125121
(
2005
).
46.
X. Y.
Zhang
,
J. S.
Chawla
,
R. P.
Deng
, and
D.
Gall
,
Phys. Rev. B
84
,
073101
(
2011
).
47.
X. Y.
Zhang
,
J. S.
Chawla
,
B. M.
Howe
, and
D.
Gall
,
Phys. Rev. B
83
,
165205
(
2011
).
48.
J. P. Neumann, Y. A. Chang, C. M. Lee
,
Acta metall.
24
,
593
(
1976
).
49.
J. T.
Okada
,
T.
Ishikawa
,
Y.
Watanabe
, and
P.-F.
Paradis
,
J. Chem. Thermodyn.
42
,
856
(
2010
).
50.
A. L.
Ortiz
,
J.-W.
Tian
,
L. L.
Shaw
, and
P. K.
Liaw
,
Scr. Mater.
62
,
129
(
2010
).
51.
P.
Zheng
,
B. D.
Ozsdolay
, and
D.
Gall
,
J. Vac. Sci. Technol. A
33
,
061505
(
2015
).
52.
K.
Young
,
M.
Young
,
S.
Chang
, and
B.
Huang
,
J. Alloys Compd.
560
,
33
(
2013
).
53.
J. M.
Purswani
,
T.
Spila
, and
D.
Gall
,
Thin Solid Films
515
,
1166
(
2006
).
54.
J. M.
Purswani
and
D.
Gall
,
Thin Solid Films
516
,
465
(
2007
).
55.
E.
Milosevic
,
S.
Kerdsongpanya
,
M. E.
McGahay
,
B.
Wang
, and
D.
Gall
,
IEEE Trans. Electron Devices
66
,
3473
(
2019
).
56.
P.
Zheng
,
T.
Zhou
, and
D.
Gall
,
Semicond. Sci. Technol.
31
,
055005
(
2016
).
57.
J. S.
Chawla
,
F.
Zahid
,
H.
Guo
, and
D.
Gall
,
Appl. Phys. Lett.
97
,
132106
(
2010
).
58.
P.Y. Zheng, R.P. Deng, and D. Gall
,
Appl. Phys. Lett.
105
,
131603
(
2014
).
59.
A.
Jog
,
P.
Zheng
,
T.
Zhou
, and
D.
Gall
,
Nanomaterials
13
,
957
(
2023
).
60.
J.S. Chawla, X. Y. Zhang and D. Gall
,
J. Appl. Phys.
110
,
043714
(
2011
).
61.
E.
Milosevic
and
D.
Gall
,
AIP Adv.
10
,
055213
(
2020
).
62.
E.
Milosevic
,
S.
Kerdsongpanya
, and
D.
Gall
, in
2018 IEEE Nanotechnology Symposium (ANTS)
(IEEE,
2018
), pp.
1
5
.
63.
P.
Zheng
and
D.
Gall
,
J. Appl. Phys.
122
,
135301
(
2017
).

Supplementary Material