Search NASA⌕ Search

SEARCH · Search NASA

Results for “weak solutions”

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.

485 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

4-Clique network minor embedding for quantum annealers

Quantum annealing is a quantum algorithm for computing solutions to combinatorial optimization problems. This study proposes a method for minor embedding optimization problems onto sparse quantum annealing hardware graphs called 4-clique network minor embedding. This method is in contrast to the standard minor embedding technique of using a path of linearly connected qubits in order to represent a logical variable state. The 4-clique minor embedding is possible on Pegasus graph connectivity, which is the native hardware graph for some of the current D-Wave quantum annealers. The Pegasus hardware graph contains many cliques of size 4, making it possible to form a graph composed entirely of paths of connected 4-cliques on which a problem can be minor-embedded. The 4-clique chains come at the cost of additional qubit usage on the hardware graph, but they allow for stronger coupling within each chain, thereby increasing chain integrity, reducing chain breaks, and allow for greater usage of the available energy scale for programming logical problem coefficients on current quantum annealers. The 4-clique minor embedding technique is compared with the standard linear path minor embedding with experiments on two D-Wave quantum annealing processors with Pegasus hardware graphs. We show proof-of-concept experiments where the 4-clique minor embeddings can use weak chain strengths while successfully carrying out the computation of minimizing random all-to-all spin glass problem instances. Published by the American Physical Society 2024

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Dynamic Modeling and Characterization of Nuclear-grade Graphite

Idaho National Labs serves as the spearhead for many innovative energy solutions to the world's energy crisis. One such solution is the INL's Microreactor which is designed to deploy to extreme/remote environments where other sources of power are either unavailable or unreliable. In order to best design these energy solutions for their operational environments, it is crucial to understand how the design, components, and materials will respond to the environmental conditions. One key material in these innovative designs is a nuclear-grade graphite known as PCEA. This study examines the behavior of PCEA graphite under dynamic loading, similar to that which may occur in extreme environments. The objective is to characterize the dynamic behavior and produce an accurate, reliable constitutive model suitable for use in simulation tools such as INL's MOOSE. Graphite specimens were tested using a Split Hopkinson Pressure Bar (SHPB) to administer the dynamic compressive load. The SHPB was charged at various pressures to produce a range of strain rates on the material in compression. Data was acquired via strain gauges on the SHPB setup, from which stress, strain, and time data were collected. Analysis revealed the stress-strain behavior of the material as well as insights into the material behavior's relationship to strain rate. Further work must continue to characterize the various other dynamic behaviors of the material which will combine to create a substantially trustworthy constitutive model for this grade of nuclear-grade graphite. Ultimately, this will allow for realistic simulation of the material in reactor designs, allowing for prediction of design weaknesses and leading to improved designs for increased resilience, security, and reliability.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Predicting Excitation Energies of Twisted Intramolecular Charge-Transfer States with the Time-Dependent Density Functional Theory: Comparison with Experimental Measurements in the Gas Phase and Solvents Ranging from Hexanes to Acetonitrile

Electronically excited states characterized by intramolecular charge transfer play an essential role in many biological processes and optical devices. The ability to make quantitative ab initio predictions of the relative energetics involved is a challenging yet desirable goal, especially for large molecules in solution. In this work, we present a data set of 61 experimental measurements of absorption and emission processes, both in the gas phase and in solvents representing a broad range of polarities, which involve intramolecular charge transfer mediated by a nonzero, "twisted" dihedral angle between one or more donor and acceptor subunits. Among a variety of density functionals investigated within the framework of linear-response theory, the "optimally tuned" LRC-ωPBE functional, which utilizes a system-specific yet nonempirical procedure to specify the range-separation parameter, emerges as the preferred choice. For the entire set of excitation energies, involving changes in dipole moment ranging from 4 to >20 Debye, the mean signed and absolute errors are 0.02 and 0.18 eV, respectively (compared, e.g., to -0.30 and 0.30 for PBE0, 0.44 and 0.47 for LRC-ωPBEh, 0.83 and 0.83 for ωB97X-V). We analyze the performance of polarizable continuum solvation models available in Q-Chem that partition the solvent response into fast and slow time scales, and clear trends emerge when measurements corresponding to the four small 4-(dimethylamino)benzonitrile (DMABN)-like molecules and a charged species are excluded. We make the case that the large errors found only for small molecules in the gas phase and weak solvents cannot be expected to improve via the optimal tuning procedure, which enforces a condition that is exact only in the well-separated donor-acceptor limit, and present empirical evidence implicating the outsized importance for small donor-acceptor systems of relaxation effects that cannot be accounted for by the linear-response time-dependent density functional theory within the adiabatic approximation. Finally, we demonstrate the utility of the optimally tuned density functional approach by targeting the charge-transfer states of a large biomimetic model system for light-harvesting structures in Photosystem II.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Revealing core-valence interactions in solution with femtosecond X-ray pump X-ray probe spectroscopy

