Search NASA⌕ Search

SEARCH · Search NASA

Results for “Charge matching”

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 163 records · Page 9

Adaptive model tuning studies for non-invasive diagnostics and feedback control of plasma wakefield acceleration at FACET-II

Beam-driven plasma wakefield acceleration (PWFA) achieves the same energy gain in a single meter, for which conventional accelerators require several kilometers, however much work is still required to match the beam quality of conventional accelerators. The PWFA processes to be studied at the FACET-II facility will utilize extremely short (down to a few fs) high energy (10 GeV), high charge (few nC) electron and eventually positron bunches. The PWFA process is extremely sensitive to the detailed longitudinal current profiles of these bunches and it would be of great benefit to have precise measurement and control of these profiles. We present an adaptive model tuning technique for FACET-II to adaptively tune online models based on real time accelerator and beam data to continuously provide a non-invasive diagnostic which predicts the longitudinal phase space (LPS) of extremely short and highly compressed electron beams, which otherwise require destructive X-band transverse deflecting cavity (XTCAV)-based measurements. Based on simulation studies, our method has the potential for: (1). The development of a non-invasive longitudinal phase space diagnostic by adaptively tuning models based on non-invasive measurements such as energy spread spectra which could be recorded at all of the bunch compressors in the facility. (2). Utilize these diagnostic to perform model-independent feedback-based to achieve desired longitudinal phase space distributions. (3). Utilize the model-independent feedback approach to maximize energy gain while minimizing emittance growth and energy gain variance of the PWFA process by tuning accelerator parameters, while monitoring what longitudinal phase space the algorithm has found, information that will be useful for further analytical and simulation studies.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Isolated [B 2 (CN) 6 ] 2– : Small Yet Exceptionally Stable Nonmetal Dianion

Here, we report the observation of a small, yet remarkably stable, metal-free hexacyanodiborate dianion [B 2 (CN) 6 ] 2– in the gas phase. Negative ion photoelectron spectroscopy (NIPES) was employed to measure its spectra at multiple laser wavelengths, yielding a 1.9 eV electron binding energy (EBE) —a remarkably high value of electronic stability and a ~2.60 eV repulsive Coulomb barrier (RCB) for electron detachment. This rationalizes the observation of this dianion, although homolytic charge-separation dissociation into two [B(CN) 3 ] •– is energetically favorable. Quantum chemical calculations demonstrate a $D_{3d}$ staggered conformation for both the dianion and radical monoanion, and the calculated EBE and RCB match the experimental values well. The simulated density of states spectrum reproduces all measured electronic transitions, while the simulated vibrational progressions for the ground state transition cover a much narrower EBE range compared to the experimental band, indicating appreciable auto-photodetachment via electronically excited dianion resonances

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A transferable classical force field to describe glyme based lithium solvate ionic liquids

A non-polarizable force field for lithium (Li + ) and bis(trifluoromethanesulfonyl)imide (TFSI – ) ions solvated in diglyme at around 0.2 mol fraction salt concentration was developed based on ab initio molecular dynamics (AIMD) simulations and a modified polymer consistent force field model. A force–torque matching based scheme, in conjunction with a genetic algorithm, was used to determine the Lennard-Jones (LJ) parameters of the ion–ion and ion–solvent interactions. This force field includes a partial charge scaling factor and a scaling factor for the 1–4 interactions. The resulting force field successfully reproduces the radial distribution function of the AIMD simulations and shows better agreement compared to the unmodified force field. The new force field was then used to simulate salt solutions with glymes of increasing chain lengths and different salt concentrations. The comparison of the MD simulations, using the new force field, with experimental data at different salt concentrations and AIMD simulations on equimolar concentrations of the triglyme system demonstrates the transferability of the force field parameters to longer glymes and higher salt concentrations. Furthermore, the force field appears to reproduce the features of the experimental x-ray structure factors, suggesting accuracy beyond the first solvation shell, for equimolar salt solutions using both triglyme and tetraglyme as the solvent. Altogether, the new force field was found to accurately reproduce the molecular descriptions of LiTFSI-glyme systems not only at various salt concentrations but also with glymes of different chain lengths. Thus, the new force field provides a useful and accurate tool to perform in silico studies of this family of systems at the atomistic level.

