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 1,009 records · Page 56

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↗

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↗

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↗

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↗

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↗

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↗

Systematically improvable multiscale solver for correlated electron systems

The development of numerical methods capable of simulating realistic materials with strongly correlated electrons, with controllable errors, is a central challenge in quantum many-body physics. Here we describe a framework for a general multiscale method based on embedding a self-energy of a strongly correlated subsystem into a self-energy generated by a method able to treat large weakly correlated systems approximately. As an example, we present the embedding of an exact diagonalization self-energy into a self-energy generated from self-consistent second-order perturbation theory. Using a quantum impurity model, generated from a cluster dynamical mean field approximation to the two-dimensional Hubbard model, as a benchmark, we illustrate that our method allows us to obtain accurate results at a fraction of the cost of typical Monte Carlo calculations. We test the method in multiple regimes of interaction strengths and dopings of the model. Here, the general embedding framework we present avoids difficulties such as double counting corrections, frequency-dependent interactions, or vertex functions. As it is solely formulated at the level of the single-particle Green's function, it provides a promising route for the simulation of realistic materials that are currently difficult to study with other methods.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

The effects of focusing on the efficiency of coherent anti-Stokes Raman scattering

The efficiency of coherent anti-Stokes Raman scattering (CARS) is calculated as a function of the input beam sizes of the Stokes and pump laser beams. A Green's function formalism is used to show that maximum efficiency occurs when the confocal parameters of the two input beams are unequal and when the phase mismatch between the anti-Stokes field and the nonlinear polarization is nonzero.

Guha, S.↗

New correction procedures for the fast field program which extend its range

A fast field program (FFP) algorithm was developed based on the method of Lee et al., for the prediction of sound pressure level from low frequency, high intensity sources. In order to permit accurate predictions at distances greater than 2 km, new correction procedures have had to be included in the algorithm. Certain functions, whose Hankel transforms can be determined analytically, are subtracted from the depth dependent Green's function. The distance response is then obtained as the sum of these transforms and the Fast Fourier Transformation (FFT) of the residual k dependent function. One procedure, which permits the elimination of most complex exponentials, has allowed significant changes in the structure of the FFP algorithm, which has resulted in a substantial reduction in computation time.

West, M.↗

Time-dependent integral transport in one-dimensional infinite media using dimensionless variables and the reduced collision formulation

A new method for calculating time-dependent integral neutron transport single collision fluxes has been developed for one-dimensional infinite media problems using dimensionless variables and the reduced collision formulation. By transforming the integral scalar flux equation into dimensionless variables, one may immediately integrate over the time domain and reduce the problem complexity to integrations strictly over the spatial domain, significantly reducing the computational work required to obtain a solution. The results have been compared to published benchmark solutions for validation. Further, reformulating the scalar flux as a Green’s function has allowed for the determination of an analytical expression for the relative error of any collision level scalar flux to the summed scalar flux. This provides the absolute number of collisions required for convergence at any given mean free time, using previous work in one-dimensional Cartesian transport and derived work herein for the one-dimensional spherical transport.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

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↗

Probing the Damage Zone at Parkfield

Rocks are heterogeneous materials that exhibit nonlinear elastic (anelastic) behavior at scales ranging from the laboratory to Earth. In the laboratory, typical, complex relationships exist between stress and strain that include hysteresis, finite relaxation times, strain rate, and history dependence. These behaviors are linked to important characteristics such as stress, porosity, permeability, material integrity, and material failure. In this work, we adopted a “pump-probe” type experiment common in laboratory studies, using solid earth tides as the low-frequency pump and empirical Green's function as the high-frequency probe. By probing the velocity at different points in the pump cycle, we constrained important information about the strain-modulus relationship. Near the San Andreas Fault, we observed strongly nonlinear elastic behavior that characterizes the damage zone. We also constrained important aspects of hysteretic behavior that are related to damage properties and possibly pore pressure. Away from the fault, the nonlinear behavior is diminished.

58 GEOSCIENCES↗

Quantized charge transport in chiral Majorana edge modes

Majorana fermions can be realized as quasiparticles in topological superconductors, with potential applications in topological quantum computing. Recently, lattices of magnetic adatoms deposited on the surface of 𝑠-wave superconductors—Shiba lattices—have been proposed as a new platform for topological superconductivity. Furthermore, these systems possess the great advantage that they are accessible via scanning-probe techniques and thus enable the local manipulation and detection of Majorana modes. Using a nonequilibrium Green's function technique we demonstrate that the topological Majorana edge modes of nanoscopic Shiba islands display universal electronic and transport properties. Most remarkably, these Majorana modes possess a quantized charge conductance that is proportional to the topological Chern number, 𝒞, and carry a supercurrent whose chirality reflects the sign of 𝒞. These results establish nanoscopic Shiba islands as promising components in future topology-based devices.

Edge states↗

Theory of subcycle time-resolved photoemission: Application to terahertz photodressing in graphene

Here, motivated by recent experimental progress we revisit the theory of pump–probe time- and angle-resolved photoemission spectroscopy (trARPES), which is one of the most powerful techniques to trace transient pump-driven modifications of the electronic properties. The pump-induced dynamics can be described in different gauges for the light–matter interaction. Standard minimal coupling leads to the velocity gauge, defined by linear coupling to the vector potential. In the context of tight-binding (TB) models, the Peierls substitution is the commonly employed scheme for single-band models. Multi-orbital extensions – including the coupling of the dipole moments to the electric field – have been introduced and tested recently. In this work, we derive the theory of time-resolved photoemission within both gauges from the perspective of nonequilibrium Green’s functions. This approach naturally incorporates the photoelectron continuum, which allows for a direct calculation of the observable photocurrent. Following this route we introduce gauge-invariant expressions for the time-resolved photoemission signal. The theory is applied to graphene pumped with short terahertz pulses, which we treat within a first-principles TB model. We investigate the gauge invariance and discuss typical effects observed in subcycle time-resolved photoemission. Our formalism is an ideal starting point for realistic trARPES simulations including scattering effects.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Many-Body State Description of Single-Molecule Electroluminescence Driven by a Scanning Tunneling Microscope

Electron transport and optical properties of a single molecule in contact with conductive materials have attracted considerable attention because of their scientific importance and potential applications. With the recent progress in experimental techniques, especially by virtue of scanning tunneling microscope (STM)-induced light emission, where the tunneling current of the STM is used as an atomic-scale source for induction of light emission from a single molecule, it has become possible to investigate single-molecule properties at subnanometer spatial resolution. Despite extensive experimental studies, the microscopic mechanism of electronic excitation of a single molecule in STM-induced light emission has yet to be clarified. Here we present a formulation of single-molecule electroluminescence driven by electron transfer between a molecule and metal electrodes based on a many-body state representation of the molecule. The effects of intramolecular Coulomb interaction on conductance and luminescence spectra are investigated using the nonequilibrium Hubbard Green’s function technique combined with first-principles calculations. We compare simulation results with experimental data and find that the intramolecular Coulomb interaction is crucial for reproducing recent experiments for a single phthalocyanine molecule. The developed theory provides a unified description of the electron transport and optical properties of a single molecule in contact with metal electrodes driven out of equilibrium, and thereby, it contributes to a microscopic understanding of optoelectronic conversion in single molecules on solid surfaces and in nanometer-scale junctions.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Ab Initio Multiplet-Plus-Cumulant Approach for Correlation Effects in X-Ray Photoelectron Spectroscopy

The treatment of electronic correlations in open-shell systems is among the most challenging problems of condensed matter theory. Current approximations are only partly successful. Ligand-field multiplet theory has been widely successful in describing intra-atomic correlation effects in x-ray spectra, but typically ignores itinerant states. The cumulant expansion for the one-electron Green’s function has been successful in describing shake-up effects but ignores atomic multiplets. More complete methods, such as dynamic mean-field theory can be computationally demanding. Here, we show that separating the dynamic Coulomb interactions into local and longer-range parts with ab initio parameters yields a combined multiplet-plus-cumulant approach that accounts for both local atomic multiplets and satellite excitations. The approach is illustrated in transition metal oxides and explains the multiplet peaks, charge-transfer satellites, and distributed background features observed in XPS experiment.

36 MATERIALS SCIENCE↗

Magnus method for electronic structure calculations at extreme conditions

We present the application of Magnus based methods to the solution of first order coupled ordinary differential equations in High Energy Density (HED) physics applications. Our focus is on the application to quantum mechanical methods, specifically on the solution of the radial Dirac equation for real and complex energies. HED applications require accurate solutions across a wide range of spatial and energy domains, including regimes where the solutions exhibit pronounced oscillatory behavior. Such cases pose significant computational challenges. We demonstrate that Magnus-based integrators can efficiently and accurately address these challenges. We discuss the implementation of the Magnus method for the solution of the radial Dirac equation, including practical considerations such as the evaluation of matrix exponentials, numerical integration, error estimation, and adaptive step size control. We also discuss the application of these methods to complex energy Green’s function techniques and the efficient approximation of integrals of the solutions relevant to HED electronic structure calculations. Here, we demonstrate the accuracy and robustness of the resulting method in applications to the free-particle case, for which analytic solutions are available for comparison, as well as the challenging case of gold at HED conditions.

general physics↗

Molecular Optoelectronics (Final Technical Report)

The project developed new theoretical tools for analysis of nanodevices subjected to external optical and bias driving. Traditional theoretical techniques for molecular optoelectronics are either based on kinetic schemes or utilize quasiparticle description. The former fails to account for quantum mechanics of the system-bath interface and may fail qualitatively when used in open interacting systems. The latter is inconvenient in (usual for molecular optical spectroscopy) many body states description of the system. We developed a new technique - the Hubbard nonequilibrium Green's functions method - which allows to avoid shortcomings of standard theoretical approaches. The Hubbard NEGF was applied to analysis of experimental measurements of photoinduced current and bias-induced electroluminescence in STM molecular junctions. We also demonstrated practical usefulness of the Hubbard NEGF within newly proposed universal theory of current-induced forces for nonadiabatic molecular dynamics which generalizes celebrated Head-Gordon and Tully expression for electronic friction. Finally, following latest trends in laser techniques we studied local responses (local fluxes and local noise spectroscopy) of biased junctions as a source of information not accessible by total responses (total fluxes and noises) of the system.

30 DIRECT ENERGY CONVERSION↗