Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “GREEN 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 1,027 records · Page 57

Molecular Optoelectronics (Final Technical Report)

The project developed new theoretical tools for analysis of nanodevices subjected to external optical and bias driving. Traditional theoretical techniques for molecular optoelectronics are either based on kinetic schemes or utilize quasiparticle description. The former fails to account for quantum mechanics of the system-bath interface and may fail qualitatively when used in open interacting systems. The latter is inconvenient in (usual for molecular optical spectroscopy) many body states description of the system. We developed a new technique - the Hubbard nonequilibrium Green's functions method - which allows to avoid shortcomings of standard theoretical approaches. The Hubbard NEGF was applied to analysis of experimental measurements of photoinduced current and bias-induced electroluminescence in STM molecular junctions. We also demonstrated practical usefulness of the Hubbard NEGF within newly proposed universal theory of current-induced forces for nonadiabatic molecular dynamics which generalizes celebrated Head-Gordon and Tully expression for electronic friction. Finally, following latest trends in laser techniques we studied local responses (local fluxes and local noise spectroscopy) of biased junctions as a source of information not accessible by total responses (total fluxes and noises) of the system.

30 DIRECT ENERGY CONVERSION↗

Interfacial contributions to spin-orbit torque and magnetoresistance in ferromagnet/heavy-metal bilayers

The thickness dependence of spin-orbit torque and magnetoresistance in ferromagnet/heavy-metal bilayers is studied here using the first-principles nonequilibrium Green's function formalism combined with the Anderson disorder model. A systematic expansion in orthogonal vector spherical harmonics is used for the angular dependence of the torque. The dampinglike torque in Co/Pt and Co/Au bilayers can be described as a sum of the spin-Hall contribution, which increases with thickness in agreement with the spin-diffusion model, and a comparable interfacial contribution. The magnetoconductance in the plane perpendicular to the current in Co/Pt bilayers is of the order of a conductance quantum per interfacial atom, exceeding the prediction of the spin-Hall model by more than an order of magnitude. This indicates that the “spin-Hall magnetoresistance,” similarly to the dampinglike torque, has a large interfacial contribution unrelated to the spin-Hall effect.

36 MATERIALS SCIENCE↗

Analytical solution for partial heating on the exterior of the pipe with application to measuring fluid flow rate

An analytical model using Green’s functions for partial external heating of a pipe is developed, which results in an exact mathematical solution for the radial and axial temperature distribution in the pipe wall. Partial heating consists of a constant heat flux function imposed over a small section of the exterior of a pipe and for a limited time duration. The solution comprises steady-state and transient parts, and an algebraic identity is used to decrease the number of summation terms in the slowly-converging steady- state part. Intrinsic verification principles are used to verify the solution. As an example application, this transient solution is applied toward the development of a simple, noninvasive method for in-field measurement of the flow rate in pipes. Furthermore, to simulate this application, a pulse of energy is imposed to the wall of the pipe, and the developed mathematical solution is used to find the flow rate of the fluid inside the pipe. An optimal experiment is designed to find the best measurement location and time.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

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↗

Strong anharmonicity dictates ultralow thermal conductivities of type-I clathrates

Type-I clathrate solids have attracted significant interest due to their ultralow thermal conductivities and sub- sequent promise for thermoelectric applications, yet the mechanisms underlying these properties are not well understood. Here, we extend the framework of vibrational dynamical mean-field theory (VDMFT) to calculate temperature-dependent thermal transport properties of solids using a many-body Green’s function approach. When applied to a coarse-grained description of 𝑋 8 Ga 16 Ge 30 , where 𝑋= Ba, Sr, we find that nonresonant scattering between cage acoustic modes and rattling modes leads to a reduction of acoustic phonon lifetimes and thus thermal conductivities. Moreover, we find that the moderate temperature dependence of conductivities above 300 K, which is consistent with experimental measurements, cannot be reproduced by textbook perturbation theory calculations, which predict a 𝑇 −1 dependence. Therefore, we suggest that nonperturbative anharmonic effects, including four- and higher-phonon scattering processes, are responsible for the ultralow thermal conductivities of type-I clathrates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Approximation of periodic Green's operator in real space using numerical integration and its use in fast Fourier transform-based micromechanical models