25 ENERGY STORAGE↗

Epitaxial (Al x Ga 1-x-y In y ) 2 O 3 Alloys Lattice Matched to Monoclinic Ga 2 O 3 Substrates

We have epitaxially stabilized a series of monoclinic (Al x Ga 1-x-y In y ) 2 O 3 alloys by careful choice of molecular beam epitaxy growth conditions, which balance alloy growth with suboxide desorption. The films are pseudomorphic to (010) β-Ga 2 O 3 substrates at thicknesses up to 150 nm with compositions ranging from (Al 0.01 Ga 0.83 In 0.16 ) 2 O 3 to (Al 0.24 Ga 0.75 In 0.03 ) 2 O 3 . The absorption edge shifts from approximately 4.62-5.14 eV with coincidently increasing Al and decreasing In mole fractions. J-V measurements reveal an increase in resistivity over four orders of magnitude with a maximum value of 4.2 x 10 5 Ω-cm for (Al 0.17 Ga 0.76 In 0.07 ) 2 O 3 , which has nearly identical lattice parameters (both in-plane and out-of-plane) to the underlying β-Ga 2 O 3 . Scanning transmission electron microscopy of this sample reveals a mostly uniform and single crystalline film, though we identify areas of non-uniform In incorporation and some γ-phase inclusions. This work demonstrates the feasibility of thick layers lattice-matched to β-Ga 2 O 3 with increased bandgap compared to phase-separation limited (Al,Ga) 2 O 3 . These alloys can enable higher bandgap epitaxial dielectrics and high sheet charge density transistors by increasing the conduction band offset with respect to β-Ga 2 O 3 .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Raman gain control in optical fibers with orbital-angular-momentum-induced chirality of light

Stimulated Raman scattering is a particularly robust nonlinearity, occurring in virtually every material because its spectral linewidth and associated frequency shift do not typically depend on phases or directions (i.e. wavevectors) of the interacting light beams. In amorphous materials such as glass fibers, Raman bandwidths are large, enabling its use as a broadband gain element. This ubiquity makes it a versatile means for achieving optical amplification or realizing lasers over a large range of pulsewidths at user-defined colors. However, this ease of deploying the effect also presents itself as a stubborn source of noise in fiber-based quantum sources or parasitic emission in fiber lasers. Here, we show that orbital angular momentum carrying light beams experiencing spin-orbit interactions yield novel phase-matching criteria for Raman scattering. This enables tailoring its spectral shape (by over half the Raman shift in a given material) as well as strength (by ∼ 100×) simply by controlling light’s topological charge – a capability of utility across the multitude of applications where modulating Raman scattering is desired.

Liu, Xiao (ORCID:000000032239762X)↗

16.52% Efficiency All-Polymer Solar Cells with High Tolerance of the Photoactive Layer Thickness

All-polymer solar cells (all-PSCs) have drawn growing attention and achieved tremendous progress recently, but their power conversion efficiency (PCE) still lags behind small-molecule-acceptor (SMA)-based PSCs due to the relative difficulty on morphology control of polymer photoactive blends. Here, in this work, low-cost PTQ10 is introduced as a second polymer donor (a third component) into the PM6:PY-IT blend to finely tune the energy-level matching and microscopic morphology of the polymer blend photoactive layer. The addition of PTQ10 decreases the π–π stacking distance, and increases the π–π stacking coherence length and the ordered face-on molecular packing orientation, which improves the charge separation and transport in the photoactive layer. Moreover, the deeper highest occupied molecular orbital energy level of the PTQ10 polymer donor than PM6 leads to higher open-circuit voltage of the ternary all-PSCs. As a result, a PCE of 16.52% is achieved for ternary all-PSCs, which is one of the highest PCEs for all-PSCs. In addition, the ternary devices exhibit a high tolerance of the photoactive layer thickness with high PCEs of 15.27% and 13.91% at photoactive layer thickness of ≈205 and ≈306 nm, respectively, which are the highest PCEs so far for all-PSCs with a thick photoactive layer.

