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 127 records · Page 7

Superconductivity from Luttinger surfaces: Emergent Sachdev-Ye-Kitaev physics with infinite-body interactions

The pairing of two electrons on a Fermi surface due to an infinitesimal attraction between them always results in a superconducting instability at zero temperature (T = 0). The equivalent question of pairing instability on a Luttinger surface (LS)—190a contour of zeros of the propagator—instead leads to a quantum critical point (QCP) that separates a non-Fermi liquid and superconductor. A surprising and little understood aspect of pair fluctuations at this QCP is that their thermodynamics maps to that of the Sachdev-Ye-Kitaev (SYK) model in the strong coupling limit. Here, we offer a simple justification for this mapping by demonstrating that (i) LS models share the reparametrization symmetry of the q → ∞ SYK model with q-body interactions close to the LS, and (ii) the enforcement of gauge invariance results in a $\frac{1}{√τ}$ (τ ~ T -1 ) behavior of the fluctuation propagator near the QCP, as is the feature of SYK conformal Green functions, but leaves the overall form of the free energy map robust.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Time Evolution of Correlation Functions in Quantum Many-Body Systems

We give rigorous analytical results on the temporal behavior of two-point correlation functions—also known as dynamical response functions or Green’s functions—in closed many-body quantum systems. We show that in a large class of translation-invariant models the correlation functions factorize at late times $\langle$A ( t ) B$\rangle$ β → $\langle$ A$\rangle$ β $\langle$ B$\rangle$ β , thus proving that dissipation emerges out of the unitary dynamics of the system. We also show that for systems with a generic spectrum the fluctuations around this late-time value are bounded by the purity of the thermal ensemble, which generally decays exponentially with system size. For autocorrelation functions we provide an upper bound on the timescale at which they reach the factorized late time value. Remarkably, this bound is only a function of local expectation values and does not increase with system size. We give numerical examples that show that this bound is a good estimate in nonintegrable models, and argue that the timescale that appears can be understood in terms of an emergent fluctuation-dissipation theorem. Our study extends to further classes of two point functions such as the symmetrized ones and the Kubo function that appears in linear response theory, for which we give analogous results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Quantum Spin-Valley Hall Kink States: From Concept to Materials Design

We propose a general and tunable platform to realize high-density arrays of quantum spin-valley Hall kink (QSVHK) states with spin-valley-momentum locking based on a two-dimensional hexagonal topological insulator. Through the analysis of Berry curvature and topological charge, the QSVHK states are found to be topologically protected by the valley-inversion and time-reversal symmetries. Remarkably, the conductance of QSVHK states remains quantized against both nonmagnetic short- and long-range and magnetic long-range disorder, verified by the Green-function calculations. Based on first-principles results and our fabricated samples, we show that QSVHK states, protected with a gap up to 287 meV, can be realized in bismuthene by alloy engineering, surface functionalization, or electric field, supporting nonvolatile applications of spin-valley filters, valves, and waveguides even at room temperature.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

tert-butoxides as precursors for atomic layer deposition of alkali metal containing thin films

Alkali metal containing materials have become increasingly attractive in a world hunting for sustainable energy materials and green functional devices. Lithium- and sodium battery technology, lead-free piezo- and ferroelectric devices, and record-breaking alkali doped tandem perovskite solar cells are among the applications where alkali metal-containing thin films get increasing attention. Atomic layer deposition (ALD) is one of the enabling thin film deposition techniques that offer chemical and geometrical versatility to realize the implementation of such thin films on an applicable scale. The drawback has until recently been a lack of available precursor chemistry that offers self-limiting growth that is fundamental to ALD. The alkali metal tert-butoxides have been shown to exhibit the necessary properties to facilitate saturating growth for Li-, Na-, K-, and Rb-containing compounds. However, the behavior of the tert-butoxides in ALD-growth has been considered difficult to unravel, with processes exhibiting limited control and low reproducibility. Very little has been reported on trends in reaction mechanisms as the mass of the alkali metal increases. Herein, we summarize the existing literature on the use of alkali metal tert-butoxides as precursors in ALD. We consider differences in the structure and behavior of the tert-butoxides as the alkali metal cation becomes heavier. In addition, we present precursor synthesis routes and key information on precursor structure, stability, and mechanistic behavior. Finally, we provide the first ever report of Cs-containing films by ALD to complement previous work on its lighter counterparts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Parallel Simulation of Beam Dynamics in Particle Accelerators [Slides]

