Next Article in Journal
Responsive Liquid Metal Droplets: From Bulk to Nano
Next Article in Special Issue
Modeling the Radial Growth of Self-Catalyzed III-V Nanowires
Previous Article in Journal
A Sustainable Strategy for Solid-Phase Extraction of Antiviral Drug from Environmental Waters by Immobilized Hydrogen Bond Acceptor
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Thermal Conductivity of GaAs Nanowire Arrays Measured by the 3ω Method

Department of Engineering Physics, McMaster University, Hamilton, Ontario, ON L8S 4L8, Canada
*
Author to whom correspondence should be addressed.
Nanomaterials 2022, 12(8), 1288; https://doi.org/10.3390/nano12081288
Submission received: 24 March 2022 / Revised: 4 April 2022 / Accepted: 7 April 2022 / Published: 10 April 2022

Abstract

:
Vertical nanowire (NW) arrays are the basis for a variety of nanoscale devices. Understanding heat transport in these devices is an important concern, especially for prospective thermoelectric applications. To facilitate thermal conductivity measurements on as-grown NW arrays, a common NW-composite device architecture was adapted for use with the 3ω method. We describe the application of this technique to obtain thermal conductivity measurements on two GaAs NW arrays featuring ~130 nm diameter NWs with a twinning superlattice (TSL) and a polytypic (zincblende/wurtzite) crystal structure, respectively. Our results indicate NW thermal conductivities of 5.2 ± 1.0 W/m-K and 8.4 ± 1.6 W/m-K in the two samples, respectively, showing a significant reduction in the former, which is the first such measurements on TSL NWs. Nearly an order of magnitude difference from the bulk thermal conductivity (~50 W/m-K) is observed for the TSL NW sample, one of the lowest values measured to date for GaAs NWs.

Graphical Abstract

1. Introduction

Interest in semiconductor nanowires (NWs) is motivated by finite size and surface effects that dominate the mechanisms of charge and phonon transport. NWs are particularly attractive as thermoelectric materials, [1,2] where efficiency is related to the dimensionless figure of merit [3].
Z T = S 2 σ T κ
Here, T is the absolute temperature, S is the Seebeck coefficient (typically a few mV/K), σ is electrical conductivity , and κ is the total thermal conductivity:
κ = κ e + κ L
The electronic component, κ e , is proportional to σ and is often much smaller than the lattice component, κ L , in semiconductors [4].
NWs stand to benefit thermoelectrics because an enhancement of the power factor ( S 2 σ ) could be achieved in low-dimensional systems, as predicted by well-known theoretical studies [5,6]; however, these confined-carrier effects have not contributed significantly to thermoelectric improvements in NWs to date [7,8,9]. Another benefit is the lower thermal conductivity of NWs compared to the bulk, which has proven more practical for increasing Z T [10]. Among the III–V compounds, measurements on InAs [11,12,13,14] and GaAs [15] NWs have shown 30–80% reductions in thermal conductivity compared to their bulk value. The influence of structural and compositional features has also been demonstrated. For example, measurements on Si [7] and Si0.96Ge0.04 [16] NWs have revealed that surface roughness can grant an additional 70% reduction compared to smooth NWs. Core-shell GaAs-AlAs NWs measured by Juntunen et al. [17] showed a non-monotonic dependence on the shell thickness, and a minimal thermal conductivity near 1 W/m-K (versus ~50 W/m-K of bulk GaAs [18,19]). Li et al. [20] measured the thermal conductivity of Si/Si0.9Ge0.1 superlattice NWs, determining values of 5 to 6 W/m-K at room temperature, much lower than Si NWs of a similar size [21].
In addition to reduced thermal conductivity, good electrical conductivity (as well as a large Seebeck coefficient) is still required for a large Z T . Interfacial and surface roughness in NWs can accomplish the former, as numerous studies have demonstrated, but the need persists for more phonon-specific mechanisms to avoid additional electron scattering. As discussed below, NWs with sharp and periodic crystallographic interfaces [22], such as twinning superlattices (TSLs), could be ideal for this purpose. Emerging growth techniques are recently enabling the controlled synthesis of TSL NWs in patterned arrays via vapor–liquid–solid (VLS) methods [23].
In a handful of III–V compounds including GaAs [24], both the zincblende (ZB) and wurtzite (WZ) phases are accessible during growth and can be selected by adjusting the growth conditions [25,26]. Polytypic (ZB/WZ) [27], phase-modulated [28], and twinning-ZB [24,26,29,30] III–V NWs can be produced in this way. Among these, disordered NW structures (Figure 1a) are associated with low thermal conductivity [27] but also low electron mobility [31], which is undesirable for devices. More ordered structures are also possible, such as twinning superlattices (TSLs; Figure 1b–d), where complementary twin segments form repeating sections that periodically rotate by 60° about the NW axis. For thermoelectric device applications, TSLs may provide a means of coherent phonon engineering [32] based on adjusting the twin period for minimal lattice thermal conductivity, as suggested by computational studies [33,34,35] and existing work on heterojunction superlattices [36]. Due to the novelty of TSL NWs and corresponding challenges in their synthesis, experimental measurements of the thermal conductivity are absent from the literature. This work provides, to the best of our knowledge, the first of such measurements on TSL NWs within the III–V material system and perhaps on TSL NWs in general.
Many of the NW thermoelectric devices proposed to date [37,38,39,40,41,42,43] employ a composite architecture, with an interstitial material introduced for mechanical support and planarization of the NW array. Among the thermal conductivity measurements reported above, the thermoreflectance-based approaches of Persson et al. [13] and Juntunen et al. [17] are applied directly to measurement devices of this type. In this case, heat flow remains highly one-dimensional, so a linear effective-medium model can be used to extract the NW thermal conductivity [13].
As an alternative to thermoreflectance measurements, we employ the AC 3ω method [44,45], which has been used previously to measure the thermal conductivity of NWs embedded in nanoporous Al2O3 [46,47,48]. In comparison to embedded (electrodeposited) NWs, free-standing VLS-grown NWs offer greater flexibility in controlling the composition and crystal structure, as mentioned earlier, with the added benefit of direct integration on Si substrates. Using a spin-on polymer for planarization, as in Refs. [13,17] and various other NW devices [49], allows for the 3ω method to be applied as a generic means of measuring the array thermal conductivity. Compared to traditional techniques, the 3ω method requires far shorter equilibration times and is insensitive to radiative losses [44]. Compared to optical techniques [13,17], the thermal penetration depth [45] is much greater in the 3ω method, because typical measurement frequencies are on the order of kHz (as opposed to MHz).
To examine the influence of crystal structure on the NW thermal conductivity, and to illustrate a novel adaptation of the 3ω method, we report measurements on GaAs NW arrays featuring polytypic ZB/WZ NWs (Figure 1a) versus twinning superlattice (TSL) NWs (Figure 1b–d). In the sections that follow, we describe the VLS growth of these NWs in dense and free-standing arrays, followed by a characterization of the crystal structure, and fabrication of the NW-composite measurement device. Next, an overview of the 3ω method is provided, with a discussion of theoretical models and the practical details of implementation. Finally, results are discussed in context with available experimental and theoretical data. Abbreviations and symbols are tabulated in Table A6 of Appendix D.

2. Materials and Methods

2.1. Nanowire Growth and Characterization

