Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Cross Flow”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 577 records · Page 32

Viewing is understanding: Graphite microstructure effects on infiltrated molten salt distribution revealed by 3D neutron tomography

Molten salt infiltration in the pore network of nuclear graphite may cause unwanted changes to graphite's local structure and mechanical and thermal properties. A detailed and comprehensive understanding of molten salt intrusion (distribution across sample cross section and penetration depth) is needed to assess its effects. Here, in this work, we report on an improved methodology for the use of neutron imaging (computed tomography) to evaluate salt penetration and distribution of a wide range of graphite grades with diverse microstructures. Neutron tomography data were acquired on the same graphite sample before and after salt intrusion; the 3D reconstructed volumes were digitally co-registered and subtracted. The difference in neutron attenuation coefficient represents direct visualization of FLiNaK (LiF–NaF–KF) salt distribution in the salt-impregnated graphite samples. This improved methodology was applied to investigate the effect of exposure times (12 h and 336 h) and of graphite microstructure when exposed to FLiNaK at 750 °C and 3 bar (gauge) pressure, starting from flowing argon at near atmospheric pressure. The results show that medium-grained and fine-grained graphites evolve to equilibrium at significantly different rates: fast salt uptake in medium-grained graphites produces salt deposits throughout the volume of graphite specimens, whereas salt infiltration in fine-grained graphites is much slower and limited to exposed surfaces.

FLiNaK infiltration↗

Mineral Scale Prevention on Electrically Conducting Membrane Distillation Membranes Using Induced Electrophoretic Mixing

The growth of mineral crystals on surfaces is a challenge across multiple industrial processes. Membrane-based desalination processes, in particular, are plagued by crystal growth (known as scaling), which restricts the flow of water through the membrane, can cause membrane wetting in membrane distillation, and can lead to the physical destruction of the membrane material. Scaling occurs when supersaturated conditions develop along the membrane surface due to the passage of water through the membrane, a process known as concentration polarization. To reduce scaling, concentration polarization is minimized by encouraging turbulent conditions and by reducing the amount of water recovered from the saline feed. In addition, antiscaling chemicals can be used to reduce the availability of cations. Here, we report on an energy-efficient electrophoretic mixing method capable of nearly eliminating CaSO4 and silicate scaling on electrically conducting membrane distillation (ECMD) membranes. The ECMD membrane material is composed of a percolating layer of carbon nanotubes deposited on porous polypropylene support and cross-linked by poly(vinyl alcohol). The application of low alternating potentials (2 Vpp,1Hz) had a dramatic impact on scale formation, with the impact highly dependent on the frequency of the applied signal, and in the case of silicate, on the pH of the solution.

54 ENVIRONMENTAL SCIENCES↗

On Black Holes as Macroscopic Quantum Objects

The relative flow of the Schwarzschild vs. the proper time during the classical evolution of a collapsing shell in the Schwarzschild coordinates practically forces us to interpret black hole formation as a highly non-local quantum process in which a shell/anti-shell pair is created within the incipient horizon, thus canceling out the original collapsing shell exactly at the horizon. By studying quantum fields in the black hole background, we reveal similar non-local effects. Among other things, the outgoing member of the Hawking pair very quickly becomes entangled with the black hole geometry (and not its partner), which is in contrast with the usual assumption that the Hawking pair is maximally entangled according to the local geometry near the horizon. Also, an infalling wave affects the black hole geometry even before it crosses the horizon. Finally, we find that a particle takes a finite amount of time to tunnel in and out of the black hole horizon, and thus avoids infinite blue and redshift in processes happening exactly at the horizon. These findings strongly support the picture of a black hole as a macroscopic quantum object.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

ThinCurr: An open-source 3D thin-wall eddy current modeling code for the analysis of large-scale systems of conducting structures

