Search NASA⌕ Search

SEARCH · Search NASA

Results for “Greens 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,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Solving quasi-free and quadratic Lindblad master equations for open fermionic and bosonic systems

The dynamics of Markovian open quantum systems are described by Lindblad master equations. For fermionic and bosonic systems that are quasi-free, i.e. with Hamiltonians that are quadratic in the ladder operators and Lindblad operators that are linear in the ladder operators, we derive the equation of motion for the covariance matrix. This determines the evolution of Gaussian initial states and the steady states, which are also Gaussian. Using ladder super-operators (a.k.a. third quantization), here we show how the Liouvillian can be transformed to a many-body Jordan normal form which also reveals the full many-body spectrum. Extending previous work by Prosen and Seligman, we treat fermionic and bosonic systems on equal footing with Majorana operators, shorten and complete some derivations, also address the odd-parity sector for fermions, give a criterion for the existence of bosonic steady states, cover non-diagonalizable Liouvillians also for bosons, and include quadratic systems. In extension of the quasi-free open systems, quadratic open systems comprise additional Hermitian Lindblad operators that are quadratic in the ladder operators. While Gaussian states may then evolve into non-Gaussian states, the Liouvillian can still be transformed to a useful block-triangular form, and the equations of motion for k-point Green’s functions form a closed hierarchy. Based on this formalism, results on criticality and dissipative phase transitions in such models are discussed in a companion paper.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Approximating accelerator impedances with resonator networks

It is common in the accelerator community to use the impedance of accelerator components to describe wake interactions in the frequency domain. However, it is often desirable to understand such wake interactions in the time domain in a general manner for excitations that are not necessarily Gaussian in nature. Here, the conventional method for doing this involves taking the inverse Fourier Transform of the component impedance, obtaining the Green's Function, and then convolving it with the desired excitation distribution. This method can prove numerically cumbersome, for a convolution integral must be evaluated for each individual point in time when the wake function is desired. An alternative to this method would be to compute the wake function analytically, which would sidestep the need for repetitive integration. Only a handful of cases, however, are simple enough for this method to be tenable. One of these cases is the case where the component in question is an RLC resonator, which has a closed-form analytical wake function solution. This means that a component which can be represented in terms of resonators can leverage this solution. As it happens, common network synthesis techniques may be used to map arbitrary impedance profiles to RLC resonator networks in a manner the accelerator community has yet to take advantage of. In this work, we will use Foster Canonical Resonator Networks and partial derivative descent optimization to develop a technique for synthesizing resonator networks that well approximate the impedances of real-world accelerator components. We will link this synthesis to the closed-form resonator wake function solution, giving rise to a powerful workflow that may be used to streamline beam dynamics simulations.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Nonequilibrium spectral moment sum rules of the Holstein–Hubbard model

Abstract We derive a general procedure for evaluating the n th derivative of a time-dependent operator in the Heisenberg representation and employ this approach to calculate the zeroth to third spectral moment sum rules of the retarded electronic Green’s function and self-energy for a system described by the Holstein–Hubbard model allowing for arbitrary spatial and time variation of all parameters (including spatially homogeneous electric fields and parameter quenches). For a translationally invariant (but time-dependent) Hamiltonian, we also provide sum rules in momentum space. The sum rules can be applied to various different phenomena like time-resolved angle-resolved photoemission spectroscopy and benchmarking the accuracy of numerical many-body calculations. This work also corrects some errors found in earlier work on simpler models.

Physics↗
DOE OSTI2019

Modeling the Kondo effect of a magnetic atom adsorbed on graphene

The Kondo effect due to a magnetic atom adsorbed on graphene is investigated using two theoretical approaches, including the non-equilibrium Green's function method within the slave-boson mean-field approximation and a newly developed tight-binding (TB) model with an effectively infinite Hubbard U term. Both methods reveal the presence of a Kondo peak in the local density of state (LDOS) of graphene near the Fermi level. A sharp Kondo peak is only observed in the odd-neighboring C sites of the C atom directly below the magnetic atom placed in top position. The peak intensity is found to decay quickly with respect to the distance between the C site and the magnetic atom. A theoretical model of scanning tunneling microscopy (STM) of graphene is presented as a means to identify the manifestation of the Kondo effect via direct topographic STM measurements. STM simulations show that in the proximity of the magnetic atom, only the C atoms of one sublattice are visible at low bias potential and thus a triangular lattice can be seen, in agreement with recent experiments. Further away from the magnetic atom, the Kondo effect dies down and eventually both sublattices are visible, leading to the recovery of the hexagonal lattice. Furthermore, this work not only provides a new TB method to study the Kondo effect in graphene, but also presents both scanning tunneling spectroscopy and microscopy fingerprints of the Kondo effect to facilitate its experimental verification.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Profiling the Quito basin (Ecuador) using seismic ambient noise