Abstract Femtosecond pump-probe spectroscopy using ultrafast optical and infrared pulses has become an essential tool to discover and understand complex electronic and structural dynamics in solvated molecular, biological, and material systems. Here we report the experimental realization of an ultrafast two-color X-ray pump X-ray probe transient absorption experiment performed in solution. A 10 fs X-ray pump pulse creates a localized excitation by removing a 1 s electron from an Fe atom in solvated ferro- and ferricyanide complexes. Following the ensuing Auger–Meitner cascade, the second X-ray pulse probes the Fe 1 s → 3 p transitions in resultant novel core-excited electronic states. Careful comparison of the experimental spectra with theory, extracts +2 eV shifts in transition energies per valence hole, providing insight into correlated interactions of valence 3 d with 3 p and deeper-lying electrons. Such information is essential for accurate modeling and predictive synthesis of transition metal complexes relevant for applications ranging from catalysis to information storage technology. This study demonstrates the experimental realization of the scientific opportunities possible with the continued development of multicolor multi-pulse X-ray spectroscopy to study electronic correlations in complex condensed phase systems.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

The “good,” the “bad,” and the “hidden” in neutron scattering and molecular dynamics of ionic aqueous solutions

We characterize a concentrated 7.3 m CaCl 2 solution, combining neutron diffraction with chloride isotopic substitution (Cl-NDIS) in null water and molecular dynamics (MD) simulations. We elucidate the solution structure, thermodynamic properties, and extent of ion pairing previously suggested as concentration-dependent and often not observed at lower concentrations. Our Cl-NDIS measurements designate the solvent-shared ion pairing as dominant and the contact ion pairing (CIP) as insignificant even under conditions close to the solubility limit. The MD models parameterized against neutron diffraction with calcium isotopic substitution (Ca-NDIS) overestimate CIP despite successfully reproducing most of the Cl-NDIS signal. This drawback originates from the fact that Ca 2+ –Cl - interactions were primarily “hidden” in the Ca-NDIS signal due to overlapping with Ca 2+ –O w and Ca 2+ –H w contributions to the total scattering. Contrary, MD models with moderate CIP and possessing generally good performance at high concentrations fail to reproduce the NDIS measurements accurately. Therefore, the electronic polarization, introduced in most of the recent MD models via scaling ionic charges, resolves some but not all parameterization drawbacks. We conclude that despite improving the quality of MD models “on average,” the question “which model is the best” has not been answered but replaced by the question “which model is better for a given research.” An overall “good” model can still be inappropriate or, in some instances, “bad” and, unfortunately, produce erroneous results. Additionally, the accurate interpretation of several NDIS datasets, complemented by MD simulations, can prevent such mistakes and help identify the strengths, weaknesses, and convenient applications for corresponding computational models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

FROSch Preconditioners for Land Ice Simulations of Greenland and Antarctica

Numerical simulations of Greenland and Antarctic ice sheets involve the solution of large-scale highly nonlinear systems of equations on complex shallow geometries. This work is concerned with the construction of Schwarz preconditioners for the solution of the associated tangent problems, which are challenging for solvers mainly because of the strong anisotropy of the meshes and wildly changing boundary conditions that can lead to poorly constrained problems on large portions of the domain. In this study, two-level generalized Dryja--Smith--Widlund (GDSW)--type Schwarz preconditioners are applied to different land ice problems, i.e., a velocity problem, a temperature problem, as well as the coupling of the former two problems. We employ the message passing interface (MPI)--parallel implementation of multilevel Schwarz preconditioners provided by the package FROSch (fast and robust Schwarz) from the Trilinos library. The strength of the proposed preconditioner is that it yields out-of-the-box scalable and robust preconditioners for the single physics problems. To the best of our knowledge, this is the first time two-level Schwarz preconditioners have been applied to the ice sheet problem and a scalable preconditioner has been used for the coupled problem. The preconditioner for the coupled problem differs from previous monolithic GDSW preconditioners in the sense that decoupled extension operators are used to compute the values in the interior of the subdomains. Additionally, several approaches for improving the performance, such as reuse strategies and shared memory OpenMP parallelization, are explored. In our numerical study we target both uniform meshes of varying resolution for the Antarctic ice sheet as well as nonuniform meshes for the Greenland ice sheet. We present several weak and strong scaling studies confirming the robustness of the approach and the parallel scalability of the FROSch implementation. Among the highlights of the numerical results are a weak scaling study for up to 32,K processor cores (8,K MPI ranks and 4 OpenMP threads) and 566,M degrees of freedom for the velocity problem as well as a strong scaling study for up to 4,K processor cores (and MPI ranks) and 68,M degrees of freedom for the coupled problem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

FROSch Preconditioners for Land Ice Simulations of Greenland and Antarctica

Numerical simulations of Greenland and Antarctic ice sheets involve the solution of large-scale highly nonlinear systems of equations on complex shallow geometries. This work is concerned with the construction of Schwarz preconditioners for the solution of the associated tangent problems, which are challenging for solvers mainly because of the strong anisotropy of the meshes and wildly changing boundary conditions that can lead to poorly constrained problems on large portions of the domain. Here, two-level GDSW (Generalized Dryja–Smith–Widlund) type Schwarz preconditioners are applied to different land ice problems, i.e., a velocity problem, a temperature problem, as well as the coupling of the former two problems. We employ the MPI-parallel implementation of multi-level Schwarz preconditioners provided by the package FROSch (Fast and Robust Schwarz)from the Trilinos library. The strength of the proposed preconditioner is that it yields out-of-the-box scalable and robust preconditioners for the single physics problems. To our knowledge, this is the first time two-level Schwarz preconditioners are applied to the ice sheet problem and a scalable preconditioner has been used for the coupled problem. The pre-conditioner for the coupled problem differs from previous monolithic GDSW preconditioners in the sense that decoupled extension operators are used to compute the values in the interior of the sub-domains. Several approaches for improving the performance, such as reuse strategies and shared memory OpenMP parallelization, are explored as well. In our numerical study we target both uniform meshes of varying resolution for the Antarctic ice sheet as well as non uniform meshes for the Greenland ice sheet are considered. We present several weak and strong scaling studies confirming the robustness of the approach and the parallel scalability of the FROSch implementation. Among the highlights of the numerical results are a weak scaling study for up to 32 K processor cores (8 K MPI-ranks and 4 OpenMP threads) and 566 M degrees of freedom for the velocity problem as well as a strong scaling study for up to 4 K processor cores (and MPI-ranks) and 68 M degrees of freedom for the coupled problem.

58 GEOSCIENCES↗
DOE OSTI2020

Scalable line and plane relaxation in a parallel structured multigrid solver

The efficient solution of sparse, linear systems that arise through the discretization of partial differential equations remains a key challenge for a range of high performance scientific simulations. One approach for reducing data movement and improving performance is by exposing and exploiting structure in a problem through the use of robust structured multilevel solvers. By choosing coarsening that preserves the structure of the problem, these methods maintain efficient structured computation and communication throughout the multigrid hierarchy. However, when coarsening is not permitted to be dependent on the operator, anisotropy must be addressed by the smoother — producing error compatible for coarse-grid correction with structured coarsening. Here, the components required in a scalable parallel structured solver are described with a focus on memory and communication efficiency of robust smoothers. While the implementation of communication and memory reduction techniques in smoothers integrated in a complete 3D solver present a significant engineering challenge, a novel approach is proposed that addresses these challenges systematically through a change to the solver’s execution model. Enabled by user-level threading paired with a set of data and communication abstractions, this approach permits seamless aggregation of communication in plane smoothers — directly reusing code for a 2D distributed multilevel cycle. Results show an effective reduction in communication costs for coarse-grid problems, and result in a speedup of 8.7x in smoothing routines shown in Fig. 12 using this approach. This produces a significant improvement to strong scalability while maintaining favorable weak scaling behavior. Finally, a parallel scaling study using a series of refined meshes is included that demonstrates the effectiveness of this approach in an application of interest.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Effect of sintering temperature on adhesion of spray-on piezoelectric transducers

Conventionally sol-gel spray-on transducers require a high-temperature (> 700 ◦C) sintering process; however, this process can affect the microstructure of the substrate material. For mechanical elbows and valves utilized for fluid transport in the energy sector, the components are designed to have a specific microstructure, and deviations from these specifications can create weak points in the system. For this reason it is important to investigate how the temperature of the deposition process affects the substrate. This paper investigates the effect of high-temperature and low-temperature (< 150 ◦C) processing conditions on the surface composition of the substrate. Furthermore, the resultant transducers from high- and low-temperature fabrication processes are compared to determine if a low-temperature processing method is feasible. For these studies a sol-gel spray-on process is employed to deposit piezoelectric ceramics onto a stainless-steel 316L substrate. Energy-dispersive X-ray spectroscopy is utilized to determine the composition of the substrate surface before and after transducer deposition. Results indicate that the high-temperature processing conditions may alter the surface composition of the metal due to a diffusion of the metal into the ceramic, which results in a metal surface that is bonded to the ceramic. Furthermore, it is shown that low-temperature processing of spray-on transducers is a viable method for transducer fabrication where the resultant transducers meet the industry minimum requirement of 30 dB signalto-noise ratio. In parallel simulation calculations, finite-element method (FEM) studies were performed to model the adhesive strength of the low-temperature processed transducer to the substrate surface. Comparisons between the simulations and experiments suggest that the bond strength is much greater than the commercial gel bonds and closer to hardened epoxy glue bonds. These results indicate that spray-on transducers fabricated under lowtemperature processing conditions are a viable solution for leave-in-place monitoring of structures.