In this paper we present a new thin-wall eddy current modeling code, ThinCurr, for studying inductively-coupled currents in 3D conducting structures -- with primary application focused on the interaction between currents flowing in coils, plasma, and conducting structures of magnetically-confined plasma devices. The code utilizes a boundary finite element method on an unstructured, triangular grid to accurately capture device structures. The new code, part of the broader Open FUSION Toolkit, is open-source and designed for ease of use without sacrificing capability and speed through a combination of Python, Fortran, and C/C++ components. Scalability to large models is enabled through use of hierarchical off-diagonal low-rank compression of the inductance matrix, which is otherwise dense. Ease of handling large models of complicated geometry is further supported by automatic determination of supplemental elements through a greedy homology approach. Here, a detailed description of the numerical methods of the code and verification of the implementation of those methods using cross-code comparisons against the VALEN code and Ansys commercial analysis software is shown.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Microsized Pore Structure Determination in EPDM Rubbers Using High-Pressure 129 Xe NMR Techniques

Microsized pore parameters, such as pore size and distance between pores in a series of model EPDM rubbers, were determined in situ under the pressure of 500 psi using 129 Xe nuclear magnetic resonance (NMR) techniques: spin–lattice (T 1 ) and spin–spin (T 2 ) relaxation measurements, pulsed-field gradient (PFG) NMR, and two-dimensional exchange spectroscopy (2D EXSY). The T 1 /T 2 (>>1) ratio for the xenon confined in the pores is larger than that for nonconfined free xenon. This suggests that almost the entire pore surface interacts with xenon atoms like a closed pore. While these pores still connect each other through very narrow diffusion/exchange channels, it is possible to observe the echo decay in PFG-NMR and cross-peaks in 2D EXSY. The results show that both diffusion (D pore ≈ 2.1 × 10 –10 m 2 /s) and exchange (exchange rate, τ exch = a few tens of milliseconds) of xenon between a pore within the material and outer surface are prolonged. The exchange distances (l), which correspond to the xenon gas penetration depth, were estimated to be 70–100 μm based on the measured diffusion coefficients and exchange rate (1/τ exch ). NMR diffraction analysis reveals that pore size (a) and pore distance (b) are on the order of magnitude of micrometers and tens of micrometers, while the diffusion coefficients of xenon gas in the diffusion channels (D eff ) are about 10 –8 m 2 /s. Overall, this study suggests that the pores with a few micrometers connected through very narrow flowing channels with the length of several tens of micrometers are developed 70 to 100 μm below the rubber surface. Furthermore, the overall steady-state diffusion of xenon is slower, approximately 2 orders of magnitudes, than the diffusion in the channel between the pores. Finally, the pore and exchange distances correlated with the composition of rubbers showed that the properties of EPDM rubber as a high-pressure gas barrier could be improved by reducing the size of cracks and the depth of gas penetration by the addition of both carbon black and silica fillers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Towards Predictive Plasma Science and Engineering through Revolutionary Multi-Scale Algorithms and Models (Final Report)

This report describes the high-level accomplishments from the Plasma Science and Engineering Grand Challenge LDRD at Sandia National Laboratories. The Laboratory has a need to demonstrate predictive capabilities to model plasma phenomena in order to rapidly accelerate engineering development in several mission areas. The purpose of this Grand Challenge LDRD was to advance the fundamental models, methods, and algorithms along with supporting electrode science foundation to enable a revolutionary shift towards predictive plasma engineering design principles. This project integrated the SNL knowledge base in computer science, plasma physics, materials science, applied mathematics, and relevant application engineering to establish new cross-laboratory collaborations on these topics. As an initial exemplar, this project focused efforts on improving multi-scale modeling capabilities that are utilized to predict the electrical power delivery on large-scale pulsed power accelerators. Specifically, this LDRD was structured into three primary research thrusts that, when integrated, enable complex simulations of these devices: (1) the exploration of multi-scale models describing the desorption of contaminants from pulsed power electrodes, (2) the development of improved algorithms and code technologies to treat the multi-physics phenomena required to predict device performance, and (3) the creation of a rigorous verification and validation infrastructure to evaluate the codes and models across a range of challenge problems. These components were integrated into initial demonstrations of the largest simulations of multi-level vacuum power flow completed to-date, executed on the leading HPC computing machines available in the NNSA complex today. These preliminary studies indicate relevant pulsed power engineering design simulations can now be completed in (of order) several days, a significant improvement over pre-LDRD levels of performance.