SUMMARY Quito, the capital of Ecuador, with more than 2.5 M inhabitants, is exposed to a high seismic hazard due to its proximity to the Pacific subduction zone and active crustal faults, both capable of generating significant earthquakes. Furthermore, the city is located in an intermontane piggy-back basin prone to seismic wave amplification. To understand the basin’s seismic response and characterize its geological structure, 20 broad and medium frequency band seismic stations were deployed in Quito’s urban area between May 2016 and July 2018 that continuously recorded ambient seismic noise. We first compute horizontal-to-vertical spectral ratios to determine the resonant frequency distribution in the entire basin. Secondly, we cross-correlate seismic stations operating simultaneously to retrieve interstations surface-wave Green’s functions in the frequency range of 0.1–2 Hz. We find that Love waves travelling in the basin’s longitudinal direction (NNE–SSW) show much clearer correlograms than those from Rayleigh waves. We then compute Love wave phase-velocity dispersion curves and invert them in conjunction with the HVSR curves to obtain shear-wave velocity profiles throughout the city. The inversions highlight a clear difference in the basin’s structure between its northern and southern parts. In the centre and northern areas, the estimated basin depth and mean shear-wave velocity are about 200 m and 1800 ms−1, respectively, showing resonance frequency values between 0.6 and 0.7 Hz. On the contrary, the basement’s depth and shear-wave velocity in the southern part are about 900 m and 2500 ms−1, having a low resonance frequency value of around 0.3 Hz. This difference in structure between the centre-north and the south of the basin explains the spatial distribution of low-frequency seismic amplifications observed during the Mw 7.8 Pedernales earthquake in April 2016 in Quito.

Pacheco, D. (ORCID:0000000282355096)↗
DOE OSTI2022

A mixed, unified forward/inverse framework for earthquake problems: fault implementation and coseismic slip estimate

SUMMARY We introduce a new finite-element (FE) based computational framework to solve forward and inverse elastic deformation problems for earthquake faulting via the adjoint method. Based on two advanced computational libraries, FEniCS and hIPPYlib for the forward and inverse problems, respectively, this framework is flexible, transparent and easily extensible. We represent a fault discontinuity through a mixed FE elasticity formulation, which approximates the stress with higher order accuracy and exposes the prescribed slip explicitly in the variational form without using conventional split node and decomposition discrete approaches. This also allows the first order optimality condition, that is the vanishing of the gradient, to be expressed in continuous form, which leads to consistent discretizations of all field variables, including the slip. We show comparisons with the standard, pure displacement formulation and a model containing an in-plane mode II crack, whose slip is prescribed via the split node technique. We demonstrate the potential of this new computational framework by performing a linear coseismic slip inversion through adjoint-based optimization methods, without requiring computation of elastic Green’s functions. Specifically, we consider a penalized least squares formulation, which in a Bayesian setting—under the assumption of Gaussian noise and prior—reflects the negative log of the posterior distribution. The comparison of the inversion results with a standard, linear inverse theory approach based on Okada’s solutions shows analogous results. Preliminary uncertainties are estimated via eigenvalue analysis of the Hessian of the penalized least squares objective function. Our implementation is fully open-source and Jupyter notebooks to reproduce our results are provided. The extension to a fully Bayesian framework for detailed uncertainty quantification and non-linear inversions, including for heterogeneous media earthquake problems, will be analysed in a forthcoming paper.

58 GEOSCIENCES↗
DOE OSTI2022

Combining different 3-D global and regional seismic wave propagation solvers towards box tomography in the deep Earth