14 SOLAR ENERGY↗

Many-body interference in kagome crystals

When electrons in metals act collectively, they enable emergent phenomena and electronic functionalities that transcend the behaviour of individual particles. Coherent collective charge motion has so far been observed primarily in superconductors, in which it arises with the formation of Cooper pairs. Here we report experimental evidence for coherent charge transport in the normal state of the kagome metal CsV 3 Sb 5 , indicative of a distinct collective electronic state. The signature is a set of magnetoresistance oscillations in mesoscopic crystalline pillars under in-plane magnetic fields, with a periodicity determined by the number of magnetic flux quanta h/e threading between adjacent kagome layers—effectively forming an interlayer Aharonov–Bohm interferometer. The cooperative nature of this phenomenon is evidenced by a non-analytic angular dependence characterized by abrupt transitions between discrete oscillation frequencies and its persistence over length scales that exceed the single-particle mean free path. Notably, the oscillation amplitude matches other anomalous electronic responses reported in CsV 3 Sb 5 , pointing to an underlying mechanism that establishes intrinsic coherence. These findings shed new light on the debated nature of correlated order in kagome metals and establish CsV 3 Sb 5 as a platform for realizing long-range coherent charge transport in the absence of superconductivity—opening new directions for coherence in correlated electron systems beyond conventional models.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

General kinetic ion-induced electron emission model for metallic walls applied to biased Z-pinch electrodes

A kinetic ion-induced electron emission (IIEE) model for general applications is developed to obtain the emitted electron energy spectrum for a distribution of ion impacts on a metallic surface. We assume an ionization cascade mechanism and use empirical models for the ion and electron stopping powers. The emission spectrum and the secondary electron yield (SEY) are validated for a variety of materials. The IIEE model is used to study the effect of IIEE on the plasma-material interactions of Z-pinch electrodes. Un-magnetized Boltzmann-Poisson simulations are performed for a Z-pinch plasma doubly bounded by two biased copper electrodes with and without IIEE at bias potentials from 0 to 9 kV. At the anode, the SEY decreases from 0 to 1 kV, but then increases at higher bias potentials. At the cathode, the SEY is much larger due to higher energy ion bombardment and grows with bias potential. As the bias potential increases, the emitted cathode electrons are accelerated to higher energies into the domain, collisionally heating the plasma. Above 1 kV, the heating is strong enough to increase the plasma potential. Despite SEY greater than 1, only a classical sheath forms as opposed to a space-charge limited or inverse sheath due to the emitted electron flux not reaching the space charge current saturation limits. Furthermore, the current in the emissionless cases saturates to a value lower than experiment. With IIEE, the current does not saturate and continues to increase with the 4 kV case, matching most closely with the experiment.

Carbon based materials↗

Augmented signal processing in Liquid Argon Time Projection Chambers with a deep neural network

The Liquid Argon Time Projection Chamber (LArTPC) is an advanced neutrino detector technology widely used in recent and upcoming accelerator neutrino experiments. It features a low energy threshold and high spatial resolution that allow for comprehensive reconstruction of event topologies. In current-generation LArTPCs, the recorded data consist of digitized waveforms on wires produced by induced signal on wires of drifting ionization electrons, which can also be viewed as two-dimensional (2D) (time versus wire) projection images of charged-particle trajectories. For such an imaging detector, one critical step is the signal processing that reconstructs the original charge projections from the recorded 2D images. For the first time, we introduce a deep neural network in LArTPC signal processing to improve the signal region of interest detection. By combining domain knowledge (e.g., matching information from multiple wire planes) and deep learning, this method shows significant improvements over traditional methods. This work details the method, software tools, and performance evaluated with realistic detector simulations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

All-Atom Molecular Dynamics Simulations of Cationic Polyelectrolyte Brushes in the Presence of Halide Counterions

