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 289 records · Page 16

Quantum Monte Carlo Method in the Steady State

We present a numerically exact steady-state inchworm Monte Carlo method for nonequilibrium quantum impurity models. Rather than propagating an initial state to long times, the method is directly formulated in the steady state. This eliminates any need to traverse the transient dynamics and grants access to a much larger range of parameter regimes at vastly reduced computational costs. We benchmark the method on equilibrium Green’s functions of quantum dots in the noninteracting limit and in the unitary limit of the Kondo regime. We then consider correlated materials described with dynamical mean field theory and driven away from equilibrium by a bias voltage. We show that the response of a correlated material to a bias voltage differs qualitatively from the splitting of the Kondo resonance observed in bias-driven quantum dots.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Cumulant methods for electron-phonon problems. I. Perturbative expansions

Here in this work, we investigate the ability of the cumulant expansion (CE) to capture one-particle spectral information in electron-phonon coupled systems at both zero and finite temperatures. In particular, we present a comprehensive study of the second- and fourth-order CEs for the one-dimensional Holstein model as compared with numerically exact methods. We investigate both finite sized systems as well as the approach to the thermodynamic limit, drawing distinctions, and connections between the behavior of systems in and away from the thermodynamic limit that enable a greater understanding of the ability of the CE to capture real-frequency information across the full range of wave vectors. We find that for zero electronic momentum, the spectral function is well described by the second-order CE at low and high temperatures. However, for nonzero electronic momenta, the CE is only accurate at high temperature. We analyze the fourth-order cumulant and find that while it improves the description of the short-time dynamics encoded in the one-particle Green's function, it can introduce divergences in the time domain as well as unphysical negative spectral weight in the spectral function. When well-behaved, the fourth-order CE does provide notable accurate corrections to the second-order CE. Finally, we use our results to comment on the use of the CE as a tool for calculating transport behavior in the realistic ab initio modeling of materials.

36 MATERIALS SCIENCE↗

Ballistic performance comparison of monolayer transition metal dichalcogenide MX{sub 2} (M = Mo, W; X = S, Se, Te) metal-oxide-semiconductor field effect transistors

We study the transport properties of monolayer MX{sub 2} (M = Mo, W; X = S, Se, Te) n- and p-channel metal-oxide-semiconductor field effect transistors (MOSFETs) using full-band ballistic non-equilibrium Green's function simulations with an atomistic tight-binding Hamiltonian with hopping potentials obtained from density functional theory. We discuss the subthreshold slope, drain-induced barrier lowering (DIBL), as well as gate-induced drain leakage (GIDL) for different monolayer MX{sub 2} MOSFETs. We also report the possibility of negative differential resistance behavior in the output characteristics of nanoscale monolayer MX{sub 2} MOSFETs.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Phonon localization in binary alloys with diagonal and off-diagonal disorder: A cluster Green's function approach

We report the development and application of a method for carrying out computational investigations of the effects of mass and force-constant (FC) disorder on phonon spectra. The method is based on the recently developed typical medium dynamical cluster approach, which is a Green's function approach. Excellent quantitative agreement with previous exact diagonalization results establishes the veracity of the method. Application of the method to a model system of binary mass and an FC-disordered system leads to several findings. A narrow resonance, significantly below the van Hove singularity, that has been termed the boson peak, is seen to emerge for low soft particle concentrations. We show, using the typical phonon spectrum, that the states constituting the boson peak cross over from being completely localized to being extended as a function of increasing soft particle concentration. In general, an interplay of mass and FC disorder is found to be cooperative in nature, enhancing phonon localization over all frequencies. However, for a certain range of frequencies, and depending on material parameters, FC disorder can delocalize the states that were localized by mass disorder, and vice versa. Modeling vacancies as weakly bonded sites with vanishing mass, we find that vacancies, even at very low concentrations, are extremely effective in localizing phonons. Thus, inducing vacancies is proposed as a promising route for efficient thermoelectrics. Lastly, we use model parameters corresponding to the alloy system, Ni 1- x Pt x , and we show that mass disorder alone is insufficient to explain the pseudogap in the phonon spectrum; the concomitant presence of FC disorder is necessary.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Theoretical insights on the electro-thermal transport properties of monolayer MoS{sub 2} with line defects

