Search NASA⌕ Search

SEARCH · Search NASA

Results for “oscillation cases”

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

SPT-3G D1: Axion early dark energy with CMB experiments and DESI observations

We present the most up-to-date constraints on axion early dark energy (AEDE) from cosmic microwave background (CMB) and baryon acoustic oscillation (BAO) measurements. In particular, we assess the impact of data from ground-based CMB experiments, the South Pole Telescope (SPT) and the Atacama Cosmology Telescope (ACT)—both with and without Planck—on constraints on AEDE. We also highlight the impact that BAO information from the Dark Energy Spectroscopic Instrument (DESI) has on these constraints. From CMB data alone, we do not find statistically significant evidence for the presence of AEDE, and we find only moderate reduction in the Hubble tension. From the latest SPT data alone, we find the maximal fractional contribution of AEDE to the cosmic energy budget is f EDE < 0.12 at 95% confidence level (CL), and the Hubble tension between the SPT and Supernovae, H0, for the Equation of State of Dark Energy (SH0ES) results is reduced to the 2.3 σ level. When combining the latest SPT, ACT, and Planck datasets, we find f EDE < 0.070 at 95% CL and the Hubble tension at the 3.6 σ level. In contrast, adding DESI data to the CMB datasets results in mild preference for AEDE and, in some cases, non-negligible reduction in the Hubble tension. From SPT + DESI , we find f EDE = 0.08 1 - 0.052 + 0.037 at 68% CL, and the Hubble tension reduces to 1.5 σ . From the combination of DESI with all three CMB experiments, we get f EDE = 0.05 5 - 0.047 + 0.024 at 68% CL and a weak preference for AEDE over Λ CDM . This data combination, in turn, reduces the Hubble tension to 2.6 σ . We highlight that this shift in parameters when adding the DESI dataset is a manifestation of the discrepancy currently present between DESI and CMB experiments in the concordance model Λ CDM .

