Search NASA⌕ Search

SEARCH · Search NASA

Results for “fixed points”

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 181 records · Page 10

Two-Level Sketching Alternating Anderson Acceleration for Complex Physics Applications

We present a novel two-level sketching extension of the Alternating Anderson–Picard (AAP) method for accelerating fixed-point iterations in challenging single- and multiphysics simulations governed by discretized PDEs. Our approach combines a static, physics-based projection that reduces the least-squares (LS) problem to the most informative field (e.g., via Schur-complement insight) with a dynamic, algebraic sketching stage driven by a backward stability analysis under Lipschitz continuity. We introduce inexpensive estimators for stability thresholds and cache-aware randomized selection strategies to balance computational cost against memory access overhead. The resulting algorithm solves reduced LS systems in place, minimizes memory footprints, and seamlessly alternates between low-cost Picard updates and Anderson mixing. Implemented in Julia, our two-level sketching AAP achieves up to 50% time-to-solution reductions compared to standard Anderson acceleration—without degrading convergence rates—on benchmark problems including Stokes, 𝑝-Laplacian, bidomain, and Navier–Stokes formulations at varying problem sizes. These results demonstrate the method’s robustness, scalability, and potential for integration into high-performance scientific computing frameworks. Our implementation is available open source in the AAP.jl library.

Barnafi, Nicolas [University of Chile, Santiago]↗

On two nucleons near unitarity with perturbative pions

We explore the impact of perturbative pions on the Unitarity Expansion in the two-nucleon S-waves of Chiral Effective Field Theory at next-to-next-to leading order (N 2 LO). Pion exchange explicitly breaks the nontrivial fixed point’s universality, i.e. invariance of S waves under both conformal and Wigner’s combined SU(4) spin–isospin transformations. On the other hand, Unitarity explicitly breaks chiral symmetry. The two seem incompatible in their respective exact-symmetry limits. χEFT with Perturbative Pions in the Unitarity Expansion resolves the apparent conflict in the Unitarity Window (phase shifts 45° ≲ δ(k) ≲ 135°), i.e. around momenta most relevant for low-energy nuclear systems. Its only LO scale is the scattering momentum; NLO adds only scattering length, effective range and non-iterated one-pion exchange (OPE); and N 2 LO only once-iterated OPE. Agreement in the 1 S 0 channel is very good.

Chiral symmetry↗

pnnl/mcl-runtime

The Minos Computing Library (MCL) is a task-based programming language and runtime for extremely heterogeneous systems. MCL facilitates writing program for heterogeneous devices and porting applications across different systems and devices. MCL supports asynchronous execution of computing tasks on all available heterogeneous devices, including GPUs, FGPAs, fixed-point accelerators, and AI accelerators. MCL provides a high-level programming interface and automatically and autonomously performs resource management, load balancing, and locality-aware scheduling.

Central, PNNL Developer↗

ORNL/Restricted-Alternating-Anderson-Picard

Implementation in the programming language Julia of the Restricted Alternating Anderson Picard (AAP) acceleration for (non)-linear fixed point iterations.

