Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “regularity”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

The Simons Observatory: science goals and forecasts for the enhanced Large Aperture Telescope

We describe updated scientific goals for the wide-field, millimeter-wave survey that will be produced by the Simons Observatory (SO). Significant upgrades to the 6-meter SO Large Aperture Telescope (LAT) are expected to be complete by 2028, and will include a doubled mapping speed with 30,000 new detectors and an automated data reduction pipeline. In addition, a new photovoltaic array will supply most of the observatory's power. The LAT survey will cover about 60% of the sky at a regular observing cadence, with five times the angular resolution and ten times the map depth of the Planck satellite. The science goals are to: (1) determine the physical conditions in the early universe and constrain the existence of new light particles; (2) measure the integrated distribution of mass, electron pressure, and electron momentum in the late-time universe, and, in combination with optical surveys, determine the neutrino mass and the effects of dark energy via tomographic measurements of the growth of structure at redshifts z ≲ 3; (3) measure the distribution of electron density and pressure around galaxy groups and clusters, and calibrate the effects of energy input from galaxy formation on the surrounding environment; (4) produce a sample of more than 30,000 galaxy clusters, and more than 100,000 extragalactic millimeter sources, including regularly sampled AGN light-curves, to study these sources and their emission physics; (5) measure the polarized emission from magnetically aligned dust grains in our Galaxy, to study the properties of dust and the role of magnetic fields in star formation; (6) constrain asteroid regoliths, search for Trans-Neptunian Objects, and either detect or eliminate large portions of the phase space in the search for Planet 9; and (7) provide a powerful new window into the transient universe on time scales of minutes to years, concurrent with observations from the Vera C. Rubin Observatory of overlapping sky.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Global stellarator coil optimization with quadratic constraints and objectives

Most present stellarator designs are produced by costly two-stage optimization: the first for an optimized equilibrium, and the second for a coil design reproducing its magnetic configuration. Few proxies for coil complexity and forces exist at the equilibrium stage. Rapid initial state finding for both stages is a topic of active research. Most present convex coil optimization codes use the least square winding surface method by Merkel (NESCOIL), with recent improvements in conditioning, regularization, sparsity, and physics objectives. While elegant, the method is limited to modeling the norms of linear functions in coil current. We present QUADCOIL, a global coil optimization method that targets combinations of linear and quadratic functions of the current. It can directly constrain and/or minimize a wide range of physics objectives unavailable in NESCOIL and REGCOIL, including the Lorentz force, magnetic energy, curvature, field-current alignment, and the maximum density of a dipole array. QUADCOIL requires no initial guess and runs nearly $10$ 2 x faster than filament optimization. Integrating it in the equilibrium optimization stage can potentially exclude equilibria with difficult-to-design coils, without significantly increasing the computation time per iteration. QUADCOIL finds the exact, global minimum in a large parameter space when possible, and otherwise finds a well-performing approximate global minimum. It supports most regularization techniques developed for NESCOIL and REGCOIL. We demonstrate QUADCOIL’s effectiveness in coil topology control, minimizing non-convex penalties, and predicting filament coil complexity with three numerical examples.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Neural posterior unfolding