Khalife, A. R. [Paris, Inst. Astrophys.] (ORCID:00↗
DOE OSTI2020

Transport, correlations, and chaos in a classical disordered anharmonic chain

Here, we explore transport properties in a disordered nonlinear chain of classical harmonic oscillators, and thereby identify a regime exhibiting behavior analogous to that seen in quantum many-body-localized systems. Through extensive numerical simulations of this system connected at its ends to heat baths at different temperatures, we computed the heat current and the temperature profile in the nonequilibrium steady state as a function of system size $\textit{N}$, disorder strength Δ, and temperature $\textit{T}$. The conductivity $κ_N$, obtained for finite length ($\textit{N}$), saturates to a value $κ_∞$ > 0 in the large $\textit{N}$ limit, for all values of disorder strength Δ and temperature $\textit{T}$ > 0. We show evidence that for any Δ > 0 the conductivity goes to zero faster than any power of $\textit{T}$ in the ($\textit{T}$/Δ) → 0 limit, and find that the form $κ_∞$ ~ $e^{{–B|\text{ln}(CΔ/T)|}^{3}}$ fits our data. This form has earlier been suggested by a theory based on the dynamics of multioscillator chaotic islands. The finite-size effect can be $κ_N < κ_∞$ due to boundary resistance when the bulk conductivity is high (the weak disorder case), or $κ_N > κ_∞$ due to direct bath-to-bath coupling through bulk localized modes when the bulk is weakly conducting (the strong disorder case). We also present results on equilibrium dynamical correlation functions and on the role of chaos on transport properties. Finally, we explore the differences in the growth and propagation of chaos in the weak and strong chaos regimes by studying the classical version of the out-of-time-ordered commutator.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Where, when, and why: Occurrence of fast-pairwise collective neutrino oscillation in three-dimensional core-collapse supernova models

Fast-pairwise collective neutrino oscillation represents a key uncertainty in the theory of core-collapse supernova (CCSN). Despite the potentially significant impact on CCSN dynamics, it is usually neglected in numerical models of CCSN because of the formidable technical difficulties of self-consistently incorporating this physics. In this paper, we investigate the prospects for the occurrence of fast flavor conversion by diagnosing electron neutrino lepton number (ELN) crossing in more than a dozen state-of-the-art three-dimensional CCSN models. ELN crossings provide a necessary condition for triggering flavor conversion. Although only zeroth and first angular moments are available from the simulations, our new method enables us to look into the angular distributions of neutrinos in momentum space and provide accurate insight into ELN crossings. Our analysis suggests that fast flavor conversion generally occurs in the postshock region of CCSNe, and that explosive models provide more favorable conditions for the flavor conversion than failed CCSNe. Here we also find that there are both common and progenitor-dependent characteristics. Classifying ELN crossings into two types, we analyze the generation mechanism of each case by scrutinizing the neutrino radiation field and matter interactions. We find key ingredients of CCSN dynamics driving the ELN crossings: proto-neutron star convection, asymmetric neutrino emission, neutrino absorptions, and scatterings. This study suggests that we need to accommodate fast flavor conversions in realistic CCSN models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Intrepid MCMC: Metropolis-Hastings with exploration

In engineering examples, one often encounters the need to sample from unnormalized distributions with complex shapes that may also be implicitly defined through a physical or numerical simulation model, making it computationally expensive to evaluate the associated density function. For such cases, MCMC has proven to be an invaluable tool. Random-walk Metropolis Methods (also known as Metropolis-Hastings (MH)), in particular, are highly popular for their simplicity, flexibility, and ease of implementation. However, most MH algorithms suffer from significant limitations when attempting to sample from distributions with multiple modes (particularly disconnected ones). Here, in this paper, we present Intrepid MCMC - a novel MH scheme that utilizes a simple coordinate transformation to significantly improve the mode-finding ability and convergence rate to the target distribution of random-walk Markov chains while retaining most of the simplicity of the vanilla MH paradigm. Through multiple examples, we showcase the improvement in the performance of Intrepid MCMC over vanilla MH for a wide variety of target distribution shapes. We also provide an analysis of the mixing behavior of the Intrepid Markov chain, as well as the efficiency of our algorithm for increasing dimensions. A thorough discussion is presented on the practical implementation of the Intrepid MCMC algorithm. Finally, its utility is highlighted through a Bayesian parameter inference problem for a two-degree-of-freedom oscillator under free vibration.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2023

The r -modes of slowly rotating, stratified neutron stars

ABSTRACT The only r-modes that exist in a globally barotropic, rotating, Newtonian star are the fundamental l = |m| solutions, where l and m are the indices of the spherical harmonic $Y_l^m$ that describe the mode’s angular dependence. This is in stark contrast to a stellar model that is non-barotropic throughout its interior, which hosts all the l ≥ |m| perturbations including radial overtones. In reality, neutron stars are stratified with locally barotropic regions. Therefore, we explore how stratification alters a star’s ability to support r-modes. We consider the globally stratified case and examine the behaviour of the modes as the star gets close to barotropicity. In this limit, we find that all but the fundamental l = |m| perturbations change character and become generic inertial modes. Restricting the analysis to l = |m| perturbations, we develop the r-mode equations in order to consider stellar models that exhibit local barotropicity. Our results for such models show that the r-mode overtones diverge and join the inertial modes. In order to see which r-modes persist and retain their character in realistic neutron stars, these calculations will need to be brought into full general relativity.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

An extended Vlasov–Fokker–Planck approach for kinetic simulations of laser plasmas

Vlasov–Fokker–Planck simulation codes occupy an important niche in modeling laser-produced plasmas, since they are well suited to studying the effect of collisions on electron kinetic phenomena, especially energy transport. One of the most important elements of energy transport is the absorption of laser light by the plasma; however, simulating this in detail requires resolving oscillations of the laser light, whose characteristic timescale is orders of magnitude shorter than the simulation time needed to study transport physics. For this reason, most Vlasov–Fokker–Planck codes used to study electron transport in laser plasmas rely on simplified models of the laser–plasma coupling. Their underlying assumptions nominally preclude their use for modeling laser light having short-scale structure in space or time, such as broadband lasers. In this work, we derive a more general computational framework suitable for arbitrarily structured laser fields. Furthermore, our approach is based on an extended set of Vlasov–Fokker–Planck equations that separately solve for the low- and high-frequency plasma response. We implement these extended Vlasov–Fokker–Planck equations in the spherical harmonic code K2 and demonstrate the performance of the method on several laser absorption test problems, with particular attention to the judicious selection of time steps, time integrators, and spherical harmonic truncation, according to the intensity and spectrum of the laser light under consideration. Comparison with the widely used Langdon absorption operator shows the Langdon operator performs remarkably well for predicting laser heating in the simple cases considered here, even in situations that would seem to violate its underlying assumptions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Critical role of structural order in bipolar redox-active molecules for organic redox flow batteries

Bipolar redox-active molecules (BRMs) have been suggested as a means to address crossover related issues in all-organic redox flow batteries (RFBs). In such species, electron donors (anolytes) and electron acceptors (catholytes) are linked by a chain, and the same chemical composition is used in both anode and cathode compartments to reduce chemical gradients and unbalanced ion transport. The resultant RFBs resemble the aqueous vanadium RFBs, but offer greater design flexibility and potentially more favorable electrochemical and physicochemical properties. Yet the complex trade-offs in these properties have been a complication. Here we developed a fundamental rubric to uncover the likely origins in the performance metrics of these BRMs that are salient for their use in RFBs. Methylene linked phenothiazine and phthalimide moieties were employed as model BRMs and variations in their properties and molecular conformations were evaluated by systematically varying the methylene linker. Our results revealed that even minimalistic changes in the linker length resulted in dramatic oscillations in the solubility and stability. Using crystallography, quantum chemistry, and molecular dynamics, we clarify that the seemingly fluctuating behavior is due to (1) the inter- and intra-molecular charge transfer between the donor and acceptor through bond or/and through space; (2) the formation of distinctive packing/clustering motifs. Both interactions are strongly dependent on the molecular conformation. As these structural factors modify the electronic structure of both electroactive functional groups and do not change monotonically, the structure–property response in BRMs is often complex. Finally, accordingly, these behaviors need to be taken into account in developing BRMs and knowledge of these structural factors will allow the more rational design of BRMs for grid-scale energy storage.

25 ENERGY STORAGE↗
DOE OSTI2025

DESI 2024 VII: cosmological constraints from the full-shape modeling of clustering measurements

We present cosmological results from the measurement of clustering of galaxy, quasar and Lyman-α forest tracers from the first year of observations with the Dark Energy Spectroscopic Instrument (DESI Data Release 1). We adopt the full-shape (FS) modeling of the power spectrum, including the effects of redshift-space distortions, in an analysis which has been thoroughly validated in a series of supporting papers as summarised in [1]. We combine the full-shape information with DESI's DR1 constraints from the baryon acoustic oscillations (BAO) of these tracers. In the flat ΛCDM cosmological model, DESI (FS+BAO), combined with a baryon density prior from Big Bang Nucleosynthesis and a weak prior on the scalar spectral index, determines matter density to Ω m = 0.2962 ± 0.0095, and the amplitude of mass fluctuations to σ 8 = 0.842 ± 0.034. The addition of the cosmic microwave background (CMB) data tightens these constraints to Ω m = 0.3056 ± 0.0049 and σ 8 = 0.8121 ± 0.0053, while further addition of the joint clustering and lensing analysis from the Dark Energy Survey Year-3 (DESY3) data further improves these measurements, and leads to a 0.4% determination of the Hubble constant, H 0 = (68.40 ± 0.27) km s -1 Mpc -1 . In models with a time-varying dark energy equation of state parametrised by w 0 and w a , combinations of DESI (FS+BAO) with CMB and type Ia supernovae continue to show the preference, previously found in the DESI DR1 BAO analysis, for w 0 > -1 and w a < 0 with similar levels of significance. DESI data, in combination with the CMB, improve the upper limits on the sum of the neutrino masses relative to the case when only the DR1 BAO was available, giving ∑m ν < 0.071 eV at 95% confidence. We finally constrain deviations from general relativity represented by two modified gravity parameters. DESI (FS+BAO) data alone measure the parameter that controls the clustering of massive particles, μ 0 = 0.11 +0.45 -0.54 , in agreement with the zero value predicted by general relativity. The combination of DESI with the CMB and the clustering and lensing analysis from DESY3 constrains both modified-gravity parameters, giving μ 0 = 0.04 ± 0.22 and Σ 0 = 0.044 ± 0.047, again in agreement with general relativity.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

The scalar chemical potential in cosmological collider physics

Non-analyticity in co-moving momenta within the non-Gaussian bispectrum is a distinctive sign of on-shell particle production during inflation, presenting a unique opportunity for the “direct detection” of particles with masses as large as the inflationary Hubble scale ( H ). However, the strength of such non-analyticity ordinarily drops exponentially by a Boltzmann-like factor as masses exceed H . In this paper, we study an exception provided by a dimension-5 derivative coupling of the inflaton to heavy-particle currents, applying it specifically to the case of two real scalars. The operator has a “chemical potential” form, which harnesses the large kinetic energy scale of the inflaton, \( {\overset{\cdot }{\phi}}_0^{1/2}\approx 60H \) , to act as an efficient source of scalar particle production. Derivative couplings of inflaton ensure radiative stability of the slow-roll potential, which in turn maintains (approximate) scale-invariance of the inflationary correlations. We show that a signal not suffering Boltzmann suppression can be obtained in the bispectrum with strength f NL ~ \( \mathcal{O} \) (0 . 01–10) for an extended range of scalar masses \( \lesssim {\overset{\cdot }{\phi}}_0^{1/2} \) , potentially as high as 10 15 GeV, within the sensitivity of upcoming LSS and more futuristic 21-cm experiments. The mechanism does not invoke any particular fine-tuning of parameters or breakdown of perturbation-theoretic control. The leading contribution appears at tree-level , which makes the calculation analytically tractable and removes the loop-suppression as compared to earlier chemical potential studies of non-zero spins. The steady particle production allows us to infer the effective mass of the heavy particles and the chemical potential from the variation in bispectrum oscillations as a function of co-moving momenta. Our analysis sets the stage for generalization to heavy bosons with non-zero spin.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Dark Energy Survey Year 3 results: Measurement of the baryon acoustic oscillations with three-dimensional clustering

The three-dimensional correlation function offers an effective way to summarize the correlation of the large-scale structure even for imaging galaxy surveys. We have applied the projected three-dimensional correlation function, ξ p to measure the baryonic acoustic oscillations (BAO) scale on the first-three years Dark Energy Survey data. The sample consists of about 7 million galaxies in the redshift range 0.6< z p < 1.1 over a footprint of 4108 deg 2 . Our theory modeling includes the impact of realistic true redshift distributions beyond Gaussian photo-z approximation. ξ p is obtained by projecting the three-dimensional correlation to the transverse direction. To increase the signal-to-noise of the measurements, we have considered a Gaussian stacking window function in place of the commonly used top-hat. ξ p is sensitive to D M (z eff )/r s , the ratio between the comoving angular diameter distance and the sound horizon. Using the full sample, D M (z eff )/r s is constrained to be19.00 ± 0.67 (top-hat) and 19.15 ± 0.58 (Gaussian) at z eff = 0.835. The constraint is weaker than the angular correlation w constraint 18.84 ± 0.50), and we trace this to the fact that the BAO signals are heterogeneous across redshift. While ξ p responds to the heterogeneous signals by enlarging the error bar, w can still give a tight bound on D M s in this case. When a homogeneous BAO-signal subsample in the range 0.7 < z p <1.0 (z eff = 0.845) is considered, ξp yields 19.80 ± 0.67 (top-hat) and 19.84 ± 0.53 (Gaussian). The latter is mildly stronger than the w constraint (19.86 ± 0.55). We find that the ξ p results are more sensitive to photo-z because ξ p keeps the three-dimensional clustering information causing it to be more prone to photo-z noise. The Gaussian window gives more robust results than the top-hat as the former is designed to suppress the low signal modes. ξ p and the angular statistics such as w have their own pros and cons, and they serve an important crosscheck with each other.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Consistent lensing and clustering in a low- S 8 Universe with BOSS, DES Year 3, HSC Year 1, and KiDS-1000