In this work, we propose an expression for the periodic first derivative of Green's function in real space. The proposed expression allows an alternative way of computing the periodic Green's operator based on periodically summing the free-space Green's operator in terms of an appropriate quadrature rule. We provide computational examples, which show the accuracy of the proposed approach, together with reduced spurious oscillations in the solution fields.

42 ENGINEERING↗

Analysis of direct and indirect noise in a next-generation aviation gas turbine combustor

A large-eddy simulation (LES) of a next-generation combustor is performed to examine effects of this combustor concepts on direct and indirect combustion noise characteristics. Direct noise is computed considering the unsteady heat release predictions while entropy fluctuations in the downstream part of the combustor are used to estimate indirect noise through the transfer function of the outlet nozzle. A low-order acoustic reconstruction technique, which utilizes the Green’s function of the configuration, is developed to compute the acoustics inside the combustor. Comparisons to experimental results are reported, showing the reasonable accuracy of the LES and combustion noise computations. The direct noise spectrum peaks at frequencies around 3 kHz because of the fast timescales associated with the chemical reactions inside the combustor. In addition, the indirect noise is dominant at low frequencies, i.e., less than 400 Hz, because of the characteristics of the entropy spectrum at the outlet nozzle. Further, compared to a conventional rich-quench-lean (RQL) combustor configuration, significant differences in the acoustics are observed, and are explained by two specific technological novelties. Firstly, the direct noise spectrum peaks at significantly higher frequencies due to the combustor’s compactness. Secondly, the entropy inhomogeneties downstream of the combustion region are an order of magnitude smaller in amplitude due to the lack of dilution with cold air. This leads to a significant reduction in the indirect noise amplitude at low frequencies, where it is the dominant source of noise. These results suggest opportunities for advanced combustion technologies to achieve substantial reductions in combustion-related noise.

42 ENGINEERING↗

Gauge constrained algorithm of variational discrete action theory at N = 3 for the multiorbital Hubbard model

The recently developed variational discrete action theory (VDAT) provides a systematic variational approach to the ground state of the quantum many-body problem, where the quality of the solution is controlled by an integer N, and increasing N monotonically approaches the exact solution. VDAT can be exactly evaluated in the d = ∞ multiorbital Hubbard model using the self-consistent canonical discrete action theory (SCDA), which requires a self-consistency condition for the integer time Green's functions. Previous work demonstrates that N = 3 accurately captures multiorbital Mott/Hund physics at a cost similar to the Gutzwiller approximation. Here we employ a gauge constraint to automatically satisfy the self-consistency condition of the SCDA at N = 3, yielding an even more efficient algorithm with enhanced numerical stability. We derive closed form expressions of the gauge constrained algorithm for the multiorbital Hubbard model with general density-density interactions, allowing VDAT at N = 3 to be straightforwardly applied to the seven-orbital Hubbard model. We present results and a performance analysis using N = 2 and N = 3 for the SU⁡(2⁢N orb ) Hubbard model in d = ∞ with N orb = 2–8, and compare to numerically exact dynamical mean-field theory solutions where available. Finally, the developments in this work will greatly facilitate the application of VDAT at N = 3 to strongly correlated electron materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Theory of molecular conductance using a modular approach

This paper probes the correlation between the conductance of a molecular wire (the property of a whole system) and its constituent backbone units (modules). By using a tight-binding Hamiltonian combined with single-particle Green’s functions, we develop an approach that enables an estimate of a conductance decay constant in terms of the Hamiltonians of molecular backbone units and the couplings between two nearest-neighbor units in the off-resonant tunneling regime. For demonstration, we examine several representative molecular systems in a framework of the Hückel model (the simplest atomistic-level model). The Hückel model can be reduced to a single-orbital-per-site formulation [A. Nitzan, Annu. Rev. Phys. Chem. 52, 681 (2001)], and each energy level in the single-orbital-per-site picture can be expressed in an explicit form including the synergistic effect of all molecular orbitals of a molecular backbone unit. Finally, based on the proposed approach, we show the correspondence between the complete destructive quantum interference and an infinite injection gap and derive the preconditions of the modified Simmons equation and the rule of intramolecular series circuits.

74 ATOMIC AND MOLECULAR PHYSICS↗

Scanning tunneling shot-noise spectroscopy in Kondo systems

Using a large- N theory in combination with the Keldysh nonequilibrium Green's function formalism, we investigate the current, differential conductance, zero-frequency shot noise, and Fano factor as measured by scanning tunneling shot-noise spectroscopy (STSNS) using a scanning tunneling microscope (STM) near single Kondo impurities and in Kondo lattices. We show that the Fano factor F exhibits a characteristic bias dependence arising from Kondo screening that is similar to the Kondo resonance observed in the differential conductance. Moreover, the line shape of F is strongly dependent on the ratio of the tunneling amplitudes for electron tunneling from the STM tip into the conduction band and electronic levels of the magnetic adatoms. We demonstrate that the Fano factor can be enhanced or suppressed due to interference effects and as such, is not only a sensitive probe for the correlation effects arising from Kondo screening, but also for quantum interference between tunneling electrons. We identify a correlation between the form of the differential conductance and the Fano factor that could be tested in future STSNS experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Anharmonic phonon behavior via irreducible derivatives: Self-consistent perturbation theory and molecular dynamics

Cubic phonon interactions are now regularly computed from first principles, and the quartic interactions have begun to receive more attention. Given this realistic anharmonic vibrational Hamiltonian, the classical phonon Green's function can be precisely measured using molecular dynamics, which can then be used to rigorously assess the range of validity for self-consistent diagrammatic approaches in the classical limit. Here we use the bundled irreducible derivative approach to efficiently and precisely compute the cubic and quartic phonon interactions of CaF 2 , systematically obtaining the vibrational Hamiltonian purely in terms of irreducible derivatives. We demonstrate that the 4-phonon sunset diagram has an important contribution to the optical phonon linewidths beyond T = 500 K. Reasonable results are obtained even at T = 900 K when performing self-consistency using the 4-phonon loop diagram and evaluating the 3-phonon bubble and 4-phonon sunset diagrams post-self-consistency. Further improvements are obtained by performing quasiparticle perturbation theory, where both the 4-phonon loop and the real part of the 3-phonon bubble are employed during self-consistency. Finally, our irreducible derivative approach to self-consistent perturbation theory is a robust tool for studying anharmonic phonons in both the quantum and classical regimes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Advanced calculations of X-ray spectroscopies with FEFF10 and Corvus

The real-space Green's function code FEFF has been extensively developed and used for calculations of X-ray and related spectra, including X-ray absorption (XAS), X-ray emission (XES), inelastic X-ray scattering, and electron energy-loss spectra. The code is particularly useful for the analysis and interpretation of the XAS fine-structure (EXAFS) and the near-edge structure (XANES) in materials throughout the periodic table. Nevertheless, many applications, such as non-equilibrium systems, and the analysis of ultra-fast pump–probe experiments, require extensions of the code including finite-temperature and auxiliary calculations of structure and vibrational properties. To enable these extensions, we have developed in tandem a new version FEFF10 and new FEFF -based workflows for the Corvus workflow manager, which allow users to easily augment the capabilities of FEFF10 via auxiliary codes. This coupling facilitates simplified input and automated calculations of spectra based on advanced theoretical techniques. The approach is illustrated with examples of high-temperature behavior, vibrational properties, many-body excitations in XAS, super-heavy materials, and fits of calculated spectra to experiment.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Effects of defects and dephasing on charge and spin currents in two-dimensional topological insulators

