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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Competing orders at higher-order Van Hove points

Van Hove points are special points in the energy dispersion, where the density of states exhibits analytic singularities. When a Van Hove point is close to the Fermi level, tendencies towards density wave orders, Pomeranchuk orders, and superconductivity can all be enhanced, often in more than one channel, leading to a competition between different orders and unconventional ground states. In this study we consider the effects from higher-order Van Hove points, around which the dispersion is flatter than near a conventional Van Hove point, and the density of states has a power-law divergence. We argue that such points are present in intercalated graphene and other materials. We use an effective low-energy model for electrons near higher-order Van Hove points and analyze the competition between different ordering tendencies using an unbiased renormalization-group approach. For purely repulsive interactions, we find that two key competitors are ferromagnetism and chiral superconductivity. For a small attractive interaction, we find an unconventional spin Pomeranchuk order, in wich the spin oder parameter winds around the Fermi surface. The supermetal state, predicted for a single higher-order Van Hove point, is an unstable fixed point in our case.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Compressed basis GMRES on high-performance graphics processing units

Krylov methods provide a fast and highly parallel numerical tool for the iterative solution of many large-scale sparse linear systems. To a large extent, the performance of practical realizations of these methods is constrained by the communication bandwidth in current computer architectures, motivating the investigation of sophisticated techniques to avoid, reduce, and/or hide the message-passing costs (in distributed platforms) and the memory accesses (in all architectures). This article leverages Ginkgo’s memory accessor in order to integrate a communication-reduction strategy into the (Krylov) GMRES solver that decouples the storage format (i.e., the data representation in memory) of the orthogonal basis from the arithmetic precision that is employed during the operations with that basis. Given that the execution time of the GMRES solver is largely determined by the memory accesses, the cost of the datatype transforms can be mostly hidden, resulting in the acceleration of the iterative step via a decrease in the volume of bits being retrieved from memory. Together with the special properties of the orthonormal basis (whose elements are all bounded by 1), this paves the road toward the aggressive customization of the storage format, which includes some floating-point as well as fixed-point formats with mild impact on the convergence of the iterative process. We develop a high-performance implementation of the “compressed basis GMRES” solver in the Ginkgo sparse linear algebra library using a large set of test problems from the SuiteSparse Matrix Collection. We demonstrate robustness and performance advantages on a modern NVIDIA V100 graphics processing unit (GPU) of up to 50% over the standard GMRES solver that stores all data in IEEE double-precision.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Analysis of Vector Particle-In-Cell (VPIC) memory usage optimizations on cutting-edge computer architectures

Vector Particle-In-Cell (VPIC) is one of the fastest plasma simulation codes in the world, with particle numbers ranging from one trillion on the first petascale system, Roadrunner, to ten trillion particles on the more recent Blue Waters supercomputer. As supercomputers continue to grow rapidly in size, so too does the gap between computing capability and memory capability. Current memory systems limit VPIC simulations greatly as the maximum number of particles that can be simulated directly depends on the available memory. In this study, we present a suite of VPIC memory optimizations (i.e., particle weight, half-precision, and fixed-point optimizations) that enable a significant increase in the number of particles in VPIC simulations. Here, we assess the optimizations’ impact on memory and runtime performance for a suite of cutting-edge computer architectures such has the NVIDIA V100 GPU, the IBM Power9, and the Fujitsu A64FX architectures. Our optimizations enable a 31.25% reduction in memory usage and up to 40% increase in the number of particles. This paper extends our work on developing particle storage format optimizations Tan et al.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Nuclear Materials Process Modeling at the Y-12 National Security Complex

The Y-12 National Security Complex (Y-12) has implemented process modeling for various accountable nuclear materials operations that are performed throughout the plant. Using a discrete, event-based dynamic simulation program, key nuclear material streams are modeled, allowing Y-12 to effectively manage numerous points of interest within the plant’s production operations. Integration of the various material processes into a single, interdependent supply and demand model is one of the ongoing focuses within Y-12’s process modeling effort. The primary purpose of using dynamic simulation modeling is to allow for analysis of the nuclear materials inventories and forecasted supplies based on future demands. Analysis of these inventories includes capacity evaluation, bottleneck mitigation, and assessments of individual pieces of equipment to inform future facility investment decisions and associated project schedules. Modeling of the nuclear materials processes throughout the complex also allows for incorporation of changes relevant to production capabilities such as the upcoming transition of specific operations to the new Uranium Processing Facility. Prior to implementation of process modeling, Y-12 forecasted supply and demand of accountable nuclear materials streams using Microsoft Excel. With deterministic models such as Microsoft Excel, the annual forecasts, generated within data input condition parameters, can only provide a fixed point of data. Fixed data cannot simulate integrated material streams and account for the possibility of occurrences and other changes that dynamic simulations take into consideration. Y-12’s dynamic process modeling allows integrated simulations of multiple accountable nuclear materials processes, including supply and demand forecasting and analysis, and is a coordinated effort involving many steps of verification and validation (V&V), site briefings, testing, reporting, data mining, planning, and documentation that spans various programs throughout the Y-12 complex.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Spontaneous breaking of U(1) symmetry in coupled complex SYK models