ABSTRACT We evaluate the consistency between lensing and clustering based on measurements from Baryon Oscillation Spectroscopic Survey combined with galaxy–galaxy lensing from Dark Energy Survey (DES) Year 3, Hyper Suprime-Cam Subaru Strategic Program (HSC) Year 1, and Kilo-Degree Survey (KiDS)-1000. We find good agreement between these lensing data sets. We model the observations using the Dark Emulator and fit the data at two fixed cosmologies: Planck (S8 = 0.83), and a Lensing cosmology (S8 = 0.76). For a joint analysis limited to large scales, we find that both cosmologies provide an acceptable fit to the data. Full utilization of the higher signal-to-noise small-scale measurements is hindered by uncertainty in the impact of baryon feedback and assembly bias, which we account for with a reasoned theoretical error budget. We incorporate a systematic inconsistency parameter for each redshift bin, A, that decouples the lensing and clustering. With a wide range of scales, we find different results for the consistency between the two cosmologies. Limiting the analysis to the bins for which the impact of the lens sample selection is expected to be minimal, for the Lensing cosmology, the measurements are consistent with A = 1; A = 0.91 ± 0.04 (A = 0.97 ± 0.06) using DES+KiDS (HSC). For the Planck case, we find a discrepancy: A = 0.79 ± 0.03 (A = 0.84 ± 0.05) using DES+KiDS (HSC). We demonstrate that a kinematic Sunyaev–Zeldovich-based estimate for baryonic effects alleviates some of the discrepancy in the Planck cosmology. This analysis demonstrates the statistical power of small-scale measurements; however, caution is still warranted given modelling uncertainties and foreground sample selection effects.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Shear-free, inhomogeneous turbulence in a stably stratified fluid

