Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “MAXED”

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

Quantum critical dynamics of a Josephson junction at the topological transition

We find the admittance Y(ω) of a Josephson junction at or near a topological transition. The dependence of the admittance on frequency and temperature at the critical point is universal and determined by the symmetries of the system. Despite the absence of a spectral gap at the transition, the dissipative response may remain weak at low energies: ReY(ω)∝max(ω,T) 2 . This behavior is strikingly different from the electromagnetic response of a normal metal. Away from the critical point, the scaling functions for the dependence of the admittance on frequency and temperature are controlled by at most two parameters.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Separation of Kondo lattice coherence from crystal electric field in CeIn 3 with Nd substitutions

The Ce m M n In 3m+2n (m=1,2;n=0,1) family has been one of the most studied families of heavy fermion compounds. This family has revealed many interesting low-temperature physics phenomena, like quantum critical points, heavy fermion superconductivity, and non-Fermi liquid behavior, when these materials are exposed to pressure, magnetic fields, and/or chemical substitution. In this work, we provide a thorough investigation of the Ce 1-x Nd x In 3 phase diagram through single crystal synthesis, x-ray diffraction, energy-dispersive spectroscopy, magnetic susceptibility, and electrical resistivity measurements. Previous electrical resistivity measurements on CeIn 3 reveal a broad maximum, T max ~50 K, which has been associated with the Kondo lattice coherence crossover and/or the crystal electric field depopulation effect as the 4f electrons condense from the high-energy quartet down to the ground state doublet. Our findings show that in the most disordered substitution region, x=0.4–0.5, these features disjoin to reveal two distinct broad humps in electrical resistivity measurements. Magnetic susceptibility and electrical resistivity data on Ce 1-x Nd x In 3 also reveal the antiferromagnetic ordering competition between CeIn 3 and NdIn 3 , where the T N of CeIn 3 is linearly suppressed to a critical concentration of x Nd ~0.6. This concentration is slightly lower than what was previously reported in nonmagnetically substituted Ce 1-x La x In 3 . Our magnetic susceptibility measurements and subsequent simulations show that in the CeIn 3 antiferromagnetic regime, x Nd ≤0.4, the Nd ions act as free paramagnets. The large magnitude of the associated paramagnetic response then masks the overlapping antiferromagnetic ordering signature of the Ce ions. Overall our study further sheds light on the underlying crystal electric field and Kondo lattice coherence interactions within the Ce m M n In 3m+2n family and could stimulate further studies of these systems via neutron diffraction or under applied pressure.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Quantum magnetism in the frustrated square lattice oxyhalides YbBi 2 IO 4 and YbBi 2 ClO 4

Square-lattice systems offer a direct route for realizing two-dimensional (2D) quantum magnetism with frustration induced by competing interactions. In this work, the square lattice materials YbBi 2 IO 4 and YbBi 2 ClO 4 were investigated using a combination of magnetization and specific-heat measurements on polycrystalline samples. Specific-heat measurements provide evidence for long-range magnetic order below $T$ N = 0.21 K (0.25 K) for YbBi 2 IO 4 (YbBi 2 ClO 4 ). On the other hand, a rather broad maximum is found in the temperature-dependent magnetic susceptibility, located at $T$ max = 0.33 K (0.38 K) in YbBi 2 IO 4 (YbBi 2 ClO 4 ), consistent with the quasi-2D magnetism expected for the large separation between the magnetic layers. Estimation of the magnetic entropy supports the expected Kramers' doublet ground state for Yb 3+ and the observed paramagnetic behavior is consistent with a well-isolated doublet. Roughly two-thirds of the entropy is consumed above $T$ N , due to a combination of the quasi-2D behavior and magnetic frustration. The impact of frustration is examined from the viewpoint of a simplified $J$ 1 -$J$ 2 square lattice model, which is frustrated for antiferromagnetic interactions. Specifically, a high-temperature series expansion analysis of the temperature-dependent specific-heat and magnetization data yields $J$ 2 /$J$ 1 = 0.30 (=0.23) for YbBi 2 IO 4 (YbBi 2 ClO 4 ). In conclusion, this simplified analysis suggests strong frustration that should promote significant quantum fluctuations in these compounds, and thus motivates future work on the static and dynamic magnetic properties of these materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Absence of a bulk signature of a charge density wave in hard x-ray measurements of UTe 2