As shown in [1], two copies of the large N Majorana SYK model can produce spontaneous breaking of a Z 2 symmetry when they are coupled by appropriate quartic terms. In this paper we similarly study two copies of the complex SYK model coupled by a quartic term preserving the U(1) x U(1) symmetry. We also present a tensor counterpart of this coupled model. When the coefficient α of the quartic term lies in a certain range, the coupled large N theory is nearly conformal. We calculate the scaling dimensions of fermion bilinear operators as functions of α. We show that the operator c$^†_{1i}$c 2i , which is charged under the axial U(1) symmetry, acquires a complex dimension outside of the line of fixed points. We derive the large N Dyson-Schwinger equations and show that, outside the fixed line, this U(1) symmetry is spontaneously broken at low temperatures because this operator acquires an expectation value. We support these findings by exact diagonalizations extrapolated to large N.

1/N expansion↗
DOE OSTI2021

Impact of Dry Cooler Air-Side Performance on a sCO2 Power Cycle for a CSP Application

Uncertainty around the design and control of the supercritical CO2 power cycle must be reduced before this technology can be implemented for large-scale grid support. To better understand the day-to-day performance of an sCO2 cycle, off-design performance calculations must be included for all power block components, and performance assumptions must be removed. This study has expanded the modeled scope to include the air-side performance for the dry cooler and has incorporated discretized heat transfer calculations for both streams through the pre-cooler to better predict off-design performance. This study considered a recompression Brayton cycle in a concentrating solar power application. The cycle model utilized fixed sCO2 turbomachinery maps for the main compressor, recompressor, and expander operating to supply approximately 10 MW gross at the design point. Fixed vendor-supplied fan curves were used to calculate the air-side performance of the dry cooler. The primary heater was modeled considering both the sCO2 and heat transfer fluid streams. Off-design performance was predicted for an ambient temperature range of 0-55℃, a HTF temperature range of 705-735℃, and a HTF mass flow range of 50-105% of the design point value. To understand the importance of modeling the air-side performance, the cycle off-design performance was also calculated using a constant CO2 outlet temperature assumption and a constant approach temperature assumption for the dry cooler. Results show that using these assumptions can significantly alter the power output and cycle efficiency predictions.

14 SOLAR ENERGY↗
DOE OSTI2020

Subcritical Multiplication with a Fixed Source

In a subcritical, multiplying medium, the system multiplication describes the expected total number of neutrons created by a single source neutron. Subcriticality plays a large role in criticality safety and thus it is vital for the subcritical multiplication factor be accurate, especially as a system approaches criticality. This work examines the accuracy of calculating the system multiplication using the MCNP6.2 ® k-eigenvalue power iteration (KCODE) method when a fixed-point source is present in a multiplying medium, for near critical systems. This work compares the standard approach for calculating system multiplication, using the fixed-source calculational approach, to a new, single k-eigenvalue power iteration approach that incorporates a fixed-source component and a fission-source component into a single calculation. For the remainder of this paper, some theoretical background and numerical results for an approximate k eigenvalue approach, an accurate fixed-source approach and a new and more accurate k-eigenvalue approach to computing system multiplication are provided.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Slow scrambling and hidden integrability in a random rotor model

In this work, we analyze the out-of-time-order correlation functions of a solvable model of a large number $\textit{N}$ of $\textit{M}$-component quantum rotors coupled by Gaussian-distributed random, infinite-range exchange interactions. We focus on the growth of commutators of operators at a temperature $\textit{T}$ above the zero temperature quantum critical point separating the spin-glass and paramagnetic phases. In the large $\textit{N, M}$ limit, the squared commutators of the rotor fields do not display any exponential growth of commutators, in spite of the absence of any sharp quasiparticlelike excitations in the disorder-averaged theory. We show that in this limit, the problem is integrable and points out interesting connections to random-matrix theory. At leading order in 1/$\textit{M}$, there are no modifications to the critical behavior but an irrelevant term in the fixed-point action leads to a small exponential growth of the squared commutator. We also introduce and comment on a generalized model involving $\textit{p}$-pair rotor interactions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Convergence analysis of single rate and multirate fixed stress split iterative coupling schemes in heterogeneous poroelastic media