High-resolution large eddy simulations are conducted of locally forced, shear-free turbulence in the presence of an initially sharp density interface. The simulations are reminiscent of oscillating grid turbulence experiments used to isolate the effect of turbulent diffusion and entrainment from background shear. By simulating such a flow we avoid common challenges of the experiments such as secondary-flow contamination due to sidewall effects and the inevitable interaction of the stratifying agent and forcing region. To address the latter concern, we add a heating term (potential energy sink) to the governing equations in the forcing layer, thereby preventing a heat flux through the source region. This modification sets up a continuous stratification in the mixed layer that is often assumed to be negligible in experiments. Despite this difference, we are able to make meaningful comparisons in terms of the overall entrainment rate, which varies as a power law with a turbulent Richardson number. Two exponents, −2 and −1, are measured depending on the definition of the Richardson number and entrainment rate used. The definition leading to −1 is consistent with most experiments, and we argue it is the superior choice if one is able to measure the relevant quantities. We also verify the self-similar scaling of turbulence velocity and length scales in the homogeneous fluid and propose ‘inner’ and ‘outer’ scalings for the stratified cases based on a local Froude number. The detailed scaling results are useful for turbulence model validation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Inverse bremsstrahlung absorption rate for super-Gaussian electron distribution functions including plasma screening

Here we provide analytic expressions for the effective Coulomb logarithm for inverse bremsstrahlung absorption which predict significant corrections to the Langdon effect and overall absorption rate compared to previous estimates. The calculation of the collisional absorption rate of laser energy in a plasma by the inverse bremsstrahlung mechanism usually makes the approximation of a constant Coulomb logarithm. We dispense with this approximation and instead take into account the velocity dependence of the Coulomb logarithm, leading to a more accurate expression for the absorption rate valid in both classical and quantum conditions. In contrast to previous work, the laser intensity enters into the Coulomb logarithm. In most laser-plasma interactions the electron distribution function is super-Gaussian [Langdon, Phys. Rev. Lett. 44, 575 (1980)], and we find the absorption rate under these conditions is increased by as much as ≈ 30% compared to previous estimates at low density. In many cases of interest the correction to Langdon's predicted reduction in absorption is large; for example at Z = 6 and Te = 400 eV the Langdon prediction for the absorption is in error by a factor of ≈ 2. However, we also account for the additional effect of plasma screening, which predicts a reduction in absorption by a similar amount (up to ≈ 30%). These two effects compete to determine the overall absorption, which may be increased or decreased, depending on the conditions. The corrections can be incorporated into radiation-hydrodynamics simulation codes by replacing the familiar Coulomb logarithm with an analytic expression which depends on the super-Gaussian order “M” and the screening length.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Once-in-a-lifetime encounter models for neutrino media. II. Quasisteady states and miscidynamic flavor evolution

We extended the once-in-a-lifetime encounter (OILE) model to stochastic interactions among neutrinos. As in the original OILE model, the new model reproduces the mean-field behavior of a dense neutrino gas for time 𝑡 ≲ (𝜇⁢𝛾) −1 , where 𝜇 measures the strength of the mean-field neutrino self-interaction potential and is proportional to the neutrino density, and the dimensionless “impact parameter” 𝛾 is a measure of the change in the flavor quantum state of a neutrino during interaction with another neutrino when the wave packets of the two neutrinos overlap. As in the mean-field case, the OILE model with random neutrino velocities experiences kinetic flavor decoherence as the flavor quantum states of the neutrinos diverge from each other. Unlike the mean-field case, however, the OILE model has a “collision term” due to the quantum entanglement among neutrinos. For 𝛾 ≪ 1, this incoherent effect can drive the neutrinos into a quasisteady state that is similar to the collective precession mode in a homogeneous and isotropic neutrino gas in the mean-field approximation. Subsequently, the collision term drives the neutrino gas adiabatically through different quasisteady states and eventually to flavor equilibration. This process is an example of miscidynamic flavor evolution, with the mixing equilibria being the quasisteady precession states.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Development, Implementation, and Verification of Partially-Ionized Collisional Multifluid Plasma Models in Drekar

