Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Green's function methods”

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 361 records · Page 20

Reexamining the relation between the binding energy of finite nuclei and the equation of state of infinite nuclear matter

The energy density is calculated in coordinate space for 12 C, 40 Ca, 48 Ca, and 208 Pb using a dispersive optical model constrained by all relevant data including the corresponding energy of the ground state. The energy density of 8 Be is also calculated using the Green’s function Monte-Carlo method employing the Argonne/Urbana two and three-body interactions. The nuclear interior minimally contributes to the total binding energy due to the 4πr 2 phase space factor. Thus, the volume contribution to the energy in the interior is not well constrained. The dispersive-optical-model energy densities are in good agreement with ab initio self-consistent Green’s function calculations of infinite nuclear matter restricted to treat only short-range and tensor correlations. These results call into question the degree to which the equation of state for nuclear matter is constrained by the empirical mass formula. In particular, the results in this work indicate that saturated nuclear matter does not require the canonical value of 16 MeV binding per particle but only about 13-14 MeV when the interior of 208 Pb is considered.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Numerically exact generalized Green's function cluster expansions for electron-phonon problems

We generalize the family of approximate momentum average methods to formulate a numerically exact, convergent hierarchy of equations whose solution provides an efficient algorithm to compute the Green's function of a particle dressed by bosons suitable in the entire parameter regime. We use this approach to extract ground-state properties and spectral functions. Our approximation-free framework, dubbed the generalized Green's function cluster expansion (GGCE), allows access to exact numerical results in the extreme adiabatic limit, where many standard methods struggle or completely fail. We showcase the performance of the method, specializing three important models of charge-boson coupling in solids and molecular complexes: the molecular Holstein model, which describes coupling between charge density and local distortions, the Peierls model, which describes modulation of charge hopping due to intersite distortions, and a more complex Holstein + Peierls system with couplings to two different phonon modes, paradigmatic of charge-lattice interactions in organic crystals. Furthermore, the GGCE serves as an efficient approach that can be systematically extended to different physical scenarios, thus providing a tool to model the frequency dependence of dressed particles in realistic settings.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Inverse scattering theory: Inverse scattering series method for one dimensional non-compact support potential

The reversion of the Born-Neumann series of the Lippmann-Schwinger equation is one of the standard ways to solve the inverse acoustic scattering problem. One limitation of the current inversion methods based on the reversion of the Born-Neumann series is that the velocity potential should have compact support. However, this assumption cannot be satisfied in certain cases, especially in seismic inversion. Based on the idea of distorted wave scattering, we explore an inverse scattering method for velocity potentials without compact support. The strategy is to decompose the actual medium as a known single interface reference medium, which has the same asymptotic form as the actual medium and a perturbative scattering potential with compact support. After introducing the method to calculate the Green’s function for the known reference potential, the inverse scattering series and Volterra inverse scattering series are derived for the perturbative potential. Analytical and numerical examples demonstrate the feasibility and effectiveness of this method. Besides, to ensure stability of the numerical computation, the Lanczos averaging method is employed as a filter to reduce the Gibbs oscillations for the truncated discrete inverse Fourier transform of each order. Our method provides a rigorous mathematical framework for inverse acoustic scattering with a non-compact support velocity potential.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Estimates for the Green's function and parameters of exponential dichotomy of a hyperbolic operator semigroup and linear relations

By applying Lyapunov's equation, the method of similar operators, and the methods of harmonic analysis, we obtain estimates for the parameters of exponential dichotomy and for the Green's function constructed for a hyperbolic operator semigroup and a hyperbolic linear relation. Estimates are obtained using quantities which are determined by the resolvent of the infinitesimal operator of the operator semigroup and of the linear relation. Bibliography: 51 titles.

97 MATHEMATICS AND COMPUTING↗

Classical and quantum computing of shear viscosity for ( 2 + 1 ) D SU(2) gauge theory

We perform a nonperturbative calculation of the shear viscosity for ( 2 + 1 )-dimensional SU(2) gauge theory by using the lattice Hamiltonian formulation. The retarded Green’s function of the stress-energy tensor is calculated from real time evolution via exact diagonalization of the lattice Hamiltonian with a local Hilbert space truncation, and the shear viscosity is obtained via the Kubo formula. When taking the continuum limit, we account for the renormalization group flow of the coupling but no additional operator renormalization. We find the ratio of the shear viscosity and the entropy density η s is consistent with a well-known holographic result 1 4 π at several temperatures on a 4 × 4 honeycomb lattice with the local electric representation truncated at j max = 1 2 . We also find the ratio of the spectral function and frequency ρ x y ( ω ) ω exhibits a peak structure when the frequency is small. Both the exact diagonalization method and simple matrix product state classical simulation method beyond j max = 1 2 on bigger lattices require exponentially growing resources. So we develop a quantum computing method to calculate the retarded Green’s function and analyze various systematics of the calculation including j max truncation and finite size effects, Trotter errors and the thermal state preparation efficiency. Our thermal state preparation method still requires resources that grow exponentially with the lattice size, but with a very small prefactor at high temperature. We test our quantum circuit on both the Quantinuum emulator and the IBM simulator for a small lattice and obtain results consistent with the classical computing ones. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Contactless electronic transport in a bio-molecular junction

Molecular electronics hold promise for next generation ultra-low power, nano-scale integrated electronics. The main challenge in molecular electronics is to make a reliable interface between molecules and metal electrodes. Interfacing metals and molecules detrimentally affects the characteristics of nano-scale molecular electronic devices. It is therefore essential to investigate alternative arrangements such as contact-less tunneling gaps wherever such configurations are feasible. We conduct ab initio density functional theory and non-equilibrium Green's functions calculations to investigate the transport properties of a biocompatible glycine molecular junction. By analyzing the localized molecular orbital energy distributions and transmission probabilities in the transport-gap, we find a glycine molecule confined between two gold electrodes, without making a contact, is energetically stable and possesses high tunneling current resembling an excellent ohmic-like interface.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

To the Theory of Inhomogeneous Electron Gas

Nonrigorous character of the density-functional theory for inhomogeneous electron gas based on the hypothesis assuming the existence of a universal density functional is demonstrated. A single-particle density matrix must be determined to calculate the ground-state energy of a finite system with a finite number of electrons. A single-particle Green function can be used to unambiguously determine the ground-state energy of an inhomogeneous electron system that satisfies the thermodynamic limit.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Ab Initio Study of Stability, Local Order, and Phase Diagram For a Series of bcc-based Transition Metal Alloys

In this work, a parameter-free electronic structure approach is applied to the study of stability and chemical order in the 15 substitutional body-centered cubic (bcc)-based alloys made of the six transition metals of groups 5 (V, Nb, Ta) and 6 (Cr, Mo, W) of the periodic table. The method is based on a Green’s function description of the electronic structure of the random alloys. Configurational order is treated within the generalized perturbation method, and temperature effects are examined with a generalized mean-field approach. In contrast to the results summarized in the assessed phase diagrams, stability and ordering trends are predicted in a broad range of alloy composition for at least seven alloys, and explanation is found in their electronic structure properties. Short-range order results, thermodynamics analysis, and bcc-based phase diagrams are also presented.

36 MATERIALS SCIENCE↗

A full-potential approach to the relativistic single-site Green's function

One major purpose of studying the single-site scattering problem is to obtain the scattering matrices and differential equation solutions indispensable to multiple scattering theory (MST) calculations. On the other hand, the single-site scattering itself is also appealing because it reveals the physical environment experienced by electrons around the scattering center. In this study, we demonstrate a new formalism to calculate the relativistic full-potential single-site Green's function. We implement this method to calculate the single-site density of states and electron charge densities. Lastly, the code is rigorously tested and with the help of Krein's theorem, the relativistic effects and full potential effects in group V elements and noble metals are thoroughly investigated.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Thermodynamic, Kinetic, and Physical Properties of Nb-U (Niobium-Uranium)

Based on an improved coupling between ab initio energetics output, phenomenological thermodynamics based on the CALPHAD (CALculation of Phase Diagrams) approach, and the use of the Thermo-Calc software, together with experimental data (whenever available), a thermodynamic assessment of the Nb-U alloy system was carried out. We made use of the full array of modern electronic-structure techniques based on density functional theory, including ab initio fully relativistic full-potential linear-muffin orbital method (FR-FP-LMTO), methods, the Scalar Relativistic Korringa-Kohn-Rostoker Atomic-Sphere Approximation (SR-KKR-ASA), and the fully relativistic or scalar relativistic Exact muffin-orbital (SR- or FR-EMTO) method, both based on the Green function formalism that has the advantage of being able to treat chemically random configuration (i.e., solid solution) with the so-called coherent potential approximation (CPA). All electronic structure calculations were performed within the Generalized Gradient Approximation (GGA) for the treatment of electron correlations. The results from this technology serve as input for CALPHAD assessment of phase diagrams

36 MATERIALS SCIENCE↗

Form-factor dependence of neutrino-nucleus cross sections using microscopic methods

To achieve its design goals, the next generation of neutrino-oscillation accelerator experiments requires percent-level predictions of neutrino-nucleus cross sections supplemented by robust estimates of the theoretical uncertainties involved. The latter arise from both approximations in solving the nuclear many-body problem and in the determination of the single- and few-nucleon quantities taken as input by many-body methods. To gauge the sensitivity of realistic nuclear many body methods to these few-nucleon quantities, we compute flux-averaged double-differential cross sections using the Green’s function Monte Carlo and spectral function methods as well as different parameterizations of the nucleon axial form factors based on either deuterium bubble-chamber data or lattice quantum chromodynamics calculations. The cross-section results are compared with available experimental data from the MiniBooNE and T2K collaborations. We also discuss the uncertainties associated with N → Δ transition form factors that enter the two-body current operator. We quantify the relations between neutrino-nucleus cross section and nucleon form factor uncertainties. These relations enable us to determine the form factor precision targets required to achieve a given cross-section precision.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Hybrid quantum-classical approach for coupled-cluster Green's function theory

The three key elements of a quantum simulation are state preparation, time evolution, and measurement. While the complexity scaling of time evolution and measurements are well known, many state preparation methods are strongly system-dependent and require prior knowledge of the system's eigenvalue spectrum. Here, we report on a quantum-classical implementation of the coupled-cluster Green's function (CCGF) method, which replaces explicit ground state preparation with the task of applying unitary operators to a simple product state. While our approach is broadly applicable to many models, we demonstrate it here for the Anderson impurity model (AIM). The method requires a number of T gates that grows as O ( N 5 ) per time step to calculate the impurity Green's function in the time domain, where N is the total number of energy levels in the AIM. Since the number of T gates is analogous to the computational time complexity of a classical simulation, we achieve an order of magnitude improvement over a classical CCGF calculation of the same order, which requires O ( N 6 ) computational resources per time step.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Calculation of the spectrum of quasiparticle electron excitations in organic molecular semiconductors

A quasiparticle electronic spectrum belongs to the characteristics of nanoobjects that are most important for applications. The following methods of calculating the electronic spectrum are analyzed: the Kohn-Sham equations of the density functional theory (DFT), the hybrid functional method, the GW approximation, and the Lehmann approximation used in the spectral representation of one-electron Green’s function. The results of these approaches are compared with the data of photoemission measurements of benzene, PTCDA, and phthalocyanine (CuPc, H{sub 2}Pc, FePc, PtPc) molecules, which are typical representatives of organic molecular semiconductors (OMS). This comparison demonstrates that the Kohn-Sham equations of DFT incorrectly reproduce the electronic spectrum of OMS. The hybrid functional method correctly describes the spectrum of the valence and conduction bands; however, the HOMO-LUMO gap width is significantly underestimated. The correct gap width is obtained in both the GW approximation and the Lehmann approach, and the total energy in this approach can be calculated in the local density approximation of DFT.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Method for Describing the Angular Distribution of Optical Radiation Scattered by a Monolayer of Ordered Spherical Particles (Normal Illumination)

We have developed a method for describing the angular distribution of intensity of radiation scattered by a monolayer of homogeneous spatially ordered monodisperse spherical particles normally illuminated by a plane circularly polarized electromagnetic wave. The method is based on the quasicrystalline approximation (QCA) of the theory of multiple scattering of waves (TMSW) using the multipole expansion of fields and the tensor Green function in vectorial spherical wavefunctions. The method is applied for analyzing the characteristics of radiation scattered by a partially ordered monolayer and a monolayer with a nonideal lattice. The results of calculations are compared with the available experimental data on the position of the first-order diffraction peak on the angular and spectral dependences of the intensity of radiation scattered by a closely packed monolayer with a nonideal triangular lattice of SiO{sub 2} particles. Good conformity of the results has been established.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Play Fairway analysis of geothermal resources across the State of Hawai‘i: 4. Updates with new groundwater chemistry, subsurface stress analysis, and focused geophysical surveys

This paper is the fourth in a series on a play fairway analysis of geothermal resources across the State of Hawai‘i. Here we describe recent exploration activities that include groundwater sampling in ten locations statewide, as well as geophysical surveys on Lana‘i, across the SW rift of Haleakala Volcano (Maui), and surrounding Mauna Kea Volcano (Hawai‘i Island). We derive a first-order method for computing topographic stress using Green’s functions, finding that topography can induce appreciable crustal stress. We develop a new method for incorporating depth information about resistivity, density, and topographic stresses into our previously published equations to calculate resource probability and confidence. We incorporate newly collected and modeled data into our calculations to update statewide maps of probability and confidence. Lana‘i Island and southeast Mauna Kea are identified as the top targets for exploratory drilling. The east rift of Haleakala Volcano, the southern-most region of Mauna Loa’s SW rift (Hawai‘i Island), and central Kauai are identified as targets for geophysical surveying.

15 GEOTHERMAL ENERGY↗

Theoretical investigation of spin-filtering in CrAs/GaAs heterostructures

The electronic structure of bulk zinc-blende GaAs, zinc-blende and tetragonal CrAs, and CrAs/GaAs supercells, computed within linear muffin-tin orbital (LMTO) local spin-density functional theory, is used to extract the band alignment for the [1,0,0] GaAs/CrAs interface in dependence of the spin orientation. With the lateral lattice constant fixed to the experimental bulk GaAs value, a local energy minimum is found for a tetragonal CrAs unit cell with a longitudinal ([1,0,0]) lattice constant reduced by ≈2%. Due to the identified spin-dependent band alignment, half-metallicity of CrAs no longer is a key requirement for spin-filtering. Based on these findings, we study the spin-dependent tunneling current in [1,0,0] GaAs/CrAs/GaAs heterostructures within the non-equilibrium Green's function approach for an effective tight-binding Hamiltonian derived from the LMTO electronic structure. Results indicate that these heterostructures are promising candidates for efficient room-temperature all-semiconductor spin-filtering devices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Fully self-consistent finite-temperature $GW$ in Gaussian Bloch orbitals for solids

In this work, we present algorithmic and implementation details for the fully self-consistent finite-temperature $GW$ method in Gaussian Bloch orbitals for solids. Our implementation is based on the finite-temperature Green's function formalism in which all equations are solved on the imaginary axis, without resorting to analytical continuation during the self-consistency. No quasiparticle approximation is employed and all matrix elements of the self-energy are explicitly evaluated. The method is tested by evaluating the band gaps of selected semiconductors and insulators. We show agreement with other, differently formulated, finite-temperature sc ⁢$GW$ implementations when finite-size corrections and basis-set errors are taken into account. By migrating computationally intensive calculations to graphics processing units, we obtain scalable results on large supercomputers with nearly optimal performance. Our work demonstrates the applicability of Gaussian orbital based sc⁢ $GW$ for ab initio correlated material simulations and provides a sound starting point for embedding methods built on top of $GW$.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