Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Renormalization and regularization”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Dynamics of a single polyampholyte chain

Polymers that feature both positive and negative charges along chains, known as polyampholytes, represent a class of materials that hold promise for a new generation of energy storage devices, the design of which will require knowledge of the underlying structure and dynamics. Here, we develop a theory based on the Rouse model for the dynamic structure factor of a single polyampholyte chain in the weak coupling regime (negligible intramolecular electrostatics) or subjected to weak external electric fields (governed by linear response). Neglecting effects of small ions, we find deviations in scaling from the classic Rouse theory and make predictions for scattering experiments performed on polyampholytes. We find that, under weak coupling with arbitrarily strong fields, the dynamics are highly dependent on the charge distribution and consequently look at two representative examples—random charge densities and periodic charge densities—with different scaling properties. Additionally, under weak fields, the dynamics are largely independent of charge distribution. Finally, we investigate the influence of hydrodynamic effects and the implications of including inertial effects in the model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Topological field theory with Haagerup symmetry

Here, we construct a (1 + 1)d topological field theory (TFT) whose topological defect lines (TDLs) realize the transparent Haagerup $\mathscr{H}$ 3 fusion category. This TFT has six vacua, and each of the three non-invertible simple TDLs hosts three defect operators, giving rise to a total of 15 point-like operators. The TFT data, including the three-point coefficients and lasso diagrams, are determined by solving all the sphere four-point crossing equations and torus one-point modular invariance equations. We further verify that the Cardy states furnish a non-negative integer matrix representation under TDL fusion. While many of the constraints we derive are not limited to this particular TFT with six vacua, we leave open the construction of TFTs with two or four vacua. Finally, TFTs realizing the Haagerup $\mathscr{H}$ 1 and $\mathscr{H}$ 2 fusion categories can be obtained by gauging algebra objects. This article makes a modest offering in our pursuit of exotica and the quest for their eventual conformity.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Sub-system self-consistency in coupled cluster theory

Here, in this article, we provide numerical evidence indicating that the single-reference coupled-cluster (CC) energies can be calculated alternatively to their copybook definition. We demonstrate that the CC energy can be reconstructed by diagonalizing the effective Hamiltonians describing correlated sub-systems of the many-body system. In the extreme case, we provide numerical evidence that the CC energy can be reproduced through the diagonalization of the effective Hamiltonian describing sub-system composed of a single electron. These properties of the CC formalism can be exploited to design protocols to define effective interactions in sub-systems used as probes to calculate the energy of the entire system and introduce a new type of self-consistency for approximate CC approaches.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Modification of Lie's transform perturbation theory for charged particle motion in a magnetic field

It is pointed out that the conventional Lie transform perturbation theory for the guiding center motion of charged particles in a magnetic field needs to be modified for ordering inconsistency. There are two reasons. First, the ordering difference between the temporal variation of gyrophase and that of the other phase space coordinates needs to be taken into account. Second, it is also important to note that the parametric limit of the derivative of a function is not equivalent to the derivative of the limit function. When these facts are taken into account, the near identity transformation rule for one form related to the Lagrangian is modified. With the modified near identity transformation rule, the drift motion of charged particles can be described in the first order, instead of the second order and beyond through a tedious expansion process as in the conventional formulation. Here, this resolves the discrepancy between the direct and Lie transform treatments in the Lagrangian perturbation theory for charged particle motion in a magnetic field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Exploring the exact limits of the real-time equation-of-motion coupled cluster cumulant Green’s functions

In this paper, we analyze the properties of the recently proposed real-time equation-of-motion coupled-cluster (RT-EOM-CC) cumulant Green’s function approach [Rehr et al., J. Chem. Phys. 152, 174113 (2020)]. We specifically focus on identifying the limitations of the original time-dependent coupled cluster (TDCC) ansatz and propose an enhanced double TDCC ansatz, ensuring the exactness in the expansion limit. In addition, we introduce a practical cluster-analysis-based approach for characterizing the peaks in the computed spectral function from the RT-EOM-CC cumulant Green’s function approach, which is particularly useful for the assignments of satellite peaks when many-body effects dominate the spectra. Our preliminary numerical tests focus on reproducing, approximating, and characterizing the exact impurity Green’s function of the three-site and four-site single impurity Anderson models using the RT-EOM-CC cumulant Green’s function approach. The numerical tests allow us to have a direct comparison between the RT-EOM-CC cumulant Green’s function approach and other Green’s function approaches in the numerical exact limit.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Renormalization of states and quasiparticles in many-body downfolding

We explore the principles of many-body Hamiltonian complexity reduction via downfolding on an effective low-dimensional representation. We show that the renormalization factor provides a unique measure of the quality of the compression as it directly represents the projection between the approximate stationary state of the many-body Hamiltonian and the full many-body wavefunction. Hence, the renormalization factor is a measure of fidelity between the effective (reduced-rank) description and the full many-body treatment for arbitrary (i.e., ground and excited) states. When the entire problem is mapped on a system of interacting quasiparticles [Romanova et al., npj Comput. Mater. 9, 126 (2023)], the effective Hamiltonians can faithfully reproduce the physics only when a clear energy scale separation exists between the subsystems and their environment. We also demonstrate that it is necessary to include quasiparticle renormalization at distinct energy scales, capturing the distinct interaction between subsystems and their surrounding environments. Numerical results from simple, exactly solvable models highlight the limitations and strengths of this approach, particularly for ground and low-lying excited states. This work lays the groundwork for applying dynamical downfolding techniques to problems concerned with (quantum) interfaces.

Green-functions technique↗
DOE OSTI2026

Lattice calculation of short-range contributions to neutrinoless double-beta decay 𝜋 − → 𝜋 + ⁢𝑒⁢𝑒 at physical pion mass

Neutrinoless double-beta (0⁢𝜈⁢𝛽⁢𝛽) decays provide an excellent probe for determining whether neutrinos are Dirac or Majorana fermions. The short-range matrix elements associated with the 𝜋 − → 𝜋 + ⁢𝑒⁢𝑒 process contribute at leading order in the 0⁢𝜈⁢𝛽⁢𝛽 decay channel 𝑛⁢𝑛 → 𝑝⁢𝑝⁢𝑒⁢𝑒 through pion exchange between nucleons. However, current lattice calculations show notable discrepancies in predicting these short-range contributions. To address this issue, we perform a lattice QCD calculation of the 𝜋 − → 𝜋 + ⁢𝑒⁢𝑒 matrix elements using domain wall fermion ensembles at the physical pion mass generated by the RBC/UKQCD Collaboration. To mitigate contamination from around-the-world effects, we develop a new method to reconstruct and subtract them directly from lattice data. We then perform nonperturbative renormalization in the regularization-independent symmetric momentum-subtraction scheme (RI/SMOM), using the (𝛾 𝜇 , 𝛾 𝜇 ) and ($\not{𝑞}$, $\not{𝑞}$) projectors. Compared with previous studies, this work reduces the uncertainties in the matrix elements and provides an independent cross-check that helps to reconcile the discrepancies among previous lattice calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Regularization scheme dependence of the counterterms in the galaxy bias expansion

In this paper we explore how different regularization prescriptions affect the counterterms in the renormalization of the galaxy bias expansion. We work in the context of primordial local non-Gaussianity including non-linear gravitational evolution. We carry out the one-loop renormalization of the field $δ$$^{ρ}_{2}$ (i.e. the square of the matter overdensity field) up to third order in gravitational evolution. Three regularization schemes are considered and their impact on the values of the counterterms is studied. We explicitly verify that the coefficients of the non-boost invariant operators are regularization scheme independent.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Second-order renormalized Hamiltonian of Yukawa theory

Using the renormalization group procedure for effective particles we calculate the effective Hamiltonians in the theory of a fermion field coupled to a scalar field via the Yukawa interaction. The theory is renormalized by the addition of counterterms. Necessary counterterms are determined by computing matrix elements of the effective Hamiltonian. All calculations are performed up to the second order in the expansion in powers of the coupling constant. Renormalized effective Hamiltonians are well-defined symmetric forms acting in the Fock space as opposed to the renormalized bare Hamiltonian, which is not well defined without regularization. We introduce computational techniques that should streamline higher-order calculations and may be of independent interest.

Ab initio calculations↗
DOE OSTI2026

Quantum critical collapse abhors a naked singularity

Classical critical collapse yields naked singularities from smooth initial data, challenging cosmic censorship, and shaping the spectrum of primordial black holes. We show that one-loop vacuum polarization near the threshold qualitatively changes this outcome by dressing the singularity with a horizon within a controlled semiclassical regime. In analytically tractable Einstein-scalar critical spacetimes, a one-loop $s$-wave treatment linearized around self-similar backgrounds shows that regularity uniquely selects an asymptotically Minkowskian, vacuum-polarization state. Its renormalized stress tensor carries a universal quantum growing mode that competes with the classical unstable mode, shifts the critical point, and generates a trapped surface along with a finite mass gap at the new threshold, thereby enforcing horizon formation even under arbitrary fine-tuning. In primordial collapse, the threshold shift enters exponentially into the formation fraction, while the mass gap truncates the low-mass tail, suggesting potentially important consequences for the predicted mass spectrum. Furthermore, these results provide a self-consistent semiclassical treatment of critical collapse and yield sharp predictions within the one-loop, near-critical, linearized regime.

Anomalies↗
DOE OSTI2023

An in-silico NMR laboratory for nuclear magnetic shieldings computed via finite fields: Exploring nucleus-specific renormalizations of MP2 and MP3

Here we developed and implemented a method-independent, fully numerical, finite difference approach to calculating nuclear magnetic resonance shieldings, using gauge-including atomic orbitals. The resulting capability can be used to explore non-standard methods, given only the energy as a function of finite-applied magnetic fields and nuclear spins. For example, standard second-order Møller-Plesset theory (MP2) has well-known efficacy for 1H and 13C shieldings and known limitations for other nuclei such as 15N and 17O. It is, therefore, interesting to seek methods that offer good accuracy for 15N and 17O shieldings without greatly increased compute costs, as well as exploring whether such methods can further improve 1H and 13C shieldings. Using a small molecule test set of 28 species, we assessed two alternatives: κ regularized MP2 (κ-MP2), which provides energy-dependent damping of large amplitudes, and MP2.X, which includes a variable fraction, X, of third-order correlation (MP3). The aug-cc-pVTZ basis was used, and coupled cluster with singles and doubles and perturbative triples [CCSD(T)] results were taken as reference values. Our κ-MP2 results reveal significant improvements over MP2 for 13C and 15N, with the optimal κ value being element-specific. κ-MP2 with κ = 2 offers a 30% rms error reduction over MP2. For 15N, κ-MP2 with κ = 1.1 provides a 90% error reduction vs MP2 and a 60% error reduction vs CCSD. On the other hand, MP2.X with a scaling factor of 0.6 outperformed CCSD for all heavy nuclei. These results can be understood as providing renormalization of doubles amplitudes to partially account for neglected triple and higher substitutions and offer promising opportunities for future applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Flavor diagonal nucleon charges using clover fermions on MILC HISQ ensembles

We present lattice results for the flavor diagonal charges of the proton from the analysis of eight ensembles generated using 2+1+1-flavors of highly improved staggered quarks by the MILC Collaboration. The calculation includes all the needed connected and disconnected contributions to nucleon three-point function. For extracting matrix elements using fits to the spectral decomposition of these correlation functions, two strategies to remove excited state contributions are employed and compared. To renormalize these charges, the 2+1-flavor mixing matrix is calculated in the regularization independent symmetric momentum subtraction intermediate scheme on the lattice. The final results are presented in the $\overline{MS}$ scheme at scale 2 GeV. The axial charges for the proton are $𝑔$$^{𝑢}_{𝐴}$ = 0.781⁢(25), $𝑔$$^{𝑑}_{𝐴}$ =−0.440⁢(39), and $𝑔$$^{𝑠}_{𝐴}$ = −0.055⁢(9); the tensor charges are $𝑔$$^{𝑢}_{𝑇}$ = 0.782⁢(28), $𝑔$$^{𝑑}_{𝑇}$ = −0.195⁢(16), and $𝑔$$^{𝑠}_{𝑇}$ = −0.0016⁢(12); and the scalar charges are $𝑔$$^{𝑢}_{𝑆}$ = 9.39⁢(88), $𝑔$$^{𝑑}_{𝑆}$ = 8.84⁢(93), and $𝑔$$^{𝑠}_{𝑆}$ = 0.37⁢(14). Results for the neutron are given by the 𝑢 ↔ 𝑑 interchange. Results for the sigma terms are 𝜎 𝜋⁢𝑁 | standard = 42⁢(6) MeV from a “standard” analysis and 𝜎 𝜋⁢𝑁 | 𝑁⁢𝜋 = 61⁢(6) MeV from an “𝑁⁢𝜋” analysis that includes the contributions of multihadron 𝑁⁢𝜋 excited states as motivated by chiral perturbation theory. Our preferred value 𝜎 𝜋⁢𝑁 | 𝑁⁢𝜋 is consistent with the phenomenological extraction from 𝜋 −𝑁 scattering data. The strangeness content of the proton, for which the standard analysis is appropriate, is 𝜎 𝑠 | standard = 35⁢(13) MeV.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1971

Fluctuation induced conductivity of superconductors above the transition temperature: Regularization of the Maki diagram

The Maki contribution to the conductivity above the superconducting transition temperature is regularized within the framework of the BCS theory. This is achieved through the renormalization of the impurity scattering vertex by inclusion of the effects of pair fluctuations. The conductivity is evaluated for a thin film. It depends only on the reduced temperature and the normal resistance per square. Fair agreement is found with Al films over a wide temperature range. Agreement is not found with experiments on Bi, Pb and Ga films, which apparently contain a strong additional pair-breaking effect. The temperature range in which interactions among fluctuations become important in the Maki conductivity is generally larger than that given by the Ginzburg criterion.

Keller, J.↗
NASA NTRS1972

Fluctuation-induced conductivity of superconductors above the transition temperature - Regularization of the Maki diagram.

The Maki contribution to the conductivity above the superconducting transition temperature is regularized within the framework of the BCS theory. This is achieved through the renormalization of the impurity-scattering vertex by inclusion of the effects of pair fluctuations. The conductivity is evaluated for a thin film. It depends only on the reduced temperature and the normal resistance per square. Fair agreement is found with Al films over a wide temperature range. Agreement is not found with experiments on Bi, Pb, and Ga films, which apparently contain a strong additional pair-breaking effect. The temperature range in which interactions among fluctuations become important in the Maki conductivity is generally larger than that given by the Ginzburg criterion.

Keller, J.↗
DOE OSTI2025

Qubit Regularization of Quantum Field Theories

To study quantum field theories on a quantum computer, we must begin with Hamiltonians defined on a finite-dimensional Hilbert space and then take appropriate limits. This approach can be seen as a new type of regularization for quantum field theories, which we refer to as qubit regularization. A related finite-dimensional regularization, known as the D-theory approach, was proposed long ago as a general framework for all quantum field theories. In this framework, the dimensionality of the local Hilbert space at each spatial point can increase as needed through an additional flavor index. To reproduce asymptotically free QFTs, most studies assume that qubit-regularized theories require extending the local Hilbert space to infinity. However, contrary to this common belief, recent discoveries in (1+1) dimensions have revealed two examples where asymptotic freedom appears to emerge within a strictly finite-dimensional local Hilbert space through a novel renormalization group (RG) flow. These findings motivate further investigation into whether asymptotically free gauge theories could also emerge within a strictly finite-dimensional local Hilbert space. To support these explorations, we propose an orthonormal basis called the monomer-dimer-tensor-network (MDTN) basis and use it to construct new types of qubit-regularized lattice gauge theories.

Chandrasekharan, Shailesh [Duke Univ., Durham, NC ↗
DOE OSTI2024

Pion valence quark distribution at physical pion mass of N f = 2 + 1 + 1 lattice QCD

Abstract We present a state-of-the-art calculation of the unpolarized pion valence-quark distribution in the framework of large-momentum effective theory (LaMET) with improved handling of systematic errors as well as two-loop perturbative matching. We use lattice ensembles generated by the MILC collaboration at lattice spacinga≈ 0.09 fm, lattice volume 64 3 × 96,N f = 2 + 1 + 1 flavors of highly-improved staggered quarks and a physical pion mass. The LaMET matrix elements are calculated with pions boosted to momentumP z ≈ 1.72 GeV with high-statistics ofO(10 6 ) measurements. We study the pion PDF in both hybrid-ratio and hybrid-regularization-independent momentum subtraction (hybrid-RI/MOM) schemes and also compare the systematic errors with and without the addition of leading-renormalon resummation (LRR) and renormalization-group resummation (RGR) in both the renormalization and lightcone matching. The final lightcone PDF results are presented in the modified minimal-subtraction scheme at renormalization scaleμ= 2.0 GeV. We show that thex-dependent PDFs are compatible between the hybrid-ratio and hybrid-RI/MOM renormalization with the same improvements. We also show that systematics are greatly reduced by the simultaneous inclusion of RGR and LRR and that these methods are necessary if improved precision is to be reached with higher-order terms in renormalization and matching.

Physics↗
DOE OSTI2022

On the closedness and geometry of tensor network state sets

Tensor network states (TNS) are a powerful approach for the study of strongly correlated quantum matter. The curse of dimensionality is addressed by parametrizing the many-body state in terms of a network of partially contracted tensors. These tensors form a substantially reduced set of effective degrees of freedom. In practical algorithms, functionals like energy expectation values or overlaps are optimized over certain sets of TNS. Concerning algorithmic stability, it is important whether the considered sets are closed because, otherwise, the algorithms may approach a boundary point that is outside the TNS set and tensor elements diverge. Here we discuss the closedness and geometries of TNS sets, and we propose regularizations for optimization problems on non-closed TNS sets. We show that sets of matrix product states (MPS) with open boundary conditions, tree tensor network states, and the multiscale entanglement renormalization ansatz are always closed, whereas sets of translation-invariant MPS with periodic boundary conditions (PBC), heterogeneous MPS with PBC, and projected entangled pair states are generally not closed. The latter is done using explicit examples like the W state, states that we call two-domain states, and fine-grained versions thereof.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Black hole horizon edge partition functions

We extend a formula for 1-loop black hole determinants by Denef, Hartnoll, and Sachdev (DHS) to spinning fields on any (d + 1)-dimensional static spherically symmetric black hole. By carefully analyzing the regularity condition imposed on the Euclidean eigenfunctions, we reveal an unambiguous bulk-edge split in the 1-loop Euclidean partition function for tensor fields of arbitrary integer spin: the bulk part captures the “renormalized” thermal canonical partition function recently discussed in [1]; the edge part is related to quasinormal modes (QNMs) that fail to analytically continue to a subset of Euclidean modes with enhanced fall-offs near the origin. Since the edge part takes the form of a path integral on S d–1 , this suggests that these are associated with degrees of freedom living on the bifurcation surface in the Lorentzian two-sided black hole geometry. For massive higher spin on static BTZ and massive vector on Nariai black holes, we find that the edge partition function is related to the QNMs with lowest overtone numbers.

79 ASTRONOMY AND ASTROPHYSICS↗