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

Gravitational Fluctuations as an Alternative to Inflation III. Numerical Results

Power spectra play an important role in the theory of inflation, and their ability to reproduce current observational data to high accuracy is often considered a triumph of inflation, largely because of a lack of credible alternatives. In previous work we introduced an alternative picture for the cosmological power spectra based on the nonperturbative features of the quantum version of Einstein’s gravity, instead of currently popular inflation models based on scalar fields. The key ingredients in this new picture are the appearance of a nontrivial gravitational vacuum condensate (directly related to the observed cosmological constant), and a calculable renormalization group running of Newton’s G on cosmological scales. More importantly, one notes the absence of any fundamental scalar fields in this approach. Results obtained previously were largely based on a semi-analytical treatment, and thus, while generally transparent in their implementation, often suffered from the limitations of various approximations and simplifying assumptions. In this work, we extend and refine our previous calculations by laying out an updated and extended analysis, which now utilizes a set of suitably modified state-of-the-art numerical programs (ISiTGR, MGCAMB and MGCLASS) developed for observational cosmology. As a result, we are able to remove some of the approximations employed in our previous studies, leading to a number of novel and detailed physical predictions. These should help in potentially distinguishing the vacuum condensate picture of quantum gravity from that of other models such as scalar field inflation. Here, besides the matter power spectrum P m (k) , we work out, in detail, predictions for what are referred to as the TT, TE, EE, BB angular spectra, as well as their closely related lensing spectra. However, the current limited precision of observational data today (especially on large angular scales) does not allow us yet to clearly prove or disprove either set of ideas. Nevertheless, by exploring in more details the relationship between gravity and cosmological matter and radiation both analytically and numerically, together with an expected future influx of increasingly accurate observational data, one can hope that the new quantum gravitational picture can be subjected to further stringent tests in the near future.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Magnetism in the two-dimensional dipolar XY model

Here, motivated by a recent experiment on a square-lattice Rydberg atom array realizing a long-range dipolar XY model [C. Chen et al., Nature (London) 616, 691 (2023)], we numerically study the model's equilibrium properties. We obtain the phase diagram, critical properties, entropies, variance of the magnetization, and site-resolved correlation functions. We consider both ferromagnetic and antiferromagnetic interactions and apply quantum Monte Carlo and pseudo-Majorana functional renormalization group techniques, generalizing the latter to a U⁡(1) symmetric setting. Our simulations perform extensive thermometry in dipolar Rydberg atom arrays and establish conditions for adiabaticity and thermodynamic equilibrium. On the ferromagnetic side of the experiment, we determine the entropy per particle S/N≈0.5, close to the one at the critical temperature, S c /N=0.585⁢(15). The simulations suggest the presence of an out-of-equilibrium plateau at large distances in the correlation function, thus motivating future studies on the nonequilibrium dynamics of the system.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Efficient simulation of moiré materials using the density matrix renormalization group

Here we present an infinite density-matrix renormalization group (DMRG) study of an interacting continuum model of twisted bilayer graphene (tBLG) near the magic angle. Because of the long-range Coulomb interaction and the large number of orbital degrees of freedom, tBLG is difficult to study with standard DMRG techniques - even constructing and storing the Hamiltonian already poses a major challenge. To overcome these difficulties, we use a recently developed compression procedure to obtain a matrix product operator representation of the interacting tBLG Hamiltonian which we show is both efficient and accurate even when including the spin, valley, and orbital degrees of freedom. To benchmark our approach, we focus mainly on the spinless, single-valley version of the problem where, at half filling, we find that the ground state is a nematic semimetal. Remarkably, we find that the ground state is essentially a k-space Slater determinant, so that Hartree-Fock and DMRG give virtually identical results for this problem. Our results show that the effects of long-range interactions in magic angle graphene can be efficiently simulated with DMRG and open up a new route for numerically studying strong correlation physics in spinful, two-valley tBLG, and other moiré materials in future work.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Entanglement properties of disordered quantum spin chains with long-range antiferromagnetic interactions

Entanglement measures are useful tools in characterizing otherwise unknown quantum phases and indicating transitions between them. Here we examine the concurrence and entanglement entropy in quantum spin chains with random long-range couplings, spatially decaying with a power-law exponent $\textit{α}$. Using the strong disorder renormalization group (SDRG) technique, we find by analytical solution of the master equation a strong disorder fixed point, characterized by a fixed point distribution of the couplings with a finite dynamical exponent, which describes the system consistently in the regime $α > \frac{1}{2}$ . A numerical implementation of the SDRG method yields a power-law spatial decay of the average concurrence, which is also confirmed by exact numerical diagonalization. However, we find that the lowest-order SDRG approach is not sufficient to obtain the typical value of the concurrence. We therefore implement a correction scheme which allows us to obtain the leading-order corrections to the random singlet state. This approach yields a power-law spatial decay of the typical value of the concurrence, which we derive both by a numerical implementation of the corrections and by analytics. Next, using numerical SDRG, the entanglement entropy (EE) is found to be logarithmically enhanced for all $\textit{α}$, corresponding to a critical behavior with an effective central charge $\textit{c}$ = ln(2), independent of $\textit{α}$. This is confirmed by an analytical derivation. Using numerical exact diagonalization (ED), we confirm the logarithmic enhancement of the EE and a weak dependence on $\textit{α}$. For a wide range of partition size $\textit{l}$, the EE fits a critical behavior with a central charge close to $\textit{c}$ = 1, which is the same as for the clean Haldane-Shastry model with a power-law-decaying interaction with $\textit{α}$ = 2. Only for small $l \ll L$, in a range which increases with the number of spins N, we find deviations which are rather consistent with the strong disorder fixed point central charge $\textit{c}$ = ln(2). Furthermore, we find using ED that the concurrence shows power-law decay, albeit with smaller power exponents than obtained by SDRG. Finally, we also present results obtained with DMRG and find agreement with ED for sufficiently small $\textit{α}$ < 2, whereas for larger $\textit{α}$ DMRG tends to underestimate the entanglement entropy and finds a faster decaying concurrence.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Triplet pair density wave superconductivity on the π -flux square lattice

Pair-density waves (PDWs) are superconducting states that spontaneously break translation symmetry in systems with time-reversal symmetry (TRS). Evidence for PDWs has been seen in several recent experiments, as well as in the pseudogap regime in cuprates. Theoretical understanding of PDWs has been largely restricted to phenomenological and numerical studies, while microscopic theories typically require strong coupling or fine tuning. In this work, we provide a symmetry-based mechanism under which PDWs emerge as a weak-coupling instability of a two-dimensional TRS metal. Combining mean-field and renormalization-group analyses, we identify a weak-coupling instability towards a triplet PDW realized in the π-flux square lattice model with on-site repulsion and moderate nearest-neighbor attraction when the Fermi level crosses Van Hove singularities at 1/4 and 3/4 fillings. Further, this PDW is protected by the magnetic translation symmetries characteristic of Hofstadter systems, of which the π-flux lattice is a special time-reversal-symmetric case.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Luther-Emery liquid and dominant singlet superconductivity in the hole-doped Haldane spin-1 chain

We investigate the pairing tendencies in the hole-doped Haldane spin-1 chain. To allow for doping, we extend the original spin chain Hamiltonian into a fermionic model involving a two-orbital Hubbard chain at intermediate or strong repulsive interaction strengths U and for degenerate orbitals. At half filling and large U, the ferromagnetic Hund's coupling, J H , generates effective spin-1 moments, with antiferromagnetic correlations between sites. Using large-scale density matrix renormalization group calculations, we accurately study the system's behavior under light hole-doping. For U = 1.6 in units of the noninteracting bandwidth and for J H /U ≳ 0.275, we find that singlet pairing dominates the long-distance physics, establishing this system as a promising platform for repulsively mediated superconductivity. We provide concrete examples of materials that could realize the physics described here. We also provide evidence that the system approaches a Luther-Emery liquid state at large system sizes, reminiscent of the behavior of doped one-orbital two-leg ladders at weak coupling, which also have superconducting tendencies. The numerically calculated central charge approaches one in the thermodynamic limit, indicating a single gapless mode as is expected for the Luther-Emery state. Exponents characterizing the power-law decays of singlet pair-pair and charge density-density correlations are determined, and found to approximately satisfy the Luther-Emery identity.

1-dimensional systems↗
DOE OSTI2024

A numerical study of bounds in the correlations of fractional quantum Hall states

We numerically compute the guiding center static structure factor \bar S(k) S ‾ ( k ) of various fractional quantum Hall (FQH) states to \mathcal{O}(k\ell)^6 𝒪 ( k ℓ ) 6 where k k is the wavenumber and \ell ℓ is the magnetic length. Employing density matrix renormalization group on an infinite cylinder of circumference L_y L y , we study the two-dimensional limit using L_y/\xi \gg 1 L y / ξ ≫ 1 , where \xi ξ is the correlation length. The main findings of our work are: 1) the ground states that deviate away from the ideal conformal block wavefunctions, do not saturate the Haldane bound, and 2) the coefficient of O(k\ell)^6 O ( k ℓ ) 6 term appears to be bounded above by a value predicted by field theories proposed in the literature. The first finding implies that the graviton mode is not maximally chiral for experimentally relevant FQH states.

