Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density correlation”

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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Strain Recovery and Defect Characterization in Mg‐Implanted Homoepitaxial GaN on High‐Quality GaN Substrates

The evolution of defects due to high‐pressure annealing of magnesium ion‐implanted epitaxial GaN grown on high‐quality GaN substrates is investigated. Changes in the implant‐induced strain are quantified as a function of annealing temperature and time. After annealing at 1300 °C for 10 min, the implant‐induced strain is fully relieved and accompanied by the presence of extended defects such as basal plane stacking faults and prismatic loops. Approximately one‐third of the original implant‐induced strain remains after annealing at 700 °C, and 5% of the original strain remains at 1000 °C for 100 min. In all cases, nearly all of the recovered strain occurs within first few minutes of annealing. A prominent increase in the asymmetric (104) triple axis X‐ray rocking curve full width at 0.01 maximum (FW0.01M) is observed after annealing at 1300 °C for 10 min. After annealing at 1300 °C for 100 min, a subsequent decrease in FW0.01M is correlated with a reduction of the extended defect density from 4 × 10 8 to 3 × 10 7 cm −2 , determined through transmission electron microscope (TEM) measurements. Further reduction in the density of the extended defects by optimizing annealing temperature and time is expected to improve the performance of GaN‐based vertical power devices.

Wang, Yekan↗
DOE OSTI2022

Correlating dislocation structures to basal and prismatic slip bands near grain boundaries in tensile-tested Mg–4Al using a multiscale electron microscopy approach

Dislocation structures where basal and prismatic slip bands meet grain boundaries in tensile-tested Mg with 4 wt% Al (Mg–4Al) were studied using a multiscale electron microscopy approach to explore the assumptions made for the dislocation pileup theory: a) slip bands consist of dislocation arrays piling up near grain boundaries, and b) stress concentration in the adjacent grain is solely caused by dislocation pileup. After post-testing electron backscatter diffraction (EBSD) scans, focused-ion-beam (FIB) lift-out specimens were prepared from regions of interest so that the specimen plane is parallel to the bulk sample surface to allow for direct correlation between mesoscale and microscale characterization. High dislocation density and the dislocation arrays in the slip bands corroborated with the first assumption. Evidence of plastic deformation in all grains, however, showed that the second assumption was false. Additionally, dislocation structures including crystal rotation boundaries and dislocations dissociating to the basal plane provide additional stress relief mechanisms near grain boundaries that are not accounted for in many crystal plasticity models, which requires auxiliary models to more accurately predict the stress and strain distribution in a microstructure.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Third-Order Møller–Plesset Theory Made More Useful? The Role of Density Functional Theory Orbitals

The practical utility of Møller–Plesset (MP) perturbation theory is severely constrained by the use of Hartree–Fock (HF) orbitals. It has recently been shown that the use of regularized orbital-optimized MP2 orbitals and scaling of MP3 energy could lead to a significant reduction in MP3 error. In this work, we examine whether density functional theory (DFT)-optimized orbitals can be similarly employed to improve the performance of MP theory at both the MP2 and MP3 levels. We find that the use of DFT orbitals leads to significantly improved performance for prediction of thermochemistry, barrier heights, noncovalent interactions, and dipole moments relative to the standard HF-based MP theory. Indeed, MP3 (with or without scaling) with DFT orbitals is found to surpass the accuracy of coupled-cluster singles and doubles (CCSD) for several data sets. We also found that the results are not particularly functional sensitive in most cases (although range-separated hybrid functionals with low delocalization error perform the best). As such, MP3 based on DFT orbitals thus appears to be an efficient, noniterative O(N 6 ) scaling wave-function approach for single-reference electronic structure computations. Scaled MP2 with DFT orbitals is also found to be quite accurate in many cases, although modern double hybrid functionals are likely to be considerably more accurate.

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

Electronic Structure of the Primary Electron Donor P 700 +• in Photosystem I Studied by Multifrequency HYSCORE Spectroscopy at X- and Q-Band

