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

Nonparametric, data-based kernel interpolation for particle-tracking simulations and kernel density estimation

Traditional interpolation techniques for particle tracking include binning and convolutional formulas that use pre-determined (i.e., closed-form, parameteric) kernels. In many instances, the particles are introduced as point sources in time and space, so the cloud of particles (either in space or time) is a discrete representation of the Green’s function of an underlying PDE. As such, each particle is a sample from the Green’s function; therefore, each particle should be distributed according to the Green’s function. In short, the kernel of a convolutional interpolation of the particle sample “cloud” should be a replica of the cloud itself. This idea gives rise to an iterative method by which the form of the kernel may be discerned in the process of interpolating the Green’s function. When the Green’s function is a density, this method is broadly applicable to interpolating a kernel density estimate based on random data drawn from a single distribution. We formulate and construct the algorithm and demonstrate its ability to perform kernel density estimation of skewed and/or heavy-tailed data including breakthrough curves.

42 ENGINEERING↗
DOE OSTI2022

Ab initio approaches to high-entropy alloys: a comparison of CPA, SQS, and supercell methods

In this work, we present a comparative study of different modeling approaches to the electronic properties of the Hf 0.05 Nb 0.05 Ta 0.8 Ti 0.05 Zr 0.05 high-entropy alloy. Common to our modeling is the methodology to compute the one-particle Green’s function in the framework of density functional theory. We demonstrate that the special quasi-random structures modeling and the supercell, i.e., the locally self-consistent multiple-scattering methods, provide very similar results for the ground state properties such as the spectral function (density of states) and the equilibrium lattice parameter. To reconcile the multiple-scattering single-site coherent potential approximation with the real space supercell methods, we included the effect of screening of the net charges of the alloy components. Based on the analysis of the total energy and spectral functions computed within the density functional theory, we found no signature for the long-range or local magnetic moments formation in the Hf 0.05 Nb 0.05 Ta 0.8 Ti 0.05 Zr 0.05 high-entropy alloy; instead, we find possible superconductivity below ~9K.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Coordinate-space solver for superfluid many-fermion systems with the shifted conjugate-orthogonal conjugate-gradient method

Self-consistent approaches to superfluid many-fermion systems in three dimensions (and their subsequent use in time-dependent studies) require a large number of diagonalizations of very large dimension Hermitian matrices, which results in enormous computational costs. We present an approach based on the shifted conjugate-orthogonal conjugate-gradient (COCG) Krylov method for the evaluation of the Green's function, from which we subsequently extract various densities (particle number, spin, current, kinetic energy, anomalous, etc.) of a nuclear system. The approach eschews the determination of the quasiparticle wave functions and their corresponding quasiparticle energies, which never explicitly appear in the construction of a single-particle Hamiltonian or needed for the calculation of various static nuclear properties, which depend only on densities. As benchmarks this study presents calculations for nuclei with axial symmetry, including the ground state of spherical (magic or semimagic) and axially deformed nuclei, the saddle point in the $^{240}Pu$ constrained fission path, and a vortex in the neutron star crust, and demonstrate the superior efficiency of the shifted COCG Krylov method over traditional approaches.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2015

Designing of spin-filtering devices in zigzag graphene nanoribbons heterojunctions by asymmetric hydrogenation and B-N doping

Using nonequilibrium Green's function in combination with the spin-polarized density functional theory, the spin-dependent transport properties of boron and nitrogen doped zigzag graphene nanoribbons (ZGNRs) heterojunctions with single or double edge-saturated hydrogen have been investigated. Our results show that the perfect spin-filtering effect (100%), rectifying behavior and negative differential resistance can be realized in the ZGNRs-based systems. And the corresponding physical analysis has been given.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Quantum calculations of the carrier mobility: Methodology, Matthiessen's rule, and comparison with semi-classical approaches

