Search NASASearch

SEARCH · Search NASA

Results for “ASYMPTOTIC FUNCTION”

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 19 records

Computation of generalised magnetic coordinates asymptotically close to the separatrix

Integrals to calculate generalised magnetic coordinates from an input magnetic flux function asymptotically close to the separatrix are presented, and implemented in the GPEC/DCON code suite. These integrals allow characterisation of the magnetic equilibrium of a diverted tokamak, in magnetic coordinates, arbitrarily close to the last closed flux surface, avoiding the numerical issues associated with calculating diverging field-line integrals near a magnetic x-point. Finally, these methods may assist ongoing efforts to develop robust asymptotic equilibrium behaviour for spectral 3D MHD codes at the separatrix.

equilibrium edge truncation

On the Sampling-Based Computation of Nash Equilibria Under Uncertainty via the Nikaido–Isoda Function

We consider the computation of an equilibrium of a stochastic Nash equilibrium problem, where the player objectives are assumed to be L 0 -Lipschitz continuous and convex, given rival decisions with convex and closed player-specific feasibility sets. To address this problem, we consider minimizing a suitably defined value function defined using the Nikaido–Isoda function. Such an avenue does not necessitate either monotonicity properties of the concatenated gradient map or potentiality requirements on the game but does require a suitable regularity requirement under which a stationary point is a Nash equilibrium. We design and analyze a sampling-enabled projected-gradient-response method, reliant on inexact resolution of a player-level best-response subproblem. Here, by deriving suitable Lipschitzian guarantees on the value function, we derive both asymptotic guarantees for the sequence of generated iterates as well as rate and complexity guarantees for computing a stationary point by appropriate choices of the sampling rate and inexactness sequence.

Nikaido-Isoda function

Linear-scaling quadruple excitations in local pair natural orbital coupled-cluster theory

Here, we present a fast, asymptotically linear-scaling implementation of the perturbative quadruples energy correction in coupled-cluster theory using local natural orbitals. Our work follows the domain-based local pair natural orbital (DLPNO) approach previously applied to lower levels of excitations in coupled-cluster theory. Our DLPNO-CCSDT(Q) algorithm uses converged doubles and triples amplitudes from a preceding DLPNO-CCSDT computation to compute the quadruples amplitude and energy in the quadruples natural orbital (QNO) basis. We demonstrate the compactness of the QNO space, showing that more than 95% of the (Q) correction can be recovered using relatively loose natural orbital cutoffs, compared to the tighter cutoffs used in pair and triples natural orbitals at lower levels of coupled-cluster theory. We also highlight the accuracy of our algorithm in the computation of relative energies, which yields deviations of sub-kJ mol −1 in relative energy compared to the canonical CCSDT(Q). Timings are conducted on a series of growing linear alkanes (up to 10 carbons and 608 basis functions) and water clusters (up to 49 water molecules and 2842 basis functions) to establish the asymptotic linear-scaling of our DLPNO-(Q) algorithm.

Auxiliary functions

Empirical investigation of nuclear correlation function distributions in lattice QCD

Two-point correlation functions of systems with baryon number 𝐵 ∈ {1,2,3,4} are investigated using lattice quantum chromodynamics (QCD). In particular, the empirical distributions of importance-sampling Monte-Carlo samples of these correlation functions are examined as a function of the spacetime separation between the two points and the baryon number. While the exact forms of these distributions are not known for QCD, recent work has determined asymptotic expressions for analogous correlation function distributions in simpler theories such as scalar field theory and the disordered phase of the 𝑂⁡(𝑁) model. The theoretical 𝑂⁡(𝑁) model distributions are found to provide an accurate description of the empirical QCD distributions at zero momentum over a wide range of temporal separations for each baryon number when assessed with a range of different statistical tests. In particular, the behavior of the baryon number 𝐵 QCD correlation function at large temporal separation is well reproduced by the 𝑂⁡(𝑁 ∼ 2/𝐵) model distribution.

Lattice field theory

A novel conditional formulation of the Vlasov–Ampère equations: a conservative, positivity, asymptotic and Gauss law preserving scheme

We propose a novel reformulation of the Vlasov–Ampère equations for plasmas that reveals discrete symmetries that enables simultaneous conservation of mass, momentum and energy; preservation of Gauss’s law; positivity of the distribution function; and consistency with quasi-neutral asymptotics. The approach employs variable and coordinate transformations to yield a coupled system comprising a modified Vlasov equation and associated moment–field equations. The modified Vlasov equation advances a conditional distribution function that excludes mass, momentum and energy densities, which are instead evolved through moment equations enforcing the relevant symmetries, conservation laws and involution constraints. This reformulation aligns naturally with a recent slow-manifold reduction technique, which separates fast electron time scales and simplifies the treatment of the quasi-neutral limit within the reduced moment–field subsystem. Using this framework, we develop a numerical method for the reduced 1D1V subsystem that, for the first time in the literature, satisfies all key physical constraints while maintaining a quasi-neutral asymptotic behaviour. The advantages of the method are demonstrated on canonical electrostatic test problems, including the multiscale ion acoustic shock wave.

