Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “KS”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Growth of Hexagonal Boron Nitride from Molten Nickel Solutions: A Reactive Molecular Dynamics Study

Metal flux methods are excellent for synthesizing high-quality hexagonal boron nitride (hBN) crystals, but the atomic mechanisms of hBN nucleation and growth in these systems are poorly understood and difficult to probe experimentally. Here, we harness classical reactive molecular dynamics (ReaxFF) to unravel the mechanisms of hBN synthesis from liquid nickel solvent over time scales up to 30 ns. These simulations mimic experimental conditions by including relatively large liquid nickel slabs containing dissolved boron and a molecular nitrogen gas phase. Overall, the reaction takes place almost exclusively on the surface of the liquid nickel, owing to the low solubility of nitrogen in bulk nickel and the intermediate species’ preference for the metal–gas interface. The formation of hBN invariably begins by reaction of dinitrogen with nickel-solvated boron atoms at the surface, forming intermediate N–N–B species, which typically evolve into B–N–B units through a short-lived intermediate where a single nitrogen atom is coordinated by one nitrogen and two boron atoms. The resulting B–N–B units, in turn, coalesce with growing hBN nuclei and carry nitrogen between hBN nanocrystals in an Ostwald ripening process. The amount of hBN produced on the tens of nanosecond time scale depends critically on the boron concentration, while having a much weaker dependence on the N 2 pressure for the regime considered (N 2 pressures of 2.5–10 MPa, Ni–B solutions with 6–12% boron by atom fraction). The highest rate of hBN formation occurs at the lowest temperature considered (1750 K, just above the melting point of nickel), while no hBN sheets are formed at 2000 K or above. An analysis of the transition pathways for nitrogen atoms shows that the final step, incorporation of small B–N motifs into larger hBN sheets, is the rate-limiting step in the regimes considered. While raising the temperature from 1750 to 2000 K has little effect on the formation of intermediates (N–N–B, B–N–B, etc.), the lack of large hBN sheets at temperatures >1900 K is explained by decreased probability of the final step and increased probability of breakup of hBN into B–N motifs.

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

Controlling rotational air lasing lineshape by carrier-envelope offset phase

The carrier-envelope offset phase (CEP) of a few-cycle optical pulse is commonly used to control electron dynamics on the attosecond timescale, whereas lasing spectra from transitions between rotational states are generally emitted over much longer durations, typically nanoseconds. Here, we demonstrate CEP control of the rotational lasing spectra corresponding to the transition from B 2 Σ$^+_u$ to X 2 Σ$^+_g$ (1) in N$^+_2$ cations, transforming its lineshape from a symmetric Lorentzian profile to an asymmetric Fano type–and vice versa. This lineshape modulation arises from the interference between the B-X coherence initiated by the main pulse and the supercontinuum (self seed) by self-phase modulation, resembling an “f-to-3f" interferometry. Additionally, for lasing lines with lower rotational quantum numbers, we observe a stronger coupling between adjacent lasing peaks, which originates from the amplification of both even- and odd-order rotational coherent emission lines. Our study presents a general framework for controlling lasing lineshapes and provides new insights into sub-optical-cycle dynamics in air lasing.

atomic and molecular interactions with photons↗
DOE OSTI2025

Light control of intramolecular nuclear dynamics by vortex electron localization

In strong-field ionization of molecules, intense light pulses are thought to have a negligible direct influence on atomic nuclei. Molecular dissociation is thus expected to be determined by the geometrical configuration of the molecular ion at the ionization instant. Contrary to this picture, we observe a counterintuitive electron-proton angular correlation and the formation of proton vortices following strong-field ionization of H 2 molecules by bicircularly polarized two-color laser fields. We explain this phenomenon by the pathway interference and localization of the residual H 2 ⁡+ electron in different angular-momentum states formed in the tail of the driving laser pulse. We validate this interpretation by combining a quantum-mechanical numerical simulation of the field-driven coupled electronic-nuclear dynamics and a semiclassical-trajectory model for the phase accumulation of the laser-driven electronic-nuclear wave packet. Our joint experimental and theoretical study reveals a general picture of vortex electron localization which can be used for controlling molecular-bond breaking with circularly polarized laser fields.

Atomic & molecular processes in external fields↗
DOE OSTI2025

Tracking the spatial and temporal links between late Cenozoic extension and magmatism in the Death Valley region (California, USA)

The relationship between late Cenozoic magmatism and extension in the central Basin and Range province (western United States) is complex, necessitating high-precision geochronology to understand its spatiotemporal connections. In the Death Valley region (California), the lack of high-precision U-Pb zircon ages has limited our understanding of the timing of pluton formation and its links to regional extension. We present new high-precision chemical abrasion-isotope dilution-thermal ionization mass spectrometry 206 Pb/ 238 U zircon ages and trace element analyses for eight Death Valley plutons. Our findings reveal three distinct phases of intrusive magmatism: (1) emplacement of shallow rapakivi granites at 13.2 Ma, (2) construction of the mid-crustal Black Mountains intrusive complex at 11.3 Ma, and (3) late emplacement of shallow, compositionally diverse intrusions at 8.2 Ma. A gap in zircon crystallization between 10 Ma and 8.2 Ma coincides with exhumation of the Black Mountains and a transition from sill to dike emplacement. The dominance of rapakivi granites in the Death Valley region, which is rare among Cenozoic granitoids, is likely a result of rapid crustal extension that induces adiabatic decompression. A comparison of the timing of volcanism, plutonism, and tectonic events in Death Valley reveals that intrusive magmatism closely tracks the locus of extension, underscoring the plutonic record as a vital link for understanding regional tectonics and changes in plate boundary dynamics during this period.

58 GEOSCIENCES↗
DOE OSTI2026

JWST Observations of Starbursts: Dust Processing in the M82 Superwind

We present JWST MIRI and NIRCam imaging of the inner ∼5 kpc of the M82 superwind at ~ $0''_{.}05 - 0''_{.}375$(∼0.9–6.5 pc) resolution. Targeted filters probe emission from polycyclic aromatic hydrocarbons (PAHs; F335M, F360M, F770W, F1130W) and continuum (F250M, F360M) The images reveal a network of cool wind filaments traced by PAHs. PAH surface brightness declines with the inverse square of distance to the midplane, suggesting that the incident radiation field from the starburst drives the observed PAH intensity out to ±2.5 kpc. The 3.3/11.3 and 3.3/7.7 μm band ratios show uniformity with distance from the starburst, though comparisons with mid-IR dust emission models indicate a modest shift toward larger PAHs. Outside the disk, 11.3/7.7 μm increases moderately, reflecting that PAHs become more neutral with distance from the starburst as they are exposed to a declining radiation field and ionization parameter. Overall, PAHs in the wind are consistent with standard-to-large sizes and standard-to-high ionization states. Including Spitzer and Herschel data, PAH abundance (q PAH ) is set at ∼1% in the starburst and remains unchanging out to ±5 kpc off the disk. This flat q PAH profile suggests that PAHs are shielded from the hot wind, perhaps residing in the surface layers of cool clouds, with possible replenishment from cloud interiors and enrichment of the halo from previous bursts. In this picture, clouds are not dense enough to promote PAH growth, and they likely undergo radiative cooling and mixing with the hot phase to survive the gauntlet for at least ∼20 Myr.