Arrays of GaAs NWs were grown on 300 μm thick p+-Si substrates ( ρ 0.005 Ω-cm) with a <111> surface orientation. GaAs NWs were grown in a 2 × 2 mm2 area on the substrate surface by the self-assisted (SA) VLS method with a Ga droplet as the seed particle, using gas source molecular beam epitaxy. NW growth details are provided in Appendix A. Two samples were grown with polytypic ZB/WZ NWs (sample A) and twinning superlattice (TSL) NWs (sample B) using identical processes, apart from a Be dopant flux introduced in the latter that induces a TSL structure due to changes in the NW sidewall surface energy [26].
After growth, the NW arrays were characterized by scanning electron microscopy (SEM), bright-field transmission electron microscopy (TEM), and high-resolution TEM (HRTEM). Samples were prepared for TEM by mechanically transferring the NWs to a Cu grid. TEM was performed using a JEOL 2010F with 200 kV accelerating voltage. Selective-area electron diffraction (SAED) was performed in the TEM to confirm crystal structure. SEM was performed using a JEOL 7000F. SEM confirmed a dense and uniform NW array (Figure 2a) from sample B with similar results to sample A with comparable NW length and diameter.
The growth procedure produced a TSL structure near the top of the NW in sample B, observed by the surface faceting in SEM (Figure 2b). A side-by-side comparison, shown in Figure 2c, highlights the structural differences between samples A and B. The crystal structure of sample A is ZB twins with WZ insertions (denoted ZB/WZ, as depicted in Figure 1a) as confirmed by HRTEM in Figure 2d. The crystal structure of sample B was identical to sample A, except near the top third of the NW where the polytypic structure was replaced with a TSL structure (depicted in Figure 1b,c), confirmed by HRTEM in Figure 2e and SAED in Figure 2f. The Ga droplet that seeded the NW growth is observable at the top of the NW in Figure 2b,c.

2.2. Device Fabrication

The 3ω measurement device is depicted in Figure 3. To fabricate the device, layers of benzocyclobutene (BCB) were applied to the NW arrays by repeated spin-coating at 7000 rpm for 77 s, followed by a 1-h cure at 250 °C in an inert N2 atmosphere after each coat. In this way, the NWs were completely submerged and insulated by an additional 3.5 μm of BCB above. A layer of microposit S1827 photoresist was then spin-coated on top of the cured BCB at 3500 rpm for 30 s, and a 1 mm2 opening was developed in the photoresist, above a suitable area of the NW array. Using the photoresist as an etch mask, a 0.5 μm deep cavity was etched into the excess BCB using reactive ion etching under 50 W power, with 35.8 sccm CF4, 5.4 sccm O2, and 1.8 sccm N2. This step removed some of the BCB above the measurement section of the NW array while allowing shorts to be avoided between the heater and longer parasitic NWs. The etch mask was then dissolved in an acetone bath.
The height of the remaining BCB excess was later measured by cross-sectional SEM, as seen in Figure 3b. In the final fabrication step, Cr (70 Å) and Pt (1500 Å) metal were deposited on the BCB by evaporation through a mask bearing the heater pattern, and the heater line was positioned diagonally across the cavity etched above the NW array. After deposition, another acetone bath was used to remove the resist and excess metal by lift-off, forming the 30 × 1000 μm heater line and four contact pads illustrated in the inset of Figure 3a.

2.3. 3ω Measurements

A Stanford Research Systems SR810 digital lock-in amplifier was used for AC measurements across a range of frequencies. The complex temperature rise ( θ ˜ 2 ω , rms ) of the heater line was then calculated from the relation [50]:
θ ˜ 2 ω , rms = [ I 1 ω , rms 2 × d R d T ] 1 ( V 3 ω , rms , re i V 3 ω , rms , im )
where I 1 ω , rms is the sinusoidal current through the heater line at the source frequency, and V 3 ω , rms , re and V 3 ω , rms , im are the in- and out-of-phase components of the third harmonic voltage across the heater line. In Equation (3), the resistance coefficient ( d R / d T ) enables the coupling of voltage and temperature via the temperature dependence of the line resistance. Resistance coefficients were obtained using a Thorlabs TED4015 temperature controller together with a Peltier element to sweep sample temperature from 15 to 30 °C, while a Keithley 2000 series multimeter was used to measure line resistances (see Figures S1–S3 in the Supplementary File, Section I). All AC measurements were done with the sample in a vacuum. The heater was contacted by copper probes in a typical four-point configuration with the substrate pressed against an aluminum block. Prior to the 3ω measurement, the electrical resistance was measured between the heater line and the aluminum block (i.e., ground) to confirm the lack of a leakage current through the sample. These measurements yielded resistances exceeding 40 MΩ, owing to the highly insulating BCB. (A small resistance would indicate significant power dissipation inside the sample, invalidating the assumed heat flow model [45].) In addition to samples A and B, a third BCB-on-Si sample (sample C) was measured to obtain the baseline thermal conductivity of the BCB. A 400 μm thick Si substrate ( ρ > 1000 Ω-cm) was used for sample C.
The standard 3 ω measurement circuit (as in Figure 2b in Ref. [51]) was used to obtain the V 3 ω , rms components and the current I 1 ω , rms = V 1 ω , rms / R ref , which was produced by a constant-current source and measured via the voltage across a precision resistor ( R ref ). Using a single lock-in amplifier, the 1 ω current measurements were conducted immediately after the 3 ω voltage measurements, over the same set of frequency points. Connections were switched externally (at the lock-in amplifier) and the contacted sample was not disturbed in any way. All measurements were conducted using the digital interface of the lock-in amplifier. We took the average reading over 10 s of continuous measurement at each frequency. Uncertainty was estimated from the standard deviation of these values, then propagated to θ ˜ 2 ω , rms via Equation (3), whereby the error on d R / d T (less than 2 %; see Supplementary File, Section I) was also incorporated.

2.4. 3ω Data Fitting

To fit the measured temperature data (from Equation (3)), the two-dimensional heat equation can be solved analytically, assuming a uniform heat flux between the line heater and the top layer [45,51]. This gives the complex temperature rise of the heater as a function of the angular frequency, ω:
T ˜ h ( ω ) = Δ T ˜ + p r th 1 + i 2 ω C h d h ( Δ T ˜ + p r th p )
where
Δ T ˜ ( ω ) = P π l κ , 1 0 1 A 1 B 1 sin 2 ( b λ ) b 2 λ 2   d λ
Here, P is the peak electrical power, while p = P / 2 b l is the heat flux, b is half the width of the heater line, and l is the length of the heater line. The integration variable, λ , is an inverse length. The parameters C h , d h , and r th represent the volumetric heat capacity, thickness, and thermal contact resistance, respectively, of the heater. In Equation (5), κ , 1 is the cross-plane thermal conductivity of the top layer (i.e., the insulating BCB). The remaining thermophysical properties are contained in the coefficients A 1 and B 1 , which are defined by
A n 1 = A n B n κ , n B n 1 κ , n 1 tanh ( B n 1 d n 1 ) 1 A n B n κ , n B n 1 κ , n 1 tanh ( B n 1 d n 1 )  
and
B n = ψ n λ 2 + i 2 ω C n κ , n  
for layer indices n = 2 , …, N , numbered in increasing order from the second layer down to the substrate (layer N ), as in Figure 3. The quantity ψ n is the anisotropy ratio, defined as ψ n = κ , n / κ , n . The coefficient A 1 of the uppermost layer is determined recursively from Equation (6) and the recurrence is terminated at the substrate layer, where A N 1 . The coefficients B n , on the other hand, are calculated directly form Equation (7). In this way, the base temperature rise in Equation (5) considers the accumulated influence of layers in the sample, while Equation (4) includes a correction accounting for the physical heater line [45]. Thermophysical properties of the sample are extracted by fitting this model (Equations (4) and (5)) to the measured temperature data (Equation (3)).
The data were fitted by minimizing the mean-square-error (MSE) defined by
ϵ ¯ ( χ ) =   1 M k = 1 M θ ˜ 2 ω k T ˜ h ( ω k ,   χ ) 2
where the vector χ contains the fitting parameters. Each layer in a sample admits four individual parameters, namely the cross-plane thermal conductivity κ , n , the volumetric heat capacity C n , the layer thickness d n , and the anisotropy ratio ψ n .
Samples A and B were modelled as three-layer structures ( N = 3 ), consisting of (1) the insulating BCB in contact with the heater, (2) the NW-BCB composite, and (3) the Si substrate, as shown in Figure 3. Sample C was modelled as a two-layer structure ( N = 2 ), consisting only of (1) a uniform BCB layer and (2) a Si substrate.
The final thermal conductivity of the NW array ( κ NW ) was calculated from the total thermal conductivity of the NW-BCB composite layer using an effective-medium model [13]:
κ NW BCB = x κ NW + ( 1 x ) κ BCB  
Here, the variable 0 x 1 represents the volume fraction of NWs in layer (2). The volume fraction must account for the NW growth yield because a fraction of oxide holes will not nucleate a NW due to parasitic effects. The yield was estimated by counting all NWs (roughly 4000 individuals) in a 30 × 30 μm2 area of each NW array versus the known density of nucleation sites. The measured thermal conductivities, κ NW BCB and κ BCB , were obtained by fitting layer (2) of samples A or B and layer (1) of sample C, respectively.

