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 289 records · Page 16

Low-temperature thermal properties of Fe-doped Co 3 Sn 2 S 2

Recent studies on Co 3 Sn 2 S 2 usually focus on its electrical and magnetic behaviors as it is a Weyl semimetal, half metal, and anisotropic ferromagnet at the same time. Only a few works looked into its low-temperature thermal characteristics. Here, in this study, we report the low-temperature heat capacity, Seebeck coefficient, and thermal conductivity of Fe-doped Co 3 Sn 2 S 2 , as well as the formation energies, band structures, and spin-polarized density of states calculated by using the spin-polarized relativistic Korringa-Kohn-Rostoker – Green’s function method (sprKKR). The heat capacity data reveal that the sample with a heavier Fe-doping level has a lower Debye temperature and a larger density of states at the Fermi level, the latter confirmed by our ab initio computational results. The Fe-doping increases the absolute value of the Seebeck coefficient at temperatures above 120 K and dramatically suppresses the thermal conductivity. Furthermore, we discuss how the external applied magnetic field affects the Seebeck coefficient and thermal conductivity of pristine Co 3 Sn 2 S 2 . Our findings provide a detailed reference and may stimulate future investigations into the low-temperature thermal transport properties of Co 3 Sn 2 S 2 and other compounds in the shandite family.

36 MATERIALS SCIENCE↗

CoSyR: A novel beam dynamics code for the modeling of synchrotron radiation effects

The self-consistent nonlinear dynamics of a relativistic charged particle beam interacting with its complete self-fields is a fundamental problem underpinning many of the accelerator design issues in high brightness beam applications, as well as the development of advanced accelerators. Particularly, synchrotron radiation induced effects in a magnetic dispersive beamline element can lead to collective beam instabilities and emittance growth. A novel beam dynamics code is developed based on a Lagrangian method for the calculation of the particles’ radiation near-fields using wavefront/wavelet meshes via the Green’s function of the Maxwell equations. These fields are then interpolated onto a moving mesh for dynamic update of the beam. This method allows radiation co-propagation and self-consistent interaction with the beam in 2D/3D simulations at greatly reduced numerical errors. Multiple levels of parallelisms are inherent in this method and implemented in our code CoSyR to enable at-scale simulations of nonlinear beam dynamics on modern computing platforms using MPI, multi-threading, and GPUs. Here, the current 2D implementation of CoSyR has been used to evaluate the transverse and longitudinal coherent radiation effects on the beam and to investigate beam optics designs proposed for mitigation of beam brightness degradation in a magnetic bunch compressor. In this paper, the design of CoSyR, as well as the benchmark with other coherent synchrotron radiation models, are described and discussed. Extension of the core algorithms to 3D is possible and planned.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Heat flow enhancement in a nanoscale plasmonic junction induced by Kondo resonances and electron-phonon coupling

Here recently, we showed that plasmon-exciton coupling can increase entropy current through a bridge coupled to plasmonic metal nanoparticles. Here we show that electron-phonon coupling can also be used to control the entropy current in similar systems. Entropy current tends to decrease due to electron-phonon coupling and to exhibit a monotonous decrease upon temperature ramping. However, an anomaly affecting the current where it is enhanced by electron-phonon coupling is indicated at around 42 times the system's Kondo temperature. We therefore report means to control heat flow by tuning the Kondo resonance through the electron-phonon coupling. We analyze the conditions that bring about these trends due to electron-phonon coupling by employing non-equilibrium Green's function formulation addressing the entropy current and the derived heat flow.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Analysis of spin frustration in an Fe III 7 cluster using a combination of computational, experimental, and magnetostructural correlation methods

