Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “branching process”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Gas Phase Synthesis of Benzene via the Propargyl Radical Self-Reaction

Polycyclic Aromatic Hydrocarbons (PAHs) have been invoked in fundamental molecular mass growth processes in our galaxy. We provide compelling evidence of the formation of the very first ringed aromatic and building block of PAHs - benzene - via the self-recombination of two resonantly stabilized propargyl (C 3 H 3 ) radicals in dilute environments using isomer selective synchrotron-based mass spectrometry coupled to theoretical calculations. Along with benzene, three other structural isomers (1,5-hexadiyne, fulvene, 2-ethynyl-1,3-butadiene) and o-benzyne are detected and their branching ratios quantified experimentally and verified with the aid of computational fluid dynamics and kinetic simulations. Furthermore, these results uncover molecular growth pathways not only in interstellar, circumstellar, and solar systems environments, but also in combustion systems gaining us a better understanding of the hydrocarbon chemistry of our universe.

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

Interactions between molecular-scale processes and hyporheic exchange for understanding Fe-S-C cycling in riparian wetlands (Final Report)

Wetlands represent some of the most productive ecosystems on the planet and critically influence global environmental health. Specifically, wetlands promote water quality by transforming nutrients and organic compounds and sequestering metals and contaminants. Riparian wetland hyporheic zones, where toxic surface water and anoxic groundwater mix, exhibit dynamic conditions that drive steep redox gradients and promote hotspots of diverse and fluctuating microbial activity. Changes in climate, water quality, and water quantity can disturb hydrologic flow and biogeochemical processing in these environments. Understanding how sulfate loading, impacted by hydrologic flux and anthropogenic inputs, influences iron and carbon cycling in wetlands will be crucial for predicting water quality issues driven by iron mineral precipitation and sorption, such as the release of heavy metals and other toxic elements. We used a fully integrated multi-scale and multi-method approach to develop a mechanistic understanding of how hydrologic flow influences coupled iron and sulfur cycles in riparian wetlands. This entailed hydrological, geochemical, and microbial observations at two locations: an anthropogenic sulfate-impacted riparian wetland in northern Minnesota, and a low-sulfate Fe-rich riparian wetland in Tims Branch at the Savannah River Site (SRS). These sites were characterized by hydrologically dynamic conditions where the stream and wetland systems oscillated between gaining (upward flow) and losing (downward flow) conditions that recharged the system occasionally with oxidants that fueled a variety of biogeochemical reactions. Aqueous geochemical measurements of surface water, groundwater, and porewater samples were made alongside solid-phase geochemical analyses of sediment gravity cores. Bulk X-ray absorption spectroscopy at the Advanced Photon Source (APS; Argonne) interrogated the speciation and distribution of Fe and S mineral phases of the sediments. Interestingly, it was discovered that, despite strongly reducing conditions, Fe(III) compounds and a variety of intermediate valence S compounds were stable in the subsurface. Indeed, compounds like thiosulfate, S(0), and intermediate valence organosulfur compounds were more prevalent than FeS and pyrite. The composition of the sediments did change with changes in hydrologic flow, showing their reactivity in changing redox conditions. The abundance of these intermediate S compounds was likely formed as a result of anaerobic oxidation by aqueous and solid-phase Fe(III) compounds, fueling a cryptic S cycle that is driving the breakdown of organic matter in the hyporheic zone. Microbiome surveys showed a core community that seemed stable across the landscape, but changed with increasing depth into the sediment. The community composition did not seem to change dramatically with changes in season or hydrologic flow, except for some organisms that have the potential to contribute to S cycling. More work is needed to confirm their functional activity. These fine process-scale analyses were placed within a dynamic field context using physical flow parameters from surface water and groundwater level measurements. These data helped shed light on sulfur-driven biogeochemical processes in hydrologically dynamic riparian wetlands, addressing a gap in our understanding about the impacts of pollution and other anthropogenic changes on ecologically sensitive environments.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Search for heavy resonances in final states with four leptons and missing transverse momentum or jets in pp collisions at $\sqrt{s}$ = 13 TeV with the ATLAS detector