3. Results and Discussion

The 3ω measurement data is shown in Figure 4a for samples A (red), B (blue), and C (green) along with curves corresponding to the fitted material properties. We considered three separate paradigms for the fit model vis a vis Equations (4) and (5); fitting with (i) the heater contribution neglected ( C h d h = 0 , r th = 0 ), (ii) the heater thermal mass neglected ( C h d h = 0 , r th 0 ), and (iii) all heater corrections included ( C h d h 0 , r th 0 ). Uncertainties on the fitted parameter values were estimated by taking the largest range within the measurement error, which is indicated by the shaded regions in Figure 4a. The best overall fit, in terms of the minimal MSE, was achieved using paradigm (ii). Paradigm (iii) produced comparable results for the NW and BCB thermal conductivities, agreeing with (ii) within uncertainty. As can be expected [45], the thermal mass of a d h = 157 nm heater line had only a small effect on the extrapolated values. Inclusion of the thermal resistance, on the other hand, greatly improved the fit. Figure S4 in Section II of the Supplementary File confirms the sensitivity of the heat model, Equations (4)–(7), to the thermal conductivities of the upper two layers vis a vis the finite penetration depth [45,51] of the temperature oscillations, shown in Figure S5, as compared to the thicknesses of layers (1) and (2) in our samples. The complete parameters are tabulated in Table A1, Table A2 and Table A3 in Appendix B. Measurements on bulk GaAs and InP substrates were used to validate the experimental setup, producing results in excellent agreement with the well-known bulk thermophysical properties of these materials [52,53,54], as shown in Figure A1 and Table A4 and Table A5 of Appendix C.
Figure 4. (a) Measured temperature amplitudes for samples A, B, and C, along with fitted curves, for the in-phase and out-of-phase (inset) components of each. (b) Thermal conductivity of GaAs NWs, with the polytypic ZB/WZ NWs (red, sample A) and the TSL NWs (blue, sample B) data points corresponding to measurements in panel (a). For comparison, the white markers indicate the experimental results of Juntunen et al. [17] and Soini et al. [15]. The dashed and dotted lines represent theoretical results from Mingo et al. [55] and Martin et al. [56], respectively, with Δ indicating surface roughness for the latter.
Figure 4. (a) Measured temperature amplitudes for samples A, B, and C, along with fitted curves, for the in-phase and out-of-phase (inset) components of each. (b) Thermal conductivity of GaAs NWs, with the polytypic ZB/WZ NWs (red, sample A) and the TSL NWs (blue, sample B) data points corresponding to measurements in panel (a). For comparison, the white markers indicate the experimental results of Juntunen et al. [17] and Soini et al. [15]. The dashed and dotted lines represent theoretical results from Mingo et al. [55] and Martin et al. [56], respectively, with Δ indicating surface roughness for the latter.
Nanomaterials 12 01288 g004
To consider thermal contact resistance between the sample layers, we used the slightly modified model of Olson et al. [57] to extend Equations (4)–(7). This approach effectively yielded zero resistance between layers (1)–(2) and between layers (2)–(3) in both samples A and B. A finite thermal resistance between the Cr/Pt heater line and the BCB of layer (1) was necessary to fit the high frequency data ( ω > 10 4 Hz) shown in Figure 4a. r th = 0.02 ± 0.01 cm2-K/W was obtained for both samples A and B.
The array thermal conductivities were determined from NW volume fractions x A = 0.056 ± 0.006 and x B = 0.084 ± 0.009 for samples A and B, respectively. Volume fractions were calculated from yields 0.50 ± 0.05 and 0.63 ± 0.05 , and NW diameters D A = 126 ± 5 nm and D B = 138 ± 5 nm, for samples A and B, respectively. The thermal conductivity of BCB was 0.19 ± 0.03 W/m-K obtained from sample C. The thermal conductivities for the disordered polytypic ZB/WZ NWs (sample A) and the TSL NWs (sample B) were κ A = 8.4 ± 1.6 W/m-K and κ B = 5.2 ± 1.0 W/m-K, as shown in Figure 4.
As a point of comparison, the bulk GaAs thermal conductivity of ~50 W/m-K [18,19] would give a thermal conductivity of ~3 W/m-K at the same volume fraction of samples A or B, where ~0.6 W/m-K was in fact measured for the NW-BCB layers. The NW thermal conductivities extracted from the samples are indeed much smaller than 50 W/m-K. The thermal conductivity of the ZB/WZ NWs from sample A was in the same range as prior theoretical and experimental results, as shown in Figure 4b (red data point). Here, the theoretical curves taken from Mingo et al. [55] (dashed) and Martin et al. [56] (dotted) indicate the approximate size dependence of the thermal conductivity for NWs with fully diffuse phonon-boundary scattering and for NWs with boundaries characterized by a root-mean-square roughness, Δ, respectively. Apart from roughness, the rate of phonon-boundary scattering is inversely proportional to the NW diameter and imparts most of the size dependence to the thermal conductivity of NWs at the diameters studied [55,56]. Based on these trends, the ZB/WZ NWs have a proportionally larger thermal conductivity compared to the 105 nm diameter NWs from Ref. [17]; however, the TSL NWs from sample B (blue data point in Figure 4b) exhibit a decrease in thermal conductivity within experimental uncertainty, despite even larger diameter. This result is among the lowest thermal conductivities achieved to date for GaAs NWs.
Lower values may yet be attainable by carefully adjusting the TSL period. As in traditional superlattices [36,58], there should exist an optimum twin period that minimizes the thermal conductivity of TSL NWs [33,34,35]. The mechanism of thermal conductivity reduction, while not fully understood [59], is usually attributed to both coherent and incoherent phonon transport in superlattices [60,61]. A fraction of phonon modes experience repeated reflections between interfaces, leading to resonances that inhibit axial propagation. Conversely, the same interfaces form resistive barriers for phonons that decohere at shorter length scales. In this view, the superlattice period selects the phonon modes (frequencies and wavelengths) affected by either scattering regime, resulting in a tunable thermal conductivity depending on the contribution of these modes to heat conduction in the periodic structure. Compared to traditional heterojunctions [20], the inherent sharpness of TSL interfaces in NWs should improve phonon coherence and provide greater tunability.
The TSL structure can be formed in NWs of various materials, provided compatible symmetry of the crystal lattice and energetic favorability for formation, as determined by growth conditions and atomic bonding. While GaAs NWs are easier to synthesize, TSLs have been observed in NWs of other III–V [22,29,30,32] and II–V compounds [62]. Because thermal conductivity is phenomenologically similar in these crystalline materials, similar behavior can be expected for TSLs. In particular, TSL NWs of materials with a lower bulk thermal conductivity than GaAs (for example, InAs and InSb [52]) merit future consideration.

4. Conclusions

In summary, we measured the thermal conductivity of GaAs NW arrays with polytypic ZB/WZ versus TSL structure, using a NW-BCB composite device structure that we adapted for the 3ω method. The thermal conductivity of GaAs TSL NWs was measured for the first time, and the result extracted from an effective-medium model indicated a thermal conductivity of κ B = 5.2 ± 1.0 W/m-K, versus κ A = 8.4 ± 1.6 W/m-K for the polytypic sample. A lower thermal conductivity was observed in the TSL sample, despite slightly larger diameter, marking a distinct reduction. Analysis of the 3 ω data indicated a significant thermal contact resistance between the BCB polymer used for planarization and the Pt/Cr heater line used for thermometry. Thermal contact resistances between sample layers, however, were found to be negligible.
These results pave the way for future systematic studies aimed at understanding the effect of crystal structure on heat flow in III–V NWs and outline an electrical characterization method with a precision comparable to thermoreflectance-based approaches. The tenfold reduction from the bulk thermal conductivity observed in our TSL NWs represents an encouraging result of NW thermoelectrics. Further reductions in the thermal conductivity may be possible by tuning the period of the TSL structure (the twin segment length), varying the NW diameter, and inducing the TSL structure in other NW materials.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/nano12081288/s1, Figures S1–S3: temperature dependence of the line resistances of sample A, B, and C; Figure S4: sensitivity of the heat model to various sample parameters; Figure S5: thermal penetration depth in BCB and NW-BCB layers.

Author Contributions

Conceptualization, A.G. and P.O.; methodology, A.G. and P.O.; software, A.G. and P.O.; validation, R.L.; formal analysis, A.G.; investigation, A.G., P.O. and N.I.G.; resources, R.L.; data curation, A.G.; writing—original draft preparation, A.G.; writing—review and editing, R.L.; visualization, A.G.; supervision, R.L.; project administration, R.L.; funding acquisition, R.L. All authors have read and agreed to the published version of the manuscript.

Funding

The authors acknowledge support of the Natural Sciences and Engineering Research Council of Canada under grant RGPIN-2018-04015.

Data Availability Statement

Data is contained within the article and Supplementary Materials.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Using plasma-enhanced chemical vapor deposition, 30 nm of SiOx was deposited on the substrate surface. To provide NW nucleation sites, arrays of holes were etched through the surface oxide in a hexagonal array pattern with a pitch of 360 nm, using electron beam lithography and reactive ion etching. The total extent of the array was confined to a 2 × 2 mm2 area on the substrate surface, as defined by the oxide template. GaAs NWs were grown by the self-assisted vapor–liquid–solid method, with a Ga droplet as the seed particle, using gas source molecular beam epitaxy. Ga was supplied from an effusion source and group V elements were supplied as dimers from a hydride gas cracker.
The polytypic zincblende/wurtzite NWs (sample A) and TSL NWs (sample B) arrays were grown using identical processes, apart from the dopant flux introduced in the latter that is used to induce a TSL structure [26]. A 250 s Ga pre-deposition and the brief growth of a short GaP base were employed to improve the NW yield [63]. GaAs NWs were then grown in three segments: (i) 50 min at a substrate temperature of 630 °C, (ii) a ramp down to 537 °C over the next 50 min, and (iii) a final 30 min of growth at 537 °C. This process ensured a high NW yield at the beginning of growth, followed by a low temperature required for TSL formation in sample B. During the growth of sample B, a Be dopant flux supplied from an elemental effusion cell was introduced in segment (i) at a NW dopant concentration of 6.0 × 10 17 cm−3, then ramped up in segment (ii), and finally held at 3.0 × 10 19 cm−3 in segment (iii). The Be dopant flux was increased gradually to avoid Ga droplet instability associated with a high dopant flux. The incorporated dopant concentration was based on previous calibrations [64,65] and no intentional dopants were introduced in sample A.

Appendix B

The thermophysical properties of all layers across the three samples, as determined by three fitting paradigms, (i)–(iii), are tabulated in Table A1, Table A2 and Table A3. Error estimates are included on the parameter values for which the model exhibits high sensitivity (see the Supplementary File, Section II). The results for the anisotropy ratio ( κ / κ ) were on the order of 10 7 for both sample A and sample B, suggesting highly one-dimensional heat flow through the NW-BCB composites. For the uniform BCB layer in sample C, we fixed ψ 1 = 1 .
Reference value ranges are provided, where available, in the rightmost columns of each table. All entries in the columns (i)–(iii) represent fitted parameter values for the corresponding fit paradigm, except for underlined values, which indicate constant parameters, and values of the MSE, ϵ ¯ . Fitted parameters exceeding their reasonable physical range are emphasized in italic text.
The inclusion of r th 0 improves the MSE by at least an order of magnitude, so paradigm (i) was discounted. Regarding paradigms (ii) and (iii), we note that thermal conductivity results for the BCB and NW-BCB layers agree within uncertainty for all three samples. Indeed, the heater thermal mass has only a small effect on the measured thermal conductivity (as per Ref. [45] (p. 2146)), considering the relative dimensions and heat capacity of the heater line versus that of the target layers. Paradigm (ii) was chosen over (iii) because it produced no unphysical values as well as the lowest overall MSE across the three samples.
Table A1. Fitted parameter values for sample A.
Table A1. Fitted parameter values for sample A.
ParameterUnits(i)
C h d h ,   r t h = 0
(ii)
r t h 0
(iii)
C h d h , r t h 0
Ref. Value
κ BCB W/m-K0.55 ± 0.030.19 ± 0.040.17 ± 0.040.18–0.29 [66,67,68,69]
C BCB J/cm3-K1.10 ± 0.031.88 ± 0.041.99 ± 0.042.19–2.29 [66,70]
κ NW BCB W/m-K0.20 ± 0.050.65 ± 0.040.70 ± 0.06
C NW BCB J/cm3-K3.70 × 10−70.414.57 × 10−7
κ p + Si W/m-K266674635–55 [71]
C p + Si J/cm3-K0.09201.220.97-
r th cm2-K/W00.01140.028-
C h J/cm3-K002.822.82 [52]
d h nm00157-
ϵ ¯ 10−3 K219.22.251.92
Table A2. Fitted parameter values for sample B.
Table A2. Fitted parameter values for sample B.
ParameterUnits(i)
C h d h ,   r t h = 0
(ii)
r t h 0
(iii)
C h d h , r t h 0
Ref. Value
κ BCB W/m-K0.47 ± 0.020.19 ± 0.040.17 ± 0.040.18–0.29 [66,67,68,69]
C BCB J/cm3-K1.18 ± 0.041.93 ± 0.042.02 ± 0.042.19–2.29 [66,70]
κ NW BCB W/m-K0.20 ± 0.050.61 ± 0.040.68 ± 0.07
C NW BCB J/cm3-K2.910.380.37
κ p + Si W/m-K8.46 × 105574835–55 [71]
C p + Si J/cm3-K4.682.022.05-
r th cm2-K/W00.0100.027-
C h J/cm3-K002.822.82 [52]
d h nm00157-
ϵ ¯ 10−3 K239.22.233.01
Table A3. Fitted parameter values for sample C.
Table A3. Fitted parameter values for sample C.
Parameter.Units(i)
C h d h ,   r t h = 0
(ii)
r t h 0
(iii)
C h d h , r t h 0
Ref. Value
κ BCB W/m-K0.20 ± 0.030.19 ± 0.030.17 ± 0.030.18–0.29 [66,67,68,69]
C BCB J/cm3-K2.22 ± 0.051.85 ± 0.052.05 ± 0.052.19–2.29 [66,70]
κ Si W/m-K2561425.03 × 104145–156 [72,73]
C Si J/cm3-K1.871.821.911.66 [74]
r th cm2-K/W00.0120.025-
C h J/cm3-K002.822.82 [52]
d h nm00157-
ϵ ¯ 10−3 K21431.3011.7