The synthesis, structure, and magnetic properties are reported for [Fe 7 O 3 (O 2 C t Bu) 9 (mda) 3 (H 2 O) 3 ] ( 1 ), where mdaH 2 is N -methyldiethanolamine. 1 was prepared from the reaction of [Fe 3 O(O 2 C t Bu) 6 (H 2 O) 3 ](NO 3 ) with mdaH 2 in a 1:~3 ratio in MeCN. The core of 1 consists of a central octahedral Fe III ion held within a non-planar Fe 6 loop by three μ 3 -O 2- and three μ 2 -RO - arms from the three mda 2- chelates. Variable-temperature dc and ac magnetic susceptibility studies revealed dominant antiferromagnetic coupling, leading to a ground state spin of S = 5 / 2 . The ground state was confirmed by a fit of magnetization data collected in the 0.1–7.0 T and 1.8–10.0 K ranges. The four Fe 2 pairwise exchange parameters ( J 1 - J 4 ) were estimated by independent methods: theoretical calculations using either broken symmetry energy differences (-46.3, -16.2, -3.9, and - 28.1 cm -1 , respectively) or Green’s function approximation methods (-41.4, -14.8, -13.2, and - 24.7 cm -1 ), and a magnetostructural correlation (MSC) previously developed for high nuclearity Fe III /O complexes (-39.5, -13.8, -6.7, and - 23.5 cm -1 ). Additionally, the J 1 - J 4 obtained from the MSC and theoretical methods were used with the program PHI to both simulate χ M T vs T as well as to serve as reasonable input values to fit the experimental data (-41.0, -11.4, -5.0, and - 27.3 cm -1 ). Analysis of the J ij led to identification of the spin frustration effects operative and the resultant spin vector alignments at each Fe III ion, thus allowing for the rationalization of the experimental ground state.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Switching limits of top-gated carbon nanotube field-effect transistors

The performance limits of carbon nanotube field-effect transistors (CNFETs) based on a recently reported process are studied by computational techniques at temperatures between 300 K and 4 K. Here, the impact of band-to-band tunneling (BTBT) and source-to-drain tunneling (SDT) is examined for devices with varying gate length through the use of simulations based on the non-equilibrium Green’s function (NEGF) formalism, and calibrated to measurements. Additionally, the case of junctionless chemical doping profiles is analyzed in contrast to electrostatic doping recently reported for test structures. The switching limits of CNFETs are further explored for devices based on carbon nanotubes (CNTs) with more favorable electronic structures.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Monte Carlo MP2-F12 for Noncovalent Interactions: The C 60 Dimer

A scalable stochastic algorithm is presented that can evaluate explicitly correlated (F12) second-order many-body perturbation (MP2) energies of weak, noncovalent, intermolecular interactions. It first transforms the formulas of the MP2 and F12 energy differences into a short sum of high-dimensional integrals of Green’s functions in real space and imaginary time. Furthermore, these integrals are then evaluated by the Monte Carlo method augmented by parallel execution, redundant-walker convergence acceleration, direct-sampling autocorrelation elimination, and control-variate error reduction. By sharing electron-pair walkers across the supermolecule and its subsystems spanned by the joint basis set, the statistical uncertainty is reduced by one to 2 orders of magnitude in the MP2 binding energy corrected for the basis-set incompleteness and superposition errors. The method predicts the MP2-F12/aug-cc-pVDZ binding energy of 19.1 ± 4.0 kcal mol –1 for the C 60 dimer at the center distance of 9.748 Å.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Anharmonic phonon dispersion in polyethylene

The second-order Green's function method for anharmonic crystals has been applied to an infinite, periodic chain of polyethylene taking into account up to quartic force constants. The frequency-independent approximation to the Dyson self-energy gives rise to numerous divergent resonances, which are fortuitous. Instead, solving the Dyson equation self-consistently with a frequency-dependent self-energy resists divergences from resonances or zero-frequency acoustic vibrations. Here, the calculated anharmonic phonon dispersion, which nonetheless displays many true resonances, and anharmonic phonon density of states furnish hitherto unknown details that explain smaller features of observed vibrational spectra.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Exciton-Defect Interaction and Optical Properties from a First-Principles T-Matrix Approach

