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 1,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Numerical evaluation of Auger recombination coefficients in relaxed and strained germanium

The potential applications of germanium and its alloys in infrared silicon-based photonics have led to a renewed interest in their optical properties. In this letter, we report on the numerical determination of Auger coefficients at T = 300 K for relaxed and biaxially strained germanium. We use a Green's function based model that takes into account all relevant direct and phonon-assisted processes and perform calculations up to a strain level corresponding to the transition from indirect to direct energy gap. We have considered excess carrier concentrations ranging from 10{sup 16} cm{sup −3} to 5 × 10{sup 19} cm{sup −3}. For use in device level simulations, we also provide fitting formulas for the calculated electron and hole Auger coefficients as functions of carrier density.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Composite electron propagator methods for calculating ionization energies

Accurate ionization energies of molecules may be determined efficiently with composite electron-propagator (CEP) techniques. These methods estimate the results of a calculation with an advanced correlation method and a large basis set by performing a series of more tractable calculations in which large basis sets are used with simpler approximations and small basis sets are paired with more demanding correlation techniques. The performance of several CEP methods, in which diagonal, second-order electron propagator results with large basis sets are combined with higher-order results obtained with smaller basis sets, has been tested for the ionization energies of closed-shell molecules from the G2 set. Useful compromises of accuracy and computational efficiency employ complete-basis-set extrapolation for second-order results and small basis sets in third-order, partial third-order, renormalized partial-third order, or outer valence Green’s function calculations. Analysis of results for vertical as well as adiabatic ionization energies leads to specific recommendations on the best use of regular and composite methods. Results for 22 organic molecules of interest in the design of photovoltaic devices, benzo[a]pyrene, Mg-octaethylporphyrin, and C{sub 60} illustrate the capabilities of CEP methods for calculations on large molecules.

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

The electronic heat capacity of YBa{sub 2}Cu{sub 3}O{sub 7−δ} superconductor

The contributions due to the point defects or disorder and anharmonicities which play deterministic role in the understanding of electronic heat capacity (EHC) of high temperature superconductors (HTS) have been investigated via electron density of states (EDOS) approach on the basis of quantum dynamical many body theory. The evaluation of EDOS has been carried out with the help of most versatile method of double time temperature dependent electron Green’s functions (GF) via a Hamiltonian (non BCS type) which includes the effects of electrons, phonons, defects, anharmonicity, and electron-phonon interactions which enables to account the effects of cubic anharmonicity besides with both the force constant changes and mass difference caused by the impurities in developing the results for EDOS and EHC. The new results reveal some striking features of EHC of HTS.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Charge transport properties of graphene: Effects of Cu-based gate electrode

Using the first-principles nonequilibrium Green's function method, we study effects of Cu and Ni@Cu used as the Cu-based gate electrode on the charge transport of graphene in the field effect transistors (FET). We find that the transmission of graphene decreases with both Cu and Ni@Cu absorbed in the scatter region. Especially, noticeable transmission gaps are present around the Femi level. The transmission gaps are still effective, and considerable cut-off regions are found under the non-equilibrium environment. The Ni@Cu depresses the transmission of graphene more seriously than the Cu and enlarges the transmission gap in armchair direction. The effects on the charge transport are attributed to the redistribution of electronic states of graphene. Both Cu and Ni@Cu induce the localization of states, so as to block the electronic transport. The Ni@Cu transforms the interaction between graphene and gate electrode from the physisorption to the chemisorption, and then induces more localized states, so that the transmission decreases further. Our results suggest that besides being used to impose gate voltage, the Cu-based gate electrode itself will have a considerable effect on the charge transport of graphene and induces noticeable transmission gap in the FET.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Effects of the molecule-electrode interface on the low-bias conductance of Cu–H{sub 2}–Cu single-molecule junctions

The atomic structure and electronic transport properties of a single hydrogen molecule connected to both symmetric and asymmetric Cu electrodes are investigated by using the non-equilibrium Green’s function formalism combined with the density functional theory. Our calculations show that in symmetric Cu–H{sub 2}–Cu junctions, the low-bias conductance drops rapidly upon stretching, while asymmetric ones present a low-bias conductance spanning the 0.2–0.3 G{sub 0} interval for a wide range of electrode separations. This is in good agreement with experiments on Cu atomic contacts in a hydrogen environment. Furthermore, the distribution of the calculated vibrational energies of the two hydrogen atoms in the asymmetric Cu–H{sub 2}–Cu junction is also consistent with experiments. These findings provide clear evidence for the formation of asymmetric Cu–H{sub 2}–Cu molecular junctions in breaking Cu atomic contacts in the presence of hydrogen and are also helpful for the design of molecular devices with Cu electrodes.

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

A many-body states picture of electronic friction: The case of multiple orbitals and multiple electronic states

We present a very general form of electronic friction as present when a molecule with multiple orbitals hybridizes with a metal electrode. To develop this picture of friction, we embed the quantum-classical Liouville equation (QCLE) within a classical master equation (CME). Thus, this article extends our previous work analyzing the case of one electronic level, as we may now treat the case of multiple levels and many electronic molecular states. We show that, in the adiabatic limit, where electron transitions are much faster than nuclear motion, the QCLE-CME reduces to a Fokker-Planck equation, such that nuclei feel an average force as well as friction and a random force—as caused by their interaction with the metallic electrons. Finally, we show numerically and analytically that our frictional results agree with other published results calculated using non-equilibrium Green’s functions. Numerical recipes for solving this QCLE-CME will be provided in a subsequent paper.

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

Controlling dynamical thermal transport of biased bilayer graphene by impurity atoms

We address the dynamical thermal conductivity of biased bilayer graphene doped with acceptor impurity atoms for AA-stacking in the context of tight binding model Hamiltonian. The effect of scattering by dilute charged impurities is discussed in terms of the self-consistent Born approximation. Green’s function approach has been exploited to find the behavior of thermal conductivity of bilayer graphene within the linear response theory. We have found the frequency dependence of thermal conductivity for different values of concentration and scattering strength of dopant impurity. Also the dependence of thermal conductivity on the impurity concentration and bias voltage has been investigated in details.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2016

An O ( N ) and parallel approach to integral problems by a kernel-independent fast multipole method: Application to polarization and magnetization of interacting particles

Large classes of materials systems in physics and engineering are governed by magnetic and electrostatic interactions. Continuum or mesoscale descriptions of such systems can be cast in terms of integral equations, whose direct computational evaluation requires O ( N 2 ) operations, where N is the number of unknowns. Such a scaling, which arises from the many-body nature of the relevant Green's function, has precluded wide-spread adoption of integral methods for solution of large-scale scientific and engineering problems. In this work, a parallel computational approach is presented that relies on using scalable open source libraries and utilizes a kernel-independent Fast Multipole Method (FMM) to evaluate the integrals in O ( N ) operations, with O ( N ) memory cost, thereby substantially improving the scalability and efficiency of computational integral methods. We demonstrate the accuracy, efficiency, and scalability of our approach in the context of two examples. In the first, we solve a boundary value problem for a ferroelectric/ferromagnetic volume in free space. In the second, we solve an electrostatic problem involving polarizable dielectric bodies in an unbounded dielectric medium. Lastly, the results from these test cases show that our proposed parallel approach, which is built on a kernel-independent FMM, can enable highly efficient and accurate simulations and allow for considerable flexibility in a broad range of applications.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

A sub k{sub B}T/q semimetal nanowire field effect transistor

The key challenge for nanoelectronics technologies is to identify the designs that work on molecular length scales, provide reduced power consumption relative to classical field effect transistors (FETs), and that can be readily integrated at low cost. To this end, a FET is introduced that relies on the quantum effects arising for semimetals patterned with critical dimensions below 5 nm, that intrinsically has lower power requirements due to its better than a “Boltzmann tyranny” limited subthreshold swing (SS) relative to classical field effect devices, eliminates the need to form heterojunctions, and mitigates against the requirement for abrupt doping profiles in the formation of nanowire tunnel FETs. This is achieved through using a nanowire comprised of a single semimetal material while providing the equivalent of a heterojunction structure based on shape engineering to avail of the quantum confinement induced semimetal-to-semiconductor transition. Ab initio calculations combined with a non-equilibrium Green's function formalism for charge transport reveals tunneling behavior in the OFF state and a resonant conduction mechanism for the ON state. A common limitation to tunnel FET (TFET) designs is related to a low current in the ON state. A discussion relating to the semimetal FET design to overcome this limitation while providing less than 60 meV/dec SS at room temperature is provided.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

State hybridization shapes the photocurrent in triple quantum dot nanojunctions

We investigated a prototype of a quantum dot based photodetector made of a dot absorber interconnected with two lateral dot filters in contact with semiconducting leads. Using the nonequilibrium Green's function technique, we found that there are two opposite evolutions of the photocurrent in the vicinity of the tunnel resonance for such a kind of nanodevice. This evolution depends on where the strongest hybridization of states happens, and hence still reveals a quantum effect. If the filter states hybridize more with the absorber states than the ones of the leads, the photocurrent shows a maximum at the tunnel resonance, while it is minimized in the opposite case.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Sub-millitesla magnetic field effects on the recombination reaction of flavin and ascorbic acid radicals

Even though the interaction of a <1 mT magnetic field with an electron spin is less than a millionth of the thermal energy at room temperature (k{sub B}T), it still can have a profound effect on the quantum yields of radical pair reactions. We present a study of the effects of sub-millitesla magnetic fields on the photoreaction of flavin mononucleotide with ascorbic acid. Direct control of the reaction pathway is achieved by varying the rate of electron transfer from ascorbic acid to the photo-excited flavin. At pH 7.0, we verify the theoretical prediction that, apart from a sign change, the form of the magnetic field effect is independent of the initial spin configuration of the radical pair. The data agree well with model calculations based on a Green’s function approach that allows multinuclear spin systems to be treated including the diffusive motion of the radicals, their spin-selective recombination reactions, and the effects of the inter-radical exchange interaction. The protonation states of the radicals are uniquely determined from the form of the magnetic field-dependence. At pH 3.0, the effects of two chemically distinct radical pair complexes combine to produce a pronounced response to ∼500 μT magnetic fields. These findings are relevant to the magnetic responses of cryptochromes (flavin-containing proteins proposed as magnetoreceptors in birds) and may aid the evaluation of effects of weak magnetic fields on other biologically relevant electron transfer processes.

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

Mesoscopic-microscopic spatial stochastic simulation with automatic system partitioning

The reaction-diffusion master equation (RDME) is a model that allows for efficient on-lattice simulation of spatially resolved stochastic chemical kinetics. Compared to off-lattice hard-sphere simulations with Brownian dynamics or Green’s function reaction dynamics, the RDME can be orders of magnitude faster if the lattice spacing can be chosen coarse enough. However, strongly diffusion-controlled reactions mandate a very fine mesh resolution for acceptable accuracy. It is common that reactions in the same model differ in their degree of diffusion control and therefore require different degrees of mesh resolution. This renders mesoscopic simulation inefficient for systems with multiscale properties. Mesoscopic-microscopic hybrid methods address this problem by resolving the most challenging reactions with a microscale, off-lattice simulation. However, all methods to date require manual partitioning of a system, effectively limiting their usefulness as “black-box” simulation codes. In this paper, we propose a hybrid simulation algorithm with automatic system partitioning based on indirect a priori error estimates. We demonstrate the accuracy and efficiency of the method on models of diffusion-controlled networks in 3D.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Similarity-transformed equation-of-motion vibrational coupled-cluster theory

A similarity-transformed equation-of-motion vibrational coupled-cluster (STEOM-XVCC) method is introduced as a one-mode theory with an effective vibrational Hamiltonian, which is similarity transformed twice so that its lower-order operators are dressed with higher-order anharmonic effects. The first transformation uses an exponential excitation operator, defining the equation-of-motion vibrational coupled-cluster (EOM-XVCC) method, and the second uses an exponential excitation-deexcitation operator. From diagonalization of this doubly similarity-transformed Hamiltonian in the small one-mode excitation space, the method simultaneously computes accurate anharmonic vibrational frequencies of all fundamentals, which have unique significance in vibrational analyses. We establish a diagrammatic method of deriving the working equations of STEOM-XVCC and prove their connectedness and thus size-consistency as well as the exact equality of its frequencies with the corresponding roots of EOM-XVCC. We furthermore elucidate the similarities and differences between electronic and vibrational STEOM methods and between STEOM-XVCC and vibrational many-body Green’s function theory based on the Dyson equation, which is also an anharmonic one-mode theory. The latter comparison inspires three approximate STEOM-XVCC methods utilizing the common approximations made in the Dyson equation: the diagonal approximation, a perturbative expansion of the Dyson self-energy, and the frequency-independent approximation. The STEOM-XVCC method including up to the simultaneous four-mode excitation operator in a quartic force field and its three approximate variants are formulated and implemented in computer codes with the aid of computer algebra, and they are applied to small test cases with varied degrees of anharmonicit

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

Algebraic motion of vertically displacing plasmas

In this paper, the vertical motion of a tokamak plasma is analytically modelled during its non-linear phase by a free-moving current-carrying rod inductively coupled to a set of fixed conducting wires or a cylindrical conducting shell. The solutions capture the leading term in a Taylor expansion of the Green's function for the interaction between the plasma column and the surrounding vacuum vessel. The plasma shape and profiles are assumed not to vary during the vertical drifting phase such that the plasma column behaves as a rigid body. In the limit of perfectly conducting structures, the plasma is prevented to come in contact with the wall due to steep effective potential barriers created by the induced Eddy currents. Resistivity in the wall allows the equilibrium point to drift towards the vessel on the slow timescale of flux penetration. The initial exponential motion of the plasma, understood as a resistive vertical instability, is succeeded by a non-linear “sinking” behaviour shown to be algebraic and decelerating. Finally, the acceleration of the plasma column often observed in experiments is thus concluded to originate from an early sharing of toroidal current between the core, the halo plasma, and the wall or from the thermal quench dynamics precipitating loss of plasma current.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

On simulation of local fluxes in molecular junctions

We present a pedagogical review of the current density simulation in molecular junction models indicating its advantages and deficiencies in analysis of local junction transport characteristics. In particular, we argue that current density is a universal tool which provides more information than traditionally simulated bond currents, especially when discussing inelastic processes. However, current density simulations are sensitive to the choice of basis and electronic structure method. We note that while discussing the local current conservation in junctions, one has to account for the source term caused by the open character of the system and intra-molecular interactions. Our considerations are illustrated with numerical simulations of a benzenedithiol molecular junction.

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

Sesame-style decomposition of KS-DFT molecular dynamics for direct interrogation of nuclear models

A common paradigm used in the construction of equations of state is to decompose the thermodynamics into a superposition of three terms: a static-lattice cold curve, a contribution from the thermal motion of the nuclei, and a contribution from the thermal excitation of the electrons. While statistical mechanical models for crystals provide tractable framework for the nuclear contribution in the solid phase, much less is understood about the nuclear contribution above the % melt temperature $(C^{(nuc)}_ν≈3R)$ and how it should transition to the high-temperature limit $(C^{(nuc)}_ν\sim\frac{3}{2}R)$. Here, we describe an algorithm for extracting both the thermal nuclear and thermal electronic contributions from quantum molecular dynamics (QMD) simulations. We then use the VASP QMD package to probe thermal nuclear behavior of liquid aluminum at normal density to compare the results to semi-empirical models – the Johnson generic model, the Chisolm high-temperature liquid model, and the CRIS model. The thermal electronic energy is also calculated from the VASP results and compared to the Kohn-Sham calculations from QE-DFT simulations, Green’s function average atom calculations, and the Thomas-Fermi atomic model.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2019

A simulation-based model for understanding the time dependent x-ray drive asymmetries and error bars in indirectly driven implosions on the National Ignition Facility

Time-dependent low-mode asymmetries are believed to play a leading role in limiting the performance of current inertial confinement fusion implosions on the National Ignition Facility (NIF) [E. I. Moses et al., Phys. Plasmas 16, 041006 (2009)]. These long wavelength modes are initiated and driven by asymmetries in the x-ray flux from the hohlraum; however, the underlying hydrodynamics of the implosion also act to modify and amplify these asymmetries. Here we present here a simulation-based model connecting the time-dependent drive asymmetry seen by the capsule to the measured inflight and hot spot symmetries. This approach is based on a Green's function analysis for which we evaluate the response of the capsule to impulses of drive asymmetry at a series of times. Our model sheds new light on the sensitivity to the drive asymmetry of an imploded capsule, giving a new tool for design. Inverting the problem and finding the drive asymmetry needed to match the experimental data allow us to tightly constrain the drive asymmetry seen by the capsule, providing an error estimate on the result. Doing so, we are able to point out when and how the complex hohlraum simulations start to deviate from what they should obtain to match the experimental data. Ultimately, we project to use this model to make some experimental recommendations to fix the time-dependent low-mode asymmetry of indirectly driven implosions and identify additional measurements to further constrain the asymmetries with a view to improving target design on the NIF.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

An experimental and theoretical study of the C 1s ionization satellites in CH 3 I

We report the C 1s ionization spectrum of CH 3 I has been studied both experimentally and theoretically. Synchrotron radiation has been employed to record polarization dependent photoelectron spectra at a photon energy of 614 eV. These spectra encompass the main-line due to the C 1s single-hole state and the peaks associated with the shake-up satellites. Vertical ionization energies and relative photoelectron intensities have been computed using the fourth-order algebraic-diagrammatic construction approximation scheme for the one-particle Green’s function and the 6-311++G** basis set. The theoretical spectrum derived from these calculations agrees qualitatively with the experimental results, thereby allowing the principal spectral features to be assigned. According to our calculations, two 2 A 1 shake-up states of the C 1s -1 σ CI → σ$*\atop{CI}$ type with singlet and triplet intermediate coupling of the electron spins (S' = 0, 1) play an important role in the spectrum and contribute significantly to the overall intensity. Both of these states are expected to have dissociative diabatic potential energy surfaces with respect to the C–I separation. Whereas the upper of these states perturbs the manifold of Rydberg states, the lower state forms a band which is characterized by a strongly increased width. Our results indicate that the lowest shake-up peak with significant spectral intensity is due to the pair (S' = 0, 1) of 2 E (C 1s -1 I 5p → σ$*\atop{CI}$) states. We predict that these 2 E states acquire photoelectron intensity due to spin-orbit interaction. Lastly, such interactions play an important role here due to the involvement of the I 5p orbitals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