Search NASASearch

SEARCH · Search NASA

Results for “Analytic linearization”

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 37 records · Page 2

Coupling to rotational manifolds to improve gas-phase pump–probe spectroscopic models

The physical picture of gas-phase optical transitions is normally presented as an isolated two-level system balanced by upward and downward processes. Isolated models assume a phenomenological treatment of collisional dephasing but do not strictly account for collisional population exchange with the rotational baths. While this assumption is valid under low-intensity conditions, where excitation is rate-limiting, isolated models can deviate from Beer’s Law at sufficient pressures and monochromatic intensities when both collisional broadening and power broadening are comparable to (or greater than) lifetime broadening, which are not uncommon conditions for cavity enhanced spectroscopies in the mid-IR spectral range. Although this problem has been addressed by rate-equation models for linear absorption measurements, a general treatment for multi-level quantum mechanical models suitable for non-linear absorption measurements (two-photon/two-color/pump–probe) is lacking. Isolated models require physical parameter inputs that disagree with expected values by at least an order of magnitude. These non-physical models undermine the ability to predict non-linear signal strengths under untested conditions and thereby limit the potential to optimize the sensitivity of non-linear spectroscopies and to expand their analytical applications (e.g., new analytes and/or buffer gases, changes in cavity free-spectral-range, changes in intracavity powers or wavelengths, and accurate investigation of physical phenomena). In this study, we derive bath-coupled models for gaseous pump–probe spectroscopy by application of the quantum Lindblad equation and detailed balance. Bath-coupled models are shown to fit data consistently across variations in intensity and agree with all physically expected values.

Cavity ring-down spectroscopy

Detection of the large-scale tidal field with galaxy multiplet alignment in the DESI Y1 spectroscopic survey

We explore correlations between the orientations of small galaxy groups, or ‘multiplets’, and the large-scale gravitational tidal field. Using data from the Dark Energy Spectroscopic Instrument (DESI) Y1 survey, we detect the intrinsic alignment (IA) of multiplets to the galaxy-traced matter field out to separations of $100\,h^{-1}$ Mpc. Unlike traditional IA measurements of individual galaxies, this estimator is not limited by imaging of galaxy shapes and allows for direct IA detection beyond redshift $z=1$. Multiplet alignment is a form of higher order clustering, for which the scale-dependence traces the underlying tidal field and amplitude is a result of small-scale ($\lt 1h^{-1}$ Mpc) dynamics. Within samples of bright galaxies, luminous red galaxies (LRG) and emission-line galaxies, we find similar scale-dependence regardless of intrinsic luminosity or colour. This is promising for measuring tidal alignment in galaxy samples that typically display no IA. DESI’s LRG mock galaxy catalogues created from the A BACUS S UMMIT N -body simulations produce a similar alignment signal, though with a 33 per cent lower amplitude at all scales. An analytic model using a non-linear power spectrum (NLA) only matches the signal down to 20 $h^{-1}$ Mpc. Our detection demonstrates that galaxy clustering in the non-linear regime of structure formation preserves an interpretable memory of the large-scale tidal field. Multiplet alignment complements traditional two-point measurements by retaining directional information imprinted by tidal forces, and contains additional line-of-sight information compared to weak lensing. This is a more effective estimator than the alignment of individual galaxies in dense, blue, or faint galaxy samples.

79 ASTRONOMY AND ASTROPHYSICS

Current-based metrology with two-terminal mesoscopic conductors

The traditional approach to quantum parameter estimation focuses on the quantum state, deriving fundamental bounds on precision through the quantum Fisher information. In most experimental settings, however, performing arbitrary quantum measurements is highly unfeasible. In open quantum systems, an alternative approach to metrology involves the measurement of stochastic currents flowing from the system to its environment. However, the present understanding of current-based metrology is mostly limited to Markovian master equations. Considering a parameter estimation problem in a two-terminal mesoscopic conductor, we identify the key elements that determine estimation precision within the Landauer-Büttiker formalism. Crucially, this approach allows us to address arbitrary coupling and temperature regimes. Furthermore, we obtain analytical results for the precision in linear-response and zero-temperature regimes. For the specific parameter estimation task that we consider, we demonstrate that the boxcar transmission function is optimal for current-based metrology in all parameter regimes.

Landauer formula

Radiation by a laser-driven flying-focus electron wave packet