Understanding exciton-defect interactions is critical for optimizing optoelectronic and quantum information applications in many materials. However, ab initio simulations of material properties with defects are often limited to high defect density. Here, we study effects of exciton-defect interactions on optical absorption and photoluminescence spectra in monolayer MoS 2 using a first-principles T-matrix approach. We demonstrate that exciton-defect bound states can be captured by the disorderaveraged Green’s function with the T-matrix approximation and further analyze their optical properties. Our approach yields photoluminescence spectra in good agreement with experiments and provides a new, computationally efficient framework for simulating optical properties of disordered 2D materials from firstprinciples.

T-matrix↗

A Nonequilibrium Thermodynamic Approach for Surface Energy Balance Closure

Abstract The surface energy imbalance, viz. that the turbulent dissipation does not fully account for available energy, has for long been an outstanding challenge in geophysical studies. In this study, we developed a novel approach based on nonequilibrium thermodynamics by representing the atmospheric boundary layer as a heat engine. In addition, an analytically tractable approach was used to estimate the ground heat flux based on Green's function approach, which in turn determines the available energy that drives the atmospheric heat engine. The proposed model was evaluated using heat fluxes measured by eight AmeriFlux eddy covariance towers with atmospheric temperature profiles recorded at adjacent radiosonde sites. The surface energy balance closure can be improved by ~11% over various landscapes, by including the estimated power production from the atmospheric heat engine.

Li, Peiyuan↗

Seismicity and Stress Associated With a Fluid-Driven Fracture: Estimating the Evolving Geometry

A coupled approach, combining the theory of rate- and state-dependent friction and methods from poroelasticity, forms the basis for a quantitative relationship between displacements and fluid leak-off from a growing fracture and changes in the rate of seismic events in the region surrounding the fracture. Poroelastic Green's functions link fracture aperture changes and fluid flow from the fracture to changes in the stress field and pore pressure in the adjacent formation. The theory of rate- and state-dependent friction provides a connection between Coulomb stress changes and variations in the rate of seismic events. Numerical modeling indicates that the Coulomb stress changes can vary significantly between formations with differing properties. The relationship between the seismicity rate changes and the changes in the formation stresses and fluid pressure is nonlinear, but a transformation produces a quantity that is linearly related to the aperture changes and fluid leak-off from the fracture. The methodology provides a means for mapping changes in seismicity into fracture aperture changes and to image an evolving fracture. An application to observed microseismicity associated with a hydrofracture reveals asymmetric fracture propagation within two main zones, with extended propagation in the upper zone. Furthermore, the time-varying volume of the fracture agrees with the injected volume, given by the integration of rate changes at the injection well, providing validation of the estimated aperture changes.

58 GEOSCIENCES↗

Regional Moment Tensor Inversion Using Rotational Observations

There are benefits from the addition of rotational motions to translational displacements for moment tensor (MT) inversions. The rotational radiation pattern is orthogonal to the shear radiation pattern, thus incorporating rotations is equivalent to gaining another observation point on the focal sphere. We demonstrated this by simulating curl and displacement wavefields for a regional distance station. Thus, one 6-C station (3-Component translational + 3-Component rotational) gathers the same information on radiation pattern as two 3-C stations at 90° azimuth from one another along the focal plane axis, which is sometimes difficult to obtain when restricted to surface sensors. We added rotational Green's functions to a regional MT inversion scheme (long-period, time-domain, and linear inversion) by computing spatial gradients from f-κ reflectivity synthetics. For rotational data, we used Array Derived Rotations (ADR) from Piñon Flats Observatory Array and Golay array deployed during IRIS Community Wavefield Demonstration Experiment. The hope is to ultimately use compact and field deployable broadband rotational seismometers instead of ADRs. Rotational motions were predicted from well-constrained deviatoric MT solutions of nine earthquakes recorded by the arrays and other seismic networks. Additionally, we formed three station sparse datasets with (3-C and 6-C) and without rotational ground motions (3-C only) for eight earthquakes testing the benefits of including rotational waveforms in MT inversions when they are equally weighted with translational waveforms. Adding rotational motions improved the double-couple components and reduced compensated-linear-vector-dipole and isotropic components in full-MT solutions using sparse datasets with poor coverage.

58 GEOSCIENCES↗

The dominant contribution of Southern Ocean heat uptake to time - evolving radiative feedback in CESM