42 ENGINEERING↗

Plug-and-play polymer microfluidic chips for hydrated, room temperature, fixed-target serial crystallography

The practice of serial X-ray crystallography (SX) depends on efficient, continuous delivery of hydrated protein crystals while minimizing background scattering. Of the two major types of sample delivery devices, fixed-target devices offer several advantages over widely adopted jet injectors, including: lower sample consumption, clog-free delivery, and the ability to control on-chip crystal density to improve hit rates. In this study we present our development of versatile, inexpensive, and robust polymer microfluidic chips for routine and reliable room temperature serial measurements at both synchrotrons and X-ray free electron lasers (XFELs). Our design includes highly X-ray-transparent enclosing thin film layers tuned to minimize scatter background, adaptable sample flow layers tuned to match crystal size, and a large sample area compatible with both raster scanning and rotation based serial data collection. The optically transparent chips can be used both for in situ protein crystallization (to eliminate crystal handling) or crystal slurry loading, with prepared samples stable for weeks in a humidified environment and for several hours in ambient conditions. Serial oscillation crystallography, using a multi-crystal rotational data collection approach, at a microfocus synchrotron beamline (SSRL, beamline 12-1) was used to benchmark the performance of the chips. Furthermore, high-resolution structures (1.3–2.7 Å) were collected from five different proteins – hen egg white lysozyme, thaumatin, bovine liver catalase, concanavalin-A (type VI), and SARS-CoV-2 nonstructural protein NSP5. Overall, our modular fabrication approach enables precise control over the cross-section of materials in the X-ray beam path and facilitates chip adaption to different sample and beamline requirements for user-friendly, straightforward diffraction measurements at room temperature.

59 BASIC BIOLOGICAL SCIENCES↗

HTGR Multiphysics Application Drivers FY26 Updates

This report summarizes FY26 progress under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program's high-temperature gas-cooled reactor (HTGR) application driver work, covering a wide range of activities such as code validation and multi-physics code assessment. 1) A detailed SAM model of the High-Temperature Engineering Test Reactor (HTTR) was developed using a unique-block grouping approach, with an extended parallel thermal network method to capture block-to-block conduction and radiation heat transfer, and applied to steady-state simulations of the HTTR 30~MW and 9~MW cases. 2) In another activity, SAM's newly implemented multi-component gas flow model was validated against the Natural convection Shutdown heat removal Test Facility (NSTF) argon ingress experiment, correctly capturing the density-driven suppression and thermal recovery of natural circulation observed when argon is introduced into the air-cooled Reactor Cavity Cooling System (RCCS) loop. 3) For the OECD/NEA High Temperature Test Facility (HTTF) benchmark, we co-led the international benchmark activities as well as the OECD/NEA final benchmark report to be released at the end of this year. 4) Finally, the coupled Griffin-SAM modeling capability for pebble-bed HTGRs was advanced by verifying the Griffin neutronics solution against Serpent Monte Carlo for a realistic non-uniform temperature distribution, resolving several deficiencies in the SAM-to-Griffin temperature transfer scheme, and enabling distinct fuel kernel, moderator, and coolant temperatures for cross section feedback. These new features were demonstrated in a PBR load-following transient.

Lee, Alvin↗

Noninvasive Imaging of a Liquid Jet

Liquid jets are found in many applications, from printing to manufacturing to entertainment. Herein this study uses three different noninvasive imaging modalities to compare resulting images of a liquid jet operating at three Reynolds numbers that cover laminar, transitional, and turbulent flow. Selected measurement quantities from each image type are also compared. High-speed backlit (BL) imaging is a simple imaging technique found in many laboratories, and this is compared to two high-speed X-ray imaging techniques, white beam (WB) imaging and focused beam (FB) radiography. BL imaging can provide a wide field of view and is easy to implement, but it only shows the presence or absence of liquid. WB imaging can show detailed contours on the surface of the liquid jet, but the imaging region is much smaller. FB radiography produces a point-source measurement and can provide the quantitative, instantaneous local liquid path length, termed the equivalent path length (EPL). All three techniques provide similar measures of jet thickness, with the FB measurements having less variation. FB measurements can also provide detailed cross sections of the average liquid jet thickness at high spatial resolutions.

