Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interactive model”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Space-time evolution of particle emission in p–Pb collisions at $\sqrt{{s}_{\text{NN}}} = 5.02$ TeV with 3D kaon femtoscopy

The measurement of three-dimensional femtoscopic correlations between identical charged kaons (K ± K ± ) produced in p–Pb collisions at center-of-mass energy per nucleon pair $\sqrt{{s}_{\text{NN}}}=5.02$ TeV with ALICE at the LHC is presented for the first time. This measurement, supplementary to those in pp and Pb–Pb collisions, allows understanding the particle-production mechanisms at different charged-particle multiplicities and provides information on the dynamics of the source of particles created in p–Pb collisions, for which a general consensus does not yet exist. It is shown that the measured source sizes increase with charged-particle multiplicity and decrease with increasing pair transverse momentum. These trends for K ± K ± are similar to the ones observed earlier in identical charged-pion and ${\text{K}}_{\text{s}}^{0}{\text{K}}_{\text{s}}^{0}$ correlations in Pb–Pb collisions at various energies and in π ± π ± correlations in p–Pb collisions at $\sqrt{{s}_{\text{NN}}}=5.02$ TeV. At comparable multiplicity, the source sizes measured in p–Pb collisions agree within uncertainties with those observed in pp collisions, and there is an indication that they are smaller than those observed in Pb–Pb collisions. The obtained results are also compared with predictions from the hadronic interaction model EPOS 3, which tends to underestimate the source size for the most central collisions and agrees with the data for semicentral and peripheral events. Furthermore, the time of maximal emission for kaons is extracted. It turns out to be comparable with the value obtained in highly peripheral Pb–Pb collisions at the same energy, indicating that the kaon emission evolution is similar to that in p–Pb collisions.

Heavy Ion Experiments↗
DOE OSTI2025

Exact Ground State of Interacting Electrons in Magic Angle Graphene

One of the most remarkable theoretical findings in magic angle twisted bilayer graphene (TBG) is the emergence of ferromagnetic Slater determinants as exact ground states for the interacting Hamiltonian at the chiral limit. This discovery provides an explanation for the correlated insulating phase which has been experimentally observed at half filling. This work is the first mathematical study of interacting models in magic angle graphene systems. These include not only TBG but also TBG-like systems featuring four flat bands per valley, and twisted trilayer graphene systems with equal twist angles. We identify symmetries of the chiral limit of the Bistritzer-MacDonald Hamiltonian that are responsible for characterizing the Hartree-Fock ground states as zero energy many-body ground states. Furthermore, for a general class of Hamiltonian, we establish criteria that the ferromagnetic Slater determinants are the unique ground states within the class of uniformly half-filled, translation invariant Slater determinants. We then demonstrate that these criteria can be explicitly verified for TBG and TBG-like systems at the chiral limit, using properties of Jacobi-θ$${\theta }$$ and Weierstrass-℘$${\wp }$$ functions.

Becker, Simon↗
DOE OSTI2026

Dihadron azimuthal asymmetry and light-quark dipole moments at the Electron-Ion Collider

Here, we propose a novel method to probe light-quark dipole moments by examining the azimuthal asymmetries between a collinear pair of hadrons in semi-inclusive deep inelastic lepton scattering off an unpolarized proton target at the Electron-Ion Collider. These asymmetries provide a means to observe transversely polarized quarks, which arise exclusively from the interference between the dipole and the Standard Model interactions, thereby depending linearly on the dipole couplings. We demonstrate that this novel approach can enhance current constraints on light-quark dipole operators by an order of magnitude, free from contamination of other new physics effects. Furthermore, it allows for a simultaneous determination of both the real and imaginary parts of the dipole couplings, offering a new avenue for investigating potential CP-violating effects at high energies.

dihadron fragmentation↗
DOE OSTI2026

Design trade-offs for residential retail tariffs and virtual power plants

Retail rate design and virtual power plants (VPPs) have the potential to shift customer electricity demand and provide economic benefits to utility customers. As the adoption of distributed energy resources (DERs) and flexible loads increases, retail tariff and program design can impact Bonbright's rate design principles including affordability, fairness, and economic efficiency. We model the effects of residential retail rates and VPP programs on power system costs in Massachusetts under a potential future system with high renewable energy and DER adoption. We model interactions among retail rate design, demand flexibility, and utility costs and identify trade-offs across different rate designs and VPP programs. We estimate that time-of-use (TOU) rates and VPP programs designed to avoid critical peak rates can lower overall system costs by 3.5 %-4.8 %. These lower costs translate to lower electricity bills for 62 %-91 % of customers, depending on the scenario. Although TOU rates with a critical peak VPP program can benefit all customer segments and are economically efficient, a VPP program with flat rates leads to the lowest overall bills for customers. We find that customers with loads that align with peak demand and who participate in critical peak VPP programs can underpay for their contribution to utility costs and shift costs to other customers. While our assumptions about mandatory TOU and/or critical peak pricing likely impact the magnitude of the results, the results highlight the trade-offs of these tariffs and programs and the importance of tariff and program design as demand becomes more flexible and responsive.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

A High-Performance Discrete-Element Framework for Simulating Flow and Jamming of Moisture Bearing Biomass Feedstocks

We developed and verified a high-performance open-source discrete element method (DEM) solver with simultaneously-supported feedstock-specific interaction models, including bonded-sphere, liquid bridge, cohesion, and non-linear contact models. Our solver uses parallel data structures on hybrid central and graphics processing unit (CPU/GPU) architectures, with favorable strong scaling performance observed for large problem sizes comprised of (100 M particles), and 4X single-node GPU speedup. The particles for corn stover feedstock were conceptualized and calibrated based on experimental measurements and results. Sensitivity analyses demonstrate that the mass flow rate from a wedge hopper is governed primarily by moisture content, friction coefficient, and cohesion energy density. The model is used to reproduce experimentally observed hopper jamming results, highlighting that the experimental no-flow trends can only be achieved by using non-spherical particles, liquid bridge and cohesion models, highlighting the importance of using concurrent feedstock specialized models for the effective representation of biomass material handling problems.

bioenergy↗
DOE OSTI2023

Measurements of the lightest hypernucleus ($_{Λ}^{3}{H}$): progress and perspective

The hyperon-nucleon (Y-N) interaction is important for the description of the equation-of-state of high baryon density matter. Hypernuclei, the cluster object of nucleons and hyperons, serve as cornerstones of a full understanding of the Y-N interaction. Recent measurements of the lightest known hypernucleus, the hypertriton’s ($_{Λ}^{3}{H}$) and anti-hypertriton’s ($\frac{H}{Λ}$$\bar{H}$) lifetime, mass and Λ separation energy have attracted interests on the subject. Its cross section and collective flow parameters have also been measured in heavy-ion collisions, which have revealed new features on its production mechanism. In this article we summarise recent measurements of $_{Λ}^{3}{H}$, focusing on the heavy-ion collisions. Here, we will discuss their implications for the $_{Λ}^{3}{H}$ properties and the constrains on the Y-N interaction models.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Interaction of Soil pH and Mineralogy Controls Soil Organic Matter Persistence through Changes in the Composition and Amount of Microbial Necromass

Microbial necromass–mineral associations are key to long-term soil organic matter (SOM) persistence. However, how soil pH and mineralogy interact to regulate SOM stability remains poorly understood. Here, we used artificial soils to test how three clay minerals (bentonite, kaolinite, and goethite), adjusted to four pH levels (5–8), affect microbial activity (respiration), microbial physiology (carbon use efficiency, CUE), microbial-derived residue material (necromass), and the formation and stability of mineral-associated organic matter (MAOM). Artificial soils were inoculated with a rhizosphere-derived microbial community cultured under the same pH conditions and on two representative simulated exudate types (organic acids and carbohydrates) and incubated for 6 weeks. In two complementary experiments, we added necromass from known microbial taxa to the same minerals across pH levels to isolate the role of necromass chemistry and loading. We found that soil pH shaped MAOM chemistry by altering microbial activity and necromass composition. In interaction with mineral type, pH also controlled MAOM thermal stability. Higher necromass loading weakened mineral-organic bonding, reducing MAOM stability, consistent with zonal mineral–organic interaction models. Our results demonstrate that microbial activity, rather than carbon use efficiency, better predicts MAOM formation and that pH-dependent necromass composition and loading govern MAOM persistence. These findings advance mechanistic understanding of SOM stabilization and have implications for predicting soil carbon dynamics under shifting environmental conditions.

carbon use efficiency↗
DOE OSTI2021

Davis Computational Spectroscopy Workflow—From Structure to Spectra

