Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Sim”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Automated Interpretation and Extraction of Topographic Information from Time of Flight Secondary Ion Mass Spectrometry Data

Time of flight secondary ion mass spectrometry (ToF SIMS) is one of the most powerful characterization tools allowing imaging of the chemical properties of various systems and materials. It allows precise studies of the chemical composition with sub-100-nm lateral and nanometer depth spatial resolution. However, comprehensive interpretation of ToF SIMS results is challengeable, because of the data volume and its multidimensionality. Furthermore, investigation of the samples with pronounced topographical features are complicated by the spectral shift. In this work we developed approach for the comprehensive ToF SIMS data interpretation based on the data analytics and automated extraction of the sample topography based on time of flight shift. We further applied this approach to investigate correlation between biological function and chemical composition in Arabidopsis roots.

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

Effects of neoclassical tearing modes and toroidal field ripple on lost alpha power in the SPARC tokamak

Here, using the SPIRAL Monte Carlo, full particle-orbit simulation code, we investigate the effects of neoclassical tearing modes (NTMs) and toroidal field (TF) ripple on alpha power losses during steady-state operation of the SPARC primary reference discharge. Model perturbations for TF ripple and the $m/n = 2/1$ and $3/2$ NTMs with exaggerated widths selected based on an H-mode plasma approaching thermal quench are added to a simulated SPARC magnetic equilibrium through which marker particles are tracked. The 3/2 and 2/1 NTMs are located at $\rho_{\mathrm{pol}} \sim 0.76$ and $\rho_\mathrm{pol} \sim 0.86$ respectively, well positioned to increase alpha particle transport into and within an outer lossy region of the plasma beyond $\rho_\mathrm{pol} \sim 0.8$ where over 95% of lost alpha particles are born. Total alpha power losses are shown to increase modestly from 1.73% lost at a minimum to 2.34% lost at a maximum, and alpha particle surface power densities form localized hotspots on the first-wall near the low-field side midplane due to NTMs and TF ripple. We establish a conservative upper limit for first-wall alpha surface power densities on a toroidally symmetric wall for typical, flattop operation and motivate the consideration of NTMs in the design of three dimensional limiter surfaces for SPARC.

alpha power losses↗
DOE OSTI2019

CHAMP cosmic rays

In this work, we study interactions of cosmological relics, $X$, of mass $m$ and electric charge $qe$ in the galaxy, including thermalization with the interstellar medium, diffusion through inhomogeneous magnetic fields and Fermi acceleration by supernova shock waves. We find that for $m \mathop{}_{\textstyle \sim}^{\textstyle <} 10^{10} q \; ~{\rm GeV}$, there is a large flux of accelerated $X$ in the disk today, with a momentum distribution $\propto 1/p^{2.5}$ extending to $(\beta p)_{max} \sim 5 \times10^4 q \; ~{\rm GeV}$. Even though acceleration in supernova shocks is efficient, ejecting $X$ from the galaxy, $X$ are continually replenished by diffusion into the disk from the halo or confinement region. For $m \mathop{}_{\textstyle \sim}^{\textstyle >} 10^{10} q \; ~{\rm GeV}$, $X$ cannot be accelerated above the escape velocity within the lifetime of the shock. The accelerated $X$ form a component of cosmic rays that can easily reach underground detectors, as well as deposit energies above thresholds, enhancing signals in various experiments. Additionally, we find that nuclear/electron recoil experiments place very stringent bounds on $X$ at low $q$; for example, $X$ as dark matter is excluded for $q$ above $10^{-9}$ for any $m$. For larger $q$, stringent bounds on the fraction of dark matter that can be $X$ are set by Cherenkov and ionization detectors. Nevertheless, very small $q$ is highly motivated by the kinetic mixing portal, and we identify regions of $(m,q)$ that can be probed by future experiments.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Reliable photometric membership (RPM) of galaxies in clusters – I. A machine learning method and its performance in the local universe

