Generalized quadrature for finite temperature Green’s function methods
Not provided.
SEARCH · Search NASA
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.
Not provided.
The Matsubara Green’s function that is used to describe temperature-dependent behavior is expressed on a numerical grid. While such a grid usually has a couple of hundred points for low-energy model systems, for realistic systems with large basis sets the size of an accurate grid can be tens of thousands of points, constituting a severe computational and memory bottleneck. In this paper, we determine efficient imaginary time grids for the temperature-dependent Matsubara Green’s function formalism that can be used for calculations on realistic systems. Here we show that, because of the use of an orthogonal polynomial transform, we can restrict the imaginary time grid to a few hundred points and reach micro-Hartree accuracy in the electronic energy evaluation. Moreover, we show that only a limited number of orthogonal polynomial expansion coefficients are necessary to preserve accuracy when working with a dual representation of the Green’s function or self-energy and transforming between the imaginary time and frequency domain.
Here, we have analyzed the electron-phonon interaction in the Periodic Anderson Model (PAM) to describe the temperature dependence of resistivity in some heavy fermion (HF) systems for finite wave vector (q) and for finite temperature (T). Since the resistivity is related to the imaginary part of the electron self energy, the expression for the same is evaluated through double time temperature dependant Green function technique of the Zubarev type. The effect of different system parameters namely the position of 4f level, E{sub 0} and the electron - phonon coupling strengths on resistivity have been studied. The results obtained give satisfactory explanations to the experimental observations.
The total energy shift of a polarized massive Dirac neutrino in an electron–positron plasma in a constant magnetic field is investigated. The calculation in the Feynman gauge is performed for the first time by using the Matsubara temperature Green functions. The dependence of the dispersion relation and the anomalous magnetic moment of the neutrino on the magnetic field strength, spin, energy, direction of motion, neutrino mass, and the plasma parameters is analyzed. The results of investigations for the massive neutrino in the limiting case are compared with those obtained earlier by other authors for the massless left neutrino.
Including finite-temperature effects from the electronic degrees of freedom in electronic structure calculations of semiconductors and metals is desired; however, in practice it remains exceedingly difficult when using zero-temperature methods, since these methods require an explicit evaluation of multiple excited states in order to account for any finite-temperature effects. Using a Matsubara Green’s function formalism remains a viable alternative, since in this formalism it is easier to include thermal effects and to connect the dynamic quantities such as the self-energy with static thermodynamic quantities such as the Helmholtz energy, entropy, and internal energy. However, despite the promising properties of this formalism, little is known about the multiple solutions of the non-linear equations present in the self-consistent Matsubara formalism and only a few cases involving a full Coulomb Hamiltonian were investigated in the past. In this study, to shed some light onto the iterative nature of the Green’s function solutions, we self-consistently evaluate the thermodynamic quantities for a one-dimensional (1D) hydrogen solid at various interatomic separations and temperatures using the self-energy approximated to second-order (GF2). At many points in the phase diagram of this system, multiple phases such as a metal and an insulator exist, and we are able to determine the most stable phase from the analysis of Helmholtz energies. Additionally, we show the evolution of the spectrum of 1D boron nitride to demonstrate that GF2 is capable of qualitatively describing the temperature effects influencing the size of the band gap.
Modeling high-temperature (tens or hundreds of eV), dense plasmas is challenging due to the multitude of non-negligible physical effects including significant partial ionization and multisite effects. These effects cause the breakdown or intractability of common methods and approximations used at low temperatures, such as pseudopotentials or plane-wave basis sets. Here we explore the Korringa-Kohn-Rostoker Green's function method at these high-temperature conditions. The method is all electron, does not rely on pseudopotentials, and uses a spherical harmonic basis set, and so avoids the aforementioned limitations. Finally, it is found to be accurate for solid density aluminum and iron plasmas when compared to a plane-wave method at low temperature, while being able to access high temperatures.
A new technique has been developed to study two dimensional heat transfer problems in gears. This technique consists of transforming the heat equation into a line integral equation with the use of Green's theorem. The equation is then expressed in terms of eigenfunctions that satisfy the Helmholtz equation, and their corresponding eigenvalues for an arbitrarily shaped region of interest. The eigenfunction are obtalned by solving an intergral equation. Once the eigenfunctions are found, the temperature is expanded in terms of the eigenfunctions with unknown time dependent coefficients that can be solved by using Runge Kutta methods. The time integration is extremely efficient. Therefore, any changes in the time dependent coefficients or source terms in the boundary conditions do not impose a great computational burden on the user. The method is demonstrated by applying it to a sample gear tooth. Temperature histories at representative surface locatons are given.
Thermodynamic properties of the interacting homogeneous electron gas are calculated using a finite-temperature cumulant Green's function approach over a broad range of densities and temperatures up to the warm dense matter regime $T ~ T_F$, where $T_F$ is the Fermi degeneracy temperature. These properties can be separated into independent particle and exchange-correlation contributions, and our focus here is on the latter. In this work our approach is based on the Galitskii-Migdal-Koltun and electron number sum rules from the finite temperature many-body Green's function formalism, together with an extension of the cumulant Green's function to finite temperature. Previously this approach yielded exchange-correlation energies and potentials in good agreement with quantum Monte-Carlo calculations. Here the method is extended for various thermodynamic quantities including the chemical potential, total energy, Helmholtz free-energy, electronic equation of state, specific heat, and isothermal compressibility, which optionally include spin dependence. We find that the exchange-correlation contributions are weakly varying at low temperature but exhibit significant temperature dependence in the WDM regime, as well as a crossover from exchange- to correlation-dominated behavior. In contrast to the $T = 0^+$ limit, we also find that renormalization effects are largely but not completely suppressed at finite temperature. Comparisons with other approaches at various levels of approximation are also discussed.
Plasmon and exciton superconductivity mechanisms are discussed. Superconductivity in a three layer metal semiconductor metal and insulator semimetal insulator sandwich structure was described in terms of the temperature dependent Green function of the longitudinal (Coulomb) field. The dependences of the superconducting transition temperature on structure parameters were obtained. In a semiconducting film, as a result of interactions of degenerate free carriers with excitons, superconductivity exists only in a certain range of parameter values, and the corresponding critical temperature is much lower than in the plasmon mechanism of superconductivity.
The theory of electron tunneling through an open nanostructure as an active element of a quantum cascade detector is developed, which takes into account the interaction of electrons with confined and interface phonons. Using the method of finite-temperature Green’s functions and the electron-phonon Hamiltonian in the representation of second quantization over all system variables, the temperature shifts and electron-level widths are calculated and the contributions of different electron-phonon-interaction mechanisms to renormalization of the spectral parameters are analyzed depending on the geometrical configuration of the nanosystem. Due to weak electron-phonon coupling in a GaAs/Al{sub 0.34}Ga{sub 0.66}As-based resonant tunneling nanostructure, the temperature shift and rf field absorption peak width are not very sensitive to the electron-phonon interaction and result from a decrease in potential barrier heights caused by a difference in the temperature dependences of the well and barrier band gaps.
In this proceeding, an idea is outlined suggesting a generic treatment on any type of input data for a numerical analytic continuation problem, which is needed when dynamical information is to be extracted from a calculationally convenient one particle imaginary time Green function. The quality of the resulting spectral function will rely only on the data to be treated, viz, data-driven. This is different from the Maximum Entropy or the Stochastic method which relies on an entropy term to guide convergence of the resulting spectral function.
Green function theory of ferromagnetism used to examine implications of cube-root behavior of magnetization near Curie point
Thermodynamic properties of gapped graphene-like structures by considering the effects of interaction between electrons and Holstein phonons have been studied. Particularly, we study the heat capacity and paramagnetic susceptibility of structures as a function of temperature within the Green’s function method with the help of Holstein model. The paramagnetic susceptibility and heat capacity can be derived by using density of states based on the Kubo formula. We have found the energy dependence of density of states for various values of gap in the presence of Holstein phonons. Finally, the temperature behaviors of specific heat and spin susceptibility of gapped graphene structure due to electron-phonon coupling have been investigated. Our results show the electron-phonon interaction leads to the appearance of a double van Hov singularity for each value of gap parameter. Also, electron-phonon coupling affects the value of heat capacity and magnetic susceptibility at low temperatures.
The class of RMn{sub 2}O{sub 5} (R = Ho, Tb, Y, Eu) compounds offers multiferroic properties where the refined magnetic zig-zag order breaks the inversion symmetry. Varying the temperature, the system undergoes a magnetic and a subsequent ferroelectric phase transition where the ferroelectricity is magnetically induced. We propose a modified anisotropic Heisenberg model that can be used as a tractable analytical model studying the properties of those antiferromagnetic zig-zag spin chains. Based on a finite temperature Green's function method, it is shown that the polarization is induced solely by different exchange couplings of the two different Mn{sup 4+} and Mn{sup 3+} magnetic ions. We calculate the excitation energy of the spin system for finite temperatures, which for its part determines the temperature dependent magnetization and polarization. The ferroelectric phase transition is manifested as a kink in the excitation energy. The variation of the polarization by an external magnetic field depends strongly on the direction of that field. Whereas, the polarization in b-direction increases with an external magnetic field as well in b-direction it can be switched for strong fields in a-direction. The results based on that modified Heisenberg model are in qualitative agreement with experimental data.
Here we present a theory of the dc Josephson effect in contacts between Fe-based and spin-singlet s -wave superconductors. The method is based on the calculation of temperature Green's function in the junction within the tight-binding model. We calculate the phase dependencies of the Josephson current for different orientations of the junction relative to the crystallographic axes of Fe-based superconductor. Further, we consider the dependence of the Josephson current on the thickness of an insulating layer and on temperature. Experimental data for PbIn / Ba 1 – x K x ( FeAs ) 2 point-contact Josephson junctions are consistent with theoretical predictions for s ± symmetry of an order parameter in this material. The proposed method can be further applied to calculations of the dc Josephson current in contacts with other new unconventional multiorbital superconductors, such as Sr 2 RuO 4 and the superconducting topological insulator Cu x Bi 2 Se 3 .
In this work, we present algorithmic and implementation details for the fully self-consistent finite-temperature $GW$ method in Gaussian Bloch orbitals for solids. Our implementation is based on the finite-temperature Green's function formalism in which all equations are solved on the imaginary axis, without resorting to analytical continuation during the self-consistency. No quasiparticle approximation is employed and all matrix elements of the self-energy are explicitly evaluated. The method is tested by evaluating the band gaps of selected semiconductors and insulators. We show agreement with other, differently formulated, finite-temperature sc $GW$ implementations when finite-size corrections and basis-set errors are taken into account. By migrating computationally intensive calculations to graphics processing units, we obtain scalable results on large supercomputers with nearly optimal performance. Our work demonstrates the applicability of Gaussian orbital based sc $GW$ for ab initio correlated material simulations and provides a sound starting point for embedding methods built on top of $GW$.
We investigated the Josephson effect between the conventional superconductor Pb0.7In0.3 and Ba0.4K0.6 (FeAs)2 single crystals by performing directional point-contact experiments. The junctions were irradiated with microwaves and the behavior of the Shapiro steps as a function of the square root of the microwave power was studied by means of the RSJ model extended to the nonautonomous case with an rf current-source term. The analysis reveals that, for current injection along the ab plane, the current-phase relation deviates from being simply described by the standard sin(Φ) (as we indeed observed along the c axis) and a dominant sin(2Φ) component is observed. The results are compared with a theory of the dc Josephson effect, formulated by using a method based on the calculation of temperature Green’s function in the junction within the tight-binding model (Burmistrova et al., Phys. Rev. B 91, 214501, 2015). The comparison shows that the experimental results here presented strongly favor an s ±-wave symmetry of the order parameter in this compound.