Differential cross section measurements are the currency of scientific exchange in particle and nuclear physics. A key challenge for these analyses is the correction for detector distortions, known as deconvolution or unfolding. Binned unfolding of cross section measurements traditionally rely on the regularized inversion of the response matrix that represents the detector response, mapping pre-detector (`particle level') observables to post-detector (`detector level') observables. In this paper we introduce Neural Posterior Unfolding, a modern, Bayesian approach that leverages normalizing flows for unfolding. By using normalizing flows for neural posterior estimation, NPU offers several key advantages including implicit regularization through the neural network architecture, fast amortized inference that eliminates the need for repeated retraining, and direct access to the full uncertainty in the unfolded result. In addition to introducing NPU, we implement a classical Bayesian unfolding method called Fully Bayesian Unfolding (FBU) in modern Python so it can also be studied. These tools are validated on simple Gaussian examples and then tested on simulated jet substructure examples from the Large Hadron Collider (LHC). We find that the Bayesian methods are effective and worth additional development to be analysis ready for cross section measurements at the LHC and beyond.

Analysis and statistical methods↗
DOE OSTI2021

The zero-emissions cost of energy: a policy concept

The energy sector generates over 70% of global greenhouse gas (GHG) emissions, but existing energy-climate policies do little to reduce GHG emissions to prevent climate change. We present a new carbon tax policy concept in which energy users would be taxed an amount equal to the cost of cleaning up the emissions that they create, wherein the tax revenues would be used to operate negative emissions technologies. The policy is based on applying the zero-emissions cost (ZEC) metric to energy prices, which includes the 'regular' cost of an energy source plus the cost of sequestering its emissions. Using the ZEC, (a) the energy sector would be emissions-neutral, (b) biofuels would be cheaper than petroleum fuels, and (c) renewable electricity would be cheaper than fossil-fuel electricity. We calculate the energy spending as a fraction of gross domestic product for 'regular' and ZEC cost scenarios. Implementing the ZEC carbon-tax policy would enable a monumental shift to a net-zero emissions energy sector, but carries a significant risk of causing an economic recession. The summary point of this commentary is to demonstrate a pathway to achieve net-zero GHG emissions from the US energy sector, and to contemplate the associated economic impacts on our society.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2020

Time-dependent relaxation of observables in complex quantum systems

Abstract We consider time-dependent relaxation of observables in quantum systems of chaotic and regular type. Using statistical arguments and exact numerical solutions we show that the spread of the initial wave function in the Hilbert space and the main characteristics of evolution of observables have certain generic features. The study compares examples of regular dynamics, a completely chaotic case of the Gaussian orthogonal ensemble, a bosonic system with random interactions, and a fully realistic case of the time evolution of various initial non-stationary states in the nuclear shell model. In the case of the Gaussian orthogonal ensemble we show that the survival probability obtained analytically also fully defines the relaxation timescale of observables. This is not the case in general. Using the realistic nuclear shell model and the quadrupole moment as an observable we demonstrate that the relaxation time is significantly longer than defined by the survival probability of the initial state. The full analysis does not show the presence of an analog of the Lyapunov exponent characteristic for examples of classical chaos.

Volya, Alexander (ORCID:0000000217656466)↗
DOE OSTI2021

Joint ptycho-tomography with deep generative priors

Abstract Joint ptycho-tomography is a powerful computational imaging framework to recover the refractive properties of a 3D object while relaxing the requirements for probe overlap that is common in conventional phase retrieval. We use an augmented Lagrangian scheme for formulating the constrained optimization problem and employ an alternating direction method of multipliers (ADMM) for the joint solution. ADMM allows the problem to be split into smaller and computationally more efficient subproblems: ptychographic phase retrieval, tomographic reconstruction, and regularization of the solution. We extend our ADMM framework with plug-and-play (PnP) denoisers by replacing the regularization subproblem with a general denoising operator based on machine learning. While the PnP framework enables integrating such learned priors as denoising operators, tuning of the denoiser prior remains challenging. To overcome this challenge, we propose a denoiser parameter to control the effect of the denoiser and to accelerate the solution. In our simulations, we demonstrate that our proposed framework with parameter tuning and learned priors generates high-quality reconstructions under limited and noisy measurement data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

The relation between quasars’ optical spectra and variability

Abstract Brightness variation is an essential feature of quasars, but its mechanism and relationship to other physical quantities are not understood well. We aimed to find the relationship between the optical variability and spectral features to reveal the regularity behind the random variation. It is known that a quasar’s Fe ii/Hβ flux ratio and equivalent width of [O iii]5007 are negatively correlated; this is called Eigenvector 1. In this work, we visualized the relationship between the position on this Eigenvector 1 (EV1) plane and how the brightness of the quasars had changed after ∼10 yr. We conducted three analyses, using a different quasar sample in each. The first analysis showed the relation between the quasars’ distributions on the EV1 plane and how much they had changed brightness, using 13438 Sloan Digital Sky Survey quasars. This result shows how brightness changes later are clearly related to the position on the EV1 plane. In the second analysis, we plotted the sources reported as “changing-look quasars” (or “changing-state quasars”) on the EV1 plane. This result shows that the position on the EV1 plane corresponds to the activity level of each source, and the bright or dim states of them are distributed on the opposite sides divided by the typical quasar distribution. In the third analysis, we examined the transition vectors on the EV1 plane using sources with multiple-epoch spectra. This result shows that the brightening and dimming sources move on a similar path and they reach a position corresponding to the opposite activity level. We also found this trend is opposite to the empirical rule that $R_{\rm {Fe\, \small {II}}}$ positively correlated with the Eddington ratio, which has been proposed based on the trends of a large number of quasars. From all these analyses, it is indicated that quasars tend to oscillate between both sides of the distribution ridge on the EV1 plane; each of them corresponds to a dim state and a bright state. This trend in optical variation suggests that significant brightness changes, such as changing-look quasars, are expected to repeat.

Astronomy & Astrophysics↗
DOE OSTI2026

Bound states of Ω baryons in light nuclei

Here, we investigate bound states of light Ω 3⁢𝑥 clusters (𝑥=𝑠,𝑐), motivated by the Ω 3⁢𝑠 ⁢𝑁 potential recently developed by the HAL QCD collaboration. To regularize this potential, we remove the deeply attractive core at 𝑟 < 0.4 fm and parametrize the long-range component (𝑟 > 0.4 fm) using a two-range Gaussian form. This procedure preserves the relevant two-body bound-state energy while having a negligible effect on the Ω 3⁢𝑠⁢ 𝑁⁢𝑁 and Ω 3⁢𝑠⁢ Ω 3⁢𝑠 ⁢𝑁 systems. An effective Ω 3⁢𝑠 ⁢𝛼 potential is then constructed by fitting a two-range Gaussian function to the long-range component of the folding potential, enabling calculations of the bound-state energies of the Ω 3⁢𝑠⁢ 𝛼, Ω 3⁢𝑠⁢ 𝛼⁢𝛼, and Ω 3⁢𝑠 ⁢Ω 3⁢𝑠 ⁢𝛼 systems. The regularization procedure leads to a substantial reduction in bound-state energies compared to those obtained with the original potential. We further extend the analysis to Ω 3⁢𝑐 -cluster systems by introducing an Ω 3⁢𝑐 ⁢𝑁 interaction, derived by comparing the existing Ω 3⁢𝑠⁢ Ω 3⁢𝑠 and Ω 3⁢𝑐 ⁢Ω 3⁢𝑐 potentials. Our results suggest that several parametrizations predict bound states in Ω 3⁢𝑐 -containing clusters. Finally, the Ω 3⁢𝑠 ⁢Ω 3⁢𝑠 interaction is described using a contactlike potential approach, motivated by the effective field theory.

binding energy & masses↗
DOE OSTI2023

Free energy and specific heat near a quantum critical point of a metal

We analyze the free energy and the specific heat for fermions interacting with a gapless boson at a quantum-critical point (QCP) in a metal. We use the Luttinger-Ward-Eliashberg formula for the free energy in the normal state, which includes contributions from bosons, fermions, and their interaction, all expressed via fully dressed fermionic and bosonic propagators. The sum of the last two contributions is the free energy F γ of an effective low-energy model of fermions with boson-mediated dynamical 4-fermion interaction V(Ω m ) ∝ 1/|Ω m | γ (the γ model). This purely electronic model has been used to analyze the interplay between non-Fermi liquid behavior and pairing near a QCP, which are both independent of the upper energy cutoff Λ. However, the specific heat C γ (T), obtained from F γ , does depend on Λ. We argue that this dependence is spurious and cancels out, once we include the contribution from bosons. We further argue that the full C(T) is the sum of the contribution from free fermions and the one from a critical boson, with the fully dressed propagator, other terms cancel out. We compare the full C(T) with the C γ (T), obtained using recently proposed regularization of F γ . We argue that for γ < 1, the full C(T) and the regularized C γ (T) differ by a γ-dependent prefactor, while for γ > 1, the full C(T) and C γ (T) differ by the positive contribution from free massless fermions (a positive constant for the electron-phonon case γ = 2). For these γ, C γ (T) is negative, but the full C(T) is positive. We argue that only the full C(T) matters as the positive and the negative contributions originate from the term in C(T) which contains the fully dressed bosonic propagator. Here, we then argue that the normal state remains stable until the pairing instability develops.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Optically induced umklapp shift currents in striped cuprates

Motivated by recent experiments that observed low-frequency second-order optical responses in doped striped superconductors, here we investigate the nonlinear electrodynamics of systems exhibiting a charge density wave (CDW) order parameter. Due to the Bragg scattering off the CDW order, an incoming spatially homogeneous electric field in addition to zero momentum current generates umklapp currents that are modulated in space at momenta of the reciprocal CDW lattice. In particular, here we predict and microscopically evaluate the umklapp shift current, a finite momentum analog of the regular shift current which represents the second-order optical process that downconverts homogeneous AC electric field into low-frequency, zero momentum current. Specifically, we evaluate real-time response functions within mean-field theory via the Keldysh technique and use the Peierls substitution to compute observables at finite momenta in lattice models. We find that systems with certain lattice symmetries (such as inversion symmetry), where the regular shift current is disallowed, may give rise to the umklapp one. We apply our framework to investigate lattice symmetries in layered materials with helical-like stripes and show that both types of shift currents provide insight into the nature of intertwined phases of matter. Finally, we discuss the relation of our findings to recent experiments in striped superconductors. Published by the American Physical Society 2024

Materials Science↗
DOE OSTI2021

Effective field theory for dilute Fermi systems at fourth order

We discuss high-order calculations in perturbative effective field theory for fermions at low energy scales. The Fermi-momentum or k F a s as expansion for the ground-state energy of the dilute Fermi gas is calculated to fourth order, both in cutoff regularization and in dimensional regularization. For the case of spin one-half fermions we find from a Bayesian analysis that the expansion is well converged at this order for |k F a s |≲ 0.5. Furthermore, we show that Padé-Borel resummations can improve the convergence for |k F a s |≲ 1. Here, our results provide important constraints for nonperturbative calculations of ultracold atoms and dilute neutron matter.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Properties of giant dipole resonances within an extended pairing model with a focus on spectral statistics

In this paper, we report on a study of the spectral features associated with dipole resonances in medium mass nuclei (E ≤ 12 MeV), as revealed in the framework of the spd-interacting boson model. The effect of pairing correlations on the theory follows from solutions obtained through an application of the Bethe Ansatz Equation. In general, calculated spectra around the critical point of the vibrational to γ-soft transitions appears to approach that of a Gaussian Orthogonal Ensemble, while near the rotational and vibrational limits of the theory the spectra show more regular behavior. Specifically, the results reveal that the statistical features of the spectra are sensitive to the vector boson pairing strength, c p , in the transition region; that is, when c p is zero, or when the system approaches one of its dynamical symmetries limits the spectrum display regular features, while for stronger c p values, or when near to the critical phase transition region, the spectral feature show more chaotic behavior. Overall, our results indicate that the statistical features are governed by the interplay between dipole resonant energies, pairing correlations, and interactions between and among the single and vector bosons modes of the theory. As part of this work we also found out that chaoticity occurs when results were fit to a Berry-Robnik distribution. Throughout our analyses, we used experimentally known information about both positive and negative parity states. Our findings suggest that dipole resonances appear to be best-described by Poisson statistics for A ≈ 32-138 nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Calculations of the $np→dγ$ reaction in chiral effective field theory

Here, we present a calculation of the radiative capture cross section of the $np→dγ$ reaction in the low-energy range, where the M⁢1 reaction channel dominates. Employing the LENPIC nucleon-nucleon interaction up to the fifth order (N4LO) that is regularized by the semilocal coordinate space regulators, we obtain the initial and final state wave functions, and evaluate the phase shifts of the scattering state and deuteron properties. We derive the transition operator from the chiral effective field theory up to the next-to-next-to leading order (N2LO), where we also regularize the transition operator using regulators consistent with those of the interactions. We compute the capture cross sections, and the results show a converging pattern with the chiral-order expansion of the nucleon-nucleon interaction, where the regulator dependence of the results is weak when higher-order nucleon-nucleon interactions are employed. We quantify the uncertainties of the cross-section results due to the chiral-order truncation. The chirally complete and consistent cross-section results are calculated up to N2LO and they compare well with the experiments and other theoretical predictions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Magnetic moments of A = 3 nuclei obtained from chiral effective field theory operators

Chiral effective field theory (χEFT) provides a framework for obtaining internucleon interactions in a systematically improvable fashion from first principles, while also providing for the derivation of consistent electroweak current operators. In this study, we apply consistently derived interactions and currents towards calculating the magnetic dipole moments of the A = 3 systems 3 H and 3 He. We focus here on LENPIC interactions obtained using semilocal coordinate-space (SCS) regularization. Starting from the momentum-space representation of the LENPIC χEFT vector current, we derive the SCS-regularized magnetic dipole operator up through next-to-next-to-leading order (N 2 LO). We then carry out no-core shell-model calculations for 3 H and 3 He systems using the SCS LENPIC interaction at N 2 LO in χEFT and evaluate the magnetic dipole moments obtained using the consistently derived one-nucleon and two-nucleon electromagnetic currents. As anticipated by prior results with χEFT currents, the current corrections through N 2 LO provide improved, but not yet complete, agreement with experiment for the 3 H and 3 He magnetic dipole moments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Octet and decuplet baryon σ terms and mass decompositions

We present a comprehensive analysis of the SU(3) octet and decuplet baryon masses and σ terms using high-precision lattice QCD data and chiral SU(3) effective theory with finite range regularization. The effects of various systematic uncertainties, including from the scale setting of the lattice data and the regularization prescriptions, are quantified. We find the pion-nucleon and strange nucleon σ terms to be σ πN = 44(3)(3) MeV and σ Ns = 50(6)(1) MeV, respectively. Furthermore, the results provide constraints on the energy-momentum tensor mass decompositions of the SU(3) octet and decuplet baryons, where we find that the trace anomaly and quark-gluon energies decrease for strange baryons due to their larger strange σ terms.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Horocycle regulator: Exact cutoff-independence in AdS/CFT

While the entanglement entropy of a single subregion in quantum field theory is formally infinite and requires regularization, certain combinations of entropies are perfectly finite in the limit that the regulator is removed, the mutual information being a common example. For generic regulator schemes, such as a holographic calculation with a uniform radial cutoff, these quantities show nontrivial dependence on the regulator at finite values of the cutoff. We investigate a holographic regularization scheme defined in three-dimensional anti-de Sitter space constructed from , curves in two-dimensional hyperbolic space perpendicular to all geodesics approaching a single point on the boundary, that leads to finite information measures that are cutoff independent, even at finite values of the regulator. We describe a broad class of such information measures, and describe how the field theory dual to the horocycle regulator is inherently nonlocal. Published by the American Physical Society 2024

Agrawal, Sristy↗
DOE OSTI2020

Two-way coupled particle-turbulence interaction: Effect of numerics and resolution on fluid and particle statistics

Euler-Lagrange point-particle simulation has emerged as a premier methodology for studying dispersed particle-laden flows. This method's popularity stems from its ability to resolve fine-scale fluid structures while also tracking individual particles at reduced cost using an appropriate particle acceleration model. However, the point-particle model has known convergence issues in that refinement of the fluid grid can lead to changes in the predicted statistics. The reasons for non-convergence are twofold: the point-particle two-way coupling force in the Navier-Stokes equations requires a numerical regularization and, without careful implementation, yields a singular force on the fluid with grid refinement. The second factor that yields grid-dependent statistics is that the point-particle force model in general depends on the undisturbed fluid velocity. When the undisturbed fluid velocity is not robustly modeled in a grid-insensitive way, the calculated force for both particles and fluid will be grid-dependent, contaminating their respective statistics. While the first issue regarding regularizing the point-particle source term has received attention in the literature, the consequences of robustly modeling the undisturbed velocity in the context of grid refinement of turbulence has received little attention. Here, we consider decaying homogeneous isotropic turbulence laden with particles at different Stokes numbers. For a given Stokes number, we systematically refine the grid and demonstrate that explicitly modeling the undisturbed fluid velocity yields relative grid insensitivity for the energy of the particle and fluid phases, as well as acceleration of the particles. We also demonstrate that an appropriately defined dissipation rate is also grid-insensitive when an undisturbed fluid velocity correction is used. In contrast, when the undisturbed fluid velocity is modeled using the conventional approach of interpolating the local fluid velocity to the particle location, we show this procedure yields divergent statistics with grid refinement. In particular, we show that higher-order interpolation of the fluid velocity in two-way coupled problems is worse than lower-order interpolation, in the absence of a correction procedure to estimate the undisturbed fluid velocity. We also examine velocity derivative statistics of the fluid phase and demonstrate that these statistics are not in general convergent even when the undisturbed fluid velocity is explicitly modeled. Collectively, the observations in this paper are used to present a philosophy on the types of questions which are answerable with point-particle methods.

42 ENGINEERING↗
DOE OSTI2021

Toward Complete Leading-Order Predictions for Neutrinoless Double β Decay

The amplitude for the neutrinoless double β (0νββ) decay of the two-neutron system nn → ppe – e – constitutes a key building block for nuclear-structure calculations of heavy nuclei employed in large-scale 0νββ searches. Assuming that the 0νββ process is mediated by a light-Majorana-neutrino exchange, a systematic analysis in chiral effective field theory shows that already at leading order a contact operator is required to ensure renormalizability. Here, we develop a method to estimate the numerical value of its coefficient (in analogy to the Cottingham formula for electromagnetic contributions to hadron masses) and validate the result by reproducing the charge-independence-breaking contribution to the nucleon-nucleon scattering lengths. Our central result, while derived in dimensional regularization, is given in terms of the renormalized amplitude Aν(|p|, |p'|) , matching to which will allow one to determine the contact-term contribution in regularization schemes employed in nuclear-structure calculations. Our results thus greatly reduce a crucial uncertainty in the interpretation of searches for 0νββ decay.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