Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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

One-Step Relativistic Driven Similarity Renormalization Group Multireference Perturbation Theory

We present an efficient implementation of a one-step relativistic second-order multireference perturbation theory based on the multireference driven similarity renormalization group (MR-DSRG) using the exact two-component (X2C) Hamiltonian, which we denote X2C-DSRG-MRPT2. We show that the X2C-DSRG-MRPT2 method can accurately capture spin–orbit coupling (SOC) effects in the electronic structure of strongly correlated systems containing elements across the periodic table. We further demonstrate that the X2C-DSRG-MRPT2 method, through its variational treatment of SOC effects, can yield spin–orbit splittings with mean absolute percentage errors consistently below 7% with respect to experimental values for systems containing up to sixth row elements. With its modest computational scaling (fourth power in system size for the perturbative step) and high accuracy, X2C-DSRG-MRPT2 provides a promising avenue for the routine treatment of relativistic effects in strongly correlated molecular systems.

Hamiltonians↗

MiniApp for Density Matrix Renormalization Group Hamiltonian Application Kernel

We present two miniapps that implement the core computational kernel of the DMRG++ application, a generic C++ code that implements the Density Matrix Renormalization Group (DMRG) algorithm. The DMRG++ core Kronecker multiplication kernel is formulated using a batched BLAS approach, with implementation that targets both multi-core CPUs using OpenMP and GPGPU using the MAGMA library. The kernel evaluates the matrix-vector multiplication of the target Hamiltonian matrix used in Lanczos algorithm for computing the lowest eigenvalue and eigenvector. The Hamiltonian matrix is expressed compactly as sums of Kronecker products of small dense matrices. We demonstrate improved performance of the miniapp on synthetic problem, and show the performance of the DMRG++ application using a plugin based on the miniapp. We also present an OpenMP miniapp that explores the use of nested parallel constructs to implement the Kronecker multiplication kernel, exploring the use of nested OpenMP worksharing and tasking abstractions to implement the multi-level parallel multiplication algorithm. The miniapp has been used as a co-design vehicle for evaluating features in the OpenMP-4.5 and upcoming OpenMP-5.0 standards.

Elwasif, Wael↗

Density matrix renormalization group with efficient dynamical electron correlation through range separation

We present a new hybrid multiconfigurational method based on the concept of range-separation that combines the density matrix renormalization group approach with density functional theory. This new method is designed for the simultaneous description of dynamical and static electron-correlation effects in multiconfigurational electronic structure problems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Spin-free formulation of the multireference driven similarity renormalization group: A benchmark study of first-row diatomic molecules and spin-crossover energetics

We report a spin-free formulation of the multireference (MR) driven similarity renormalization group (DSRG) based on the ensemble normal ordering of Mukherjee and Kutzelnigg (1997). This ensemble averages over all microstates of a given total spin quantum number, and therefore, it is invariant with respect to SU(2) transformations. As such, all equations may be reformulated in terms of spin-free quantities and they closely resemble those of spin-adapted closed-shell coupled cluster (CC) theory. The current implementation is used to assess the accuracy of various truncated MR-DSRG methods (perturbation theory up to third order and iterative methods with single and double excitations) in computing the constants of 33 first-row diatomic molecules. The accuracy trends for these first-row diatomics are consistent with our previous benchmark on a small subset of closed-shell diatomic molecules. We then present the first MR-DSRG application on transition-metal complexes by computing the spin splittings of the [Fe(H 2 O) 6 ] 2+ and [Fe(NH 3 ) 6 ] 2+ molecules. A focal point analysis (FPA) shows that third-order perturbative corrections are essential to achieve reasonably converged energetics. The FPA based on the linearized MR-DSRG theory with one- and two-body operators and up to a quintuple- ζ basis set predicts the spin splittings of [Fe(H 2 O) 6 ] 2+ and [Fe(NH 3 ) 6 ] 2+ to be –35.7 and –17.1 kcal mol –1 , respectively, showing good agreement with the results of local CC theory with singles, doubles, and perturbative triples.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A general non-Abelian density matrix renormalization group algorithm with application to the C{sub 2} dimer

We extend our previous work [S. Sharma and G. K.-L. Chan, J. Chem. Phys. 136, 124121 (2012)], which described a spin-adapted (SU(2) symmetry) density matrix renormalization group algorithm, to additionally utilize general non-Abelian point group symmetries. A key strength of the present formulation is that the requisite tensor operators are not hard-coded for each symmetry group, but are instead generated on the fly using the appropriate Clebsch-Gordan coefficients. This allows our single implementation to easily enable (or disable) any non-Abelian point group symmetry (including SU(2) spin symmetry). We use our implementation to compute the ground state potential energy curve of the C{sub 2} dimer in the cc-pVQZ basis set (with a frozen-core), corresponding to a Hilbert space dimension of 10{sup 12} many-body states. While our calculated energy lies within the 0.3 mE{sub h} error bound of previous initiator full configuration interaction quantum Monte Carlo and correlation energy extrapolation by intrinsic scaling calculations, our estimated residual error is only 0.01 mE{sub h}, much more accurate than these previous estimates. Due to the additional efficiency afforded by the algorithm, the excitation energies (T{sub e}) of eight lowest lying excited states: a{sup 3}Π{sub u}, b{sup 3}Σ{sub g}{sup −}, A{sup 1}Π{sub u}, c{sup 3}Σ{sub u}{sup +}, B{sup 1}Δ{sub g}, B{sup ′1}Σ{sub g}{sup +}, d{sup 3}Π{sub g}, and C{sup 1}Π{sub g} are calculated, which agree with experimentally derived values to better than 0.06 eV. In addition, we also compute the potential energy curves of twelve states: the three lowest levels for each of the irreducible representations {sup 1}Σ{sub g}{sup +}, {sup 1}Σ{sub u}{sup +}, {sup 1}Σ{sub g}{sup −}, and {sup 1}Σ{sub u}{sup −}, to an estimated accuracy of 0.1 mE{sub h} of the exact result in this basis.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Renormalization group flow of the Aharonov–Bohm scattering amplitude

The Aharonov–Bohm elastic scattering with incident particles described by plane waves is revisited by using the phase-shifts method. The formal equivalence between the cylindrical Schrödinger equation and the one-dimensional Calogero problem allows us to show that up to two scattering phase-shifts modes in the cylindrical waves expansion must be renormalized. The renormalization procedure introduces new length scales giving rise to spontaneous breaking of the conformal symmetry. The knowledge of the exact beta function permits us to describe the renormalization group flows within the two-parametric family of renormalized Aharonov–Bohm scattering amplitudes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Renormalization group evolution of multi-gluon correlators in high energy QCD

Many-body QCD in leading high energy Regge asymptotics is described by the Balitsky-JIMWLK hierarchy of renormalization group equations for the x evolution of multi-point Wilson line correlators. These correlators are universal and ubiquitous in final states in deeply inelastic scattering and hadronic collisions. For instance, recently measured di-hadron correlations at forward rapidity in deuteron-gold collisions at the Relativistic Heavy Ion Collider (RHIC) are sensitive to four and six point correlators of Wilson lines in the small x color fields of the dense nuclear target. We evaluate these correlators numerically by solving the functional Langevin equation that describes the Balitsky-JIMWLK hierarchy. We compare the results to mean-field Gaussian and large Nc approximations used in previous phenomenological studies. We comment on the implications of our results for quantitative studies of multi-gluon final states in high energy QCD.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Renormalization group analysis of the Reynolds stress transport equation