ABSTRACT We introduce a new method to determine galaxy cluster membership based solely on photometric properties. We adopt a machine learning approach to recover a cluster membership probability from galaxy photometric parameters and finally derive a membership classification. After testing several machine learning techniques (such as stochastic gradient boosting, model averaged neural network and k-nearest neighbours), we found the support vector machine algorithm to perform better when applied to our data. Our training and validation data are from the Sloan Digital Sky Survey main sample. Hence, to be complete to $M_r^* + 3$, we limit our work to 30 clusters with $z$phot-cl ≤ 0.045. Masses (M200) are larger than $\sim 0.6\times 10^{14} \, \mathrm{M}_{\odot }$ (most above $3\times 10^{14} \, \mathrm{M}_{\odot }$). Our results are derived taking in account all galaxies in the line of sight of each cluster, with no photometric redshift cuts or background corrections. Our method is non-parametric, making no assumptions on the number density or luminosity profiles of galaxies in clusters. Our approach delivers extremely accurate results (completeness, C $\sim 92{\rm{ per\ cent}}$ and purity, P $\sim 87{\rm{ per\ cent}}$) within R200, so that we named our code reliable photometric membership. We discuss possible dependencies on magnitude, colour, and cluster mass. Finally, we present some applications of our method, stressing its impact to galaxy evolution and cosmological studies based on future large-scale surveys, such as eROSITA, EUCLID, and LSST.

Lopes, Paulo A. A.↗
DOE OSTI2021

Evidence for sub-Chandrasekhar Type Ia supernovae from the last major merger

ABSTRACT We investigate the contribution of sub-Chandrasekhar mass Type Ia supernovae to the chemical enrichment of the Gaia Sausage galaxy, the progenitor of a significant merger event in the early life of the Milky Way. Using a combination of data from Nissen & Schuster, the GALactic Archaeology with HERMES (GALAH) Data Release 3 [with 1D non-local thermal equilibrium (NLTE) abundance corrections], and the Apache Point Observatory Galactic Evolution Experiment (APOGEE) Data Release 16, we fit analytic chemical evolution models to a nine-dimensional chemical abundance space (Fe, Mg, Si, Ca, Cr, Mn, Ni, Cu, and Zn) in particular focusing on the iron-peak elements, Mn and Ni. We find that low [Mn/Fe] $\sim -0.15\, \mathrm{dex}$ and low [Ni/Fe] $\sim -0.3\, \mathrm{dex}$ Type Ia yields are required to explain the observed trends beyond the [α/Fe] knee of the Gaia Sausage (approximately at [Fe/H] $=-1.4\, \mathrm{dex}$). Comparison to theoretical yield calculations indicates a significant contribution from sub-Chandrasekhar mass Type Ia supernovae in this system (from ${\sim} 60\, \mathrm{per\, cent}$ to $100\, \mathrm{per\, cent}$ depending on the theoretical model with an additional ${\pm} 10\, \mathrm{per\, cent}$ systematic from NLTE corrections). We compare to results from other Local Group environments including dwarf spheroidal galaxies, the Magellanic Clouds, and the Milky Way’s bulge, finding the Type Ia [Mn/Fe] yield must be metallicity dependent. Our results suggest that sub-Chandrasekhar mass channels are a significant, perhaps even dominant, contribution to Type Ia supernovae in metal-poor systems, whilst more metal-rich systems could be explained by metallicity-dependent sub-Chandrasekhar mass yields, possibly with additional progenitor mass variation related to star formation history, or an increased contribution from Chandrasekhar mass channels at higher metallicity.

Sanders, Jason L.↗
DOE OSTI2021

Formation and evolution of protostellar accretion discs – II. From 3D simulation to a simple semi-analytic model of Class 0/I discs

