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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Density functional study of relaxation of adsorbate vibration modes: Dominance of anharmonic interaction

Formulation and density functional workflow for calculating the lifetime of vibrational modes of molecular adsorbates on solid surfaces due to vibration–phonon coupling are presented. The anharmonic coupling is invoked to give the correct description of the origin of temperature dependence. Using pyrrolidine (C 4 H 9 N) absorbed on the Cu(001) surface as a concrete example, here we show that the anharmonic coupling can be one to two orders more significant than the harmonic interaction for the broadening of vibrational spectra, especially as temperature increases. These results challenge the common assumption that the anharmonic interaction is weak and call for attention of considering its effect in quantum relaxation and related problems.

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

Roadmap on methods and software for electronic structure based simulations in chemistry and materials

This Roadmap article provides a succinct, comprehensive overview of the state of electronic structure methods and software for molecular and materials simulations. Seventeen distinct sections collect insights by 51 leading scientists in the field. Each contribution addresses the status of a particular area, as well as current challenges and anticipated future advances, with a particular eye towards software related aspects and providing key references for further reading. Foundational sections cover density functional theory and its implementation in real-world simulation frameworks, Green's function based many-body perturbation theory, wave-function based and stochastic electronic structure approaches, relativistic effects and semiempirical electronic structure theory approaches. Subsequent sections cover nuclear quantum effects, real-time propagation of the electronic structure, challenges for computational spectroscopy simulations, and exploration of complex potential energy surfaces. The final sections summarize practical aspects, including computational workflows for complex simulation tasks, the impact of current and future high-performance computing architectures, software engineering practices, education and training to maintain and broaden the community, as well as the status of and needs for electronic structure based modeling from the vantage point of industry environments. Overall, the field of electronic structure software and method development continues to unlock immense opportunities for future scientific discovery, based on the growing ability of computations to reveal complex phenomena, processes and properties that are determined by the make-up of matter at the atomic scale, with high precision.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Floquet dynamics of disordered bands with isolated critical energies

We investigate localization properties of driven models which exhibit a subextensive number of extended states in the static setting. We consider instances where the extended modes are or are not protected by topological considerations. To this end, we contrast the strongly driven, disordered, lowest Landau level, which we refer to as the random Landau model (RLM), with the random dimer model (RDM); the RDM also has a subextensive set of delocalized modes in the middle of the spectrum whose origin is not topological. We map the driven models on to a higher-dimensional effective model and numerically compute the localization length as a function of disorder strength, drive amplitude, and frequency using the recursive Green's function method. Our numerical results indicate that, in the presence of a strong drive (low frequency and/or large drive amplitude), the topologically protected RLM continues to exhibit a spectrum with both localized and delocalized (or critical) modes, but the spectral range of delocalized modes is enhanced by the driving. This occurs due to an admixture of the localized modes with extended modes arising due to the topologically protected critical energy in the middle of the spectrum. On the other hand, in the RDM, a weak drive immediately localizes the entire spectrum. This occurs in contrast to the naive expectation from perturbation theory that mixing between localized and delocalized modes generically enhances the delocalization of all modes. Our work highlights the importance of the origin of the delocalized modes in the localization properties of the corresponding Floquet model.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Unusual Exchange Couplings and Intermediate Temperature Weyl State in Co 3 Sn 2 S 2

Understanding magnetism and its possible correlations to topological properties has emerged to the forefront as a difficult topic in studying magnetic Weyl semimetals. Co 3 Sn 2 S 2 is a newly discovered magnetic Weyl semimetal with a kagome lattice of cobalt ions and has triggered intense interest for rich fantastic phenomena. In this study, we report the magnetic exchange couplings of Co 3 Sn 2 S 2 using inelastic neutron scattering and two density functional theory (DFT) based methods: constrained magnetism and multiple-scattering Green’s function methods. Co 3 Sn 2 S 2 exhibits highly anisotropic magnon dispersions and linewidths below T C , and paramagnetic excitations above T C . The spin-wave spectra in the ferromagnetic ground state is well described by the dominant third-neighbor “across-hexagon” J d model. Our density functional theory calculations reveal that both the symmetry-allowed 120° antiferromagnetic orders support Weyl points in the intermediate temperature region, with distinct numbers and the locations of Weyl points. Our study highlights the important role Co 3 Sn 2 S 2 can play in advancing our understanding of kagome physics and exploring the interplay between magnetism and band topology.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Identifying topological superconductivity in two-dimensional transition-metal dichalcogenides

