SEARCH · Search NASA

Results for “Harmonic numbers”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 181 records · Page 10

All sources NASA NTRS DOE OSTI

DOE OSTI2022

E2 rotational invariants of 0$^{+}_{1}$ and 2$^{+}_{1}$ states for 106 Cd: The emergence of collective rotation

The collective structure of 106 Cd is elucidated by multi-step Coulomb excitation of a 3.849 MeV/A beam of 106 Cd on a 1.1 mg/cm 2 208 Pb target using GRETINA-CHICO2 at ATLAS. Fourteen E2 matrix elements were obtained. The nucleus 106 Cd is a prime example of emergent collectivity that possesses a simple structure: it is free of complexity caused by shape coexistence and has a small, but collectively active number of valence nucleons. This work follows in a long and currently active quest to answer the fundamental question of the origin of nuclear collectivity and deformation, notably in the cadmium isotopes. The results are discussed in terms of phenomenological models, the shell model, and Kumar-Cline sums of E2 matrix elements. The < 0$^{+}_{2}$ ||E2||2$^{+}_{1}$ > matrix element is determined for the first time, providing a total, converged measure of the electric quadrupole strength, < Q 2 >, of the first-excited 2$^{+}_{1}$ level relative to the 0$^{+}_{1}$ ground state, which does not show an increase as expected of harmonic and anharmonic vibrations. Strong evidence for triaxial shapes in weakly collective nuclei is indicated; collective vibrations are excluded. This is contrary to the only other cadmium result of this kind in 114 Cd by C. Fahlander et al., Nucl. Phys. A485, 327 (1988), which is complicated by low-lying shape coexistence near midshell.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2020

Density-Matrix Based Extended Lagrangian Born–Oppenheimer Molecular Dynamics

Extended Lagrangian Born–Oppenheimer molecular dynamics [ Phys. Rev. Lett. 2008, 100, 123004] is presented for Hartree–Fock theory, where the extended electronic degrees of freedom are represented by a density matrix, including fractional occupation numbers at elevated electronic temperatures. In contrast to regular direct Born–Oppenheimer molecular dynamics simulations, no iterative self-consistent field optimization is required prior to the force evaluations. To sample regions of the potential energy landscape where the gap is small or vanishing, which leads to particular convergence problems in regular direct Born–Oppenheimer molecular dynamics simulations, an adaptive integration scheme for the extended electronic degrees of freedom is presented. The integration scheme is based on a tunable, low-rank approximation of a fourth-order kernel, K, that determines the metric tensor, T ≡ K T K, used in the extended harmonic oscillator of the Lagrangian that generates the dynamics of the electronic degrees of freedom. Here, the formulation and algorithms provide a general guide to implement extended Lagrangian Born–Oppenheimer molecular dynamics for quantum chemistry, density functional theory, and semiempirical methods using a density matrix formalism.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Protection and Restoration Solutions to Reliable and Resilient Integration of Grid-connected PV Installations and Distributed Energy Resources: Design, Testbed, Proof of Work and Impact Studies (Final Report)