The long-sought pair density wave (PDW) is an exotic phase of matter in which charge density wave (CDW) order is intertwined with the amplitude or phase of coexisting, superconducting order. Originally predicted to exist in copper-oxides, circumstantial evidence for PDW order now exists in a variety of materials. Recently, scanning tunneling microscopy (STM) studies have reported evidence for a three-component CDW at the surface of the heavy-fermion superconductor, UTe 2 , persisting below its superconducting transition temperature [Aishwarya et al., Nature 618, 928933 (2023); Gu et al., Nature 618, 921927 (2023); LaFleur et al., Nat. Commun. 15, 4456 (2024)]. Here, in this study, we use hard x-ray diffraction measurements on crystals of UTe 2 at 𝑇 = 1.9 and 12 K to search for a bulk signature of this CDW. Using STM measurements as a constraint, we calculate the expected locations of CDW superlattice peaks and sweep a large volume of reciprocal space in search of a signature. We fail to find any evidence for a CDW near any of the expected superlattice positions in many Brillouin zones. We estimate an upper bound on the CDW lattice distortion of 𝑢 max ≲ 4×10 −3 Å. Our results suggest that the CDW observed in STM is either purely electronic, somehow lacking a signature in the lattice, or is restricted to the material surface.

Kengle, Caitlin S.↗
DOE OSTI2020

Hadronic interaction model sibyll 2.3d and extensive air showers

We present a new version of the hadron interaction event generator SIBYLL. While the core ideas of the model have been preserved, the new version handles the production of baryon pairs and leading particles in a new way. In addition, production of charmed hadrons is included. Updates to the model are informed by high-precision measurements of the total and inelastic cross sections with the forward detectors at the LHC that constrain the extrapolation to ultrahigh energy. Minimum-bias measurements of particle spectra and multiplicities support the tuning of fragmentation parameters. This paper demonstrates the impact of these changes on air-shower observables such as X max and N μ , drawing comparisons with other contemporary cosmic-ray interaction models.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Modifying PyUltraLight to model scalar dark matter with self-interactions

Here we introduce a modification of the pysiultralight code that models the dynamical evolution of ultralight axionlike scalar dark matter fields. Our modified code, pysiultralight, adds a quartic, self-interaction term to reflect the one which arises naturally in axionlike particle models. Using a particle mass of 10−22 eV/c 2 , we show that pysiultralight produces spatially oscillating solitons, exploding solitons, and collapsing solitons which prior analytic work shows will occur with attractive self-interactions. Using our code we calculate the oscillation frequency as a function of soliton mass and equilibrium radius in the presence of attractive self-interactions. We show that when the soliton mass is below the critical mass ($M_c$ = $\frac{\sqrt{3}}{2} M_{max}$) described by Chavanis and the initial radius is within a specific range, solitons are unstable and explode. We test the maximum mass criteria described by Chavanis and Chavanis and Delfini for a soliton to collapse when attractive self-interactions are included. We also analyze both binary soliton collisions and a soliton rotating around a central mass with attractive and repulsive self-interactions. We find that when attractive self-interactions are included, the density profiles get distorted after a binary collision. We also find that a soliton is less susceptible to tidal stripping when attractive self-interactions are included. We find that the opposite is true for repulsive self-interactions in that solitons would be more easily tidally stripped. Including self-interactions might therefore influence the survival timescales of infalling solitons.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Perturbation theory models for LSST-era galaxy clustering: Tests with subpercent mock catalog measurements in Fourier and configuration space

Here, we analyze the clustering of galaxies using the z = 1.006 snapshot of the CosmoDC2 simulation, a high-fidelity synthetic galaxy catalog designed to validate analysis methods for the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). We present sub-percent measurements of the galaxy auto-correlation and galaxy-dark matter cross correlations in Fourier space and configuration space for a magnitude-limited galaxy sample. At these levels of precision, the statistical errors of the measurement are comparable to the systematic effects present in the simulation and measurement procedure; nevertheless, using a hybrid-PT model, we are able to model non-linear galaxy bias with 0.5% precision up to scales of k max = 0.5 h/Mpc and r min = 4 Mpc/h. While the linear bias parameter is measured with 0.01% precision, other bias parameters are determined with considerably weaker constraints and sometimes bimodal posterior distributions. We compare our fiducial model with lower dimensional models where the higher-order bias parameters are fixed at their co-evolution values and find that leaving these parameters free provides significant improvements in our ability to model small scale information. We also compare bias parameters for galaxy samples defined using different magnitude bands and find agreement between samples containing equal numbers of galaxies. Finally, we compare bias parameters between Fourier space and configuration space and find moderate to significant tension between the two approaches. Although our model is often unable to fit the CosmoDC2 galaxy samples within the 0.1% precision of our measurements, our results suggest that the hybrid-PT model used in this analysis is capable of modeling non-linear galaxy bias within the percent level precision needed for upcoming galaxy surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Binary neutron star mergers as a probe of quark-hadron crossover equations of state