Kumar, Prashant (ORCID:0000000258002768)↗
DOE OSTI2020

Many-Body Effects in FeN4 Center Embedded in Graphene

We introduce a computational approach to study porphyrin-like transition metal complexes, bridging density functional theory and exact many-body techniques, such as the density matrix renormalization group (DMRG). We first derive a multi-orbital Anderson impurity Hamiltonian starting from first principles considerations that qualitatively reproduce generalized gradient approximation (GGA)+U results when ignoring inter-orbital Coulomb repulsion U ′ and Hund exchange J. An exact canonical transformation is used to reduce the dimensionality of the problem and make it amenable to DMRG calculations, including all many-body terms (both intra- and inter-orbital), which are treated in a numerically exact way. We apply this technique to FeN 4 centers in graphene and show that the inclusion of these terms has dramatic effects: as the iron orbitals become single occupied due to the Coulomb repulsion, the inter-orbital interaction further reduces the occupation, yielding a non-monotonic behavior of the magnetic moment as a function of the interactions, with maximum polarization only in a small window at intermediate values of the parameters. Furthermore, U ′ changes the relative position of the peaks in the density of states, particularly on the iron d z 2 orbital, which is expected to affect the binding of ligands greatly.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Anomalous dimensions from gauge couplings in SMEFT with right-handed neutrinos

Standard Model Neutrino Effective Field Theory (SMNEFT) is an effective theory with Standard Model (SM) gauge-invariant operators constructed only from SM and right-handed neutrino fields. For the full set of dimension-six SMNEFT operators, we present the gauge coupling terms of the one-loop anomalous dimension matrix for renormalization group evolution (RGE) of the Wilson coefficients between a new physics scale and the electroweak scale. We find that the SMNEFT operators can be divided into five subsets which are closed under RGE. Our results apply for both Dirac and Majorana neutrinos. We also discuss the operator mixing pattern numerically and comment on some interesting phenomenological implications.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Fragility of the nematic spin liquid induced by diagonal couplings in the square-lattice SU(3) model

Here, we present a large-scale density matrix renormalization group (DMRG) study of the spin-1 SU(3) bilinear-biquadratic model on the square lattice, which was suggested to host a nematic spin liquid state in recent DMRG calculations. We report that this spin liquid appears to strongly compete with a three-sublattice magnetic order. To further study the competition between the two states, and the reason of the emergent nematic spin liquid, we included an additional next-nearest-neighbor SU(3) symmetric interactions along one of the two plaquette diagonal directions. This allows to tune the square-lattice model to the triangular-lattice model. By computing spin correlation functions and various order parameters, we find that the three-sublattice order may develop at infinitesimal additional new coupling, at least within the precision of our study. Compared with the previous findings that the nematic spin liquid is stable in extended parameter regions with additional couplings that respect the lattice symmetries of the square lattice, we argue that here the diagonal couplings, which frustrate the bipartite-lattice structure, rapidly suppress the two-sublattice fluctuations and the three-sublattice order thus wins. This numerical result is consistent with the conjecture that the nematic spin liquid emerges from the competition between two- and three-sublattice fluctuations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2004

Complexity Induced Anisotropic Bimodal Intermittent Turbulence in Space Plasmas

The "physics of complexity" in space plasmas is the central theme of this exposition. It is demonstrated that the sporadic and localized interactions of magnetic coherent structures arising from the plasma resonances can be the source for the coexistence of nonpropagating spatiotemporal fluctuations and propagating modes. Non-Gaussian probability distribution functions of the intermittent fluctuations from direct numerical simulations are obtained and discussed. Power spectra and local intermittency measures using the wavelet analyses are presented to display the spottiness of the small-scale turbulent fluctuations and the non-uniformity of coarse-grained dissipation that can lead to magnetic topological reconfigurations. The technique of the dynamic renormalization group is applied to the study of the scaling properties of such type of multiscale fluctuations. Charged particle interactions with both the propagating and nonpropagating portions of the intermittent turbulence are also described.

Chang, Tom↗
DOE OSTI2023

High precision tests of QCD without scale or scheme ambiguities: The 40th anniversary of the Brodsky–Lepage–Mackenzie method

A key issue in making precise predictions in QCD is the uncertainty in setting the renormalization scale μ r and thus determining the correct values of the QCD running coupling α s (μ r ) at each order in the perturbative expansion of a QCD observable. It has often been conventional to simply set the renormalization scale to the typical scale of the process Q and vary it in the range μ r $\in$ [Q/2, 2Q] in order to estimate the theoretical error. This is the practice of Conventional Scale Setting (CSS). The resulting CSS prediction will however depend on the theorist’s choice of renormalization scheme and the resulting pQCD series will diverge factorially. It will also disagree with renormalization scale setting used in QED and electroweak theory thus precluding grand unification. A solution to the renormalization scale-setting problem is offered by the Principle of Maximum Conformality (PMC), which provides a systematic way to eliminate the renormalization scale-and-scheme dependence in perturbative calculations. The PMC method has rigorous theoretical foundations, it satisfies Renormalization Group Invariance (RGI) and preserves all self-consistency conditions derived from the renormalization group. The PMC cancels the renormalon growth, reduces to the Gell-Mann–Low scheme in the N c → 0 Abelian limit and leads to scale- and scheme-invariant results. The PMC has now been successfully applied to many high-energy processes. In this article we summarize recent developments and results in solving the renormalization scale and scheme ambiguities in perturbative QCD. In particular, we present a recently developed method the PMC ∞ and its applications, comparing the results with CSS. The method preserves the property of renormalizable SU(N)/U(1) gauge theories defined as Intrinsic Conformality (iCF). This property underlies the scale invariance of physical observables and leads to a remarkably efficient method to solve the conventional renormalization scale ambiguity at every order in pQCD. This new method reflects the underlying conformal properties displayed by pQCD at NNLO, eliminates the scheme dependence of pQCD predictions and is consistent with the general properties of the PMC. A new method to identify conformal and β-terms, which can be applied either to numerical or to theoretical calculations is also shown. We present results for the thrust and C-parameter distributions in e + e - annihilation showing errors and comparison with the CSS. We also show results for a recent innovative comparison between the CSS and the PMC ∞ applied to the thrust distribution investigating both the QCD conformal window and the QED N c → 0 limit. In order to determine the thrust distribution along the entire renormalization group flow from the highest energies to zero energy, we consider the number of flavors near the upper boundary of the conformal window. In this flavor-number regime the theory develops a perturbative infrared interacting fixed point. These results show that PMC ∞ leads to higher precision and introduces new interesting features in the PMC. In fact, this method preserves with continuity the position of the peak, showing perfect agreement with the experimental data already at NNLO. We also show a detailed comparison of the PMC ∞ with the other PMC approaches: the multi-scale-setting approach (PMCm) and the single-scale-setting approach (PMCs) by comparing their predictions for three important fully integrated quantities R e+e- , R$_Τ$ and Γ (H→$b\bar{b}$) up to the four-loop accuracy.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Accurate localization of Kosterlitz-Thouless-type quantum phase transitions for one-dimensional spinless fermions

We investigate the charge-density wave (CDW) transition for one-dimensional spinless fermions at half band filling with nearest-neighbor electron transfer amplitude t and interaction V. The model is equivalent to the anisotropic XXZ Heisenberg model for which the Bethe Ansatz provides an exact solution. For V>V c =2t, the CDW order parameter and the single-particle gap are finite but exponentially small, as is characteristic for a Kosterlitz-Thouless transition. It is notoriously difficult to locate such infinite-order phase transitions in the phase diagram using approximate analytical and numerical approaches. Second-order Hartree-Fock theory is qualitatively applicable for all interaction strengths, and predicts the CDW transition to occur at V$^{(2)}_{c,2}$≈1.5t. Second-order Hartree Fock theory is almost variational because the density of quasiparticle excitations is small. We apply the density-matrix renormalization group (DMRG) for periodic boundary conditions for system sizes up to 514 sites, which permits a reliable extrapolation of all physical quantities to the thermodynamic limit, apart from the critical region. We investigate the ground-state energy, the gap, the order parameter, the momentum distribution, the quasiparticle density, and the density-density correlation function to locate V c from the DMRG data. In conclusion, tracing the breakdown of the Luttinger liquid and the peak in the quasiparticle density at the band edge permits us to reproduce V c with an accuracy of one percent.

1-dimensional spin chains↗
DOE OSTI2023

Triangular lattice Hubbard model physics at intermediate temperatures