42 ENGINEERING↗

Complexing Agent-Assisted Membraneless Zinc–Iodine Aqueous Batteries

A membrane is required for conventional zinc–iodine aqueous batteries, since soluble polyiodides cross over to the anode side and react with zinc metal spontaneously. Making the battery membraneless increases ion transport and reduces its cost and overall footprint. In this paper, a membraneless Zn–I 2 aqueous battery is demonstrated, employing a complexing agent, 1-butyl-1-methylpyrrolidinium iodide (MBPI), to promote the formation of I 5 – -containing, phase-separated polyiodides upon charging, to minimize self-discharge and suppress Zn dendrite growth. With an additional 0.3 M MBPI in 4 M ZnI 2 electrolyte, the membraneless battery achieved 65 cycles with >85% Coulombic efficiency, whereas the MBPI-free control failed immediately. Additionally, a volumetric capacity of 14.3 Ah L –1 was achieved, surpassing those of most membraneless batteries reported to date regardless of redox chemistry, and underscores the potential of complexing agents in simplifying the architecture of conventional Zn–I 2 flow batteries.

25 ENERGY STORAGE↗

Mesoscale fractal whey protein particles derived from microscale linear-shaped protein assemblies (Part 1): Manufacturing method and particle characteristics

Whey protein isolates (WPI) are widely used in processed foods for their versatile functional properties. Modifying the structural properties of proteins by assembling them into mesoscale or microscale particles may improve their functionality and broaden their applications. This study aims to manufacture and characterize mesoscale whey protein particles (WPP) derived from WPI. Two types of WPP, WPP1 (0.05 mL/min) and WPP2 (0.25 mL/min), were prepared through a multistep approach involving liquid antisolvent (LAS) precipitation, heat treatment, and microfluidization. Liquid antisolvent precipitation was performed by injecting a 20% (wt/vol) WPI dispersion (pH 7) into an ethanol-glycerol mixture (75:25, vol/vol) under laminar flow, followed by heat treatment at 80°C for 20 min as a particle hardening step. This process produced stable fiber- and ribbon-shaped whey protein assemblies (WPA), which served as precursors to WPP. Subsequent microfluidization (150 MPa, 6 passages) reduced the size of WPA, yielding mesoscale WPP with irregular morphologies and a more uniform size distribution, as revealed by microscopy and dynamic light scattering. ζ-Potential and fluorescence labeling indicated higher surface charge and surface hydrophobicity of WPP compared with untreated WPI. The WPP showed internal mass fractal and surface fractal structures at larger length scales, analyzed using small-angle X-ray scattering. Fourier transform infrared spectroscopy demonstrated an increased fraction of intermolecular β-sheets in WPP, suggesting that hydrogen bonding contributed to their formation. Gel electrophoresis confirmed that disulfide bonds served as the primary cross-links stabilizing the WPP structure. Furthermore, turbidity measurements showed that WPP exhibited superior colloidal phase stability compared with untreated WPI and maintained high colloidal stability under both acidic and neutral pH conditions.

Antisolvent precipitation↗

Reviews and syntheses: Variable inundation across Earth's terrestrial ecosystems