Here, we describe an automated workflow that connects a series of atomic simulation tools to investigate the relationship between atomic structure, lattice dynamics, materials properties, and inelastic neutron scattering (INS) spectra. Starting from the atomic simulation environment (ASE) as an interface, we demonstrate the use of a selection of calculators, including density functional theory (DFT) and density functional tight binding (DFTB), to optimize the structures and calculate interatomic force constants. We present the use of our workflow to compute the phonon frequencies and eigenvectors, which are required to accurately simulate the INS spectra in crystalline solids like diamond and graphite as well as molecular solids like rubrene. We have also implemented a machine-learning force field based on Chebyshev polynomials called the Chebyshev interaction model for efficient simulation (ChIMES) to improve the accuracy of the DFTB simulations. We then explore the transferability of our DFTB/ChIMES models by comparing simulations derived from different training sets. We show that DFTB/ChIMES demonstrates ~100× reduction in computational expense while retaining most of the accuracy of DFT as well as yielding high accuracy for different materials outside of our training sets. The DFTB/ChIMES method within the workflow expands the possibilities to use simulations to accurately predict materials properties of increasingly complex structures that would be unfeasible with ab initio methods.

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

Global Simulation of Electron Cyclotron Harmonic Wave Instability in a Storm-Time Magnetosphere

Electron cyclotron harmonic (ECH) waves are electrostatic emissions between the ECHs and play a dominant role for precipitating energetic electrons in the magnetotail. Statistically, the ECH wave intensity is stronger at nightside and dawnside than at dayside and duskside. In this study, we, for the first time, simulate the global ECH wave evolution during a geomagnetic storm event using Ring current Atmosphere interactions Model with Self-Consistent Magnetic field (RAM-SCB) combined with a linear growth rate solver.We find that the simulation results are generally consistent with the statistical and real-time observations. The ECH wave instability is much stronger at nightside and dawnside, compared to the instability at dayside and duskside. Before a geomagnetic storm (quiet time), the unstable regions of the ECH waves lie beyond $5R_E$ with a weak instability level. During the main phase of a geomagnetic storm, the unstable regions can extend to a lower altitude $(3R_E)$ with a strong instability level. Finally, during the recovery phase, the unstable regions return to $5$ ~ $6R_E$. We also find that the inner boundary of unstable ECH wave regions is coincident with the plasmapause location during the whole geomagnetic storm event.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

RAM-SCB Runs on Request at CCMC

The Ring current–Atmosphere interactions Model (RAM) with Self-Consistent magnetic (B) field (SCB) combines a large-scale kinetic model of ring current plasma with a three-dimensional (3-D) force-balanced model of the terrestrial magnetic field. RAM-SCB simulates the evolution of major ion species (H + , O + , and He + by default) and electrons as a function of azimuth, radial distance, energy, and pitch angle. Simulation outputs include the Dst index, and pressure and differential flux of the modeled species, thus providing benefit as a science code and to inform surface charging hazard within the model domain. Version 2.2 of this open-source simulation code has now—as of 8 August 2023—been made available for Runs-On-Request via the Community Coordinated Modeling Center.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Strongly correlated Chern insulators in magic-angle twisted bilayer graphene

Interactions between electrons and the topology of their energy bands can create unusual quantum phases of matter. Most topological electronic phases appear in systems with weak electron–electron interactions. The instances in which topological phases emerge only as a result of strong interactions are rare and mostly limited to those realized in intense magnetic fields. The discovery of flat electronic bands with topological character in magic-angle twisted bilayer graphene (MATBG) has created a unique opportunity to search for strongly correlated topological phases. Here we introduce a local spectroscopic technique using a scanning tunnelling microscope to detect a sequence of topological insulators in MATBG with Chern numbers C = ±1, ±2 and ±3, which form near filling factors of ±3, ±2 and ±1 electrons per moiré unit cell, respectively, and are stabilized by modest magnetic fields. One of the phases detected here (C = +1) was previously observed when the sublattice symmetry of MATBG was intentionally broken by a hexagonal boron nitride substrate, with interactions having a secondary role. We also demonstrate that strong electron–electron interactions alone can produce not only the previously observed phase, but also other unexpected Chern insulating phases in MATBG. The full sequence of phases that we observe can be understood by postulating that strong correlations favour breaking time-reversal symmetry to form Chern insulators that are stabilized by weak magnetic fields. Our findings illustrate that many-body correlations can create topological phases in moiré systems beyond those anticipated from weakly interacting models.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Measurement of the axial vector form factor from antineutrino–proton scattering

