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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

A theoretical and numerical study of polarimetric scattering and emission from random rough surfaces with anisotropic directional spectrum

In this paper, theoretical and numerical results of the polarimetric scattering and emission from random rough surfaces with anisotropic directional spectrum are presented for the remote sensing of ocean and soil surfaces. The polarimetric scattered field for rough dielectric surfaces is derived to the second order by the small perturbation method (SPM). It is found that the second-order scattered field is coherent in nature, and its coefficients for different polarizations present the lowest-order corrections to the Fresnel reflection coefficients of the surfaces. In addition, the cross-polarized (HV and VH) components of the coherent fields are reciprocal and not zero for surfaces with anisotropic directional spectrum when the azimuth angle of the incident direction is not aligned with the symmetry directions of surfaces. In order to verify the energy conservation condition of the theoretical results, which is important if the theory is to be applied to the passive polarimetry of rough surfaces, a Monte Carlo simulation is performed to numerically calculate the polarimetric reflectivities of one-dimensional random rough surfaces which are generated with a prescribed power-law spectrum in the spectral domain and transformed to the spatial domain by the FFT. The surfaces simulated by this approach are periodic with the period corresponding to the low-wavenumber cutoff. To calculate the scattering from periodic dielectric surfaces, the authors present a new numerical technique which applies the Floquet theorem to reduce the problem to one period and does not require the evaluation of one-dimensional periodic Green's function used in the conventional method of moment formulation. Once the scattering coefficients are obtained, the polarimetric Stokes vectors for the emission from the random surfaces are then calculated according to the Kirchhoff's law and are illustrated as functions of relative azimuth observation and row directions. The second-order SPM is also reduced to this one-dimensional case, and theoretical results for all four Stokes parameters are in excellent agreement with the numerical results obtained from the Monte Carlo simulation. Moreover, the second order coherent fields are found to be necessary for the theoretical evaluation of the third and fourth Stokes parameters. Without that, the reflectivities of random rough surfaces are significantly over-estimated, and the signs of the third and the fourth Stokes parameters are found to be incorrect, if calculated by using the first-order incoherent fields alone.

Yueh, S. H.↗
DOE OSTI2014

Rectification induced in N{sub 2}{sup AA}-doped armchair graphene nanoribbon device

By using non-equilibrium Green function formalism in combination with density functional theory, we investigated the electronic transport properties of armchair graphene nanoribbon devices in which one lead is undoped and the other is N{sub 2}{sup AA}-doped with two quasi-adjacent substitutional nitrogen atoms incorporating pairs of neighboring carbon atoms in the same sublattice A. Two kinds of N{sub 2}{sup AA}-doped style are considered, for N dopants substitute the center or the edge carbon atoms. Our results show that the rectification behavior with a large rectifying ratio can be found in these devices and the rectifying characteristics can be modulated by changing the width of graphene nanoribbons or the position of the N{sub 2}{sup AA} dopant. The mechanisms are revealed to explain the rectifying behaviors.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2014

Negative differential resistance devices by using N-doped graphene nanoribbons

Recently, extensive efforts have been devoted to the investigations of negative differential resistance (NDR) behavior in graphene. Here, by performing fully self-consistent density functional theory calculations combined with non-equilibrium Green's function technique, we investigate the transport properties of three molecules from conjugated molecule, one-dimension alkane chain, and single molecule magnet, which are sandwiched between two N-doped zigzag and armchair graphene nanoribbons (GNRs). We observe robust NDR effect in all examined molecular junctions including benzene, alkane, and planar four-coordinated Fe complex. Through the analyses of the calculated electronic structures and the bias-dependent transmission coefficients, we find that the narrow density of states of N-doped GNRs and the bias-dependent effective coupling between the discrete frontier molecular orbitals and the subbands of N-doped GNRs are responsible for the observed NDR phenomenon. These theoretical findings imply that N-doped GNRs hold great potential for building NDR devices based on various molecules.

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

A WENO-solver combined with adaptive momentum discretization for the Wigner transport equation and its application to resonant tunneling diodes

We present a novel numerical scheme for the deterministic solution of the Wigner transport equation, especially suited to deal with situations in which strong quantum effects are present. The unique feature of the algorithm is the expansion of the Wigner function in local basis functions, similar to finite element or finite volume methods. This procedure yields a discretization of the pseudo-differential operator that conserves the particle density on arbitrarily chosen grids. The high flexibility in refining the grid spacing together with the weighted essentially non-oscillatory (WENO) scheme for the advection term allows for an accurate and well-resolved simulation of the phase space dynamics. A resonant tunneling diode is considered as test case and a detailed convergence study is given by comparing the results to a non-equilibrium Green's functions calculation. The impact of the considered domain size and of the grid spacing is analyzed. The obtained convergence of the results towards a quasi-exact agreement of the steady state Wigner and Green's functions computations demonstrates the accuracy of the scheme, as well as the high flexibility to adjust to different physical situations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