LW3D was used to explore the expected CSR from an undulator-based diagnostic using a first-principles approach. It was also used for comparison with an analytic Green-function approach. LW3D was extended to model transport systems involving dipoles and quadrupoles. It was used to explore the radiation produced in a 4-dipole bunch compressor. Radiation ahead of the electron bunch at the beamline exit was shown to be associated with wavefronts produced in upstream dipoles overtaking the bunch.

43 PARTICLE ACCELERATORS↗

Variational Monte Carlo Calculations of A ≤ 4 Nuclei with an Artificial Neural-Network Correlator Ansatz

Here, the complexity of many-body quantum wave functions is a central aspect of several fields of physics and chemistry where nonperturbative interactions are prominent. Artificial neural networks (ANNs) have proven to be a flexible tool to approximate quantum many-body states in condensed matter and chemistry problems. In this work we introduce a neural-network quantum state ansatz to model the ground-state wave function of light nuclei, and approximately solve the nuclear many-body Schrodinger equation. Using efficient stochastic sampling and optimization schemes, our approach extends pioneering applications of ANNs in the field, which present exponentially scaling algorithmic complexity. We compute the binding energies and point-nucleon densities of A ≤ 4 nuclei as emerging from a leading-order pionless effective field theory Hamiltonian. We successfully benchmark the ANN wave function against more conventional parametrizations based on two- and three-body Jastrow functions, and virtually exact Green's function Monte Carlo results.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Nonparametric, data-based kernel interpolation for particle-tracking simulations and kernel density estimation

Traditional interpolation techniques for particle tracking include binning and convolutional formulas that use pre-determined (i.e., closed-form, parameteric) kernels. In many instances, the particles are introduced as point sources in time and space, so the cloud of particles (either in space or time) is a discrete representation of the Green’s function of an underlying PDE. As such, each particle is a sample from the Green’s function; therefore, each particle should be distributed according to the Green’s function. In short, the kernel of a convolutional interpolation of the particle sample “cloud” should be a replica of the cloud itself. This idea gives rise to an iterative method by which the form of the kernel may be discerned in the process of interpolating the Green’s function. When the Green’s function is a density, this method is broadly applicable to interpolating a kernel density estimate based on random data drawn from a single distribution. We formulate and construct the algorithm and demonstrate its ability to perform kernel density estimation of skewed and/or heavy-tailed data including breakthrough curves.

42 ENGINEERING↗

Quantifying the effect of 3D models on moment tensor results using synthetic data

Moment tensors provide vital information on seismic source properties for seismic events. Moment tensors require seismic wavespeed models to compute Green’s functions, which measure the impulse response between a given source and receiver. Traditionally, researchers have used one-dimensional velocity models to calculate Green’s functions since 1D Green’s functions are computationally cheap to compute. Local 1D velocity models can also accurately model waveforms at short distances (< 500 km). However, 1D velocity models do not account for lateral heterogeneity, which can cause significant misfit in tectonically complex regions such as the Middle East (Covellone and Savage, 2012). Green’s functions calculated using 3D seismic wavespeed models have been shown to perform better in tectonically complex regions (Covellone and Savage, 2012; Kintner and Modrak, 2022), so we are interested in quantifying the effect of considering 3D structure on moment tensor inversion results. The Middle East is an ideal study area for a synthetic moment tensor test for two reasons. Firstly, the Middle East is a tectonically complex region that has been heavily studied. Secondly, the Middle East has significant seismic activity throughout the region, but imperfect station coverage due to limited open data through large swaths of the domain. The tectonic complexity and uneven station coverage will test real-world performance even in a synthetic experiment.

58 GEOSCIENCES↗

Lepton–Nucleus Interactions within Microscopic Approaches

This review paper emphasizes the significance of microscopic calculations with quantified theoretical error estimates in studying lepton–nucleus interactions and their implications for electron scattering and accelerator neutrino oscillation measurements. We investigate two approaches: Green’s Function Monte Carlo and the extended factorization scheme, utilizing realistic nuclear target spectral functions. In our study, we include relativistic effects in Green’s Function Monte Carlo and validate the inclusive electron scattering cross section on carbon using available data. We compare the flux-folded cross sections for neutrino-carbon scattering with T2K and MINER$v$ A experiments, noting the substantial impact of relativistic effects in reducing the theoretical curve strength when compared to MINER$v$ A data. Additionally, we demonstrate that quantum Monte Carlo-based spectral functions accurately reproduce the quasi-elastic region in electron scattering data and T2K flux-folded cross sections. By comparing results from Green’s Function Monte Carlo and the spectral function approach, which share a similar initial target state description, we quantify errors associated with approximations in the factorization scheme and the relativistic treatment of kinematics in Green’s Function Monte Carlo.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Atomistic modeling of energy band alignment in CdTe(1 0 0) and CdTe(1 1 1) surfaces

An atomic-scale perspective of energy band alignment in CdTe surfaces has not been spatially studied despite the major role surfaces play in forming interfaces within CdTe-based thin film photovoltaic devices. Atomistic modeling based on density functional theory coupled with surface Green’s function is used for calculating energy band alignment of CdTe surfaces. The CdTe(1 0 0) ((1×1) and c(2×2) reconstruction) and CdTe(1 1 1) ((1×1) and (2×2) reconstruction) facets without and with surface relaxation provide insightful band bending characteristics that influence charge carrier transport. Results show that unrelaxed (1×1) CdTe(1 0 0) and CdTe(1 1 1) surfaces bend the valence band downward with surface polarity dictating the surface potential magnitude. The reconstructed CdTe(1 0 0) c(2×2) and CdTe(1 1 1) (2×2) surfaces result in favorable surface electronic features in relation to their unreconstructed variants. In addition, the structurally relaxed CdTe(1 1 1) surfaces develop an internal energy cusp potential that may enhance hole charge transport toward the back of CdTe solar cell devices. In conclusion, energy band alignments calculated within the study lead to a detailed understanding of how CdTe surfaces may affect CdTe-based thin film photovoltaic applications.

36 MATERIALS SCIENCE↗

Pre-equilibrium evolution of conserved charges with initial conditions in the ICCING Monte Carlo event generator

Heavy-ion collisions can be well described through relativistic viscous hydrodynamics, but questions still remain when hydrodynamics is applicable because the initial state may begin very far from equilibrium. Thus, a pre-equilibrium evolution phase is used to bridge the gap between the initial state and hydrodynamics. KøMPøST is one such pre-equilibrium model that propagates the energy-momentum tensor by decomposing it into the background and fluctuations around that background, whose evolution is captured by Green's functions. We extend this formalism to include conserved charges and calculate the corresponding nonequilibrium Green's functions in the relaxation-time approximation. The ICCING algorithm initializes conserved charges in the initial state by sampling $g$ → $q$$\overline{q}$ splitting probabilities and is, thus, perfectly positioned to implement Green's functions for charge propagation. As a result, we show that this method alters the initial-state charge geometries and is applicable in central to mid-central collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Method of Green’s Potentials for Elliptic PDEs in Domains with Random Apertures

Problems with topological uncertainties appear in many fields ranging from nano-device engineering to the design of bridges. In many of such problems, a part of the domains boundaries is subjected to random perturbations making inefficient conventional schemes that rely on discretization of the whole domain. Here, we study elliptic PDEs in domains with boundaries comprised of a deterministic part and random apertures, and apply the method of modified potentials with Green’s kernels defined on the deterministic part of the domain. This approach allows to reduce the dimension of the original differential problem by reformulating it as a boundary integral equation posed on the random apertures only. The multilevel Monte Carlo method is then applied to this modified integral equation and its optimal ϵ -2 asymptotical complexity is shown. Finally, we provide the qualitative analysis of the proposed technique and support it with numerical results.

97 MATHEMATICS AND COMPUTING↗

The Effects of Atmospheric Models on the Estimation of Infrasonic Source Functions at the Source Physics Experiment

We invert infrasound signals for an equivalent seismoacoustic source function using different atmospheric models to produce the necessary Green’s functions. The infrasound signals were produced by a series of underground chemical explosions as part of the Source Physics Experiment (SPE). In a previous study, we inverted the infrasound data using so-called predictive atmospheric models, which were based on historic, regional-scaled, publicly available weather observations interpolated onto a 3D grid. For the work presented here, we invert the same infrasound data, but using atmospheric models based on weather data collected in a time window that includes the approximate time of the explosion experiments, which we term postdictive models. We build two versions of the postdictive models for each SPE event: one that is based solely on the regional scaled observations, and one that is based on both regional scaled observations combined with on-site observations obtained by a weather sonde released at the time of the SPE. We then invert the observed data set three times, once for each atmospheric model type. We find that the estimated seismoacoustic source functions are relatively similar in waveform shape regardless of which atmospheric model that we used to construct the Green’s functions. However, we find that the amplitude of the estimated source functions is systematically dependent on the atmospheric model type: using the predictive atmospheric models to invert the data generally yields estimated source functions that are larger in amplitude than those estimated using the postdictive models.

54 ENVIRONMENTAL SCIENCES↗

NESSi : The N on- E quilibrium S ystems S imulation package