It is anticipated that the gravitational radiation detected in future gravitational wave (GW) detectors from binary neutron star (NS) mergers can probe the high-density equation of state (EOS). In this work, we perform the first simulations of binary NS mergers which adopt various parametrizations of the quark-hadron crossover (QHC) EOS. These are constructed from combinations of a hadronic EOS ($n_b$<2 $n_0$) and a quark-matter EOS ($n_b$>5 $n_0$), where $n_b$ and $n_0$ are the baryon number density and the nuclear saturation density, respectively. At the crossover densities (2 $n_0$ < $n_b$ < 5 $n_0$) the QHC EOSs continuously soften, while remaining stiffer than hadronic and first-order phase transition EOSs, achieving the stiffness of strongly correlated quark matter. This enhanced stiffness leads to significantly longer lifetimes of the postmerger NS than that for a pure hadronic EOS. We find a dual nature of these EOSs such that their maximum chirp GW frequencies $f_{\text{max}}$ fall into the category of a soft EOS while the dominant peak frequencies ($f_{\text{peak}}$) of the postmerger stage fall in between that of a soft and stiff hadronic EOS. An observation of this kind of dual nature in the characteristic GW frequencies will provide crucial evidence for the existence of strongly interacting quark matter at the crossover densities for QCD.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Simple Hamiltonian for quantum simulation of strongly coupled $(2+1)D$ SU(2) lattice gauge theory on a honeycomb lattice

Here, we find a simple spin Hamiltonian to describe physical states of $(2+1)$-dimensional SU(2) lattice gauge theory on a honeycomb lattice with a truncation of the electric field representation at $j_{max}=\frac{1}{2}$. The simple spin Hamiltonian contains only local products of Pauli matrices, even though Gauss’s law has been completely integrated out.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Quantum computing for energy correlators

In recent years, energy correlators have emerged as powerful observables for probing the fragmentation dynamics of high-energy collisions. We introduce the first numerical strategy for calculating energy correlators using the Hamiltonian lattice approach, providing access to the intriguing nonperturbative dynamics of these observables. Furthermore, motivated by rapid advances in quantum computing hardware and algorithms, we propose a quantum algorithm for calculating energy correlators in quantum field theories. This algorithm includes ground state preparation, the application of source, sink, energy flux and real-time evolution operators, and the Hadamard test. We validate our approach by applying it to the SU(2) pure gauge theory in 2 + 1 dimensions on 3 × 3 and 5 × 5 honeycomb lattices with 𝑗 max = $\frac{1}{2}$ at various couplings, utilizing both classical methods and the quantum algorithm, the latter tested using the IBM emulator for specific configurations. The results are consistent with the expected behavior of the strong coupling regime and motivate a more comprehensive study to probe the confinement dynamics across the weak and strong coupling regimes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Emergent hydrodynamic mode on SU(2) plaquette chains and quantum simulation

We search for emergent hydrodynamic modes in real-time Hamiltonian dynamics of 2+1-dimensional SU(2) lattice gauge theory on a quasi-one-dimensional plaquette chain, by numerically computing symmetric correlation functions of energy densities on lattice sizes of about 20 with the local Hilbert space truncated at 𝑗 max = $\frac{1}{2}$. Because of the Umklapp processes, we only find a mode for energy diffusion. The symmetric correlator exhibits transport peak near zero frequency with a width approximately proportional to momentum squared at small momentum, when the system is fully quantum ergodic, as indicated by the eigenenergy level statistics. This transport peak leads to a power-law 𝑡 −$\frac{1}{2}$ decay of the symmetric correlator at late time, also known as the long-time tail, as well as diffusionlike spreading in position space. We also introduce a quantum algorithm for computing the symmetric correlator on a quantum computer and find it gives results consistent with exact diagonalization when tested on the IBM emulator. Finally we discuss the future prospect of searching for the sound modes.

Hamiltonian systems↗
DOE OSTI2023

Generic two-phase coexistence in a type-2 Schloegl model for autocatalysis on a square lattice: Analysis via heterogeneous master equations

