Moment tensor catalogue of California and Nevada using 3D Green's functions
Explore the source record for details and available documents.
Geosciences↗
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.
Spin-orbit-proximitized ferromagnetic metal by monolayer transition metal dichalcogenide: Atlas of spectral functions, spin textures, and spin-orbit torques in Co/MoSe 2 , Co/WSe 2 and Co/TaSe 2 heterostructures
In this work, the heterostructures composed of an ultrathin ferromagnetic metal (FM) and a material hosting strong spin-orbit (SO) coupling are the principal resource for SO torque and spin-to-charge conversion nonequilibrium effects in spintronics. The key quantity in theoretical description of these effects is nonequilibrium spin density, which can appear on any monolayer of the heterostructure through which the current is flowing as long as the monolayer bands are affected by the native or proximity induced SO coupling. Here we demonstrate how hybridization of wave functions of Co layer and a monolayer of transition metal dichalcogenides (TMDs)—such as semiconducting MoSe 2 and WSe 2 or metallic TaSe 2 —can lead to dramatic transmutation of electronic and spin structure of Co within some distance away from its interface with TMD, when compared to the bulk of Co or its surface in contact with a vacuum. This is due to proximity induced SO splitting of Co bands encoded in the spectral functions and spin textures on its monolayers, which we obtain using noncollinear density functional theory (ncDFT) combined with equilibrium Green function (GF) calculations. In fact, SO splitting is present due to structural inversion asymmetry of the bilayer—i.e., just the presence of the interface—even if SO coupling within TMD monolayer is artificially switched off in ncDFT calculations. However, switching it on makes the effects associated with proximity SO coupling within Co layer about five times larger. Injecting spin-unpolarized charge current through SO-proximitized monolayers of Co generates nonequilibrium spin density over them, so that its cross product with the magnetization of Co determines SO torque. The SO torque computed via first-principles quantum transport methodology, which combines ncDFT with nonequilibrium GF calculations, can be used as the screening parameter to identify optimal combination of materials and their interfaces for applications in spintronics. In particular, we identify heterostructure two-monolayer-Co/monolayer-WSe 2 as the most optimal one, at least in the clean limit.
2-dimensional systems↗
Foundations of variational discrete action theory
Variational wave functions and Green's functions are two important paradigms for solving quantum Hamiltonians, each having their own advantages. Here we detail the variational discrete action theory (VDAT), which exploits the advantages of both paradigms in order to approximately solve the ground state of quantum Hamiltonians. VDAT consists of two central components: the sequential product density matrix (SPD) ansatz and a discrete action associated with the SPD. The SPD is a variational ansatz inspired by the Trotter decomposition and characterized by an integer $\mathscr{N}$, recovering many well-known variational wave functions, in addition to the exact solution for $\mathscr{N}$ = ∞. The discrete action describes all dynamical information of an effective integer time evolution with respect to the SPD. We generalize the path integral to our integer time formalism, which converts a dynamic correlation function in integer time to a static correlation function in a compound space. We also generalize the usual many-body Green's function formalism to integer time, which results in analogous but distinct mathematical structures, yielding integer time versions of the generating functional, Dyson equation, and Bethe-Salpeter equation. We prove that the SPD can be exactly evaluated in the multiband Anderson impurity model (AIM) by summing a finite number of diagrams. For the multiband Hubbard model, we prove that the self-consistent canonical discrete action approximation (SCDA), which is the integer time analog of the dynamical mean-field theory, exactly evaluates the SPD for d = ∞. VDAT within the SCDA provides an efficient yet reliable method for capturing the local physics of quantum lattice models, which will have broad applications for strongly correlated electron materials. More generally, VDAT should find applications in various many-body problems in physics.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Broadband impedance modeling and single bunch instabilities estimations of the advanced light source upgrade project
We report the Advanced Light Source Upgrade (ALS-U) is a 4th generation diffraction-limited soft x-ray radiation source, consisting of a new accumulator ring (AR) and a new storage ring (SR). In both rings coupling-impedance driven instabilities need careful evaluation to ensure meeting the machine’s high-performance goals. This paper presents the workflow followed in building the impedance models and the beam-stability analysis based on those models. We follow the commonly accepted approach of separating the resistive-wall and the geometric parts of the impedance; the former is obtained by analytical formulas, the latter by numerical electro-magnetic codes (primarily CST Studio software) with perfectly-conducting boundary conditions. Impedance budgets are established and pseudo-Green functions calculated to be used in beam dynamics studies. We also present various ways to cross-check simulation results for reliable impedance modeling. Finally, the crucial single-bunch instability current thresholds for various operation modes are determined and discussed.
07 ISOTOPE AND RADIATION SOURCES↗
A computational approach to complex junctions and interfaces
Molecular- and nano- junctions, surfaces, interfaces, and interfacial processes have been the subjects of this DOE supported program for about 20 years. We have investigated electronic and magnetic structure, spin-dependent charge transport, and emerging phenomena at surfaces and interfaces in a number of physical systems ranging from nano-clusters, to nanowires, to two-dimensional (2D) crystals and their interactions, to bulk matter. Eighty-four refereed papers were published [1-84], with a majority effort from the group led by the PI, and over a hundred presentations were made at international conferences. We use first-principles methods in the framework of density functional theory in conjunction with nano-equilibrium Green function techniques, semi-empirical Boltzmann transport theory, tight-binding models, and the Hubbard model and many-body wavefunction approaches. In pursuing our scientific goals, we have also developed a number of computational algorithms and computational packages.
2D materials↗
Approximating and incorporating model uncertainty in an inversion for seismic source functions: Preliminary results
We present preliminary work on propagating model uncertainty into the estimation of the time domain source time functions of the seismic source. Our method is based on an estimated model covariance function, which we estimate from the data. The model covariance function is then used to construct a suite of surrogate Greens functions which we use in a Monte Carlo type inversion scheme. The result is a probability density function of the six independent source time functions, each of which corresponds to an individual component of the seismic moment tensor. We compare the results of our method with those obtained using a computationally expensive finite difference Monte Carlo method and find that our new method produces results that are deficient in low frequencies. The advantage of our new method, which we term the Karhunen-Loeve Monte Carlo (KLMC) method, is that is several orders of magnitude faster than our current method, which uses a finite difference scheme to produce the suite of forward models.
42 ENGINEERING↗
Dynamic mode decomposition for extrapolating nonequilibrium Green's-function dynamics
The Hartree-Fock generalized Kadanoff-Baym ansatz (HF-GKBA) offers an approximate numerical procedure for propagating the two-time nonequilibrium Green's function (NEGF). Here, using the GW self-energy, we compare the HF-GKBA to exact results for a variety of systems with long- and short-range interactions, different two-body interaction strengths, and various nonequilibrium preparations. We find excellent agreement between the HF-GKBA and exact time evolution in models when more realistic long-range exponentially decaying interactions are considered. Furthermore, this agreement persists for long times and for intermediate to strong interaction strengths. In large systems, HF-GKBA becomes prohibitively expensive for long-time evolutions. For this reason, we look at the use of dynamical mode decomposition (DMD) to reconstruct long-time NEGF trajectories from a sample of the initial trajectory. Using no more than 16% of the total time evolution, we reconstruct the total trajectory with high fidelity. Our results show the potential for DMD to be used in conjunction with HF-GKBA to calculate long-time trajectories in large-scale systems.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
GW with hybrid functionals for large molecular systems
A low-cost approach for stochastically sampling static exchange during time-dependent Hartree–Fock-type propagation is presented. This enables the use of an excellent hybrid density functional theory (DFT) starting point for stochastic GW quasiparticle energy calculations. Generalized Kohn–Sham molecular orbitals and energies, rather than those of a local-DFT calculation, are used for building the Green function and effective Coulomb interaction. The use of an optimally tuned hybrid diminishes the starting point dependency in one-shot stochastic GW, effectively avoiding the need for self-consistent GW iterations.
Chemistry↗
A numerical Poisson solver with improved radial solutions for a self-consistent locally scaled self-interaction correction method
Abstract The universal applicability of density functional approximations is limited by self-interaction error made by these functionals. Recently, a novel one-electron self-interaction-correction (SIC) method that uses an iso-orbital indicator to apply the SIC at each point in space by scaling the exchange-correlation and Coulomb energy densities was proposed. The locally scaled SIC (LSIC) method is exact for the one-electron densities, and unlike the well-known Perdew–Zunger SIC (PZSIC) method recovers the uniform electron gas limit of the uncorrected density functional approximation, and reduces to PZSIC method as a special case when isoorbital indicator is set to the unity. Here, we present a numerical scheme that we have adopted to evaluate the Coulomb potential of the electron density scaled by the iso-orbital indicator required for the self-consistent LSIC calculations. After analyzing the behavior of the finite difference method (FDM) and the green function solution to the radial part of the Poisson equation, we adopt a hybrid approach that uses the FDM for the Coulomb potential due to the monopole and the GF for all higher-order terms. The performance of the resultant hybrid method is assessed using a variety of systems. The results show improved accuracy than earlier numerical schemes. We also find that, even with a generic set of radial grid parameters, accurate energy differences can be obtained using a numerical Coulomb solver in standard density functional studies.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
First-principles study of L -shell iron and chromium opacity at stellar interior temperatures
Recently developed free-energy density functional theory (DFT)-based methodology for optical property calculations of warm dense matter has been applied for studying L-shell opacity of iron and chromium at T = 182 eV. We use Mermin–Kohn–Sham density functional theory with a ground-state and a fully temperature-dependent generalized gradient approximation exchange-correlation (XC) functionals. It is demonstrated that the role of XC at such a high-T is negligible due to the total free-energy of interacting system being dominated by the noninteracting free-energy term in agreement with estimations for the homogeneous electron gas. Furthermore, our DFT predictions are compared to the radiative emissivity and opacity of dense plasmas model, to the real-space Greens function method, and to experimental measurements. Good agreement is found between all three theoretical methods, and in the bound–continuum region for Cr when compared to the experiment, while the discrepancy between direct DFT calculations and the experiment for Fe remains essentially the same as for plasma-physics models.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Efficient separation and coprecipitation for simplified cathode recycling
Hydrometallurgical recycling of spent lithium-ion batteries is among the most promising recovery approaches. The current hydrometallurgical method for battery recycling faces challenges such as complex separation and precipitation processes, and environmental concerns from the use of caustic inorganic acid and hazardous hydrogen peroxide. Our proposed modified method, namely polyol-metallurgy, uses citric acid in ethylene glycol as dual-function green solution to overcome these obstacles. This bifunctional solution leaches valuable metal ions from cathode materials (e.g., LiCoO 2 ), and then acts as chelating agents to selectively precipitate cobalt through an esterification reaction, without the need for additional precipitation agents. The leaching efficiency of cobalt and lithium reaches 99.55% and 97.65%, respectively, and more than 96% of cobalt could be directly self-precipitated and recovered. Further, the unique characteristics of the dual-function solution also avoids impurities from Al foils and PVDF/carbon black films, enabling the simple separation process. Therefore, the process can be completed in one-pot system with efficient leaching, separation and coprecipitation, making it more efficient and feasible for operation than existing alternatives.
25 ENERGY STORAGE↗
Chemical structure–property relationships in nanocelluloses
Abstract Nanocelluloses, nanoscale cellulosic particles or fibrils, are an appealing class of nanomaterials with enormous potential for sustainable development. They have great promise to target some of the 17 Sustainable Development Goals outlined by the United Nations. Nanocelluloses are derived from the most abundant biopolymer in the world, cellulose, which have multiple hydroxyl groups on their surface, enabling a vast range of chemical modifications. The interplay between the chemical structure and physicochemical properties of nanocelluloses is crucial to engineering the next generation of green, functional nanomaterials. In this review paper, we provide a comprehensive account of the structure–property relationships in nanocelluloses, which may establish the basis for the design of precision bio‐based materials. Chemical modifications of nanocelluloses, including hydrolysis, oxidation, esterification, amidation, and polymer grafting, provide an array of chemical and physical properties that open opportunities in a diverse spectrum of applications. We review how chemical modifications affect the physicochemical properties of nanocelluloses, such as thermal stability, hydrophobicity, and dispersibility in nonpolar media. Understanding nanocellulose chemical structure–property relationships is integral to the development of sustainable material platforms based on the most renewable biopolymer on earth.
59 BASIC BIOLOGICAL SCIENCES↗
Notes on the complex Sachdev-Ye-Kitaev model
We describe numerous properties of the Sachdev-Ye-Kitaev model for complex fermions with N $\gg$ 1 flavors and a global U(1) charge. We provide a general definition of the charge in the ( G, Σ) formalism, and compute its universal relation to the infrared asymmetry of the Green function. The same relation is obtained by a renormalization theory. The conserved charge contributes a compact scalar field to the effective action, from which we derive the many-body density of states and extract the charge compressibility. We compute the latter via three distinct numerical methods and obtain consistent results. Finally, we present a two dimensional bulk picture with free Dirac fermions for the zero temperature entropy.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Monte Carlo simulations for understanding the transport properties of metallic nanowires
We report the conductance calculation of metallic nanowires of different metals under the application of a tensile force. The elongation induced by the tensile force was simulated by means of the canonical Monte Carlo method combined with the embedded atom method to compute the energy of the system. We studied nanowires composed of Ag, Au, Pt and Cu. The conductance was calculated by means of the Landauer formula and the non-equilibrium Green function formalism. As previously verified experimentally and using other simulation techniques, the conductance is quantized. In addition to this, we have found that there is a strong correlation between the magnitude of the conductance and the width of the nanowire at its thinnest point. Conductance histograms were analyzed for all the different metals.
36 MATERIALS SCIENCE↗
Two-dimensional shaping of Solov'ev equilibrium with vacuum using external coils
In this work, we demonstrate a method for constructing the Solov'ev equilibrium with any given 2D shape surrounded by a vacuum region using external poloidal field coils, based on the method previously developed by T. Xu and R. Fitzpatrick [Nuclear Fusion, 59, 064002 (2019)]. The computational domain consists of two parts: the plasma region, where the solution is the same as the Solov'ev solution, and the vacuum region, where the magnetic field generated by external coils as well as plasma current is determined using the Green function method through a matching condition near the separatrix. However, the method is not limited to the Solov'ev equilibrium in particular. In conclusion, the accuracy, efficiency, and robustness of such a scheme suggest that this method may be applied to the 2D shaping of a tokamak plasma with a vacuum region using external coils in general.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
The effect of curvature induced broken potential vorticity conservation on drift wave turbulences
Abstract The correlation theory of turbulence suppression (Zhang and Mahajan 1993 Phys. Fluids B5 2000) by velocity shear was constructed by invoking the ansatz of potential vorticity conservation (PVC) that holds for relatively simple (slab) models of drift wave turbulence. It is, therefore, surprising that a detailed modern simulation of the H mode pedestal, using the gyrokinetic code GENE (Hatch et al 2018 Plasma Phys. Control. Fusion 60 084003), found ‘striking agreement’ with the predictions of the analytic model. To understand the reasons for this remarkable agreement, an extended theory that contains finite (magnetic) curvature, and which does not conserve potential vorticity, is developed and ‘solved’ by calculating an inhomogeneous Green function reflecting the fact that the new system has a potential vorticity source. It is, then, demonstrated that the effect of the broken PVC is insignificant for the normal operation parameters in tokamaks; the correction due to curvature is at the order of 2 L 0 / R , where L 0 is the scale length of the local gradient and R is the major radius. The excellent agreement between simulation and slab model 1993 theory is therefore not accidental; the latter can be applied with confidence to the tokamak pedestal.
Physics↗
Imaging the crust and uppermost mantle structure of Portugal (West Iberia) with seismic ambient noise
SUMMARY We present a new high-resolution 3-D shear wave velocity (Vs) model of the crust and uppermost mantle beneath Portugal, inferred from ambient seismic noise tomography. We use broad-band seismic data from a dense temporary deployment covering the entire Portuguese mainland between 2010 and 2012 in the scope of the WILAS project. Vertical component data are processed using phase correlation and phase weighted stack to obtain empirical Green functions (EGFs) for 2016 station pairs. Further, we use a random sampling and subset stacking strategy to measure robust Rayleigh-wave group velocities in the period range 7–30 s and associated uncertainties. The tomographic inversion is performed in two steps: First, we determine group-velocity lateral variations for each period. Next, we invert them at each grid point using a new trans-dimensional inversion scheme to obtain the 3-D shear wave velocity model. The final 3-D model extends from the upper crust (5 km) down to the uppermost mantle (60 km) and has a lateral resolution of ∼50 km. In the upper and middle crusts, the Vs anomaly pattern matches the tectonic units of the Variscan Massif and Alpine basins. The transition between the Lusitanian Basin and the Ossa Morena Zone is marked by a contrast between moderate- and high-velocity anomalies, in addition to two arched earthquake lineations. Some faults, namely, the Manteigas–Vilariça–Bragança fault and the Porto–Tomar–Ferreira do Alentejo fault, have a clear signature from the upper crust down to the uppermost mantle (60 km). Our 3-D shear wave velocity model offers new insights into the continuation of the main tectonic units at depth and contributes to better understanding the seismicity of Portugal.
Silveira, Graça (ORCID:0000000221102554)↗
Pairing instability on a Luttinger surface: A non-Fermi liquid to superconductor transition and its Sachdev-Ye-Kitaev dual
Superconductivity results from an instability of the Fermi surface—contour of poles of the single-particle propagator—to an infinitesimally small attraction between electrons. Here, we instead discuss the analogous problem on a model Luttinger surface, or contour of zeros of the Green function. At zero temperature ($β → \infty$) and a critical interaction strength ($u_{c∞}$) characterized by the residue of self-energy pole, we find that the pair susceptibility diverges leading to a superconducting instability. We evaluate the pair fluctuation partition function and find that the spectral density in the normal state has an interaction-driven, power-law $\frac{1}{\sqrt{ω}}$ type, van-Hove singularity (vHS), indicating non-Fermi liquid (NFL) physics. Crucially, in the strong coupling limit ($βu \gg 1$), the leading-order fluctuation free-energy terms in the normal state of this NFL-SC transition resembles the equivalent [$\textit{O}(1)$] terms of the Sachdev-Ye-Kitaev (SYK) model. This free-energy contribution takes a simple form $–βF = βu_{c∞} – γ \text{ln}(βu_{c∞}$), where $\gamma$ is a constant equal to $\frac{1}{2}$. Weak impurity scattering ($τ \gg β^{–1}$) leaves the low-energy spectral density unaffected but leads to an interaction-driven enhancement of superconductivity. Our results in this work shed light on the role played by order-parameter fluctuations in providing the key missing link between Mott physics and strongly coupled toy models exhibiting gravity duals.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