Appendix C

The thermophysical properties of the bulk GaAs and InP substrates, which were measured to validate the experimental setup, are shown in Table A4 and Table A5 in this section, respectively. Figure A1 shows the measured temperature amplitudes and corresponding fit lines for both samples.
Figure A1. Measured and fitted temperature amplitudes for bulk GaAs and InP samples. Shaded regions around the data points indicate the measurement uncertainty.
Figure A1. Measured and fitted temperature amplitudes for bulk GaAs and InP samples. Shaded regions around the data points indicate the measurement uncertainty.
Nanomaterials 12 01288 g0a1
Table A4. Fitted parameter values for bulk GaAs.
Table A4. Fitted parameter values for bulk GaAs.
ParameterUnits Ref. Value
κ GaAs W/m-K49 ± 249–56 [19,52]
C GaAs J/cm3-K1.9 ± 0.11.8 [54]
r th cm2-K/W0.00011-
C h J/cm3-K2.822.82 [52]
d h nm157-
ϵ ¯ 10−3 K20.14
Table A5. Fitted parameter values for bulk InP.
Table A5. Fitted parameter values for bulk InP.
ParameterUnits Ref. Value
κ InP W/m-K77 ± 270–80 [52,75]
C InP J/cm3-K1.5 ± 0.11.5 [76]
r th cm2-K/W0.00011-
C h J/cm3-K2.822.82 [52]
d h nm157-
ϵ ¯ 10−3 K20.15

Appendix D

Table A6. Abbreviations and symbols.
Table A6. Abbreviations and symbols.
AbbreviationsSymbolsSymbols Units
BCBbenzocyclobutene Z T thermoelectric figure of merit1
MSEmean squared error S Seebeck coefficientV/K
NWnanowire σ electrical conductivityS/m
SAself-assisted T absolute temperatureK
SAEDselective-area electron diffraction κ total thermal conductivityW/m-K
SEMscanning electron microscopy κ e electronic thermal conductivityW/m-K
TEMtransmission electron microscopy κ L lattice thermal conductivityW/m-K
TSLtwinning superlattice ρ electrical resistivityΩ-m
VLSvapor–liquid–solid θ ˜ 2 ω , rms measured temperature riseK
WZwurtzite I 1 ω , rms measured source currentA
ZBzincblende d R / d T heater line resistance coefficientΩ/K
V 3 ω , rms , re measured in-phase 3 ω voltageV
V 3 ω , rms , im measured out-of-phase 3 ω voltageV
ω frequencyHz
T ˜ h corrected heater temperature riseK
Δ T ˜ base heater temperature riseK
p heat fluxW/m2
r th heater line thermal contact resistancem2-K/W
C h heater line heat capacityJ/m3-K
d h heater line thicknessm
P peak electrical powerW
l heater line lengthm
κ , n cross-plane thermal conductivity of layer nW/m-K
C n heat capacity of layer nJ/m3-K
d n thickness of layer nm
ψ n thermal conductivity anisotropy ratio of layer n1
A n recursive coefficient Equation (5)1
B n coefficient Equation (5)m−1
b heater line half-widthm
λ integration variable Equation (5)m−1
ϵ ¯ mean squared errorK2
χ sample parameters vector-
κ NW BCB measured NW-BCB thermal conductivityW/m-K
κ NW measured NW array thermal conductivityW/m-K
κ BCB measured BCB thermal conductivityW/m-K
x A ( B ) NW volume fraction in sample A(B)1
D A ( B ) NW diameter in sample A(B)m
κ A ( B ) NW array thermal conductivity in sample A(B)W/m-K