Moiré systems offer an exciting playground to study the many-body effects of strongly correlated electrons in regimes that are not easily accessible in conventional material settings. Motivated by a recent experiment on WSe 2 /WS 2 moiré bilayers [Y. Tang et al., Nature (London) 579, 353 (2020)], which realized a triangular superlattice with a small hopping (of approximately 10 Kelvin), with tunable density of holes, we explore the Hubbard model on the triangular lattice for intermediate temperatures t ≲ T < U. Employing finite-temperature Lanczos calculations and closely following the fitting protocols used in the experiment, we recover the observed trends in the reported Curie-Weiss temperature Θ with filling, using the reported interaction strength U/t=20. Here, we focus on the large increase of |Θ| on decreasing the density below half filling and the sign change of Θ at higher fillings, which signals the onset of ferromagnetism. The increase in |Θ| is also seen in the t–J model (the low-energy limit of the Hubbard model) in the intermediate temperature range, which we clarify is opposite to the trend in its high temperature limit. Differences between the low, intermediate, and high temperature regimes are discussed. Our numerical calculations also capture the crossover between short-range antiferromagnetic to ferromagnetic order in the intermediate temperature regime, a result broadly consistent with the experimental findings. We find that this behavior is a finite-temperature remnant of the underlying zero-temperature phase transition, which we explore with ground-state density matrix renormalization group calculations. We provide evidence of ferromagnetism characterized by weak (but robust) correlations that explain the small Θ seen in the experiment.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Non-linear dynamics of jet quenching

We develop a comprehensive analytic framework for jet quenching in QCD media, based on a medium-induced parton cascade sourced by collinear virtual splittings. We show that the energy flow out of the jet cone, driven by turbulent gluon cascades, is governed by a non-linear rate equation that resums gluon splittings at arbitrary angles and is enhanced by the medium length, L. The solution of this equation sets the initial condition for a non-linear DGLAP-like evolution equation, which describes the collinear early vacuum cascade resolved by the medium at angles exceeding the medium resolution angle, θ c . For asymptotic jet energies, the medium-induced cascade displays an exponential behavior that generalizes the Poisson-like distribution of parton energy loss. This formulation enables the resummation of leading contributions in α s ln(1/R), and α s ln(R/θ c ), and powers of α s L. We briefly explore the limit of strong quenching, where analytic treatments are feasible, offering insights into the impact of parton cascades on jet quenching. These results provide guidance for future numerical simulations and analytical investigations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1995

Evaluation of Turbulence Models for Unsteady Flows of an Oscillating Airfoil

Unsteady flowfields of a two-dimensional oscillating airfoil are calculated using an implicit, finite-difference, Navier Stokes numerical scheme. Five widely used turbulence models are used with the numerical scheme to assess the accuracy and suitability of the models for simulating the retreating blade stall of helicopter rotor in forward flight. Three unsteady flow conditions corresponding to an essentially attached flow, light-stall, and deep-stall cases of an oscillating NACA 0015 wing experiment were chosen as test cases for computations. Results of unsteady airloads hysteresis curves, harmonics of unsteady pressures, and instantaneous flowfield patterns are presented. Some effects of grid density, time-step size, and numerical dissipation on the unsteady solutions relevant to the evaluation of turbulence models are examined. Comparison of unsteady airloads with experimental data show that all models tested are deficient in some sense and no single model predicts airloads consistently and in agreement with experiment for the three flow regimes. The chief findings are that the simple algebraic model based on the renormalization group theory (RNG) offers some improvement over the Baldwin Lomax model in all flow regimes with nearly same computational cost. The one-equation models provide significant improvement over the algebraic and the half-equation models but have their own limitations. The Baldwin-Barth model overpredicts separation and underpredicts reattachment. In contrast, the Spalart-Allmaras model underpredicts separation and overpredicts reattachment.

Srinivasan, G. R.↗
DOE OSTI2024

Relevance of one-loop SMEFT matching in the 2HDM

The two-Higgs doublet model (2HDM) is a well understood alternative to the Standard Model of particle physics. If the new particles included in the 2HDM are at an energy scale much greater than the weak scale, the theory can be matched to the Standard Model effective field theory (SMEFT). We compute for the first time the complete one-loop matching at dimension-6. We compare its numerical impact with that of tree-level matching at dimension-8 by performing a global fit to single Higgs and precision electroweak measurements, and we emphasize the importance of comparing one-loop SMEFT results with corresponding one-loop results in the full 2HDM model. In the SMEFT, we consider the relative importance of both one-loop matching and the inclusion of renormalization group evolution. Our results demonstrate the necessity of studying the impact of various expansions to quantify the uncertainties of the SMEFT matching. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