We present vanadium dioxide (VO2) thin films having high resistivity contrast with silicon substrates through use of crystallized alumina (Al2O3) buffer layers, engineered for this purpose. We first optimized the process by depositing VO2 onto C-plane sapphire substrates prior to alumina thin films. The latter of which were grown via atomic layer deposition on silicon substrates. We then applied rapid thermal annealing (RTA) to crystallize the alumina films. Scanning electron microscopy results indicated a thickness of 107 nm for each VO2 film, which yielded hot–cold resistivity contrast ratios of 9.76 × 104, 1.46 × 104, and 3.66 × 103, when deposited on the C-plane sapphire, the annealed buffers, and the as-deposited alumina buffers, respectively. Atomic force microscopy of the film surface roughness of the VO2 films indicated root mean squared roughness (Rq) of 4.56 nm, 6.79 nm, and 3.30 nm, respectively, for the films grown on the C-plane sapphire, annealed buffers, and as-deposited buffers. Finally, x-ray diffraction (XRD) of the VO2 films indicated the desired composition and strong (0h0)/(00h) texturing, when deposited on both the C-plane sapphire and the annealed alumina buffer layers. XRD results indicated a series of peaks corresponding to the α-Al2O3/C-plane sapphire, and an XRD analysis of the buffers alone confirmed crystallization of the buffer layer via RTA. The process defined in this paper produced a series of highly textured VO2 films making them most valuable for the integration of VO2 with silicon-based devices.

The versatility and quasi-passive attributes of phase change materials (PCMs) in terms of reacting directly or “instinctually” to external stimuli (e.g., heat, pressure, light, and strain) make them particularly useful in sensing and switching applications. Such factors can mean replacing specialized traditional integrated circuits, or improving sensor performance by either exploiting nonlinear behavior, or by creating reconfiguration by activating multiple states of geometry through the strategic actuation of different PCM regions. Vanadium dioxide (VO2) is particularly appealing for this exploitation/reconfiguration in that it has a low metal-insulator transition (MIT) temperature (68 °C), which exhibits a reversible, sharp, nonlinear change. These properties are advantageous in terms of lower power and faster switching devices, unlike other PCMs (e.g., germanium-antimony-telluride), that require a high temperature excitation (i.e., ∼618 °C) to revert to a dielectric state.1 Although the temperature and hysteresis width of the VO2 phase transition is tunable, as is the MIT itself on a domain level,2 this work focuses on growing VO2 films using silicon (Si) substrates that exhibit similar properties to that grown on a C-plane sapphire.

The C-plane sapphire was chosen because it yields VO2 films whose electrical properties closely resemble those of bulk VO2: particularly, a transition temperature of 68 °C.3 Other works have presented comparisons of high-quality VO2 films on C-plane with R-3,4 and M-plane4 sapphire, but both alternate sapphire substrates show a shift in transition temperature. All experimental VO2 films in this work maintain the transition temperature of 68 °C so the most direct control comparison is the C-plane sapphire.

VO2 applications include radio frequency (RF) switches,5,6 infrared asymmetric transmission limiters,7 and thermal imaging sensors.8 Its sharp change in resistivity makes it an ideal material for a shunt RF switch, where the room temperature VO2 is RF-transparent in the passing state and heated to a conducting state to short the signal line to the ground plane, causing high reflection.6,9 Such switches and similar mechanisms are useful in a multitude of RF applications, including millimeter wave phase shifters10 and reactive power limiters,6 reconfigurable antennas,11 and metasurfaces.12 As Si is widely used in a multitude of such research and industry concepts, additional applications of VO2 such as micro-actuators,2,13,14 thermochromics,2,15,16 optical modulators,2,17 field-effect transistors,18 power meters,2 electrostatic discharge protection,19 and gas sensors2 give merit to the integration of VO2 on Si substrates. The key novelty of this work is the deposition of high-contrast VO2 films on Si substrates through use of annealed Al2O3 buffer layers, thus paving the way for the easy integration of such devices in the future.

Unfortunately, the limited deposition of state of the art VO2 thin films to single crystal sapphire substrates,4 where there is a strong lattice match between the C-plane sapphire and VO2, currently hinders the pragmatic use of VO2. Specifically, the extremely low etch rates20 and lack of semiconductor properties prevent the use of sapphire for bulk etch/micromachining and in standard device processes. Indeed, silicon substrates are most commonly used in device processing, for both traditional integrated circuit applications and techniques used in micro-electromechanical systems. Unfortunately, the large lattice mismatch between Si and VO2 inhibits the efficacy of direct deposition of VO2. Such films exhibit comparatively small change in the resistivity (<1 × 103) and large hysteresis (∼8 °C) in the phase transition region.21 In cases where larger hysteresis is beneficial, like memory applications, a recent overview discusses intentional manipulation of hysteresis.2 Intuitively, it is better to start with VO2 having bulk properties and modify from there than to accept a VO2 film with weak orientation and low phase purity. In Table I, we provide a comparison of theoretical lattice mismatch of VO2 on both the C-plane sapphire and Si.

TABLE I.

Lattice mismatch of tetragonal rutile VO2 as deposited on a C-plane sapphire and (100) silicon.

C-plane sapphire (α-Al2O3)Silicon-d (Fd3¯m)
VO2 (P42/mnm4.19% 16.1% 
C-plane sapphire (α-Al2O3)Silicon-d (Fd3¯m)
VO2 (P42/mnm4.19% 16.1% 

To improve the performance of VO2 thin films on Si substrates, several materials for buffer layers such as yttria-stabilized zirconia,22 HfO2,23 and amorphous Al2O321,24 were proposed. The advantage of this work over previous studies using Al2O3 buffer layers is the new comparison of VO2 on amorphous alumina with VO2 on annealed, polycrystalline alumina buffers. Another distinguishing factor is that this work also uses DC sputtering rather than RF sputtering21 and yields a higher resistivity contrast ratio.21,24 Although thermally oxidized silicon dioxide buffer layers25 have shown resistivity contrast ratios of 1.3 × 104, this fabrication technique is limited to Si and cannot be adapted to other substrates. It is, however, possible to use the thin film atomic layer deposition (ALD) technique with various substrates such as GaAs or GaN.26,27

VO2 films deposited on single crystal sapphire substrates have been described as epitaxial, since their sapphire lattice determines their orientation.28–30 The VO2 films presented in this work have clear epitaxial relation to their underlying sapphire substrates and annealed buffers, but we refrain from using the term epitaxial for two reasons; the sputtering process described here is not a standard epitaxy process, and in the buffer layer case, the underlying growth surface is polycrystalline, not monocrystalline. However, we also abstain from using the term polycrystalline to describe the VO2, as it discounts the substantial dependence of our VO2 preferred texturing on the underlying crystal substrate and buffer layer structures. Therefore, we use the modified term textured-epitaxial to convey the dependence of the VO2 films' structure on epitaxial interactions with both the single crystal sapphire and local sites in the polycrystalline alumina buffer. As such, we sought to create a polycrystalline α-Al2O3 buffer layer to create similar VO2 growth conditions on the bulk sapphire. The alumina buffer is not categorized as, nor intended to be epitaxial nor textured-epitaxial in itself, but rather it is intended to eliminate the effect of the Si diamond cubic lattice on the VO2. The VO2, Si, and Al2O3 lattices are illustrated conceptually as cross sections in Fig. 1.

FIG. 1.

(a) Crystallographic schematic of rutile VO2 on diamond silicon illustrating the large lattice parameter mismatch that leads to high strain vs (b) rutile VO2 on an alumina buffer layer on silicon eases the lattice parameter mismatch lowering the overall strain and increasing the electrical resistivity contrast.

FIG. 1.

(a) Crystallographic schematic of rutile VO2 on diamond silicon illustrating the large lattice parameter mismatch that leads to high strain vs (b) rutile VO2 on an alumina buffer layer on silicon eases the lattice parameter mismatch lowering the overall strain and increasing the electrical resistivity contrast.

Close modal

Our proposed solution is to implement alumina buffer layers deposited using ALD, which we crystallize using rapid thermal annealing (RTA). X-ray diffraction (XRD) illustrates no crystallographic differences of the VO2 structure when deposited on annealed polycrystalline alumina buffer vs C-plane sapphire, whereas VO2 deposited on amorphous alumina buffer indicates significantly reduced crystallinity. Furthermore, XRD analysis of alumina thin films deposited on diamond silicon without VO2 confirms the emergence of polycrystallinity upon annealing. The VO2 films used here were deposited with DC reactive sputtering, which is advantageous over a number of techniques such as pulsed laser deposition, solgel, ALD, chemical vapor deposition, molecular beam epitaxy,16 metal organic chemical vapor deposition,30 and hydrothermal growth31 making it possible to extend the sputtering to a larger batch processing for efficient manufacturing.

1. Vanadium dioxide deposition

This VO2 deposition process is an extension of previous efforts to create reconfigurable terahertz filters with VO2 using reactive sputtering.32 An incrementing of the O2 flow during annealing and controlling for a variation in the deposition rate between trials increased the resistivity contrast ratio from 7 × 104 to 9.76 × 104 and decreased the hysteresis width from 6 °C to 4 °C. As the measurement process was also improved, the 6 °C hysteresis figure refers to early trials within this work. The C-plane sapphire wafers are cleaned with organic solvents and dried. A dummy Si wafer is then used to eliminate contaminants in the sputtering chamber and coat the chamber in VO2 prior to deposition on the experimental sample.

First, the sample is cleaned with plasma for 7 min using 20 sccm Ar+ plasma and 50 W RF power at 4.0 Pa. Next, the chamber is heated at vacuum to 650 °C prior to flowing 1.62 sccm O2 and 20 sccm Ar (7.5% O2: total gas ratio) at 1.33 Pa pressure. The DC power is increased to 250 W and the sputtering target shutter left open for the time required to ensure the deposition of 107 nm of VO2, based upon the rate measured. Once the deposition is complete, the sample is annealed, in situ, at the same temperature in 1.33 Pa, 10 sccm O2 environment for 15 min. We use an experimental process, the iterative results of which are detailed in Sec. III A, to establish these parameters. Overall, we established them through five key iterations, described in Table II in Sec. III A.

TABLE II.

VO2 deposition trials. Rather than listing the constant process variables in the table, they are described in the text for clarity. Trials deemed that process failures are also omitted and key process changes marked with *.

Sample IDAr flow (sccm)O2 flow (sccm)O2:total gas ratio (%)Deposition time (min:s)Annealing O2 flow (sccm)Resistivity contrast ratio (unitless)Hysteresis width (°C)Transition range (°C)
19.96 1.58 7.34 16:40 1.58 3.11 × 104 12 
19.96 *1.62 *7.51 16:40 *1.62 1.99 × 104 18 
19.96 *1.58 *7.34 16:40 *4.96 6.29 × 104 12 
19.96 1.58 7.34 16:40 *9.97 6.88 × 104 12 
19.96 1.58 7.34 *20:00 9.97 9.76 × 104 12 
Sample IDAr flow (sccm)O2 flow (sccm)O2:total gas ratio (%)Deposition time (min:s)Annealing O2 flow (sccm)Resistivity contrast ratio (unitless)Hysteresis width (°C)Transition range (°C)
19.96 1.58 7.34 16:40 1.58 3.11 × 104 12 
19.96 *1.62 *7.51 16:40 *1.62 1.99 × 104 18 
19.96 *1.58 *7.34 16:40 *4.96 6.29 × 104 12 
19.96 1.58 7.34 16:40 *9.97 6.88 × 104 12 
19.96 1.58 7.34 *20:00 9.97 9.76 × 104 12 

We discovered that the deposition rate can vary significantly between trials when the same process was employed multiple times with varying results. However, controlling for the variable deposition rate (by measuring before each deposition) returned the expected results, thus confirming repeatability of this process. Therefore, the final iteration of the process is also used to deposit on the experimental buffer layers as well.

2. Alumina deposition and annealing

To fabricate the alumina thin films, Al2O3 with a thickness of 45 nm is deposited onto the Si substrates. The precursors are trimethylaluminum and H2O heated to 300 °C with N2 as the purging gas, followed by the RTA process at 950 °C for 60 s in N2 ambient. This process is essential for the deposition of high-contrast, highly textured VO2 films since the crystallization of alumina creates a stronger basis for VO2 growth than do the amorphous films.

Electrical resistivity, crystallinity, thickness, and surface morphology serve as the criteria for characterizing the VO2 films grown in this study. To capture the full heating and cooling hysteresis loop for each sample, the resistivity was measured using a 4-point probe, hot plate, and surface probe thermometer. Parallel-beam x-ray diffraction (XRD) is next used to determine the crystal structure of the VO2 and for the Al2O3 buffers. All XRD patterns are collected on a Rigaku SmartLab diffractometer with Cu-Kα radiation, in a 2θ/θ collection mode, with a step size of 0.01°, a speed of 12 s/step, and all slits set to 1.00 mm. Scanning electron microscopy (SEM) is used to measure thickness and observe crystal grain structure and surface morphology, and atomic force microscopy (AFM) used to measure the surface roughness.

Resistivity contrast ratio, defined as the resistivity at room temperature divided by that above the MIT, is a key figure of merit for VO2 thin films. A hot plate, surface probe thermometer, and 4-point probe were used to measure the sheet resistance over the phase change of the VO2 samples. We measured the samples from room temperature through and past the transition region and again during the cool down to capture the full hysteretic loop of the phase change behavior. We then multiplied the sheet resistance values by the thickness to calculate resistivity.

As indicated in the resistivity hysteresis loops in Fig. 2, on the C-plane sapphire, our films exhibit state of the art performance, with nearly five orders of magnitude resistivity contrast ratio at 9.76 × 104. The iterations of the fabrication process are detailed in Table II with their respective contrast ratio and hysteresis width results. A significantly decreased resistivity contrast ratio at 3.66 × 103, characterized the VO2 thin films deposited on as-deposited alumina buffer layers. Such findings are expected in that ALD alone yields amorphous films that cannot create a lattice match for VO2 growth. Based upon our hypothesis that RTA is suitable in transforming the alumina buffer layer from amorphous to polycrystalline, we deposited VO2 on these annealed buffer layers. The result was an increase in the resistivity contrast ratio by a factor of 4, which yielded a new ratio of 1.46 × 104.

FIG. 2.

Resistivity vs temperature measurements comparing (a) the electrical performance of the VO2 thin films on ALD thin films, both amorphous (as-deposited) and annealed, as well as bulk crystalline C-plane sapphire wafers and (b) each deposition trial iteration as detailed in Table I. Sample IDs are labeled A through E, corresponding to the process details in Table I. Sample A is the first trial, and Sample E is from the finalized process.

FIG. 2.

Resistivity vs temperature measurements comparing (a) the electrical performance of the VO2 thin films on ALD thin films, both amorphous (as-deposited) and annealed, as well as bulk crystalline C-plane sapphire wafers and (b) each deposition trial iteration as detailed in Table I. Sample IDs are labeled A through E, corresponding to the process details in Table I. Sample A is the first trial, and Sample E is from the finalized process.

Close modal

Parallel-beam XRD analysis confirmed that depositing VO2 on annealed polycrystalline alumina does not cause a change in the crystal structure compared to when it is deposited on the C-plane sapphire, with both yielding monoclinic (P21/c) VO2 at room temperature. Figure 3 illustrates low-crystallinity monoclinic VO2 when deposited on amorphous alumina with a low intensity (020)/(002) peak around 39.8°. However, when deposited on annealed polycrystalline alumina and C-plane sapphire, there are strong, clear diffraction peaks at 39.8°, corresponding to VO2 texturing in the (020) or (002) planes supported by the peaks at approximately 86°, corresponding to the (040) or (004) plane. Because the d-spacings of the (002), and (020) peaks as well as the (004), and (040) peaks are almost exactly coincidental, they cannot be differentiated in this circumstance.3,33–41 These 2θ ranges were chosen because no other diffraction peaks were shown during previous XRD scans on similar samples, which can be found in figures in the supplementary material. Other literature corroborates these two peaks and shows that for (020)/(002) textured VO2, the rest of the spectrum does not show any peaks of interest.3,4

FIG. 3.

XRD of (a) monoclinic VO2 deposited on a C-plane sapphire, 45nm annealed alumina, and 45nm as-deposited alumina, illustrating higher crystallinity VO2 and (020)/(002) texture when deposited on the annealed alumina (b) 25nm alumina before and after annealing illustrating polycrystallinity emerging on top of the diamond silicon substrate with a silicon intermediate, Ia3¯.

FIG. 3.

XRD of (a) monoclinic VO2 deposited on a C-plane sapphire, 45nm annealed alumina, and 45nm as-deposited alumina, illustrating higher crystallinity VO2 and (020)/(002) texture when deposited on the annealed alumina (b) 25nm alumina before and after annealing illustrating polycrystallinity emerging on top of the diamond silicon substrate with a silicon intermediate, Ia3¯.

Close modal

VO2 on the C-plane sapphire has a slight (040)/(004) peak shift of 0.5° toward a higher 2θ from simulated bulk diffraction patterns; whereas, VO2 on the annealed polycrystalline alumina (040)/(004) peak shift of 0.9° toward a higher 2θ angle. However, both (040)/(004) peaks are broadened. This shift is most likely due to the highly non-uniform but overall compressive strain on the system, which is also supported by the (040)/(004) peak broadening. The peak broadening depends on the Bragg angle, indicating that it is most likely due to a combination of the increase in strain, which could be different in individual grains, and crystallite size. This broadening originates from the lattice mismatch between the Si-d substrate (Fd3¯m) and the α-Al2O3 buffer (R3¯cH), shown in Fig. 3, which strains the buffer placing extra strain on the VO2 that is absent when deposited on the single-crystalline C-plane sapphire. Additionally, there is strain inherent in the VO2 films; since they are deposited above the MIT, they are grown in the tetragonal rutile structure and subsequently transition to monoclinic when cooled. Figure 3 illustrates the initial crystallization of VO2 with this preferred (020) growth texture even without the alumina anneal. However, the low signal-to-noise ratio and absence of the (040) peak indicates the presence of very little crystalline VO2 without annealing the buffer.

The small peaks at 38.1° in both the C-plane sapphire and annealed alumina and at 41.7° in the C-plane sapphire correspond to the (110) and (006) planes, respectively, in the α-Al2O3 substrate and buffer. However, when Al2O3 is deposited and annealed on the silicon substrate by itself, it crystallizes into γ-Al2O3(Fd3¯m) as shown in Fig. 2(b). The 19.6°, 32.3°, and 39.3° correspond to the (111), (220), and (222) planes in γ-Al2O3, respectively. As with VO2 deposition, there are slight peak shifts of 0.2° and 0.4° for the (211) and (220) peaks as compared to simulated bulk diffraction patterns associated with the lattice mismatch placing strain on the system. Since the RTA performed on the alumina buffer creates only γ-Al2O3, it is likely that some local topotaxial conversion28,42 in the alumina occurs during the VO2 deposition process from the γ-phase to α-phase.

In the XRD scans of alumina buffers on Si, the highest intensity peak coincides with the (100) single-crystalline Si-d (Fd3¯m) substrate. However, the second highest intensity peak at the 33.2° corresponds to the (121) plane of Si-HP (Ia3¯), which is supported by the presence of two peaks at approximately 91° indexing to the (532) plane of Si-HP. The presence of this phase most likely manifests as a lattice transition layer between the Si-d substrate and the crystallized γ-Al2O3 due to a large lattice mismatch of 46.2%. Similar to the VO2 crystallized on annealed Al2O3, the slight peak shifts in the XRD patterns are due to the strain, with peak broadening due to a combination of larger crystallite size and the non-uniformity of the strain from individual grains in the annealed polycrystalline buffer.

In Fig. 4, we show the theoretical lattice mismatch values of the deposited materials with the substrates and buffers. The low lattice mismatch of the tetragonal rutile VO2 (R), which is associated with the deposition temperature, is the central consideration for forming highly textured VO2. We perform the deposition process at 650 °C, which is well above the MIT, to ensure the creation of a VO2 (R) film during the growth and annealing process. This VO2 (R) film evolves into VO2 (M) upon cooling. From these observations, we conclude that the annealed alumina ultimately mimics to a significant degree the growth conditions leading to a highly crystalline VO2 (M), which the C-plane sapphire provides.

FIG. 4.

Absolute lattice mismatch of the deposited compounds (left) on the substrates or buffers (top). VO2 (M) is included, since that is the phase observed in XRD. However, it is not the growth lattice, and therefore its mismatch is not necessarily considered.

FIG. 4.

Absolute lattice mismatch of the deposited compounds (left) on the substrates or buffers (top). VO2 (M) is included, since that is the phase observed in XRD. However, it is not the growth lattice, and therefore its mismatch is not necessarily considered.

Close modal

With the similarities in mind between VO2 as deposited on the C-plane sapphire and annealed alumina, we classify the latter as textured-epitaxial as well. The polycrystalline nature of the buffer and the sputtering/in situ annealing process preclude using the term epitaxial, but the data strongly suggest local epitaxial interactions as the cause of the similarity in these films.

Each VO2 film is 107 nm thick, with negligible variation based on the crystallinity of the buffer layer. To confirm the thickness of the films, we observed a cross section of each sample deposited on an alumina buffer layer using SEM. As shown in Fig. 5, the grains of the VO2 thin films are evident in both the samples with annealed and as-deposited buffer layers, with the film thickness not dependent on the buffer layer properties. Note also the crystallization of the buffer layer as hypothesized and observed using XRD by the visual inspection of the micrographs compared in Fig. 5. This observation of the crystal formation further corroborates the hypothesis that the application of RTA to ALD alumina thin films creates a basis for higher contrast, more highly textured VO2 than as-deposited alumina.

FIG. 5.

SEM micrographs of VO2 deposited on [(a) and (b)] annealed Al2O3 and [(c) and (d)] as-deposited Al2O3 buffer layers. Each buffer layer is 45nm thick, and each VO2 film is 107nm thick. The left column [(a) and (c)] shows cross sections of the films, and the right column [(b) and (d)] shows the same samples on a 45° tilt.

FIG. 5.

SEM micrographs of VO2 deposited on [(a) and (b)] annealed Al2O3 and [(c) and (d)] as-deposited Al2O3 buffer layers. Each buffer layer is 45nm thick, and each VO2 film is 107nm thick. The left column [(a) and (c)] shows cross sections of the films, and the right column [(b) and (d)] shows the same samples on a 45° tilt.

Close modal

Additionally, we performed AFM to study the VO2 films' surface roughness. The VO2 grown on as-deposited alumina showed the lowest root mean squared roughness (Rq), of 3.30 nm, the results of which were expected as it was the most amorphous film as shown by XRD and SEM. The Rq roughness of the VO2 on annealed alumina and C-plane sapphire were 6.79 nm and 4.56 nm, respectively. Although VO2 on the C-plane sapphire exhibits a more predominant crystalline arrangement, as shown by XRD, it was smoother, as indicated by the difference in roughness between the sapphire wafer and the annealed alumina buffer. Sapphire is polished to an arithmetic mean (Ra) roughness of <0.3 nm according to the manufacturer, while the Ra roughness of the annealed alumina thin film is 0.68 nm, measured using AFM. Note also in Fig. 6 that VO2 on as-deposited alumina has a rounded grain structure and the smallest size, while VO2 on annealed alumina has larger grains with some variance in size and shape, and VO2 on the C-plane sapphire has consistently sized and shaped grains.

FIG. 6.

AFM images of 1 × 1 μm (top row) and 5 × 5 μm (bottom row) of VO2 on the C-plane sapphire [(a) and (d)]; the annealed alumina buffer [(b) and (e)]; and the as-deposited alumina [(c) and (f)].

FIG. 6.

AFM images of 1 × 1 μm (top row) and 5 × 5 μm (bottom row) of VO2 on the C-plane sapphire [(a) and (d)]; the annealed alumina buffer [(b) and (e)]; and the as-deposited alumina [(c) and (f)].

Close modal

We developed a repeatable process for depositing high-resistivity contrast, textured-epitaxial VO2 on C-plane sapphire and on alumina thin film buffer layers on silicon substrates, using DC sputtering. The textured-epitaxial VO2 on sapphire exhibited state of the art resistivity contrast of 9.76 × 104, the polycrystalline VO2 on as-deposited alumina showed resistivity contrast of 3.66 × 103, and annealing the alumina increased that result by a factor of 4 to 1.46 × 104 in the respective textured-epitaxial VO2. XRD, SEM, and AFM analysis of our films shows that the annealed alumina buffer layers produce textured-epitaxial VO2 films similar to those grown directly on the bulk C-plane sapphire, as well as crystallization of the alumina, particularly in the γ-phase. Our VO2 deposition process, combined with silicon substrates, creates high-contrast, textured-epitaxial films; thus, it is more easily integrated with a much wider range of devices than VO2 strictly grown on sapphire wafers.

See the supplementary material for a full 2θ angle range (10°–70°) x-ray diffraction scan of VO2 on a C-plane sapphire substrate.

The authors wish to acknowledge the support of the National Science Foundation (NSF) Grant No. 1845370, the Air Force Research Lab (AFRL)/Defense Associated Graduate Student Innovators (DAGSI) Award No. RY6-OSU-19-2-AFRL2, and the Ohio State University Institute for Materials Research. The Center for Electron Microscopy and Analysis (CEMAS) at the Ohio State University was the location for the x-ray diffraction performed in this effort. C.E.W. and V.D.N. acknowledge financial support from the Smart Vehicles Concepts Center, an NSF Industry–University Cooperative Research Center (NSF IIP 1738723) for seed funding.

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

1.
A.
Sarangan
,
J.
Duran
,
V.
Vasilyev
,
N.
Limberopoulos
,
I.
Vitebskiy
, and
I.
Anisimov
,
IEEE Photonics J.
10
,
2200409
(
2018
).
2.
R.
Shi
,
N.
Shen
,
J.
Wang
,
W.
Wang
,
A.
Amini
,
N.
Wang
, and
C.
Cheng
,
Appl. Phys. Rev.
6
,
011312
(
2019
).
3.
T.-H.
Yang
,
R.
Aggarwal
,
A.
Gupta
,
H.
Zhou
,
R. J.
Narayan
, and
J.
Narayan
,
J. Appl. Phys.
107
,
053514
(
2010
).
4.
Y.
Zhao
,
J.
Hwan Lee
,
Y.
Zhu
,
M.
Nazari
,
C.
Chen
,
H.
Wang
,
A.
Bernussi
,
M.
Holtz
, and
Z.
Fan
,
J. Appl. Phys.
111
,
053533
(
2012
).
5.
F.
Dumas-Bouchiat
,
C.
Champeaux
,
A.
Catherinot
,
A.
Crunteanu
, and
P.
Blondy
,
Appl. Phys. Lett.
91
,
223505
(
2007
).
6.
A.
Crunteanu
,
J.
Givernaud
,
J.
Leroy
,
D.
Mardivirin
,
C.
Champeaux
,
J.-C.
Orlianges
,
A.
Catherinot
, and
P.
Blondy
,
Sci. Technol. Adv. Mater.
11
,
065002
(
2010
).
7.
R.
Thomas
,
E.
Makri
,
T.
Kottos
,
B.
Shapiro
, and
I.
Vitebskiy
,
Phys. Rev. A
98
,
053806
(
2018
).
8.
C.
Chen
,
X.
Yi
,
X.
Zhao
, and
B.
Xiong
,
Sens. Actuators A
90
,
212
(
2001
).
9.
J.
Lee
,
D.
Lee
,
S. J.
Cho
,
J.-H.
Seo
,
D.
Liu
,
C.-B.
Eom
, and
Z.
Ma
,
Appl. Phys. Lett.
111
,
063110
(
2017
).
10.
C.
Hillman
,
P. A.
Stupar
,
J. B.
Hacker
,
Z.
Griffith
,
M.
Field
, and
M.
Rodwell
, in IEEE MTT-S International Microwave Symposium Digest, Tampa, FL (2014).
11.
D.
Anagnostou
,
D.
Torres
,
T.
Teeslink
, and
N.
Sepulveda
,
IEEE Antennas Propag. Mag.
(published online 2020).
12.
M. R. M.
Hashemi
,
S.-H.
Yang
,
T.
Wang
,
N.
Sepulveda
, and
M.
Jarrahi
,
Sci. Rep.
6
,
35439
(
2016
).
13.
R.
Cabrera
,
E.
Merced
, and
N.
Sepulveda
,
J. Microelectromech. Syst.
23
,
243
(
2014
).
14.
E.
Merced
,
N.
Dávila
,
D.
Torres
,
R.
Cabrera
,
F. E.
Fernández
, and
N.
Sepúlveda
,
Smart Mater. Struct.
21
,
105009
(
2012
).
15.
Y.
Choi
,
Y.
Jung
, and
H.
Kim
,
Thin Solid Films
615
,
437
(
2016
).
16.
H.-C.
Ho
,
Y.-C.
Lai
,
K.
Chen
,
T. D.
Dao
,
C.-H.
Hsueh
, and
T.
Nagao
,
Appl. Surf. Sci.
495
,
143436
(
2019
).
17.
L.
Sánchez
,
S.
Lechago
, and
P.
Sanchis
,
Opt. Lett.
40
,
1452
(
2015
).
18.
K.
Abbas
,
J.
Hwang
,
G.
Bae
,
H.
Choi
, and
D. J.
Kang
,
ACS Appl. Mater. Interfaces
9
,
13571
(
2017
).
19.
S. M.
Bohaichuk
,
M. M.
Pelella
,
Y.
Sun
,
Z.
Zhang
,
S.
Ramanathan
, and
E.
Pop
,
IEEE Electron Device Lett.
41
,
292
(
2020
).
20.
K. R.
Williams
,
K.
Gupta
, and
M.
Wasilik
,
J. Microelectromech. Syst.
12
,
761
(
2003
).
21.
Y.
Xiong
,
Q.-Y.
Wen
,
Z.
Chen
,
W.
Tian
,
T.-L.
Wen
,
Y.-L.
Jing
,
Q.-H.
Yang
, and
H.-W.
Zhang
,
J. Phys. D Appl. Phys.
47
,
455304
(
2014
).
22.
A.
Gupta
,
R.
Aggarwal
,
P.
Gupta
,
T.
Dutta
,
R. J.
Narayan
, and
J.
Narayan
,
Appl. Phys. Lett.
95
,
111915
(
2009
).
23.
B.
Viswanath
,
C.
Ko
,
Z.
Yang
, and
S.
Ramanathan
,
J. Appl. Phys.
109
,
063512
(
2011
).
24.
D.
Zhang
,
T.
Wen
,
Y.
Xiong
,
D.
Qiu
, and
Q.
Wen
,
Nano Micro Lett.
9
,
29
(
2017
).
25.
S. J.
Yun
,
J. W.
Lim
,
J.-S.
Noh
,
B.-G.
Chae
, and
H.-T.
Kim
,
Jpn. J. Appl. Phys.
47
,
3067
(
2008
).
26.
S.
Jiang
,
G.
He
,
M.
Liu
,
L.
Zhu
,
S.
Liang
,
W.
Li
,
Z.
Sun
, and
M.
Tian
,
Adv. Electr. Mater.
4
,
1700543
(
2018
).
27.
A.
Hiraiwa
,
T.
Sasaki
,
S.
Okubo
,
K.
Horikawa
, and
H.
Kawarada
,
J. Appl. Phys.
123
,
155303
(
2018
).
28.
A.
Matsuda
,
Y.
Nozawa
,
S.
Kaneko
, and
M.
Yoshimoto
,
Appl. Surf. Sci.
480
,
956
(
2019
).
29.
K.
Okimura
and
Y.
Suzuki
,
Jpn. J. Appl. Phys.
50
,
065803
(
2011
).
30.
A. M.
Makarevich
,
I. I.
Sadykov
,
D. I.
Sharovarov
,
V. A.
Amelichev
,
A. A.
Adamenkov
,
D. M.
Tsymbarenko
,
A. V.
Plokhih
,
M. N.
Esaulkov
,
P. M.
Solyankin
, and
A. R.
Kaul
,
J. Mater. Chem. C
3
,
9197
(
2015
).
31.
R.
Wang
,
H.
Wang
,
Z.
An
,
J.
He
,
C.
Zhang
,
G.
Pan
, and
X.
Li
,
J. Appl. Phys.
125
,
163104
(
2019
).
32.
V.
Sanphuang
,
N.
Ghalichechian
,
N. K.
Nahar
, and
J. L.
Volakis
,
IEEE Trans. Terahertz Sci. Technol.
6
,
583
(
2016
).
33.
A.
Gupta
,
J.
Narayan
, and
T.
Dutta
,
Appl. Phys. Lett.
97
,
151912
(
2010
).
34.
J.
Nag
,
E. A.
Payzant
,
K. L.
More
, and
R. F.
Haglund
,
Appl. Phys. Lett.
98
,
251916
(
2011
).
35.
J.
Narayan
and
V. M.
Bhosle
,
J. Appl. Phys.
100
,
103524
(
2006
).
36.
Y.
Zhou
and
S.
Ramanathan
,
J. Appl. Phys.
112
,
074114
(
2012
).
37.
L. L.
Fan
,
Y. F.
Wu
,
C.
Si
,
C. W.
Zou
,
Z. M.
Qi
,
L. B.
Li
,
G. Q.
Pan
, and
Z. Y.
Wu
,
Thin Solid Films
520
,
6124
(
2012
).
38.
L. L.
Fan
,
Y. F.
Wu
,
C.
Si
,
G. Q.
Pan
,
C. W.
Zou
, and
Z. Y.
Wu
,
Appl. Phys. Lett.
102
,
011604
(
2013
).
39.
Y. X.
Guo
,
Y. F.
Liu
,
C. W.
Zou
,
Z. M.
Qi
,
Y. Y.
Wang
,
Y. Q.
Xu
,
X. L.
Wang
,
F.
Zhang
, and
R.
Zhou
,
Appl. Phys. A
115
,
1245
(
2014
).
40.
P.
Mandal
,
A.
Speck
,
C.
Ko
, and
S.
Ramanathan
,
Opt. Lett.
36
,
1927
(
2011
).
41.
Z. P.
Wu
,
S.
Yamamoto
,
A.
Miyashita
,
Z. J.
Zhang
,
K.
Narum
, and
H.
Naramoto
,
J. Phys. Condens. Matter
10
,
L765
(
1998
).
42.
R. D.
Shannon
and
R. C.
Rossi
, “The Definition of Topotaxy,” eScholarship, Lawrence Berkeley National Laboratory, Recent Work (
1963
).

Supplementary Material