The pressure velocity correlation and return to isotropy term in the Reynolds stress transport equation are analyzed using the Yakhot-Orszag renormalization group. The perturbation series for the relevant correlations, evaluated to lowest order in the epsilon-expansion of the Yakhot-Orszag theory, are infinite series in tensor product powers of the mean velocity gradient and its transpose. Formal lowest order Pade approximations to the sums of these series produce a fast pressure strain model of the form proposed by Launder, Reece, and Rodi, and a return to isotropy model of the form proposed by Rotta. In both cases, the model constant are computed theoretically. The predicted Reynolds stress ratios in simple shear flows are evaluated and compared with experimental data. The possibility is discussed of driving higher order nonlinear models by approximating the sums more accurately.

Rubinstein, R.↗

Spin filtering in a double quantum dot device: Numerical renormalization group study of the internal structure of the Kondo state

A double quantum dot device, connected to two channels that only interact through interdot Coulomb repulsion, is analyzed using the numerical renormalization group technique. Using a two-impurity Anderson model, and realistic parameter values [S. Amasha, A. J. Keller, I. G. Rau, A. Carmi, J. A. Katine, H. Shtrikman, Y. Oreg, and D. Goldhaber-Gordon, Phys. Rev. Lett. 110, 046604 (2013)], it is shown that, by applying a moderate magnetic field and independently adjusting the gate potential of each quantum dot at half-filling, a spin-orbital SU(2) Kondo state can be achieved where the Kondo resonance originates from spatially separated parts of the device. Our results clearly link this spatial separation effect to currents with opposing spin polarizations in each channel, i.e., the device acts as a spin filter. In addition, an experimental probe of this polarization effect is suggested, pointing to the exciting possibility of experimentally probing the internal structure of an SU(2) Kondo state.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Spectral functions with the density matrix renormalization group: Krylov-space approach for correction vectors

Frequency-dependent correlations, such as the spectral function and the dynamical structure factor, help illustrate condensed matter experiments. Within the density matrix renormalization group (DMRG) framework, an accurate method for calculating spectral functions directly in frequency is the correction-vector method. The correction vector can be computed by solving a linear equation or by minimizing a functional. Our paper proposes an alternative to calculate the correction vector: to use the Krylov-space approach. Further, this paper also studies the accuracy and performance of the Krylov-space approach, when applied to the Heisenberg, the t-J, and the Hubbard models. The cases we studied indicate that the Krylov-space approach can be more accurate and efficient than the conjugate gradient, and that the error of the former integrates best when a Krylov-space decomposition is also used for ground state DMRG.

36 MATERIALS SCIENCE↗

Turbulent transport of a passive-scalar field by using a renormalization-group method

A passive-scalar field is considered to evolve under the influence of a turbulent fluid governed by the Navier-Stokes equation. Turbulent-transport coefficients are calculated by small-scale elimination using a renormalization-group method. Turbulent processes couple both the viscosity and the diffusivity. In the absence of any correlation between the passive-scalar fluctuations and any component of the fluid velocity, the renormalized diffusivity is essentially the same as if the fluid velocity were frozen, although the renormalized equation does contain higher-order nonlinear terms involving viscosity. This arises due to the nonlinear interaction of the velocity with itself. In the presence of a finite correlation, the turbulent diffusivity becomes coupled with both the velocity field and the viscosity. There is then a dependence of the turbulent decay of the passive scalar on the turbulent Prandtl number.

Hossain, Murshed↗

Generator-coordinate reference states for spectra and 0 ν β β decay in the in-medium similarity renormalization group

We use a reference state based on symmetry-restored states from deformed mean-field or generator-coordinate-method (GCM) calculations in conjunction with the in-medium similarity- renormalization group (IMSRG) to compute spectra and matrix elements for neutrinoless double- beta (0$v$ββ) decay. Because the decay involves ground states from two nuclei, we use evolved operators from the IMSRG in one nucleus in a subsequent GCM calculation in the other. We benchmark the resulting IMSRG+GCM method against complete shell-model diagonalization for both the energies of low-lying states in 48 Ca and 48 Ti and the 0$v$ββ matrix element for the decay of 48 Ca, all in a single valence shell. Our approach produces better spectra than either the IMSRG with a spherical-mean-field reference or GCM calculations with unevolved operators. For the 0$v$ββ matrix element the improvement is slight.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Nonlinear relativistic plasma resonance: Renormalization group approach

An analytical solution to the nonlinear set of equations describing the electron dynamics and electric field structure in the vicinity of the critical density in a nonuniform plasma is constructed using the renormalization group approach with allowance for relativistic effects of electron motion. It is demonstrated that the obtained solution describes two regimes of plasma oscillations in the vicinity of the plasma resonance— stationary and nonstationary. For the stationary regime, the spatiotemporal and spectral characteristics of the resonantly enhanced electric field are investigated in detail and the effect of the relativistic nonlinearity on the spatial localization of the energy of the plasma relativistic field is considered. The applicability limits of the obtained solution, which are determined by the conditions of plasma wave breaking in the vicinity of the resonance, are established and analyzed in detail for typical laser and plasma parameters. The applicability limits of the earlier developed nonrelativistic theories are refined.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Random Crystal Field in a Mixed Spin S = 1/2 and S = 3/2 Ising Model by Renormalization Group Theory

In this work, we study the Ising model with mixed spins S = 1/2 and S = 3/2 on the hypercubic lattice and the random crystal field at the sublattice with S = 3/2 described by a two peaks law. To achieve this, we use an approximation of position space renormalization group (PSRG) namely Migdal-Kadanoff in which we use both the free energy derivative and the flow in the parameter space of the Hamiltonian. For all values of the random probability, the critical behavior is determined via the critical exponents at the second-order fixed points while at low temperatures; the discontinuities of the of the free energy derivative provide the positions of the first-order transitions. The introduction of a minimal amount of disorder causes a change in the phase diagram showing the relevance of disorder for d = 2 and d = 3. The second-order transition remains always with the same critical exponents as those of the pure model, and a new first-order transition appears at very low temperature. Also, a comparison with other similar works is given.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Functional renormalization group study of nuclear and neutron matter

A chiral model based on nucleons interacting via boson exchange is investigated. Fluctuation effects are included consistently beyond the mean-field approximation in the framework of the functional renormalization group. The liquid-gas phase transition of symmetric nuclear matter is studied in detail. No sign of a chiral restoration transition is found up to temperatures of about 100 MeV and densities of at least three times the density of normal nuclear matter. Moreover, the model is extended to asymmetric nuclear matter and the constraints from neutron star observations are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Renormalization Group Tutorial

Complex physical systems sometimes have statistical behavior characterized by power- law dependence on the parameters of the system and spatial variability with no particular characteristic scale as the parameters approach critical values. The renormalization group (RG) approach was developed in the fields of statistical mechanics and quantum field theory to derive quantitative predictions of such behavior in cases where conventional methods of analysis fail. Techniques based on these ideas have since been extended to treat problems in many different fields, and in particular, the behavior of turbulent fluids. This lecture will describe a relatively simple but nontrivial example of the RG approach applied to the diffusion of photons out of a stellar medium when the photons have wavelengths near that of an emission line of atoms in the medium.

Bell, Thomas L.↗

Renormalization group analysis for the quasi-1D superconductor BaFe 2 S 3

Motivated by the discovery of superconductivity in the two-leg, quasi-one dimensional ladder compound, BaFe 2 S 3 we present a renormalization group study of electrons moving on a two leg, two orbital ladder, subjected to Hubbard repulsion U and Hund’s coupling J . In our calculations, we adopt tightbinding parameters obtained from ab-initio studies on this material. At incommensurate filling, the long wavelength analysis displays four phases as a function of 0 ≤ J/U < 1. We show that a fully gapped superconductor is stabilized at sufficiently large Hund’s coupling, the relative phases at the three Fermi points are “+,-,-”. By contrast, when the system is tuned to half filling, Umklapp scattering gives rise to Mott insulating phases. Furthermore, we discuss the general implications of our study for the broad class of iron-based superconductors.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