Abstract. The structure, function, and dynamics of Earth's terrestrial ecosystems are profoundly influenced by how often (frequency) and how long (duration) they are inundated with water. A diverse array of natural and human-engineered systems experience temporally variable inundation whereby they fluctuate between inundated and non-inundated states. Variable inundation spans extreme events to predictable sub-daily cycles. Variably inundated ecosystems (VIEs) include hillslopes, non-perennial streams, wetlands, floodplains, temporary ponds, tidal systems, storm-impacted coastal zones, and human-engineered systems. VIEs are diverse in terms of inundation regimes, water chemistry and flow velocity, soil and sediment properties, vegetation, and many other properties. The spatial and temporal scales of variable inundation are vast, ranging from sub-meter to whole landscapes and from sub-hourly to multi-decadal. The broad range of system types and scales makes it challenging to predict the hydrology, biogeochemistry, ecology, and physical evolution of VIEs. Despite all experiencing the loss and gain of an overlying water column, VIEs are rarely considered together in conceptual, theoretical, modeling, or measurement frameworks and approaches. Studying VIEs together has the potential to generate mechanistic understanding that is transferable across a much broader range of environmental conditions, relative to knowledge generated by studying any one VIE type. We postulate that enhanced transferability will be important for predicting changes in VIE function in response to global change. Here we aim to catalyze cross-VIE science that studies drivers and impacts of variable inundation across Earth's VIEs. To this end, we complement expert mini-reviews of eight major VIE systems with overviews of VIE-relevant methods and challenges associated with scale. We conclude with perspectives on how cross-VIE science can derive transferable understanding via unifying conceptual models in which the impacts of variable inundation are studied across multi-dimensional environmental space.

54 ENVIRONMENTAL SCIENCES↗

Elucidating the effect of preheating temperature on melt pool morphology variation in Inconel 718 laser powder bed fusion via simulation and experiment

In laser powder bed fusion (L-PBF) additive manufacturing, the mechanical performance, microstructure and defects of fabricated parts are closely associated with the melt pool morphology, e.g., its dimension and shape through the building process. Past studies have largely focused on how the process parameters such as laser power and scan speed affect melt pool characteristics. In this study, the melt pool morphology variation as a function of preheating temperature in the conduction, transition, and keyhole regimes and the underlying mechanisms in each regime are investigated through ex-situ sample characterization and computation thermal fluid dynamics (CtFD) simulation. Single tracks with different combinations of laser power and scan speed are deposited on an Inconel 718 bare plate preheated to a temperature range of 100–500 °C in the experiment. Significant changes are observed in melt pool morphology as a function of preheating temperature from optical measurements of melt track cross sections. The depth of melt pool in the three regimes increases monotonically with preheating temperature, e.g., at 500 °C, the experimental melt pool depth is increased by 49% in conduction regime, 34% in transition regime and 33% in keyhole regime, respectively, while the variation of melt pool width in each regime does not all follow an increasing trend but depends on the melt pool regimes. Melt pool width variation in the conduction and transition regimes is found to depend on the enhanced heat conduction directly related to temperature dependent thermal properties. Through validated CtFD simulations, it is found that in the keyhole regime the evaporation mass, recoil pressure, and laser drilling effect is enhanced with higher preheating temperature, which gives rise to a deeper melt pool. The simulations also reveal that preheating temperature significantly elongates the melt track length due to the increased flow rate and strong recoil pressure that accelerates the backward flow.

36 MATERIALS SCIENCE↗

The influence of thermomechanical treatment pathways on texture and mechanical properties in ARB Cu/Nb nanolaminates

Accumulative roll-bonded Cu/Nb nanolaminates (ARB Cu/Nb) possess high strength, thermal stability, and radiation tolerance arising from a high content of heterophase interfaces at fine layer thicknesses. These properties can be tailored by processing parameters used in the ARB Cu/Nb fabrication process, in which layer thickness, thermal history, and strain pathway determine the interface types and resultant properties found in the material. In this work, we subject ARB Cu/Nb to annealing, and then two different rolling pathways – one where rolling direction is held constant (longitudinal rolling, or LR), and one where rolling direction is rotated by 90º and held constant thereafter (cross rolling, or CR). Rolling is performed on ARB Cu/Nb over a targeted range of layer thicknesses from 193 to 25 nm and resultant bulk textures measured by neutron diffraction are correlated with mechanical properties measured by miniaturized tensile tests. The annealing procedure sharpens texture in both phases. We find that Cu mostly develops the same texture in LR and CR. In contrast, Nb develops a distinct texture along the CR pathway that is distinct from the LR texture. The composite texture of Cu/Nb is thus distinct between LR and CR pathways. This difference in texture development between Cu and Nb as a function of strain after change in rolling direction demonstrates the viability for deliberate pairing of Cu LR and Nb CR textures at a desired layer thickness. For mechanical properties, we find that differences in texture do not result in yield or flow stress differences above a layer thickness of 25 nm. Below a layer thickness of 25 nm, despite similar Taylor factors, yield and flow stress and are significantly different. In conclusion, this indicates texture only influences mechanical behavior at low layer thickness, where interface structure dominates mechanical properties.