We consider the development of multifluid models for partially ionized multispecies plasmas. The models are composed of a standard set of five-moment fluid equations for each species plus a description of electromagnetics. The most general model considered utilizes a full set of fluid equations for each charge state of each atomic species, plus a set of fluid equations for electrons. The fluid equations are coupled through source terms describing electromagnetic coupling, ionization, recombination, charge exchange, and elastic scattering collisions in the low-density coronal limit. The form of each of these source terms is described in detail, and references for required rate coefficients are identified for a diverse range of atomic species. Initial efforts have been made to extend these models to incorporate some higher-density collisional effects, including ionization potential depression and three- body recombination. Some reductions of the general multifluid model are considered. First, a reduced multifluid model is derived which averages over all of the charge states (including neutrals) of each atomic species in the general multifluid model. The resulting model maintains full consistency with the general multifluid model from which it is derived by leveraging a quasi-steady-state collisional ionization equilibrium assumption to recover the ionization fractions required to make use of the general collision models. Further reductions are briefly considered to derive certain components of a single-fluid magnetohydrodynamics (MHD) model. In this case, a generalized Ohm's law is obtained, and the standard MHD resistivity is expressed in terms of the collisional models used in the general multifluid model. A number of numerical considerations required to obtain robust implementations of these multifluid models are discussed. First, an algebraic flux correction (AFC) stabilization approach for a continuous Galerkin finite element discretization of the multifluid system is described in which the characteristic speeds used in the stabilization of the fluid systems are synchronized across all species in the model. It is demonstrated that this synchronization is crucial in order to obtain a robust discretization of the multifluid system. Additionally, several different formulations are considered for describing the electromagnetics portion of the multifluid system using nodal continuous Galerkin finite element discretizations. The formulations considered include a parabolic divergence cleaning method and an implicit projection method for the traditional curl formulation of Maxwell's equations, a purely- hyperbolic potential-based formulation of Maxwell's equations, and a mixed hyperbolic-elliptic potential-based formulation of Maxwell's equations. Some advantages and disadvantages of each formulation are explored to compare solution robustness and the ease of use of each formulation. Numerical results are presented to demonstrate the accuracy and robustness of various components of our implementation. Analytic solutions for a spatially homogeneous damped plasma oscillation are derived in order to verify the implementation of the source terms for electromagnetic coupling and elastic collisions between fluid species. Ionization balance as a function of electron temperature is evaluated for several atomic species of interest by comparing to steady-state calculations using various sets of ionization and recombination rate coefficients. Several test problems in one and two spatial dimensions are used to demonstrate the accuracy and robustness of the discretization and stabilization approach for the fluid components of the multifluid system. This includes standard test problems for electrostatic and electromagnetic shock tubes in the two-fluid and ideal shock-MHD limits, a cylindrical diocotron instability, and the GEM challenge magnetic reconnection problem. A one-dimensional simplified prototype of an argon gas puff configuration as deployed on Sandia's Z-machine is used as a demonstration to exercise the full range of capabilities associated with the general multifluid model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Chaotic and integrable magnetic fields in one-dimensional hybrid Vlasov–Maxwell equilibria