Recently, the accurate modeling of flow–structure interactions has gained more attention and importance for both petroleum and environmental engineering applications. Of particular interest is the coupling between subsurface flow and reservoir geomechanics. Different single rate and multirate iterative and explicit coupling schemes have been proposed and analyzed in the past. In addition, Banach fixed point contraction results were obtained for iterative coupling schemes, and conditionally stable results were obtained for explicit coupling schemes. In this work, we will consider the mathematical analysis of the single rate and multirate fixed stress split iterative coupling schemes for spatially heterogeneous poroelastic media. We will re–establish the contractivity for both schemes in the localized case, and we will show that heterogeneities come at the expense of imposing more restricted conditions on the number of fine flow time steps that can be taken within one coarse mechanics time step in the multirate case. Our mathematical analysis is supplemented by numerical simulations validating our derived upper bounds. Finally, to the best of our knowledge, this is the first rigorous mathematical analysis of the multirate fixed–stress split iterative coupling scheme in heterogeneous poroelastic media.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Superadditivity at large charge

The weak gravity conjecture has been invoked to conjecture that the dimensions of charged operators in a CFT should obey a superadditivity relation (sometimes referred to as convexity). In this paper, we study superadditivity of the operator spectrum in theories expanded about the semi-classical saddle point that dominates correlators of large charge operators. We explore this in two contexts. The first is a model with two scalar fields that carry different charges, at a non-trivial Wilson-Fisher fixed point. A careful analysis of the semi-classics for this two field model demonstrates that ‘quantum’ violations of superadditivity (those not forbidden by the conjecture) persist in the large charge regime. We then turn to study the general properties of CFTs at large charge as bottom-up EFTs. By a trial and error procedure we come up with a seemingly consistent family of examples violating the conjecture. In so doing the presence of a genuine dilaton field appears necessary. On the one hand our result demonstrates that the superadditivity conjecture cannot be proven purely on the basis of a bottom-up analysis. On the other hand, the need for a dilaton, with the corresponding infinite fine tuning, indicates the conjecture-violating EFTs are unlikely to be UV completable.

effective field theories↗
DOE OSTI2020

Interaction-Induced Metallicity in a Two-Dimensional Disordered Non-Fermi Liquid

The interplay of interactions and disorder in two-dimensional (2D) electron systems has actively been studied for decades. The paradigmatic approach involves starting with a clean Fermi liquid and perturbing the system with both disorder and interactions. Here, we start instead with a clean non-Fermi liquid near a 2D ferromagnetic quantum critical point and consider the effects of disorder. In contrast with the disordered Fermi liquid, we find that our model does not suffer from runaway flows to strong coupling and the system has a marginally stable fixed point with perfect conduction.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Nonequilibrium universality of the nonreciprocally coupled 𝑂⁡(𝑛 1 ) × 𝑂⁡(𝑛 2 ) model

Nonequilibrium dynamics play an important role in all contexts of physics, both classical and quantum as well as living and nonliving, so it is crucial to develop a foundational understanding of nonequilibrium phase transitions. In this work we investigate an important class of nonequilibrium dynamics in the form of nonreciprocal interactions. In particular we study how nonreciprocal coupling between two 𝑂⁡(𝑛𝑖) order parameters (with 𝑖 = 1,2) affects the universality at a multicritical point, extending the analysis of J. T. Young et al. [Phys. Rev. X 10, 011039 (2020)], which considered the case 𝑛 1 = 𝑛 2 = 1, i.e., a ℤ 2 × ℤ 2 model. We show that nonequilibrium fixed points (NEFPs) emerge for a broad range of 𝑛 1 ,𝑛 2 and exhibit intrinsically nonequilibrium critical phenomena, namely a violation of fluctuation-dissipation relations at all scales and underdamped oscillations near criticality in contrast to the overdamped relaxational dynamics of the corresponding equilibrium models. Furthermore, the NEFPs exhibit an emergent discrete scale invariance in certain physically relevant regimes of 𝑛 1 ,𝑛 2 , but not others, depending on whether the critical exponent 𝜈 is real or complex. The boundary between these two regions is described by an exceptional point in the renormalization group (RG) flow, leading to distinctive features in correlation functions and the phase diagram. Another contrast with the previous work is the number and stability of the NEFPs as well as the underlying topology of the RG flow. Lastly, we investigate an extreme form of nonreciprocity where one order parameter is independent of the other order parameter but not vice versa. Unlike the ℤ 2 × ℤ 2 model, which becomes nonperturbative in this case, we identify a distinct nonequilibrium universality class whose dependent field similarly violates fluctuation-dissipation relations but does not exhibit discrete scale invariance or underdamped oscillations near criticality.

Critical phenomena↗
DOE OSTI2023

The Gross-Neveu-Yukawa archipelago