An exact solution of the Dirac equation in the presence of an arbitrary electromagnetic plane wave is found, which corresponds to a focused electron wave packet, with the focus of the wave packet moving at the speed of light in the opposite direction of the average momentum of the electron wave packet (unless the plane wave is so intense to reflect the electron). The photon spectrum emitted by such an electron wave packet in the presence of a linearly polarized plane wave is studied both analytically and numerically. The spectrum is also compared with the one emitted by a single-momentum, plane-wave electron in the case of the electron being initially counterpropagating (on average for the flying-focus case) with the plane wave and within the locally constant field approximation. It is found that if the electron flying-focus wave packet is focused beyond a Compton wavelength, the angular distribution of the emitted radiation along the magnetic field of the electromagnetic plane wave is broader than for an electron with definite momentum. Corresponding the maximum value of the photon yield on the transverse plane is smaller in the flying-focus electron case. This could represent an experimental signature of a laser-driven flying-focus electron wave packet.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS

Integrable symplectic maps with a polygon tessellation

Identifying integrable dynamics remains a formidable challenge, and despite centuries of research, only a handful of examples are known to date. In this article, we explore a distinct form of area-preserving (symplectic) mappings derived from the stroboscopic Poincaré cross section of a kicked rotator—an oscillator subjected to an external force periodically switched on in short pulses. The significance of this class of problems extends to various applications in physics and mathematics, including particle accelerators, crystallography, and studies of chaos. Notably, Suris's theorem constrains the integrability within this category of mappings, outlining potential scenarios with analytic invariants of motion. In this paper, we challenge the assumption of the analyticity of the invariant by exploring piecewise linear transformations on a torus ( T 2 ) and associated systems on the plane ( R 2 ), incorporating arithmetic quasiperiodicity and discontinuities. Introducing a new automated technique, we discovered previously unknown scenarios featuring polygonal invariants that form perfect tessellations and, moreover, fibrations of the plane or torus. This work reveals a novel category of planar tilings characterized by discrete symmetries that emerge from the invertibility of transformations and are intrinsically linked to the presence of integrability. Our algorithm relies on the analysis of the Poincaré rotation number and its piecewise monotonic nature for integrable cases, contrasting with the noisy behavior in the case of chaos, thereby allowing for clear separation. Some of the newly discovered systems exhibit the peculiar behavior of “integrable diffusion,” characterized by infinite and quasirandom hopping between tiles while being confined to a set of invariant segments. Finally, through the implementation of a smoothening procedure, all mappings can be generalized to quasi-integrable scenarios with suppressed volume occupied by chaotic trajectories, thereby opening doors to potential practical applications. Published by the American Physical Society 2024

43 PARTICLE ACCELERATORS

The 3D clustering of Lyman Alpha Emitters measured with DESI

We present a clustering analysis of Lyman-$α$ emitters (LAEs) using spectroscopic observations from the Dark Energy Spectroscopic Instrument (DESI) of candidates selected from the Blanco/DECam Intermediate-Band Imaging Survey (IBIS). We measure the two-point correlation function and the power spectrum, including cross-correlations with DESI quasars. Using both analytical and halo occupation distribution (HOD) simulation-based modeling, we find a linear bias of $b \sim 2.31$--$2.62$ for LAEs over the redshift range $2.26 < z < 3.41$. The analytical modeling also provides constraints on the strength of radiative transfer effects, while the HOD analysis characterizes the LAE-halo connection across multiple models. Finally, we quantify the magnitude of non-perturbative clustering effects such as Fingers of God in the LAE population, providing essential input for the accurate modeling of LAE-based cosmological analyses in forthcoming high-redshift surveys such as DESI-II.