Understanding the response of the charged polyelectrolyte (PE) brushes and brush-supported water and ions to the changes in the nature of screening counterions is significant in developing strategies for utilizing PE brushes in various applications. In this paper, we employ all-atom molecular dynamics (MD) simulations for studying the behavior of cationic [poly(2-(methacryloyloxy)ethyl) trimethylammonium] (PMETA) brushes and brush-supported water and ions in the presence of different halide (X: I – , Br – , Cl – , and F – ) screening counterions. We find that despite the F– ion having the largest charge density, the extent of binding of the counterions on the PMETAX brushes varies as I – > Br – > Cl – > F – , leading to PMETAX brush height being least with I– counterions and greatest with F – counterions. This trend in the binding of the halide ions matches the previous experimental result and can be explained by identifying the chaotropic nature of the I – and Br – ions that promote a disruption of water structure and a more favorable binding of the ions to the polymer chains. Such a binding trend also ensures that the order of the water molecules around the PMETAX chains or the counterions as well as the number of water–water hydrogen bonds inside the brush layer increases in the following order of the counterion-specific PMETAX brushes: F – > Cl – > Br – > I – . Furthermore, halide-ion-PMETAX-chain binding takes place via both interchain and intrachain bridging: intrachain bridging dominates for the case of counterions that show enhanced binding (I – and Br – ), while interchain bridging is more favored for the case of counterions that show weakened binding (F – and Cl – ). Also, a greater degree of intrachain bridging leads to greater compressibility and flexibility of the brush layer. Lastly, we show that the mobility of the halide ions follows a nonmonotonic trend with the charge density: the mobility decreases as I – < Br – < Cl – (as their binding to the PMETAX chains varies as I – > Br – > Cl – ), but the mobility of F – ions is in between that of I – and Br – ions. We argue that the strongly attached hydration layer and the ensuing friction inside the brush layer lead to such a reduced mobility of the F – ions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Magnetic excitations and absence of charge order in the van der Waals ferromagnet Fe 4.75 GeTe 2

Understanding the ground state of van der Waals (vdW) magnets is crucial for designing devices leveraging these platforms. Here, we investigate the magnetic excitations and charge order in Fe 4.75 GeTe 2 , a vdW ferromagnet with ≈ 315 K Curie temperature. Using Fe L 3 -edge resonant inelastic X-ray scattering, we observe a dual nature of magnetic excitations, comprising a coherent magnon and a broad non-dispersive continuum extending up to 150 meV, 50% higher than in Fe 2.72 GeTe 2 . The continuum intensity is sinusoidally modulated along the stacking direction L, with a period matching the inter-slab distance. Our results indicate that while the dual character of the magnetic excitations is generic to Fe-Ge-Te vdW magnets, Fe 4.75 GeTe 2 exhibits a longer out-of-plane magnetic correlation length, suggesting enhanced 3D magnetic character. Furthermore, resonant X-ray diffraction reveals that previously reported ±(1/3, 1/3, L) peaks originate from crystal structure rather than from charge order.

Ferromagnetism↗

Probing the Electronic Structure of [B 10 H 10 ] 2– Dianion Encapsulated by an Octamethylcalix[4]pyrrole Molecule

Despite being an important closo-borate in the condensed phase boron chemistry, isolated B 10 H 10 2- is electronically unstable and has never been detected in the gas phase. Herein, we report a successful capture of this fleeting species through binding with an octamethylcalix[4]pyrrole (omC4P) molecule to form a stable gaseous omC4P·[B 10 H 10 ] 2- complex and its characterizations utilizing negative ion photoelectron spectroscopy (NIPES). The recorded NIPE spectrum, contributed from both omC4P and [B 10 H 10 ] 2- , is deconvoluted by subtracting the omC4P contribution to yield a [B 10 H 10 ] 2- spectrum. The obtained [B 10 H 10 ] 2- spectrum consists of four major bands spanning electron binding energy (EBE) range from 1 to 5 eV with the EBE gaps matching excellently with the energy intervals of computed highly lying occupied molecular orbitals of the B 10 H 10 2- dianion. Finally, this study showcases a generic method to utilize omC4P to capture unstable multiply charged anions in the gas phase for experimental determination of their electronic structures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Vanadate Retention by Iron and Manganese Oxides