The primary electron donor P 700 of the photosystem I ( PS I) is a heterodimer consisting of two chlorophyll molecules. A series of electron-transfer events immediately following the initial light excitation leads to a stabilization of the positive charge by its cation radical form, P 700 +• . Here, the electronic structure of P 700 +• and, in particular, its asymmetry with respect to the two chlorophyll monomers is of fundamental interest and is not fully understood up to this date. Here, we apply multifrequency X- (9 GHz) and Q-band (35 GHz) hyperfine sublevel correlation (HYSCORE) spectroscopy to investigate the electron spin density distribution in the cation radical P 700 +• of PS I from a thermophilic cyanobacterium Thermosynechococcus elongatus . Six 14 N and two 1 H distinct nuclei have been resolved in the HYSCORE spectra and parameters of the corresponding nuclear hyperfine and quadrupolar hyperfine interactions were obtained by combining the analysis of HYSCORE spectral features with direct numerical simulations. Based on a close similarity of the nuclear quadrupole tensor parameters, all of the resolved 14 N nuclei were assigned to six out of total eight available pyrrole ring nitrogen atoms (i.e., four in each of the chlorophylls), providing direct evidence of spin density delocalization over the both monomers in the heterodimer. Using the obtained experimental values of the 14 N electron-nuclear hyperfine interaction parameters, the upper limit of the electron spin density asymmetry parameter is estimated as R A/B upper = 7.7 ± 0.5, while a tentative assignment of 14 N observed in the HYSCORE spectra yields R B/A = 3.1 ± 0.5.

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

Insights into Nonelectroactive C–C Bond Formation on Cu(100) during Electrochemical CO 2 Reduction from Multiconfigurational Wavefunction Theory

Carbon–carbon (C–C) bond formation is necessary for hydrocarbon (and oxygenate) synthesis beyond methane (and formate/formic acid) during electrochemical CO and CO 2 reduction (ECOR and ECO 2 R). Cu has notable ability to form hydrocarbons compared to other pure metals. In particular, the (100) facet of face-centered cubic Cu forms ethylene competitively with H 2 and methane during both ECOR and ECO 2 R. Past simulations based on density functional theory (DFT) with standard exchange-correlation functional approximations predict fast nonelectroactive C–C bond formation channels involving adsorbed (*) CO together with another *CO, formyl (*CHO), or hydroxymethylidyne (*COH), forming OC*–*CO, OC*–CHO*, and OC*–*COH, respectively. Such simulations support the prevailing hypothesis that emergence of C 2 products is kinetically determined at the early stages of the reduction chemistry. Here we show, via simulations with more accurate many-body, i.e., “correlated”, wavefunction theory (enabled by an embedding scheme), that the coupling of *CO with a *CO or a *COH (previously predicted at the same level of theory to kinetically dominate over *CHO as the one-electron reduction product of *CO) is highly activated (kinetically impeded), with free energy barriers >1 eV, in contradiction to previous DFT-based simulations. Intriguingly, we find that the coupling of two adjacent *COHs incurs only a small barrier (<0.3 eV) and is exoergic (< –1 eV); however, given the predicted low surface mobility of *COH, the emergence of HOC*–*COH is also improbable, at least at low *COH coverages. We therefore conclude that it is highly unlikely for *CO to participate in nonelectroactive C–C bond formation on pristine Cu(100), contrary to conventional wisdom, and that the energetically favorable *COH dimerization may occur only after substantial buildup of *COH on the surface.

adsorption↗
DOE OSTI2026

Photocatalytic Ammonia Synthesis using Fe-Based MOFs: The Role of Ligand Functionalization