SUMMARY In previous publications, we presented a general framework, which we called ‘box tomography’, that allows the coupling of any two different numerical seismic wave propagation solvers, respectively outside and inside a target region, or ‘box’. The goal of such hybrid wavefield computations is to reduce the cost of computations in the context of full-waveform inversion for structure within the target region, when sources and/or receivers are located at large distances from the box. Previously, we had demonstrated this approach with sources and receivers outside the target region in a 2-D acoustic spherical earth model, and demonstrated and applied this methodology in the 3-D spherical elastic Earth in a continental scale inversion in which all stations were inside the target region. Here we extend the implementation of the approach to the case of a 3-D global elastic earth model in the case where both sources and stations are outside the box. We couple a global 3-D solver, SPECFEM3D_GLOBE, for the computation of the wavefield and Green’s functions in a reference 3-D model, with a regional 3-D solver, RegSEM, for the computation of the wavefield within the box, by means of time-reversal mirrors. We briefly review key theoretical aspects, showing in particular how only the displacement is needed to be stored at the boundary of the box. We provide details of the practical implementation, including the geometrical design of the mirrors, how we deal with different sizes of meshes in the two solvers, and how we address memory-saving through the use of B-spline compression of the recorded wavefield on the mirror. The proposed approach is numerically efficient but also versatile, since adapting it to other solvers is straightforward and does not require any changes in the solver codes themselves, as long as the displacement can be recovered at any point in time and space. We present benchmarks of the hybrid computations against direct computations of the wavefield between a source and an array of stations in a realistic geometry centred in the Yellowstone region, with and without a hypothetical plume within the ‘box’, and with a 1-D or a 3-D background model, down to a period of 20 s. The ultimate goal of this development is for applications in the context of imaging of remote target regions in the deep mantle, such as, for example, Ultra Low Velocity Zones.

Geochemistry & Geophysics↗
DOE OSTI2025

MTUQ: a framework for estimating moment tensors, point forces, and their uncertainties

SUMMARY We introduce MTUQ, an open-source Python package for seismic source estimation and uncertainty quantification, emphasizing flexibility and operational scalability. MTUQ provides MPI-parallelized grid search and global optimization capabilities, compatibility with 1-D and 3-D Green’s function database formats, customizable data processing, C-accelerated waveform and first-motion polarity misfit functions, and utilities for plotting seismic waveforms and visualizing misfit and likelihood surfaces. Applicability to a range of full- and constrained-moment tensor, point force, and centroid inversion problems is possible via a documented application programming interface, accompanied by example scripts and integration tests. We demonstrate the software using three different types of seismic events: (1) a 2009 intraslab earthquake near Anchorage, Alaska; (2) an episode of the 2021 Barry Arm landslide in Alaska; and (3) the 2017 Democratic People’s Republic of Korea underground nuclear test. With these events, we illustrate the well-known complementary character of body waves, surface waves, and polarities for constraining source parameters. We also convey the distinct misfit patterns that arise from each individual data type, the importance of uncertainty quantification for detecting multimodal or otherwise poorly constrained solutions, and the software’s flexible, modular design.

58 GEOSCIENCES↗
DOE OSTI2020

Brilliant angle-independent structural colours preserved in weevil scales from the Swiss Pleistocene

Extant weevils exhibit a remarkable colour palette that ranges from muted monochromatic tones to rainbow-like iridescence, with the most vibrant colours produced by three-dimensional photonic nanostructures housed within cuticular scales. Although the optical properties of these nanostructures are well understood, their evolutionary history is not fully resolved, in part due to a poor knowledge of their fossil record. Here, we report three-dimensional photonic nanostructures preserved in brightly coloured scales of two weevils, belonging to the genus Phyllobius or Polydrusus , from the Pleistocene (16–10 ka) of Switzerland. The scales display vibrant blue, green and yellow hues that resemble those of extant Phyllobius/Polydrusus . Scanning electron microscopy and small-angle X-ray scattering analyses reveal that the subfossil scales possess a single-diamond photonic crystal nanostructure. In extant Phyllobius/Polydrusus , the near-angle-independent blue and green hues function primarily in crypsis. The preservation of far-field, angle-independent structural colours in the Swiss subfossil weevils and their likely function in substrate matching confirm the importance of investigating fossil and subfossil photonic nanostructures to understand the evolutionary origins and diversification of colours and associated behaviours (e.g. crypsis) in insects.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Cooperative effects in thin dielectric layers: Long-range Dicke superradiance

The realization and control of collective quantum effects so far have predominantly focused on cold atomic ensembles. Quantum photonic platforms, with their engineered Green's functions and integration capability of advanced solid-state quantum emitters, provide opportunities to explore regimes of light-matter interaction beyond the scope of atomic systems. In this work, we demonstrate that embedding quantum emitters within a thin dielectric layer fundamentally alters their collective radiative behavior. The optical modes in the dielectric layer mediate long-range dipole-dipole interactions between emitters, enabling both total and directional superradiance between emitters separated by several wavelengths. Crucially, this mechanism supports Dicke superradiance even in parameter regimes where standard settings fail to support an interaction, unveiling a dimensionality-driven enhancement of cooperative effects. By bridging many-body quantum optics and photonic engineering, our work reveals a distinct interplay between surrounding dimensionality and collective quantum dynamics. Experimental realization of these predictions, readily achievable in solid-state quantum optics platforms, paves the way for scalable, directional quantum light sources and frontiers in many-body quantum optics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Non-Hermitian quantum mechanics approach for extracting and emulating continuum physics based on bound-state-like calculations: Detailed description