Here we use a 3D radiative non-ideal magnetohydrodynamic simulation to investigate the formation and evolution of a young protostellar disc from a magnetized pre-stellar core. The simulation covers the first ${\sim }10\, {\rm kyr}$ after protostar formation and shows a massive, weakly magnetized disc with radius that initially grows and then saturates at ${\sim }30\, {\rm au}$. The disc is gravitationally unstable with prominent large-amplitude spiral arms. We use our simulation results and a series of physical arguments to construct a predictive and quantitative physical picture of Class 0/I protostellar disc evolution from several aspects, including (i) the angular-momentum redistribution in the disc, self-regulated by gravitational instability to make most of the disc marginally unstable; (ii) the thermal profile of the disc, well-approximated by a balance between radiative cooling and accretion heating; and (iii) the magnetic-field strength and magnetic-braking rate inside the disc, regulated by non-ideal magnetic diffusion. Using these physical insights, we build a simple 1D semi-analytic model of disc evolution. We show that this 1D model, when coupled to a computationally inexpensive simulation for the evolution of the surrounding pseudo-disc, can be used reliably to predict disc evolution in the Class 0/I phase. The predicted long-term evolution of disc size, which saturates at ${\sim }30\, {\rm au}$ and eventually shrinks, is consistent with a recent observational survey of Class 0/I discs. Such hierarchical modelling of disc evolution circumvents the computational difficulty of tracing disc evolution through Class 0/I phase with direct, numerically converged simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Stability and pulsation of the first dark stars

ABSTRACT The first bright objects to form in the Universe might not have been ‘ordinary’ fusion-powered stars, but ‘dark stars’ (DSs) powered by the annihilation of dark matter (DM) in the form of weakly interacting massive particles (WIMPs). If discovered, DSs can provide a unique laboratory to test DM models. DSs are born with a mass of the order of M⊙ and may grow to a few million solar masses; in this work we investigate the properties of early DSs with masses up to $\sim \! 1000 \, \mathrm{ M}_\odot$, fueled by WIMPS weighing 100 GeV. We improve the previous implementation of the DM energy source into the stellar evolution code mesa. We show that the growth of DSs is not limited by astrophysical effects: DSs up to $\sim \!1000 \, \mathrm{ M}_{\odot }$ exhibit no dynamical instabilities; DSs are not subject to mass-loss driven by super-Eddington winds. We test the assumption of previous work that the injected energy per WIMP annihilation is constant throughout the star; relaxing this assumption does not change the properties of the DSs. Furthermore, we study DS pulsations, for the first time investigating non-adiabatic pulsation modes, using the linear pulsation code gyre. We find that acoustic modes in DSs of masses smaller than $\sim \! 200 \, \mathrm{ M}_\odot$ are excited by the κ − γ and γ mechanism in layers where hydrogen or helium is (partially) ionized. Moreover, we show that the mass-loss rates potentially induced by pulsations are negligible compared to the accretion rates.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Production efficiencies of sound waves in the intracluster medium driven by AGN jets

Feedback from active galactic nuclei (AGNs) is believed to be the most promising solution to the cooling flow problem in cool-core clusters. Dissipation of sound waves is considered as one of the possible heating mechanisms; however, its relative contribution to heating remains unclear. To estimate the energy budget for heating, we perform 3D hydrodynamic simulations of AGN jet injections in a Perseus-like cluster and quantify the amount of energy stored in the forms of weak shocks and waves. For this work, e find that, for a single jet injection with typical parameters in cool-core clusters, ${\sim}9{{\ \rm per\ cent}}$ of the total jet energy is stored in compressional waves (including both shocks and waves). However, due to the destructive effects among randomly phased waves as well as the dissipation of shock energies, in our simulations including self-regulated AGN feedback, no more than $3{{\ \rm per\ cent}}$ of the total injected energy goes into compressional waves. We further separate the energy contribution from shocks and waves and find that, for a single outburst, the shocks can only contribute to ${\sim}20{-}30{{\ \rm per\ cent}}$ of the total compressional energy in the inner radii and quickly dissipate away. In the self-regulated case where shocks are repeatedly generated, shocks completely dominate over sound waves in the inner region and can still provide ${\sim}40{-}50{{\ \rm per\ cent}}$ of the total compressional energy at outer radii. Our results suggest that the production of sound waves is not as efficient as what was previously found, and thus sound wave dissipation may be a subdominant source of heating in cool-core clusters.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