Schloegl's second model (also known as the quadratic contact process) on a square lattice involves spontaneous annihilation of particles at lattice sites at rate p, and their autocatalytic creation at unoccupied sites with n ≥ 2 occupied neighbors at rate k n . Here, kinetic Monte Carlo (KMC) simulation reveals that these models exhibit a nonequilibrium discontinuous phase transition with generic two-phase coexistence: the p value for equistability of coexisting populated and vacuum states, p eq (S), depends on the orientation or slope, S, of a planar interface separating those phases. The vacuum state displaces the populated state for p > eq (S), and the opposite applies for p < p eq (S) for 0 < S < ∞. The special “combinatorial” rate choice k n = n(n–1)/12 facilitates an appealing simplification of the exact master equations for the evolution of spatially heterogeneous states in the model, which aids analytic investigation of these equations via hierarchical truncation approximations. Truncation produces coupled sets of lattice differential equations which can describe orientation-dependent interface propagation and equistability. The pair approximation predicts that p eq (max) = p eq (S = 1) = 0.09645 and p eq (min) = p eq (S → ∞) = 0.08827, values deviating less than 15% from KMC predictions. In the pair approximation, a perfect vertical interface is stationary for all p < p eq (S = ∞) = 0.08907, a value exceeding p eq (S → ∞). One can regard an interface for large S → ∞ as a vertical interface decorated with isolated kinks. For p < p eq (S = ∞), the kink can move in either direction along this otherwise stationary interface depending upon p, but for p = p eq (min) the kink is also stationary.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Partially averaged Navier-Stokes closure modeling for variable-density turbulent flow

We report this work extends the framework of the partially averaged Navier-Stokes (PANS) equations to variable-density flow, i.e., multimaterial and/or compressible mixing problems with density variations and production of turbulence kinetic energy by both shear and buoyancy mechanisms. The proposed methodology is utilized to derive the PANS BHR-LEVM closure. This includes a priori testing to analyze and develop guidelines toward the efficient selection of the parameters controlling the physical resolution and, consequently, the range of resolved scales of PANS. Two archetypal test-cases involving transient turbulence, hydrodynamic instabilities, and coherent structures are used to illustrate the accuracy and potential of the method: the Taylor-Green vortex at Reynolds number Re = 3000, and the Rayleigh-Taylor flow at Atwood number 0.5 and (Re) max ≈ 500. These representative problems, for which turbulence is generated by shear and buoyancy processes, constitute the initial validation space of the new model, and their results are comprehensively discussed in two subsequent studies. The computations indicate that PANS can accurately predict the selected flow problems, resolving only a fraction of the scales of large-eddy simulation and direct numerical simulation strategies. The results also reiterate that the physical resolution of the PANS model must guarantee that the key instabilities and coherent structures of the flow are resolved. The remaining scales can be modeled through an adequate turbulence scale-dependent closure.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Destructive Error Interference in Product-Formula Lattice Simulation

Quantum computers can efficiently simulate the dynamics of quantum systems. Here, we study the cost of digitally simulating the dynamics of several physically relevant systems using the first-order product-formula algorithm. We show that the errors from different Trotterization steps in the algorithm can interfere destructively, yielding a much smaller error than previously estimated. In particular, we prove that the total error in simulating a nearest-neighbor interacting system of n sites for time t using the first-order product formula with r time slices is O ( n t / r + n t 3 / r 2 ) when n t 2 / r is less than a small constant. Given an error tolerance ϵ , the error bound yields an estimate of max { O ( n 2 t / ϵ ) , O ( n 2 t 3 / 2 / ϵ 1 / 2 ) } for the total gate count of the simulation. The estimate is tighter than previous bounds and matches the empirical performance observed in Childs et al. [ Proc. Natl. Acad. Sci. U.S.A. 115 , 9456 (2018) ]. We also provide numerical evidence for potential improvements and conjecture an even tighter estimate for the gate count.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Bending Magnet Synchrotron Radiation Imaging with Large Orbital Collection Angles

Synchrotron radiation (SR) from bending magnets, wigglers, and undulators is now extensively produced for users at storage ring based light sources, with unique properties in terms of average brightness and stability. Here, we present a profound study of bending magnet SR intensity distribution in the image plane of a focusing optical system. Measurements of this intensity distribution at the MAX-IV low emittance storage ring are compared to theoretical predictions, and found to be in excellent agreement. This work shows upon the possibility of performing high resolution emittance diagnostics with visible or near-visible SR on upcoming low-emittance storage ring based light sources. As a byproduct of our study, we derive a closed analytical expression for the intensity distribution from a zero-emittance beam, in the limiting case of wide orbital collection angles. This expression finally allows us to demonstrate the meeting between classical electrodynamics applied to SR emission and focusing, and the Landau and Lifshitz prediction of radiation intensity distribution nearby a caustic.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Optimal Twirling Depth for Classical Shadows in the Presence of Noise