Anthropogenic emissions of vanadium (V) into terrestrial and aquatic surface systems now match those of geogenic processes, and yet, the geochemistry of vanadium is poorly described in comparison to other comparable contaminants like arsenic. In oxic systems, V is present as an oxyanion with a +5 formal charge on the V center, typically described as H x VO 4 (3–x)– , but also here as V(V). Iron (Fe) and manganese (Mn) (oxy)hydroxides represent key mineral phases in the cycling of V(V) at the solid–solution interface, and yet, fundamental descriptions of these surface-processes are not available. Here, we utilize extended X-ray absorption fine structure (EXAFS) and thermodynamic calculations to compare the surface complexation of V(V) by the common Fe and Mn mineral phases ferrihydrite, hematite, goethite, birnessite, and pyrolusite at pH 7. Inner-sphere V(V) complexes were detected on all phases, with mononuclear V(V) species dominating the adsorbed species distribution. Our results demonstrate that V(V) adsorption is exergonic for a variety of surfaces with differing amounts of terminal –OH groups and metal–O bond saturations, implicating the conjunctive role of varied mineral surfaces in controlling the mobility and fate of V(V) in terrestrial and aquatic systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Characterizing W sources in the all-W wall, all-RF WEST tokamak environment

In this work, experimental data, together with interpretive modeling tools, are examined to study trends in the tungsten (W) source in the all-W environment of the WEST tokamak, both from the divertor and from the main chamber. In particular, a poloidal limiter protecting an ion cyclotron resonance heating (ICRH) antenna is used as proxy for main chamber sourcing. The key study is carried out by stepping up lower hybrid current drive (LHCD) power, as the only auxiliary power source. Limiter and divertor W sources exhibit a qualitatively similar proportionality to the total power crossing the separatrix, P SEP , although the main chamber source remains substantially lower than the divertor source, for the range of P SEP accessible in the experiments. Intepretive modeling of the limiter source is carried out with a particle-in-cell (PIC) sheath model coupled to a surface sputtering model. Oxygen is used as a proxy for all light impurity species allowing for characterization of the critical W erosion regions. To get a good quantitative match to the data, it is necessary to assume that the oxygen arrives at the surface mostly at high ionization stages (4+ and above). A separate simulation with SOLEDGE-EIRENE, constrained to measured upstream scrape-off-layer plasma profiles, gives oxygen fractional abundances that are compatible with the PIC simulation result. This is understood to arise from transport processes that dominate over recombination. Substituting the LHCD by ICRH, in an equivalent experiment, the local W source exhibits a 3× enhancement. This can be matched by the simulation, by assuming local RF electric field rectification, based on ~100 eV peak-to-peak, near-antennna electric field. This work has highlighted the particular importance of understanding the ion charge state balance of light impurities as these are most likely the dominant sputtering species in fusion devices with high-Z walls.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

QED with the electron as a very small black-hole-like object

We pursue a simple semi-classical particle-based picture of QED in which the electron is assumed to contain a point massless and chargeless core that interacts with the electromagnetic vacuum in ways that are not dissimilar to the manner in which small black holes may interact with photons. Modifications are made to obtain known properties of the electron including the presence of rest energy, Lamb shift, magnetic moments (both free and bound), charge, Larmor emission, and vacuum polarization. The present study does not provide proofs but is rather an attempt to reverse-engineer the photon-electron and black-hole-electron interaction properties needed to match experimental observations. The ease by which many QED observables can be reproduced suggests further study of the presented framework may be warranted. Future detailed analysis of the exchange of long wavelength Hawking radiation between a pair of small black holes will be of particular interest.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Multi-modal Energy-optimal Trip Scheduling in Real-time (METS-R) for Transportation Hubs (Final Report)