Radiative feedbacks are found to vary with time in both historical records and future warming projections. Previous studies proposed two factors that determine the variation of radiative feedbacks: (i) the evolution of tropical sea surface warming patterns and (ii) the tropical-extratropical contrast of ocean heat uptake. Our results bridge the two factors by evaluating the remote impact from the extratropical ocean on tropical temperature patterns, accounting for the changes in radiative feedbacks. Based on the Green’s Function approach that quantifies the non-local contributions of regional ocean heat uptake, we show that the net radiative feedback evolution in CESM can be mostly attributed to the heat uptake variations in the Southern Ocean. The enhanced surface warming associated with the weakened heat uptake decades after quadrupling CO2 is not confined over the Southern Ocean, but extends to tropical Southeastern Pacific, which leads to decreasing tropospheric stability and more positive cloud feedbacks.

Lin, Yuan-Jen↗

Synthetic Evaluation of Infrasonic Multipole Waveform Inversion

Acoustic source inversions estimate the mass flow rate of volcanic explosions or yield of chemical explosions and provide insight into potential source directionality. However, the limitations of applying these methods to complex sources and their ability to resolve a stable solution have not been investigated in detail. We perform synthetic infrasound waveform inversions that use 3-D Green’s functions for a variety of idealized and realistic deployment scenarios using both a flat plane and Yasur volcano, Vanuatu as examples. We investigate the ability of various scenarios to retrieve the input source functions and relative amplitudes for monopole and multipole (monopole and dipole) inversions. Infrasound waveform inversions appear to be a robust method to quantify mass flow rates from simple sources (monopole) using deployments of infrasound sensors placed around a source, but care should be taken when analyzing and interpreting results from more complex acoustic sources (multipole) that have significant directional components. In the examples we consider the solution is stable for monopole inversions with a signal-to-noise ratio greater than five and the dipole component is small. For most scenarios investigated, the vertical dipole component of the multipole explosion source is poorly constrained and can impact the ability to recover the other source term components. Because multipole inversions are ill-posed for many deployments, a low residual does not necessarily mean the proper source vector has been recovered. Synthetic studies are able to help investigate the limitations and place bounds on information that may be missing using monopole and multipole inversions for potentially directional sources.

58 GEOSCIENCES↗

Neutral Mode Dominates the Forced Global and Regional Surface Temperature Response in the Past and Future

Abstract Using a large suite of Green's function perturbation experiments, we construct the emergent dynamical operator for surface temperature and extract the most excitable mode. The leading mode turns out to be the most dominant mode excited by the climate forcing of doubling CO 2 , alone capturing 56% of the total spatial variances of the surface temperature response. The pattern of the leading mode is partly organized by the atmospheric annular modes in both hemispheres, and further modulated by the feedbacks from ocean circulations. Though derived from a single model, the leading mode can capture the spatio‐temporal variabilities of the global mean surface temperature over the past century in both the CMIP6 model simulations and observations. Moreover, the leading mode is most efficiently excited by radiative forcings from the midlatitude bands centered around 45°N and 25°S, where most human activities take place, suggesting the potency of human influence on the climate change.

54 ENVIRONMENTAL SCIENCES↗

Explaining Forcing Efficacy With Pattern Effect and State Dependence

The magnitude of global surface temperature change in response to unit radiative forcing depends on the type and magnitude of forcing agent—a concept known as a “forcing efficacy.” However, the mechanisms behind the forcing efficacy are still unclear. In this study, we perform a set of simulations using CESM1 to calculate the efficacy of 10 different forcing agents defined in terms of fixed-SST effective radiative forcing, and then use a Green's function approach to show that each forcing efficacy can be largely understood in terms of the radiative feedbacks associated with the different surface temperature patterns induced by the forcing agents (a pattern effect). We also quantify how the state dependence of feedbacks on global mean surface temperature anomalies impacts forcing efficacies. The results show that the forcing efficacy can be well reconstructed with a combination of pattern effect and state dependence.

58 GEOSCIENCES↗

A Comprehensive Stress Drop Map From Trench to Depth in the Northern Chilean Subduction Zone