The classical shadows protocol is an efficient strategy for estimating properties of an unknown state p using a small number of state copies and measurements. In its original form, it involves twirling the state with unitaries from some ensemble and measuring the twirled state in a fixed basis. It was recently shown that for computing local properties, optimal sample complexity (copies of the state required) is remarkably achieved for unitaries drawn from shallow depth circuits composed of local entangling gates, as opposed to purely local (zero depth) or global twirling (infinite depth) ensembles. Here, we consider the sample complexity as a function of the depth of the circuit, in the presence of noise. We find that this noise has important implications for determining the optimal twirling ensemble. Under fairly general conditions, we (i) show that any single-site noise can be accounted for using a depolarizing noise channel with an appropriate damping parameter f, (ii) compute thresholds f th at which optimal twirling reduces to local twirling for Pauli operators, (iii) nth order Renyi entropies (n ≥2), and (iv) provide a meaningful upper bound t max on the optimal circuit depth for any finite noise strength f, which applies to observables and entanglement entropy measurements. In conclusion, these thresholds strongly constrain the search for optimal strategies to implement shadow tomography and are easily tailored to the experimental system at hand.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

BioXTAS RAW 2 : new developments for a free open-source program for small-angle scattering data reduction and analysis

BioXTAS RAW is a free open-source program for reduction, analysis and modelling of biological small-angle scattering data. Here, the new developments in RAW version 2 are described. These include improved data reduction using pyFAI ; updated automated Guinier fitting and D max finding algorithms; automated series ( e.g. size-exclusion chromatography coupled small-angle X-ray scattering or SEC-SAXS) buffer- and sample-region finding algorithms; linear and integral baseline correction for series; deconvolution of series data using regularized alternating least squares ( REGALS ); creation of electron-density reconstructions using electron density via solution scattering ( DENSS ); a comparison window showing residuals, ratios and statistical comparisons between profiles; and generation of PDF reports with summary plots and tables for all analysis. Furthermore, there is now a RAW API, which can be used without the graphical user interface (GUI), providing full access to all of the functionality found in the GUI. In addition to these new capabilities, RAW has undergone significant technical updates, such as adding Python 3 compatibility, and has entirely new documentation available both online and in the program.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

DS-GL: Advancing Graph Learning via Harnessing the Power of Nature within Dynamic Systems

With the rapid digitization of the world, an increasing number of real-world applications are turning to nonEuclidean data, modeled as graphs. Due to their intrinsic high complexity and irregularity, learning from graph data demands tremendous computational power. Recently, CMOS-compatible Ising machines, i.e., dynamic systems composed of CMOS components, have emerged as a new approach that harnesses the inherent power of natural annealing within dynamic systems to efficiently resolve binary optimization problems and have been adopted for traditional graph computation, such as max-cut. However, when performing complex Graph Learning (GL) tasks, Ising machines face significant hurdles: (i) they are inherently binary and thus ill-suited for real-valued problems; (ii) their expensive all-to-all coupling network that guarantees effective natural annealing poses daunting scalability concerns. To address these challenges, this paper proposes a nature-powered graph learning framework dubbed DS-GL, which is the first effort to transform the process of solving graph learning problems into the natural annealing process within a parameterized dynamic system embodied as a CMOS chip. To tackle the two major hurdles, DS-GL first augments the Ising machine architecture to modify the self-reaction term of its Hamiltonian function from linear to quadratic, effectively serving as an energy regulator. This adjustment maintains the system’s original physical interpretation while enabling it to process continuous, real-valued data. Second, to address the scaling issue, DS-GL further upgrades the real-valued dense Ising machine by decomposing it into a mesh-based multi-PE dynamic system that supports efficient distributed spatial-temporal co-annealing across different PEs through sparse interconnects. By exploiting the inherent sparsity and component structures in real-world graphs, DS-GL is able to map complex graph learning tasks onto the scalable dynamic system while maintaining high accuracy. Evaluations with three diverse GL applications across six real-world datasets, including traffic flow and COVID-19 prediction, show that DS-GL can deliver from 102× to 106× speedups and 500× energy reduction over Graph Neural Networks on GPUs, with 5% - 20% accuracy enhancement.

Song, Ruibing↗