36 MATERIALS SCIENCE↗

Divergence Reduction in Monte Carlo Neutron Transport with On-GPU Asynchronous Scheduling

While Monte Carlo Neutron Transport (MCNT) is near-embarrasingly parallel, the effectively unpredictable lifetime of neutrons can lead to divergence when MCNT is evaluated on GPUs. Divergence is the phenomenon of adjacent threads in a warp executing different control flow paths; on GPUS, it reduces performance because each work group may only execute one path at a time. The process of Thread Data Remapping (TDR) resolves these discrepancies by moving data across hardware such that data in the same warp will be processed through similar paths. A common issue among prior implementations of TDR is the synchronous nature of its remapping and processing cycles, which exhaustively sort data produced by prior processing passes and exhaustively evaluate the sorted data. In another work, we defined a method of remapping data through an asynchronous scheduler which allows for work to be stored in shared memory and deferred arbitrarily until that work is a viable option for low-divergence evaluation. This article surveys a wider set of cases, with the goal of characterizing performance trends across a more comprehensive set of parameters. These parameters include cross sections of scattering/capturing/fission, use of implicit capture, source neutron counts, simulation time spans, and tuned memory allocations. Across these cases, we have recorded minimum and average execution times, as well as a heuristically tuned near-optimal memory allocation size for both synchronous and asynchronous scheduling. Across the collected data, it is shown that the asynchronous method is faster and more memory efficient in the majority of cases, and that it requires less tuning to achieve competitive performance.

Computer Science↗

NSTX-U Machine Core Vacuum Seals Upgrade Design

The vacuum in the National Spherical Torus Experiment Upgrade (NSTX-U) is in the range of 2.66x10⁻⁶ Pa (2x 10⁻⁸ torr). To enhance vacuum seals and lower gas permeation in the NSTX-U and therefore ensuring operation safety and quality, the machine core vacuum seal design has been upgraded using double elastomeric seals and pumped interspaces (guard channels). In NSTX-U machine core, there are six major vacuum seals, including upper and lower vacuum vessel main flanges, upper and lower bellows flanges, and two in the ceramic isolation ring assembly that are mounted to the upper vacuum vessel main flange. The elastomeric seals are commercially available and customized fluorocarbon (Viton) O-rings. To ensure appropriate O-ring compression ratios for the vacuum seal, the O-ring grooves are customized with appropriate tolerances considering specific situations for each seal. In this upgrade, the vacuum in each interspace is designed to reach 13.33 Pa (0.1 torr), thus the overall pressure difference crossing each double O-ring seal is expected to be reduced by at least three orders of magnitude, and accordingly the overall rates of leakage and permeation will be lowered by three orders of magnitude. The corresponding pump-down times for each seal, at room temperature, are estimated at about 60 s. Additionally, the vacuum pumped interspaces will be able to function as a real-time leak monitor.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

CHARACTERIZATION AND DEVELOPMENT OF THE ALL–AQUEOUS COPPER THERMALLY REGENERATIVE AMMONIA BATTERY

Low temperature thermal energy is an abundant and untapped energy source for large scale electrical power and energy production around the world. Many different stationary sources such as industrial and thermal power plants produce significant amounts of low temperature heat that goes unused. The primary reason for the energy loss is because the temperature of the heat is often too low for traditional systems to economically generate power from this waste heat. Many different technologies are being researched to produce grid–scale power from low temperature heat to increase grid efficiency, advance decarbonization efforts, and progress electrical power production to a more distributed system. A thermally regenerative battery (TRB) is one promising technology in this space, as it has large power and energy densities relative to competitive technologies and therefore has the potential to disrupt the current power generation mix. TRBs combine two unit operations (a redox flow battery and a distillation column) which are easily scalable and mature, meaning that if a viable TRB chemistry can be identified, its development to pilot and commercial scale can be streamlined. The all–aqueous, copper thermally regenerative ammonia battery (Cuaq–TRAB) was the first TRB that did not use an electrodeposition/dissolution reaction by stabilizing Cu(I) and Cu(II) in the electrolyte, allowing all electroactive species to be stable in an aqueous phase. Through these electrolyte modifications, coulombic efficiency, open circuit potential, and copper solubility all increased, resulting in large gains to power density, energy density, and energy efficiency.