The rate of extreme coronal line emitters in the Baryon Oscillation Spectroscopic Survey LOWZ sample

Extreme coronal line emitters (ECLEs) are a rare class of galaxy that exhibit strong, high-ionization iron coronal emission lines in their spectra. In some cases, these lines are transient and may be the result of tidal disruption event (TDEs). To test this connection, we calculate the rate of variable ECLEs (vECLEs) at redshift $\sim 0.3$. We search for ECLEs in the Baryon Oscillation Spectroscopic Survey (BOSS) LOWZ sample and discover two candidate ECLEs. Using follow-up spectra from the Dark Energy Spectroscopic Instrument and Gemini Multi-Object Spectrograph, and mid-infrared observations from the Wide-field Infrared Survey Explorer, we determine that one of these galaxies is a vECLE. Using this galaxy, we calculate the galaxy-normalized vECLE rate at redshift $\sim 0.3$ to be $R_\mathrm{G}=1.6~^{+3.8}_{-1.4}\times 10^{-6}~\mathrm{galaxy}^{-1}~\mathrm{yr}^{-1}$and the mass-normalized rate to be $R_\mathrm{M}=7~^{+16}_{-6}\times 10^{-18}~\mathrm{M_\odot ^{-1}}~\mathrm{yr}^{-1}$. This is then converted to a volumetric rate of $R_\mathrm{V}=1.8~^{+4.5}_{-1.5}\times 10^{-9}~\mathrm{Mpc}^{-3}~\mathrm{yr}^{-1}$. Formally, the LOWZ vECLE rates are $2 \!-\! 4$ times lower than the rates calculated from the Sloan Digital Sky Survey Legacy sample at redshift $\sim 0.1$. However, given the large uncertainties on both measurements, they are consistent with each other at $1\sigma$. Both the galaxy-normalized and volumetric rates are one to two orders of magnitude lower than TDE rates from the literature, consistent with vECLEs being caused by $5 \!-\! 20$ per cent of all TDEs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Statistical estimates of the binary properties of rotational variables

ABSTRACT We present a model to estimate the average primary masses, companion mass ranges, the inclination limit for recognizing a rotational variable, and the primary mass spreads for populations of binary stars. The model fits a population’s binary mass function distribution and allows for a probability that some mass functions are incorrectly estimated. Using tests with synthetic data, we assess the model’s sensitivity to each parameter, finding that we are most sensitive to the average primary mass and the minimum companion mass, with less sensitivity to the inclination limit and little to no sensitivity to the primary mass spread. We apply the model to five populations of binary spotted rotational variables identified in ASAS-SN, computing their binary mass functions using RV data from APOGEE. Their average primary mass estimates are consistent with our expectations based on their CMD locations ($\sim 0.75 \, {\rm M}_{\odot }$ for lower main sequence primaries and $\sim 0.9$–$1.2 \, {\rm M}_{\odot }$ for RS CVn and sub-subgiants). Their companion mass range estimates allow companion masses down to $M_2/M_1\simeq 0.1$, although the main sequence population may have a higher minimum mass fraction ($\sim 0.4$). We see weak evidence of an inclination limit $\gtrsim 50^{\circ }$ for the main sequence and sub-subgiant groups and no evidence of an inclination limit in the other groups. No groups show strong evidence for a preferred primary mass spread. We conclude by demonstrating that the approach will provide significantly better estimates of the primary mass and the minimum mass ratio and reasonable sensitivity to the inclination limit with 10 times as many systems.

Phillips, Anya (ORCID:000900051914974X)↗
DOE OSTI2018

The formation of supermassive black holes from Population III.1 seeds. I. Cosmic formation histories and clustering properties