Here we study the superconducting pairing instabilities and gap functions for prototypical two-dimensional (2D) transition-metal dichalcogenides (TMDCs) WS 2 , MoTe 2 , and MoS 2 in the 2H phase under both hole and electron doping at 10 K. Our first-principles quantum many-body Green's function approach allows us to treat the full d and p manifold of orbitals with strong spin-orbit coupling, yielding pairing predictions with material-specific detail. The resulting gap functions exhibit a variety of mixed-parity superconducting states, including s, p, d, f, d ± id, and p ± ip pairing modes. In particular, we predict 3% and 4% hole-doped WS 2 to be a chiral p ± ip topological superconductor. For 1% hole-doped MoS 2 , we find a competition between three doubly degenerate chiral and nonchiral instabilities. Overall, the relative pairing strengths are found to follow the Fermi-surface topology, due to nesting between the Fermi-surface sheets. Finally, we discuss our predictions in relation to available experimental data and classify the topology of the predicted superconducting pairing symmetries.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Cumulant methods for electron-phonon problems. II. The self-consistent cumulant expansion

Here in this paper we present a self-consistent cumulant expansion (SC-CE) and investigate its accuracy for the one-dimensional Holstein model with and without phonon dispersion. We show that for finite lattice sizes, the numerical integration of the SC-CE equations becomes unstable at long times. This defect is partially ameliorated when studying systems in the thermodynamic limit, enabling the demonstration that the SC-CE corrects many deficits of the standard perturbative CE in the (nondispersive) Holstein model. The natural phonon damping that arises in the more realistic dispersed Holstein model renders the SC-CE stable, allowing for a complete assessment of the method. Here, we find that self-consistency dramatically corrects many of the failures found in the perturbative CE, but also introduces some unphysical features. Finally, we comment on the potential use of SC-CE as a tool for calculating Green's functions in generic many-body problems.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A theory for colors of strongly correlated electronic systems

Many strongly correlated transition metal insulators are colored, even though they have band gaps much larger than the highest energy photons from the visible light. An adequate explanation for the color requires a theoretical approach able to compute subgap excitons in periodic crystals, reliably and without free parameters—a formidable challenge. The literature often fails to disentangle two important factors: what makes excitons form and what makes them optically bright. We pick two archetypal cases as examples: NiO with green color and MnF 2 with pink color, and employ two kinds of ab initio many body Green’s function theories; the first, a perturbative theory based on low-order extensions of the $GW$ approximation, is able to explain the color in NiO, while the same theory is unable to explain why MnF 2 is pink. We show its color originates from higher order spin-flip transitions that modify the optical response, which is contained in dynamical mean-field theory (DMFT). We show that symmetry lowering mechanisms may determine how ‘bright’ these excitons are, but they are not fundamental to their existence.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Effect of surface functional groups on MXene conductivity

In this work, we report the in-plane electron transport in the MXenes (i.e., within the MXene layers) as a function of composition using the density-functional tight-binding method, in conjunction with the non-equilibrium Green’s functions technique. Our study reveals that all MXene compositions have a linear relationship between current and voltage at lower potentials, indicating their metallic character. However, the magnitude of the current at a given voltage (conductivity) has different trends among different compositions. For example, MXenes without any surface terminations (Ti 3 C 2 ) exhibit higher conductivity compared to MXenes with surface functionalization. Among the MXenes with –O and –OH termination, those with –O surface termination have lower conductivity than the ones with –OH surface terminations. Interestingly, conductivity changes with the ratio of –O and –OH on the MXene surface. Our calculated I–V curves and their conductivities correlate well with transmission functions and the electronic density of states around the Fermi level. The surface composition-dependent conductivity of the MXenes provides a path to tune the in-plane conductivity for enhanced pseudocapacitive performance.

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

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↗
DOE OSTI2021

Fluctuations analysis of spin susceptibility: Néel ordering revisited in dynamical mean field theory