The construction of kinetic equilibrium states is important for studying stability and wave propagation in collisionless plasmas. Thus, many studies over the past decades have been focused on calculating Vlasov–Maxwell equilibria using analytical and numerical methods. However, the problem of kinetic equilibrium of hybrid models is less studied, and self-consistent treatments often adopt restrictive assumptions ruling out cases with irregular and chaotic behaviour, although such behaviour is observed in spacecraft observations of space plasmas. In this paper, we develop a one-dimensional (1-D), quasineutral, hybrid Vlasov–Maxwell equilibrium model with kinetic ions and massless fluid electrons and derive associated solutions. The model allows for an electrostatic potential that is expressed in terms of the vector potential components through the quasineutrality condition. The equilibrium states are calculated upon solving an inhomogeneous Beltrami equation that determines the magnetic field, where the inhomogeneous term is the current density of the kinetic ions and the homogeneous term represents the electron current density. We show that the corresponding 1-D system is Hamiltonian, with position playing the role of time, and its trajectories have a regular, periodic behaviour for ion distribution functions that are symmetric in the two conserved particle canonical momenta. For asymmetric distribution functions, the system is nonintegrable, resulting in irregular and chaotic behaviour of the fields. The electron current density can modify the magnetic field phase space structure, inducing orbit trapping and the organization of orbits into large islands of stability. Thus, the electron contribution can be responsible for the emergence of localized electric field structures that induce ion trapping. We also provide a paradigm for the analytical construction of hybrid equilibria using a rotating two-dimensional harmonic oscillator Hamiltonian, enabling the calculation of analytic magnetic fields and the construction of the corresponding distribution functions in terms of Hermite polynomials.

Kaltsas, Dimitrios A. (ORCID:0000000300769015)↗
DOE OSTI2022

C2QA - Bosonic Qiskit

The practical benefits of hybrid quantum information processing hardware that contains continuous-variable objects (bosonic modes such as mechanical or electromagnetic oscillators) in addition to traditional (discrete-variable) qubits have recently been demonstrated by experiments with bosonic codes that reach the break-even point for quantum error correction [1]–[5] and by efficient Gaussian boson sampling simulation of the Franck-Condon spectra of triatomic molecules [6] that is well beyond the capabilities of current qubit-only hardware. The goal of this Co-design Center for Quantum Advantage (C2QA) project is to develop an instruction set architecture (ISA) for hybrid qubit/bosonic mode systems that contains an inventory of the fundamental operations and measurements that are possible in such hardware. The corresponding abstract machine model (AMM) would also contain a description of the appropriate error models associated with the gates, measurements and time evolution of the hardware. This information has been implemented as an extension of IBM Qiskit. IBM Qiskit is an open-source software development toolkit (SDK) for simulating the quantum state of a quantum circuit and for running the same circuits on prototype hardware within the IBM Quantum Experience. We introduce the Bosonic Qiskit software to enable the simulation of hybrid qubit/bosonic systems using the existing Qiskit software development kit [7]. This implementation can be used for simulating new hybrid systems, verifying proposed physical systems, and modeling systems larger than can currently be constructed. We also cover tutorials and example use cases included within the software to study JaynesCummings models, bosonic Hubbard models, plotting Wigner functions and animations, and calculating maximum likelihood estimations using Wigner functions.

quantum, boson, Qiskit↗
DOE OSTI2020

The impact of the fiducial cosmology assumption on BAO distance scale measurements

ABSTRACT Standard analysis pipelines for measurements of Baryon Acoustic Oscillations (BAO) in galaxy surveys make use of a fiducial cosmological model to guide the data compression required to transform from observed redshifts and angles to the measured angular and radial BAO peak positions. In order to remove any dependence on the fiducial cosmology from the results, all models compared to the data should mimic the compression and its dependence on the fiducial model. In practice, approximations are made when testing models: (1) There is assumed to be no residual dependence on the fiducial cosmology after reconstruction, (2) differences in the distance–redshift relationship are assumed to match a linear scaling, and (3) differences in clustering between true and fiducial models are assumed to be removed by the free parameters used to null the non-BAO signal. We test these approximations using the current standard measurement procedure with a set of halo catalogues from the aemulus suite of N-body simulations, which span a range of wCDM cosmological models. We focus on reconstruction of the primordial BAO and locating the BAO. For the range of wCDM cosmologies covered by the aemulus suite, we find no evidence for systematic errors in the measured BAO shift parameters α∥ and α⊥ to $\lt 0.1\%$. However, the measured errors $\sigma _{\alpha _{\parallel }}$ and $\sigma _{\alpha _{\bot }}$ show a notable absolute increase by up to +0.001 and +0.002, respectively, in the case that the fiducial cosmology does not match the truth. These effects on the inferred BAO scale will be important, given the precision of measurements expected from future surveys including DESI, Euclid, and WFIRST.

79 ASTRONOMY AND ASTROPHYSICS↗