Barnafi, NicolásAlejandro [Pontifica Universidad C↗

Vertical nanoscale strain-induced electronic localization in epitaxial La2/3Sr1/3MnO3 films with ZrO2 nanopillar inclusions

Abstract Unusual electrical transport properties associated with weak or strong localization are sometimes found in disordered electronic materials. Here, we report experimental observation of a crossover of electronic behavior from weak localization to enhanced weak localization due to the spatial influence of disorder induced by ZrO 2 nanopillars in (La 2/3 Sr 1/3 MnO 3 ) 1− x :(ZrO 2 ) x ( x = 0, 0.2, and 0.3) nanocomposite films. The spatial strain regions, identified by scanning transmission electron microscopy and high-resolution x-ray diffraction, induce a coexistence of two-dimentional (2D) and three-dimentional (3D) localization and switches to typical 2D localization with increasing density of ZrO 2 pillars due to length scale confinement, which interestingly accords with enhancing vertically interfacial strain. Based on the excellent agreement of our experimental results with one-parameter scaling theory of localization, the enhanced weak localization exists in metal range close to the fixed point. These films provide a tunable experimental model for studying localization in particular the transition regime by appropriate choice of the second epitaxial phase. Graphical Abstract

36 MATERIALS SCIENCE↗

Influence of Lake Ice Biases in Reanalysis Data on Downscaled Climate Simulations over the Great Lakes Region

This data package contains observation-based and model-simulated datasets (all provided in NetCDF format) for evaluating how wintertime lake-ice representation affects regional weather and climate over the Laurentian Great Lakes (freshwater lake ecosystem) during the high–ice-cover winter of 2009. The observational component includes: (1) Stage IV gridded precipitation at 4 km, hourly resolution for January–February 2009 over the Great Lakes region (radar–gauge multisensor precipitation analyses); (2) Great Lakes Surface Environmental Analysis (GLSEA) satellite-derived lake-ice coverage at 1.3 km, daily resolution for the 2009 winter months, providing ice coverage over Lakes Superior, Michigan, Huron, Erie, and Ontario; and (3) in situ measurements at the Standard Rock site on Lake Superior from the Great Lakes Evaporation Network (GLEN) at hourly resolution, including near-surface atmospheric variables and sensible and latent heat fluxes (air–lake exchange) at a fixed point location. The modeling component provides corresponding fields from two simulations, both archived at 4 km, hourly resolution: a standalone Weather Research Forecasting model (WRF) run driven by the European Centre for Medium-Range Weather Forecasts (ECMWF) Reanalysis 5 (ERA5), and a two-way coupled model using WRF and the Finite Volume Community Ocean Model (WRF-FVCOM, a 3-D hydrodynamic lake model). These outputs include variables relevant to air–lake interaction and lake-effect processes (e.g., near-surface temperature, humidity, wind, precipitation, and surface turbulent fluxes), enabling direct comparison with the observational datasets. Users can analyze and visualize these NetCDF files with common tools such as Python (e.g., xarray, netCDF4, numpy, pandas), NCO/CDO, Panoply, or ncview; NetCDF variables can also be converted to other formats (e.g., CSV, GeoTIFF) using these utilities.

EARTH SCIENCE > ATMOSPHERE > ATMOSPHERIC TEMPERATU↗

Relocation of the Seismicity of the Caucasus Region

Our objective is to improve the view of the seismicity in the Caucasus region using instrumental data between 1951 and 2019. To create a comprehensive catalog, we combine the bulletins of local agencies and the International Seismological Centre, and use an advanced single-event location algorithm, iLoc, to obtain better locations. We show that relocations with iLoc, using travel-time predictions from the 3D upper mantle velocity model, Regional Seismic Travel Time, improve the locations. Then, using the iLoc results as initial locations and the ground-truth events identified in the iLoc results as fix points, we apply Bayesloc, a multiple-event location algorithm, to simultaneously relocate the entire seismicity of the Caucasus region. Here we demonstrate that the simultaneous relocation of the seismicity with Bayesloc clarifies the location and geometry of major active structures accommodating ongoing convergence between the Arabian and Eurasian continents between the Black and Caspian Seas. Among our major findings is the confirmation of widespread seismicity in the mantle beneath the northern flank of the Greater Caucasus and central Caspian, resulting from north-dipping subduction of the Kura and South Caspian basins and the identification of a discrete band of crustal seismicity beneath the southern flank of the Greater Caucasus.

58 GEOSCIENCES↗

A model of persistent breaking of continuous symmetry

We consider a UV-complete field-theoretic model in general dimensions, including d=2+1, which consists of two copies of the long-range vector models, with O(m) and O(N-m) global symmetry groups, perturbed by double-trace operators. Using conformal perturbation theory we find weakly-coupled IR fixed points for N\geq 6 N ≥ 6 that reveal a spontaneous breaking of global symmetry. Namely, at finite temperature the lower rank group is broken, with the pattern persisting at all temperatures due to scale-invariance. We provide evidence that the models in question are unitary and invariant under full conformal symmetry. Furthermore, we show that this model exhibits a continuous family of weakly interacting field theories at finite N.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Stress Tensor flows, birefringence in non-linear electrodynamics and supersymmetry

We identify the unique stress tensor deformation which preserves zero-birefringence conditions in non-linear electrodynamics, which is a 4d 4 d version of the T\overline{T} T T ¯ operator. We study the flows driven by this operator in the three Lagrangian theories without birefringence - Born-Infeld, Plebanski, and reverse Born-Infeld - all of which admit ModMax-like generalizations using a root- T\overline{T} T T ¯ -like flow that we analyse in our paper. We demonstrate one way of making this root- T\overline{T} T T ¯ -like flow manifestly supersymmetric by writing the deforming operator in \mathcal{N} = 1 𝒩 = 1 superspace and exhibit two examples of superspace flows. We present scalar analogues in d = 2 d = 2 with similar properties as these theories of electrodynamics in d = 4 d = 4 . Surprisingly, the Plebanski-type theories are fixed points of the classical T\overline{T} T T ¯ -like flows, while the Born-Infeld-type examples satisfy new flow equations driven by relevant operators constructed from the stress tensor. Finally, we prove that any theory obtained from a classical stress-tensor-squared deformation of a conformal field theory gives rise to a related “subtracted” theory for which the stress-tensor-squared operator is a constant.

Physics↗

Prethermalization in coupled one-dimensional quantum gases

We consider the problem of the development of steady states in one-dimensional Bose gas tubes that are weakly coupled to one another through a density-density interaction. We analyze this development through a Boltzmann collision integral approach. We argue that when the leading order of the collision integral, where single particle-hole excitations are created in individual gases, is dominant, the state of the gas evolves first to a non-thermal fixed point, i.e. a prethermalization plateau. This order is dominant when a pair of tubes are inequivalent with, say, different temperatures or different effective interaction parameters, \gamma γ . When both tubes are in the strongly interacting regime we additionally characterize this non-thermal prethermalization plateau by constructing the quasi-conserved quantities that control the existence of this plateau as well as the associated generalized Gibbs ensemble.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

(SPT-)LSM theorems from projective non-invertible symmetries

Projective symmetries are ubiquitous in quantum lattice models and can be leveraged to constrain their phase diagram and entanglement structure. In this paper, we investigate the consequences of projective algebras formed by non-invertible symmetries and lattice translations in a generalized 1+1 1 + 1 D quantum XY model based on group-valued qudits. This model is specified by a finite group G G and enjoys a projective \mathsf{Rep}(G)× Z(G) 𝖱 𝖾 𝗉 ( G ) × Z ( G ) and translation symmetry, where symmetry operators obey a projective algebra in the presence of symmetry defects. For invertible symmetries, such projective algebras imply Lieb-Schultz-Mattis (LSM) anomalies. However, this is not generally true for non-invertible symmetries, and we derive a condition on G G for the existence of an LSM anomaly. When this condition is not met, we prove an SPT-LSM theorem: any unique and gapped ground state is necessarily a non-invertible weak symmetry protected topological (SPT) state with non-trivial entanglement, for which we construct an example fixed-point Hamiltonian. The projectivity also affects the dual symmetries after gauging \mathsf{Rep}(G)× Z(G) 𝖱 𝖾 𝗉 ( G ) × Z ( G ) sub-symmetries, giving rise to non-Abelian and non-invertible dipole symmetries, as well as non-invertible translations. We complement our analysis with the SymTFT, where the projectivity causes it to be a topological order non-trivially enriched by translations. Throughout the paper, we develop techniques for gauging \mathsf{Rep}(G) 𝖱 𝖾 𝗉 ( G ) symmetry and inserting its symmetry defects on the lattice, which are applicable to other non-invertible symmetries.

Pace, Salvatore D. (ORCID:0000000306093335)↗

Conformal geometry from entanglement

In a physical system with conformal symmetry, observables depend on cross-ratios, measures of distance invariant under global conformal transformations (conformal geometry for short). We identify a quantum information-theoretic mechanism by which the conformal geometry emerges at the gapless edge of a 2+1D quantum many-body system with a bulk energy gap. We introduce a novel pair of information-theoretic quantities (\mathfrak{c}_{\textrm{tot}}, \eta) ( 𝔠 tot , η ) that can be defined locally on the edge from the wavefunction of the many-body system, without prior knowledge of any distance measure. We posit that, for a topological groundstate, the quantity \mathfrak{c}_{\textrm{tot}} 𝔠 tot is stationary under arbitrary variations of the quantum state, and study the logical consequences. We show that stationarity, modulo an entanglement-based assumption about the bulk, implies (i) \mathfrak{c}_{\textrm{tot}} 𝔠 tot is a non-negative constant that can be interpreted as the total central charge of the edge theory. (ii) \eta η is a cross-ratio, obeying the full set of mathematical consistency rules, which further indicates the existence of a distance measure of the edge with global conformal invariance. Thus, the conformal geometry emerges from a simple assumption on groundstate entanglement. We show that stationarity of \mathfrak{c}_{\textrm{tot}} 𝔠 tot is equivalent to a vector fixed-point equation involving \eta η , making our assumption locally checkable. We also derive similar results for 1+1D systems under a suitable set of assumptions.

Kim, Isaac H.↗

Equations of Motion for the Vertical Rigid-Body Rotor: Linear and Nonlinear Cases

Centuries ago, the prolific mathematician Leonhard Euler (1707–1783) wrote down the equations of motion (EOM) for the heavy symmetrical top with one point fixed. The resulting set of equations turned out to be nonlinear and had a limited number of closed-form solutions. Today, tools such as transfer matrix and finite elements enable the calculation of the rotor dynamic properties for rotor-bearing systems. Some of these tools rely on the “linearized” version of the EOM to calculate the eigenvalues, unbalance response, or transients in these systems. In fact, industry standards mandate that rotors be precisely balanced to have safe operational characteristics. However, in some cases, the nonlinear aspect of the EOM should be considered. The purpose of this paper is to show examples of how the linear vs. nonlinear formulations differ. This paper will also show how excessive unbalance is capable of dramatically altering the behavior of the system and can produce chaotic motions associated with the “jump” phenomenon.

42 ENGINEERING↗

3D Equilibrium Reconstruction and MHD Studies in Stellarators (Final Report)

This is the final technical report for the grant Three Dimensional Equilibrium Reconstruction and MHD Studies in Stellarators. The major goals of this project were to: further develop V3FIT code capabilities, continue to support experiments using the V3FIT code, continue theoretical support for the CTH experiment at Auburn, particularly in the area of NIMROD simulations of CTH discharges. Dr. Lopez Ortiz has made significant progress on his analysis of NIMROD simulations of CTH MHD phenomena. In particular, his study has focused on the fixed point dynamics and energy flow amongst toroidal Fourier modes in sawtooth oscillations simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Electricity Market Dynamics as Smooth Games: Organizing Gradient-Based Learners and Identifying Instabilities

Wide-spread adoption of machine learning and distributed energy production has the potential to fundamentally change behaviours of energy markets. In practice, along with the benefits, previously unseen instabilities have begun emerging in markets currently embracing these technologies. In this report we study electricity markets as dynamical systems arising in smooth games, where energy consumers are modelled as gradient-based learners who interact on a graph through distributed energy producers. We give simple examples of how instabilities can arise even when the consumers’ local objectives seem reasonable and how these instabilities can be remedied. Drawing insights from these examples, we begin developing a general model of optimization-driven electricity market dynamics and give conditions for ensuring local Nash equilibria of the system are stable fixed points of the dynamics. Our results demonstrate how models of electricity market dynamics might be used for deriving formal guarantees in order to stabilize markets and enact regulatory policies in the face of market shocks.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

The Connection of Adiabaticity to the Fourier Transform

This is probably the first thing I’d suggest anyone learns about adiabatic transitions because it tells you how a linear system reacts to an arbitrary forcing function. Not all systems are linear but they are usually locally linear near their stable fixed point. I think Scott Berg pointed out the relation to the Fourier transform to me and there’s probably literature about this under generic dynamical systems.

43 PARTICLE ACCELERATORS↗

Tokamak disruption simulation

UT contribution in this Tokamak Disruption Simulation (TDS) is to develop a parallel high‐order hybridized Discontinuous Galerkin (HDG) methods for large‐s cale MHD simulations. The following are the major goals: 1) Construction of HDG methods for linearized MHD, 2) Rigorous analysis for the HDG formulations for linearized MHD, 3) 2D and 3D simulations to verify the convergent of the HDG methods for linearized MHD, 4) multigrid solvers/preconditioners for HDG formulations, 5) HDG for reconnection problems, 6) Divergence cleaning with HDG, 7) HDG formulations for nonlinear MHD, 8) Picard fixed point HDG approach, 9) IMEX HDG‐DG for nonlinear HDG; 10) parallel large‐scale HDG for linear and nonlinear MHD simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Flavor Fragmentation Function Factorization

A definition of partonic jet flavor that is both theoretically well-defined and experimentally robust would have profound implications for measurements and predictions especially for heavy flavor applications. Recently, a definition of jet flavor was introduced as the net flavor flowing along the direction of the Winner-Take-All axis of a jet which is soft safe to all orders, but not collinear safe. Here, we exploit the lack of collinear safety and propose a factorization theorem of perturbative flavor fragmentation functions that resum collinear divergences and describe the evolution of flavor from the short distance of jet production to the long distance at which hadronization occurs. Collinear flavor evolution is governed by a small modification of the DGLAP equations. We present a detailed all-orders analysis and identify exact relations that must hold amongst the various anomalous dimensions by probability conservation and the existence of fixed points of the renormalization group flow. We explicitly validate the factorization theorem at one-loop order, and demonstrate its consistency at two loops in particular flavor channels. Starting at two-loops, constraints on phase space imposed by flavor measurements potentially allow for non-trivial soft contributions, but we demonstrate that they are scaleless and so explicitly vanish, ensuring that soft particles are summed inclusively and all divergences are exclusively collinear in nature. This factorization theorem opens the door to precision calculations with identified flavor in the infrared.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