We perform a bootstrap analysis of a mixed system of four-point functions of bosonic and fermionic operators in parity-preserving 3d CFTs with O(N) global symmetry. Our results provide rigorous bounds on the scaling dimensions of the O(N)-symmetric Gross-Neveu-Yukawa (GNY) fixed points, constraining these theories to live in isolated islands in the space of CFT data. We focus on the cases N = 1, 2, 4, 8, which have applications to phase transitions in condensed matter systems, and compare our bounds to previous analytical and numerical results.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

The myth of URANS

Since the 1990s, RANS practitioners have observed spontaneous unsteadiness in RANS simulations. Some have suggested deliberately using this as a method of resolving large turbulent structures. However, to date, no one has produced a theoretical justification for this unsteady RANS (URANS) approach. Here, we extend the dynamical system fixed point analysis to create a theoretical model for URANS dynamics. The results are compared to URANS simulations for homogeneous isotropic decaying turbulence. The model shows that URANS can predict incorrect decay rates and that the solution tends towards steady RANS over time. Similar analysis for forced turbulence shows a fixed modelled energy of about 30% of total energy, regardless of the model parameters. The same analysis can be used to show how hybrid type models can begin to address these issues.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Properties of the new $ \mathcal{N} $ = 1 AdS 4 vacuum of maximal supergravity

The recent comprehensive numerical study of critical points of the scalar potential of four-dimensional N = 8, SO(8) gauged supergravity using Machine Learning software in [1] has led to a discovery of a new N = 1 vacuum with a triality-invariant SO(3) symmetry. Guided by the numerical data for that point, we obtain a consistent SO(3) × Z 2 -invariant truncation of the N = 8 theory to an N = 1 supergravity with three chiral multiplets. Critical points of the truncated scalar potential include both the N = 1 point as well as two new non-supersymmetric and perturbatively unstable points not found by previous searches. Studying the structure of the submanifold of SO(3) × Z 2 -invariant supergravity scalars, we find that it has a simple interpretation as a submanifold of the 14-dimensional Z $3\atop{2}$ -invariant scalar manifold (SU(1, 1)/U(1)) 7 , for which we find a rather remarkable superpotential whose structure matches the single bit error correcting (7, 4) Hamming code. This 14-dimensional scalar manifold contains approximately one quarter of the known critical points. We also show that there exists a smooth supersymmetric domain wall which interpolates between the new N = 1 AdS4 solution and the maximally supersymmetric AdS4 vacuum. Using holography, this result indicates the existence of an N = 1 RG flow from the ABJM SCFT to a new strongly interacting conformal fixed point in the IR.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Many-body localization near the critical point

Here, we examine the many-body localization (MBL) phase transition in one-dimensional quantum systems with quenched randomness and short-range interactions. Following recent works, we use a strong-randomness renormalization group (RG) approach where the phase transition is due to the so-called avalanche instability of the MBL phase. We show that the critical behavior can be determined analytically within this RG. On a rough qualitative level the RG flow near the critical fixed point is similar to the Kosterlitz-Thouless (KT) flow as previously shown, but there are important differences in the critical behavior. Thus, we show that this MBL transition is in a universality class that is different from KT. The divergence of the correlation length corresponds to critical exponent $\textit{ν → ∞}$, but the divergence is weaker than for the KT transition.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Interaction-induced velocity renormalization in magic-angle twisted multilayer graphene

Abstract Twistronics heterostructures provide a novel route to control the electronic single particle velocity and thereby to engineer strong effective interactions. Here we show that the reverse may also hold, i.e. that these interactions strongly renormalize the band structure. We demonstrate this mechanism for alternating-twist magic-angle three- and four-layer graphene at charge neutrality and in the vicinity of a phase transition which can be described by an Ising Gross-Neveu critical point corresponding, e.g. to the onset of valley Hall or Hall order. While the non-interacting model displays massless Dirac excitations with strongly different velocities, we show that interaction corrections make them equal in the infrared. However, the renormalization group flow of the velocities and of the coupling to the critical bosonic mode is strongly non-monotonic and dominated by the vicinity of a repulsive fixed point. We predict experimental consequences of this theory for tunneling and transport experiments and discuss the expected behavior at other quantum critical points, including those corresponding to intervalley coherent ordering.

Classen, Laura↗
DOE OSTI2020

Quantum Criticality in the 2D Quasiperiodic Potts Model

Quantum critical points in quasiperiodic magnets can realize new universality classes, with critical properties distinct from those of clean or disordered systems. Here, we study quantum phase transitions separating ferromagnetic and paramagnetic phases in the quasiperiodic q-state Potts model in 2+1D. Using a controlled real-space renormalization group approach, we find that the critical behavior is largely independent of q, and is controlled by an infinite-quasiperiodicity fixed point. In conclusion, the correlation length exponent is found to be ν = 1, saturating a modified version of the Harris-Luck criterion.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