Here in this paper, we revisit the antiferromagnetic (AF) phase diagram of the single-band three-dimensional half-filled Hubbard model on a simple cubic lattice studied within the dynamical mean field theory (DMFT). Although this problem has been investigated extensively in the literature, a comprehensive understanding of the impact of the different one- and, in particular, two-particle local correlation functions of DMFT on the AF transition temperature is still missing. We have, hence, performed a fluctuation analysis of T N with respect to different local bosonic fluctuations (charge, spin, particle-particle) contained in the two-particle vertex of DMFT. Our results indicate that, beyond weak coupling, the screening of the DMFT vertex by local fluctuations leads to an enhancement of T N with respect to a random phase approximation (RPA) like calculation where this vertex is replaced by the bare interaction. The overall suppression of T N in DMFT with respect to RPA is then solely due to the incoherence introduced by the DMFT self-energy in the one-particle Green's functions. This illustrates the Janus-faced role of the local moment formation in the DMFT solution of the Hubbard model, which leads to completely opposite effects in the one- and two-particle correlation functions.

36 MATERIALS SCIENCE↗
DOE OSTI2024

An asymptotic-preserving semi-Lagrangian algorithm for the anisotropic heat transport equation with arbitrary magnetic fields

Here, we extend the recently proposed semi-Lagrangian algorithm for the extremely anisotropic heat transport equation [Chacón et al., J. Comput. Phys ., 272 (2014)] to deal with arbitrary magnetic field topologies. The original scheme (which showed remarkable numerical properties) was valid for the so-called tokamak-ordering regime, in which the magnetic field magnitude was not allowed to vary much along field lines. The proposed extension maintains the attractive features of the original scheme (including the analytical Green's function, which is critical for tractability) with minor modifications, while allowing for completely general magnetic fields. The accuracy and generality of the approach are demonstrated by numerical experiment with an analytical manufactured solution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Minimal two-body quantum absorption refrigerator

Here, we study the phenomenon of absorption refrigeration, where refrigeration is achieved by heating instead of work, in two different setups: a minimal set up based on coupled qubits and two nonlinearly coupled resonators. Considering ZZ interaction between the two qubits, we outline the basic ingredients required to achieve cooling. Using local as well as global master equations, we observe that inclusion of an XX type term in the qubit-qubit coupling is detrimental to cooling. We compare the cooling effect obtained in the qubit case with that of nonlinearly coupled resonators (multilevel system) where the ZZ interaction translates to a Kerr-type nonlinearity. For small to intermediate strengths of nonlinearity, we observe that multilevel quantum systems, for example qutrits, give better cooling effect compared to the qubits. Using Keldysh nonequilibrium Green's function formalism, we go beyond first order sequential tunneling processes and study the effect of higher order processes on refrigeration. We find a reduced cooling effect compared to the master equation calculations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Full-potential KKR within the removed-sphere method: A practical and accurate solution to the Poisson equation

Here an efficient and accurate generalization of the removed-sphere method (RSM) to solve the Poisson equation for total charge density in a solid with space-filling convex Voronoi polyhedra (VPs) and any symmetry is presented. The generalized RSM avoids the use of multipoles and VP shape functions for cellular integrals, which have associated ill-convergent large, double-internal L sums in spherical-harmonic expansions, so that fast convergence in single-L sums is reached. Our RSM adopts full Ewald formulation to work for all configurations or when symmetry breaking occurs, such as for atomic displacements or elastic constant calculations. The structure-dependent coefficients A L that define RSM can be calculated once for a fixed structure and speed up the whole self-consistent-field procedure. The accuracy and rapid convergence properties are confirmed using two analytic models, including the Coulomb potential and energy. We then implement the full-potential RSM using the Green's function Korringa-Kohn-Rostoker (KKR) method for real applications and compare the results with other first-principle methods and experimental data, showing that they are equally as accurate.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Variational Discrete Action Theory

In this work, we propose the variational discrete action theory (VDAT) to study the ground state properties of quantum many-body Hamiltonians. VDAT is a variational theory based on the sequential product density matrix (SPD) ansatz, characterized by an integer $\mathscr{N}$, which monotonically approaches the exact solution with increasing $\mathscr{N}$. To evaluate the SPD, we introduce a discrete action and a corresponding integer time Green’s function. We use VDAT to exactly evaluate the SPD in two canonical models of interacting electrons: the Anderson impurity model and the d = ∞ Hubbard model. For the latter, we evaluate $\mathscr{N}$ = 2 – 4, where $\mathscr{N}$ = 2 recovers the Gutzwiller approximation (GA), and we show that $\mathscr{N}$ = 3, which exactly evaluates the Gutzwiller-Baeriswyl wave function, provides a truly minimal yet precise description of Mott physics with a cost similar to that of the GA. VDAT is a flexible theory for studying quantum Hamiltonians, competing both with state-of-the-art methods and simple, efficient approaches all within a single framework.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Ab initio approaches to high-entropy alloys: a comparison of CPA, SQS, and supercell methods

In this work, we present a comparative study of different modeling approaches to the electronic properties of the Hf 0.05 Nb 0.05 Ta 0.8 Ti 0.05 Zr 0.05 high-entropy alloy. Common to our modeling is the methodology to compute the one-particle Green’s function in the framework of density functional theory. We demonstrate that the special quasi-random structures modeling and the supercell, i.e., the locally self-consistent multiple-scattering methods, provide very similar results for the ground state properties such as the spectral function (density of states) and the equilibrium lattice parameter. To reconcile the multiple-scattering single-site coherent potential approximation with the real space supercell methods, we included the effect of screening of the net charges of the alloy components. Based on the analysis of the total energy and spectral functions computed within the density functional theory, we found no signature for the long-range or local magnetic moments formation in the Hf 0.05 Nb 0.05 Ta 0.8 Ti 0.05 Zr 0.05 high-entropy alloy; instead, we find possible superconductivity below ~9K.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Sparsity of Radiating Characteristic Modes on Infinite Periodic Structures

Characteristic modes on infinite periodic structures are studied using spectral dyadic Green’s functions. This formulation demonstrates that, in contrast to the modal analysis of finite structures, the number of radiating characteristic modes is limited by unit cell size and incident wave vector (i.e., scan angle or phase shift per unit cell). Here, the reflection tensor is decomposed into modal contributions from radiating modes, indicating that characteristic modes are a predictably sparse basis in which to study reflection phenomena.

42 ENGINEERING↗
DOE OSTI2020

DFTB+, a software package for efficient approximate density functional theory based atomistic simulations

DFTB+ is a versatile community developed open source software package offering fast and efficient methods for carrying out atomistic quantum mechanical simulations. By implementing various methods approximating density functional theory (DFT), such as the density functional based tight binding (DFTB) and the extended tight binding method, it enables simulations of large systems and long timescales with reasonable accuracy while being considerably faster for typical simulations than the respective ab initio methods. Based on the DFTB framework, it additionally offers approximated versions of various DFT extensions including hybrid functionals, time dependent formalism for treating excited systems, electron transport using non-equilibrium Green’s functions, and many more. DFTB+ can be used as a user-friendly standalone application in addition to being embedded into other software packages as a library or acting as a calculation-server accessed by socket communication. We give an overview of the recently developed capabilities of the DFTB+ code, demonstrating with a few use case examples, discuss the strengths and weaknesses of the various features, and also discuss on-going developments and possible future perspectives.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Tensor train continuous time solver for quantum impurity models

The simulation of strongly correlated quantum impurity models is a significant challenge in modern condensed matter physics that has multiple important applications. Thus far, the most successful methods for approaching this challenge involve Monte Carlo techniques that accurately and reliably sample perturbative expansions to any order. However, the cost of obtaining high precision through these methods is high. Recently, tensor train decomposition techniques have been developed as an alternative to Monte Carlo integration. In this study, we apply these techniques to the single-impurity Anderson model at equilibrium by calculating the systematic expansion in power of the hybridization of the impurity with the bath. Furthermore, we demonstrate the performance of the method in a paradigmatic application, examining the first-order phase transition on the infinite-dimensional Bethe lattice, which can be mapped to an impurity model through dynamical mean field theory. Our results indicate that using tensor train decomposition schemes allows the calculation of finite-temperature Green's functions and thermodynamic observables with unprecedented accuracy. The methodology holds promise for future applications to frustrated multiorbital systems, using a combination of partially summed series with other techniques pioneered in diagrammatic and continuous time quantum Monte Carlo.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