We discuss carrier mobilities in the quantum Non-Equilibrium Green's Functions (NEGF) framework. We introduce a method for the extraction of the mobility that is free from contact resistance contamination and with minimal needs for ensemble averages. We focus on silicon thin films as an illustration, although the method can be applied to various materials such as semiconductor nanowires or carbon nanostructures. We then introduce a new paradigm for the definition of the partial mobility μ{sub M} associated with a given elastic scattering mechanism “M,” taking phonons (PH) as a reference (μ{sub M}{sup −1}=μ{sub PH+M}{sup −1}−μ{sub PH}{sup −1}). We argue that this definition makes better sense in a quantum transport framework as it is free from long range interference effects that can appear in purely ballistic calculations. As a matter of fact, these mobilities satisfy Matthiessen's rule for three mechanisms [e.g., surface roughness (SR), remote Coulomb scattering (RCS) and phonons] much better than the usual, single mechanism calculations. We also discuss the problems raised by the long range spatial correlations in the RCS disorder. Finally, we compare semi-classical Kubo-Greenwood (KG) and quantum NEGF calculations. We show that KG and NEGF are in reasonable agreement for phonon and RCS, yet not for SR. We discuss the reasons for these discrepancies.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Coupled Cluster Green's function formulations based on the effective Hamiltonians

In this work, we demonstrate that the effective Hamiltonians obtained with the downfolding procedure based on double unitary coupled cluster (DUCC) ansatz can be used in the context of Green’s function coupled cluster (GFCC) formalism to calculate spectral functions of molecular systems. This combined approach (DUCC-GFCC) provides a significant reduction of numerical effort and good agreement with the corresponding all-orbital GFCC methods in energy windows that are consistent with the choice of active space. These features are demonstrated on the example of two benchmark systems: H 2 O and N 2 , where DUCC-GFCC calculations were performed for active spaces of various sizes.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2018

Sensitivity of Surface Temperature to Oceanic Forcing via q -Flux Green’s Function Experiments. Part I: Linear Response Function

This paper explores the use of the linear response function (LRF) to relate the mean sea surface temperature (SST) response to prescribed ocean heat convergence ( q flux) forcings. Two methods for constructing the LRF based on the fluctuation–dissipation theorem (FDT) and Green’s function (GRF) are examined. A 900-yr preindustrial simulation by the Community Earth System Model coupled with a slab ocean model (CESM–SOM) is used to estimate the LRF using FDT. For GRF, 106 pairs of CESM–SOM simulations with warm and cold q-flux patches are performed. FDT is found to have some skill in estimating the SST response to a q-flux forcing when the local SST response is strong, but it fails in inverse estimation of the q-flux forcing for a given SST pattern. In contrast, GRF is shown to be reasonably accurate in estimating both SST response and q-flux forcing. Possible degradation in FDT may be attributed to insufficient data sampling, significant departure of the SST distribution from Gaussianity, and the nonnormality of the constructed operator. The GRF-based LRF is then used to (i) generate a global surface temperature sensitivity map that shows the q-flux forcing in higher latitudes to be 3–4 times more effective than low latitudes in producing global surface warming, and (ii) identify the most excitable SST mode (neutral vector) that shows marked resemblance to the interdecadal Pacific oscillation (IPO). The latter discovery suggests that the IPO-like fluctuation exists in the absence of the coupling to the ocean dynamics. Coupling to the ocean dynamics in CESM, on the other hand, only enhances the spectral power of the IPO at interannual time scales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Tuning negative differential resistance in single-atomic layer boron-silicon sheets

Using density functional theory and nonequilibrium Green's function formalism for quantum transport calculation, we have quantified the ballistic transport properties along different directions in two-dimensional boron-silicon (B-Si) compounds, as well as the current response to bias voltage. The conductance of the most B-Si devices is higher than the conductance of one-atom-thick boron and silicene. Furthermore, the negative differential resistance phenomenon can be found at certain B-Si stoichiometric composition, and it occurs at various bias voltages. Also, the peak-to-valley ratio is sensitive to the B-Si composition and dependent of the direction considered for B-Si monolayers. The present findings could be helpful for applications of the single-atomic layer B-Si sheets in the field of semiconductor devices or low-dimensional electronic devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Fully-relativistic full-potential multiple scattering theory: A pathology-free scheme

The Green function plays an essential role in the Korringa–Kohn–Rostoker(KKR) multiple scattering method. In practice, it is constructed from the regular and irregular solutions of the local Kohn–Sham equation and robust methods exist for spherical potentials. However, when applied to a non-spherical potential, numerical errors from the irregular solutions give rise to pathological behaviors of the charge density at small radius. Here we present a full-potential implementation of the fully-relativistic KKR method to perform ab initio self-consistent calculation by directly solving the Dirac differential equations using the generalized variable phase (sine and cosine matrices) formalism Liu et al. (2016). The pathology around the origin is completely eliminated by carrying out the energy integration of the single-site Green function along the real axis. Here, by using an efficient pole-searching technique to identify the zeros of the well-behaved Jost matrices, we demonstrated that this scheme is numerically stable and computationally efficient, with speed comparable to the conventional contour energy integration method, while free of the pathology problem of the charge density. As an application, this method is utilized to investigate the crystal structures of polonium and their bulk properties, which is challenging for a conventional real-energy scheme. The noble metals are also calculated, both as a test of our method and to study the relativistic effects.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Carbon doping induced peculiar transport properties of boron nitride nanoribbons p-n junctions