References

  1. Zebarjadi, M.; Esfarjani, K.; Dresselhaus, M.S.; Ren, Z.F.; Chen, G. Perspectives on Thermoelectrics: From Fundamentals to Device Applications. Energy Environ. Sci. 2012, 5, 5147–5162. [Google Scholar] [CrossRef] [Green Version]
  2. Vineis, C.J.; Shakouri, A.; Majumdar, A.; Kanatzidis, M.G. Nanostructured Thermoelectrics: Big Efficiency Gains from Small Features. Adv. Mater. 2010, 22, 3970–3980. [Google Scholar] [CrossRef]
  3. Rowe, D.M. CRC Handbook of Thermoelectrics; CRC Press: Boca Raton, FL, USA, 2018; ISBN 978-0-429-95667-6. [Google Scholar]
  4. Pichanusakorn, P.; Bandaru, P. Nanostructured Thermoelectrics. Mater. Sci. Eng. R Rep. 2010, 67, 19–63. [Google Scholar] [CrossRef]
  5. Hicks, L.D.; Dresselhaus, M.S. Thermoelectric Figure of Merit of a One-Dimensional Conductor. Phys. Rev. B 1993, 47, 16631–16634. [Google Scholar] [CrossRef]
  6. Hicks, L.D.; Dresselhaus, M.S. Effect of Quantum-Well Structures on the Thermoelectric Figure of Merit. Phys. Rev. B 1993, 47, 12727–12731. [Google Scholar] [CrossRef]
  7. Hochbaum, A.I.; Chen, R.; Delgado, R.D.; Liang, W.; Garnett, E.C.; Najarian, M.; Majumdar, A.; Yang, P. Enhanced Thermoelectric Performance of Rough Silicon Nanowires. Nature 2008, 451, 163–167. [Google Scholar] [CrossRef]
  8. Seol, J.H.; Moore, A.L.; Saha, S.K.; Zhou, F.; Shi, L.; Ye, Q.L.; Scheffler, R.; Mingo, N.; Yamada, T. Measurement and Analysis of Thermopower and Electrical Conductivity of an Indium Antimonide Nanowire from a Vapor-Liquid-Solid Method. J. Appl. Phys. 2007, 101, 023706. [Google Scholar] [CrossRef]
  9. Schmidt, V.; Mensch, P.F.J.; Karg, S.F.; Gotsmann, B.; Das Kanungo, P.; Schmid, H.; Riel, H. Using the Seebeck Coefficient to Determine Charge Carrier Concentration, Mobility, and Relaxation Time in InAs Nanowires. Appl. Phys. Lett. 2014, 104, 012113. [Google Scholar] [CrossRef]
  10. Hasan, N.; Wahid, H.; Nayan, N.; Ali, M.S.M. Inorganic Thermoelectric Materials: A Review. Int. J. Energy Res. 2020, 44, 6170–6222. [Google Scholar] [CrossRef]
  11. Zhou, F.; Moore, A.L.; Bolinsson, J.; Persson, A.; Fröberg, L.; Pettes, M.T.; Kong, H.; Rabenberg, L.; Caroff, P.; Stewart, D.A.; et al. Thermal Conductivity of Indium Arsenide Nanowires with Wurtzite and Zinc Blende Phases. Phys. Rev. B 2011, 83, 205416. [Google Scholar] [CrossRef] [Green Version]
  12. Ren, P.; Zhu, X.; Han, J.; Xu, J.; Ma, L.; Li, H.; Zhuang, X.; Zhou, H.; Zhang, Q.; Xia, M.; et al. Synthesis and Diameter-Dependent Thermal Conductivity of InAs Nanowires. Nano-Micro Lett. 2014, 6, 301–306. [Google Scholar] [CrossRef] [Green Version]
  13. Persson, A.I.; Koh, Y.K.; Cahill, D.G.; Samuelson, L.; Linke, H. Thermal Conductance of InAs Nanowire Composites. Nano Lett. 2009, 9, 4484–4488. [Google Scholar] [CrossRef] [PubMed]
  14. Swinkels, M.; van Delft, M.; Oliveira, D.; Cavalli, A.; Zardo, I.; van der Heijden, R.; Bakkers, E. Diameter Dependence of the Thermal Conductivity of InAs Nanowires. Nanotechnology 2015, 26, 385401. [Google Scholar] [CrossRef] [PubMed]
  15. Soini, M.; Zardo, I.; Uccelli, E.; Funk, S.; Koblmüller, G.; Abstreiter, G. Thermal Conductivity of GaAs Nanowires Studied by Micro-Raman Spectroscopy Combined with Laser Heating. Appl. Phys. Lett. 2010, 97, 263107. [Google Scholar] [CrossRef]
  16. Kim, H.; Park, Y.-H.; Kim, I.; Kim, J.; Choi, H.-J.; Kim, W. Effect of Surface Roughness on Thermal Conductivity of VLS-Grown Rough Si1−xGex Nanowires. Appl. Phys. A 2011, 104, 23–28. [Google Scholar] [CrossRef]
  17. Juntunen, T.; Koskinen, T.; Khayrudinov, V.; Haggren, T.; Jiang, H.; Lipsanen, H.; Tittonen, I. Thermal Conductivity Suppression in GaAs–AlAs Core–Shell Nanowire Arrays. Nanoscale 2019, 11, 20507–20513. [Google Scholar] [CrossRef] [Green Version]
  18. Holland, M.G. Analysis of Lattice Thermal Conductivity. Phys. Rev. 1963, 132, 2461–2471. [Google Scholar] [CrossRef]
  19. Carlson, R.O.; Slack, G.A.; Silverman, S.J. Thermal Conductivity of GaAs and GaAs1− xPx Laser Semiconductors. J. Appl. Phys. 1965, 36, 505–507. [Google Scholar] [CrossRef]
  20. Li, D.; Wu, Y.; Fan, R.; Yang, P.; Majumdar, A. Thermal Conductivity of Si/SiGe Superlattice Nanowires. Appl. Phys. Lett. 2003, 83, 3186–3188. [Google Scholar] [CrossRef] [Green Version]
  21. Li, D.; Wu, Y.; Kim, P.; Shi, L.; Yang, P.; Majumdar, A. Thermal Conductivity of Individual Silicon Nanowires. Appl. Phys. Lett. 2003, 83, 2934–2936. [Google Scholar] [CrossRef]
  22. Caroff, P.; Dick, K.A.; Johansson, J.; Messing, M.E.; Deppert, K.; Samuelson, L. Controlled Polytypic and Twin-Plane Superlattices in III–V Nanowires. Nat. Nanotechnol. 2009, 4, 50–55. [Google Scholar] [CrossRef] [PubMed]
  23. McIntyre, P.C.; Fontcuberta i Morral, A. Semiconductor Nanowires: To Grow or Not to Grow? Mater. Today Nano 2020, 9, 100058. [Google Scholar] [CrossRef]
  24. Burgess, T.; Breuer, S.; Caroff, P.; Wong-Leung, J.; Gao, Q.; Tan, H.H.; Jagadish, C. Twinning Superlattice Formation in GaAs Nanowires. ACS Nano 2013, 7, 8105–8114. [Google Scholar] [CrossRef] [PubMed]
  25. Panciera, F.; Baraissov, Z.; Patriarche, G.; Dubrovskii, V.G.; Glas, F.; Travers, L.; Mirsaidov, U.; Harmand, J.-C. Phase Selection in Self-Catalyzed GaAs Nanowires. Nano Lett. 2020, 20, 1669–1675. [Google Scholar] [CrossRef] [PubMed]
  26. Goktas, N.I.; Sokolovskii, A.; Dubrovskii, V.G.; LaPierre, R.R. Formation Mechanism of Twinning Superlattices in Doped GaAs Nanowires. Nano Lett. 2020, 20, 3344–3351. [Google Scholar] [CrossRef] [PubMed]
  27. Weathers, A.; Moore, A.L.; Pettes, M.T.; Salta, D.; Kim, J.; Dick, K.; Samuelson, L.; Linke, H.; Caroff, P.; Shi, L. Phonon Transport and Thermoelectricity in Defect-Engineered InAs Nanowires. In MRS Proceedings; Cambridge University Press: Cambridge, UK, 2012; Volume 1404. [Google Scholar] [CrossRef]
  28. Fedorov, V.V.; Dvoretckaia, L.N.; Kirilenko, D.A.; Mukhin, I.S.; Dubrovskii, V.G. Formation of Wurtzite Sections in Self-Catalyzed GaP Nanowires by Droplet Consumption. Nanotechnology 2021, 32, 495601. [Google Scholar] [CrossRef]
  29. Algra, R.E.; Verheijen, M.A.; Borgström, M.T.; Feiner, L.-F.; Immink, G.; van Enckevort, W.J.P.; Vlieg, E.; Bakkers, E.P.A.M. Twinning Superlattices in Indium Phosphide Nanowires. Nature 2008, 456, 369–372. [Google Scholar] [CrossRef] [Green Version]
  30. Yuan, X.; Guo, Y.; Caroff, P.; He, J.; Tan, H.H.; Jagadish, C. Dopant-Free Twinning Superlattice Formation in InSb and InP Nanowires. Phys. Status Solidi RRL 2017, 11, 1700310. [Google Scholar] [CrossRef]
  31. Karg, S.F. Full Thermoelectric Characterization of InAs Nanowires Using MEMS Heater/Sensors. Nanotechnology 2014, 25, 305702. [Google Scholar] [CrossRef]
  32. De Luca, M.; Fasolato, C.; Verheijen, M.A.; Ren, Y.; Swinkels, M.Y.; Kölling, S.; Bakkers, E.P.A.M.; Rurali, R.; Cartoixà, X.; Zardo, I. Phonon Engineering in Twinning Superlattice Nanowires. Nano Lett. 2019, 19, 4702–4711. [Google Scholar] [CrossRef] [Green Version]
  33. Hu, M.; Poulikakos, D. Si/Ge Superlattice Nanowires with Ultralow Thermal Conductivity. Nano Lett. 2012, 12, 5487–5494. [Google Scholar] [CrossRef] [PubMed]
  34. Porter, A.; Tran, C.; Sansoz, F. Intrinsic Nanotwin Effect on Thermal Boundary Conductance in Bulk and Single-Nanowire Twinning Superlattices. Phys. Rev. B 2016, 93, 195431. [Google Scholar] [CrossRef] [Green Version]
  35. Xiong, S.; Kosevich, Y.A.; Sääskilahti, K.; Ni, Y.; Volz, S. Tunable Thermal Conductivity in Silicon Twinning Superlattice Nanowires. Phys. Rev. B 2014, 90, 195439. [Google Scholar] [CrossRef]
  36. Venkatasubramanian, R. Lattice Thermal Conductivity Reduction and Phonon Localizationlike Behavior in Superlattice Structures. Phys. Rev. B 2000, 61, 3091–3097. [Google Scholar] [CrossRef]
  37. Elyamny, S.; Dimaggio, E.; Magagna, S.; Narducci, D.; Pennelli, G. High Power Thermoelectric Generator Based on Vertical Silicon Nanowires. Nano Lett. 2020, 20, 4748–4753. [Google Scholar] [CrossRef] [PubMed]
  38. Xu, B.; Khouri, W.; Fobelets, K. Two-Sided Silicon Nanowire Array/Bulk Thermoelectric Power Generator. IEEE Electron Device Lett. 2014, 35, 596–598. [Google Scholar] [CrossRef]
  39. Curtin, B.M.; Fang, E.W.; Bowers, J.E. Highly Ordered Vertical Silicon Nanowire Array Composite Thin Films for Thermoelectric Devices. J. Electron. Mater. 2012, 41, 887–894. [Google Scholar] [CrossRef]
  40. Li, Y.; Buddharaju, K.; Singh, N.; Lo, G.Q.; Lee, S.J. Chip-Level Thermoelectric Power Generators Based on High-Density Silicon Nanowire Array Prepared With Top-Down CMOS Technology. IEEE Electron Device Lett. 2011, 32, 674–676. [Google Scholar] [CrossRef]
  41. Keyani, J.; Stacy, A.M.; Sharp, J. Assembly and Measurement of a Hybrid Nanowire-Bulk Thermoelectric Device. Appl. Phys. Lett. 2006, 89, 233106. [Google Scholar] [CrossRef]
  42. Wang, W.; Jia, F.; Huang, Q.; Zhang, J. A New Type of Low Power Thermoelectric Micro-Generator Fabricated by Nanowire Array Thermoelectric Material. Microelectron. Eng. 2005, 77, 223–229. [Google Scholar] [CrossRef]
  43. Abramson, A.R.; Kim, W.C.; Huxtable, S.T.; Yan, H.; Wu, Y.; Majumdar, A.; Tien, C.-L.; Yang, P. Fabrication and Characterization of a Nanowire/Polymer-Based Nanocomposite for a Prototype Thermoelectric Device. J. Microelectromech. Syst. 2004, 13, 505–513. [Google Scholar] [CrossRef]
  44. Cahill, D.G. Thermal Conductivity Measurement from 30 to 750 K: The 3ω Method. Rev. Sci. Instrum. 1990, 61, 802–808. [Google Scholar] [CrossRef]
  45. Borca-Tasciuc, T.; Kumar, A.R.; Chen, G. Data Reduction in 3ω Method for Thin-Film Thermal Conductivity Determination. Rev. Sci. Instrum. 2001, 72, 2139–2147. [Google Scholar] [CrossRef]
  46. Paterson, J.; Singhal, D.; Tainoff, D.; Richard, J.; Bourgeois, O. Thermal Conductivity and Thermal Boundary Resistance of Amorphous Al2O3 Thin Films on Germanium and Sapphire. J. Appl. Phys. 2020, 127, 245105. [Google Scholar] [CrossRef]
  47. Singhal, D.; Paterson, J.; Ben-Khedim, M.; Tainoff, D.; Cagnon, L.; Richard, J.; Chavez-Angel, E.; Fernandez, J.J.; Sotomayor-Torres, C.M.; Lacroix, D.; et al. Nanowire Forest of Pnictogen–Chalcogenide Alloys for Thermoelectricity. Nanoscale 2019, 11, 13423–13430. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Singhal, D.; Paterson, J.; Tainoff, D.; Richard, J.; Ben-Khedim, M.; Gentile, P.; Cagnon, L.; Bourgault, D.; Buttard, D.; Bourgeois, O. Measurement of Anisotropic Thermal Conductivity of a Dense Forest of Nanowires Using the 3ω Method. Rev. Sci. Instrum. 2018, 89, 084902. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Goktas, N.I.; Wilson, P.; Ghukasyan, A.; Wagner, D.; McNamee, S.; LaPierre, R.R. Nanowires for Energy: A Review. Appl. Phys. Rev. 2018, 5, 041205. [Google Scholar] [CrossRef]
  50. Dames, C. Measuring the Thermal Conductivity of Thin Films: 3 Omega and Related Electrothermal Methods. Annu. Rev. Heat Transf. 2013, 16, 43. [Google Scholar] [CrossRef]
  51. Dames, C.; Chen, G. 1ω, 2ω, and 3ω Methods for Measurements of Thermal Properties. Rev. Sci. Instrum. 2005, 76, 124902. [Google Scholar] [CrossRef]
  52. Cardarelli, F. Materials Handbook: A Concise Desktop Reference; Springer International Publishing: Cham, Switzerland, 2018; ISBN 978-3-319-38923-3. [Google Scholar]
  53. Lide, D.R.; Baysinger, G.; Chemistry, S.; Berger, L.I.; Goldberg, R.N.; Kehiaian, H.V. CRC Handbook of Chemistry and Physics; CRC Press: Boca Raton, FL, USA, 2004; p. 2661. ISBN 0-8493-0485-7. [Google Scholar]
  54. Blakemore, J.S. Semiconducting and Other Major Properties of Gallium Arsenide. J. Appl. Phys. 1982, 53, R123–R181. [Google Scholar] [CrossRef]
  55. Mingo, N.; Broido, D.A. Lattice Thermal Conductivity Crossovers in Semiconductor Nanowires. Phys. Rev. Lett. 2004, 93, 246106. [Google Scholar] [CrossRef] [PubMed]
  56. Martin, P.N.; Aksamija, Z.; Pop, E.; Ravaioli, U. Reduced Thermal Conductivity in Nanoengineered Rough Ge and GaAs Nanowires. Nano Lett. 2010, 10, 1120–1124. [Google Scholar] [CrossRef] [PubMed]
  57. Olson, B.W.; Graham, S.; Chen, K. A Practical Extension of the 3ω Method to Multilayer Structures. Rev. Sci. Instrum. 2005, 76, 053901. [Google Scholar] [CrossRef]
  58. Venkatasubramanian, R.; Siivola, E.; Colpitts, T.; O’Quinn, B. Thin-Film Thermoelectric Devices with High Room-Temperature Figures of Merit. Nature 2001, 413, 597–602. [Google Scholar] [CrossRef]
  59. Sheng, P. Introduction to Wave Scattering, Localization, and Mesoscopic Phenomena, 2nd ed.; Springer Series in Materials Science; Springer: Berlin/Heidelberg, Germany; New York, NY, USA, 2006; ISBN 978-3-540-29155-8. [Google Scholar]
  60. Tian, Z.; Esfarjani, K.; Chen, G. Green’s Function Studies of Phonon Transport across Si/Ge Superlattices. Phys. Rev. B 2014, 89, 235307. [Google Scholar] [CrossRef] [Green Version]
  61. Mu, X. Ultra-Low Thermal Conductivity in Si/Ge Hierarchical Superlattice Nanowire. Sci. Rep. 2015, 5, 16697. [Google Scholar] [CrossRef] [Green Version]
  62. Choi, S.B.; Song, M.S.; Kim, Y. Zn3P2 Twinning Superlattice Nanowires Grown on Fluorine-Doped Tin Oxide Glass Substrates. J. Phys. Chem. C 2019, 123, 4597–4604. [Google Scholar] [CrossRef]
  63. Wilson, D.P.; Dubrovskii, V.G.; LaPierre, R.R. Improving the Yield of GaAs Nanowires on Silicon by Ga Pre-Deposition. Nanotechnology 2021, 32, 265301. [Google Scholar] [CrossRef] [PubMed]
  64. Goktas, N.I.; Fiordaliso, E.M.; LaPierre, R.R. Doping Assessment in GaAs Nanowires. Nanotechnology 2018, 29, 234001. [Google Scholar] [CrossRef]
  65. Dubrovskii, V.G.; Hijazi, H.; Goktas, N.I.; LaPierre, R.R. Be, Te, and Si Doping of GaAs Nanowires: Theory and Experiment. J. Phys. Chem. C 2020, 124, 17299–17307. [Google Scholar] [CrossRef]
  66. Pinel, S.; Marty, A.; Tasselli, J.; Bailbe, J.-P.; Beyne, E.; Van Hoof, R.; Marco, S.; Morante, J.R.; Vendier, O.; Huan, M. Thermal Modeling and Management in Ultrathin Chip Stack Technology. IEEE Trans. Comp. Packag. Technol. 2002, 25, 244–253. [Google Scholar] [CrossRef]
  67. Gracias, A.; Tokranova, N.; Thelen, B.C.M.; Castracane, J. Influence of Diamond Nanoparticles on the Thermal Properties of Benzocyclobutene (BCB). Phys. Status Solidi A 2011, 208, 684–690. [Google Scholar] [CrossRef]
  68. Moghaddam, S.; Kiger, K.T.; Modafe, A.; Ghodssi, R. A Novel Benzocyclobutene-Based Device for Studying the Dynamics of Heat Transfer during the Nucleation Process. J. Microelectromech. Syst. 2007, 16, 1355–1366. [Google Scholar] [CrossRef]
  69. Yang, J. Processing Procedures for CYCLOTENE 3000 Series Dry Etch Resins; DOW Chemical Co.: Midland, MI, USA, 2005. [Google Scholar]
  70. Modafe, A.; Ghalichechian, N.; Powers, M.; Khbeis, M.; Ghodssi, R. Embedded Benzocyclobutene in Silicon: An Integrated Fabrication Process for Electrical and Thermal Isolation in MEMS. Microelectron. Eng. 2005, 82, 154–167. [Google Scholar] [CrossRef]
  71. Slack, G.A. Thermal Conductivity of Pure and Impure Silicon, Silicon Carbide, and Diamond. J. Appl. Phys. 1964, 35, 3460–3466. [Google Scholar] [CrossRef]
  72. Maycock, P.D. Thermal Conductivity of Silicon, Germanium, III-V Compounds and III-V Alloys. Solid-State Electron. 1967, 10, 161–168. [Google Scholar] [CrossRef]
  73. Glassbrenner, C.J.; Slack, G.A. Thermal Conductivity of Silicon and Germanium from 3°K to the Melting Point. Phys. Rev. 1964, 134, A1058–A1069. [Google Scholar] [CrossRef]
  74. Flubacher, P.; Leadbetter, A.J.; Morrison, J.A. The Heat Capacity of Pure Silicon and Germanium and Properties of Their Vibrational Frequency Spectra. Philos. Mag. 1959, 4, 273–294. [Google Scholar] [CrossRef]
  75. Childs, G.E.; Ericks, L.J.; Powell, R.L. Thermal Conductivity of Solids at Room Temperature and Below: A Review and Compilation of the Literature; National Bureau of Standards: Boulder, CO, USA, 1973.
  76. Piesbergen, U. Die Durchschnittlichen Atomwärmen der AIIIBv-Halbleiter AlSb, GaAs, GaSb, InP, InAs, LnSb und die Atomwärme des Elements Germanium Zwischen 12 und 273°K. Z. Naturforsch. A 1963, 18, 141–147. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Three-dimensional structures for (a) polytypic zincblende/wurtzite (ZB/WZ) and (b) twinning superlattice (TSL) nanowires, showing [111]A and [111]B surface faceting in the latter. (c) Bilayer stacking in the TSL, seen in an orthogonal projection along [ 11 2 ¯ ], exhibiting reversal of the normal ABC stacking sequence across the twin plane. (d) Twinning is equivalent to a rotation of the crystal structure about the NW axis by 60°, illustrated by tetrahedral primitives of bulk GaAs. The indicated growth direction applies to the entire figure.
Figure 1. Three-dimensional structures for (a) polytypic zincblende/wurtzite (ZB/WZ) and (b) twinning superlattice (TSL) nanowires, showing [111]A and [111]B surface faceting in the latter. (c) Bilayer stacking in the TSL, seen in an orthogonal projection along [ 11 2 ¯ ], exhibiting reversal of the normal ABC stacking sequence across the twin plane. (d) Twinning is equivalent to a rotation of the crystal structure about the NW axis by 60°, illustrated by tetrahedral primitives of bulk GaAs. The indicated growth direction applies to the entire figure.
Nanomaterials 12 01288 g001
Figure 2. (a) SEM image of the NW array from sample B; (b) SEM image of a NW top, showing surface faceting due to a TSL; (c) Side-by-side TEM images of a polytypic NW from sample A (left) and a TSL NW from sample B (right); (d) HRTEM image near the center of the same NW from sample A, showing a polytypic WZ/ZB NW structure; (e) HRTEM image near the center of the TSL in sample B; (f) Selective-area electron diffraction pattern confirming ZB twins in the TSL region of sample B.
Figure 2. (a) SEM image of the NW array from sample B; (b) SEM image of a NW top, showing surface faceting due to a TSL; (c) Side-by-side TEM images of a polytypic NW from sample A (left) and a TSL NW from sample B (right); (d) HRTEM image near the center of the same NW from sample A, showing a polytypic WZ/ZB NW structure; (e) HRTEM image near the center of the TSL in sample B; (f) Selective-area electron diffraction pattern confirming ZB twins in the TSL region of sample B.
Nanomaterials 12 01288 g002
Figure 3. Diagrams illustrating the device layers used to model heat conduction: (1) the electrically insulating BCB excess, (2) the target NW-BCB composite layer, and (3) the silicon substrate. (a) 3D schematic of the device, with an inset showing the heater line and four contact pads. (b) Cross-sectional SEM image of a device from sample B.
Figure 3. Diagrams illustrating the device layers used to model heat conduction: (1) the electrically insulating BCB excess, (2) the target NW-BCB composite layer, and (3) the silicon substrate. (a) 3D schematic of the device, with an inset showing the heater line and four contact pads. (b) Cross-sectional SEM image of a device from sample B.
Nanomaterials 12 01288 g003
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ghukasyan, A.; Oliveira, P.; Goktas, N.I.; LaPierre, R. Thermal Conductivity of GaAs Nanowire Arrays Measured by the 3ω Method. Nanomaterials 2022, 12, 1288. https://doi.org/10.3390/nano12081288

AMA Style

Ghukasyan A, Oliveira P, Goktas NI, LaPierre R. Thermal Conductivity of GaAs Nanowire Arrays Measured by the 3ω Method. Nanomaterials. 2022; 12(8):1288. https://doi.org/10.3390/nano12081288

Chicago/Turabian Style

Ghukasyan, Ara, Pedro Oliveira, Nebile Isik Goktas, and Ray LaPierre. 2022. "Thermal Conductivity of GaAs Nanowire Arrays Measured by the 3ω Method" Nanomaterials 12, no. 8: 1288. https://doi.org/10.3390/nano12081288

APA Style

Ghukasyan, A., Oliveira, P., Goktas, N. I., & LaPierre, R. (2022). Thermal Conductivity of GaAs Nanowire Arrays Measured by the 3ω Method. Nanomaterials, 12(8), 1288. https://doi.org/10.3390/nano12081288

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop