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 109 records · Page 6

Green’s functions on a renormalized lattice: An improved method for the integer quantum Hall transition

Highlights: • Speed up of numerical calculations by exact renormalization group steps. • Practical approach to handle delicate scaling analysis. • Absence of evidence to prefer marginal corrections over power-law corrections. We introduce a performance-optimized method to simulate localization problems on bipartite tight-binding lattices. It combines an exact renormalization group step to reduce the sparseness of the original problem with the recursive Green’s function method. We apply this framework to investigate the critical behavior of the integer quantum Hall transition of a tight-binding Hamiltonian defined on a simple square lattice. In addition, we employ an improved scaling analysis that includes two irrelevant exponents to characterize the shift of the critical energy as well as the corrections to the dimensionless Lyapunov exponent. We compare our findings with the results of a conventional implementation of the recursive Green’s function method, and we put them into broader perspective in view of recent development in this field.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Nonequilibrium Green's function formulation of intersubband absorption for nonparabolic single-band effective mass Hamiltonian

The formulas are derived that enable calculations of intersubband absorption coefficient within nonequilibrium Green's function method applied to a single-band effective-mass Hamiltonian with the energy dependent effective mass. The derivation provides also the formulas for the virtual valence band components of the two-band Green's functions which can be used for more exact estimation of the density of states and electrons and more reliable treatment of electronic transport in unipolar n-type heterostructure semiconductor devices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Intersublattice magnetocrystalline anisotropy using a realistic tight-binding method based on maximally localized Wannier functions

Using a realistic tight-binding Hamiltonian based on maximally localized Wannier functions, we investigate the two-ion magnetocrystalline anisotropy energy (MAE) in L 1 0 transition metal compounds. MAE contributions from throughout the Brillouin zone are obtained using magnetic force theorem calculations with and without perturbation theory. The results from both methods agree with each other, and both reflect features of the Fermi surface. The intrasublattice and intersublattice contributions to MAE are evaluated using a Green's function method. We find that the sign of the intersublattice contribution varies among compounds, and that its amplitude may be significant, suggesting MAE can not be resolved accurately in a single-ion manner. The results are further validated by scaling spin-orbit-coupling strength in density functional theory. Altogether, this realistic tight-binding method provides an effective approach to evaluate and analyze MAE while retaining the accuracy of corresponding first-principles methods.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Ab initio simulations of metal contacts for graphene-based devices

The precise atomic structure of a metal contact significantly affects the performance of nanoscale electronic devices. We use an accurate, DFT-based non-equilibrium Green’s function method to evaluate various metal contacts with graphene or graphene nanoribbons. Furthermore, for surface metal contacts not chemically bound to graphene, Ti contacts have lower resistance than those of Au, Ca, Ir, Pt, and Sr. However, as an edge contact, Ti has larger resistance than Au. Bridging O atoms at Ti and Au edge contacts lowers the transmission by over 30%.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Methods for calculating X-ray diffuse scattering from a crystalline medium with spheroidal quantum dots

Two independent approaches to calculate the angular distribution of X-ray diffusion scattering from a crystalline medium with spheroidal quantum dots (QDs) have been proposed. The first method is based on the analytical solution involving the multipole expansion of elastic strain fields beyond QDs. The second approach is based on calculations of atomic displacements near QDs by the Green’s function method. An analysis of the diffuse scattering intensity distribution in the reciprocal space within these two approaches shows that both methods yield similar results for the chosen models of QD spatial distribution.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Subsonic aerodynamic and flutter characteristics of several wings calculated by the SOUSSA P1.1 panel method

The SOUSSA (steady, oscillatory, and unsteady subsonic and supersonic aerodynamics) program is the computational implementation of a general potential flow analysis (by the Green's function method) that can generate pressure distributions on complete aircraft having arbitrary shapes, motions and deformations. Some applications of the initial release version of this program to several wings in steady and oscillatory motion, including flutter are presented. The results are validated by comparisons with other calculations and experiments. Experiences in using the program as well as some recent improvements are described.

Yates, E. C., Jr.↗

Radiative interaction of atmosphere and surface: write up with elements of code

In passive satellite remote sensing of the Earth, separation of the path radiance (atmosphere-only contribution) from the surface reflection remains a “significant challenge”. Recent literature names it among the gaps in radiative transfer (RT) topics that “require continued research in the near future”. The challenge comes from multiple reflections (bouncing) between the atmosphere and surface – radiative interaction. In this paper we use a known RT technique, the matrix-operator method (MOM), and a new modification of the monochromatic vector RT (vRT) code IPOL (Intensity and POLarization) to simulate the interaction of a plane-parallel atmosphere and a few widely used surface reflection models. Following the idea of the Green’s function method, IPOL no longer takes the surface model parameters on input. Instead, it provides the path radiance, and the atmospheric reflection and transmission matrices as output. Despite many RT codes use the MOM formalism, this output does not seem common. The surface reflection matrix is computed externally. Therefore, this paper extends the Green’s function atmospheric correction technique to the case of polarized light. Aiming clarity rather than performance, we explain in Python the structure of the surface matrices for the isotropic (Lambertian), directional unpolarized, and polarized ocean reflection models. We then combine these surface matrices and the precomputed IPOL output to get numerically accurate signal at the top of atmosphere (TOA) and test it vs. published benchmarks. Then, for each benchmark scenario we show how to get the surface from the TOA signal, i.e. perform the RT-based atmospheric correction.

radiative transfer↗

Development of Hexagonal Nodal Diffusion Codes at National Tsing-Hua University

The development of nodal methods for hexagonal geometry reactors at National Tsing-Hua University is reviewed. In particular, the development of the Hybrid Nodal Green's Function Method in hexagonal geometry is re-cast into the Direct-Coarse Mesh Finite Difference Method for efficient implementation with residual solution techniques such as GMRES. The method is further modified to include a rigorous treatment of non-physical terms that arise from the transverse integration procedure. The treatment follows the conformal mapping approach first developed by Y-A Chao and his co-workers. The overall development is implemented into a code written in C++ and tested on a number of well known benchmarks. The new code is shown to be very highly accurate while extremely efficient computationally.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Adjoint method and runaway electron avalanche

The adjoint method for the study of runaway electron dynamics in momentum space Liu et al (2016 Phys. Plasmas 23 010702) is rederived using the Green's function method, for both the runaway probability function (RPF) and the expected loss time (ELT). The RPF and ELT obtained using the adjoint method are presented, both with and without the synchrotron radiation reaction force. In conclusion, the adjoint method is then applied to study the runaway electron avalanche. Both the critical electric field and the growth rate for the avalanche are calculated using this fast and novel approach.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Vibrational properties of TaW alloy using modified embedded atom method potential

Force-constants up to second neighbours of pure transition metal Ta and TaW alloy are determined using the modified embedded atom method (MEAM) potential. The obtained force-constants are used to calculate the phonon dispersion of pure Ta and TaW alloy. As a further application of MEAM potential, the force-constants are used to calculate the local vibrational density of states and mean square thermal displacements of pure Ta and W impurity atoms with Green’s function method. The calculated results are found to be in agreement with the experimental measurements.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Absence of Coulomb Blockade in the Anderson Impurity Model at the Symmetric Point

In this study, we investigate the characteristics of the electric current in the so-called symmetric Anderson impurity model. We study the nonequilibrium model using two complementary approximate methods, the perturbative quantum master equation approach to the reduced density matrix and a self-consistent equation of motion approach to the nonequilibrium Green's function. We find that, at a particular symmetry point, an interacting Anderson impurity model recovers the same steady-state current as an equivalent noninteracting model, akin a two-band resonant level model. We show this in the Coulomb blockade regime for both high and low temperatures, where either the approximate master equation approach or the Green's function method provides accurate results for the current. We conclude that the steady-state current in the symmetric Anderson model at this regime does not encode characteristics of a many-body interacting system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Metamaterial-based theoretical description of light scattering by metallic nano-hole array structures

We have experimentally and theoretically investigated the light-matter interaction in metallic nano-hole array structures. The scattering cross section spectrum was measured for three samples each having a unique nano-hole array radius and periodicity. Each measured spectrum had several peaks due to surface plasmon polaritons. The dispersion relation and the effective dielectric constant of the structure were calculated using transmission line theory and Bloch's theorem. Using the effective dielectric constant and the transfer matrix method, the surface plasmon polariton energies were calculated and found to be quantized. Using these quantized energies, a Hamiltonian for the surface plasmon polaritons was written in the second quantized form. Working with the Hamiltonian, a theory of scattering cross section was developed based on the quantum scattering theory and Green's function method. For both theory and experiment, the location of the surface plasmon polariton spectral peaks was dependant on the array periodicity and radii of the nano-holes. Good agreement was observed between the experimental and theoretical results. It is proposed that the newly developed theory can be used to facilitate optimization of nanosensors for medical and engineering applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

The Determination of the Large-Scale Circulation of the Pacific Ocean from Satellite Altimetry using Model Green's Functions

A Green's function method for obtaining an estimate of the ocean circulation using both a general circulation model and altimetric data is demonstrated. The fundamental assumption is that the model is so accurate that the differences between the observations and the model-estimated fields obey a linear dynamics. In the present case, the calculations are demonstrated for model/data differences occurring on very a large scale, where the linearization hypothesis appears to be a good one. A semi-automatic linearization of the Bryan/Cox general circulation model is effected by calculating the model response to a series of isolated (in both space and time) geostrophically balanced vortices. These resulting impulse responses or 'Green's functions' then provide the kernels for a linear inverse problem. The method is first demonstrated with a set of 'twin experiments' and then with real data spanning the entire model domain and a year of TOPEX/POSEIDON observations. Our present focus is on the estimate of the time-mean and annual cycle of the model. Residuals of the inversion/assimilation are largest in the western tropical Pacific, and are believed to reflect primarily geoid error. Vertical resolution diminishes with depth with 1 year of data. The model mean is modified such that the subtropical gyre is weakened by about 1 cm/s and the center of the gyre shifted southward by about 10 deg. Corrections to the flow field at the annual cycle suggest that the dynamical response is weak except in the tropics, where the estimated seasonal cycle of the low-latitude current system is of the order of 2 cm/s. The underestimation of observed fluctuations can be related to the inversion on the coarse spatial grid, which does not permit full resolution of the tropical physics. The methodology is easily extended to higher resolution, to use of spatially correlated errors, and to other data types.

Stammer, Detlef↗

Chemical Source Inversion using Assimilated Constituent Observations in an Idealized Two-dimensional System

We present a source inversion technique for chemical constituents that uses assimilated constituent observations rather than directly using the observations. The method is tested with a simple model problem, which is a two-dimensional Fourier-Galerkin transport model combined with a Kalman filter for data assimilation. Inversion is carried out using a Green's function method and observations are simulated from a true state with added Gaussian noise. The forecast state uses the same spectral spectral model, but differs by an unbiased Gaussian model error, and emissions models with constant errors. The numerical experiments employ both simulated in situ and satellite observation networks. Source inversion was carried out by either direct use of synthetically generated observations with added noise, or by first assimilating the observations and using the analyses to extract observations. We have conducted 20 identical twin experiments for each set of source and observation configurations, and find that in the limiting cases of a very few localized observations, or an extremely large observation network there is little advantage to carrying out assimilation first. However, in intermediate observation densities, there decreases in source inversion error standard deviation using the Kalman filter algorithm followed by Green's function inversion by 50% to 95%.

Tangborn, Andrew↗

Curie temperature of ultrathin ferromagnetic layer with Dzyaloshinskii-Moriya interaction

We investigate the effect of the Dzyaloshinskii-Moriya interaction (DMI) on the Curie temperature of the ultrathin ferromagnetic layers. It has been known that the Curie temperature of the ferromagnet depends on spin wave excitation energies, and they are affected by DMI. Therefore, the ferromagnetic transition temperature of the ultrathin ferromagnetic layer must be sensitive on the DMI. We find that the Curie temperature depends on the DMI by using the double time Green's function method. Since the DMI is arisen by the inversion symmetry breaking structure, the DMI is always important in the inversion symmetry breaking ultrathin ferromagnetic layers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Single particle density of trapped interacting quantum gases

An expression for single particle density for trapped interacting gases has been obtained in first order of interaction using Green’s function method. Results are easily simplified for homogeneous quantum gases and are found to agree with famous results obtained by Huang-Yang-Luttinger and Lee-Yang.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Phonon dispersion and local density of states in NiPd alloy using modified embedded atom method potential

A modified embedded atom method (MEAM) potential model up to second neighbours has been used to calculate the phonon dispersions for Ni{sub 0.55}Pd{sub 0.45} alloy in which Pd is introduced as substitutional impurity. Using the force-constants obtained from MEAM potential, the local vibrational density of states in host Ni and substitutional Pd atoms using Green’s function method has been calculated. The calculation of phonon dispersions of NiPd alloy shows a good agreement with the experimental results. Condition of resonance mode has also been investigated and resonance mode in the frequency spectrum of impurity atom at low frequency is observed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