A search for a new heavy boson produced via gluon-fusion in the four-lepton channel with missing transverse momentum or jets is performed. The search uses proton-proton collision data equivalent to an integrated luminosity of 139 fb -1 at a centre-of-mass energy of 13 TeV collected by the ATLAS detector between 2015 and 2018 at the Large Hadron Collider. This study explores the decays of heavy bosons: R → SH and A → ZH, where R is a CP-even boson, A is a CP-odd boson, H is a CP-even boson, and S is considered to decay into invisible particles that are candidates for dark matter. In these processes, S → invisible and H → ZZ. The Z boson associated with the heavy scalar boson H decays into all decay channels of the Z boson. The mass range under consideration is 390–1300 (320–1300) GeV for the R (A) boson and 220–1000 GeV for the H boson. No significant deviation from the Standard Model backgrounds is observed. The results are interpreted as upper limits at a 95% confidence level on the cross-section times the branching ratio of the heavy resonances.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Impacts of molecular architecture on the radiation-induced degradation and reaction kinetics of hydrophobic diglycolamides with the solvated electron and the dodecane radical cation

Given their proposed use as trivalent actinide–lanthanide separation ligands, the role of molecular architecture on the radiation robustness of diglycolamide (DGA) molecules has been investigated. This study examined three prototypical molecules with differences in their aliphatic chain architecture: N,N,N′,N′-tetra(n-octyl)diglycolamide (TODGA), N,N,N′,N′-tetra(2-ethylhexyl)diglycolamide (T2EHDGA), and N,N′-dimethyl-N,N′-dioctyldiglycolamide (DMDODGA). Rate coefficients and activation parameters are reported for the reactivity of each DGA with the solvated electron (e solv − ) and the corresponding dodecane radical cation (RH˙ + ) over the temperature range of 10.0 to 44.1 °C. These measurements indicate that DMDODGA is the most chemically reactive with both transient radicals, which may be attributed to this molecule's more accessible backbone. Complementary gamma dose accumulation studies (≤ 600 kGy) under envisioned process conditions—50 mM DGA in n-dodecane solvent—afforded dose constants for the loss of DGA of d = (3.41 ± 0.07) × 10 −3 , (4.19 ± 0.09) × 10 −3 , and (4.65 ± 0.09) × 10 −3 kGy −1 for T2EHDGA, DMDODGA, and TODGA, respectively. These dose constants indicate that varying DGA architecture affords subtle differences in chemical reactivity, leading to varying rates of radiolytic degradation under envisioned actinide–lanthanide separation conditions. However, more ambitious DGA frameworks, such as modifying the backbone, branching of the aliphatic chains, and/or changing the size of the chain may be required for larger gains in radiolytic longevity while optimizing actinide–lanthanide selectivity.

Arrhenius parameters↗
DOE OSTI2023

Efficient and selective dual-pathway polyolefin hydro-conversion over unexpectedly bifunctional M/TiO 2 -anatase catalysts

The huge scale of plastic waste generation and its environmental consequences drive the demands for catalytic plastic upcycling processes. Here, we report efficient polypropylene (PP) hydro-conversion over Ru and Ni supported on a sol-gel anatase TiO 2 (TiO 2 -A-SG). A small number of Brønsted acid sites on TiO 2 -A-SG enables hydrocracking, improving efficiency, steering the selectivity towards high-valued C 4-20 , and allowing more isomerization, compared to Ru-based monofunctional hydrogenolysis catalysts. The relative contribution of hydrocracking increases with lower hydrogen partial pressure, more branched substrates, and lower Ru loading. Ni/TiO 2 -A-SG exhibits superior hydrocracking activity and selectivity to Ni on conventional Brønsted-acidic supports, particularly zeolites. Further, mechanistic studies shows fast isomerization and sequential ß-scissions caused by strong polymer-catalyst interaction and fast cracking with abundant 3 C. This favors C 4-12 formation and prevents secondary reactions. This work demonstrated a novel, highly efficient noble-metal-free catalyst for plastic waste upcycling, while advancing the mechanistic understanding of polyolefin hydro-conversion.

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

Kinetics and Mechanism of the Singlet Oxygen Atom Reaction with Dimethyl Ether

Here, we combine in situ laser spectroscopy, quantum chemistry, and kinetic calculations to study the reaction of a singlet oxygen atom with dimethyl ether. Infrared laser absorption spectroscopy and Faraday rotation spectroscopy are used for the detection and quantification of the reaction products OH, H 2 O, HO 2 , and CH 2 O on submillisecond time scales. Fitting temporal profiles of products with simulations using an in-house reaction mechanism allows product branching to be quantified at 30, 60, and 150 Torr. The experimentally determined product branching agrees well with master equation calculations based on electronic structure data and transition state theory. The calculations demonstrate that the dimethyl peroxide (CH 3 OOCH 3 ) generated via O-insertion into the C–O bond undergoes subsequent dissociation to CH 3 O + CH 3 O through energetically favored reactions without an intrinsic barrier. This O-insertion mechanism can be important for understanding the fate of biofuels leaking into the atmosphere and for plasma-based biofuel processing technologies.

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

Promoting multiexciton interactions in singlet fission and triplet fusion upconversion dendrimers

Abstract Singlet fission and triplet-triplet annihilation upconversion are two multiexciton processes intimately related to the dynamic interaction between one high-lying energy singlet and two low-lying energy triplet excitons. Here, we introduce a series of dendritic macromolecules that serve as platform to study the effect of interchromophore interactions on the dynamics of multiexciton generation and decay as a function of dendrimer generation. The dendrimers (generations 1–4) consist of trimethylolpropane core and 2,2-bis(methylol)propionic acid (bis-MPA) dendrons that provide exponential growth of the branches, leading to a corona decorated with pentacenes for SF or anthracenes for TTA-UC. The findings reveal a trend where a few highly ordered sites emerge as the dendrimer generation grows, dominating the multiexciton dynamics, as deduced from optical spectra, and transient absorption spectroscopy. While the dendritic structures enhance TTA-UC at low annihilator concentrations in the largest dendrimers, the paired chromophore interactions induce a broadened and red-shifted excimer emission. In SF dendrimers of higher generations, the triplet dynamics become increasingly dominated by pairwise sites exhibiting strong coupling (Type II), which can be readily distinguished from sites with weaker coupling (Type I) by their spectral dynamics and decay kinetics.

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

ETG turbulence in a tokamak pedestal

This paper explores the fundamental characteristics of electron-temperature-gradient (ETG)-driven turbulence in the tokamak pedestal. The extreme gradients in the pedestal produce linear instabilities and nonlinear turbulence that are distinct from the corresponding ETG phenomenology in the core plasma. The linear system exhibits multiple (greater than ten) unstable eigenmodes at each perpendicular wave vector, representing different toroidal and slab branches of the ETG instability. Proper orthogonal decomposition of the nonlinear fluctuations reveals no clear one-to-one correspondence between the linear and nonlinear modes for most wave vectors. Moreover, nonlinear frequencies deviate strongly from those of the linear instabilities, with spectra peaking at positive frequencies, which is opposite the sign of the ETG instability. The picture that emerges is one in which the linear properties are preserved only in a narrow range of k-space. Outside this range, nonlinear processes produce strong deviations from both the linear frequencies and eigenmode structures. This is interpreted in the context of critical balance, which enforces alignment between the parallel scales and fluctuation frequencies. We also investigate the nonlinear saturation processes. We observe a direct energy cascade from the injection scale to smaller scales in both perpendicular directions. However, in the bi-normal direction, there is also nonlocal inverse energy transfer to larger scales. Neither streamers nor zonal flows dominate the saturation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Compare linear-system solver and preconditioner stacks with emphasis on GPU performance and propose phase-2 NGP solver development pathway

The goal of the ExaWind project is to enable predictive simulations of wind farms comprised of many megawatt-scale turbines situated in complex terrain. Predictive simulations will require computational fluid dynamics (CFD) simulations for which the mesh resolves the geometry of the turbines and captures the rotation and large deflections of blades. Whereas such simulations for a single turbine are arguably petascale class, multi-turbine wind farm simulations will require exascale-class resources. The primary physics codes in the ExaWind project are Nalu-Wind, which is an unstructured-grid solver for the acoustically incompressible Navier-Stokes equations, and OpenFAST, which is a whole-turbine simulation code. The Nalu-Wind model consists of the mass-continuity Poisson-type equation for pressure and a momentum equation for the velocity. For such modeling approaches, simulation times are dominated by linear-system setup and solution for the continuity and momentum systems. For the ExaWind challenge problem, the moving meshes greatly affect overall solver costs as reinitialization of matrices and recomputation of preconditioners is required at every time step. In this report we evaluated GPU-performance baselines for the linear solvers in the Trilinos and hypre solver stacks using two representative Nalu-Wind simulations: an atmospheric boundary layer precursor simulation on a structured mesh, and a fixed-wing simulation using unstructured overset meshes. Both strong-scaling and weak-scaling experiments were conducted on the OLCF supercomputer Summit and similar proxy clusters. We focused on the performance of multi-threaded Gauss-Seidel and two-stage Gauss-Seidel that are extensions of classical Gauss-Seidel; of one-reduce GMRES, a communication-reducing variant of the Krylov GMRES; and algebraic multigrid methods that incorporate the afore-mentioned methods. The team has established that AMG methods are capable of solving linear systems arising from the fixed-wing overset meshes on CPU, a critical intermediate result for ExaWind FY20 Q3 and Q4 milestones. For the fixed-wing strong-scaling study (model with 3M grid-points), the team identified that Nalu-Wind simulations with the new Trilinos and hypre solvers scale to modest GPU counts, maintaining above 70% efficiency up to 6 GPUs. However, there still remain significant bottlenecks to performance: matrix assembly (hypre), AMG setup (hypre and Trilinos) In the weak-scaling experiments (going from 0.4M to 211M gridpoints), it's shown that the solver apply phases are faster on GPUs, but that Nalu-Wind simulation times grow, primarily due to the multigrid-setup process. Finally, based on the report outcomes, we propose a linear solver path-forward for the remainder of the ExaWind project. Near term, the NREL team will continue their work on GPU-based linear-system assembly. They will also investigate how the use of alternatives to the NVIDIA UVM (unified virtual memory) paradigm affects performance. Longer term, the NREL team will evaluate algorithmic performance on other types of accelerators and merge their improvements back to the main hypre repository branch. Near term, the Trilinos team will address performance bottlenecks identified in this milestone, such as implementing a GPU-based segregated momentum solve and reusing matrix graphs across linear-system assembly phases. Longer term, the Trilinos team will do detailed analysis and optimization of multigrid setup.

17 WIND ENERGY↗
DOE OSTI2026

Process-oriented evaluation of quasi-stationary Rossby waves and their impact on surface air temperature extremes in dynamical downscaling over North America

Quasi-stationary Rossby waves are a crucial component of the general circulation and play a significant role in regional water and energy cycles, as well as in extreme events. However, process-oriented evaluation for Rossby waves is rarely performed for dynamical downscaling simulations. To close this gap, we evaluate three classes of dynamical downscaling approaches, with a focus on quasi-stationary Rossby waves and their impact on surface air temperature over North America during Northern Hemisphere summer. The three classes of models differ in the way large-scale forcing is provided: a limited-area model (LAM) constrained only by lateral boundary conditions, represented by RegCM4 from the North American branch of the Coordinated Regional Downscaling Experiment (NA-CORDEX), a LAM with spectral nudging to maintain consistency in large-scale dynamics with the forcing data, represented by the Weather Research and Forecasting (WRF) model simulation in NA-CORDEX, and a global variable-resolution model with smoothly varying grid spacings, represented by the Community Atmosphere Model version 5.4, with the Model for Prediction Across Scales (MPAS) as its dynamical core (CAM-MPAS). With no constraints on the atmospheric dynamics, CAM-MPAS exhibits several mean biases in the upper-level circulations over the Pacific Coast region: a weaker subtropical jet, a northward-shifted mid-latitude jet, and an overestimated southerly flow. With the lateral boundary constraint alone, RegCM4 also exhibits weaker jets and overestimated southerly winds off the West Coast. Rossby ray theory reveals that those wind biases direct incoming Rossby waves northward. The erroneously routed Rossby waves distort the relationship between the accumulation of wave activity over the US West Coast and surface temperature anomalies over the Southern Great Plains, which emerges approximately 4 d after the convergence of wave-activity flux in the ERA-Interim reanalysis. Furthermore, the response of heatwaves to the extreme wave activity flux is not reproduced by the two models, a serious drawback as a dynamical downscaling framework is expected to connect large-scale forcing to local-scale phenomena. The WRF model employing spectral nudging is largely free from the aforementioned problems. A pair of sensitivity simulations suggests that spectral nudging is the key to improving the dynamics of quasi-stationary Rossby waves and their impact on surface air temperature. Our results also demonstrate the effectiveness of Rossby wave diagnostics that allow for realistic background flows for assessing the credibility of dynamical downscaling over North America, where incoming Rossby waves propagate through complex circulation patterns before traveling across the continent.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Search for nonresonant Higgs boson pair production in final states with two bottom quarks and two photons in proton-proton collisions at $ \sqrt{s} $ = 13 TeV

A search for nonresonant production of Higgs boson pairs via gluon-gluon and vector boson fusion processes in final states with two bottom quarks and two photons is presented. The search uses data from proton-proton collisions at a center-of-mass energy of $ \sqrt{s} $ = 13 TeV recorded with the CMS detector at the LHC, corresponding to an integrated luminosity of 137 fb$^{−1}$. No significant deviation from the background-only hypothesis is observed. An upper limit at 95% confidence level is set on the product of the Higgs boson pair production cross section and branching fraction into $ \gamma \gamma \mathrm{b}\overline{\mathrm{b}} $. The observed (expected) upper limit is determined to be 0.67 (0.45) fb, which corresponds to 7.7 (5.2) times the standard model prediction. This search has the highest sensitivity to Higgs boson pair production to date. Assuming all other Higgs boson couplings are equal to their values in the standard model, the observed coupling modifiers of the trilinear Higgs boson self-coupling κ$_{λ}$ and the coupling between a pair of Higgs bosons and a pair of vector bosons c$_{2V}$ are constrained within the ranges −3.3 < κ$_{λ}$< 8.5 and −1.3 < c$_{2V}$< 3.5 at 95% confidence level. Constraints on κ$_{λ}$ are also set by combining this analysis with a search for single Higgs bosons decaying to two photons, produced in association with top quark-antiquark pairs, and by performing a simultaneous fit of κ$_{λ}$ and the top quark Yukawa coupling modifier κ$_{t}$.[graphic not available: see fulltext]

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Chemical insights into the multi-regime low-temperature oxidation of di-n-propyl ether: Jet-stirred reactor experiments and kinetic modeling

To further understand the combustion characteristics and the reaction pathways of acyclic ethers, in this work the oxidation of di-n-propyl ether (DPE) was investigated in a jet-stirred reactor (JSR) combined with a photoionization molecular-beam mass spectrometer. The experiments were carried out at near-atmospheric pressure (700 Torr) and over a temperature range of 425–850 K. Based on the experimental data and previous studies on ether oxidation, a new kinetic model was constructed and used to interpret the oxidation chemistry of DPE. In DPE oxidation, a high reactivity at low temperatures and two negative temperature coefficient (NTC) zones were observed. These behaviors are explained in this work by taking advantage of the obtained species information and the modeling analyses: the two NTC zones are caused by the competition of chain branching and termination reactions of the fuel itself and specific oxidation intermediates, respectively. Furthermore, the general requirements to have double-NTC behavior are discussed. A variety of crucial fuel-specific C 6 species, such as ketohydroperoxides and diones, were detected in the species pool of DPE oxidation. Their formation pathways are illuminated based on rate-of-production (ROP) analyses. Propanal was identified as the most abundant small molecule intermediate, and its related reactions have an important impact on the oxidation process of DPE. Both acetic acid and propionic acid were detected in high concentrations. A new formation pathway of propionic acid is proposed and incorporated into the kinetic model to achieve a more accurate prediction for propionic acid mole fractions.

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

S -factor and scattering-parameter extractions from ${}^{3}\mathrm{He}+{}^{4}\mathrm{He}{ \rightarrow }^{7}\mathrm{Be}+\gamma $

Previous studies of the reaction ${}^{3}\mathrm{He}+{}^{4}\mathrm{He}\to {}^{7}\mathrm{Be}+\gamma $ have mainly focused on providing the best central value and error bar for the S factor at solar energies. Experimental measurements of this capture reaction at higher energies, the ${}^{3}\mathrm{He}$–${}^{4}\mathrm{He}$ scattering phase shifts, as well as properties of ${}^{7}\mathrm{Be}$ and its excited state, have been used to constrain the theoretical models employed for this purpose. Here we show that much more information than was previously appreciated can be extracted from angle-integrated capture data alone. We use the next-to-leading-order (NLO) amplitude in an effective field theory (EFT) for ${}^{3}\mathrm{He}+{}^{4}\mathrm{He}\to {}^{7}\mathrm{Be}+\gamma $ to perform the extrapolation. At this order the EFT describes the capture process using an s-wave scattering length and effective range, the asymptotic properties of ${}^{7}\mathrm{Be}$ and its excited state, and short-distance contributions to the E1 capture amplitude. We extract the multi-dimensional posterior of all these parameters via a Bayesian analysis that uses capture data below 2 MeV. We find that properties of the ${}^{7}\mathrm{Be}$ ground and excited states are well constrained. The total S factor $S(0)\,={0.577}_{-0.016}^{+0.015}$ keV b, while the branching ratio for excited- to ground-state capture at zero energy, ${Br}(0)={0.406}_{-0.011}^{+0.013}$, both at 68% degree of belief. This S(0) is broadly consistent with other recent evaluations, and agrees with the previously recommended value $S(0)=0.56\pm 0.03\,\,\mathrm{keV}$ b, but has a smaller error bar. We also find significant constraints on ${}^{3}\mathrm{He}$–${}^{4}\mathrm{He}$ scattering parameters, and we obtain constraints on the angular distribution of capture gamma rays, which is important for interpreting experiments. The path forward for this reaction seems to lie with better measurements of the scattering phase shifts and S(E)'s angular dependence away from zero energy, together with better understanding of the asymptotic normalization coefficients of the ${}^{7}$Be bound states' wave functions. Lastly, data on these could further reduce the uncertainty on S(0).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Transitory sensitivity in automatic chemical kinetic mechanism analysis

Abstract Detailed chemical kinetic mechanisms are necessary for resolving many important chemical processes. As the chemistry of smaller molecules has become better grounded and quantum chemistry calculations have become cheaper, kineticists have become interested in constructing progressively larger kinetic mechanisms to model increasingly complex chemical processes. These large kinetic mechanisms prove incredibly difficult to refine and time‐consuming to interpret. Traditional sensitivity analysis on a large mechanism can range from inconvenient to practically impossible without special techniques to reduce the computational cost. We first present a new time‐local sensitivity analysis we term transitory sensitivity analysis. Transitory sensitivity analysis is demonstrated in an example to accurately identify traditionally sensitive reactions at an 18,000x speed up over traditional sensitivities. By fusing transitory sensitivity analysis with more traditional time‐local branching, pathway, and cluster analyses, we develop an algorithm for efficient automatic mechanism analysis. This automatic mechanism analysis at a time point is able to identify the reactions a target is most sensitive to using transitory sensitivity analysis and then propose hypotheses why the reaction might be sensitive using branching, pathway, and cluster analyses. We implement these algorithms within the reaction mechanism simulator (RMS) package, which enables us to report the automatic mechanism analysis results in highly readable text formats and in molecular flux diagrams.

Johnson, Matthew S.↗
DOE OSTI2023

Simulating the Excited-State Dynamics of Polaritons with Ab Initio Multiple Spawning

Over the past decade, there has been a growth of interest in polaritonic chemistry, where the formation of hybrid light–matter states (polaritons) can alter the course of photochemical reactions. These hybrid states are created by strong coupling between molecules and photons in resonant optical cavities and can even occur in the absence of light when the molecule is strongly coupled with the electromagnetic fluctuations of the vacuum field. Here, we present a first-principles model to simulate nonadiabatic dynamics of such polaritonic states inside optical cavities by leveraging graphical processing units (GPUs). Our first implementation of this model is specialized for a single molecule coupled to a single-photon mode confined inside the optical cavity but with any number of excited states computed using complete active space configuration interaction (CASCI) and a Jaynes–Cummings-type Hamiltonian. Using this model, we have simulated the excited-state dynamics of a single salicylideneaniline (SA) molecule strongly coupled to a cavity photon with the ab initio multiple spawning (AIMS) method. We demonstrate how the branching ratios of the photodeactivation pathways for this molecule can be manipulated by coupling to the cavity. We also show how one can stop the photoreaction from happening inside of an optical cavity. Finally, we also investigate cavity-based control of the ordering of two excited states (one optically bright and the other optically dark) inside a cavity for a set of molecules, where the dark and bright states are close in energy.

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

Search for Baryon-Number-Violating Processes in B − Decays to the Ξ ¯ c 0 Λ ¯ c − Final State

We report the results of the first search for B − decays to the Ξ ¯ c 0 Λ ¯ c − final state using 711 fb − 1 of data collected at the ϒ ( 4 S ) resonance with the Belle detector at the KEKB asymmetric-energy e + e − collider. The results are interpreted in terms of both direct baryon-number-violating B − decay and Ξ c 0 − Ξ ¯ c 0 oscillations which follow the standard model decay B − → Ξ c 0 Λ ¯ c − . We observe no evidence for baryon number violation and set the 95% confidence-level upper limits on the ratio of baryon-number-violating and standard model branching fractions B ( B − → Ξ ¯ c 0 Λ ¯ c − ) / B ( B − → Ξ c 0 Λ ¯ c − ) to be < 2.7 % and on the effective angular frequency of mixing ω in Ξ c 0 − Ξ ¯ c 0 oscillations to be < 0.76 ps − 1 (equivalent to τ mix > 1.3 ps ). Published by the American Physical Society 2024

Physics↗
DOE OSTI2020

Investigation of Hierarchical Structure Formation in Merocyanine Photovoltaics

Merocyanines (MCs) are a versatile class of small-molecule dyes. Their optoelectronic properties are easily tunable by chemically controlling their donor–acceptor strength, and their structural properties can be tuned by simple side-chain substitution. This paper demonstrates a novel series of MCs featuring an indoline donor with varying hydrocarbon side-chain length (from 6 to 12 carbons) and a tert-butyl-thiazole acceptor, labeled InTBT. Bulk heterojunction organic photovoltaics are fabricated with a [6,6]-phenyl-C 61 -butyric acid methyl ester (PCBM) acceptor and characterized. Films composed of I8TBT:PCBM and I9TBT:PCBM produced the highest power conversion efficiency of 4.5%, which suggests that the morphology is optimized by controlling the side-chain length. Hierarchical structure formation in InTBT:PCBM films is studied using grazing incidence X-ray diffraction (GIXRD), small-angle neutron scattering (SANS), and atomic force microscopy (AFM). When mixed with PCBM, InTBTs with ≤8 side-chain carbons form pure crystalline domains, while InTBTs with ≥9 side-chain carbons mix well with PCBM. SANS demonstrates that increasing side-chain length increases the InTBT-rich domain size. In addition, a branched hexyl–dodecyl side-chain IHDTBT:PCBM film was studied and found to exhibit the worst-performance organic photovoltaic (OPV) device. The large-branched side chain inhibited mixing between IHDTBT and PCBM resulting in large segregated phases.

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

Search for heavy neutral Higgs bosons A and H in the $t\bar{t}$Z channel in proton-proton collisions at 13 TeV

A direct search for new heavy neutral Higgs bosons A and H in the $t\bar{t}$Z channel is presented, targeting the process pp → A → ZH with H → $t\bar{t}$. For the first time, the channel with decays of the Z boson to muons or electrons in association with all-hadronic decays of the $t\bar{t}$ system is targeted. The analysis uses proton-proton collision data collected at the CERN LHC with the CMS experiment at $\sqrt{𝑠}$ = 13 TeV, which correspond to an integrated luminosity of 138 fb −1 . No signal is observed. Upper limits on the product of the cross section and branching fractions are derived for narrow resonances A and H with masses up to 2100 and 2000 GeV, respectively, assuming A boson production through gluon fusion. The results are also interpreted within two-Higgs-doublet models, where A and H are CP-odd and CP-even states, respectively, complementing and substantially extending the reach of previous searches.

2HDM↗