Ebina, H. [UC, Berkeley; LBL, Berkeley] (ORCID:000

A Tutorial on Bayesian analysis of linear shock compression data

Gas gun and other shock compression experiments often produce shock wave velocity measurements that are linearly associated with particle velocity. Traditionally, this empirical relationship is quantified with a single Hugoniot curve that is estimated using least squares regression. However, for downstream modeling and simulation tasks, it is often more useful to have multiple Hugoniot curves in the pressure–volume plane that are consistent with the data. We employ Bayesian uncertainty quantification methods as a framework for propagating measurement uncertainty through to model parameters and predictions. Specifically, this Tutorial shows how to sample multiple Hugoniot curves in the pressure–volume plane that are consistent with the shock wave-particle velocity measurements in a two-step Bayesian approach. First, we obtain an analytical expression for the posterior distribution of the linear model parameters using Bayesian linear regression. Second, we propagate samples from the posterior distribution through the Rankine–Hugoniot equations to yield Hugoniot curves in the pressure–volume plane. The procedure is demonstrated with publicly available data on argon, copper, and nickel, and compared against bootstrapping and linear regression. The Bayesian procedure is shown to be interpretable, computationally inexpensive, and less sensitive than an alternative bootstrapping approach to the removal of the point in the copper dataset that has the largest particle velocity. As a Tutorial on Bayesian methodology for the shock compression community, we provide several derivations and explanations that make this paper self-contained, and make all code and data available at github.com/llnl/BALSCD.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC

High-order limiting methods using maximum principle bounds derived from the Boltzmann equation I: Euler equations

The use of limiting methods for high-order numerical approximations of hyperbolic conservation laws generally requires defining an admissible region/bounds for the solution. In this work, we present a novel approach for computing solution bounds and limiting for the Euler equations through the kinetic representation provided by the Boltzmann equation, which allows for extending limiters designed for linear advection directly to the Euler equations. Given an arbitrary set of solution values to compute bounds over (e.g., numerical stencil) and a desired linear advection limiter, the proposed approach yields an analytic expression for the admissible region of particle distribution function values, which may be numerically integrated to yield a set of bounds for the density, momentum, and total energy. Further, these solution bounds are shown to preserve positivity of density/pressure/internal energy and, when paired with a limiting technique, can robustly resolve strong discontinuities while recovering high-order accuracy in smooth regions without any ad hoc corrections (e.g., relaxing the bounds). This approach is demonstrated in the context of an explicit unstructured high-order discontinuous Galerkin/flux reconstruction scheme for a variety of difficult problems in gas dynamics, including cases with extreme shocks and shock-vortex interactions. Furthermore, this work presents a foundation for limiting techniques for more complex macroscopic governing equations that can be derived from an underlying kinetic representation for which admissible solution bounds are not well-understood.

42 ENGINEERING

Laser polarization effects on crossed-beam energy transfer in inertial confinement fusion

Crossed-beam energy transfer (CBET) is a consequential process in inertial confinement fusion (ICF) experiments. CBET depends on the polarization state of the interacting lasers and can, in turn, alter their polarization. Here, we derive analytical expressions for two-beam CBET with linearly and circularly polarized light, and present simulations of CBET that include polarization effects for ICF conditions relevant to experiments on the National Ignition Facility. The results show that CBET mixes the polarization of the beams as they propagate through the target and leads to azimuthal power imbalances between beams within a given cone (i.e., at the same polar angle in the target chamber). These variations are shown to correlate with variations in measured backscatter between beams of a same quadruplet. Using circularly polarized, instead of linearly polarized, light produces similar overall levels of CBET within a cone of beams but tends to reduce the azimuthal variations. This is expected to slightly improve irradiation symmetry and reduce the risk of backscatter from stimulated Brillouin scattering. Circular polarization may therefore be beneficial for the next generation of laser fusion drivers.

Physics - Plasma physics

Machine Learning–Augmented Laser-Induced Breakdown Spectroscopy for Spectral Discrimination of Iron Oxalates

Enhanced characterization and phase identification of post-PUREX Pu Oxalates (PuOXA) are pivotal for nonproliferation and pre-detonation nuclear forensics. Despite significant advances in the characterization of PuO 2 samples, little is known about the impact of both the chemical structure and oxidation states of PuOXA (i.e., Pu(III) and Pu(IV)) have on optical emission signatures. Here, we demonstrate the analytical capabilities of laser-induced breakdown spectroscopy (LIBS) applied to Fe(II) and Fe(III) oxalate samples as surrogates for PuOXA, highlighting the discriminating features in the LIBS emission spectra arising from differences in the oxidation states within mixed FeOXA samples. We report the enhancement of spectral feature selection using Principal Component Analysis (PCA), which enables the analytical superiority of machine learning algorithms such as Linear Discriminant Analysis (LDA), Quadratic Discriminant Analysis (QDA), Partial Least Squares Regression (PLSR), Support Vector Regression (SVR), and Random Forest Regression (RFR) over conventional univariate techniques for phase discrimination and chemometric analysis. Cluster analysis revealed how both matrix effects and laser ablation influence cluster separability by introducing spectral artifacts that misdirect the maximization of variance. PCA-selected emission lines were used in the regression models, demonstrating that both univariate and multivariate linear regression models (i.e., PLSR and SVR) can achieve acceptable performance, with machine learning models outperforming conventional calibration regressions. Furthermore, the application of non-linearly activated PCA-selected emission lines illustrates how simplifying the data while retaining captured variance enables the use of less complex and more computationally efficient models. Furthermore, this is particularly evident in the underperformance of RFR, which suffers from increased computational costs and overfitting owing to its high complexity.

Oxalates

Taylor-Expansion-Based Robust Power Flow in Unbalanced Distribution Systems: A Hybrid Data-Aided Method

Traditional power flow methods often adopt certain assumptions designed for passive balanced distribution systems, thus lacking practicality for unbalanced operation. moreover, their computation accuracy and efficiency are heavily subject to unknown errors and bad data in measurements or prediction data of distributed energy resources (ders). to address these issues, this paper proposes a hybrid data-aided robust power flow algorithm in unbalanced distribution systems, which combines taylor series expansion knowledge with a data-driven regression technique. the proposed method initiates a linearization power flow model to derive an explicitly analytical solution by modified taylor expansion. to mitigate the approximation loss that surges due to the der integration and bad data, we further develop a data-aided robust support vector regression approach to estimate the errors efficiently. comparative analysis in the 13-bus and 123-bus ieee unbalanced feeders shows that the proposed hybrid algorithm achieves superior computational efficiency, with guaranteed accuracy and robustness against outliers.

data-driven

Modification of ion-temperature-gradient turbulence by impurities in stellarator plasmas

Recent nonlinear gyrokinetic simulations have shown that impurities can strongly modify the turbulent heat flux in stellarator plasmas. Here, the ion-temperature-gradient (ITG) dispersion relation in a plasma containing impurities is analytically solved in certain limits and an expression for the modification of the ITG growth rate by impurities is derived. The analytical expression is the sum of three terms corresponding to three different physical causes (impurity density gradient, impurity temperature gradient and dilution) of the change in the growth rate. The scalings predicted analytically for the modification of the growth rate are shown to be reproduced by linear gyrokinetic simulations. The conditions for reduction or increase of the ITG growth by impurities are also correctly predicted by the analytical solution to the dispersion relation. Finally, a remarkable correlation is found between the analytical expression for the modification of the growth rate and the modification of the turbulent heat flux obtained from nonlinear gyrokinetic simulations.

Calvo, Iván [Research Centre for Energy, Environme

Electron– K -phonon interaction in twisted bilayer graphene

We develop an analytic theory to describe the interaction between electrons and K phonons and study its influence on superconductivity in the bare bands of twisted bilayer graphene (TBG). We find that, due to symmetry and the two-center approximation, only one optical K phonon (160meV) of graphene is responsible for the intervalley electron-phonon interaction. Here, this phonon has recently been found in angular-resolved photoemission spectroscopy to be responsible for replicas of the TBG flat bands. By projecting the interaction to the TBG flat bands, we perform the full symmetry analysis of the phonon-mediated attractive interaction and pairing channels in the Chern basis, and show that several channels are guaranteed to have gapless order parameters. From the linearized gap equations, we find that the highest T c pairing induced by this phonon is a singlet gapped s-wave inter-Chern-band order parameter, followed closely by a gapless nematic d-wave intra-Chern-band order parameter. We justify these results analytically, using the topological heavy-fermion mapping of TBG which has allowed us to obtain an analytic form of a phonon-mediated attractive interaction and to analytically solve the linearized and T = 0 gap equations. For the intra-Chern-band channel, the nematic state with nodes is shown to be stabilized in the chiral flat-band limit. While the flat-band Coulomb interaction can be screened sufficiently enough—around the Van Hove singularities—to allow for electron-phonon based superconductivity, it is unlikely that this effect can be maintained in the lower density of states excitation bands around the correlated insulator states.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND

Microstructurally Strained Pyrochlore–Perovskite Biphasic Electrocatalysts for the Oxygen Evolution Reaction

Efficiency of water splitting for hydrogen production is often limited by the sluggish kinetics of multiple electronic transfers required in the heterogeneous oxygen evolution reaction (OER). Catalyst design for reducing the high OER overpotential remains a major scientific challenge. Lattice-strain engineering, a method for tuning the electronic structure and surface geometric configuration of active sites, may greatly affect the interaction between adsorbates and catalytic surfaces for high activity and stability. Here, in this study, we present the synthesis of biphasic oxides of YPrSrRuMnO x , which consists of distinct phases of Y 2 Ru 2 O 7 pyrochlore and (Pr 0.7 Sr 0.3 )MnO 3 perovskite, and the development of a suitable analytical approach to study the strain–catalytic property relationship. Linear sweep voltammetry results reveal that the biphasic oxide exhibits approximately 3.1 times greater mass activity and 2.4 times larger turnover frequency (TOF) than single-phase Y 2 Ru 2 O 7 in the 0.1 M HClO 4 electrolyte. The biphasic catalyst is also about 3 times more stable than the single-phase oxide under acidic conditions. X-ray photoelectron spectroscopy, nitrogen isotherm, and electrochemical surface area analyses indicate that the oxidation state, specific surface area, and electrochemical surface area do not cause enough difference in the observed enhancement of OER performance. We examined the effects of microstrain on electrocatalysis, originating from lattice mismatch between different phases, using three different structural models. Specifically, we compared the Williamson–Hall method, standard stress–strain analysis, and Rietveld refinement in analyzing the structure–property relationship. Strain mapping using geometric phase analysis (GPA) further revealed significant microstrain and lattice dislocations localized near phase boundaries in the biphasic oxide, in contrast to the uniform strain in single-phase materials. The results reveal that the increased microstrain correlates well with the improved OER performance, as the biphasic oxide catalyst exhibits 2–3 times greater microstrain than Y 2 Ru 2 O 7 pyrochlore.

electrocatalysts

Towards large-scale quantum optimization solvers with few qubits

Quantum computers hold the promise of more efficient combinatorial optimization solvers, which could be game-changing for a broad range of applications. However, a bottleneck for materializing such advantages is that, in order to challenge classical algorithms in practice, mainstream approaches require a number of qubits prohibitively large for near-term hardware. Here we introduce a variational solver for MaxCut problems over $m={{\mathcal{O}}}({n}^{k})$ binary variables using only n qubits, with tunable k > 1. The number of parameters and circuit depth display mild linear and sublinear scalings in m , respectively. Moreover, we analytically prove that the specific qubit-efficient encoding brings in a super-polynomial mitigation of barren plateaus as a built-in feature. Altogether, this leads to high quantum-solver performances. For instance, for m = 7000, numerical simulations produce solutions competitive in quality with state-of-the-art classical solvers. In turn, for m = 2000, experiments with n = 17 trapped-ion qubits feature MaxCut approximation ratios estimated to be beyond the hardness threshold 0.941. Our findings offer an interesting heuristics for quantum-inspired solvers as well as a promising route towards solving commercially-relevant problems on near-term quantum devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC

Linear complexity

We present factorization and solution phases for a new linear complexity direct solver designed for concurrent batch operations on fine-grained parallel architectures, for matrices amenable to hierarchical representation. We focus on the strong-admissibility-based $\mathscr{H}^{2}$ format, where strong recursive skeletonization factorization compresses remote interactions. We build upon previous implementations of $\mathscr{H}^{2}$ matrix construction for efficient factorization and solution algorithm design, which are illustrated graphically in stepwise detail. The algorithms are ‘blackbox’ in the sense that the only inputs are the matrix and right-hand side, without analytical or geometrical information about the origin of the system. We demonstrate linear complexity scaling in both time and memory on four representative families of dense matrices up to one million in size. Parallel scaling up to 16 threads is enabled by a multi-level matrix graph coloring and avoidance of dynamic memory allocations thanks to prefix-sum memory management. An experimental backward error analysis is included. We break down the timings of different phases, identify phases that are memory-bandwidth limited, and discuss alternatives for phases that may be sensitive to the trend to employ lower precisions for performance.

Boukaram, Wajih

Final Technical Report: Testing the robustness of CMB-S4 inflation measurements to foreground models

This report summarizes the research activities for DOE award SC0024462, which focused on enhancing the robustness of CMB-S4 inflation measurements by characterizing and mitigating Galactic foreground interference. The project grounded foreground models in the physical processes of the Interstellar Medium with MHD and semi-analytic filament models, and developed analysis pipelines, such as Needlet Internal Linear Combination (NILC) techniques, to isolate faint cosmological signals.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS

Analytic model for neutral penetration and plasma fueling

Neutral atoms recycled from wall interaction interact with confined plasma, thereby refuelling it, most strongly in the region closest to the wall. This occurs near the X-point in diverted configurations, or else near the wall itself in limited configurations. A progression of analytic models is developed for neutral density in the vicinity of a planar or linear source in an ionising domain. First-principles neutral transport simulations with DEGAS2 are used throughout to test the validity and limits of the model when using equivalent sources. The model is further generalised for strong plasma gradients or the inclusion of charge exchange. An important part of the problem of neutral fuelling from recycling is thereby isolated and solved with a closed-form analytic model. A key finding is that charge exchange with the confined plasma can be significantly simplified with a reasonable sacrifice of accuracy by treating it as a loss. The several assumptions inherent to the model (and the simulations with which it is compared) can be adapted according to the particular behaviour of neutrals in the divertor and the manner in which they cross the separatrix.

fusion plasma