M. Sinding, Kyle↗
DOE OSTI2021

Improving the astrometric solution of the Hyper Suprime-Cam with anisotropic Gaussian processes

Context. We study astrometric residuals from a simultaneous fit of Hyper Suprime-Cam images. Aims. We aim to characterize these residuals and study the extent to which they are dominated by atmospheric contributions for bright sources. Methods. We used Gaussian process interpolation with a correlation function (kernel) measured from the data to smooth and correct the observed astrometric residual field. Results. We find that a Gaussian process interpolation with a von Kármán kernel allows us to reduce the covariances of astrometric residuals for nearby sources by about one order of magnitude, from 30 mas 2 to 3 mas 2 at angular scales of ~1 arcmin. This also allows us to halve the rms residuals. Those reductions using Gaussian process interpolation are similar to recent result published with the Dark Energy Survey dataset. We are then able to detect the small static astrometric residuals due to the Hyper Suprime-Cam sensors effects. We discuss how the Gaussian process interpolation of astrometric residuals impacts galaxy shape measurements, particularly in the context of cosmic shear analyses at the Rubin Observatory Legacy Survey of Space and Time.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Bipolar Membranes for Ion Management in (Photo)Electrochemical Energy Conversion

Conspectus (Photo)electrochemical energy conversion is important in the development of a carbon-neutral energy economy because it can provide a pathway for mitigating the intermittency of renewable energy sources such as wind and solar. In order to operate efficiently, these technologies, which include photoelectrochemical cells, water and CO 2 electrolyzers, fuel cells, and redox flow batteries, require fast charge transfer kinetics at the electrode/electrolyte interface as well as robust ion management in the electrolyte. In conventional electrolyzers and fuel cells, the electrolyte is strongly acidic or basic and ionic current is carried by H + or OH – ions. In contrast, photoelectrodes and electrocatalysts for water splitting are often studied in buffered solutions. The question of ion balance in these systems led us to analyze the polarization losses due to ion concentration gradients in cells that employed various buffer–membrane combinations. Continuously driving the buffer ions across an ionomer membrane not only lowers the buffer capacity of an aqueous electrolyte but also introduces pH gradients that result in significant energy losses. To address the problem, we and other groups have studied the use of reverse-biased bipolar membranes (BPMs) in (photo)electrolytic cells. BPMs consist of an anion exchange layer (AEL) laminated with a cation exchange layer (CEL) and are usually equipped with a catalytic layer in between to accelerate the water dissociation reaction. At the AEL/CEL interface, water dissociates into protons and hydroxide ions, which replenish those consumed at the cathode and anode. Compared to conventional water electrolyzers with proton/anion exchange membranes (PEM/AEM), BPM electrolyzers provide the unique advantage of continuously operating the cathode and anode under different pH conditions, which is desirable when the two electrode reactions have different pH requirements. BPMs also enable the use of buffered electrolytes at pH values that are optimized for electrode stability and product selectivity in applications such as CO 2 electrolysis. Product crossover losses and CO 2 pumping can be dramatically reduced in BPM-based CO 2 electrolyzers, relative to conventional alkaline membranes, by electrostatic repulsion (of anionic products) and electroosmotic drag (for neutral products). BPM-based gas fed CO 2 electrolyzers can achieve high current density, but they suffer from low Faradaic efficiency (FE) due to the acidic local environment of the CEL. This problem can be mitigated by adding an aqueous buffering layer or by creating a weak acid cation exchange film on the CEL face of the membrane. The use of BPMs in fuel cells and redox flow batteries offers some interesting advantages. Configurations with both reverse and forward bias have been studied, but forward bias has been favored due to material compatibility, reaction kinetics, and thermodynamic considerations. The net reaction at the AEL/CEL interface is the acid–base neutralization reaction, which has a high inherent reaction rate constant, but in the BPM is limited to a nanometric space-charge layer and requires efficient catalysis to achieve high current density. Understanding the mechanism of the acid–base neutralization and the opposite process, the water dissociation reaction, will be essential for improving the performance of forward-biased BPMs. In conclusion, this Account reviews our current understanding of the working mechanisms of BPMs and discusses how we can use them to effectively manage ions for various (photo)electrochemical applications.