To gain a better understanding of the complex transients in a utility grid with a large number of solar PV installations coupled by PWM inverters during the protection and restoration period of the power grid, in this project a kW-level experimental system with dominating inverter-based resources and a hardware-in-loop simulation system with transmission and distribution models, as well as several SEL protective relays have been built and used. The major research findings of the project are summarized below in two perspectives: Experimental research: A self-organized ultra-high frequency solitary waveform is discovered and demonstrated in the hardware experimental system. Despite the familiarities, such a solitary waveform is distinct from any harmonics, transient, or resonance waves in that, it is 1) non-dispersive over time and space, 2) not associated with any active source or the linear superposition of sources, 3) half-cycle asymmetric, 4) not responding to filters or change of system characteristic resonance frequency, 5) ubiquitous as it occurs simultaneously everywhere in the system from DC supply, power lines, and the utility grid, and 6) explosive through tripping the protection or damaging susceptible devices or circuits. Analytical study: The nature of such a new waveform and an analytical explanation of its formation are studied based on related physics and non-linear science principles, with the following highlights: 1) such a new waveform follows the solutions to the Non-linear Schrödinger equation, so the classic linear perturbation theory is unable to explain or predict such a unique waveform as confirmed with the research team of RTDS; 2) the critical condition of occurrence of such a waveform is derived, which indicates that the breakings of solitary wave is a system synchronous issue between the utility grid in the 60Hz phasor domain and duty-ratio modulation of DC sources in the inverter switching frequency domain; 3) such a waveform carries energy mass so it could be detrimental; 4) such a waveform is deceiving as it is not readily detectable in the energy propagation direction by the primary protection equipment, while it is detrimental in the perpendicular direction of energy propagation, i.e. voltage direction. Thus, it has more likely challenges to the second primary equipment and devices, particularly at the weakest link and point such as aging insulation and inappropriate setting of susceptible devices. Therefore, such a waveform can be easily ignored in the aftermath investigation. The intellectual merits: The experimental discovery reveals certain unfamiliar transient phenomena that could challenge the integration of large-scale grid-connected solar PV installation during the group ride-through period of solar PV installations, commissioning of large-solar farms, or dramatic change of solar radiation conditions. Particularly, the research findings suggest that the occurrence of the unfamiliar transient phenomena is uniquely associated with the inverter-based solar PV installation, which is less likely to happen for rotary energy systems. The physics and non-linear science-based research work laid down an analytical path to the challenging transient stability problems including those that have been observed currently and future calls. Both experimental and analytical research results explain the limitation of many current research effects including some DoE research undertakings as well as possible solutions. The broader impacts: The research outcome of the project advances the understanding of the possible transient problems for the integration of inverter-based solar PV installation. It also highlights the theoretical and technical barriers for applying the conventional theory and technology to design and implement countermeasures against their adverse impacts of the large-scale solar PV installation and operation on grid reliability and security. More broadly, the research outcomes have shed some light on the open, fundamental challenges in the integration of large-scale solar photovoltaic energy into current and next-generation power grids across the country, and worldwide. The advanced physics and non-linear science-based analysis open up a path to harmonization of the renewables for societal energy needs via building a resilient and sustainable electric power infrastructure.

14 SOLAR ENERGY↗

DOE OSTI2020

Impacts of swinging door lossy compression of synchrophasor data

As the number of operational Phasor Measurement Units (PMUs) is increasing, compression is becoming crucial for efficiently transferring and storing enormous amount of synchrophasor data continuously measured from a power system. Some utilities in the U.S. have already adopted a nonlinear compression technique, called Swinging Door Algorithm (SDA), assuming compressed data meets technical and regulatory requirements without resulting in any degradation of the measurements. This paper studies impacts of this compression algorithm on field measured and simulated PMU data, and presents results from spectral analyses of raw and decompressed measurements. The results show blind application of the SDA compression alters frequency content in the measurements and introduces false harmonics and aliasing. These can lead to false conclusions being drawn from the measurements about the power system condition.

42 ENGINEERING↗

DOE OSTI2026

Leveraging large language models to address data scarcity in machine learning for graphene synthesis