Using the nonequilibrium Keldysh Green's function formalism, we investigate the effect of defects on the electronic structure and transport properties of two-dimensional topological insulators (TI). We demonstrate how the spatial flow of charge changes between the topologically protected edge and bulk states and show that elastically and inelastically scattering defects that preserve the time-reversal symmetry of the TI lead to qualitatively different effects on the TI's local electronic structure and its transport properties. Moreover, we show that the recently predicted ability to create highly spin-polarized currents by breaking the time-reversal symmetry of the TI via magnetic defects [J. S. Van Dyke and D. K. Morr, Phys. Rev. B 93, 081401 (2016)] is robust against the inclusion of a Rashba spin-orbit interaction and the effects of dephasing, and remains unaffected by changes over a wide range of the TI's parameters. We discuss how the sign of the induced spin currents changes under symmetry operations, such as reversal of bias and gate voltages, or spatial reflections. Finally, we show that the insight into the interplay between topology and symmetry of the magnetic defects can be employed for the creation of intriguing quantum phenomena, such as highly localized magnetic fields inside the TI.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Ab initio self-consistent many-body theory of polarons at all couplings

We present a theoretical framework to describe polarons from first principles within a many-body Green’s function formalism. Starting from a general electron-phonon Hamiltonian, we derive a self-consistent Dyson equation in which the phonon-mediated self-energy is composed by two distinct terms. One term is the Fan-Migdal self-energy and describes dynamic electron-phonon processes, the other term is a contribution to the self-energy originating from the static displacements of the atomic nuclei in the polaronic ground state. The lowest-order approximation to the present theory yields the standard many-body perturbation theory approach to electron-phonon interactions in the limit of large polarons, and the ab initio polaron equations introduced in the limit of small polarons. Here, a practical recipe to implement the present unifying formalism in first-principles calculations is outlined. We apply our method to the Fröhlich model, and obtain remarkably accurate polaron energies at all couplings, in line with Feynman’s polaron theory and diagrammatic Monte Carlo calculations. We also recover the well-known results of Fröhlich and Pekar at weak and strong coupling, respectively. The present approach enables predictive many-body calculations of polarons in real materials at all couplings.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Direct and indirect spin current generation and spin-orbit torques in ferromagnet/nonmagnet/ferromagnet trilayers

Spin-orbit torques in ferromagnet/nonmagnet/ferromagnet trilayers are studied using a combination of symmetry analysis, circuit theory, semiclassical simulations, and first-principles calculations using the nonequilibrium Green's function method with supercell disorder averaging. Here, we focus on unconventional processes involving the interplay between the two ferromagnetic layers, which are classified into direct and indirect mechanisms. The direct mechanism involves spin current generation by one ferromagnetic layer and its subsequent absorption by the other. In the indirect mechanism, the in-plane spin-polarized current from one ferromagnetic layer “leaks” into the other layer, where it is converted into an out-of-plane spin current and reabsorbed by the original layer. The direct mechanism results in a predominantly dampinglike torque, which damps the magnetization towards a certain direction 𝐬 𝑑 . The indirect mechanism results in a predominantly fieldlike torque with respect to a generally different direction 𝐬 𝑓 . Similarly to the current-in-plane giant magnetoresistance, the indirect mechanism is only active if the thickness of the nonmagnetic spacer is smaller than or comparable to the mean free path. Numerical calculations for a semiclassical model based on the Boltzmann equation confirm the presence of both direct and indirect mechanisms of spin current generation. First-principles calculations reveal sizable unconventional spin-orbit torques in Co/Cu/Co, Py/Cu/Py, and Co/Pt/Co trilayers and provide strong evidence of indirect spin current generation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Quantifying EGS Reservoir Complexity with an Integrated Geophysical Approach-Improved Resolution Ambient Seismic Noise Interferometry

The primary goal of the project was met, as improved resolution (from kilometers to tens/hundreds of meters and to a depth of at least 3 km) was added to the seismic interferometry-estimated parameters relative to the work reported in AltaRock Energy Inc. (2014) for the Dixie Valley Geothermal Wellfield. A procedure for initial fault identification, using seismic interferometry, was developed. The project team further improved, and applied a new, cost-effective, and non-invasive exploration method, developed and validated since 2010, during two previous experiments, at Dixie Valley and Soda Lake, Nevada. The method, initially designed for hydrothermal exploration, was adapted to EGS exploration conditions, which require deeper crustal information. Regions of elevated energy for both the 2-7 Hz and 5-15 Hz frequency bands are located above intersections of multiple faults. Thus, an unprecedented, innovative seismic field deployment and ambient seismic noise and signal processing were successfully used to estimate crustal structure and to identify fault locations at depths up to 3 km, at a cost per kilometer of survey that in some cases may be as much as ten times less than the cost of a controlled-source survey. Based on valid and useful seismic data collected and analyzed during Phase I, in Phase II the team has further quantified statistically the advantage of using seismic analysis, in combination with other geological and geophysical techniques, for EGS target identification. A power spectral density analysis of the passive-seismic empirical Green’s functions reveals significantly higher power recorded over the exploration area at the power plant, and in the drilling area, near the BB station 10201, especially on the horizontal components. In this study, the most important relationship observed during the prior report, namely the one between velocity and temperature was explored. The areas that show anomalously high predicted temperature in the geostatistics occur where multiple faults intersect. Not surprisingly, such areas span between the Dixie Valley range front fault and the basin-ward Piedmont faults, in the vicinity of numerous wells. These locations also fall in the area of high EGS favorability from the prior 2014 study. The anomalies also coincide with the high power spectral density values. The result of altering all favorability criteria for the new data reached a maximum of 1/8th of a point and were not large enough to change any favorability value reported in the 2014 study. The new data reinforce the results obtained from the previous 2014 study. The current high-resolution ambient seismic noise project provided excellent results in defining faults in the wellfield. As such the identification of known faults by the project demonstrated the utility of this method for identifying faults at an exploration stage in geothermal development. Such fault identification, while extremely useful in exploration investigation, has not added sufficient new data for the generation of significantly different DVGW favorability maps in the current project. We believe the cost savings of our passive-seismic technique for EGS, over common industry active-source survey techniques, may have a minimum of 30%, and a maximum of 90%. There are many significant non-cost benefits of passive seismic surveys. Our passive-source method provides an environmentally-friendly seismic survey with minimal to negligible impacts. The economic benefit is greater EGS project feasibility.

15 GEOTHERMAL ENERGY↗

Atomistic modeling of the metallic-to-semiconducting phase boundaries in monolayer MoS{sub 2}

Recent experimental demonstration on the coexistence of metallic and semiconducting phases in the same monolayer MoS{sub 2} crystal has attracted much attention for its use in ultra-low contact resistance-MoS{sub 2} transistors. However, the electronic structures of the metallic-to-semiconducting phase boundaries, which appear to dictate the carrier injection in such transistors, are not yet well understood. In this letter, interfacing the 2H and 1T′ polytypes appropriately, we first model the “beta” (β) and the “gamma” (γ) phase boundaries, and demonstrate good agreement with experiential results. We then apply first-principles based density functional theory to calculate the electronic structures for those optimized geometries. We further employ non equilibrium Green's function formalism to evaluate the transmission spectra and the local density of states (LDOS) in order to assess the Schottky barrier nature of the phase boundaries. Our study reveals that while the γ boundary yields p-type Schottky barrier, the β boundary leads to the distinct symmetric Schottky barrier with an atomically sharp transition region. This understanding could be useful for designing high performance transistors using phase-engineered MoS{sub 2} crystals.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Toward ab-initio nuclear theory calculations of 𝛿 𝐶

We propose a theory framework to study the isospin-symmetry breaking correction 𝛿 C in superallowed nuclear 𝛽 decays, crucial for the precise determination of |𝑉 𝑢⁢𝑑 |. Based on a general assumptions of the isovector dominance in isospin-symmetry breaking interactions, we construct a set of functions 𝐹 𝑇 𝑧 which involve nuclear matrix elements of isovector monopole operators and the nuclear Green's function. Via the functions 𝐹 𝑇 𝑧 , a connection of 𝛿 C to measurable electroweak nuclear radii is established, providing an experimental gauge of the theory accuracy of 𝛿 C . Here, we outline a strategy to perform ab initio calculations of 𝐹 𝑇 𝑧 based on the Lanczos algorithm, and discuss its similarity with other nuclear-structure-dependent inputs in nuclear 𝛽 decays.

Ab initio calculations↗