β-phase gallium oxide ( β-Ga2O3) has drawn significant attention due to its large critical electric field strength and the availability of low-cost high-quality melt-grown substrates. Both aspects are advantages over gallium nitride (GaN) and silicon carbide (SiC) based power switching devices. However, because of the poor thermal conductivity of β-Ga2O3, device-level thermal management is critical to avoid performance degradation and component failure due to overheating. In addition, for high-frequency operation, the low thermal diffusivity of β-Ga2O3 results in a long thermal time constant, which hinders the use of previously developed thermal solutions for devices based on relatively high thermal conductivity materials (e.g., GaN transistors). This work investigates a double-side diamond-cooled β-Ga2O3 device architecture and provides guidelines to maximize the device’s thermal performance under both direct current (dc) and high-frequency switching operation. Under high-frequency operation, the use of a β-Ga2O3 composite substrate (bottom-side cooling) must be augmented by a diamond passivation overlayer (top-side cooling) because of the low thermal diffusivity of β-Ga2O3.
more »
« less
High throughput substrate screening for interfacial thermal management of β-Ga2O3 by deep convolutional neural network
Electronic devices get smaller and smaller in every generation. In micro-/nano-electronic devices such as high electron mobility transistors, heat dissipation has become a crucial design consideration due to the ultrahigh heat flux that has a negative effect on devices' performance and their lifetime. Therefore, thermal transport performance enhancement is required to adapt to the device size reduction. β-Ga2O3 has recently gained significant scientific interest for future power devices because of its inherent material properties such as extremely wide bandgap, outstanding Baliga's figure of merit, large critical electric field, etc. This work aims to use a machine learning approach to search promising substrates or heat sinks for cooling β-Ga2O3, in terms of high interfacial thermal conductance (ITC), from large-scale potential structures taken from existing material databases. With the ITC dataset of 1633 various substrates for β-Ga2O3 calculated by full density functional theory, we trained our recently developed convolutional neural network (CNN) model that utilizes the fused orbital field matrix (OFM) and composition descriptors. Our model proved to be superior in performance to traditional machine learning algorithms such as random forest and gradient boosting. We then deployed the CNN model to predict the ITC of 32 716 structures in contact with β-Ga2O3. The CNN model predicted the top 20 cubic and noncubic substrates with ITC on the same level as density functional theory (DFT) results on β-Ga2O3/YN and β-Ga2O3/MgO interfaces, which has the highest ITC of 1224 and 1211 MW/m2K, respectively, among the DFT-ITC datasets. Phonon density of states, group velocity, and scattering effect on high heat flux transport and consequently increased ITC are also investigated. Moderate to high phonon density of states overlap, high group velocity, and low phonon scattering are required to achieve high ITC. We also found three Magpie descriptors with strong Pearson correlation with ITC, namely, mean atomic number, mean atomic weight, and mean ground state volume per atom. Calculations of such descriptors are computationally efficient, and therefore, these descriptors provide a new route for quickly screening potential substrates from large-scale material pools for high-performance interfacial thermal management of high-electron mobility transistor devices.
more »
« less
- PAR ID:
- 10541004
- Publisher / Repository:
- APS
- Date Published:
- Journal Name:
- Journal of Applied Physics
- Volume:
- 135
- Issue:
- 20
- ISSN:
- 0021-8979
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
β-Ga2O3 is an emerging ultra-wide bandgap semiconductor, holding a tremendous potential for power-switching devices for next-generation high power electronics. The performance of such devices strongly relies on the precise control of electrical properties of β-Ga2O3, which can be achieved by implantation of dopant ions. However, a detailed understanding of the impact of ion implantation on the structure of β-Ga2O3 remains elusive. Here, using aberration-corrected scanning transmission electron microscopy, we investigate the nature of structural damage in ion-implanted β-Ga2O3 and its recovery upon heat treatment with the atomic-scale spatial resolution. We reveal that upon Sn ion implantation, Ga2O3 films undergo a phase transformation from the monoclinic β-phase to the defective cubic spinel γ-phase, which contains high-density antiphase boundaries. Using the planar defect models proposed for the γ-Al2O3, which has the same space group as β-Ga2O3, and atomic-resolution microscopy images, we identify that the observed antiphase boundaries are the {100}1/4 ⟨110⟩ type in cubic structure. We show that post-implantation annealing at 1100 °C under the N2 atmosphere effectively recovers the β-phase; however, nano-sized voids retained within the β-phase structure and a γ-phase surface layer are identified as remanent damage. Our results offer an atomic-scale insight into the structural evolution of β-Ga2O3 under ion implantation and high-temperature annealing, which is key to the optimization of semiconductor processing conditions for relevant device design and the theoretical understanding of defect formation and phase stability.more » « less
-
Teherani, Ferechteh H.; Rogers, David J. (Ed.)We demonstrated a metal-organic chemical vapor deposition (MOCVD) of smooth, thick, and monoclinic phase-pure gallium oxide (Ga2O3) on c-plane sapphire using silicon-oxygen bonding (SiOx) as a phase stabilizer. The corundum (α), monoclinic (β), and orthorhombic (ε) phases of Ga2O3 with a bandgap in the 4.4 – 5.1 eV range, are promising materials for power semiconductor devices and deep ultraviolet (UV) solar-blind photodetectors. The MOCVD systems are extensively used for homoepitaxial growth of β-Ga2O3 on (001), (100), (010), and (¯2 01) β-Ga2O3 substrates. These substrates are rare/expensive and have very low thermal conductivity; thus, are not suitable for high-power semiconductor devices. The c-plane sapphire is typically used as a substrate for high-power devices. The β-Ga2O3 grows in the (¯2 01) direction on sapphire. In this direction, the presence of high-density oxygen dangling bonds, frequent stacking faults, twinning, and other phases and planes impede the heteroepitaxy of thick β-Ga2O3. Previously phase stabilizations with SiOx have been reported for tetragonal and monoclinic hafnia. We were able to grow ~580nm thick β-Ga2O3 on sapphire by MOCVD at 750 oC through phase stabilization using silane. The samples grown with silane have a reduction in the surface roughness and resistivity from 10.7 nm to 4.4 nm and from 371.75 Ω.cm to 135.64 Ω.cm, respectively. These samples show a pure-monoclinic phase determined by x-ray diffraction (XRD); have tensile strain determined by Raman strain mapping. These results show that a thick, phase-pure -Ga2O3 can be grown on c-plane sapphire which can be suitable for creating power devices with better thermal management.more » « less
-
Abstract Among various models for estimating interfacial thermal conductance (ITC) across different material interfaces, the diffuse mismatch model (DMM) has been generally evaluated as a reliable approach for material interfaces at high temperatures. The previous works by DMM have indicated the correct order of magnitude of ITC in isotropic material interfaces. However, it cannot accurately reproduce the ITC for low-dimensional anisotropic layered materials that are desired for many potential applications. Also, the inappropriate mode matching process approximation of the phonon dispersion curve tends to overestimate the ITC. In this paper, we revisited and updated the numerical method in our previous work that utilizes a mode-to-mode comparison within the DMM framework to predict ITC with the first-principles accuracy. We employed this model to calculate ITCs between layered materials such as MoS2 and graphite and metals such as Al, Au, and Cr. We then compared our values with previous literature data from calculations of phonon dispersion curve and experimental data from time-domain thermoreflectance measurements. With a better mode matching algorithm, the updated numerical method can predict the ITCs with improved accuracy. Further analysis also confirmed that counting only the three acoustic modes and neglecting the low-frequency optical modes lead to significant underestimation of the ITC using DMM.more » « less
-
Lattice thermal conductivity (κL) is a crucial characteristic of crystalline solids with significant implications for thermal management, energy conversion, and thermal barrier coating. The advancement of computational tools based on density functional theory (DFT) has enabled the effective utilization of phonon quasi-particle-based approaches to unravel the underlying physics of various crystalline systems. While the higher order of anharmonicity is commonly used for explaining extraordinary heat transfer behaviors in crystals, the impact of exchange-correlation (XC) functionals in DFT on describing anharmonicity has been largely overlooked. The XC functional is essential for determining the accuracy of DFT in describing interactions among electrons/ions in solids and molecules. However, most XC functionals in solid-state physics are primarily focused on computing the properties that only require small atomic displacements from the equilibrium (within the harmonic approximation), such as harmonic phonons and elastic constants, while anharmonicity involves larger atomic displacements. Therefore, it is more challenging for XC functionals to accurately describe atomic interactions at the anharmonicity level. In this study, we systematically investigate the room-temperature κL of 16 binary compounds with rocksalt and zincblende structures using var- ious XC functionals such as local density approximation (LDA), Perdew-Burke-Ernzerhof (PBE), revised PBE for solid and surface (PBEsol), optimized B86b functional (optB86b), revised Tao-Perdew-Staroverov-Scuseria (revTPSS), strongly constrained and appropriately normed functional (SCAN), regularized SCAN (rSCAN) and regularized-restored SCAN (r2SCAN) in combination with different perturbation orders, including phonon within harmonic approximation (HA) plus three- phonon scattering (HA+3ph), phonon calculated using self-consistent phonon theory (SCPH) plus three-phonon scattering (SCPH+3ph), and SCPH phonon plus three- and four-phonon scattering (SCPH+3,4ph). Our results show that the XC functional exhibits strong entanglement with perturbation order and the mean relative absolute error (MRAE) of the computed κL is strongly influenced by both the XC functional and perturbation order, leading to error cancellation or amplification. The minimal (maximal) MRAE is achieved with revTPSS (rSCAN) at the HA+3ph level, SCAN (r2SCAN) at the SCPH+3ph level, and PBEsol (rSCAN) at the SCPH+3,4ph level. Among these functionals, PBEsol exhibits the highest accuracy at the highest perturbation order. The SCAN- related functionals demonstrate moderate accuracy but are suffer from numerical instability and high computational costs. Furthermore, the different impacts of quartic anharmonicity on κL in rocksalt and zincblende structures are identified by all XC functionals, attributed to the distinct lattice anharmonicity in these two structures. These findings serve as a valuable reference for selecting appropriate functionals for describing anharmonic phonons and offer insights into high-order force constant calculations that could facilitate the development of more accurate XC functionals for solid materials.more » « less
An official website of the United States government