Machine learning in experimental materials science faces significant challenges due to the scarcity of data, which are costly and time-consuming to generate, particularly when relying on in-house experiments. Literature data mining offers a potential solution but introduces issues like mixed data quality, inconsistent formats, and non-uniform reporting of synthesis parameters, resulting in partially missing and heterogeneous features across the dataset. Here, we propose data imputation and feature engineering methods that employ pre-trained large language models (LLMs) to enhance machine learning performance on scarce, heterogeneous datasets, demonstrated on graphene CVD synthesis data and the ML-HydPARK hydrogen storage dataset. GPT models perform data imputation via tailored prompting and semantic normalization of inconsistently reported features through embeddings, for example, to harmonize the complex nomenclature of CVD substrates. Beyond yielding more diverse and richer feature representations than traditional methods such as K-nearest neighbors (KNN) and Multivariate Imputation by Chained Equations (MICE), LLM-based data imputation is evaluated against dataset characteristics and prompting strategies. We vary the level of autonomy granted to the LLM, from generic prompting that leverages pre-trained knowledge for autonomous data generation to data-informed prompting that constrains outputs using target-specific information, and demonstrate which level of autonomy yields superior imputation performance across datasets and feature types. The proposed data engineering methods markedly improve downstream performance; for example, in graphene layer number classification using a support vector machine (SVM), binary accuracy increases from 39% to 65% and ternary accuracy from 52% to 72%. Fine-tuning experiments on both datasets show that combining our proposed LLM-based data imputation and feature encoding methods with numerical machine learning predictors outperforms standalone fine-tuned LLM predictors in data-scarce settings. The proposed strategies emphasize data enhancement techniques rather than refining learning architectures or regularizing loss functions, offering a broadly applicable framework for improving machine learning performance on scarce, inhomogeneous datasets.

Chemical vapor deposition↗

DOE OSTI2023

Magnon orbital angular momentum of ferromagnetic honeycomb and zigzag lattice models

By expanding the gauge 𝜆 𝑛 ⁡(𝐤) for magnon band 𝑛 in harmonics of momentum 𝐤=(𝑘,𝜙), we demonstrate that the only observable component of the magnon orbital angular momentum 𝑂 𝑛 ⁡(𝐤) is its angular average over all angles 𝜙, denoted by 𝐹 𝑛 ⁡(𝑘). Although 𝐹 𝑛 ⁡(𝑘) vanishes for antiferromagnetic honeycomb and zigzag (0<𝐽 1 <𝐽 2 ) lattices, it is nonzero for the ferromagnetic (FM) versions of those lattices in the presence of Dzyaloshinskii-Moriya interactions. For a FM zigzag model with equal exchange interactions 𝐽 1⁢𝑥 and 𝐽 1⁢𝑦 along the 𝑥 and 𝑦 axes, the magnon bands are degenerate along the boundaries of the Brillouin zone with 𝑘 𝑥 −𝑘 𝑦 =±𝜋/𝑎 and the Chern numbers 𝐶 𝑛 are not well defined. However, a revised model with 𝐽 1⁢𝑦 ≠𝐽 1⁢𝑥 lifts those degeneracies and produces well-defined Chern numbers of 𝐶𝑛=±1 for the two magnon bands. When 𝐽 1⁢𝑦 =𝐽 1⁢𝑥 , the thermal conductivity 𝜅 𝑥⁢𝑦⁡ (𝑇) of the FM zigzag lattice is largest for 𝐽2/𝐽1>6 but is still about four times smaller than that of the FM honeycomb lattice at high temperatures. Due to the removal of band degeneracies, 𝜅𝑥⁢𝑦⁡(𝑇) is slightly enhanced when 𝐽 1⁢𝑦 ≠𝐽 1⁢𝑥 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

Equilibrium carbon isotope fractionation factors of hydrocarbons: Semi-empirical force-field method

Here, we have calculated the reduced partition function ratios for carbon isotopes (β-factor) of 67 hydrocarbons (alkanes, alkenes, alkynes, cycloalkanes, and aromatics), including their 267 single-substituted isotopomers. The calculations were performed using the harmonic oscillator – rigid rotator model and the Urey / Bigeleisen-Mayer approach. Normal frequencies of molecular vibrations of the isotopologues were calculated from the molecular structures data and valence field force constants stored in the Light-handled Elucidation of Vibrations (LEV) database elaborated by Dr. Gribov and his colleagues in Russia. The LEV database is constructed from experimental spectroscopic and structural data by solving inverse problems of molecular vibrations. The LEV is internally consistent and the most comprehensive database available to date. The β-factor were calculated in the temperature range of 200–800 K with a 10 K step. Our calculations predict that the β-factors increase with an increasing number of C atoms within the same groups of hydrocarbons (e.g., alkanes). Our calculations also show a general descending order of 13 C enrichments among the different groups of hydrocarbons: cycloalkanes, aromatics, alkenes (double bonds) and isoalkanes, alkanes, alkynes (triple bonds). Position-specific, intramolecular isotope effects within hydrocarbons are determined by the β-factors of C in the different functional groups in the order: quaternary (C), tertiary (methine, CH), secondary (methylene, CH 2 ), primary (methyl, CH 3 ) and double bond (C=C), saturated bond (C - C), triple bond (HC≡). Our calculations on bulk and position-specific carbon isotope β-factors of the hydrocarbons, which are generally consistent with very limited ab initio calculations in the literature, are internally consistent and the most comprehensive to date for future applications to position-specific isotope geochemistry of hydrocarbons.

58 GEOSCIENCES↗

DOE OSTI2023

Measurement of 2D density profiles using a second-harmonic, dispersion interferometer

A second-harmonic, dispersion interferometer is used to image large-area (≃5 cm 2 ) plasma-jet and gas-jet density profiles. Achromatic telescopes magnify the diameters of the primary-laser beam (1064 nm) and its second-harmonic (532 nm) before probing the sample and de-magnify the beam diameters after the sample, where the primary beam transfers its phase change to a second, second-harmonic beam, allowing the sample’s dispersive-phase change to be measured between two, orthogonally polarized second harmonic beams. The telescopes produce an azimuthally symmetric, dispersive-phase shift in the sample + background phase-change image and in the background phase-change image, which is removed by digital subtraction. The interferometer’s performance was verified using standard-optical components as dispersive elements (BK7 lenses and wedge plates), resolving a minimum, phase-change sensitivity of ΔΦ min ≳ 15 mrad and spatial resolution of Δx res ≃ 100 μm. The phase change produced by unknown-density objects (a pulsed-plasma-jet and a pulsed-gas-jet) was measured, and their data were used to recover the original, 2D density profiles using an inverse Abel transform: peak-number density, N gas ≃ 6 × 10 20 cm –3 and N e ≃ 5 × 10 16 cm –3 ; line-integrated density, ∫N gas dl ~ 2 × 10 19 cm –2 and ∫N e dl ~ 1 × 10 16 cm –2 . The techniques and methods developed here are scalable to even larger probe-beam diameters and frame-capture rates, leading to a diagnostic capability that is well-suited for applications involving the real-time measurement of density.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Cybersecurity Certification Recommendations for Interconnected Grid Edge Devices and Inverter Based Resources

Escalating deployment of PV and grid-edge devices on the distribution grid has increased the sustainability and efficiency of the electric grid. However, the increasing number of distributed energy resources (DERs) deployed creates a heightened cyber-physical interdependency on the distribution grid and thus creates more vectors for cyber-attacks to exploit through information and communication technology (ICT) systems and networks. For example, control signal packets can be modified, intercepted, or corrupted due to vulnerabilities in communication protocols used by microgrid controllers and grid edge devices for power control. Therefore, to mitigate and prevent cyber-attacks on grid edge devices and the inverter-based resources connected to the distribution grid, the U.S. Department of Solar Energy Technologies Office (SETO) awarded funding to the National Renewable Energy Laboratory and Sandia National Laboratory (SNL) to research, develop, and harmonize cybersecurity standards for Photovoltaic (PV) systems and for other kinds of DERs. To help develop a standard for DER cybersecurity, NREL established certification recommendations and test cases, in consensus with the solar industry and UL, for ensuring intrinsic design security for DERs. These recommendations were developed to bolster the cybersecure functionalities such as TLS, MAC, CRL, session resumption/renegotiation, and password, system, and service security management within the DER devices. The proposed test cases verify authentication, authorization, confidentiality, and data integrity for data and communications of DERs that use Transmission Control Protocol/Internet Protocol (TCP/IP). They were also developed to protect DER communications from eavesdropping, replay, man-in-the-middle, denial of service (DoS), spoofing through security certificates, least-privilege violation, and brute-force credentials. This report, which has been validated and reviewed by UL, expands upon those test cases to provide DER cybersecurity certification recommendations which increase DER resiliency and help to mitigate cyber-attacks. UL's collaboration with NREL and approval of this document will accelerate the adoption of a UL standard for DER cybersecurity.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2021

Correlations of azimuthal anisotropy Fourier harmonics with subevent cumulants in p Pb collisions at s N N = 8.16 TeV

Event-by-event long-range correlations of azimuthal anisotropy Fourier coefficients ($v_n$) in 8.16 TeV p Pb data, collected by the CMS experiment at the CERN Large Hadron Collider, are extracted using a subevent four-particle cumulant technique applied to very low multiplicity events. Each combination of four charged particles is selected from either two, three, or four distinct subevent regions of a pseudorapidity range from -2.4 to 2.4 of the CMS tracker, and with transverse momentum between 0.3 and 3.0 GeV. Using the subevent cumulant technique, correlations between $v_n$ of different orders are measured as functions of particle multiplicity and compared to the standard cumulant method without subevents over a wide event multiplicity range. At high multiplicities, the $v_2$ and $v_3$ coefficients exhibit an anticorrelation; this behavior is observed consistently using various methods. The $v_2$ and $v_4$ correlation strength is found to depend on the number of subevents used in the calculation. As the event multiplicity decreases, the results from different subevent methods diverge because of different contributions of noncollective or few-particle correlations. Correlations extracted with the four-subevent method exhibit a tendency to diminish monotonically toward the lowest multiplicity region (about 20 charged tracks) investigated. These findings extend previous studies to a significantly lower event multiplicity range and establish the evidence for the onset of long-range collective multiparticle correlations in small system collisions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Ion cyclotron emission in Maxwellian plasmas

Destabilization of ion cyclotron waves—waves with frequencies close to ion cyclotron harmonics—in inhomogeneous plasmas with Maxwell velocity distribution is considered. A new mechanism of destabilization of these waves is found, whereas the known interaction between cyclotron waves and drift waves is shown to be hardly able to lead to instabilities in realistic tokamak plasmas. Our finding is that the resonance wave–particle interaction in the presence of temperature gradient can change the diamagnetic drift frequency in such a way that the destabilizing influence of the diamagnetic drift exceeds Landau damping. This occurs when the energy of resonance particles is sufficiently high, which is always the case due to the infinite tail of the Maxwell distribution. Particles with higher energies provide a larger ratio of drive to damping, but the number of resonance particles and, thus, instability growth rates are exponentially small when particle energy is very high. Therefore, only moderately suprathermal particles can lead to observable instabilities. A condition for instabilities driven by these particles is obtained. Destabilization of electrostatic cyclotron waves and ordinary cyclotron waves is studied.

Kolesnichenko, Ya. I. (ORCID:0000000323219681)↗

DOE OSTI2020

Modal analysis of the behavior of inertial particles in turbulence subjected to Stokes drag

The clustering of small heavy inertial particles subjected to Stokes drag in turbulence, which is the concern of present study, is known to be minimal at small and large Stokes number and substantial at St = $\textit{O}$(1). This nonmonotonic trend, which has been shown computationally and experimentally, is yet to be explained analytically. In this study, we obtain an analytical expression for the Lyapunov exponents that quantitatively predicts this trend. The sum of the exponents, which is the normalized rate of change of the signed volume of a small cloud of particles, is correctly predicted to be negative and positive at small and large Stokes numbers, respectively, asymptoting to $\textit{τQ}$ as $\textit{τ}$ → 0 and $τ^{–1/2}|Q|^{1/4}$ as $\textit{τ} → ∞$, where $\textit{τ}$ is the particle relaxation time and $\textit{Q(τ)}$ is the difference between the norm of the rotation- and strain-rate tensors computed along the particle trajectory, which remains negative at all $\textit{τ}$ in turbulent flows. Additionally, the trajectory crossing is predicted only in hyperbolic flows where Q < 0 for sufficiently inertial particles with a $\textit{τ}$ that scales with $|Q|^{–1/2}$. Following the onset of crossovers, a transition from clustering to dispersion is predicted correctly. We show these behaviors are not unique to three-dimensional isotropic turbulence and can be reproduced closely by a one-dimensional mono-harmonic flow, which appears as a fundamental canonical problem in the study of particle clustering. Analysis of this one-dimensional canonical flow shows that the rate of clustering, quantified as the product of the Lyapunov exponent and particle relaxation time, is bounded by –1/2, behaving with extreme nonlinearity in the hyperbolic flows and always remaining positive in the elliptic flows. Finally, this model problem also confirms that dispersion in a hyperbolic flow will occur only if particles' trajectory cross. These findings, which stem from our analysis, are corroborated by the direct numerical simulations.

42 ENGINEERING↗

DOE OSTI2023

A statistical mechanics framework for polymer chain scission, based on the concepts of distorted bond potential and asymptotic matching

To design increasingly tough, resilient, and fatigue-resistant elastomers and hydrogels, the relationship between controllable network parameters at the molecular level (bond type, non-uniform chain length, entanglement density, etc.) to macroscopic quantities that govern damage and failure must be established. Many of the most successful constitutive models for elastomers have been rooted in statistical mechanical treatments of polymer chains. Typically, such constitutive models have used variants of the freely jointed chain model with rigid links. However, since the free energy state of a polymer chain is dominated by enthalpic bond distortion effects as the chain approaches its rupture point, bond extensibility ought to be accounted for if the model is intended to capture chain rupture. To that end, a new bond potential is supplemented to the freely jointed chain model (as derived in the u FJC framework of Buche and Silberstein (2021) and Buche et al. (2022)), which we have extended to yield a tractable, closed-form model of single chain behavior that should be amenable to continuum-level constitutive model development. Inspired by the asymptotically matched u FJC model response in both the low/intermediate chain force and high chain force regimes, a simple, quasi-polynomial bond potential energy function is derived. This bond potential exhibits harmonic behavior near the equilibrium state and anharmonic behavior for large bond stretches tending to a characteristic energy plateau (akin to the Lennard-Jones and Morse bond potentials). Using this bond potential, approximate yet highly-accurate analytical functions for bond stretch and chain force dependent upon chain stretch are established. Then, using this polymer chain model, a stochastic thermal fluctuation-driven chain rupture framework is developed. This framework is based upon a force-modified tilted bond potential that accounts for distortional bond potential energy, allowing for the derivation and subsequent calculation of the dissipated chain scission energy. Here, the cases of rate-dependent and rate-independent scission are accounted for throughout the rupture framework. The impact of Kuhn segment number on chain rupture behavior is also investigated. The model is fit to single chain mechanical response data collected from atomic force microscopy tensile tests for validation and to glean deeper insight into the molecular physics taking place. Due to their analytical nature, this polymer chain model and the associated rupture framework can in the future be implemented in finite element models accounting for fracture and fatigue in polydisperse elastomer networks.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Universal Jaynes-Cummings Control of an Oscillator

The Jaynes-Cummings (JC) interaction-the coherent exchange of excitations between a two-level system and a harmonic oscillator-is one of the fundamental interactions of quantum optics, realized across platforms such as cavity quantum electrodynamics, trapped ions, mechanical resonators, and superconducting circuits. Although JC interactions and qubit rotations form a universal gate set for oscillator control, practical implementations have not been demonstrated. Here we develop and experimentally demonstrate universal JC-based oscillator control by compiling arbitrary unitary gates into sequences of JC interactions and qubit rotations. In our experiment, the oscillator is realized using a mode of a high quality factor microwave cavity and the ancilla qubit using a superconducting transmon circuit, with the JC interaction implemented by a sideband interaction enabled by the Josephson nonlinearity. The native gates are constructed to be closed below a chosen cutoff photon number, encoding a qudit with suppressed leakage errors, while ancilla relaxation errors are detectable. We further find that the dispersive shift serves as a compilation resource that reduces circuit depths. We demonstrate universal qudit control and implement a single-qutrit gate set with a mean post-selected process fidelity of 96%, as well as ququart and ququint shift gates. These results establish Jaynes-Cummings control as a practical route to universal oscillator control, enabling programmable bosonic processors across a variety of quantum platforms.

Huang, Jordan [Rutgers U., Piscataway]↗

DOE OSTI2021

Baseline predictions of elliptic flow and fluctuations for the RHIC Beam Energy Scan using response coefficients

Currently the RHIC Beam Energy Scan is exploring a new region of the Quantum Chromodynamic phase diagram at large baryon densities that approaches nuclear astrophysics regimes. This provides an opportunity to study relativistic hydrodynamics in a regime where the net conserved charges of baryon number, strangeness, and electric charge play a role, which will significantly change the theoretical approach to simulating the baryon-dense Quark-Gluon Plasma. In this work, we detail many of the important changes needed to adapt both initial conditions and the medium to baryon-rich matter. Then, we make baseline predictions for the elliptical flow and fluctuations based on extrapolating the physics at LHC and top RHIC energies to support future analyses of where and how the new baryon-dense physics causes these extrapolations to break down. First we compare eccentricities across beam energies, exploring their underlying assumptions; we find the the extrapolated initial state is predicted to be nearly identical to that at AuAu $\sqrt{s_{NN}}$ = 200 GeV. Then the final flow harmonic predictions are based on linear+cubic response. We discuss preliminary STAR results in order to determine the implications that they have for linear+cubic response coefficients at the lowest beam energy of AuAu $\sqrt{s_{NN}}$ = 7 GeV.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

Third international challenge to model the medium- to long-range transport of radioxenon to four Comprehensive Nuclear-Test-Ban Treaty monitoring stations

In 2015 and 2016, atmospheric transport modeling challenges were conducted in the context of the Comprehensive Nuclear-Test-Ban Treaty (CTBT) verification, however, with a more limited scope with respect to emission inventories, simulation period and number of relevant samples (i.e., those above the Minimum Detectable Concentration (MDC)) involved. Therefore, a more comprehensive atmospheric transport modeling challenge was organized in 2019. Stack release data of Xe-133 were provided by the Institut National des Radioéléments/IRE (Belgium) and the Canadian Nuclear Laboratories/CNL (Canada) and accounted for in the simulations over a three (mandatory) or six (optional) months period. Best estimate emissions of additional facilities (radiopharmaceutical production and nuclear research facilities, commercial reactors or relevant research reactors) of the Northern Hemisphere were included as well. Model results were compared with observed atmospheric activity concentrations at four International Monitoring System (IMS) stations located in Europe and North America with overall considerable influence of IRE and/or CNL emissions for evaluation of the participants’ runs. Participants were prompted to work with controlled and harmonized model set-ups to make runs more comparable, but also to increase diversity. It was found that using the stack emissions of IRE and CNL with daily resolution does not lead to better results than disaggregating annual emissions of these two facilities taken from the literature if an overall score for all stations covering all valid observed samples is considered. A moderate benefit of roughly 10% is visible in statistical scores for samples influenced by IRE and/or CNL to at least 50% and there can be considerable benefit for individual samples. Effects of transport errors, not properly characterized remaining emitters and long IMS sampling times (12–24 h) undoubtedly are in contrast to and reduce the benefit of high-quality IRE and CNL stack data. Complementary best estimates for remaining emitters push the scores up by 18% compared to just considering IRE and CNL emissions alone. Despite the efforts undertaken the full multi-model ensemble built is highly redundant. An ensemble based on a few arbitrary runs is sufficient to model the Xe-133 background at the stations investigated. The effective ensemble size is below five. An optimized ensemble at each station has on average slightly higher skill compared to the full ensemble. However, the improvement (maximum of 20% and minimum of 3% in RMSE) in skill is likely being too small for being exploited for an independent period.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

MARS-Q modeling of kink-peeling instabilities in DIII-D QH-mode plasma

In quiescent H-mode (QH-mode) regime, edge harmonic oscillations (EHOs) are believed to provide necessary radial transport to prevent occurrence of large edge localized modes. A systematic modeling study is performed here on the low- n EHOs in a DIII-D QH-mode plasma (Chen et al 2016 Nucl. Fusion 56 076011), by utilizing the MARS-Q code (Liu et al 2013 Phys. Plasmas 20 042503). Both the n = 1 and n = 2 instabilities are found to be strongly localized near the plasma edge, exhibiting the edge-peeling characteristics. The DIII-D resistive wall is found to have minor effects on these instabilities. The plasma resistivity is found to strongly modify the mode growth rate. Assuming the Spitzer model for the plasma resistivity, the computed mode growth rate scales as S sup>-1/3 with S being the Lundquist number. Toroidal flow of the plasma slightly stabilizes these edge localized kink-peeling modes. Drift kinetic effects all have a destabilization effect on these modes. Non-perturbative magneto-hydrodynamic-kinetic hybrid computations find that the drift kinetic effects associated with thermal particle species push the peak location of the eigenmode radially inward but still in the pedestal region. The modeled plasma temperature and density fluctuations in the plasma edge region, as well as the poloidal magnetic field perturbations along both the low and high field sides of the plasma surface, are in good agreement with experimental measurements. Finally, the quasi-linear initial value simulations find a strong non-linear interplay between the kink-peeling instability and the toroidal flow near the plasma edge. The combined effect of the damping of the flow amplitude and change of the edge flow shear is found to be the stabilizing factor for the kink-peeling mode, leading to the mode saturation and thus EHOs.

Physics↗

DOE OSTI2025

Extreme Longitudinal Compression of Optimized Beams for MEV Ultrafast Electron Diffraction (Final Technical Report)

We worked out the design of a high repetition rate MeV energy ultrafast electron diffraction instrument based on the existing Cornell photoinjector, which can readily be applied to the presented findings. This example is a blueprint of other similarly arranged UED setups. Using particle tracking simulations in conjunction with multiobjective genetic algorithm optimization, we explored the smallest bunch lengths, emittance, and probe spot sizes achievable. As two limits, we defined stroboscopic conditions (with single electrons per pulse) and operation with 10 5 electrons per bunch which may be suitable for single-shot diffraction images. In the stroboscopic case, the flexibility provided by the many cavity bunching and acceleration allows for longitudinal phase space linearization without a higher harmonic field, providing sub-fs bunch lengths at the sample. Given low emittance photoemission conditions, these small bunch lengths can be maintained with probe transverse sizes at the single micron (1 μm) scale and below. In the case of 10 5 electrons per pulse, we simulated state-of-the-art 5D brightness conditions: rms bunch lengths of 10 fs with 3-nm normalized emittances, while permitting repetition rates as high as 1.3 GHz. We showed that in conjunction with collimating apertures, a novel focusing scheme achieves very high-quality emittance compensation for the central core of the beam composing 40% of particles, for a resulting beam size of 5 μm (rms). Finally, to aid in the design of new SRF-based ultrafast electron diffraction machines, we simulated the trade-off between the number of cavities used and achievable bunch length and emittance. In the longitudinal dimension, we made use of the fact that MeV UED requires much lower energy than the 15-MeV maxi mum energy of Cornell’s CBETA injector, and we may therefore use several of the SRF cavities for bunch length compression. In practice, we used a genetic optimization algorithm to choose the phases and amplitudes of the cavities appropriately for optimal bunching. In the zero space charge case, we found that bunching and acceleration are distributed across the six cavities in a way that produces a linearizing effect. And we showed that the ultimate bunch length can be limited by time-of-flight differences arising from transverse size and transverse momentum spread. The space charge code developed and used for this development is now permanent part of the Bmad accelerator simulation code and has already contributed to other developments, e.g., for the EIC electron cooler design.

43 PARTICLE ACCELERATORS↗