We calculate cosmic distributions in space and time of the formation sites of the first, ‘Pop III.1’ stars, exploring a model in which these are the progenitors of all supermassive black holes (SMBHs), seen in the centres of most large galaxies. Pop III.1 stars are defined to form from primordial composition gas in dark matter minihaloes with |${\sim }10^6\, \mathrm{M}_\odot$| that are isolated from neighbouring astrophysical sources by a given isolation distance, d_iso. We assume Pop III.1 sources are seeds of SMBHs, based on protostellar support by dark matter annihilation heating that allows them to accrete a large fraction of their minihalo gas, i.e. |${\sim }10^5\, \mathrm{M}_\odot$|⁠. Exploring d_iso from 10 to |$100\, \rm {kpc}$| (proper distances), we predict the redshift evolution of Pop III.1 source and SMBH remnant number densities. The local, |$z$| = 0 density of SMBHs constrains |$d_{\rm {iso}}\lesssim 100\, \rm {kpc}$| (i.e. |$3\, \rm {Mpc}$| comoving distance at |$z$| ≃ 30). In our simulated (⁠|${\sim }60\, \rm {Mpc})^3$| comoving volume, Pop III.1 stars start forming just after |$z$| = 40. Their formation is largely complete by |$z$| ≃ 25–20 for d_iso = 100–|$50\, \rm {kpc}$|⁠. We follow source evolution to |$z$| = 10, by which point most SMBHs reside in haloes with |${\gtrsim }10^8\, \mathrm{M}_\odot$|⁠. Over this period, there is relatively limited merging of SMBHs for these values of d_iso. We also predict SMBH clustering properties at |$z$| = 10: feedback suppression of neighbouring sources leads to relatively flat angular correlation functions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Cosmological constraints from galaxy–lensing cross-correlations using BOSS galaxies with SDSS and CMB lensing

ABSTRACT We present cosmological parameter constraints based on a joint modelling of galaxy–lensing cross-correlations and galaxy clustering measurements in the SDSS, marginalizing over small-scale modelling uncertainties using mock galaxy catalogues, without explicit modelling of galaxy bias. We show that our modelling method is robust to the impact of different choices for how galaxies occupy dark matter haloes and to the impact of baryonic physics (at the $\sim 2{{\ \rm per\ cent}}$ level in cosmological parameters) and test for the impact of covariance on the likelihood analysis and of the survey window function on the theory computations. Applying our results to the measurements using galaxy samples from BOSS and lensing measurements using shear from SDSS galaxies and CMB lensing from Planck, with conservative scale cuts, we obtain $S_8\equiv \left(\frac{\sigma _8}{0.8228}\right)^{0.8}\left(\frac{\Omega _\mathrm{ m}}{0.307}\right)^{0.6}=0.85\pm 0.05$ (stat.) using LOWZ × SDSS galaxy lensing, and S8 = 0.91 ± 0.1 (stat.) using combination of LOWZ and CMASS × Planck CMB lensing. We estimate the systematic uncertainty in the galaxy–galaxy lensing measurements to be $\sim 6{{\ \rm per\ cent}}$ (dominated by photometric redshift uncertainties) and in the galaxy–CMB lensing measurements to be $\sim 3{{\ \rm per\ cent}}$, from small-scale modelling uncertainties including baryonic physics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

The XRISM first-light observation: Velocity structure and thermal properties of the supernova remnant N 132D

Abstract We present an initial analysis of the X-Ray Imaging and Spectroscopy Mission (XRISM) first-light observation of the supernova remnant (SNR) N 132D in the Large Magellanic Cloud. The Resolve microcalorimeter has obtained the first high-resolution spectrum in the 1.6–10 keV band, which contains K-shell emission lines of Si, S, Ar, Ca, and Fe. We find that the Si and S lines are relatively narrow, with a broadening represented by a Gaussian-like velocity dispersion of $\sigma _v \sim 450$ km s$^{-1}$. However, the Fe He$\alpha$ lines are substantially broadened with $\sigma _v \sim 1670$ km s$^{-1}$. This broadening can be explained by a combination of the thermal Doppler effect due to the high ion temperature and the kinematic Doppler effect due to the SNR expansion. Assuming that the Fe He$\alpha$ emission originates predominantly from the supernova ejecta, we estimate the reverse shock velocity at the time when the bulk of the Fe ejecta were shock heated to be $-1000 \lesssim V_{\rm rs}$ (km s$^{-1}$) $\lesssim 3300$ (in the observer frame). We also find that Fe Ly$\alpha$ emission is redshifted with a bulk velocity of $\sim 890$ km s$^{-1}$, substantially larger than the radial velocity of the local interstellar medium surrounding N 132D. These results demonstrate that high-resolution X-ray spectroscopy is capable of providing constraints on the evolutionary stage, geometry, and velocity distribution of SNRs.