1D1V

An asymptotic-preserving semi-Lagrangian algorithm for the anisotropic heat transport equation with arbitrary magnetic fields

Here, we extend the recently proposed semi-Lagrangian algorithm for the extremely anisotropic heat transport equation [Chacón et al., J. Comput. Phys ., 272 (2014)] to deal with arbitrary magnetic field topologies. The original scheme (which showed remarkable numerical properties) was valid for the so-called tokamak-ordering regime, in which the magnetic field magnitude was not allowed to vary much along field lines. The proposed extension maintains the attractive features of the original scheme (including the analytical Green's function, which is critical for tractability) with minor modifications, while allowing for completely general magnetic fields. The accuracy and generality of the approach are demonstrated by numerical experiment with an analytical manufactured solution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY

Casimir-Polder potential on an excited atom near an atomic array

We develop a microscopic description of the fluctuation-mediated Casimir-Polder (CP) shifts on a 'test' two-level atom placed near a two-dimensional atomic array of two-level atoms. We derive the resonant and off-resonant CP potentials experienced by the excited test atom using fourth-order perturbation theory, under the assumption that the test atom resonance is far detuned from those of the array atoms. The total potential on the test atom can be described as the sum of the pairwise resonant and off-resonant potentials resulting from its interaction with the individual atoms of the array. We analyze the asymptotic scaling of CP shifts as a function of the test atom-array separation, and its dependence on various system parameters: array spacing and size, and dipole orientation of the array atoms. Our results bridge the description of CP potential across two distinct regimes: (i) from a single-atom limit where we recover the well-known two-atom Van der Waals potential, (ii) to a macroscopic boundary limit, where we demonstrate new asymptotic scaling laws. We demonstrate that these scaling laws can be tuned via the microscopic parameters of the atomic array, establishing atomically-controlled arrays as a versatile platform for tailoring fluctuation-induced QED phenomena.

FOS: Physical sciences

Quantum Hamilton-Jacobi theory, spectral path integrals, and exact WKB analysis

We propose a new way to perform path integrals in quantum mechanics by using a quantum version of Hamilton-Jacobi (HJ) theory. In classical mechanics, Hamilton-Jacobi theory is a powerful formalism, however, its utility is not explored in quantum theory beyond approximation schemes. The canonical transformation enables one to set the new Hamiltonian to constant or zero, but keeps the information about solution in Hamilton’s characteristic function. To benefit from this in quantum theory, one must work with a formulation in which classical Hamiltonian is used. This uniquely points to phase space path integral. However, the main variable in HJ formalism is energy, not time. Thus, we are led to consider the Fourier transform of the path integral, the spectral path integral Z ˜ ( E ) . The evaluation of path integrals reduces to determining the quantum Hamilton characteristic functions (which can be achieved via an asymptotic analysis) and a discrete sum over the quantum period lattice, generalizing Gutzwiller’s sum. Published by the American Physical Society 2025

Türe, Mustafa (ORCID:0009000975968618)

Generative Vulnerability Assessment for Cyber-Physical Systems

Cyber-physical systems (CPS) are highly susceptible to malicious attacks due to their complex dynamics and interconnectivity. A comprehensive understanding of their vulnerabilities is essential for designing effective resilience measures. This paper presents a data-driven attack generative system for evaluating the vulnerability of CPS. The proposed approach formulates the vulnerability assessment problem as determining the feasibility of a specific attack set based on two boundary functions that represent the effectiveness and stealthiness of attacks. The attack generative model is trained using a custom loss function, with two universal approximators designed to learn the effectiveness and stealthiness functions simultaneously. Theoretical results for successful generation and asymptotic convergence of the resulting training algorithm are given. As a result, the proposed approach is evaluated via numerical simulation of an IEEE 14-bus system and gas pipeline systems, demonstrating its viability in learning how to attack nonlinear CPS and identify potential vulnerabilities.

Computer systems organization

Response tailoring of elasto-plastic trusses

Abstract In this work we tailor the response of trusses loaded beyond their yield limit. The truss structures are modeled using finite strain theory and rate-independent elasto-plasticity. We design trusses with a tailored mechanical response that is between “elastic” and “elastic-ideal-plastic” subject to the volume constraint. The design updates are generated by the gradient-based Method of Moving Asymptotes (MMA) solver and the sensitivities of the response functions are computed using a path-dependent adjoint sensitivity analysis. The computations are performed in Matlab.

42 ENGINEERING

Quasifragmentation functions in the massive Schwinger model

We introduce the concept of the quark quasifragmentation function (qFF) using an equal-time and spatially boosted form of the Collins-Soper fragmentation function where the out-meson fragment is replaced by the current asymptotic condition. We derive the qFF for a fermion in two-dimensional quantum electrodynamics (QED2) using the Kogut-Susskind Hamiltonian after a mapping onto spin qubits in a spatial lattice with open boundary conditions. This form is suitable for quantum computations. We compute the qFF by exact diagonalization of the spin Hamiltonian. The results are compared to the qFF following from the Drell-Levy-Yan result for QED2, both at strong and weak coupling, and to two-dimensional quantum chromodynamics in the lowest Fock approximation.

Fragmentation functions

Small-𝑥 asymptotics of the leading-twist flavor-singlet quark TMDs

In this paper, we investigate the small-𝑥 behavior of the flavor-singlet, leading-twist quark transverse-momentum-dependent parton distribution functions (TMDs) using the light-cone operator treatment. This formalism allows us to express TMD operators at small 𝑥 in terms of polarized dipole amplitudes, enabling a systematic approach to their small-𝑥 evolution. We derive the evolution equations for these TMDs and solve them within the large-𝑁 𝑐 approximation under the linearized, double-logarithmic approximation, where 𝑁 𝑐 represents the number of quark colors. Expanding on previous work on unpolarized and helicity TMDs, we present the small-𝑥 asymptotics for a comprehensive set of TMDs, including the Sivers function, helicity worm-gear, transversity, pretzelosity, Boer-Mulders, and transversity worm-gear distributions. Our results provide a complete picture of the small-𝑥 asymptotic behavior for all leading-twist flavor-singlet quark TMDs. We also discuss the implications of our findings for phenomenological applications and outline potential avenues for further research, particularly in understanding nonlinear effects and extending beyond the double-logarithmic approximation and large-𝑁 𝑐 approximations.

Adamiak, Daniel [Thomas Jefferson National Acceler

The effective number of parameters in kernel density estimation

We devise a new formula for measuring the effective degrees of freedom (EDoF) in kernel density estimation (KDE). Starting from the orthogonal polynomial sequence (OPS) expansion for the ratio of the empirical to the oracle density, we show how convolution with the kernel leads to a new OPS with respect to which one may express the resulting KDE. The expansion coefficients of the two OPS systems can then be related via a kernel sensitivity matrix, which leads to a natural oracle definition of EDoF through the trace operator. Asymptotic properties of the (empirical) plug-in EDoF are worked out through influence functions, and connections with other empirical EDoFs are established. Minimization of Kullback-Leibler divergence is investigated as an alternative to integrated squared error based bandwidth selection rules, yielding a new normal scale rule. The methodology, which arises from a proper oracle formulation and is not restricted to convolution kernels, suggests the possibility of a new bandwidth selection rule based on an information criterion such as AIC.

bandwidth selection

Asymptotic inconsistency of the cumulative algorithm for laser-induced damage probability analysis

The “cumulative algorithm” is a data analysis method that has been proposed to provide an objective, nonparametric determination of laser-induced damage probability as a function of fluence from experimental data that contain both damaged sites and undamaged sites (i.e., 1-on-1 or S-on-1 testing protocols). In this work, the limitations of this approach are explored by considering the asymptotic limit of a large number of test sites. It is shown that the cumulative algorithm does not converge to the true probability distribution and significantly underestimates the damage probability near the damage onset. Here, based on the results of this work, the cumulative algorithm is not recommended for accurate estimation of damage probability.

Computational methods

Theoretical Assessment of the Transition Between Electron Emission Mechanisms for Nonplanar Diodes

Theoretically and computationally describing the operation of nanodiodes requires characterizing the transitions between multiple electron emission mechanisms for nanodiodes with complicated geometries. This motivates our development of techniques to determine when simplified theories for individual mechanisms suffice compared to more complete, but more computationally expensive, models. Leveraging recent theories that define a canonical gap distance to translate planar theory to nonplanar diodes, we derive the conditions for the transitions among thermal emission, field emission, and space-charge-limited current density (SCLCD) in vacuum and with collisions for non-Cartesian coordinate systems, including spherical, cylindrical, and prolate spheroidal coordinate systems. Particle-in-cell (PIC) simulations of the current density as a function of applied voltage for a tip-to-plate geometry in vacuum agreed qualitatively with the asymptotes for thermal emission at low voltage and SCLCD at higher voltage using the canonical gap distance. As a result, this demonstrates the utility of this approach for guiding system design and suggests future extensions to save simulation time for more realistic geometries that are more computationally expensive.

Conformal mapping

An Exploration of the Equation of State Dependence of Core-collapse Supernova Explosion Outcomes and Signatures

We explore, using a state-of-the-art simulation code in 3D and late-enough times to witness final observables, the dependence of core-collapse supernova explosions on the nuclear equation of state (EOS). Going beyond questions of explodability, we compare final explosion energies, nucleosynthetic yields, recoil kicks, and gravitational-wave and neutrino signatures using the SFHo and DD2 nuclear EOSs for a 9 M ⊙ /solar-metallicity progenitor star. The DD2 EOS is stiffer and has a lower effective nucleon mass. The result is a more extended protoneutron star (PNS) and lower central densities. As a consequence, the mean neutrino energies, final explosion energy, and recoil kick speed are lower. Moreover, the evolution of PNS convection differs between the two EOS models in significant ways. This translates in part into interestingly altered neutrino “light” curves and noticeably altered gravitational-wave signal strengths and frequency characteristics that may be diagnostic. The faster exploding model (SFHo) yields slightly more neutron-rich ejecta and more species with atomic weights between 60 and 90 and a weak r-process. However, this is merely a preliminary study. The next step is a more comprehensive and multiprogenitor set of 3D supernova simulations for various EOSs to late times when the observables have asymptoted. Such a future investigation will have a direct bearing on the neutron star and black hole birth mass functions and the quest toward a fully quantitative theory of supernova observables.

Rusakov, Aleksandr [Princeton University, NJ (Unit

Dynamically Learning Incentives for Load Control

As electrical generation becomes more distributed and volatile, and loads become more uncertain, controllability of distributed energy resources (DERs), regardless of their ownership status, will be necessary for grid reliability. Grid operators lack direct control over end-users' grid interactions, such as energy usage, but incentives can influence behavior -- for example, an end-user that receives a grid-driven incentive may adjust their consumption or expose relevant control variables in response. A key challenge in studying such incentives is the lack of data about human behavior, which usually motivates strong assumptions, such as distributional assumptions on compliance or rational utility-maximization. In this paper, we propose a general incentive mechanism in the form of a constrained optimization problem -- our approach is distinguished from prior work by modeling human behavior (e.g., reactions to an incentive) as an arbitrary unknown function. We propose feedback-based optimization algorithms to solve this problem that each leverage different amounts of information and/or measurements. We show that each converges to an asymptotically stable incentive with (near)-optimality guarantees given mild assumptions on the problem. Finally, we evaluate our proposed techniques in voltage regulation simulations on standard test beds. We test a variety of settings, including those that break assumptions required for theoretical convergence (e.g., convexity, smoothness) to capture realistic settings. In this evaluation, our proposed algorithms are able to find near-optimal incentives even when the reaction to an incentive is modeled by a theoretically difficult (yet realistic) function.

demand response

Data-Driven Recommendation of Optimal Tuning Scheme for Range-Separated Hybrid Functionals in Solution-Phase UV/Vis Absorption Energy Prediction

Time-dependent density functional theory (TDDFT) combined with range-separated hybrid (RSH) functionals and a tuned range-separation parameter γ offers a computationally economical approach for high-throughput excited- state property predictions. The γ-tuning procedure in the gas phase is well established. However, no agreement on the best γ- tuning procedure has been made when considering the solvent effect with implicit solvent models like the polarizable continuum model (PCM). To answer that question, this study created a diverse dataset with 937 molecules with experimental solutionphase UV/vis absorption spectra. Three γ-tuning methods, the gasphase γ-tuning (GPγT), the partial vertical γ-tuning (PVγT), and the strict vertical γ-tuning (SVγT), were evaluated for the ωPBEh functional over the entire dataset. Additional benchmarks are done for the optimally tuned screened range-separated hybrid combined with the PCM approach (SRSH-PCM) and the solvation-mediated tuning procedure (sol-med-OT). Our findings revealed that the optimal γ-values obtained by the PVγT and the SVγT are significantly smaller than the GPγT. This trend holds consistently across all molecules in our dataset, and we explained the origin of this phenomenon. TDDFT calculations with PVγTand SVγT-tuned γ-values and default global Fock exchange fraction achieve superior performance compared to those using GPγTtuned or default γ and slightly outperform SRSH-PCM and sol-med-OT with similar or lesser computational cost. Furthermore, we found that the smaller γ-values from SVγT captured the expected 1/(εR) asymptotic behavior in the solution phase, resulting in accurate prediction of solution-phase CT excitations, consistent with the screened asymptote behavior encoded in SRSH-PCM. These results show that SVγT is the best scheme for high-throughput UV/vis absorption spectrum calculations using the ωPBEh functional from a data-driven perspective.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH