Search NASA⌕ Search

SEARCH · Search NASA

Results for “Greens function”

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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Competing conduction mechanisms in high performance carbon nanotube fibers

The performance of carbon nanotube (CNT) cables, a contender for copper-wire replacement, is tied to its metallic and semi-conducting-like conductivity responses with temperature; the origin of the semi-conducting-like response however is an underappreciated incongruity in literature. With controlled aspect-ratio and doping-degree, over 61 unique cryogenic experiments including anisotropy and Hall measurements, CNT cable performance is explored at extreme temperatures (65 mK) and magnetic-fields (60 T). A semi-conducting-like conductivity response with temperature becomes temperature-independent approaching absolute-zero, uniquely demonstrating for the first time the necessity of heterogeneous fluctuation induced tunneling; complete de-doping leads to localized hopping, contrasting graphite's pure metallic-like response. High-field magneto-resistance (including novel +22 % longitudinal magneto-resistance near room-temperature) is analyzed with hopping and classical two-band models, both yielding a similar parameter useful for conductor development. Varying field-orientation angle uncovers significant two- and four-fold symmetries that are shown to be from Aharonov-Bohm-like corrections to the curvature-induced bandgap, a first for macroscale CNT fibers. Tight-binding calculations using Green's Function formalism model the largest, coherent transport to-date in commensurate CNT bundles in magnetic-field, revealing non-uniform transmission across bundle cross-sections with doping restoring uniformity; independent of doping, transport in bundle-junction-bundle systems are predominantly from CNTs adjacent to the other bundle—demonstrating that smaller bundles are more efficient for electronic transport. The final impact is predicting the ultimate conductivity of heterogeneous CNT cables using temperature and field-dependent transport, surpassing conductivity of traditional metals.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Universal scaling of electron transmission for nearly ballistic and quantum dragon nanodevices

Here, we predict two different universal scaling regimes for the quantum transmission of metallic nanodevices following the addition of a small amount of uncorrelated disorder. A nanodevice is connected to two thin semi-infinite uniform leads, and the Non-Equilibrium Green’s Function (NEGF) methodology yields the electron transmission $\mathscr{T}(E)$ as a function of the injected electron energy $E$. Ballistic nanodevices have no disorder and have $\mathscr{T}(E)=1$for all $E$ that allow electron propagation in the leads. Quantum dragon nanodevices can have extremely strong properly correlated disorder, and still have $\mathscr{T}(E)=1$for all $E$. Additional uncorrelated site disorder leads to Fano resonances in $\mathscr{T}(E)$. Averaging over the uncorrelated disorder we predict using perturbation theory two universal scaling regimes for $\mathscr{T}_{ave}(E)$. The functional form of both universal scaling regimes depend on the device length and width, energy, and variance of the uncorrelated disorder. The second scaling regime, valid for small but somewhat larger uncorrelated disorder than the first scaling regime, also has the form dependent on the density of states of the system. These two scaling regimes are demonstrated to be valid via large scale computer calculations.

Ballistic transport↗
DOE OSTI2021

Full spectrum optical constant interface to the Materials Project

Optical constants characterize the interaction of materials with light and are important properties in material design. Here we present a Python-based Corvus workflow for simulations of full spectrum optical constants from the visible and ultraviolet to hard x-ray wavelengths based on the real-space Green’s function code FEFF10 and structural data from the Materials Project (MP). The Corvus workflow manager and its associated tools provide an interface to FEFF10 and the MP database. The workflow parallelizes the FEFF computations of optical constants over all absorption edges for each material in the MP database specified by a unique MP-ID. The workflow tools determine the distribution of computational resources needed for that case. Similarly, the optical constants for selected sets of materials can be computed in a single-shot. Additionally, to illustrate the approach, we present results for several elemental solids in the periodic table, as well as a sample compound, and compare our predictions with experimental results. In addition, we provide a database of calculated results for all elements for which there is a stable elemental solid at standard conditions available in the Materials Project database. As in x-ray absorption spectra, these results are interpreted in terms of an atomic-like background and fine-structure contributions.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Mobilities of polydisperse hard spheres near a no-slip wall

In the article, we have derived an analytical formulation for far-field hydrodynamic interactions among unequal size hard spheres near a no-slip wall and have implemented the formulation in a Stokesian dynamics model to simulate a suspension of polydisperse particles in a semi-bounded domain. The formulation is based on the multipole expansion of the boundary integral of Stokes flow, and the mobility tensors are deduced from Fáxen’s law together with Green’s function for Stokes flow near a no-slip wall. Lubrication approximation is incorporated to account for the close-distance interactions between any two particles and also between the particles and the wall. The implementation is validated against previous formulations for equal size particles and against a boundary-element code for unequal size particles. The code can be used to simulate the Stokesian or Brownian interaction of unequal particles in presence of a no-slip wall. In the current study, we applied the Stokesian dynamics model to investigate the trajectories of two unequal hard spheres and their redistribution during sedimentation parallel to a wall. We further used it to demonstrate the cases of many-particle sedimentation toward or parallel to a wall.

catalysis (heterogeneous)↗
DOE OSTI2018

Exact diagonalization library for quantum electron models

We present an exact diagonalization C++ template library (EDLib) for solving quantum electron models,including the single-band finite Hubbard cluster and the multi-orbital impurity Anderson model. Theobservables that can be computed using EDLib are single particle Green’s functions and spin–spincorrelation functions. This code provides three different types of Hamiltonian matrix storage that canbe chosen based on the model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Reproducibility in G 0 W 0 calculations for solids

Ab initio many-body perturbation theory within the GW approximation is a Green's function formalism widely used in the calculation of quasiparticle excitation energies of solids. In what has become an increasingly standard approach, Kohn–Sham eigenenergies, generated from a DFT calculation with a strategically-chosen exchange–correlation functional “starting point”, are used to construct G and W, and then perturbatively corrected by the resultant GW self-energy. In practice, there are several ways to construct the GW self-energy, and these can lead to variations in predicted quasiparticle energies. For example, for ZnO and TiO 2 , the GW fundamental gaps reported in the literature can vary by more than 1 eV depending on the GW code used. In this work, we calculate and analyze GW quasiparticle (QP) energies of these and other systems with three different GW codes: BERKELEYGW, ABINIT and YAMBO. Through a systematic analysis of the GW implementation of these three codes, we identify the primary origin of major discrepancies between codes reported in prior literature to be the different implementations the Coulomb divergence in the Fock exchange term and the frequency integration scheme of the GW self-energy. We then eliminate these discrepancies by using common numerical methods and algorithms, demonstrating that the same quasiparticle energies for a given material can be obtained with different codes, within numerical differences ascribable to the technical details of the underling implementations. This work will be important for users and developers in assessing the precision of future GW applications and methods.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Characterization and integration of the singular test integrals in the method-of-moments implementation of the electric-field integral equation

In this paper, we characterize the logarithmic singularities arising in the method of moments from the Green’s function in integrals over the test domain, and we use two approaches for designing geometrically symmetric quadrature rules to integrate these singular integrands. These rules exhibit better convergence properties than quadrature rules for polynomials and, in general, lead to better accuracy with a lower number of quadrature points. In this work, we demonstrate their effectiveness for several examples encountered in both the scalar and vector potentials of the electric-field integral equation (singular, near-singular, and far interactions) as compared to the commonly employed polynomial scheme and the double Ma–Rokhlin–Wandzura (DMRW) rules, whose sample points are located asymmetrically within triangles.

42 ENGINEERING↗
DOE OSTI2020

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↗
DOE OSTI2023

Analytical solution and parameter estimation for heat of wetting and vapor adsorption during spontaneous imbibition in tuff

An analytical expression is derived for the thermal response observed during spontaneous imbibition of water into a dry core of zeolitic tuff. Sample tortuosity, thermal conductivity, and thermal source strength are estimated from fitting an analytical solution to temperature observations during a single laboratory test. Here, the closed-form analytical solution is derived using Green's functions for heat conduction in the limit of “slow” water movement; that is, when advection of thermal energy with the wetting front is negligible. The solution has four free fitting parameters and is efficient for parameter estimation. Laboratory imbibition data used to constrain the model include a time series of the mass of water imbibed, visual location of the wetting front through time, and temperature time series at six locations. The thermal front reached the end of the core hours before the visible wetting front. Thus, the predominant form of heating during imbibition in this zeolitic tuff is due to vapor adsorption in dry zeolitic rock ahead of the wetting front. The separation of the wetting front and thermal front in this zeolitic tuff is significant, compared to wetting front behavior of most materials reported in the literature. This work is the first interpretation of a thermal imbibition response to estimate transport (tortuosity) and thermal properties (including thermal conductivity) from a single laboratory test.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Single-cone finite difference scheme for the (2+1)D Dirac von Neumann equation

An explicit finite difference scheme is presented for the von Neumann equation for (2+1)D Dirac fermions. It is founded upon a staggered space–time grid which ensures a single-cone energy dispersion and performs the time-derivative in one sweep using a three-step leap-frog procedure. It enables a space–time-resolved numerical treatment of the mixed-state dynamics of Dirac fermions within the effective single-particle density matrix formalism. Energy-momentum dispersion, stability and convergence properties are derived. Elementary numerical tests to demonstrate stability properties use parameters which pertain to topological insulator surface states. A method for the simulation of charge injection from an electric contact is presented and tested numerically. Potential extensions of the scheme to a Dirac–Lindblad equation, real-space–time Green's function formulations, and higher-order finite-difference schemes are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

An asymptotic-preserving 2D-2P relativistic Drift-Kinetic-Equation solver for runaway electron simulations in axisymmetric tokamaks

We propose an asymptotic-preserving (AP), uniformly convergent numerical scheme for the relativistic collisional Drift-Kinetic Equation (rDKE) to simulate runaway electrons in axisymmetric toroidal magnetic field geometries typical of tokamak devices. The approach is derived from an exact Green's function solution with numerical approximations of quantifiable impact, and results in a simple, two-step operator-split algorithm, consisting of a collisional Eulerian step, and a Lagrangian orbit-integration step with analytically prescribed kernels. The AP character of the approach is demonstrated by analysis of the dominant numerical errors, as well as by numerical experiments. We demonstrate the ability of the algorithm to provide accurate answers regardless of plasma collisionality on a circular axisymmetric tokamak geometry.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Topological Hall effect in magnetic topological insulator films

Geometric Berry phase can be induced either by spin–orbit coupling, giving rise to the anomalous Hall effect in ferromagnetic materials, or by chiral spin texture, such as skyrmions, leading to the topological Hall effect. Recent experiments have revealed that both phenomena can occur in topological insulator films with magnetic doping, thus providing us with an intriguing platform to study the interplay between these two phenomena. In this work, we report on a numerical simulation of the anomalous Hall and topological Hall effects in a four-band model that can properly describe the quantum well states in the magnetic topological insulator films by combining Landauer–Büttiker formula and the iterative Green’s function method. Our numerical results suggest that spin–orbit coupling in this model plays a different role in the quantum transport in the clean and disordered limits. In the clean limit, spin–orbit coupling mainly influences the longitudinal transport but does not have much effect on topological Hall conductance. In the disordered limit, the longitudinal transport is determined by disorder scattering and spin–orbit coupling is found to affect strongly the topological Hall conductance. Here, this sharp contrast unveils a dramatic interplay between spin–orbit coupling and disorder effect in topological Hall effect in magnetic topological insulator systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Phase-field modeling of dislocation–interstitial interactions

The mechanical behavior of body-centered cubic (BCC) materials can be dramatically affected by the presence of interstitial solute atoms. Here we present a new phase-field dislocation dynamics formulation to include the diffusion of interstitials. Short-range interactions are accounted for by a concentration-dependent lattice energy, and long-range interactions are accounted for by modifications to the elastic energy. The interstitial diffusion law introduces gradients that require methods for minimizing Gibbs oscillations, which is done via a modified Green’s function. The formulation is general to any solute-solvent system and is applied here to Nb-O as a model system, whose interstitial parameters are obtained from ab initio calculations. The effect of O on the core structures of Nb edge and screw dislocations is calculated. The diffusion of O to form interstitial atmospheres around dislocation cores is simulated, as well as the critical stresses required for dislocations to break away or cross slip from these atmospheres. Future applications of the method to simulate complex interstitial embrittlement mechanisms are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Experimental confirmation of first-principles thermal conductivity in Zirconium-doped ThO 2

The degradation of thermal conductivity in advanced nuclear fuels due to the accumulation of fission products and irradiation-induced defects is inevitable, and must be considered as part of safety and efficiency analyses of nuclear reactors. Here, this study examines the thermal conductivity of a zirconium-doped ThO 2 crystal, synthesized via the hydrothermal method using a spatial domain thermoreflectance technique. Zirconium is one of the soluble fission products in oxide fuels that can effectively scatter heat-carrying phonons in the crystalline lattice of fuel. Thus, thermal property measurements of zirconium-doped ThO 2 single crystals provide insights into the effects of substitutional zirconium doping, isolated from extrinsic factors such as grain boundary scattering. The experimental results are compared with first-principles calculations of the lattice thermal conductivity of ThO 2 , employing an iterative solution of the Peierls-Boltzmann transport equation. Additionally, the non-perturbative Green's function methodology is utilized to compute phonon-point defect scattering rates, accounting for local distortions around point defects, including mass difference changes, interatomic force constants, and structural relaxation. The congruence between the predicted results from first-principles calculations and the measured temperature-dependent thermal conductivity validates the computational methodology. Furthermore, the methodologies employed in this study enable systematic investigations of thermal conductivity reduction by fission products, potentially leading to the development of more accurate fuel performance codes.

36 - MATERIALS SCIENCE↗
DOE OSTI2022

Low-temperature thermal properties of Fe-doped Co 3 Sn 2 S 2

Recent studies on Co 3 Sn 2 S 2 usually focus on its electrical and magnetic behaviors as it is a Weyl semimetal, half metal, and anisotropic ferromagnet at the same time. Only a few works looked into its low-temperature thermal characteristics. Here, in this study, we report the low-temperature heat capacity, Seebeck coefficient, and thermal conductivity of Fe-doped Co 3 Sn 2 S 2 , as well as the formation energies, band structures, and spin-polarized density of states calculated by using the spin-polarized relativistic Korringa-Kohn-Rostoker – Green’s function method (sprKKR). The heat capacity data reveal that the sample with a heavier Fe-doping level has a lower Debye temperature and a larger density of states at the Fermi level, the latter confirmed by our ab initio computational results. The Fe-doping increases the absolute value of the Seebeck coefficient at temperatures above 120 K and dramatically suppresses the thermal conductivity. Furthermore, we discuss how the external applied magnetic field affects the Seebeck coefficient and thermal conductivity of pristine Co 3 Sn 2 S 2 . Our findings provide a detailed reference and may stimulate future investigations into the low-temperature thermal transport properties of Co 3 Sn 2 S 2 and other compounds in the shandite family.

36 MATERIALS SCIENCE↗
DOE OSTI2022

CoSyR: A novel beam dynamics code for the modeling of synchrotron radiation effects