individual objects↗
DOE OSTI2025

In-orbit performance of the soft X-ray imaging telescope Xtend aboard XRISM

Here, we present a summary of the in-orbit performance of the soft X-ray imaging telescope Xtend onboard the X-Ray Imaging and Spectroscopy Mission (XRISM), based on in-flight observation data, including first-light celestial objects, calibration sources, and results from the cross-calibration campaign with other currently operating X-ray observatories. XRISM/Xtend has a large field of view of ${38{^{\prime }_{.}}5}$ $\times$ ${38{^{\prime }_{.}}5}$, covering an energy range of 0.4–13 keV, as demonstrated by the first-light observation of the galaxy cluster Abell 2319. It also features an energy resolution of 170–180 eV at 6 keV, which meets the mission requirement and enables us to resolve He-like and H-like Fe K$\alpha$ lines. Throughout the observation during the performance verification phase, we confirm that two issues identified in the Soft X-ray Imager (SXI) onboard the previous Hitomi mission—light leakage and crosstalk events—are addressed and suppressed in the case of Xtend. A joint cross-calibration observation of the bright quasar 3C 273 results in an effective area measured to be $\sim$420 cm$^{2}$ at1.5 keV and $\sim$310 cm$^{2}$ at 6.0 keV, which matches values obtained in ground tests. We also continuously monitor the health of Xtend by analyzing overclocking data, calibration source spectra, and day-Earth observations; the readout noise is stable and low, and contamination is negligible even one year after launch. A low background level compared with other major X-ray instruments onboard satellites, combined with the largest grasp ($\Omega _{\rm eff}\sim 60$ cm$^2$ deg$^2$) of Xtend, will not only support Resolve analysis, but also enable significant scientific results on its own. This includes near-future follow-up observations and transient searches in the context of time-domain and multi-messenger astrophysics.

instrumentation: detectors: Xtend↗
DOE OSTI2025

XRISM view of a stellar flare: High-resolution Fe K spectra of HR 1099, an RS CVn-type star

A high-resolution X-ray spectroscopic observation was made of the RS CVn-type binary star HR 1099 using the Resolve instrument onboard XRISM for its calibration purposes. During the $\sim$400 ks telescope time covering 1.5 binary orbit, a flare lasting for $\sim$100 ks was observed with a released X-ray radiation energy of ${\sim }10^{34}$ erg, making it the first stellar flare ever observed with an X-ray microcalorimeter spectrometer. The flare peak count rate is 6.4 times higher than that in quiescence and is distinguished clearly in time thanks to the long telescope time. Many emission lines were detected in the 1.7–10 keV range both in the flare and quiescent phases. Using the high spectral resolution of Resolve in the Fe K band (6.5–7.0 keV), we resolved the inner-shell lines of Fe XIX-XXVI as well as the outer-shell lines of Fe XXV-XXVI . These lines have peaks in the contribution functions at different temperatures over a wide range, allowing us to construct the differential emission measure (DEM) distribution over the electron temperature of 1–10 keV (roughly 10–100 MK) based only on Fe lines, thus without an assumption of the elemental abundance. The reconstructed DEM has a bimodal distribution, and only the hotter component increased during the flare. The elemental abundance was derived based on the DEM distribution thus constructed. A significant abundance increase was observed during the flare for Ca and Fe, which are some of the elements with the lowest first ionization potential among those analyzed, but not for Si, S, and Ar. This behavior is seen in some giant solar flares and the present result is a clear example in stellar flares.

Astronomy and AstroPhysics↗
DOE OSTI2025