Cross, Nicholas↗

Metalorganic chemical vapor deposition of ZnGeN 2 films on GaN: effects of cation stoichiometry on surface morphology and crystallinity

Novel optoelectronic device designs based on the heterostructures of III-N and II-IV-N 2 are promising to advance device performance significantly. For example, by utilizing InGaN-ZnGeN 2 quantum well (QW) structures instead of pure InGaN, the band structure engineering in the QW active region can lead to improved electron-hole wavefunction overlap, and thus enhance the radiative efficiency for photon generation. These novel heterostructures have great potential to address the current challenge of low quantum efficiency in InGaN QW based light emitting diodes emitting in green and beyond. The materials development of ZnGeN 2 is still at an early stage as compared to the much matured GaN material system. In an ideal octet rule preserving ordered structure of ZnGeN 2 , every N atom is coordinated by exactly two Zn and two Ge atoms. However, local violation of octet rule can be caused by non-ideal coordination of N by the cations. The disordered structure is thermodynamically less favorable but can still be achieved, for example, in kinetics-limited growth regime. The ordered ZnGeN 2 has a bandgap very close to that of GaN (~3.4 eV) and a lattice mismatch of <0.1% with GaN. Interestingly, the valence band of ZnGeN 2 has been predicted to be ~1 eV above that of GaN, which has inspired novel designs for high efficiency light emitters. In this work, we investigated the metalorganic chemical vapor deposition (MOCVD) of ZnGeN 2 films on GaN/c-sapphire templates. Diethylzinc (DEZn), germane (GeH4) and ammonia were used as the precursors for Zn, Ge and N, respectively. A systematic study was conducted to investigate the cation stoichiometry as a function of growth temperature (TG), total reactor pressure (P) and DEZn/GeH 4 molar flow rate ratio (RII/IV). Under the investigated growth window, the Zn/(Zn+Ge) composition in the films, determined from energy dispersive X-ray spectroscopy, decreased monotonically with increase in TG but increased with increase in P and RII/IV. Atom probe tomography data did not indicate the presence of any secondary phases such as Zn 3 N 2 or Ge 3 N 4 . The surface morphology and crystallinity of the grown films had strong correlation with the Zn/(Zn+Ge) composition. The scanning electron microscopy images showed that the near-stoichiometric films have planar surfaces whereas Zn-rich films had crystallites on their surface and the Zn-poor films had faceted surface. Scanning transmission electron microscopy (STEM) imaging revealed that the Zn-rich and Zn-poor films have columnar and filament-like morphology, respectively, whereas the near-stoichiometric films have continuous film-like cross-sectional morphology. TEM nano-diffraction patterns as well as X-ray diffraction 2θ-ω scan profiles indicate that the near stoichiometric films are single crystalline. Nano-diffraction pattern of the stoichiometric films resembled that of a disordered ZnGeN 2 structure. Room temperature Raman spectra of near-stoichiometric films showed only the phonon density of states like features of a cation disordered ZnGeN 2 . Cathodoluminescence and photoluminescence spectra measured at different temperatures had similar features with peak emission wavelength at ~ 2 eV. In conclusion, the stoichiometry of ZnGeN 2 films can be widely tuned by tuning the MOCVD growth parameters. The surface morphology and the crystallinity of the films were found to have strong correlation with the Zn/(Zn+Ge) composition. The stoichiometric ZnGeN 2 films grown on GaN were demonstrated with uniform surface morphology and high crystalline quality. The results from this work will provide pathway to implement ZnGeN 2 in device structures.

Karim, Md Rezaul↗