Yan, Zhifei↗
DOE OSTI2026

Epitaxy of Beta-Ga2O3 on Highly-Offcut (Greater than 10 Degrees) Substrates

Beta-Ga2O3 has emerged as a leading candidate for next-generation power electronics, radio frequency (RF) switches, and extreme environment electronics due to a wide band gap (4.6 - 4.9 eV), high dopability (approximately 40 meV activation energy for an isolated silicon donor), and melt growth characteristics resulting in commercially available 4-inch substrates and commercial demonstrations of 6-inch substrates by multiple techniques. The (100) surface of Ga2O3 is highly desirable from a device and epitaxy standpoint - bulk growth of (100) material is more scalable than (010), the surface is nearly lattice-matched to p-type partner NiO, and Al2O3 incorporates at higher concentrations without phase separation. However, the epitaxial growth rate on (100) surfaces is less than 10% of other faces due to weak bonding and favorable desorption. Recent demonstrations have shown growth rate improvements from 0.4 nm/min to 1.5 nm/min by growing on (100) wafers that are offcut 6 degrees in the -c direction.1 These films show step-flow growth from (-201) step-edges and high electron mobility. Despite these exciting results, offcuts greater than 6 degrees have not been explored due to the waste associated with grinding and polishing large offcuts. In this talk we will discuss the molecular beam epitaxy (MBE) growth and properties of Beta-Ga2O3 grown on (100) substrates offcut in the -c direction up to 13.4 degrees. These large offcuts are enabled by edge-fed film-defined growth (EFG) where the offcut is grown into the surface by pulling the crystal through the EFG die with the seed crystal rotated by the desired offcut angle. We will demonstrate that 13.4 degrees offcut substrates still exhibit a terraced (100) surface, and that a >10x increase (4.8 nm/min) in growth rate is achieved. As previously reported on lower offcuts, we observe 100% reversal of substrate twin domains around the (001) direction at the substrate-epilayer interface. We will discuss electrical properties including record-low (by MBE) unintentional doping densities of < 5E15 cm-3.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Ferrotorryweiserite, Rh 5 Fe 10 S 16 , a New Mineral Species from the Sisim Placer Zone, Eastern Sayans, Russia, and the Torryweiserite–Ferrotorryweiserite Series

Ferrotorryweiserite, Rh 5 Fe 10 S 16 , occurs as small grains (≤20 µm) among droplet-like inclusions (up to 50 μm in diameter) of platinum-group minerals (PGM), in association with oberthürite or Rh-bearing pentlandite, laurite, and a Pt-Pd-Fe alloy (likely isoferroplatinum and Fe-Pd-enriched platinum), hosted by placer grains of Os-Ir alloy (≤0.5 mm) in the River Ko deposit. The latter is a part of the Sisim placer zone, which is likely derived from ultramafic units of the Lysanskiy layered complex, southern Krasnoyarskiy kray, Russia. The mineral is opaque, gray to brownish gray in reflected light, very weakly bireflectant, not pleochroic to weakly pleochroic (grayish to light brown tints), and weakly anisotropic. The calculated density is 5.93 g·cm –3 . Mean results (and ranges) of four WDS analyses are: Ir 18.68 (15.55–21.96), Rh 18.34 (16.32–20.32), Pt 0.64 (0.19–1.14), Ru 0.03 (0.00–0.13), Os 0.07 (0.02–0.17), Fe 14.14 (13.63–14.64), Ni 13.63 (12.58–14.66), Cu 4.97 (3.42–6.41), Co 0.09 (0.07–0.11), S 29.06 (28.48–29.44), and total 99.66 wt.%. They correspond to the following formula calculated for a total of 31 atoms per formula unit: (Rh 3.16 Ir 1.72 Pt 0.06 Ru 0.01 Os 0.01 ) Σ 4.95 (Fe 4.48 Ni 4.11 Cu 1.38 Co 0.03 ) Σ10.00 S 16.05 . The results of synchrotron micro-Laue diffraction studies indicate that ferrotorryweiserite is trigonal; its probable space group is $R\bar{3}m$ (#166) based on its Ni-analog, torryweiserite. The unit-cell parameters refined from 177 reflections are a = 7.069 (2) Å, c = 34.286 (11) Å, V = 1484 (1) Å 3 , and Z = 3. The c:a ratio is 4.8502. The strongest eight peaks in the X-ray diffraction pattern derived from results of micro-Laue diffraction study [d in Å(hkil)(I)] are 2.7950 ($20\bar{2}5$) (100); 5.7143 (0006) (60); 1.7671 ($22\bar{4}0$) (44.4); 3.0486 ($20\bar{2}1$) (39.4); 5.7650 ($10\bar{1}2$) (38.6); 2.5956 ($20\bar{2}7$) (37.8); 3.0058 ($11\bar{2}6$) (36.5); and 1.5029 ($4\bar{2}\bar{2}12$) (35.3). Ferrotorryweiserite and the associated PGM crystallized from microvolumes of residual melt at late stages of crystallization of grains of Os- and Ir-dominant alloys occurred in lode zones of chromitites of the Lysanskiy layered complex. In a particular case, the residual melt is disposed peripherally around a core containing a disequilibrium association of magnesian olivine (Fo 72.9–75.6 ) and albite (Ab 81.6–86.4 ), with the development of skeletal crystals of titaniferous augite: Wo 40.8–43.2 En 26.5–29.3 Fs 20.3–22.6 Aeg 6.9–9.5 (2.82–3.12 wt.% TiO 2 ). Ferrotorryweiserite represents the Fe-dominant analog of torryweiserite. We also report occurrences of ferrotorryweiserite in the Marathon deposit, Coldwell Complex, Ontario, Canada, and infer the existence of the torryweiserite–ferrotorryweiserite solid solution in other deposits and complexes.

58 GEOSCIENCES↗
DOE OSTI2026

Fast Growth of Ga2O3 on Highly Offcut (100)-Oriented Substrates

Beta-Ga2O3 has emerged as a leading candidate for next-generation power electronics, radio frequency (RF) switches, and extreme environment electronics due to a wide band gap (4.6 - 4.9 eV), high dopability (approximately 40 meV activation energy for an isolated silicon donor), and melt growth characteristics resulting in commercially available 4-inch substrates and commercial demonstrations of 6-inch substrates by multiple techniques. The (100) surface of Ga2O3 is highly desirable from a device and epitaxy standpoint - bulk growth of (100) material is more scalable than (010), the surface is nearly lattice-matched to p-type partner NiO, and Al2O3 incorporates at higher concentrations without phase separation. More importantly, the impact ionization coefficients along the [100] direction are low, leading to the highest possible critical fields. This is advantageous compared to the current state of the art, (001), due to reduced surface defects and increased possible breakdown voltage. However, the epitaxial growth rate on (100) surfaces is less than 10% of other faces due to weak bonding and favorable desorption, and on-axis (100) growth easily forms twin domains. Recent demonstrations have shown growth rate improvements from 0.4 nm/min to 1.5 nm/min by growing on (100) wafers that are offcut 6 degrees in the -c direction. These films show step-flow growth from (20-1) step-edges and high electron mobility due to suppressed twins. Despite these exciting results, offcuts greater than 6 degrees have not been explored due to the waste associated with grinding and polishing large offcuts. In this talk we will discuss the molecular beam epitaxy (MBE) growth and properties of Beta-Ga2O3 grown on (100) substrates offcut in the -c direction up to 13.4 degrees. These large offcuts are enabled by edge-fed film-defined growth (EFG) where the offcut is grown into the surface by pulling the crystal through the EFG die with the seed crystal rotated by the desired offcut angle. We will demonstrate that 13.4 degrees offcut substrates still exhibit a terraced (100) surface, and that a >10x increase (>5 nm/min) in growth rate is achieved. As previously reported on lower offcuts, we observe reversal of substrate twin domains around the (001) direction at the substrate-epilayer interface. We will discuss electrical properties including record-low (by MBE) unintentional doping densities of < 5E15 cm-3 and critical breakdown field in Schottky barrier diodes comparable with state-of-the-art (001) Ga2O3 without edge termination.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Adaptive Protection and Validated Models to Enable Deployment of High Penetrations of Solar PV (PV-MOD)

The availability and validation of various PV models in commercial tools differ, with some models not yet thoroughly validated for advanced inverter functionalities and reliable performance under weak system conditions. Many existing models do not fully incorporate new inverter control functions, which can affect system stability. The increasing deployment of solar PV and other inverter-based resources (IBRs), including distributed energy resources (DERs), is influencing the reliable operation of protection schemes in distribution systems and microgrids. Emerging adaptive protection schemes (APS) offer new opportunities for protecting these systems during varying configurations and DER operating conditions, though their demonstration and validation remain limited. Adaptive protection schemes face similar challenges, as they are typically designed for specific configurations. There is a growing need for tools and methodologies to streamline the deployment of adaptive protection for safe and reliable DER integration. The project main objective was to develop and validate high-fidelity generic models of solar PV facilities for stability, protection, EMT, and QSTS analyses. This objective was achieved, and these models can now be integrated into commercial software tools, enabling utilities, vendors, and developers to study high-penetration PV systems more confidently. The project also demonstrated advanced applications of these models, including the design and deployment of adaptive protection schemes in high-penetration field applications and microgrids, supporting grid safety and reliability. Several milestones were reached by the end of the project. A sophisticated inverter test plan was developed, and inverters representative of the North American marketplace were selected. EPRI and NREL tested various inverters, conforming to IEEE standards. Improvements were made to existing generic models of IBR units, IBR plants, and aggregated feeders for various analyses. The first generic electromagnetic transient (EMT) model for a solar PV plant was developed, conforming to IEEE Std 2800™-2022 and validated against laboratory measurements of a 2.2 MVA large-scale battery energy storage system (BESS) inverter. That model was then used to produce reference responses illustrating examples of validated and verified IBR plant models that pass or fail tests for technical minimum capability and performance as specified in the IEEE standard. The developed, tested, and validated generic models can be used for transmission planning, stability assessments, expansion planning, and evaluating potential future IBR interconnection requirements. They can also support interconnection screens and conformity assessments of IBR plants, including solar PV. The project significantly contributed to the ongoing standardization and model-based representation and verification of IBR responses. The project further addressed challenges of common distribution protection schemes with increasing deployment of DER by developing, validating, and demonstrating adaptive protection schemes (APS) that can improve the reliable and safe integration of DER into distribution systems. New APS were designed using improved DER models for three common distribution systems: a radial feeder, a meshed network, and a microgrid. Modeling and hardware-in-the-loop (HIL) testing of the APS were conducted, successfully showing their effectiveness and selectivity. Proof-of-concept field demonstration was achieved for two APS, i.e., one on a radial feeder and another one in a microgrid. Field demonstration could not be achieved for the APS on a meshed network, primarily due apprehension of one utility partner and also due to limited access to the protective algorithms in the network protectors. Guidelines developed from the lessons learned in the project lay out the general process followed in the design, installation, and commissioning of APS for various distribution systems. Distribution utility partners’ apprehension about field demonstration of the new APS were addressed—with varying success—by taking a stepped risk-management approach of modeling of a wide range of sensitivities first, performing in-depth proof-of-concept testing in the laboratory including HIL next, and finally deliberately implementing and commissioning the actual protection equipment and algorithms into parts of—or in parallel operation to—the three real distribution systems. Future work should include pilot projects that further show the acceptable performance of the developed APS before these schemes be rolled out more widely. Inclusion of both utility and original equipment manufacturers (OEMs) in future projects could increase chances of successful field demonstration. Despite challenges in achieving the field demonstration goal of the project for all three APS, the research significantly contributed to the innovation of adaptive protection solutions for scalable and reliable DER integration into distribution systems. This project significantly enhances the understanding of the impact of using appropriate inverter models on distribution and transmission (T&D) systems. By addressing the limitations of existing generic models, the project introduces high-fidelity models for stability, protection, electromagnetic transient (EMT), and quasi-static time series (QSTS) analyses. These models, integrated into commercial software tools, enable utilities, vendors, and developers to confidently study high-penetration PV systems. The project also demonstrates advanced applications, including adaptive protection schemes (APS) for distribution systems and microgrids, ensuring grid safety and reliability. The technical effectiveness and economic feasibility of the methods are evident through the development and validation of sophisticated inverter test plans and the selection of representative inverters. Testing by EPRI and NREL on retail, commercial, and utility-scale inverters, conforming to IEEE standards, underscores the robustness of the models. Improvements to existing generic models for various analyses further enhance their validity and applicability. The project also identifies gaps in common distribution protection schemes and designed new APS using improved DER models, demonstrating their effectiveness through modeling and hardware-in-the-loop (HIL) testing. The project’s benefits to the public are manifold. By advancing the standardization and model-based representation of IBR response, it supports transmission planning, stability assessments, and future IBR interconnection requirements. The generic models can facilitate better communication between transmission planners and developers, supporting expected IBR plant capability and performance. Additionally, the development of APS for radial feeders, meshed networks, and microgrids supports the integration of distributed energy resources (DERs) into distribution systems, enhancing grid reliability and safety. The project’s emphasis on thorough testing and simplicity in design ensures practical and scalable solutions for DER integration.

14 SOLAR ENERGY↗
DOE OSTI2025

Probing the Sustainable Reduction of CO 2 , N 2 and NO 3 to Fuels and Chemicals using Non-Traditional Porphyrinoid Catalysts

For the last several decades, numerous strategies have been proposed to be able to interconvert electricity and commodity fuels for transportation and other applications. Even as solar and wind energy become more widely available, these renewable energy sources are not always available where the population is most dense or during peak times of energy demand. One plausible solution to the issue of electricity storage and transport is to use electricity to drive the formation of high-energy compounds and chemical fuels. As a result, the electrochemical conversion of CO 2 into chemical fuels has received major attention as a multi-pronged approach for electricity coversion, storage, and transport. Similar strategies are also attractive for conversion of N 2 and NO x into value added compounds such as ammonia. Metalloporphyrin and metallocorrole complexes are comprised of aromatic tetrapyrrole scaffolds and have been extensively studied as homogeneous catalysts for critical catalytic processes like the electrochemical CO 2 reduction reaction (eCO 2 RR) to generate CO as feedstock for the Fischer–Tropsch process to generate a variety of hydrocarbons as chemical fuels. Despite their significance and impact, efficiently driving such catalytic reduction processes is challenging for multiple reasons. Such processes require multiple-electron transfer events, as opposed to single-electron redox chemistry that forms high-energy singly-reduced intermediates. In addition, these electron transfer events must be coupled to proton-transfer steps and avoid application of high overpotentials where proton (H + ) reduction to generate H 2 (as well as other side-reactions) can lower the selectivity and energy efficiency of the eCO 2 RR process. Certain Fe(III) aromatic tetrapyrroles (e.g., porphyrins and corroles) have been reported to overcome these challenges, but are often difficult to synthesize and modify, and have a tendency to support single-electron redox chemistry at the metal as opposed to multielectron redox involving the ligand and metal in concert. Other families of tetrapyrroles in which all four meso-carbons are reduced (i.e., sp 3 -hybridized) are also known in the literature. These non-aromatic tetrapyrroles are known as porphyrinogens and support multi-electron redox chemistry which has been elaborated by several groups, prompting researchers to consider whether these scaffolds may improve the kinetics and efficiencies of multi-electron steps attendant to activation of thermodynamically stable small-molecule substrates such as CO 2 , N 2 , NO 3 − , and NO 2 − . Although porphyrinogens support multi-electron redox properties, the four sp 3 -hybridized meso-carbons inherent to the tetrapyrrole core completely disrupt π-conjugation between pyrrolic units, which compromises the photochemical properties of porphyrinogens in comparison to that of traditional porphyrinoids. Moreover, the highly reducing nature of porphyrinogens predisposes these platforms to extreme air and water sensitivity. Accordingly, such metalated porphyrinogens are often pyrophoric and/or incompatible with many common organic solvents (CH 2 Cl 2 , CHCl 3 , CH 3 CN, EtOAc, etc.). Hence, reports of efficient small-molecule activation or catalysis supported by porphyrinogens have been limited. To overcome the instability of porphyrinogens while retaining the scaffold’s attractive multi-electron redox properties, we have directed attention to establishing less-well studied groups of non-aromatic tetrapyrroles (i.e., phlorins, biladienes, and isocorroles). These non-traditional tetrapyrroles each contain just a single sp 3 -hybridized meso-carbon, which ablates the platforms aromaticity but still provides an extended π-framework. In addition to developing innovative platforms that may be used to sustainably generate value-added chemicals, fuels, and ammonia via either photo- or electrocatalytic approaches our work also has significantly broadened our understanding of how to prepare and modulate the properties of non-aromatic tetrapyrroles for other applications. Our efforts entailed a synergistic partnership and active collaborations with Drs. David C. Grills and Mehmed Z. Ertem (both of Brookhaven National Lab) to better characterize isocorroles (and related non-aromatic tetrapyrroles) using a variety of advanced spectroscopic and theoretical methods. Our combined efforts have shown that the combination of properties that isocorroles (and related tetrapyrroles) provide results in good chemical stability paired with unique redox and photochemical characteristics that are not typically supported by simple and/or unadorned aromatic tetrapyrroles. Beside redox chemistry, the photochemistry of isocorroles and other non-aromatic tetrapyrroles containing one sp 3 -hybridized meso-carbon is also appealing due to their absorption in the long-visible to near-IR regions, which are essential for photocatalysis and related applications that benefit from direct excitation at these wavelengths. The results reported vastly improve our understanding of non-aromatic tetrapyrrole synthesis, properties and utility for activation of small molecule substrates such as O 2 and CO 2 in the presence of weak cationic organic acids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