Scattering of high energy particles from nucleons probes their structure, as was done in the experiments that established the non-zero size of the proton using electron beams. The use of charged leptons as scattering probes enables measuring the distribution of electric charges, which is encoded in the vector form factors of the nucleon. Scattering weakly interacting neutrinos gives the opportunity to measure both vector and axial vector form factors of the nucleon, providing an additional, complementary probe of their structure. The nucleon transition axial form factor, F A , can be measured from neutrino scattering from free nucleons, ν μ n → μ – p and ν¯ μ p → μ + n , as a function of the negative four-momentum transfer squared (Q 2 ). Up to now, FA(Q 2 ) has been extracted from the bound nucleons in neutrino–deuterium scattering, which requires uncertain nuclear corrections. Here we report the first high-statistics measurement, to our knowledge, of the ν¯ μ p → μ + n cross-section from the hydrogen atom, using the plastic scintillator target of the MINERvA experiment, extracting F A from free proton targets and measuring the nucleon axial charge radius, r A , to be 0.73 ± 0.17 fm. The antineutrino–hydrogen scattering presented here can access the axial form factor without the need for nuclear theory corrections, and enables direct comparisons with the increasingly precise lattice quantum chromodynamics computations. Finally, the tools developed for this analysis and the result presented are substantial advancements in our capabilities to understand the nucleon structure in the weak sector, and also help the current and future neutrino oscillation experiments to better constrain neutrino interaction models.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Many-body reactive force field development for carbon condensation in C/O systems under extreme conditions

In this paper, we describe the development of a reactive force field for C/O systems under extreme temperatures and pressures, based on the many-body Chebyshev Interaction Model for Efficient Simulation (ChIMES). The resulting model, which targets carbon condensation under thermodynamic conditions of 6500 K and 2.5 g cm –3 , affords a balance between model accuracy, complexity, and training set generation expense. We show that the model recovers much of the accuracy of density functional theory for the prediction of structure, dynamics, and chemistry when applied to dissociative condensed phase systems at 1:1 and 1:2 C:O ratios, as well as molten carbon. Our C/O modeling approach exhibits a 104 increase in efficiency for the same system size (i.e., 128 atoms) and a linear system size scalability over standard quantum molecular dynamics methods, allowing the simulation of significantly larger systems than previously possible. We find that the model captures the condensed-phase reaction-coupled formation of carbon clusters implied by recent experiments, and that this process is susceptible to strong finite size effects. Overall, we find the present ChIMES model to be well suited for studying chemical processes and cluster formation at pressures and temperatures typical of shock waves. We expect that the present C/O modeling paradigm can serve as a template for the development of a broader high pressure–high temperature force-field for condensed phase chemistry in organic materials.

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

Calculation of the detonation state of HN 3 with quantum accuracy

HN 3 is a unique liquid energetic material that exhibits ultrafast detonation chemistry and a transition to metallic states during detonation. We combine the Chebyshev interaction model for efficient simulation (ChIMES) many-body reactive force field and the extended-Lagrangian multiscale shock technique molecular dynamics method to calculate the detonation properties of HN 3 with the accuracy of Kohn–Sham density-functional theory. ChIMES is based on a Chebyshev polynomial expansion and can accurately reproduce density-functional theory molecular dynamics (DFT-MD) simulations for a wide range of unreactive and decomposition conditions of liquid HN 3 . We show that addition of random displacement configurations and the energies of gas-phase equilibrium products in the training set allows ChIMES to efficiently explore the complex potential energy surface. Schemes for selecting force field parameters and the inclusion of stress tensor and energy data in the training set are examined. Structural and dynamical properties and chemistry predictions for the resulting models are benchmarked against DFT-MD. We demonstrate that the inclusion of explicit four-body energy terms is necessary to capture the potential energy surface across a wide range of conditions. Our results generally retain the accuracy of DFT-MD while yielding a high degree of computational efficiency, allowing simulations to approach orders of magnitude larger time and spatial scales. The techniques and recipes for MD model creation we present allow for direct simulation of nanosecond shock compression experiments and calculation of the detonation properties of materials with the accuracy of Kohn–Sham density-functional theory.

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

