Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “initial state”

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 1,153 records · Page 64
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Improved Benchmarking of Cohesive Elements in Abaqus Standard for Predicting Disbond and Delamination in Composite Structures

Traditional approaches for aircraft certification require the assumption of an initial flaw condition, either represented as barely visible impact damage (BVID) or through inclusion of a Teflon insert to serve as surrogate damage. Based on the initial composite damage state, the structure must be shown to demonstrate structural durability and damage tolerance (DaDT) according to the following criteria: a. Damage displays no detrimental growth under cyclic loading b. The structure is able to sustain design limit load (DLL) Currently, the only available manner for validating structural performance is through test. Since damage can occur over a wide variety of areas within a structure, this approach has proven to be increasingly expensive and time consuming for composite airframes and acreage structure within the design-test-certification building block. A further complicating factor is the requirement to accurately capture the most critical damage morphologies as a starting condition. To understand the severity of the damage, it is either required to experimentally determine the most critical areas at tremendous expense or rely on legacy data of similar structural testing, which limits design space expansion. A preferred solution is to use advanced analysis to provide improved understanding of load margins for critical locations based on a wide variety of potential starting damage conditions. The standard industry approach for DaDT certification adheres to the use of the traditional virtual crack closure technique (VCCT) method. VCCT is generally a preferred method because it conforms to the current certification principles of damage from a known flaw, and when used correctly, can be effective at predicting delamination propagation under static and cyclic loading. The VCCT method requires the inclusion of an initial flaw in the finite element (FE) model requiring a-priori knowledge of the flaw location. This in turn requires a plethora of analysis cases to be examined to cover a reasonable span of potential damage states. Additionally, the VCCT approach requires node-to-node connectivity rendering it incompatible with the best practices and approaches for using continuum damage mechanics (CDM) based progressive damage and failure analysis (PDFA) tools within a typical FE solver. Alternatives to VCCT have emerged in the form of cohesive elements which utilize the cohesive zone model (CZM). Unlike VCCT which models linear elastic fracture mechanics, cohesive elements couples continuum and fracture based responses through the use of bilinear traction separation laws. These laws are defined based on a penalty stiffness, a cohesive strength, and a strain energy release rate. The approach can be mesh regularized with native cohesive elements within many FE solvers such as Abaqus and LS-DYNA. In Phase I of the NASA Advanced Composites Consortium (ACC) post-buckled stiffened panel with BVID, Strength and Life [1], the performance of cohesive elements were benchmarked in comparison to VCCT and LEFM solutions and showed good agreement using Abaqus explicit [2]. To realize savings on current and future programs, it is still necessary to close technical gaps related to the use of cohesive elements with Abaqus Standard. Within a program environment, standard finite element analysis is the preferred analytical capability for quasi-static loading as it eliminates uncertainty due to oscillatory behavior commonly seen with explicit analysis. This oscillatory behavior creates difficulties in writing margins of safety based on the analysis. The use of negative tangent stiffness material models complicates convergence which typically requires the use of numerical controls such as viscous damping to overcome. To date, there has not been a comprehensive study on how to establish best practices for cohesive element convergence for predictive capability within the Abaqus implicit solver. In pursuit of these goals, under the NASA ACC program, several numerical benchmark problems were proposed including pure mode I (double cantilevered beam – DCB), pure mode II (end notch flexure – ENF), and symmetric/unsymmetric evolving mixed mode (single leg bend – SLB). This paper focuses on the use of cohesive elements to model the delamination through the use of CZM. Specifically, finite element models for the DCB, ENF, symmetric SLB, and unsymmetric SLB, are developed and various solution controls for convergence are studied to develop a best practice. Once the best practice has been developed, the predictive capability of the objective CZM model is used to analyze the hat pull-off strength of a standard hat stiffened configuration under various loading conditions.

Abaqus↗
DOE OSTI2026

Metal-Centered Photoredox Catalysis using d 6 Transition Metal-based Chromophores