Here, this work applies a reduced basis method to study the continuum physics of a finite quantum system—either few or many-body. Specifically, I develop reduced-order models, or emulators, for the underlying inhomogeneous Schrödinger equation and train the emulators against the equation's bound-state-like solutions at complex energies. The emulators rapidly and accurately interpolate and extrapolate the matrix elements of the Hamiltonian resolvent operator (Green's function) across a parameter space that includes both complex energy and other real-valued physical inputs in the Schrödinger equation. The spectra, discretized and compressed as the result of emulation, and the associated resolvent matrix elements (or amplitudes), have the defining characteristics of non-Hermitian quantum mechanics calculations, featuring complex eigenenergies with negative imaginary parts and branch cuts moved below the real axis in the complex energy plane. Therefore, one now has a method that extracts continuum physics from bound-state-like calculations and emulates those extractions in the input parameter space. Building on a prior Letter [Zhang, Phys. Rev. Lett. 135, 242501 (2025)], this article provides the full theoretical details, a comprehensive analysis of the method's performance, and a brief discussion of how it can be coupled with existing continuum approaches to perform emulations in their input parameter spaces.

ab initio calculations↗
DOE OSTI2025

Exact statistics of helical wormlike chains with twist-bend coupling

Here, we present a solution for the Green's function for the general case of a helical wormlike chain with twist-bend coupling and demonstrate the applicability of our solution for evaluating general structural and mechanical chain properties. We find that twist-bend coupling renormalizes the persistence length and the force-extension curves relative to wormlike chains. Analysis of intrinsically twisted polymers shows that incorporation of twist-bend coupling results in the oscillatory behavior in principal tangent correlations that are observed in some studies of synthetic polymers. The exact nature of our solution provides a framework to evaluate the role of twist-bend coupling on polymer properties, such as DNA bound to histones and motivates the reinterpretation of existing biopolymer experimental data.

Ghosh, Ashesh [University of California, Berkeley,↗
DOE OSTI2025

Anisotropic light-tailored RKKY interaction in two-dimensional 𝑑-wave altermagnets

Altermagnets are known in spintronics for their intrinsic spin-splitting and unconventional magnetic responses, particularly to magnetic impurities. However, effectively controlling the magnetic exchange interactions in altermagnets is challenging for practical applications. Here, in this work, we propose using circularly polarized light to tune the Ruderman-Kittel-Kasuya-Yosida (RKKY) interaction in two-dimensional 𝑑-wave altermagnets. Using the real-space retarded Green's functions approach, our results show that while the Heisenberg and Ising exchanges dominate, a notable Dzyaloshinskii–Moriya (DM) interaction also plays a key role. Furthermore, the inherent strength of altermagnetism imprints chirp-like signatures into the magnetic responses, which can be dynamically tuned via light. We mainly demonstrate that gate-induced Rashba spin-orbit coupling is essential in response to light—light selectively and anisotropically adjusts the DM interaction without affecting the other exchanges. Our findings further indicate that rotating the altermagnet by 45° relative to the light's polarization direction generates a Dirac-like dispersion and different DM interactions. We finally extract critical thresholds where light reverses DM interactions along one axis or balances both in-plane components. The anisotropic light-driven control of RKKY interactions in altermagnets not only highlights their unique properties but also opens new avenues for engineering tailored magnetic characteristics in spintronic applications.

altermagnetism↗
DOE OSTI2026

Transformer-based operator learning framework for self-energy in strongly correlated systems

We introduce Σ-Attention, a transformer-based operator-learning framework for approximating the self-energy operator of strongly correlated electronic systems. By creating a batched dataset that combines results from three complementary approaches, i.e., many-body perturbation theory, strong-coupling expansion, and exact diagonalization, each effective in specific parameter regimes, Σ-Attention is applied to learn an accurate approximation for the self-energy operator that is valid across a wide range of parameter regimes. This hybrid strategy leverages the strengths of existing methods while relying on the transformer's ability to generalize beyond individual limitations. More importantly, the scalability of the transformer architecture allows the learned self-energy to be extended to systems with larger sizes, leading to much improved computational scaling. Using the one-dimensional Hubbard model, we demonstrate that Σ-Attention can accurately predict the Matsubara Green's function of large systems with a wide range of coupling strength. Our framework offers a promising and scalable pathway for studying strongly correlated systems with many possible generalizations.

Zhu, Yuanran↗
DOE OSTI2025

Analytical expression of a finite, long, conical canted-cosine-theta coil for particle collider interaction regions

Magnets in the accelerator interaction region (IR) present significant challenges because of high field requirements and limited available space. Conical-shaped magnets offer advantages in these environments by allowing closer placement to the interaction point while maintaining clearance from synchrotron radiation. Interestingly, numerical studies have shown that conical canted-cosine-theta (CCT) designs produce a constant field distribution along the axial direction in the IR quadrupoles for the Electron-Ion Collider (EIC) at Brookhaven National Laboratory. However, the field harmonics generated by conical CCT windings are not yet fully understood. This paper presents an analytical approach to describe the magnetic field produced by a conical surface current and proposes a method for designing conical CCT magnets for accelerator applications. First, we begin with a surface current sheet having a general cosine-theta distribution in spherical coordinates and solve the vector potential using the Green’s function. The magnetic fields generated by the conical current sheet are expressed using associated Legendre polynomials. These results are then related to circular field harmonics and integral field harmonics for designing a coil that produces a pure multipole field. Next, a single layer of the conical CCT winding path is produced based on the cosine-theta current distribution. Finally, the magnetic field quality of dipole and quadrupole conical CCT coils with multiple layers is verified using the Biot-Savart law.

Yang, Ye↗
DOE OSTI2019

Quantum two-photon emission in a photonic cavity

In this paper, we derive an expression for the two-photon spontaneous emission (TPSE) rate of an excited quantum emitter in the presence of arbitrary bodies in its vicinities. After investigating the influence of a perfectly conducting plate on the TPSE spectral distribution (Purcell effect), we demonstrate the equivalence of our expression with the more usual formula written in terms of the corresponding dyadic Green's function. We establish a convenient relation between the TPSE spectral distribution and the corresponding Purcell factors of the system. Next, we consider an emitter close to a dielectric medium and show that, in the near-field regime, the TPSE spectral distribution is substantially enhanced and changes abruptly at the resonance frequencies. Finally, motivated by the suppression that may occur in the one-photon spontaneous emission of an excited atom between two parallel conducting plates, we discuss the TPSE for this same situation and show that complete suppression can never occur for s → s transitions.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Quantum membrane phases in synthetic lattices of cold molecules or Rydberg atoms

Here, we calculate properties of dipolar interacting ultracold molecules or Rydberg atoms in a semisynthetic three-dimensional configuration—one synthetic dimension plus a two-dimensional real-space optical lattice or periodic microtrap array—using the stochastic Green's function quantum Monte Carlo method. Through a calculation of thermodynamic quantities and appropriate correlation functions, along with their finite-size scalings, we show that there is a second-order transition to a low-temperature phase in which two-dimensional “sheets” form in the synthetic dimension of internal rotational or electronic states of the molecules or Rydberg atoms, respectively. Simulations for different values of the interaction V, which acts between atoms or molecules that are adjacent both in real and synthetic space, allow us to compute a phase diagram. We find a finite-temperature transition at sufficiently large V as well as a quantum phase transition—a critical value V c below which the transition temperature vanishes.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Three-dimensional analysis of microbunched electron cooling

Microbunched electron cooling (MBEC, Ratner, ), a special category of coherent electron cooling, has been under consideration as a candidate method for the strong hadron cooling component of the Electron Ion Collider (EIC) at BNL. Most of the theoretical treatments of this scheme so far have been based on simplified models for the key space charge effect in which electrons and/or hadrons are approximated by charged disks. Although such disk-based models provide a flexible tool that captures a good part of the underlying physics, they lack the fidelity and robustness of a rigorous analysis that explicitly treats the hadron and cooler electron beams as collections of point charges. In this paper, we present such a rigorous, three-dimensional (3D) model, both from a theory-based and a simulation perspective. Specifically, we calculate the generalized wakefield (or Green’s function) of the cooling system—a crucial quantity that determines the overall performance—paying special attention to the dependence of the wake on the hadron transverse positions. Additionally, our treatment takes into account previously neglected 3D effects that can cause damping of the plasma oscillations in the amplification section of the cooler, such as focusing and the angular spread of the electron beam. Published by the American Physical Society 2025

43 PARTICLE ACCELERATORS↗