Two dimensional (2D) materials demonstrate several novel electrical, mechanical, and thermal properties which are quite distinctive to those of their bulk form. Among many others, one important potential application of the 2D material is its use in the field of energy harvesting. Owing to that, here we present a detailed study on electrical as well as thermal transport of monolayer MoS{sub 2}, in quasi ballistic regime. Besides the perfect monolayer in its pristine form, we also consider various line defects which have been experimentally observed in mechanically exfoliated MoS{sub 2} samples. For calculating various parameters related to the electrical transmission, we employ the non-equilibrium Green's function-density functional theory combination. However, to obtain the phonon transmission, we take help of the parametrized Stillinger-Weber potential which can accurately delineate the inter-atomic interactions for the monolayer MoS{sub 2}. Due to the presence of line defects, we observed significant reductions in both the charge carrier and the phonon transmissions through a monolayer MoS{sub 2} flake. Moreover, we also report a comparative analysis showing the temperature dependency of the thermoelectric figure of merit values, as obtained for the perfect as well as the other defective 2D samples.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Modeling Sound Propagation Through Non-Axisymmetric Jets

A method for computing the far-field adjoint Green's function of the generalized acoustic analogy equations under a locally parallel mean flow approximation is presented. The method is based on expanding the mean-flow-dependent coefficients in the governing equation and the scalar Green's function in truncated Fourier series in the azimuthal direction and a finite difference approximation in the radial direction in circular cylindrical coordinates. The combined spectral/finite difference method yields a highly banded system of algebraic equations that can be efficiently solved using a standard sparse system solver. The method is applied to test cases, with mean flow specified by analytical functions, corresponding to two noise reduction concepts of current interest: the offset jet and the fluid shield. Sample results for the Green's function are given for these two test cases and recommendations made as to the use of the method as part of a RANS-based jet noise prediction code.

Leib, Stewart J.↗

Understanding the Effect of Twisting Graphene Sheet on Its Magnetoresistance and Spin Filtration Properties

Spin-dependent quantum transport properties in twisted graphene is calculated using density functional theory (DFT) and non-equilibrium Green’s function (NEGF) formulation. Twisting increases the band gap and reduces the spin transport in twisted graphene. At low bias, twisting restricts any spin-up current in antiparallel configuration (APC) which results in higher magnetoresistance (MR), and at higher bias, twisting increases the spin-up current slightly, resulting in lower MR. The MR obtained in twisted graphene at low biases is greater than that of pristine graphene. High spin filtration is observed in parallel configuration and APC for twisted graphene at higher biases which is greater than that of pristine graphene, but at low voltages, pristine gives higher spin filtration than twisted graphene.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Enhanced Magnetoresistance in In-Plane Monolayer MoS{sub 2} with CrO{sub 2} Electrodes

Magnetoresistance of monolayer MoS{sub 2} is reported to increase when used in in-plane configuration with CrO{sub 2} half-metal ferromagnet (HMF) electrodes. Density functional theory (DFT)- and non-equilibrium Green’s function (NEGF)-based simulations show that high magnetoresistance (MR) values of ∼860% can be achieved in in-plane monolayer MoS{sub 2} with CrO{sub 2} electrodes, which is much higher than the MR value of ∼300% for nine-layer and the MR value of ∼70% for single-layer out-of-plane MoS{sub 2} reported in Dolui et al., Phys. Rev. B 90(R), 041401 (2014) past by other researchers. High spin-injection efficiency ∼100% is also obtained at high bias voltages. High MR and perfect spin filtration suggests the importance of this configuration in spintronics applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

The electronic transport behavior of hybridized zigzag graphene and boron nitride nanoribbons

In this present work, we have investigated the electronic transport properties of the hybridized structure constructed by the zigzag graphene and boron-nitride (BN) nanoribbons (Z-B{sub n}N{sub m}C{sub p}, n + m + p = 16) through employing nonequilibrium Green's functions in combination with the density-functional theory. The results demonstrate that the electronic transport properties of the hybridized Z-B{sub n}N{sub m}C{sub p} nanoribbons are strongly dependent on the width of boron-nitride or graphene nanoribbons. When the numbers of n and m are not equal, the negative differential resistance behavior is observed, which can be modulated by varying the width of BN nanoribbons. The conductance of the hybridized Z-B{sub n}N{sub m}C{sub p} nanoribbons with odd numbers of zigzag carbon chains also increases by the width of BN nanoribbons.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Impact of incomplete metal coverage on the electrical properties of metal-CNT contacts: A large-scale ab initio study

Using a dedicated combination of the non-equilibrium Green function formalism and large-scale density functional theory calculations, we investigated how incomplete metal coverage influences two of the most important electrical properties of carbon nanotube (CNT)-based transistors: contact resistance and its scaling with contact length, and maximum current. These quantities have been derived from parameter-free simulations of atomic systems that are as close as possible to experimental geometries. Physical mechanisms that govern these dependences have been identified for various metals, representing different CNT-metal interaction strengths from chemisorption to physisorption. Our results pave the way for an application-oriented design of CNT-metal contacts.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Spectral function of the Holstein polaron at finite temperature