The self-consistent nonlinear dynamics of a relativistic charged particle beam interacting with its complete self-fields is a fundamental problem underpinning many of the accelerator design issues in high brightness beam applications, as well as the development of advanced accelerators. Particularly, synchrotron radiation induced effects in a magnetic dispersive beamline element can lead to collective beam instabilities and emittance growth. A novel beam dynamics code is developed based on a Lagrangian method for the calculation of the particles’ radiation near-fields using wavefront/wavelet meshes via the Green’s function of the Maxwell equations. These fields are then interpolated onto a moving mesh for dynamic update of the beam. This method allows radiation co-propagation and self-consistent interaction with the beam in 2D/3D simulations at greatly reduced numerical errors. Multiple levels of parallelisms are inherent in this method and implemented in our code CoSyR to enable at-scale simulations of nonlinear beam dynamics on modern computing platforms using MPI, multi-threading, and GPUs. Here, the current 2D implementation of CoSyR has been used to evaluate the transverse and longitudinal coherent radiation effects on the beam and to investigate beam optics designs proposed for mitigation of beam brightness degradation in a magnetic bunch compressor. In this paper, the design of CoSyR, as well as the benchmark with other coherent synchrotron radiation models, are described and discussed. Extension of the core algorithms to 3D is possible and planned.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Heat flow enhancement in a nanoscale plasmonic junction induced by Kondo resonances and electron-phonon coupling

Here recently, we showed that plasmon-exciton coupling can increase entropy current through a bridge coupled to plasmonic metal nanoparticles. Here we show that electron-phonon coupling can also be used to control the entropy current in similar systems. Entropy current tends to decrease due to electron-phonon coupling and to exhibit a monotonous decrease upon temperature ramping. However, an anomaly affecting the current where it is enhanced by electron-phonon coupling is indicated at around 42 times the system's Kondo temperature. We therefore report means to control heat flow by tuning the Kondo resonance through the electron-phonon coupling. We analyze the conditions that bring about these trends due to electron-phonon coupling by employing non-equilibrium Green's function formulation addressing the entropy current and the derived heat flow.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Analysis of spin frustration in an Fe III 7 cluster using a combination of computational, experimental, and magnetostructural correlation methods

The synthesis, structure, and magnetic properties are reported for [Fe 7 O 3 (O 2 C t Bu) 9 (mda) 3 (H 2 O) 3 ] ( 1 ), where mdaH 2 is N -methyldiethanolamine. 1 was prepared from the reaction of [Fe 3 O(O 2 C t Bu) 6 (H 2 O) 3 ](NO 3 ) with mdaH 2 in a 1:~3 ratio in MeCN. The core of 1 consists of a central octahedral Fe III ion held within a non-planar Fe 6 loop by three μ 3 -O 2- and three μ 2 -RO - arms from the three mda 2- chelates. Variable-temperature dc and ac magnetic susceptibility studies revealed dominant antiferromagnetic coupling, leading to a ground state spin of S = 5 / 2 . The ground state was confirmed by a fit of magnetization data collected in the 0.1–7.0 T and 1.8–10.0 K ranges. The four Fe 2 pairwise exchange parameters ( J 1 - J 4 ) were estimated by independent methods: theoretical calculations using either broken symmetry energy differences (-46.3, -16.2, -3.9, and - 28.1 cm -1 , respectively) or Green’s function approximation methods (-41.4, -14.8, -13.2, and - 24.7 cm -1 ), and a magnetostructural correlation (MSC) previously developed for high nuclearity Fe III /O complexes (-39.5, -13.8, -6.7, and - 23.5 cm -1 ). Additionally, the J 1 - J 4 obtained from the MSC and theoretical methods were used with the program PHI to both simulate χ M T vs T as well as to serve as reasonable input values to fit the experimental data (-41.0, -11.4, -5.0, and - 27.3 cm -1 ). Analysis of the J ij led to identification of the spin frustration effects operative and the resultant spin vector alignments at each Fe III ion, thus allowing for the rationalization of the experimental ground state.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