The nonequilibrium dynamics of correlated many-particle systems is of interest in connection with pump–probe experiments on molecular systems and solids, as well as theoretical investigations of transport properties and relaxation processes. Nonequilibrium Green’s functions are a powerful tool to study interaction effects in quantum many-particle systems out of equilibrium, and to extract physically relevant information for the interpretation of experiments. Here, we present the open-source software package NESSi (The Non-Equilibrium Systems Simulation package) which allows to perform many-body dynamics simulations based on Green’s functions on the L-shaped Kadanoff–Baym contour. NESSi contains the library libcntr which implements tools for basic operations on these nonequilibrium Green’s functions, for constructing Feynman diagrams, and for the solution of integral and integro-differential equations involving contour Green’s functions. The library employs a discretization of the Kadanoff–Baym contour into time points and a high-order implementation of integration routines. The total integrated error scales up to $\mathcal{O}(N^{-7})$, which is important since the numerical effort increases at least cubically with the simulation time. A distributed-memory parallelization over reciprocal space allows large-scale simulations of lattice systems. We provide a collection of example programs ranging from dynamics in simple two-level systems to problems relevant in contemporary condensed matter physics, including Hubbard clusters and Hubbard or Holstein lattice models. The libcntr library is the basis of a follow-up software package for nonequilibrium dynamical mean-field theory calculations based on strong-coupling perturbative impurity solvers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Accurate numerical, integral methods for computing drift-kinetic Trubnikov-Rosenbluth potentials

A novel numerical method is employed to compute the integral form of the axi-symmetric Trubnikov-Rosenbluth potentials. Two methods for quadrature in pitch-angle are described and their convergence properties are studied. Careful attention is given to quadrature over a singular Green's function. Here it is shown that an infinite series representation of the Green's function can be used more efficiently than its closed form involving complete elliptic integrals. Then a collocation method in speed, with its associated quadrature scheme, is laid out and its convergence properties are studied. Using the proposed scheme, accurate low-order moments of the field collision operator are obtained using relatively few velocity space degrees of freedom. The scheme is showcased by solving for the equilibrium, axi-symmetric bootstrap current in tokamaks. A C 0 Gauss-Lobatto-Legendre finite element pitch-angle basis with vertex nodes at the trapped/passing boundary is shown, in the context of the integral methods used, to be much more efficient than the more common Legendre polynomial expansion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Insights into negative differential resistance in MoS 2 Esaki diodes: A first-principles perspective

MoS 2 is a two-dimensional material with a band gap depending on the number of layers and tunable by an external electric field. The experimentally observed intralayer band-to-band tunneling and interlayer band-to-band tunneling in this material present an opportunity for new electronic applications in tunnel field-effect transistors. However, such a widely accepted concept has yet to be been supported by theoretical investigations based on first principles. In this paper, using density functional theory, in conjunction with nonequilibrium Green's function techniques and our electric field gating method, enabled by a large-scale computational approach, we study the relation between band alignment and transmission in planar and side-stack MoS 2 p–i–n junction configurations. Here, we demonstrate the presence of negative differential resistance for both in-plane and interlayer current, a staple characteristic of tunnel diode junctions, and analyze the physical origin of such an effect. Electrostatic potentials, the van der Waals barrier, and a complex band analysis are also examined for a thorough understanding of Esaki diodes.

74 ATOMIC AND MOLECULAR PHYSICS↗

Real-space representation of the quasiparticle self-consistent GW self-energy and its application to defect calculations

The quasiparticle self-consistent (QS) GW (G for Green's function, W for screened Coulomb interaction) approach incorporates the corrections of the quasiparticle energies from their Kohn-Sham density functional theory (DFT) eigenvalues by means of an energy-independent and Hermitian self-energy matrix usually given in the basis set of the DFT eigenstates. By expanding these into an atom-centered basis set (specifically here the linearized muffin-tin orbitals) a real space representation of the self-energy corrections becomes possible. In this work, We show that this representation is relatively short-ranged. This offers opportunities to construct the self-energy of a complex system from parts of the system by a cut-and-paste method. Specifically for a point defect, represented in a large supercell, the self-energy can be constructed from those of the host and a smaller defect-containing cell. The self-energy of the periodic host can be constructed simply from a GW calculation for the primitive cell. We show for the case of the As Ga in GaAs that the defect part can already be well represented by a minimal eight-atom cell and allows us to construct the self-energy for a 64-atom cell in good agreement with direct QSGW calculations for the large cell. Using this approach to an even larger 216-atom cell shows the defect band approaches an isolated defect level. The calculations also allow us to identify a second defect band which appears as a resonance near the conduction band minimum. The results on the extracted defect levels agree well with Green's function calculations for an isolated defect and with experimental data.

36 MATERIALS SCIENCE↗