Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical renormalization group”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Disentangling observable dependence in SCET I and SCET II anomalous dimensions: angularities at two loops

The resummation of radiative corrections to collider jet observables using soft collinear effective theory is encoded in differential renormalization group equations (RGEs), with anomalous dimensions depending on the observable under consideration. This observable dependence arises from the ultraviolet (UV) singular structure of real phase space integrals in the effective field theory. We show that the observable dependence of anomalous dimensions in SCET I problems can be disentangled by introducing a suitable UV regulator in real radiation integrals. Resummation in the presence of the new regulator can be performed by solving a two-dimensional system of RGEs in the collinear and soft sectors, and resembles many features of resummation in SCET II theories by means of the rapidity renormalization group. We study the properties of SCET I with the additional regulator and explore the connection with the system of RGEs in SCET II theories, highlighting some universal patterns that can be exploited in perturbative calculations. As an application, we compute the two-loop soft and jet anomalous dimensions for a family of recoil-free angularities and give new analytic results. This allows us to study the relations between the SCET I and SCET II limits for these observables. We also discuss how the extra UV regulator can be exploited to calculate anomalous dimensions numerically, and the prospects for numerical resummation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Dihadron correlations in small- x DIS at NLO: transverse momentum dependent fragmentation

We compute the inclusive dihadron cross-section in Deep Inelastic Scattering at next-to-leading order (NLO) and small x in the Color Glass Condensate. We focus on the kinematic limit where the hadrons are produced at forward rapidities (in the direction of the virtual photon) and back-to-back in the transverse plane. Our calculation demonstrates that the coefficient of the Sudakov double logarithm for this process is –$\frac{αs}{2π}$ [C F + $\frac{Nc}{2}$] instead of –$\frac{αsNc}{4π}$ when back-to-back jets are measured in the final state. To preserve the universality of the Sudakov soft factor associated with the Weizsäcker-Williams transverse momentum dependent (TMD) gluon distribution, we promote the collinear fragmentation functions into TMD fragmentation functions. We then perform the resummation of the Sudakov logarithms through Collins-Soper-Sterman evolution of the TMD fragmentation functions and the Weizsäcker-Williams TMD gluon distribution. Finally, analytic expressions are obtained for the NLO coefficient functions in the $\overline{MS}$-scheme. These results pave the way towards numerically calculating dihadron correlations at small x at the future Electron-Ion Collider with full NLO accuracy.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Time-dependent approach to inelastic scattering spectroscopies in and away from equilibrium: Beyond perturbation theory

We propose a nonperturbative numerical approach to calculate the spectrum of a many-body Hamiltonian with time and momentum resolution by exactly recreating a scattering event using the time-dependent Schrödinger equation. Akin to an actual inelastic scattering experiment, we explicitly account for the incident and scattered particles (e.g., photons, neutrons, electrons, etc.) in the Hamiltonian and obtain the spectrum by measuring the energy and momentum lost by the particle after interacting with the sample. We illustrate the method by calculating the spin excitations of a Mott-insulating Hubbard chain after a sudden quench with the aid of the time-dependent density matrix renormalization group method. Our formalism can be applied to different forms of spectroscopies, such as neutron and Compton scattering and electron-energy-loss spectroscopy, as an example.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Entanglement Renormalization for Quantum Field Theories with Discrete Wavelet Transforms

We propose an adaptation of Entanglement Renormalization for quantum field theories that, through the use of discrete wavelet transforms, strongly parallels the tensor network architecture of the Multiscale Entanglement Renormalization Ansatz (a.k.a. MERA). Our approach, called wMERA, has several advantages of over previous attempts to adapt MERA to continuum systems. In particular, (i) wMERA is formulated directly in position space, hence preserving the quasi-locality and sparsity of entanglers; and (ii) it enables a built-in RG flow in the implementation of real-time evolution and in computations of correlation functions, which is key for efficient numerical implementations. As examples, we describe in detail two concrete implementations of our wMERA algorithm for free scalar and fermionic theories in (1+1) spacetime dimensions. Possible avenues for constructing wMERAs for interacting field theories are also discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

One-dimensional Holstein model revisited

Here, we analyze the global ground-state (quantum) phase diagram of the one-dimensional spinful Holstein model at half filling as a function of the strength of the electron-phonon coupling (represented by the strength of the phonon-induced attraction, U) and the phonon frequency ω 0 . In addition to reanalyzing the various asymptotic regimes, we carry out density-matrix renormalization group simulations to correct previous inferences concerning the antiadiabatic (large ω 0 ) and strong-coupling (large U) regimes. There are two distinct phases—a fully gapped commensurate charge-density-wave phase and a spin-gapped Luther-Emery phase with a gapless charge mode—separated by a phase boundary, with a shape that reflects different microscopic physics in the weak- and strong-coupling limits.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Spinon continuum in the Heisenberg quantum chain compound Sr 2 ⁢V 3 ⁢O 9

Magnetic excitations in the spin chain candidate Sr 2 ⁢V 3 ⁢O 9 have been investigated by inelastic neutron scattering on a single crystal sample. A spinon continuum with a bandwidth of ~22 meV is observed along the chain formed by alternating magnetic V 4+ and nonmagnetic V 5+ ions, which reveals the importance of the orbital degree of freedom in determining the chain axis as identified by prior electronic structure calculations. Incipient magnetic Bragg peaks due to weak ferromagnetic interchain couplings emerge when approaching the magnetic transition at T N ~ 5.3 K, while the excitations remain gapless within the instrumental resolution. Comparisons to the Bethe ansatz, density matrix renormalization group calculations, and effective field theories confirm Sr 2 ⁢V 3 ⁢O 9 as a host of weakly coupled S = 1/2 chains dominated by antiferromagnetic intrachain interactions of ~7.1(1) meV.

1-dimensional spin chains↗
DOE OSTI2022

Perturbative calculations of gravitational form factors at large momentum transfer

We perform a perturbative QCD analysis of the gravitational form factors (GFFs) of nucleon at large momentum transfer. We derive the explicit factorization formula of the GFFs in terms of twist-3 and twist-4 light-cone distribution amplitudes of nucleon. Power behaviors for these GFFs are obtained from the leading order calculations. Numeric results of the quark and gluon contributions to various GFFs are presented with model assumptions for the distribution amplitudes in the literature. We also present the perturbative calculations of the scalar form factor $\langle P$' $|F^2|P \rangle$ for pion and proton at large momentum transfer.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

A multilayer multi-configurational approach to efficiently simulate large-scale circuit-based quantum computers on classical machines

Here, a multilayer multi-configurational theory framework is adapted to simulate circuit-based quantum computers. Quantum addition of superpositions of an exponential number of summands is performed in polynomial time with high accuracy. We demonstrate numerically accurate calculations including up to one million qubits for entangling benchmarks. Simulation cost can be assessed by entropy-based entanglement measures. For the considered systems, we show that the entanglement only grows weakly with the system size. The present simulations demonstrate how quantum algorithms in low-entropy regimes can be used efficiently on classically simulated quantum computers.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Renormalization group effects in dark matter interactions

We present a renormalization-group (RG) analysis of dark matter interactions with the standard model, where dark matter is allowed to be a component of an electroweak multiplet, and has a mass at or below the electroweak scale. We consider, in addition to the gauge interactions, the complete set of effective operators for dark matter interactions with the standard model above the weak scale, up to and including mass dimension six. We calculate the RG evolution of these operators from the high scale Λ down to the weak scale, and perform the matching to the tower of effective theories below the weak scale. We also summarize the RG evolution below the weak scale and the matching to the nonrelativistic nuclear interactions. We present several numerical examples and show that in certain cases the dark matter — nucleus scattering rate can change by orders of magnitude when the electroweak running is included.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Next-to-next-to-leading power corrections to unpolarized Semi-Inclusive Deep Inelastic Scattering

Semi-Inclusive Deep Inelastic Scattering (SIDIS) is a key tool for exploring the three-dimensional structure of the nucleon through Transverse Momentum Dependent parton distributions and fragmentation functions. While leading-power contributions to the SIDIS cross-section are well established, next-to-leading power (NLP) corrections of order 1/Q and next-to-next-to-leading power (NNLP) corrections of order 1/Q 2 to the hadronic tensor have only recently begun to be systematically investigated. These corrections are essential for reliable phenomenology and interpretation of modern high-precision data. In recent papers by one of the authors, NNLP corrections to the Drell-Yan process were derived using the rapidity factorization formalism. In the present work, we extend this approach to SIDIS and obtain analytic expressions for the unpolarized structure functions. We derive NNLP corrections that include convolutions of unpolarized distributions, f 1 , with unpolarized fragmentation functions, D 1 , and Boer-Mulders functions, ${h}_1^{\perp }$, with Collins fragmentation functions, ${H}_1^{\perp }$. We compare our results with previous formulations, provide numerical studies, confront our predictions with HERMES and COMPASS measurements, and present predictions for future experiments at Jefferson Lab and the Electron-Ion Collider.

deep inelastic scattering↗
DOE OSTI2024

Efficient simulation of low-temperature physics in one-dimensional gapless systems

Here, we discuss the computational efficiency of the finite-temperature simulation with minimally entangled typical thermal states (METTS). To argue that METTS can be efficiently represented as matrix product states, we present an analytic upper bound for the average entanglement Rényi entropy of METTS for a Rényi index 0 < q ≤ 1. In particular, for one-dimensional (1D) gapless systems described by conformal field theories, the upper bound scales as O⁡(cN 0 ⁢log⁡β) where c is the central charge and N is the system size. Furthermore, we numerically find that the average Rényi entropy exhibits a universal behavior characterized by the central charge and is roughly given by half of the analytic upper bound. Based on these results, we show that METTS can provide a speedup compared to employing the purification method to analyze thermal equilibrium states at low temperatures in 1D gapless systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Phase diagram of a distorted kagome antiferromagnet and application to Y-kapellasite

Abstract We investigate the magnetism of a previously unexplored distorted spin-1/2 kagome model consisting of three symmetry-inequivalent nearest-neighbor antiferromagnetic Heisenberg couplings J ⬡ , J , and $$J^{\prime}$$ J ′ , and uncover a rich ground state phase diagram even at the classical level. Using analytical arguments and numerical techniques we identify a collinear $$\overrightarrow{Q}=0$$ Q → = 0 magnetic phase, two unusual non-collinear coplanar $$\overrightarrow{Q}=(1/3,1/3)$$ Q → = ( 1 / 3 , 1 / 3 ) phases and a classical spin liquid phase with a degenerate manifold of non-coplanar ground states, resembling the jammed spin liquid phase found in the context of a bond-disordered kagome antiferromagnet. We further show with density functional theory calculations that the recently synthesized Y-kapellasite Y 3 Cu 9 (OH) 19 Cl 8 is a realization of this model and predict its ground state to lie in the region of $$\overrightarrow{Q}=(1/3,1/3)$$ Q → = ( 1 / 3 , 1 / 3 ) order, which remains stable even after the inclusion of quantum fluctuation effects within variational Monte Carlo and pseudofermion functional renormalization group. The presented model opens a new direction in the study of kagome antiferromagnets.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Superspin renormalization and slow relaxation in random spin systems

We develop an excited-state real-space renormalization group (RSRG-X) formalism to describe the dynamics of conserved densities in randomly interacting spin-12 systems. Our formalism is suitable for systems with U(1) and Z2 symmetries, and we apply it to chains of randomly positioned spins with dipolar XX+YY interactions, as arise in Rydberg quantum simulators and other platforms. The formalism generates a sequence of effective Hamiltonians that provide approximate descriptions for dynamics on successively smaller energy scales. These effective Hamiltonians involve “superspins”: two-level collective degrees of freedom constructed from (anti)aligned microscopic spins. Conserved densities can then be understood as relaxing via coherent collective spin flips. For the well-studied simpler case of randomly interacting nearest-neighbor XX+YY chains, the superspins reduce to single spins. Our formalism also leads to a numerical method capable of simulating the dynamics up to an otherwise inaccessible combination of large system size and late time. Focusing on disorder-averaged infinite-temperature autocorrelation functions, in particular the spin survival probability Sp¯(t), we demonstrate quantitative agreement between our algorithm and exact diagonalization (ED) at low but nonzero frequencies. Such agreement holds for chains with nearest-neighbor, next-nearest-neighbor, and long-range dipolar interactions. Our results indicate decay of Sp¯(t) slower than any power law and feature no significant deviation from the ∼1/ln2(t) asymptote expected from the infinite-randomness fixed-point of the nearest-neighbor model. We also apply the RSRG-X formalism to two-dimensional long-range systems of moderate size and find slow late-time decay of Sp¯(t).

Zhao, Yi J↗
DOE OSTI2021

Quantum phase transitions in long-range interacting hyperuniform spin chains in a transverse field

Hyperuniform states of matter are characterized by anomalous suppression of long-wavelength density fluctuations. While most of the interesting cases of disordered hyperuniformity are provided by complex many-body systems such as liquids or amorphous solids, classical spin chains with certain long-range interactions have been shown to demonstrate the same phenomenon. Such spin-chain systems are ideal models for exploring the effects of quantum mechanics on hyperuniformity. It is well-known that the transverse field Ising model shows a quantum phase transition (QPT) at zero temperature. Under the quantum effects of a transverse magnetic field, classical hyperuniform spin chains are expected to lose their hyperuniformity. High-precision simulations of these cases are complicated because of the presence of highly nontrivial long-range interactions. We perform an extensive analysis of these systems using density matrix renormalization group simulations to study the possibilities of phase transitions and the mechanism by which they lose hyperuniformity. Even for a spin chain of length 30, we see discontinuous changes in properties like the “τ order metric” of the ground state, the measure of hyperuniformity, and the second cumulant of the total magnetization along the x-direction, all suggestive of first-order QPTs. An interesting feature of the phase transitions in these disordered hyperuniform spin chains is that, depending on the parameter values, the presence of a transverse magnetic field may lead remarkably to an increase in the order of the ground state as measured by the “τ order metric,” even if hyperuniformity is lost. Therefore, it would be possible to design materials to target specific novel quantum behaviors in the presence of a transverse magnetic field. Our numerical investigations suggest that these spin chains can show no more than two QPTs. We further analyze the long-range interacting spin chains via the Jordan-Wigner mapping onto a system of spinless fermions, showing that under the pairwise-interaction approximation and a mean-field treatment, there can be at most two quantum phase transitions. Here, based on these numerical and theoretical explorations, we conjecture that for spin chains with long-range pair interactions that have convergent cosine transforms, there can be a maximum of two quantum phase transitions at zero temperature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Conservation laws in coupled cluster dynamics at finite temperature

We extend the finite-temperature Keldysh non-equilibrium coupled cluster theory (Keldysh-CC) [A. F. White and G. K.-L. Chan, J. Chem. Theory Comput. 15, 6137–6253 (2019)] to include a time-dependent orbital basis. When chosen to minimize the action, such a basis restores local and global conservation laws (Ehrenfest’s theorem) for all one-particle properties while remaining energy conserving for time independent Hamiltonians. We present the time-dependent Keldysh orbital-optimized coupled cluster doubles method in analogy with the formalism for zero-temperature dynamics, extended to finite temperatures through the time-dependent action on the Keldysh contour. To demonstrate the conservation property and understand the numerical performance of the method, we apply it to several problems of nonequilibrium finite-temperature dynamics: a 1D Hubbard model with a time-dependent Peierls phase, laser driving of molecular H2, driven dynamics in warm-dense silicon, and transport in the single impurity Anderson model.

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

Accelerating resonant spectroscopy simulations using multishifted biconjugate gradient

Resonant spectroscopies, which involve intermediate states with finite lifetimes, provide important insights into collective excitations in quantum materials that are otherwise inaccessible. However, theoretical understanding in this area is often limited by the numerical challenges of solving Kramers-Heisenberg-type response functions for large-scale systems. To address this, we introduce a multishifted biconjugate gradient algorithm that exploits the shared structure of Krylov subspaces across spectra with varying incident energies, effectively reducing the computational complexity to that of linear spectroscopies. Both mathematical proofs and numerical benchmarks confirm that this algorithm substantially accelerates spectral simulations, achieving constant complexity independent of the number of incident energies, while ensuring accuracy and stability. This development provides a scalable, versatile framework for simulating advanced spectroscopies in quantum materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

A time-dependent momentum-resolved scattering approach to core-level spectroscopies

While new light sources allow for unprecedented resolution in experiments with X-rays, a theoretical understanding of the scattering cross-section is lacking. In the particular case of strongly correlated electron systems, numerical techniques are quite limited, since conventional approaches rely on calculating a response function (Kramers-Heisenberg formula) that is obtained from a perturbative analysis of scattering processes in the frequency domain. This requires a knowledge of a full set of eigenstates in order to account for all intermediate processes away from equilibrium, limiting the applicability to small tractable systems. In this work, we present an alternative paradigm, recasting the problem in the time domain and explicitly solving the time-dependent Schrödinger equation without the limitations of perturbation theory: a faithful simulation of the scattering processes taking place in actual experiments, including photons and core electrons. We show how this approach can yield the full time and momentum resolved Resonant Inelastic X-Ray Scattering (RIXS) spectrum of strongly interacting many-body systems. We demonstrate the formalism with an application to Mott insulating Hubbard chains using the time-dependent density matrix renormalization group method, which does not require a priory knowledge of the eigenstates and can solve very large systems with dozens of orbitals. This approach can readily be applied to systems out of equilibrium without modification and generalized to other spectroscopies.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Super-leading logarithms in top-quark pair production at hadron colliders

To date, the appearance and resummation of “super-leading” logarithms in hadron-hadron collisions has been studied only for massless parton states. We extend the formalism to include an arbitrary number of massive final states. We derive the corresponding anomalous dimension and identify an additional Coulomb phase that gives rise to a new source of super-leading logarithms. We then perform a systematic leading-logarithmic resummation of these contributions for 2 → M processes. Finally, we analyze the numerical impact in partonic scattering processes for $t\bar{t}$ production, including a treatment of the Sommerfeld enhancement observed near threshold.

Effective Field Theories of QCD↗