By applying nonequilibrium Green's function combined with density functional theory, we investigate the electronic transport properties of carbon-doped p-n nanojunction based on hexagonal boron nitride armchair nanoribbons. The calculated I-V curves show that both the center and edge doping systems present obvious negative differential resistance (NDR) behavior and excellent rectifying effect. At low positive bias, the edge doping systems possess better NDR performance with larger peak-to-valley ratio (∼10{sup 5}), while at negative bias, the obtained peak-to-valley ratio for both of the edge and center doping systems can reach the order of 10{sup 7}. Meanwhile, center doping systems present better rectifying performance than the edge doping ones, and giant rectification ratio up to 10{sup 6} can be obtained in a wide bias range. These outstanding transport properties are explained by the evolution of the transmission spectra and band structures with applied bias, together with molecular projected self-consistent Hamiltonian eigenvalues and eigenstates.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Process-Independent Effective Coupling from QCD Green's Functions to Phenomenology

This article reports on a very recent proposal for a new type of process-independent QCD effective charge (Binosi et al. in Phys Rev D 96(5):054026, 2017. https://doi.org/10.1103/PhysRevD.96.054026) defined, as an anologue of the Gell-Mann-Low effective charge in QCD, on the ground of nothing but the knowledge of the gauge-field two-point Green's function, albeit modified within a particular computational framework; namely, the combination of pinch technique and background field method which makes possible a systematic rearranging of classes of diagrams in order to redefine the Green's function and have them obey linear QED-like Slavnov-Taylor identities. We have here calculated that effective charge, shown how strikingly well it compares to a process-dependent effective charge based on the Bjorken sum rule; and, finally, employed it in an exploratory calculation of the proton electromagnetic form factor in the hard scattering regime.

confinement↗
NASA NTRS1995

The Schwinger Variational Method

Variational methods have proven invaluable in theoretical physics and chemistry, both for bound state problems and for the study of collision phenomena. The application of the Schwinger variational (SV) method to e-molecule collisions and molecular photoionization has been reviewed previously. The present chapter discusses the implementation of the SV method as applied to e-molecule collisions. Since this is not a review of cross section data, cross sections are presented only to server as illustrative examples. In the SV method, the correct boundary condition is automatically incorporated through the use of Green's function. Thus SV calculations can employ basis functions with arbitrary boundary conditions. The iterative Schwinger method has been used extensively to study molecular photoionization. For e-molecule collisions, it is used at the static exchange level to study elastic scattering and coupled with the distorted wave approximation to study electronically inelastic scattering.

Huo, Winifred M.↗
DOE OSTI2021

Nuclear energy density functionals grounded in ab initio calculations

Here, we discuss the construction of a nuclear energy density functional (EDF) from ab initio computations and advocate the need for a methodical approach that is free from ad hoc assumptions. The equations of state (EoSs) of symmetric nuclear and pure neutron matter are computed using the chiral NNLO sat and the phenomenological AV4' + UIX c Hamiltonians as inputs to self-consistent Green's function (SCGF) and auxiliary field diffusion Monte Carlo (AFDMC) methods. We propose a convenient parametrization of the EoS as a function of the Fermi momentum and fit it on the SCGF and AFDMC calculations. We apply the ab initio based EDF to carry out an analysis of the binding energies and charge radii of different nuclei in the local density approximation. The NNLO sat -based EDF produces encouraging results, whereas the AV4' + UIX c -based one is farther from experiment. Possible explanations of these different behaviors are suggested, and the importance of gradient and spin-orbit terms is analyzed. Our paper paves the way for a practical and systematic way to merge ab initio nuclear theory and density functional theory, while shedding light on some critical aspects of this procedure.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2015

An ab initio study of multiple phonon scattering resonances in silicon germanium alloys

We have computed here phonon scattering rates and density of states in silicon germanium alloys using Green's function calculations and density functional theory. This method contrasts with the virtual crystal approximation (VCA) used in conjunction with Fermi's golden rule, which cannot capture resonance states occurring through the interaction of substitutional impurities with the host lattice. These resonances are demonstrated by density of states and scattering rate calculations in the dilute limit and show broadening as the concentration increases. Although these deviations become significant from the VCA at high frequencies, the relaxation times obtained for these phonon modes are small in both the full scattering theory and the VCA, resulting in their negligible contribution to thermal transport.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Excitons and Davydov splitting in sexithiophene from first-principles many-body Green’s function theory

Organic semiconductors have promising and broad applications in optoelectronics. Understanding their electronic excited states is important to help us control their spectroscopic properties and performance of devices. There have been a large amount of experimental investigations on spectroscopies of organic semiconductors, but theoretical calculation from first principles on this respect is still limited. Here, we use density functional theory (DFT) and many-body Green’s function theory, which includes the GW method and Bethe-Salpeter equation, to study the electronic excited-state properties and spectroscopies of one prototypical organic semiconductor, sexithiophene. The exciton energies of sexithiophene in both the gas and bulk crystalline phases are very sensitive to the exchange-correlation functionals used in DFT for ground-state structure relaxation. We investigated the influence of dynamical screening in the electron-hole interaction on exciton energies, which is found to be very pronounced for triplet excitons and has to be taken into account in first principles calculations. In the sexithiophene single crystal, the energy of the lowest triplet exciton is close to half the energy of the lowest singlet one. While lower-energy singlet and triplet excitons are intramolecular Frenkel excitons, higher-energy excitons are of intermolecular charge-transfer type. The calculated optical absorption spectra and Davydov splitting are in good agreement with experiments.

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

First-principles calculations of thermoelectric properties of TiN/MgO superlattices: The route for an enhancement of thermoelectric effects in artificial nanostructures

We present the thermoelectric properties of TiN/MgO superlattices employing first-principles calculation techniques. The Seebeck coefficients, the electrical conductances, the thermal conductances, and the figure of merit are investigated employing electrical and thermal transport calculations based on density functional theory combined with the nonequilibrium Green's function and nonequilibrium molecular dynamics simulation methods. The TiN/MgO superlattices with a small lattice mismatch at the interfaces are ideal systems to study the way for an enhancement of thermoelectric properties in artificial nanostructures. We find that the interfacial scattering between the two materials in the metal/insulator superlattices causes the electrical conductance to change rapidly, which enhances the Seebeck coefficient significantly. We show that the figure of merit for the artificial superlattice nanostructures has a much larger value compared with that of the bulk material and changes drastically with the superlattice configurations at the atomistic level.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

A mixed, unified forward/inverse framework for earthquake problems: fault implementation and coseismic slip estimate

SUMMARY We introduce a new finite-element (FE) based computational framework to solve forward and inverse elastic deformation problems for earthquake faulting via the adjoint method. Based on two advanced computational libraries, FEniCS and hIPPYlib for the forward and inverse problems, respectively, this framework is flexible, transparent and easily extensible. We represent a fault discontinuity through a mixed FE elasticity formulation, which approximates the stress with higher order accuracy and exposes the prescribed slip explicitly in the variational form without using conventional split node and decomposition discrete approaches. This also allows the first order optimality condition, that is the vanishing of the gradient, to be expressed in continuous form, which leads to consistent discretizations of all field variables, including the slip. We show comparisons with the standard, pure displacement formulation and a model containing an in-plane mode II crack, whose slip is prescribed via the split node technique. We demonstrate the potential of this new computational framework by performing a linear coseismic slip inversion through adjoint-based optimization methods, without requiring computation of elastic Green’s functions. Specifically, we consider a penalized least squares formulation, which in a Bayesian setting—under the assumption of Gaussian noise and prior—reflects the negative log of the posterior distribution. The comparison of the inversion results with a standard, linear inverse theory approach based on Okada’s solutions shows analogous results. Preliminary uncertainties are estimated via eigenvalue analysis of the Hessian of the penalized least squares objective function. Our implementation is fully open-source and Jupyter notebooks to reproduce our results are provided. The extension to a fully Bayesian framework for detailed uncertainty quantification and non-linear inversions, including for heterogeneous media earthquake problems, will be analysed in a forthcoming paper.

58 GEOSCIENCES↗