Photocatalytic ammonia (NH 3 ) synthesis offers a carbon-neutral alternative to the Haber−Bosch process, which generates 42 million metric tons of CO 2 equivalent emissions annually. However, solar-to-ammonia conversion with contemporary photocatalysts remains far from practical requirements, and understanding the limiting factors in systems with well-defined active sites is crucial. Here, we show how the μ 3 -oxo-centered trinuclear Fe cluster in MIL-101(Fe) functions as the catalytic motif for N 2 -to-NH 3 conversion through combined experimental and computational investigations. Comparative studies with a molecular analogue demonstrate that the cluster is stabilized within the MOF framework, sustaining redox cycling and maintaining high catalytic activity. We systematically functionalized the dicarboxylate ligands of MIL-101(Fe) with −NH 2 , −Br, −NO 2 , −F, and −CF 3 to probe how ligand chemistry modulates Fe electron density, N 2 adsorption capacity, and proton availability, correlating these properties with catalytic performance using spectroscopic and surface characterization techniques alongside timeresolved infrared to assess excited-state lifetimes. F-functionalization optimally balances N 2 activation, proton availability at Fe active sites, and excited-state lifetimes, boosting NH 3 production by ∼ 60% relative to unmodified MIL-101(Fe). This study of ligandfunctionalized MIL-101(Fe) MOFs uncovers the underlying structure-activity relationships and advances design principles for solardriven NH 3 synthesis.

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

Tuning Solid Electrolyte Interphase Formation before Plating Onset in Anode-Free Sodium Batteries

Sodium (Na) batteries are of growing interest due to the higher earth abundance of sodium than lithium, as well as their promising theoretical energy density when metallic Na anodes are used. However, Na plating and stripping are heavily influenced by the physicochemical properties of the solid electrolyte interphase (SEI), which is directly influenced by the solvent and salt used for the electrolyte. While most studies focus on the SEI that forms on the surface of Na metal after plating, we expand this analysis by identifying a nanoscale “pre-plating” SEI that forms on the current collector (CC) prior to the onset of Na plating. Here, we systematically investigate an array of Na salt and glyme solvents in the electrolyte and determine the associated impacts on pre-plating SEI formation on aluminum CCs. By combining analytical electrochemistry approaches with a multimodal suite of spectroscopy techniques (X-ray, infrared, and Raman), supported by density functional theory calculations, we reveal a direct correlation between the Na + coordination environment and pre-plating SEI composition. We find that longer-chain glymes produce larger proportions of organic alkoxide products in the interphase, consistent with increased Na + −glyme interactions, while the fraction of salt-derived inorganic products (e.g., NaF) correlates with Na + −anion coordination. These insights highlight the critical influence of electrolyte composition particularly solvent identity and Na + coordination on the initial SEI formation in anode-free Na batteries.

anode-free batteries↗
DOE OSTI2021

Accurate simulation of surfaces and interfaces of ten FCC metals and steel using Lennard–Jones potentials

Abstract The earlier integration of validated Lennard–Jones (LJ) potentials for 8 fcc metals into materials and biomolecular force fields has advanced multiple research fields, for example, metal–electrolyte interfaces, recognition of biomolecules, colloidal assembly of metal nanostructures, alloys, and catalysis. Here we introduce 12-6 and 9-6 LJ parameters for classical all-atom simulations of 10 further fcc metals (Ac, Ca (α), Ce (γ), Es (β), Fe (γ), Ir, Rh, Sr (α), Th (α), Yb (β)) and stainless steel. The parameters reproduce lattice constants, surface energies, water interfacial energies, and interactions with (bio)organic molecules in 0.1 to 5% agreement with experiment, as well as qualitative mechanical properties under standard conditions. Deviations are reduced up to a factor of one hundred in comparison to earlier Lennard–Jones parameters, embedded atom models, and density functional theory. We also explain a quantitative correlation between atomization energies from experiments and surface energies that supports parameter development. The models are computationally very efficient and applicable to an exponential space of alloys. Compatibility with a wide range of force fields such as the Interface force field (IFF), AMBER, CHARMM, COMPASS, CVFF, DREIDING, OPLS-AA, and PCFF enables reliable simulations of nanostructures up to millions of atoms and microsecond time scales. User-friendly model building and input generation are available in the CHARMM-GUI Nanomaterial Modeler. As a limitation, deviations in mechanical properties vary and are comparable to DFT methods. We discuss the incorporation of reactivity and features of the electronic structure to expand the range of applications and further increase the accuracy.