Here we compute the Holstein polaron spectral function on a one-dimensional ring using the finite-temperature T Lanczos method. With increasing T additional features in the spectral function emerge even at temperatures below the phonon frequency. We observe a substantial spread of the spectral weight towards lower frequencies and the broadening of the quasiparticle (QP) peak. In the weak-coupling regime the QP peak merges with the continuum in the high-T limit. In the strong-coupling regime the main features of the low-T spectral function remain detectable up to the highest T used in our calculations. The effective polaron mass shows a nonmonotonic behavior as a function of T at small phonon frequency but increases with T at larger frequencies. The self-energy remains k independent even at elevated T in the frequency range corresponding to the polaron band, while at higher frequencies it develops a distinguishable k dependence. Analytical expressions for the first few frequency moments are derived, and they agree well with those extracted from numerical calculations in a wide-T regime.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Out-of-equilibrium Majorana zero modes in interacting Kitaev chains

Here, we employ a time-dependent real-space local density-of-states method to study the movement and fusion of Majorana zero modes in the one-dimensional interacting Kitaev model, based on the time evolution of many-body states. We analyze the dynamics and both fusion channels of Majoranas using time-dependent potentials, either creating walls or wells. For fast moving Majoranas, we unveil nonequilibrium signatures of the “strong-zero-mode” operator (quasiparity degeneracy in the full spectrum) and its breakdown in the presence of repulsive Coulomb interactions. Focusing on forming a full electron after fusion, we also discuss the upper and lower limits on the Majorana speed needed to reduce nonadiabatic effects and to avoid poisoning due to decoherence.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Modeling ion sensing in molecular electronics

We examine the ability of molecules to sense ions by measuring the change in molecular conductance in the presence of such charged species. The detection of protons (H{sup +}), alkali metal cations (M{sup +}), calcium ions (Ca{sup 2+}), and hydronium ions (H{sub 3}O{sup +}) is considered. Density functional theory (DFT) is used within the Keldysh non-equilibrium Green's function framework (NEGF) to model electron transport properties of quinolinedithiol (QDT, C{sub 9}H{sub 7}NS{sub 2}), bridging Al electrodes. The geometry of the transport region is relaxed with DFT. The transport properties of the device are modeled with NEGF-DFT to determine if this device can distinguish among the M{sup +} + QDT species containing monovalent cations, where M{sup +} = H{sup +}, Li{sup +}, Na{sup +}, or K{sup +}. Because of the asymmetry of QDT in between the two electrodes, both positive and negative biases are considered. The electron transmission function and conductance properties are simulated for electrode biases in the range from −0.5 V to 0.5 V at increments of 0.1 V. Scattering state analysis is used to determine the molecular orbitals that are the main contributors to the peaks in the transmission function near the Fermi level of the electrodes, and current-voltage relationships are obtained. The results show that QDT can be used as a proton detector by measuring transport through it and can conceivably act as a pH sensor in solutions. In addition, QDT may be able to distinguish among different monovalent species. This work suggests an approach to design modern molecular electronic conductance sensors with high sensitivity and specificity using well-established quantum chemistry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Fractional charge and spin errors in self-consistent Green’s function theory

We examine fractional charge and spin errors in self-consistent Green’s function theory within a second-order approximation (GF2). For GF2, it is known that the summation of diagrams resulting from the self-consistent solution of the Dyson equation removes the divergences pathological to second-order Møller-Plesset (MP2) theory for strong correlations. In the language often used in density functional theory contexts, this means GF2 has a greatly reduced fractional spin error relative to MP2. The natural question then is what effect, if any, does the Dyson summation have on the fractional charge error in GF2? To this end, we generalize our previous implementation of GF2 to open-shell systems and analyze its fractional spin and charge errors. We find that like MP2, GF2 possesses only a very small fractional charge error, and consequently minimal many electron self-interaction error. This shows that GF2 improves on the critical failings of MP2, but without altering the positive features that make it desirable. Furthermore, we find that GF2 has both less fractional charge and fractional spin errors than typical hybrid density functionals as well as random phase approximation with exchange.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Transition voltages of vacuum-spaced and molecular junctions with Ag and Pt electrodes

The transition voltage of vacuum-spaced and molecular junctions constructed with Ag and Pt electrodes is investigated by non-equilibrium Green's function formalism combined with density functional theory. Our calculations show that, similarly to the case of Au-vacuum-Au previously studied, the transition voltages of Ag and Pt metal-vacuum-metal junctions with atomic protrusions on the electrode surface are determined by the local density of states of the p-type atomic orbitals of the protrusion. Since the energy position of the Pt 6p atomic orbitals is higher than that of the 5p/6p of Ag and Au, the transition voltage of Pt-vacuum-Pt junctions is larger than that of both Ag-vacuum-Ag and Au-vacuum-Au junctions. When one moves to analyzing asymmetric molecular junctions constructed with biphenyl thiol as central molecule, then the transition voltage is found to depend on the specific bonding site for the sulfur atom in the thiol group. In particular agreement with experiments, where the largest transition voltage is found for Ag and the smallest for Pt, is obtained when one assumes S binding at the hollow-bridge site on the Ag/Au(111) surface and at the adatom site on the Pt(111) one. This demonstrates the critical role played by the linker-electrode binding geometry in determining the transition voltage of devices made of conjugated thiol molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Theory of tunneling between two-dimensional electron layers driven by spin pumping: Adiabatic regime and beyond

Tunneling spectroscopy between parallel two-dimensional (2D) electronic systems provides a powerful method to probe the underlying electronic properties by measuring tunneling conductance. Here, in this work, we present a theoretical framework for spin transport in 2D-to-2D tunneling systems, driven by spin pumping. This theory applies to a vertical heterostructure where two layers of metallic 2D electron systems are separated by an insulating barrier, with one layer exchange coupled to a magnetic layer driven at resonance. Utilizing a nonperturbative Floquet-Keldysh formalism, we derive general expressions for the tunneling spin and charge currents across a broad range of driving frequencies, extending beyond the traditional adiabatic pumping regime. At low frequencies, we obtain analytical results that recover the known behaviors in the adiabatic regime. However, at higher frequencies, our numerical findings reveal significant deviations in the dependence of spin and charge currents on both frequency and precession angle. This work offers fresh insights into the role of magnetization dynamics in tunneling transport, opening up new avenues for exploring nonadiabatic spin pumping phenomena.

Green's function methods↗

Helical Organic and Inorganic Polymers

Despite being a staple of synthetic plastics and biomolecules, helical polymers are scarcely studied with Gaussian-basis-set ab initio electron-correlated methods on an equal footing with molecules. This article introduces an ab initio second-order many-body Green’s-function [MBGF(2)] method with nondiagonal, frequency-dependent Dyson self-energy for infinite helical polymers using screw-axis-symmetry-adapted Gaussian-spherical-harmonics basis functions. Together with the Gaussian-basis-set density-functional theory for energies, analytical atomic forces, translational-period force, and helical-angle force, it can compute correlated energy, quasiparticle energy bands, structures, and vibrational frequencies of an infinite helical polymer, which smoothly converge at the corresponding oligomer results. These methods can handle incommensurable structures, which have an infinite translational period and are hard to characterize by any other method, just as efficiently as commensurable structures. We apply them to polyethylene (2/1 helix), polyacetylene (Peierls’ system), and polytetrafluoroethylene (13/6 helix) to establish the quantitative accuracy of MBGF(2)/cc-pVDZ in simulating their (angle-resolved) ultraviolet photoelectron spectra, and of B3LYP/cc-pVDZ or 6-31G** in reproducing their structures, infrared and Raman band positions, phonon dispersions, and (coherent and incoherent) inelastic neutron scattering spectra. We then predict the same properties for infinitely catenated chains of nitrogen or oxygen and discuss their possible metastable existence under ambient conditions. They include planar zigzag polyazene (N2)x (Peierls’ system), 11/3-helical isotactic polyazane (NH)x, 9/4-helical isotactic polyfluoroazane (NF)x, and 7/2-helical polyoxane (O)x as potential high-energy-density materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Controlling Dzyaloshinskii-Moriya interaction in a centrosymmetric nonsymmorphic crystal

Presence of the Dzyaloshinskii-Moriya (DM) interaction in limited noncentrosymmet ric materials leads to novel spin textures and exotic chiral physics. Emergence of DM interaction in centrosymmetric crystals could greatly enrich material realization. Here we show that an itinerant centrosymmetric crystal respecting a nonsymmorphic space group is a new platform for the DM interaction. Taking P4/nmm space group as an example, we demonstrate that the Ruderman-Kittel-Kasuya-Yosida (RKKY) in teraction induces the DM interactions, in addition to the Heisenberg exchange and the Kaplan-Shekhtman-Entin-wohlman-Aharony (KSEA) interaction. Here, the direction of DM vector depends on the positions of magnetic atoms in the real space, and the amplitude depends on the location of the Fermi surface in the reciprocal space. The diversity stems from the position-dependent site groups and the momentum-dependent electronic structures guaranteed by the nonsymmorphic symmetries. Our study unveils the role of the nonsymmorphic symmetries in affecting magnetism, and suggests that the nonsymmorphic crystals can be promising platforms to design magnetic interactions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