Lie-algebraic classical simulations for quantum computing

The classical simulation of quantum dynamics plays an important role in our understanding of quantum complexity and in the development of quantum technologies. Efficient techniques such as those based on the Gottesman-Knill theorem for Clifford circuits, tensor networks for low entanglement-generating circuits, or Wick's theorem for fermionic Gaussian states have become central tools in quantum computing. In this work, we contribute to this body of knowledge by presenting a framework for classical simulations, dubbed “𝔤-sim”, which is based on the underlying Lie algebraic structure of the dynamical process. When the dimension of the algebra grows at most polynomially in the system size, there exist observables for which the simulation is efficient. Indeed, we show that 𝔤-sim enables new regimes for classical simulations, is able to deal with certain forms of noise in the evolution, as well as can be used to tackle several paradigmatic variational and nonvariational quantum computing tasks. For the former, we perform Lie-algebraic simulations to train and optimize parametrized quantum circuits (thus effectively showing that some variational models can be dequantized), design enhanced parameter initialization strategies, solve tasks of quantum circuit synthesis, and train a quantum-phase classifier. For the latter, we report large-scale noiseless and noisy simulations on benchmark problems. By comparing the limitations of 𝔤-sim and certain Wick's theorem-based simulations, we find that the two methods become inefficient for different types of states or observables, hinting at the existence of distinct, nonequivalent resources for classical simulation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Effect of crab cavity high order modes on the coupled-bunch stability of High-Luminosity Large Hadron Collider

High frequency High Order Modes can significantly affect the coupled-bunch stability in a circular accelerator. With a large enough shunt impedance they may drive a coupled-bunch transverse instability. We have developed an analytical model and implemented it in the NHT Vlasov solver. The results show that the instability is characterized by the excitation of many azimuthal intra-bunch modes, which have similar growth rates, that make the traditional remedies such as a flat resistive feedback and chromaticity inefficient in suppressing it. For the High Luminosity Large Hadron Collider this could result in a significant increase of the stabilizing Landau octupole current, up to $\sim 100$ A ($\sim 20\%$ of the maximum available current). In order to limit the increase below 10 A ($\sim 2\%$), the transverse shunt impedance has to be kept below 1 M$\Omega$/m.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Beyond universal behavior in the one-dimensional chain with random nearest-neighbor hopping

Here, we study the one-dimensional nearest-neighbor tight-binding model of electrons with independently distributed random hopping and no on-site potential (i.e., off-diagonal disorder with particle-hole symmetry, leading to sublattice symmetry, for each realization). For nonsingular distributions of the hopping, it is known that the model exhibits a universal, singular behavior of the density of states $ρ(E) ~ 1/|E\: \text{ln}^3 |E||$ and of the localization length $ξ(E) ~ |\text{ln}|E||$, near the band center $\textit{E}$ = 0 . (This singular behavior is also applicable to random $\textit{XY}$ and Heisenberg spin chains; it was first obtained by Dyson for a specific random harmonic oscillator chain.) Simultaneously, the state at$\textit{E}$ = 0 shows a universal, subexponential decay at large distances $\sim \text{exp}[–\sqrt{r/r_0}]$. In this study, we consider singular, but normalizable, distributions of hopping, whose behavior at small $\textit{t}$ is of the form $\sim 1/[t \text{ln}^{λ+1}(1/t)]$, characterized by a single, continuously tunable parameter λ > 0. We find, using a combination of analytic and numerical methods, that while the universal result applies for λ > 2, it no longer holds in the interval 0 < λ < 2. In particular, we find that the form of the density of states singularity is enhanced (relative to the Dyson result) in a continuous manner depending on the nonuniversal parameter λ; simultaneously, the localization length shows a less divergent form at low energies and ceases to diverge below λ = 1. For λ < 2, the fall-off of the $\textit{E}$ = 0 state at large distances also deviates from the universal result and is of the form $\sim\text{exp}[–(r/r_0)^{1/λ}]$ , which decays faster than an exponential for λ < 1.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