Controlling excited-state reactivity in transition metal chromophores is central to light-driven chemistry, with applications spanning from solar energy conversion to photoredox catalysis. Ru(II) and Ir(III) chromophores currently dominate applications within this field due to their long-lived charge-transfer (CT) excited states, which support bimolecular electron transfer chemistry. However, their elemental scarcity and relative lack of selectivity for differentiating oxidative versus reductive pathways has in part motivated the exploration of alternative platforms that could, in principle, offer distinct mechanistic opportunities. Earth-abundant, valence-isoelectronic 3d 6 complexes of Fe(II) and Co(III) have attracted significant attention due to their electronic similarity to Ru(II) and Ir(III), suggesting the potential for complementary photophysical behavior. In contrast to their heavier congeners, however, the weaker ligand fields of first-row transition metals promote rapid relaxation from initially populated CT states to low-lying metal-centered (MC) excited states, typically of 5 T 2 or 3 T 1 character. Because these MC states involve redistribution of electron density within MC d-orbitals rather than charge separation between metal and ligand orbitals, they were historically considered poor candidates for use in photochemical transformations.

Charge transfer↗
DOE OSTI2020

Measurement of nuclear effects on ψ(2S) production in p-Pb collisions at $ \sqrt{s_{\mathrm{NN}}}$ = 8.16 TeV

Inclusive ψ (2S) production is measured in p-Pb collisions at the centre-of-mass energy per nucleon-nucleon pair \( \sqrt{s_{\mathrm{NN}}} \) = 8.16 TeV, using the ALICE detector at the CERN LHC. The production of ψ (2S) is studied at forward (203 < y cms < 3 . 53) and backward (-4.46 < y cms < -2.96) centre-of-mass rapidity and for transverse momentum p T < 12 GeV/ c via the decay to muon pairs. In this paper, we report the integrated as well as the y cms - and p T -differential inclusive production cross sections. Nuclear effects on ψ (2S) production are studied via the determination of the nuclear modification factor that shows a strong suppression at both forward and backward centre-of-mass rapidities. Comparisons with corresponding results for inclusive J/ ψ show a similar suppression for the two states at forward rapidity (p-going direction), but a stronger suppression for ψ (2S) at backward rapidity (Pb-going direction). As a function of p T , no clear dependence of the nuclear modification factor is found. The relative size of nuclear effects on ψ (2S) production compared to J/ ψ is also studied via the double ratio of production cross sections [ σ ψ (2S) /σ J /ψ ]pPb/[ σ ψ (2S) /σ J /ψ ] pp between p-Pb and pp collisions. The results are compared with theoretical models that include various effects related to the initial and final state of the collision system and also with previous measurements at \( \sqrt{s_{\mathrm{NN}}} \) = 5.02 TeV.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

An NV − center in magnesium oxide as a spin qubit for hybrid quantum technologies

Recent predictions suggest that oxides, such as MgO and CaO, could serve as hosts of spin defects with long coherence times and thus be promising materials for quantum applications. However, in most cases, specific defects have not yet been identified. Here, by using a high-throughput first-principles framework and advanced electronic structure methods, we identify a negatively charged complex between a nitrogen interstitial and a magnesium vacancy in MgO with favorable electronic and optical properties for hybrid quantum technologies. We show that this NV − center has stable triplet ground and excited states, with singlet shelving states enabling optical initialization and spin-dependent readout. We predict several properties, including absorption, emission, and zero-phonon line energies, as well as zero-field splitting tensor, and hyperfine interaction parameters, which can aid in the experimental identification of this defect. Our calculations show that due to a strong pseudo-Jahn Teller effect and low-frequency phonon modes, the NV − center in MgO is subject to a substantial vibronic coupling. We discuss design strategies to reduce such coupling and increase the Debye-Waller factor, including the effect of strain and the localization of the defect states. We propose that the favorable properties of the NV − defect, along with the technological maturity of MgO, could enable hybrid classical-quantum applications, such as spintronic quantum sensors and single qubit gates.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Impact of CuInP 2 S 6 –metal interfaces on the stabilization of polar phases and polarization switching

The multifunctionality of two-dimensional ferroelectric CuInP 2 S 6 (CIPS) arises from the existence of multiple polar phases combined with a high ionic conductivity that facilitates polarization switching in unusual ways. The van der Waals (vdW) layered structure provides ultrathin flakes and ideal interfaces to integrate with other materials for microelectronics and neuromorphic elements. However, device integration necessitates metal contacts to read, write, or transmit signals. In this work, we find that different types of metal–CIPS interfaces strongly impact the stabilization of specific polar phases and the field-induced transitions between the polarization states. Cu electrodes initially suppress the piezoresponse, whereas, at CIPS–Ag interfaces, the electromechanical signal is increased. Under electric fields, the Cu electrodes, Ag electrodes and surrounding CIPS surfaces can show distinct switching behavior as different phases and polarization orientations are stabilized. These findings highlight that metal–CIPS interfaces provide the opportunity to optimize functional material properties.

Ferroelectrics and multiferroics↗
NASA NTRS1987

A numerical study of the nonlinear thermal stability of solar loops

A time-dependent numerical model is used to investigate the nonlinear thermal stability of static loops of various heights. Simulations show that the instability of a hot state with loop heights of less than about 1000 km is physically significant, with an initially hot atmosphere in low-lying compact loops evolving to an extended atmosphere with temperatures far below 100,000 K. Results also show that high-lying loops are stable to all reasonable perturbations, including those of large initial amplitude and long wavelength. The simulation results suggest that low-lying compact loops should not be common to the sun, and that cool loops with temperatures near 100,000 K must be formed in the cool state initially and cannot evolve from preexisiting loops.

Klimchuk, J. A.↗
DOE OSTI2020

Direct observation of minibands in a twisted graphene/WS 2 bilayer

Stacking two-dimensional (2D) van der Waals materials with different interlayer atomic registry in a heterobilayer causes the formation of a long-range periodic superlattice that may bestow the heterostructure with properties such as new quantum fractal states or superconductivity. Recent optical measurements of transition metal dichalcogenide (TMD) heterobilayers have revealed the presence of hybridized interlayer electron-hole pair excitations at energies defined by the superlattice potential. The corresponding quasiparticle band structures, so-called minibands, have remained elusive, and no such features have been reported for heterobilayers composed of a TMD and another type of 2D material. We introduce a new x-ray capillary technology for performing microfocused angle-resolved photoemission spectroscopy with a spatial resolution of ~1 μm, and directly observe minibands at certain twist angles in mini Brillouin zones (mBZs). We discuss their origin in terms of initial and final state effects by analyzing their dispersion in distinct mBZs.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Low Cost (CAPEX and variable): Tool design for cell and module fabrication with thin, free-standing silicon wafers

This project aimed to develop technologies that can potentially enable free-standing thin (<80 μm) wafer in today’s manufacturing lines with high production yield, and thereby reduce capex barriers of silicon photovoltaics cells and modules. One of the major benefits is that thin wafer dramatically reduces the amount of polysilicon required. As a result, it can lead to reduction in the capital expenditures associated with polysilicon refining and wafer fabrication, which together are more than half of the total capital expenditure to manufacture Si PV module. We focused our efforts on developing the tools needed to enable high yield fabrication of wafer, cell, and module with thin silicon wafers. However, as the wafer thickness reduces, the major challenge is that wafer breakage increases significantly. Three technological areas were explored in this project to improve the production yield of the silicon wafer, namely detection of edge cracks via dark-field near-infrared (NIR) scattering; (2) wafer handling using controlled temperature profiles; (3) manufacturable low-stress cell interconnection for multiwire. First, the formation of wafer cracks in submillimeter length is one of the reasons that causes wafer breakages. Crack detection tools are needed to enable us to locate and track the wafer crack during manufacturing, so that we can improve the process to reduce initialization. The state-of-the-art crack detection technique cannot fulfill the need for measuring submillimeter edge cracks detrimental for thin wafers. The prototype developed in this project demonstrates the scanning of microcracks near wafer edges. With a semi-automatic laboratory setup, the submillimeter cracks were reliably detected near the edges in multi-Si wafers. The smallest detectable crack is 200 µm in length in slow scans; and submillimeter cracks are detected in high-throughput scans at the scan speed of >0.5 m/s, which is compatible with the inline detection of a manufacturing line at least 1 sec/wafer. This detection limit is a significant advancement in comparison to the benchmarked industrial tool. Second, wafer handling with the well-controlled temperature profile was thought to be the solution to reduce crack initiation and propagation during the manufacturing. However, without crack detection being widely adopted in production line, we did not find a strong industrial pull toward this technology. We did an initial literature survey and then diverted our efforts to the other tasks. Third, the innovation on low-stress multi-wire interconnection tackles a fundamental problem in the standard interconnection scheme. The standard over-under “zig-zag” interconnection induces a significant amount of stress into the soldering point on the cell whenever PV module is under stress, e.g., thermal cycling, transportation, and installation. Therefore, the interconnection process was re-designed in this project to allow for significant movement between adjacent solar cells, e.g., due to thermal expansion, without building up stresses on the solar cells or solder joints. The new interconnection method with the cross-connect wire also simplifies the tabbing and stringing process by replacing the standard over-under method with an off-cell interconnect from top to bottom. A manual tabbing and stringing tool for this new process was prototyped in the lab to demonstrate the fabrication of this new interconnection design. During the test with brass sheets as our “testing cells”, it was found that the mechanical cycling test only broke the interconnection after more than 50,000 cycles, which is equivalent to more than 130 years of the day-and-night thermal cycles. Lastly, throughout the project, we continuously analyzed the PV market with techno-economic analysis to identify the opportunity for thin Si adoption. Even though the drastic cost reduction has already happened in the past five years, our analysis results indicated that we can still save quite significantly in both manufacturing cost and factory capex, Particularly, in order to grow the PV manufacturing capacity to multi-terawatt level, reducing the thickness of silicon wafer is one of the most effective ways to quickly reduce factory capex, and sustain the high growth rate. The technologies developed in this project are readily available to provide some assistances in tackling the production yield problem.

14 SOLAR ENERGY↗
NASA NTRS2012

Theoretical Studies of N2-broadened Half-widths of H2O Lines Involving High j States

Based on the properties of the energy levels and wave functions of H2O states, one can categorize H2O lines into individually defined groups such that within the same group, the energy levels and the wave functions associated with two paired lines have an identity property while those associated with different pairs have a similarity property. Meanwhile, by thoroughly analyzing processes used to calculate N2-broadened half-widths, it was found that the 'Fourier series' of W(sup a)(sub L(sub 1))(sub K(sub 1))(sub K(sub 1)) (t; j(sub f) T(sub f) and W(sup a)(sub L(sub 1))(sub K(sub 1))(sub K(sub 1)) (t; j(sub i) T(sub i), and a factor P(sub 222) (j(sub f) T(sub f) j(sub i) T(sub i)) are the key items in the Robert-Bonamy formalism to distinguish contributions to ReS2(r(sub c)) among different transitions of j(sub f) T(sub f) - j(sub i). However, these items are completely determined by the energy levels and the wave functions associated with their initial and final states and they must bear the latter's features as well. Thus, it becomes obvious that for two paired lines in the same group, their calculated half-widths must be almost identical and the values associated with different pairs must vary smoothly as their ji values vary. Thus, the pair identity and the smooth variation rules are established within individual groups of lines. One can use these rules to screen half-width data listed in HITRAN and to improve the data accuracies.

water↗
DOE OSTI2024

Initial energy-momentum to final flow: A general framework for heavy-ion collisions

The evolution of a relativistic heavy-ion collision is typically understood as a process that transmutes the initial geometry of the system into the final momentum distribution of observed hadrons, which can be described via a cumulant expansion of the initial distribution of energy density and is represented at leading order as the well-known eccentricity scaling of anisotropic flow. We extend this framework to include the contribution from initial momentum-space properties, as encoded in other components of the energy-momentum tensor. Furthermore, we confirm the validity of the framework in state-of-the-art hydrodynamic simulations of large and small systems. With this framework, it is possible to separate the effects of early time dynamics from those of final-state evolution, even in the case when the distribution of energy does not fully determine subsequent evolution, as for example, in small systems. Specifically, we answer the question of when and how azimuthal correlations from the initial state survive to the final state. In very small systems such as 𝑝−𝑝, for example, initial momentum degrees of freedom dominate over energy. Thus, even if the system forms a quark-gluon plasma that is well described by hydrodynamics, the usual hydrodynamic picture of the transmutation of initial geometry to final momentum anisotropy is broken. Nevertheless, we show that the hydrodynamic response to the full energy-momentum tensor can be well understood in a similar manner as larger systems. Additionally, this framework elucidates the generic features of the system's evolution that are responsible for the impressive success of hydrodynamic simulations, but which may still hold even in cases where hydrodynamics is not applicable.

Relativistic heavy-ion collisions↗
DOE OSTI2023

High harmonic generation governed by edge states in triangular graphene quantum dots

Here, we study theoretically interaction of short optical pulses with triangular graphene quantum dots. If such quantum dots have zigzag edges then there are in-gap degenerate edge states, the number of which depends on the size of the system. We show that the nonlinear optical response, such as high harmonic generation, of triangular quantum dots is sensitive to the initial electron population of their edge states. In general, the emission spectra of quantum dots have weak dependence on the number of occupied edge states, but if half of the edge states are initially occupied, which can be realized only in the quantum dots with even number of edge states, then the even high harmonics are strongly suppressed. The suppression is the strongest when the frequency of the pulse is well below the band gap and it is weak when the pulse frequency becomes comparable to the band gap of the system.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

U.S. Distribution Transformer Demand Phase III - Key Drivers and Managing Demand [Slides]

This presentation demonstrates a significant analysis, on forecasting the demand for distribution transformers. The analysis is conducted for the United States, estimating the initial in-service capacity of these assets, and forecasting demand for these assets through 2050 with several sensitivities conducted. It examines not only demand for these assets, but importantly, how utility planning practices can impact the demand for theses assets in time. Under load growth scenarios, whether utilities practice like-for-like replacement strategies as failures occur, or whether they practice proactive up-sizing, anticipating electric load growth, can have major impacts on future demand. It also examines several other growth factors, such as the increasing demand for step-up transformers, which share many of the same characteristics as distribution transformers, and the demand for specific transformers for large project growth from data centers.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Evidence for Dissociation and Ionization in Shock Compressed Nitrogen to 800 GPa

Triple bonding in the nitrogen molecule (N 2 ) is among the strongest chemical bonds with a dissociation enthalpy of 9.8 eV/molecule. Nitrogen is therefore an excellent test bed for theoretical and numerical methods aimed at understanding how bonding evolves under the influence of the extreme pressures and temperatures of the warm dense matter regime. Here, we report laser-driven shock experiments on fluid molecular nitrogen up to 800 GPa and 4.0 g/cm 3 . Line-imaging velocimetry measurements and impedance matching method with a quartz reference yield shock equation of state data of initially precompressed nitrogen. Comparison with numerical simulations using path integral Monte Carlo and density functional theory molecular dynamics reveals clear signatures of chemical dissociation and the onset of L-shell ionization. Furthermore, combining data along multiple shock Hugoniot curves starting from densities between 0.76 and 1.29 g/cm 3 , our study documents how pressure and density affect these changes in chemical bonding and provides benchmarks for future theoretical developments in this regime, with applications for planetary interior modeling, high energy density science, and inertial confinement fusion research.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Pseudo-equilibrium theory for extrinsic doping control of the topological semimetal Cd 3 As 2

The standard approach for predicting defect equilibria from first principles assumes that the solid-state system is initially in a thermodynamic equilibrium with the external atomic reservoirs. This “growth step” is then often followed by a temperature quench in a “pseudo-equilibrium” in which some or all defect concentrations are frozen in until only the Fermi level E F remains to be equilibrated. However, this protocol does not account for the possibility of site exchanges which can create important defect redistributions as long as short-range defect migration is kinetically permissible. To model this redistribution, we developed an approach to solve for the non-equilibrium chemical potentials as a function of temperature while maintaining the overall defect stoichiometry. We then apply this approach to the Dirac semimetal Cd 3 As 2 to model extrinsic doping with group 1/11 and 14 elements. Undoped Cd 3 As 2 exhibits an undesirable mismatch between E F and the Dirac point. This unintentional electron doping originates from intrinsic defects and is difficult to overcome through adjustment of synthesis conditions alone. Employing our pseudo-equilibrium modeling, we identify extrinsic doping strategies for realizing doping-balanced Cd 3 As 2 at the relatively low temperatures accessible in thin-film growth of this material.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Iron(III)–Oxo Cluster Chemistry with Dimethylarsinate Ligands: Structures, Magnetic Properties, and Computational Studies

A program has been initiated to develop Fe III /oxo cluster chemistry with the ‘pseudo-carboxylate’ ligand dimethylarsinate (Me 2 AsO 2 - ) for comparison with the well investigated Fe III /oxo/carboxylate cluster area. The synthesis and characterization of three polynuclear Fe III complexes are reported, [Fe 12 O 4 (O 2 C t Bu) 8 (O 2 AsMe 2 ) 17 (H 2 O) 3 ]Cl 3 (1), Na 2 [Fe 12 Na 2 O 4 (O 2 AsMe 2 ) 20 (NO 3 ) 6 (Me 2 AsO 2 H) 2 (H 2 O) 4 ](NO 3 ) 6 (2) and [Fe 3 (O 2 AsMe 2 ) 6 (Me 2 AsO 2 H) 2 (hqn) 2 ](NO 3 ) (3), where hqn is 8-hydroxyquinoline. The Fe 12 core of 1 is a type never previously encountered in Fe III carboxylate chemistry, consisting of two Fe 6 units each of which comprises two {Fe 3 (μ 3 -O 2- )} units bridged by three Me 2 AsO 2 - groups and linked into an Fe 12 loop structure by two anti-anti η 1 :η 1 :μ Me 2 AsO 2 - groups, a bridging mode extremely rare with carboxylates. 2 also consists of two Fe 6 units, differing in their ligation from those in 1, and this time linked together into a linear structure by a central {Na 2 (NO 3 ) 2 } bridging unit. 3 is a linear Fe 3 complex with no monoatomic bridges between Fe III ions, a very rare situation in Fe III chemistry with any ligands, and unprecedented in Fe carboxylate chemistry. The distinct differences observed in arsinate vs carboxylate ligation modes are rationalized largely based on the greater basicity of the former vs the latter. Variable-temperature dc and ac magnetic susceptibility data reveal all Fe 2 pairwise interactions to be antiferromagnetic. For 1 and 2, the different J ij couplings were estimated by use of a magnetostructural correlation for high nuclearity Fe III -oxo clusters and by density functional theory calculations using broken symmetry methods, allowing identification of their relative spin vector alignments and thus rationalization of their S = 0 ground states. The J ij values were then used as input values to give excellent fits of the experimental χM T vs T data. For 3, the fits of the experimental χM T vs T data to the Van Vleck equation or with PHI gave a very weak J 12 = -0.8(1) cm -1 (H = –2JŠ i ·Š j convention) between adjacent Fe III ions, and an S = 5/2 ground state. Furthermore, these initial Fe III arsinate complexes also provide structural parameters that help validate literature assignments of arsinate binding modes to iron oxide/hydroxide minerals as part of environmental concerns of using arsenic-containing herbicides in agriculture.

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

An accurate and efficient Ehrenfest dynamics approach for calculating linear and nonlinear electronic spectra

Linear and nonlinear electronic spectra provide an important tool to probe the absorption and transfer of electronic energy. Here, we introduce a pure state Ehrenfest approach to obtain accurate linear and nonlinear spectra that is applicable to systems with large numbers of excited states and complex chemical environments. We achieve this by representing the initial conditions as sums of pure states and unfolding multi-time correlation functions into the Schrödinger picture. By doing this, we show that one can obtain significant improvements in accuracy over the previously used projected Ehrenfest approach and that these benefits are particularly pronounced in cases where the initial condition is a coherence between excited states. While such initial conditions do not arise when calculating linear electronic spectra, they play a vital role in capturing multidimensional spectroscopies. We demonstrate the performance of our method by showing that it is able to quantitatively capture the exact linear, 2D electronic spectroscopy, and pump–probe spectra for a Frenkel exciton model in slow bath regimes and is even able to reproduce the main spectral features in fast bath regimes.

Chemistry↗
DOE OSTI2024

A Pseudoreversible Normalizing Flow for Stochastic Dynamical Systems with Various Initial Distributions

Here, we present a pseudoreversible normalizing flow method for efficiently generating samples of the state of a stochastic differential equation (SDE) with various initial distributions. The primary objective is to construct an accurate and efficient sampler that can be used as a surrogate model for computationally expensive numerical integration of SDEs, such as those employed in particle simulation. After training, the normalizing flow model can directly generate samples of the SDE’s final state without simulating trajectories. The existing normalizing flow model for SDEs depends on the initial distribution, meaning the model needs to be retrained when the initial distribution changes. The main novelty of our normalizing flow model is that it can learn the conditional distribution of the state, i.e., the distribution of the final state conditional on any initial state, such that the model only needs to be trained once and the trained model can be used to handle various initial distributions. This feature can provide a significant computational saving in studies of how the final state varies with the initial distribution. Additionally, we propose to use a pseudoreversible network architecture to define the normalizing flow model, which has sufficient expressive power and training efficiency for a variety of SDEs in science and engineering, e.g., in particle physics. We provide a rigorous convergence analysis of the pseudoreversible normalizing flow model to the target probability density function in the Kullback–Leibler divergence metric. Numerical experiments are provided to demonstrate the effectiveness of the proposed normalizing flow model.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Scalar Co-SIMP dark matter: models and sensitivities

In this work, we present UV completions of the recently proposed number-changing Co-SIMP freeze-out mechanism. In contrast to the standard cannibalistic-type dark matter picture that occurs entirely in the dark sector, the 3 → 2 process setting the relic abundance in this case requires one Standard Model particle in the initial and final states. This prevents the dark sector from overheating and leads to rich experimental signatures. We generate the Co-SIMP interaction with a dark sector consisting of two scalars, with the mediator coupling to either nucleons or electrons. In either case, the dark matter candidate is naturally light : nucleophilic interactions favor the sub-GeV mass range and leptophilic interactions favor the sub-MeV mass range. Viable thermal models in these lighter mass regimes are particularly intriguing to study at this time, as new developments in low-threshold detector technologies will begin probing this region of parameter space. While particles in the sub-MeV regime can potentially impact light element formation and CMB decoupling, we show that a late-time phase transition opens up large fractions of parameter space. These thermal light dark matter models can instead be tested with dedicated experiments. We discuss the viable parameter space in each scenario in light of the current sensitivity of various experimental probes and projected future reach.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