Kanhaiya, Krishan (ORCID:0000000236222655)↗
DOE OSTI2025

Radiation damages the silicates present in polyphasic mineral aggregates causing concrete’s degradation

While many U.S. nuclear power plants have submitted Subsequent License Renewal Application to operate beyond 60 years, others are already considering Operations Beyond Eighty years. In such cases, concrete biological shields are exposed to neutron and gamma radiation exceeding prescribed thresholds. Radiation-induced volumetric expansion (RIVE), extensively studied in single crystals, may also contribute to the degradation of polycrystalline aggregates. Since minerals differ in atomic structure and chemistry, radiation can affect them in distinct ways. This study examines quartzite, marble, and limestone to evaluate how irradiation affects their physical attributes and chemical reactivity. Results show crystalline silicates experience significant RIVE damage and enhanced reactivity in alkaline solutions compared to non-irradiated phases. Enhanced intra- and inter-granular dissolution could compromise aggregate integrity. An empirical correlation links silicate dissolution rate to atomic constraints (density, rigidity) and radiation dose, providing a predictive framework for estimating changes in silicate aggregate properties within radiation-exposed concrete.

Bouissonnié, Arnaud [Univ. of California, Los Ange↗
DOE OSTI2020

Controlling protein adsorption modes electrostatically

Protein adsorption on surfaces is ubiquitous in biology and in biotechnology. There are various forces required for controlling protein adsorption. Here, we introduce an explicit ion coarse-grained molecular dynamics simulation approach for studying the effects of electrostatics on protein adsorption, and 2D protein assembly on charged surfaces. Our model accounts for the spatial distribution of protein charges. We use catalase as our model protein. We find that the preferential adsorption mode of proteins at low protein concentration on a charged surface is “standing up”. When the protein concentration in a solution increases to reach a critical density on the surface, the adsorption mode switches from “standing up” to a mixed state “flat on” and “standing up”, which increases the lateral correlations among the adsorbed proteins. As such, the changes in the adsorption mode arise from the protein adsorption that cancel the surface charge and the protein–protein repulsion. This correlated surface structure melts as the salt concentration increases because the charged surface is cancelled by the salt ions and the proteins de-adsorb. For the case of strongly charged surfaces the “standing up” conformation remains more favorable even at high protein adsorption at low salt concentrations since in that conformation the surface charge is cancelled more effectively, generating an even more laterally correlated structure. We elucidate the effects of parameters such as surface charge density, salt concentration, and protein charges on the different adsorption modes and the structure and organization of proteins on the charged surfaces. Finally, this study provides a guide for controlling protein assembly on surfaces.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Excited states in variational Monte Carlo using a penalty method

Here, in this article, the authors present a technique using variational Monte Carlo to solve for excited states of electronic systems. This technique is based on enforcing orthogonality to lower energy states, which results in a simple variational principle for the excited states. Energy optimization is then used to solve for the excited states. This technique is applied to the well-characterized benzene molecule, in which ~10 000 parameters are optimized for the first 12 excited states. Agreement within ~0.2 eV is obtained with higher scaling coupled cluster methods; small disagreements with experiment are likely due to vibrational effects.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Infrared spectroscopy of RG–Co + (H 2 O) complexes (RG = Ar, Ne, He): The role of rare gas “tag” atoms

RG n –Co + (H 2 O) cation complexes (RG = Ar, Ne, He) are generated in a supersonic expansion by pulsed laser vaporization. Here, complexes are mass-selected using a time-of-flight spectrometer and studied with infrared laser photodissociation spectroscopy, measuring the respective mass channels corresponding to the elimination of the rare gas “tag” atom. Spectral patterns and theory indicate that the structures of the ions with a single rare gas atom have this bound to the cobalt cation opposite the water moiety in a near-C 2v arrangement. The O–H stretch vibrations of the complex are shifted compared to those of water because of the metal cation charge-transfer interaction; these frequencies also vary systematically with the rare gas atom attached. The efficiencies of photodissociation also vary with the rare gas atoms because of their widely different binding energies to the cobalt cation. The spectrum of the argon complex could only be measured when at least three argon atoms were attached. In the case of the helium complex, the low binding energy allows the spectra to be measured for the low-frequency H–O–H scissors bending mode and for the O–D stretches of the deuterated analog. The partially resolved rotational structure for the antisymmetric O–H and O–D stretches reveals the temperature of these complexes (6 K) and establishes the electronic ground state. The helium complex has the same 3 B 1 ground state as the tag-free complex studied previously by Metz and co-workers (2013), but the A rotational constant is contaminated by vibrational averaging from the bending motion of the helium.

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

Structures of glasses created by multiple kinetic arrests

X-ray scattering has been used to characterize glassy itraconazole (ITZ) prepared by cooling at different rates. Faster cooling produces ITZ glasses with lower (or zero) smectic order with more sinusoidal density modulation, larger molecular spacing, and shorter lateral correlation between the rod-like molecules. We find that each glass is characterized by not one, but two fictive temperatures T f (the temperature at which a chosen order parameter is frozen in the equilibrium liquid). The higher T f is associated with the regularity of smectic layers and lateral packing, while the lower T f with the molecular spacings between and within smectic layers. This indicates that different structural features are frozen on different timescales. The two timescales for ITZ correspond to its two relaxation modes observed by dielectric spectroscopy: the slower δ mode (end-over-end rotation) is associated with the freezing of the regularity of molecular packing and the faster α mode (rotation about the long axis) with the freezing of the spacing between molecules. Finally, our finding suggests a way to selectively control the structural features of glasses.

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

High-throughput ab initio reaction mechanism exploration in the cloud with automated multi-reference validation

Quantum chemical calculations on atomistic systems have evolved into a standard approach to studying molecular matter. These calculations often involve a significant amount of manual input and expertise, although most of this effort could be automated, which would alleviate the need for expertise in software and hardware accessibility. Here, we present the AutoRXN workflow, an automated workflow for exploratory high-throughput electronic structure calculations of molecular systems, in which (i) density functional theory methods are exploited to deliver minimum and transition-state structures and corresponding energies and properties, (ii) coupled cluster calculations are then launched for optimized structures to provide more accurate energy and property estimates, and (iii) multi-reference diagnostics are evaluated to back check the coupled cluster results and subject them to automated multi-configurational calculations for potential multi-configurational cases. All calculations are carried out in a cloud environment and support massive computational campaigns. Key features of all components of the AutoRXN workflow are autonomy, stability, and minimum operator interference. We highlight the AutoRXN workflow with the example of an autonomous reaction mechanism exploration of the mode of action of a homogeneous catalyst for the asymmetric reduction of ketones.

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

The surface chemistry of colloidal lead halide perovskite nanowires

This study explored the interplay between the ligand–surface chemistry of colloidal CsPbBr 3 nanowires (NWs) and their optical properties. The ligand equilibrium was probed using nuclear magnetic resonance spectroscopy, and by perturbing the equilibrium via dilution, the gradual removal of ligands from the CsPbBr 3 surface was observed. This removal was correlated with an increase in the surface defect density, as suggested by a broadening of the photoluminescence (PL) spectrum, a decrease in the PL quantum yield (PLQY), and quenching of the PL decay. These results highlight similar surface binding between the traditional CsPbBr 3 quantum dots and our NWs, thereby expanding the scope of well-established ligand chemistry to a relatively unexplored nanocrystal morphology. In conclusion, by controlling the dilution factor, it was revealed that CsPbBr 3 NWs achieve a PLQY of 72% ± 2% and a relatively long average PL lifetime of 400 ± 10 ns, without relying on additional surface passivation techniques, such as ligand exchange.

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

Variation in supercurrent connectivity and vortex pinning of state-of-the-art pulsed laser deposited RE Ba 2 Cu 3 O 7- x coated conductors

REBa 2 Cu 3 O 7-x (REBCO, RE: rare earth) coated conductors (CCs) suffer from great critical current I c differences between different manufacturers, I c variations within individual manufacturers, and often significant lengthwise fluctuations. The understanding of such variations is complicated by the lack of a direct correlation between I c and the critical current density J c . In fact, although J c is the fundamental property determined by the local vortex pinning landscape, I c is often limited by variable current blocking mechanisms. An important practical complexity is commercial practices of performing J c and I c evaluations based on I c at 77 K and self-field (sf), where connectivity variation dominates over vortex pinning variations. However, at higher fields and lower temperatures, vortex pinning becomes more complex and highly variable, making predictions of I c and J c at arbitrary temperature T, magnetic field H, and field orientation θ, quite uncertain. To address some aspects of this problem, we conducted detailed spool-to-spool performance characterization on recently manufactured REBCO CCs. Despite J c (77 K, sf) varies by only ∼11%, J c (77 K, 1 T) of its minimum and maximum (for H//ab-plane) show variations of ∼21% and ∼32%, respectively. This emphasizes the shortcoming in using J c (77 K, sf) as parameter for evaluating even the low field performance. An even more remarkable spool-to-spool J c variation of ∼68% was observed at 20 K and 15 T for H//ab-plane. To identify the origin of such lack of reproducibility we performed microstructural characterizations, which revealed, within the REBCO layer, large variation in the density of copper oxide (CuO x ) particles ranging from 0.1 to 2 μm in size. We believe that they play a decisive role in reducing the effective cross-section of the REBCO layer by not simply blocking current themselves, but also by nucleating off-axis REBCO grains, whose misoriented grain boundaries adversely impact REBCO grain-to-grain connectivity. The REBCO growth associated with high density CuO x particles also leads to the more disordered spatial arrangement of BaHfO 3 precipitate arrays, which, when self-aligning along the ab-planes, generate stronger pinning enhancing J c (H//ab) at all temperatures. In this way, we established that variations of both connectivity and vortex pinning are thus directly coupled. Our results also explain why the so-called ‘lift-factor’, typically defined by the ratio between I c (T, H) and I c (77 K, sf), frequently turns out to be unreliable.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Adaptive machine learning for time-varying systems: low dimensional latent space tuning

Machine learning (ML) tools such as encoder-decoder convolutional neural networks (CNN) can represent incredibly complex nonlinear functions which map between combinations of images and scalars. For example, CNNs can be used to map combinations of accelerator parameters and images which are 2D projections of the 6D phase space distributions of charged particle beams as they are transported between various particle accelerator locations. Despite their strengths, applying ML to time-varying systems, or systems with shifting distributions, is an open problem, especially for large systems for which collecting new data for re-training is impractical or interrupts operations. Particle accelerators are one example of large time-varying systems for which collecting detailed training data requires lengthy dedicated beam measurements which may no longer be available during regular operations. We present a novel method of adaptive ML for time-varying systems. Our approach is to map very high (N ≈ 100k) dimensional inputs (a combination of scalar parameters and images) into the low dimensional (N ≈ 2) latent space at the output of the encoder section of an encoder-decoder CNN. We then actively tune the low dimensional latent space-based representation of complex system dynamics by the addition of an adaptively tuned feedback vector directly before the decoder sections builds back up to our image-based high-dimensional phase space density representations. This method allows us to learn correlations within and to quickly tune the characteristics of incredibly large parameter space systems and to track their evolution in real time based on feedback without massive new data sets for re-training. We demonstrate that our method can accurately predict and track the phase space of charged particle beams at various locations in a particle accelerator by adaptively adjusting in real-time while the unknown input beam distribution of the accelerator is changing in shape, charge, and offset and while the RF system of the accelerator itself is also changing in an unpredictable way. For FACET-II we demonstrate that such an approach has the potential to use transverse deflecting cavity and energy spread spectrum beam measurements to accurately predict 2D projections of the 6D phase space of the electron beam at the plasma wakefield acceleration interaction point where such diagnostics are unavailable.

47 OTHER INSTRUMENTATION↗