Cronin, Serena A. [University of Maryland, College↗
DOE OSTI2026

Comparison of structurally diverse simulation models for prediction of epidemic outcomes caused by a long-distance dispersed pathogen

Long-distance dispersal (LDD) pathogens pose substantial challenges for epidemic control due to their ability to generate new infection foci at great distances. While various modeling approaches have been developed to understand and manage such outbreaks, little work has compared how models of different structures behave under shared conditions. Here, in this study, we compare four structurally distinct epidemiological models — EPIMUL, GEMF, PoPS, and Warwick — each adapted to simulate the spread of wheat stripe rust (WSR), a wind-dispersed LDD pathogen, under identical epidemiological parameters and dispersal kernel. Using data from a controlled field experiment, we evaluate the ability of each model to replicate disease prevalence under nine intervention scenarios that vary in timing and culling area. While the models differ substantially in design — ranging from spatial grid-based to network-based and raster-based frameworks — the shared dispersal kernel allowed for close alignment in their predictions. All models accurately captured general epidemic trends, particularly the strong effect of early intervention on disease suppression. We qualitatively compared their behavioral responses across scenarios and also evaluated an ensemble prediction by averaging across model outputs. Our findings highlight how integrating shared epidemiological components into distinct modeling frameworks can improve consistency and accuracy, while reinforcing the importance of early culling in managing LDD pathogen outbreaks.

Dispersal kernel↗
DOE OSTI2025

Thermal performance and energy consumption validation of an occupied local government office building outfitted with ceiling tile phase change materials

Buildings present an opportunity for energy conservation and the modulation of peak energy demand through controlled Heating, Ventilation, and Air Conditioning (HVAC) energy use. The administrative and office building stock in the United States holds potential to achieve energy and demand savings through retrofits such as insulation, weatherization, and thermal energy storage. Specifically, there is a need to validate passive phase change material (PCM) applications in full scale in aging administrative buildings in the US to evaluate the energy benefits. Aim of this study was to conduct a whole building level thermal and energy validation of an operational building and explore an alternative method for evaluating energy efficiency. To accomplish this, the study employed PCMs in the drop ceiling and carry out an energy audit and on-site measurement of HVAC systems' energy demand and consumption. A full-scale EnergyPlus energy model, modeled by the authors, served as a baseline for evaluation. The results show that calibrated model's envelope temperature measures fall within the accepted errors. HVAC energy simulation results also fall within the accepted errors for monthly and hourly pre- and post- PCM retrofit electricity and natural gas data. The novelty of this study is that it employees energy scales per Heating Degree Hour and Cooling Degree Hour, in contrast to the commonly used Heating Degree Days and Cooling Degree Days as reported in the literature to analyze energy savings. These findings underscore the pivotal role of a calibrated model in assessing the efficacy of a singular energy measure, like a PCM-retrofitted ceiling, in an occupied office building.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

ED-cPSD: Fast Phase-Size Distribution via Sequential Erosion-Dilation

The Erosion-Dilation continuous Phase-Size Distribution, ED-cPSD, is an application for calculating continuous pore and particle-size distribution from digital reconstructions and/or image-based structural data. It is based on the erosion-dilation continuous phase-size distribution method. A continuous size distribution is a measure of the probability density of finding a particle or pore of a certain size. These distributions are of interest in any field of study involving porous media, including but not limited to electrochemistry, petroleum engineering, geology, and food science. The algorithm behind the software provides a computationally efficient way to calculate phase-size distributions for large domains. For a 3D battery electrode reconstruction with 1.3 x 10 8 voxels, the particle size distribution is derived in under 2 min on a desktop, while also retaining flexibility and computational efficiency for HPC-scale multi-threading. The software can handle structures with over 10 9 voxels. The algorithm is roughly 280 times faster than a previous version on the same task.

Characterization↗
DOE OSTI2025

Autocorrelation and Multifractal Detrended Fluctuation Analyses Reveal Superdiffusive Mass Transport in Solvent-Filled Nanoporous Media

Fluorescence fluctuation spectroscopy experiments were conducted to better understand the complex mass transport dynamics of organic molecules in liquid-filled nanoporous media. Anodic aluminum oxide (AAO) membranes incorporating 10 and 20 nm diameter cylindrical pores were employed as model materials. Nile red (NR) dye was used as a fluorescent tracer. The dye was dissolved separately in ethanol and toluene at a concentration of 20 nM and used to fill the membrane nanopores. Confocal fluorescence microscopy was employed to capture photon intensity time series data reflecting apparent diffusion of the dye within the pores. Autocorrelation of these data revealed that NR diffusion within the membranes occurred over a broad range of time scales. The autocorrelation decays were fit to a model for one-dimensional diffusion incorporating both fast and slow components having apparent diffusion coefficients, D f and D s , differing by a factor of ∼100. The fast mechanism was attributed to hindered bulk-like diffusion in the central pore cavity, while slow diffusion likely involved absorption of the dye to the pore surfaces. Unfortunately, important evidence of diffusion anomalies is lost in the broad autocorrelation decays obtained. The method of multifractal detrended fluctuation analysis (MF-DFA) was applied to the same data as a means to overcome this limitation. MF-DFA revealed that time series acquired from within the nanopores were multifractal and exhibited evidence of anomalous superdiffusion, likely resulting from the participation of a desorption-mediated diffusion mechanism. Monte Carlo simulations of time series modeling desorption-mediated diffusion in cylindrical nanopores provided support for this assignment. Here, the new knowledge gained affords an improved understanding of hydrocarbon dynamics within nanoporous oil and gas shales.

Diffusion↗
DOE OSTI2025

Distinguishing Desirable and Undesirable Reactions in Multicomponent Systems for Redox Activation of the Uranyl Ion

Although it has been established that covalent functionalization of the U–O bonds in the uranyl dication (UO 2 2+ ) generally requires use of strong reductants and electrophiles, little work has examined how interactions between the individual reaction components could affect final outcomes in solution. Here, the patterns of such reactivity have been studied in a UO 2 2+ -containing model system supported by a workhorse pentadentate ligand, 2,2′-[(methylimino)bis(2,1-ethanediylnitrilomethylidyne)]bis-phenol. Oxo activation and functionalization have been tested with (i) electrochemical and chemical reduction, and (ii) coordinating and noncoordinating solvents. In acetonitrile, uranyl reduction was achieved cleanly, but treatment of the reduced species with tris(pentafluorophenyl)borane (BCF) resulted in a mixture of products arising from direct electron transfer to BCF. In dichloromethane (CH 2 Cl 2 ), electrochemical reduction of uranyl was achieved cleanly, but clean chemical reactivity was inaccessible. Despite these challenges, one trinuclear and oxo-deficient uranium-containing product was crystallized from CH 2 Cl 2 solution and characterized; thus, desirable electrophilic reactivity can proceed to some degree in CH 2 Cl 2 with BCF. Computational studies were used to investigate the properties of the trinuclear uranium product and the changes that could be inducible by further reduction. Here, taken together, the reactivity patterns identified here could inform design of improved systems for actinyl oxo functionalization.

cations↗
DOE OSTI2025

Accurate Force Field for Carbon Dioxide–Silica Interactions Based on Density Functional Theory

Fluid–silica interfaces are ubiquitous in chemistry, occurring in both natural geochemical environments and practical applications ranging from separations to catalysis. Simulations of these interfaces have been, and continue to be, a significant avenue for understanding their behavior. A constraining factor, however, is the availability of accurate force fields. Most simulations use traditional “mixing rules” to determine nonbonded dispersion interactions, an approach that has not been critically examined. Here, in this study, we present Lennard-Jones parameters for the interaction of carbon dioxide with silica interfaces that are optimized to reproduce density functional theory (DFT)-based binding energies. The modeling is based on the recently developed silica-DDEC force field, whose atomic charges are consistent with DFT calculations. Standard mixing rules are found to predict weaker CO 2 binding to silica than that obtained from DFT, an effect corrected by the optimized parameters given here. This behavior extends to other silica force fields (Clayff and Gulmen-Thompson), and the present Lennard-Jones parameters improve their performance as well. The effects of improved Lennard-Jones parameters on the structural and dynamical properties of condensed CO 2 in silica slit pores are also examined.

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

Ab Initio Molecular Dynamics Study of the Reduction of Acetone by the Hydrated Electron

We have investigated the reaction dynamics of the reduction of acetone by the solvated electron in water using ab initio molecular dynamics simulations at 298 and 373 K. The rate constants derived from the simulations are consistent with experimental observations that the reaction has a low activation energy. Detailed analyses carried out to shed light on the mechanism of the reaction show that solvent reorganization plays a key role as a reaction coordinate, as expected from the Marcus electron transfer theory. Furthermore, constrained density functional theory calculations indicate that the electronic coupling is large, placing the reaction in the adiabatic limit. Indeed, the activation energies and rate constants of the simulations are in accordance with the predictions of adiabatic Marcus theory.

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

Imaging Three-Dimensional Molecular Structure and Dynamics with Multiparticle Covariance and Cumulant Coulomb Explosion Analysis

Coulomb explosion imaging (CEI) provides a direct means of imaging molecular geometry by correlating fragment ion momenta following the fragmentation of a molecular polycation. Here, we demonstrate the use of three-body covariance and four-body cumulant analysis to extract three-dimensional (3D) structural information from the X-ray-induced Coulomb explosion of tert-butyl iodide (C4H9I). Site-selective ionization at the iodine 4d edge with intense femtosecond soft X-ray pulses from an X-ray free-electron laser (XFEL) enables rapid charge buildup and molecular breakup. By correlating ionic fragments in the molecular frame, we isolate complete dissociation channels and reveal subtle structural changes, such as umbrella-type motion of the branched alkyl chain, during the ionization process. Comparison with point-charge simulations of the Coulomb explosion shows close agreement, validating the approach. Furthermore, these results establish covariance/cumulant mapping as a powerful strategy for imaging complex three-dimensional molecular structures and point the way toward time-resolved CEI using both XFEL and tabletop sources for capturing ultrafast structural dynamics.

Imaging↗
DOE OSTI2025

Ultrafast and Highly Mobile Photocarriers in Monolayer WSe 2 Doped by α-RuCl 3

We report an experimental investigation of photocarrier dynamics in a heterostructure formed by a WSe 2 monolayer stacked on a multilayer α-RuCl 3 flake. The samples were fabricated using mechanical exfoliation and dry transfer techniques. Photoluminescence measurements showed that the excitonic photoluminescence of WSe 2 is quenched by more than 3 orders of magnitude upon contact with α-RuCl 3 . Transient absorption measurements revealed an ultrashort photocarrier lifetime of 0.5 ps in the heterostructure. Spatially resolved transient absorption microscopy measurements were employed to probe the transport properties of these photocarriers, resulting in a room-temperature diffusion coefficient of 370 cm 2 s –1 , which is higher than the exciton diffusion coefficients of most previously studied 2D semiconductors. Furthermore, these results suggest that α-RuCl 3 -doped WSe 2 could be utilized in high-speed optoelectronic devices.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Attosecond control and measurement of chiral photoionization dynamics

Many chirality-sensitive light–matter interactions are governed by chiral electron dynamics. Therefore, the development of advanced technologies making use of chiral phenomena would critically benefit from measuring and controlling chiral electron dynamics on their natural attosecond timescales. Such endeavours have so far been hampered by the lack of characterized circularly polarized attosecond pulses, an obstacle that has recently been overcome. Here we introduce chiroptical spectroscopy with attosecond pulses and demonstrate attosecond coherent control over photoelectron circular dichroism (PECD), as well as the measurement of chiral asymmetries in the forward–backward and angle-resolved photoionization delays of chiral molecules. We show that co-rotating attosecond and near-infrared (IR) pulses can nearly double the PECD and even change its sign compared with single-photon ionization. We demonstrate that chiral photoionization delays depend on both polar and azimuthal angles of photoemission in the light-propagation frame, requiring 3D momentum resolution. We measure forward–backward chiral-sensitive delays of up to 60 as and polar-angle-resolved photoionization delays of up to 240 as, which include an asymmetry of about 60 as originating from chirality in the continuum–continuum transitions. Attosecond chiroptical spectroscopy opens the door to quantitatively understanding and controlling the dynamics of chiral molecules on the electronic timescale.

Atomic and molecular interactions with photons↗
DOE OSTI2025

Tracking oxygen vacancy migration in memristor devices using operando hard X-ray photoelectron spectroscopy

Fabrication of ultrathin (sub-2 nm) oxide semiconductor memristors poses the fundamental challenge of achieving oxide growth with atomic precision in terms of electronic structure and defect formation. Recently, ultrathin memristors consisting of bilayers of mixed Al 2 O 3 and MgO atomic layers were fabricated using an in vacuo atomic layer deposition process. This approach offers a unique platform for precise atomic control of oxygen vacancies in the device in which the vacancies are introduced by placing MgO atomic layers between pristine Al 2 O 3 layers. In this work, we present a systematic operando Hard X-ray PhotoElectron Spectroscopy (HAXPES) study of the switching of such memristors, combined with complementary current-voltage and capacitance-voltage (C-V) measurements. We used a memristor stack of Al 2 O 3 /MgO atomic layers, with the MgO-containing oxide deposited on the bottom Al metal electrode and a pure Al 2 O 3 layer below the top Pd electrode. HAXPES analysis shows a substantial change in the chemical shift of the Aluminum oxide when switching between the ”OFF” and ”ON” states indicative of a redistribution of oxygen vacancies in the device active layer. Interestingly, subsequent switching to the OFF state shows hysteretic behavior indicating the retention of some oxygen vacancies in the top Al 2 O 3 layer. This vacancy retention can be correlated with the stochastic behavior of the switching voltage observed in these devices. C–V measurements show a clear frequency-dependent response in the OFF state, consistent with enhanced polarization and vacancy trapping at low frequencies.

47 OTHER INSTRUMENTATION↗
DOE OSTI2026

Effects of confinement and pressure on the structure and dynamics of carbon dioxide in silica slit pores

An understanding of carbon dioxide fluid properties within geological mesopores is important in applications ranging from carbon sequestration to shale oil recovery. Here, a molecular dynamics study is presented that aims to shed light on these systems by simulation of CO 2 fluid confined in β-cristobalite silica slit pores with different pore widths and pressure conditions. The weakly associating nature of carbon dioxide leads to little difference in structural and interfacial dynamical properties for different pore sizes and pressures. Rather, the behavior is found to be dominated by the entropic effects, namely, how the CO 2 organizes next to the silica surface. The CO 2 self-diffusion coefficient shows the strongest pore size dependence. It is strongly diminished in small pores and does not reach the bulk fluid value even in 6 nm pores. It also decreases with pressure, yielding an activation volume that increases with pore size. These results provide new insight into the behavior of a compressible fluid in nanoscale confinement.

Godahewa, Sahan M. [Univ. of Kansas, Lawrence, KS ↗