A new Green's function Monte Carlo algorithm for the solution of the two-dimensional nonlinear Poisson–Boltzmann equation: Application to the modeling of the communication breakdown problem in space vehicles during re-entry

The objective of this paper is the exposition of a recently-developed, novel Green's function Monte Carlo (GFMC) algorithm for the solution of nonlinear partial differential equations and its application to the modeling of the plasma sheath region around a cylindrical conducting object, carrying a potential and moving at low speeds through an otherwise neutral medium. The plasma sheath is modeled in equilibrium through the GFMC solution of the nonlinear Poisson–Boltzmann (NPB) equation. The traditional Monte Carlo based approaches for the solution of nonlinear equations are iterative in nature, involving branching stochastic processes which are used to calculate linear functionals of the solution of nonlinear integral equations. Over the last several years, one of the authors of this paper, K. Chatterjee has been developing a philosophically-different approach, where the linearization of the equation of interest is not required and hence there is no need for iteration and the simulation of branching processes. Instead, an approximate expression for the Green's function is obtained using perturbation theory, which is used to formulate the random walk equations within the problem sub-domains where the random walker makes its walks. However, as a trade-off, the dimensions of these sub-domains have to be restricted by the limitations imposed by perturbation theory. The greatest advantage of this approach is the ease and simplicity of parallelization stemming from the lack of the need for iteration, as a result of which the parallelization procedure is identical to the parallelization procedure for the GFMC solution of a linear problem. The application area of interest is in the modeling of the communication breakdown problem during a space vehicle's re-entry into the atmosphere. However, additional application areas are being explored in the modeling of electromagnetic propagation through the atmosphere/ionosphere in UHF/GPS applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Improving Spectral Resolution from Real-time Evolution for Correlated Systems

Abstract The quality of numerically simulated spectra using real-time evolution methods for strongly correlated systems is affected by both the length of simulation time and the system size, limiting resolution in both frequency and momentum. In this work, we propose a computationally cheap, linear autoregressive machine learning-based framework to extend short-time and short-distance results over a wider range. We use the proposed method to extend the lesser Green’s function for both the Hubbard model and the much more computationally challenging Hubbard-extended Holstein model. This technique significantly improves both the frequency and momentum resolution of the single-particle removal spectrum $${\mathcal{A}}(k,\omega )$$ A ( k , ω ) , allowing the observation of otherwise obscured spectral features due to electron-phonon coupling.

Tang, Ta↗
DOE OSTI2021

Estimation of the orientation of stress in the Earth’s crust without earthquake or borehole data

Abstract Mechanical stress acting in the Earth’s crust is a fundamental property that is important for a wide range of scientific and engineering applications. The orientation of maximum horizontal compressive stress can be estimated by inverting earthquake source mechanisms and measured directly from borehole-based measurements, but large regions of the continents have few or no observations. Here we present an approach to determine the orientation of maximum horizontal compressive stress by measuring stress-induced anisotropy of nonlinear susceptibility, which is the derivative of elastic modulus with respect to strain. Laboratory and Earth experiments show that nonlinear susceptibility is azimuthally dependent in an anisotropic stress field and is maximum in the orientation of maximum horizontal compressive stress. We observe this behavior in the Earth—in Oklahoma and New Mexico, U.S.A, where maximum nonlinear susceptibility coincides with the orientation of maximum horizontal compressive stress measured using traditional methods. Our measurements use empirical Green’s functions and solid-earth tides and can be applied at different temporal and spatial scales.

58 GEOSCIENCES↗
DOE OSTI2013

Microsegregation effects on the thermal conductivity of silicon-germanium alloys

A silicon-germanium (SiGe) alloy is a promising candidate for thermoelectric materials; while it shows a significantly reduced thermal conductivity (κ) as compared to pure Si and Ge, the κ values obtained from previous experiments and computations tend to be widely scattered. We present here a computational analysis of thermal transport in SiGe, particularly the effects of the local segregation (microsegregation) of alloying elements. Our nonequilibrium molecular dynamics simulations confirm the strong dependence of κ on the Si:Ge ratio and the occurrence of the minimum κ around Si{sub 0.8}Ge{sub 0.2}, consistent with existing experimental observations. Moreover, our study clearly demonstrates that the κ of Si{sub 0.8}Ge{sub 0.2} increases substantially and monotonically as Ge atoms undergo segregation; that is, the magnitude of alloy scattering is found to be sensitive to homogeneity in the distribution of alloying elements. Nonequilibrium Green's function analysis also shows that such microsegregation enhances phonon transmission due to the reduced number of scattering centers. The findings highlight that distribution homogeneity, along with composition, can be a critical factor in determining the κ of SiGe alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2015

First-order exchange and self-energy corrections to static density correlation function of a spin-polarized two-dimensional quantum electron fluid

Static density-density correlation function has been calculated for a spin-polarized two-dimensional quantum electron fluid by including the first-order exchange and self-energy corrections to the random-phase approximation (RPA). This is achieved by determining these corrections to the RPA linear density-density response function, obtained by solving the equation of motion for the single-particle Green’s function. Resulting infinite hierarchy of equations (involving higher-order Green’s functions) is truncated by factorizing the two-particle Green’s function as a product of the single-particle Green’s function and one-particle distribution function. Numerical results of correlation function are compared directly against the quantum Monte Carlo simulation data due to Tanatar and Ceperley for different coupling parameter (r{sub s}) values. We find almost exact agreement for r{sub s}=1, with a noticeable improvement over the RPA. Its quality, however, deteriorates with increasing r{sub s}, but correction to RPA is quite significant.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Application of a self-consistent NEGF procedure to study the coherent transport with phase breaking scattering in low dimensional systems

We have studied Quantum Transport with dephasing in Low Dimensional systems. Here, we apply a self-consistent NEGF procedure to study the transport mechanism in low-dimensional systems with phase breaking scatterers. Under this we have determined the transmission coefficient of a very small Multi-Moded Nanowire which is under a small bias potential of few meV. We have calculated the transmission of this device first with no scatterers. Then we have introduced scatterers in the device and calculated the transmission for the device.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Valley-dependent Lorentz force and Aharonov-Bohm phase in strained graphene p - n junction

Veselago lens focusing in graphene p-n junction is promising for realizations of new generation electron optics devices. However, the effect of the strain-induced Aharonov-Bohm interference in a p-n junction has not been discussed before. We provide an experimentally feasible setup based on the Veselago lens in which the presence of strain can result in both the valley-dependent Lorentz force and Aharonov-Bohm interference. Inparticular, by employing the Green’s function and tight-binding methods, we study the strain induced by dislocations and line defects in a p-n junction and show how the resulting Aharonov-Bohm phase and interference can be detected. Furthermore, for a different strain configuration, e.g., corresponding to corrugated graphene, we find strong signatures of valley splitting induced by the fictitious magnetic field. Our proposal can be useful for mapping elastic deformations and defects, and for studying valley-dependent effects in graphene.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Toward a predictive theory of correlated materials

Correlated electron materials display a rich variety of notable properties ranging from unconventional superconductivity to metal-insulator transitions. These properties are of interest from the point of view of applications but are hard to treat theoretically, as they result from multiple competing energy scales. Although possible in more weakly correlated materials, theoretical design and spectroscopy of strongly correlated electron materials have been a difficult challenge for many years. In this paper, by treating all the relevant energy scales with sufficient accuracy, complementary advances in Green’s functions and quantum Monte Carlo methods open a path to first-principles computational property predictions in this class of materials.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Electron-phonon interaction in three-barrier nanosystems as active elements of quantum cascade detectors

The theory of electron tunneling through an open nanostructure as an active element of a quantum cascade detector is developed, which takes into account the interaction of electrons with confined and interface phonons. Using the method of finite-temperature Green’s functions and the electron-phonon Hamiltonian in the representation of second quantization over all system variables, the temperature shifts and electron-level widths are calculated and the contributions of different electron-phonon-interaction mechanisms to renormalization of the spectral parameters are analyzed depending on the geometrical configuration of the nanosystem. Due to weak electron-phonon coupling in a GaAs/Al{sub 0.34}Ga{sub 0.66}As-based resonant tunneling nanostructure, the temperature shift and rf field absorption peak width are not very sensitive to the electron-phonon interaction and result from a decrease in potential barrier heights caused by a difference in the temperature dependences of the well and barrier band gaps.

36 MATERIALS SCIENCE↗
NASA NTRS2003

Analysis of Photothermal Characterization of Layered Materials: Design of Optimal Experiments

In this paper numerical calculations are presented for the steady-periodic temperature in layered materials and functionally-graded materials to simulate photothermal methods for the measurement of thermal properties. No laboratory experiments were performed. The temperature is found from a new Green s function formulation which is particularly well-suited to machine calculation. The simulation method is verified by comparison with literature data for a layered material. The method is applied to a class of two-component functionally-graded materials and results for temperature and sensitivity coefficients are presented. An optimality criterion, based on the sensitivity coefficients, is used for choosing what experimental conditions will be needed for photothermal measurements to determine the spatial distribution of thermal properties. This method for optimal experiment design is completely general and may be applied to any photothermal technique and to any functionally-graded material.

Cole, Kevin D.↗
DOE OSTI2014

Expansion techniques for collisionless stellar dynamical simulations

We present graphics processing unit (GPU) implementations of two fast force calculation methods based on series expansions of the Poisson equation. One method is the self-consistent field (SCF) method, which is a Fourier-like expansion of the density field in some basis set; the other method is the multipole expansion (MEX) method, which is a Taylor-like expansion of the Green's function. MEX, which has been advocated in the past, has not gained as much popularity as SCF. Both are particle-field methods and optimized for collisionless galactic dynamics, but while SCF is a 'pure' expansion, MEX is an expansion in just the angular part; thus, MEX is capable of capturing radial structure easily, while SCF needs a large number of radial terms. We show that despite the expansion bias, these methods are more accurate than direct techniques for the same number of particles. The performance of our GPU code, which we call ETICS, is profiled and compared to a CPU implementation. On the tested GPU hardware, a full force calculation for one million particles took ∼0.1 s (depending on expansion cutoff), making simulations with as many as 10{sup 8} particles fast for a comparatively small number of nodes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Rigorous Ab Initio Quantum Embedding for Quantum Chemistry Using Green’s Function Theory: Screened Interaction, Nonlocal Self-Energy Relaxation, Orbital Basis, and Chemical Accuracy

We present a detailed discussion of the self-energy embedding theory (SEET), which is a quantum embedding scheme allowing us to describe a chosen subsystem very accurately while keeping the description of the environment at a lower level. We apply SEET to molecular examples where our chosen subsystem is made out of a set of strongly correlated orbitals while the weakly correlated orbitals constitute an environment. Consequently, a highly accurate method is used to calculate the self-energy for the system, while a lower-level method is employed to find the self-energy for the environment. Such a self-energy separation is very general, and to make the SEET procedure applicable to multiple systems, a detailed and practical procedure for the evaluation of the system and environment self-energy is necessary. We list all of the intricacies for one of the possible procedures while focusing our discussion on many practical implementation aspects such as the choice of best orbital basis, impurity solver, and many steps necessary to reach chemical accuracy. Finally, on a set of carefully chosen molecular examples, we demonstrate that SEET, which is a controlled, systematically improvable Green’s function method, can be as accurate as established wave function quantum chemistry methods.

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

A meshless stochastic method for Poisson–Nernst–Planck equations

A plethora of biological, physical, and chemical phenomena involve transport of charged particles (ions). Its continuum-scale description relies on the Poisson–Nernst–Planck (PNP) system, which encapsulates the conservation of mass and charge. The numerical solution of these coupled partial differential equations is challenging and suffers from both the curse of dimensionality and difficulty in efficiently parallelizing. We present a novel particle-based framework to solve the full PNP system by simulating a drift–diffusion process with time- and space-varying drift. We leverage Green’s functions, kernel-independent fast multipole methods, and kernel density estimation to solve the PNP system in a meshless manner, capable of handling discontinuous initial states. The method is embarrassingly parallel, and the computational cost scales linearly with the number of particles and dimension. We use a series of numerical experiments to demonstrate both the method’s convergence with respect to the number of particles and computational cost vis-à-vis a traditional partial differential equation solver.

Chemistry↗
DOE OSTI2017

Propagator quantum chemical study of S-cis-(Z)-2-(2-formylethenyl)pyrrole: electronic structure and aspects of intramolecular hydrogen bond manifestation in ionization spectra

For the purpose of investigation of the electronic structures of functionalized pyrroles with potential biological activity the electronic structures and ionization spectra of S-cis-(Z)-2-(2-formylethenyl)pyrrole (FP) were calculated by the propagator quantum chemical method. The calculations were performed using the third-order algebraic diagrammatic construction method (ADC(3)) for one-particle Green´s function (electronic propagator) and the 6–31G** basis set. Going from FP (possessing the intramolecular hydrogen bond H⋯O) to its conformation FPR (without H⋯O bonding), the O1s-ionization energy and the ionization energy of the σ-type lone electron pair orbital of the O atom decrease by ~0.2 eV, which is a consequence of stopping the electron density transfer from the O atom. A strong electron density transfer through the hydrogen bond from the O atom to the NH group occurs in the nitrogen core level ionization spectrum as evidenced by a lower N1s-ionization energy of FP (by ~0.7 eV ) compared to that of FPR. The valence shell ionization spectra of FP and FPR calculated using the ADC(3) method are characterized by a high density of the satellite states. The results obtained indicate that the electronic structures of the compounds of the considered class are characterized by pronounced effects of electron correlation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