We compute stress drops for earthquakes in Northern Chile recorded between 2007 and 2021. By applying two analysis techniques, (a) the spectral ratio (SR) method and (b) the spectral decomposition (SDC) method, a stress drop map for the subduction zone consisting of 51,510 stress drop values is produced. We build an extended set of empirical Green’s functions (EGF) for the SR method by systematic template matching. Outputs are used to compare with results from the SDC approach, where we apply cell-wise obtained global EGF's to compensate for the structural heterogeneity of the subduction zone. We find a good consistency of results of the two methods. The increased spatial coverage and quantity of stress drop estimates from the SDC method facilitate a consistent stress drop mapping of the different seismotectonic domains. Albeit only small differences of median stress drop, strike-perpendicular depth sections clearly reveal systematic variations, with earthquakes at different seismotectonic locations exhibiting distinct values. In particular, interface seismicity is characterized by the lowest observed median value, whereas upper plate earthquakes show noticeably higher stress drop values. Intermediate depth earthquakes show comparatively high average stress drop and a rather strong depth-dependent increase of median stress drop. Additionally, we observe spatio-temporal variability of stress drops related to the occurrence of the two megathrust earthquakes in the study region. The presented study is the first coherent large scale 3D stress drop mapping of the Northern Chilean subduction zone. It provides an important component for further detailed analysis of the physics of earthquake ruptures.

58 GEOSCIENCES↗

CANVAS: An Adjoint Waveform Tomography Model of California and Nevada

Abstract We present the California‐Nevada Adjoint Simulations (CANVAS) model, an adjoint waveform tomography model of the crust and uppermost mantle of the states of California and Nevada. We used WUS256 (Rodgers et al., 2022, https://doi.org/10.1029/2022jb024549 ) as the starting model and iteratively decreased the minimum period of CANVAS from 30 to 12 s. CANVAS was iterated in two distinct stages: the first stage with source mechanisms from the Global Centroid Moment Tensor (GCMT) catalog and the second stage with inverted moment tensors (MT) using the CANV_WUS model (Doody et al., 2023, https://doi.org/10.1029/2023jb026463 ). We show that updating the MTs with 3D Green's functions improved waveform fits and azimuthal coverage of windowed data used to calculate the gradients. As for the model itself, we improved waveform fits over WUS256, particularly in the dispersed surface waves. CANVAS resolved tectonic features seen in other models and accurately defined the depth to basement of major basins, including the Central Valley and the Ventura Basin. We propose CANVAS as a starting model for crustal tomography models on smaller scales.

58 GEOSCIENCES↗

A Frequency Domain Methodology for Quantitative Evaluation of Diffuse Wavefield With Applications to Seismic Imaging

Abstract Ambient Noise Imaging (ANI) of subsurface structures relies on seismic interferometry of diffuse seismic wavefields. However, the lack of effective methods to quantify and identify highly diffuse waves hampers applications of ANI, particularly in evaluating seismic attenuation and monitoring structural changes with high temporal resolution. Conventional ANI approaches require data normalization, which effectively suppresses the non‐diffuse component with large amplitude but also results in significant loss of amplitude and phase information in the continuous seismic records. In this study, we propose a frequency domain method to quantitatively evaluate the degree of diffuseness of seismic wavefields by analyzing their statistical characteristics of modal amplitudes for stationarity and randomness. Tests on synthetic waveform and field nodal records show that the proposed method can effectively distinguish between diffuse and non‐diffuse waveforms for either single‐ or three‐component data. As an application, we identify a 60‐s‐long diffuse coda of a local M 2.2 earthquake recorded by a dense nodal array on the San Jacinto Fault Zone, and successfully extract high‐quality dispersion curve andQ‐value without performing data normalization. These results are consistent with those obtained by conventional methods that assess the correlation between coherency and the Green's function, and by modeling ballistic waves generated by road traffic. Our proposed method can advance the imaging of subsurface velocity and attenuation structures as well as monitoring temporal changes for scientific studies and engineering applications.

Geochemistry & Geophysics↗