Ultra-wide bandgap β-Ga 2 O 3 films: Optical, phonon, and temperature response properties

Optical and phonon interactions of Ga 2 O 3 thin films with nanocrystalline morphology were studied at extreme temperatures. The films were grown using a sputtering technique and analyzed via temperature response transmission, Raman scattering, and high-resolution deep-UV photoluminescence (PL). Raman modes indicated that the structure corresponds to the β-phase. The optical-gap at the range of 77–620 K exhibited a redshift of ~200 meV, with a temperature coefficient of ~0.4 meV/K. The optical-gap at room-temperature is 4.85 eV. The electron–phonon interaction model at that temperature range pointed to a low energy phonon, ~31 meV, that is involved in the thermal properties of the optical-gap. Detailed Urbach energy analysis indicated that defects are the dominant mechanism controlling the band-edge characteristics even at an elevated temperature regime where phonon dominance is usually expected. Defects are attributed to the disordered forms of graphite that were detected via Raman scattering and to the granular morphology of the film. A deep-UV laser with an above bandgap exaction line of 5.1 eV was employed to map the PL of the films. The highly resolved spectra, even at room-temperature, show a strong emission of ~3.56 eV attributed to self-trapped holes (STHs). The STH is discussed and modeled in terms of the self-trapped exciton. Moreover, a very distinct but low-intensity emission was found at 4.85 eV that agrees with the value of the optical-gap and is attributed to bandgap recombination. The intensity ratio between the STH and that of the bandgap was found to be 6:1.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Resonance three-wave interactions and strange attractor

We report it is shown that the incorporation of linear sink/source terms in the three-wave resonance interaction model results in the time dependence of the wave amplitudes, which could exhibit the properties of a strange attractor. This finding demonstrates that the transition to turbulent dynamics of the waves could be related not only to the coupling of wave triads but also to the establishing of the strange attractor-like dynamics within individual wave triads.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Preliminary benchmarks and analysis of boundary conditions in a trenched microstructured silicon radiation detector

Microstructured neutron detectors have the benefit of enhanced neutron detection efficiency as compared to planar devices, achieved by etching 6 LiF-filled trenches on the top surface of a silicon PIN diode. This sensor geometry results in a complex electric field distribution and depletion characteristics within the diode under reverse bias. For the first time on record, the effects of a fixed oxide charge on the microstructured device depletion characteristics and mobile carrier transport is investigated. Prototype detectors were fabricated with non-conformal surface doping. Capacitance voltage and current voltage measurements were performed for these prototypes and compared with COMSOL Multiphysics simulations. A spectral response from an 241Am alpha particle source was acquired and analyzed. It was found that monoenergetic alpha particles produce three prominent peaks in the pulse height spectrum output by the device. The peaks were confirmed by simulations to correlate with dead layers and incident trajectories into the microstructure. It was also found that significant differences in pulse rise time result, corresponding with events arriving in a low-field region in the fins and a high-field region in the bulk. Geant4 was utilized for radiation transport, interaction modeling, and benchmarking the spectral data. The results of this simulation work provide confidence in the ability to attain and benchmark electrical characteristics and spectral data for semiconductor radiation detectors employing complex microstructures.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

The Gaussian nature of the band-edge of ZnO microcrystalline thin films

A straightforward analytical approach based on the derivative of the absorption coefficient is presented, which enables probing the nature of the band edge (BE) of ZnO microcrystalline films. The study was conducted via transmission experiments at temperatures of 77–532 K and repeated for samples annealed up to 1073 K. It was found that the derivative of the natural log of the BE absorption coefficient resulted in a Gaussian function. The Gaussian linewidth is used in the electron–phonon (e–p) interaction model to characterize the defect-state of the films. The BE of the as-grown film was found to exhibit no thermal dependence and no e–p coupling, indicative of a disordered crystal. Upon annealing and improvement of the film quality, the thermal phonons became more activated, but only above room temperature with a phonon energy of ~75 meV, while up to room temperature, the impact of phonons on the BE is insignificant. A disorder–order transition was determined to take place at an annealing temperature of ~673 K. X-ray diffraction concurs with these results. The study indicates that the prevalent defects are of structural nature due to the inherent granular morphology of the films. This defect was found to dominate the behavior of the BE even at the elevated temperature regime rather than thermal phonons.

36 MATERIALS SCIENCE↗