This report summarizes the work performed under the award number EE0008524. The project develops the Multi-modal Energy-optimal Trip Scheduling in Real-time (METS-R) platform as the next-generation transportation solution based on autonomous electric vehicles (AEV) serving passenger trips from and to urban transportation hubs, to substantially reduce transportation energy consumption. Extensive data collection and analyses were first conducted to understand the demand patterns and energy consumption of hub-based on-road trips. Then, a data-driven framework that consists of an analytical module and a simulation module was proposed. For the analytical module, five planning + operation tools were developed to support the planning and energy-efficient operations of urban AEV services: the charging station planning that robotically allocates charging supplies based on the stationary charging demand distribution; the transit planning and demand adaptive scheduling model that efficiently generates\ candidate transit routes from hubs to other places and dynamically adjusts the transit time table to fit the current demand; the online energy-efficient routing that learns the energy-optimal paths from observations of link-level energy consumption in real-time; the hub-based ridesharing that matches trip requests together with account for the uncertainty of future trip demand and vehicle supply; and finally, the integrated demand prediction and anomaly detection pipeline that leverages the flight/train time table and support other planning/operation tools. To demonstrate the performance of these tools, a scalable high-performance agent-based simulator was built. We divided the urban space into multiple service zones where each zone was considered as an agent for passenger generation and vehicle charging. Two types of AEV agents were coded to model two types of mobility services: AEV taxi and AEV transit. For the AEV taxi, the team implemented the functions of pickup/drop-off passengers, energy-efficient routing, ridesharing, fleet rebalancing, and recharging. For the AEV bus, the team implemented the functions of demand-adaptive route scheduling, passenger boarding, and recharging. A high-performance computing framework was introduced to receive various profiling information (such as link energy updates, vehicle speed) from the simulator instances and communicate the operational commands back to the instances. The numerical experiments show that each of the proposed operational algorithms can reduce energy consumption and improve system efficiency. Furthermore, there exists the need to collectively consider multiple planning + operational strategies as multiple strategies can influence each other in terms of performance impacts. Recommendations for future work related to AEV planning and simulation are discussed.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

DFT simulation of the X-ray diffraction pattern of aluminum-ion-intercalated graphite used as the cathode material of the aluminum-ion battery

Many breakthroughs have been achieved in rechargeable aluminum-ion battery technologies in recent years. Most recently, operando X-ray diffraction (XRD) combined with density functional theory (DFT) calculations was reported to study the chloroaluminate anion (AlCl 4 -)-intercalated graphite cathode of the battery. However, there are quite a few discrepancies between the measured and simulated XRD patterns. This work is focused on the simulation of XRD patterns of graphite intercalation compounds (GICs) with DFT calculations. Our results reveal that both the geometry of AlCl 4 - in graphite and the gallery height of GICs are dependent on the intercalant density. At low intercalant density, the gallery height keeps constant, but at high intercalant densities, the gallery height is linearly related to the intercalant density. Our simulated XRD patterns are highly consistent with the measured operando XRD patterns. Not only do the angles of the peaks match very well, but also the relative intensities and the corresponding electrode capacities show reasonable agreement with the experimental results. The DFT simulation of the XRD pattern provides significant information on the stage index and the charge capacity of the GIC electrode.

25 ENERGY STORAGE↗

Charge Transfer Plasmonics in Bespoke Graphene/α-RuCl 3 Cavities

Surface plasmon polaritons (SPPs) provide a window into the nano-optical, electrodynamic response of their host material and its dielectric environment. Graphene/α-RuCl 3 serves as an ideal model system for imaging SPPs since the large work function difference between these two layers facilitates charge transfer that hole dopes graphene with n ~ 10 13 cm –2 free carriers. In this work, we study the emergent THz response of graphene/α-RuCl 3 heterostructures using our home-built cryogenic scanning near-field optical microscope. Using phase-resolved imaging, we clearly observe long wavelength, heavily damped THz SPPs in a series of variable-size graphene cavities. From this, we extract the plasmonic wavelength and scattering rate in the graphene/α-RuCl 3 heterostructures. We determine that the measured plasmon wavelength and electronic scattering rate match our heterostructures’ theoretically predicted values. Here, our results demonstrate that shaping graphene into bespoke cavity structures enables observation and quantification of SPPs in heavily doped graphene that are largely not addressable with other experimental techniques. Moreover, the manifest lack of metallicity observed in the adjacent doped α-RuCl 3 layer provides significant constraints on the nature of the interfacial charge transfer in this 2D heterostructure.

36 MATERIALS SCIENCE↗