Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Structural variation”

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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Large change of interlayer vibrational coupling with stacking in Mo 1-x W x Te 2

Stacking variations in quasi-two-dimensional materials can have an important influence on material properties, such as changing the topology of the band structure. Unfortunately, the weakness of van der Waals (vdW) interactions makes it difficult to compute the stacking dependence of properties, and even in a material as simple as graphite the stacking energetics remain unclear. Mo 1 - x W x Te 2 is a material in which three differently stacked phases are conveniently accessible by temperature changes: 1 T ' , T d * , and the reported Weyl semimetal phase T d . Additionally, the transitions proceed via layer sliding, and the corresponding interlayer shear mode (ISM) is relevant not just for the stacking energetics but also for understanding the relationship between Weyl physics and structural changes. However, the interlayer interactions of Mo 1 - x W x Te 2 are not well understood, with wide variation in computed properties. We report inelastic neutron scattering of the ISM in a Mo 0.91 W 0.09 Te 2 crystal. The ISM energies are generally consistent with the linear chain model, as expected given the weak interlayer interaction, though there are some discrepancies from predicted intensities. However, the interlayer force constants K x in the T d * and 1 T ' phases are substantially weaker than that of T d at 75(3) and 83(3)%, respectively. Considering that the relative positioning of atoms in neighboring layers is approximately the same regardless of overall stacking, our results suggest that longer-range influences, such as stacking-induced electronic band-structure changes, may be responsible for the substantial change in the interlayer vibrational coupling and thus the C 55 elastic constant. These findings should elucidate the stacking energetics of Mo 1 - x W x Te 2 and other vdW layered materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Probing individual single atom electrocatalyst sites by advanced analytical scanning transmission electron microscopy

Single atom electrocatalysts (SAEs) are promising next-generation materials for promoting a variety of important reactions, such as the oxygen reduction, nitrogen reduction, and CO 2 reduction reactions. While bulk characterization techniques such as X-ray absorption spectroscopy and Mössbauer spectroscopy have significantly enhanced our understanding of these catalysts, direct probing of individual single metal atom sites at the atomic scale is necessary to understand local variations in the properties of these sites and accelerate design and synthesis of improved SAEs. Aberration-corrected scanning transmission electron microscopy (STEM) has become a powerful tool for providing this type of atomic-scale information about SAE metal sites. These sites are typically unstable under the electron beam, however, which, in combination with conventional acquisition methods and detectors, has limited the type and quantity of information obtainable by spectroscopic STEM techniques. Here, we map multiple individual SAE metal sites in a nitrogen-doped carbon containing atomically dispersed Fe and Re (FeReNC) at the atomic scale by direct electron detection electron energy-loss spectroscopy (EELS). Direct electron detection provides an improved signal-to-noise ratio over conventional scintillator-based detectors and enables detection and real space localization of weak signals. In addition, we demonstrate an automated method for identification of metal atom positions, placement of the probe on these sites, and simultaneous EELS and energy dispersive X-ray spectroscopic (EDS) signal acquisition. This simultaneous acquisition of EELS and EDS provides access to the composition and bonding of a wide range of SAE metal sites. In this study, focusing the probe directly on the metal sites also increases the relevant data acquisition rate by more than an order of magnitude over two-dimensional mapping, enabling improved statistical measurements of site properties. The versatility, sensitivity, and speed that these techniques provide enhances our ability to probe the local elemental and chemical environment of a large number of individual SAE metal site structures at the atomic scale, enabling an improved understanding of the variations in the local properties of these electrocatalysts to be gained. As a result, significantly increased information about individual metal sites will be available to future electrochemical studies through these techniques, accelerating the development of advanced SAEs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Distortion of a Relativistic Jet Echoing a Magnetic Flux Eruption

Magnetized accretion onto spinning black holes (BHs) can accumulate a large magnetic flux across the event horizon and launch a pair of relativistic jets via the Blandford–Znajek mechanism. In the magnetically saturated (arrested) state, excess magnetic flux is ejected from the BH in episodic magnetic flux eruptions, which result in a significant yet temporary reduction of jet power. We analyze results of a high-resolution 3D general-relativistic magnetohydrodynamic numerical simulation of geometrically thick magnetically saturated accretion onto a high-spin Kerr BH for a single cycle of magnetic flux eruption and accumulation. We show that following an eruption, a weakened jet develops a strong helical distortion with a distinct structure of magnetic fields—the poloidal field along the jet core is unaffected by the eruption; while toroidal field lines, ejected from the BH during the eruption and later readvected onto it, form poloidal “bypasses” along the inner jet sheath. Such a distortion may appear in sources fed by geometrically thick accretion flows as an asymmetric superluminal knot, strongly interacting with the jet sheath along an oblique working surface. The jet section repowered by magnetic flux reaccumulated on the BH is tilted by a few degrees, implying significant variations in radiation boost towards observers of BL Lac blazars. The intrinsic structure of the jet spine is consistent with axisymmetric semi-analytical models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Electron-Proton Scattering Event Generation using Structured Tokenization

Recent work such as Omnijet-$\alpha$ has demonstrated that effective tokenization combined with transformer-based architectures can produce effective foundation models for jet physics. While tokenization may help models capture generalizable event characteristics, it also introduces discretization errors that may compromise the precision required for downstream physics analyses. As the number and complexity of the particle features grow, these errors are likely to grow proportionally. In this study, we investigate new tokenization strategies to improve the application of generative transformer models to \textsc{Pythia8} simulations of electron-proton scattering at the Electron-Ion Collider. Specifically, we propose a feature-based structured tokenization approach that utilizes multiple tokens per particle, improving expressivity, while reducing the total number of unique tokens needed. We evaluate this method against grid-based binning, K-means clustering, and vector-quantized variational auto-encoders on the event simulations. Our results show that feature-based structured tokenization reduces discretization error, leading to more accurate generative modeling of particle-level events.

Goldenberg, Steven [Thomas Jefferson National Acce↗
DOE OSTI2025

Electron-Proton Scattering Event Generation using Structured Tokenization

Recent work such as Omnijet-$\alpha$ has demonstrated that effective tokenization combined with transformer-based architectures can produce effective foundation models for jet physics. While tokenization may help models capture generalizable event characteristics, it also introduces discretization errors that may compromise the precision required for downstream physics analyses. As the number and complexity of the particle features grow, these errors are likely to grow proportionally. In this study, we investigate new tokenization strategies to improve the application of generative transformer models to \textsc{Pythia8} simulations of electron-proton scattering at the Electron-Ion Collider. Specifically, we propose a feature-based structured tokenization approach that utilizes multiple tokens per particle, improving expressivity, while reducing the total number of unique tokens needed. We evaluate this method against grid-based binning, K-means clustering, and vector-quantized variational auto-encoders on the event simulations. Our results show that feature-based structured tokenization reduces discretization error, leading to more accurate generative modeling of particle-level events.

Goldenberg, Steven [Thomas Jefferson National Acce↗
DOE OSTI2023

Strained Superlattice photocathode development using CBE and MBE

High polarization electron beams for accelerators are generated using strained superlattice GaAs-based photocathodes. A collaboration of researchers at University of California Santa Barbara and Jefferson Lab has been investigating growth of SSL photocathodes using either chemical- or molecular-beam epitaxy (CBE or MBE) to re-establish a source of high polarization photocathode material. While calibrating growth parameters for the now standard GaAs/GaAsP strained superlattice structure, the UCSB personnel encountered the drawbacks of growing structures with phosphorus in a MBE or CBE system, including high vapor pressure byproducts, phosphorus contamination in both the GaAs layer and all subsequent materials grown in that chamber, and highly flammable and toxic chamber residues. For these reasons, the UCSB team began additionally investigating strained superlattices of InAlGaAs/AlGaAs, which has been successfully demonstrated in the literature1, to be a high QE, high polarization photocathode material. For GaAs/GaAsP, varying phosphorus content affects both the band gap and polarization. However for InAlGaAs/AlGaAs, In and Al content can be varied independently and allows optimization of both polarization and bandgap. Several variations on InAlGaAs/AlGaAs superlattice photocathodes, including those with distributed Bragg reflector (DBR) structures to enhance QE, have been grown at UCSB and sent to Jefferson Lab for testing. The initial results are quite promising, with polarization over 80% and QE about 0.3%, and ongoing test results will be presented.

Stutzman, Marcy↗
DOE OSTI2022

Evolution of the Heliotail Lobes over a Solar Cycle as Measured by IBEX

Energetic neutral atom (ENA) measurements by IBEX reveal that the heliotail comprises an energy-dependent multilobe structure. We examine the heliotail evolution over 11 yr of IBEX observations covering a full solar cycle (SC). We find the following: (1) The heliotail structure persists over the entire SC, comprising three ENA-enhanced and two ENA-suppressed lobes. (2) Lobe sizes and locations are generally stable but exhibit variations in ENA fluxes driven by the SC. (3) Lobe centers follow a cyclic behavior over multiple SC phases, indicating direct signatures of slow and fast solar wind (SW) interactions in the inner heliosheath (IHS). (4) The tilted plane passing through the port–starboard lobes' centers oscillates in latitude but maintains its tilt from the ecliptic plane, likely a consequence of the interstellar magnetic field draping around the heliosphere. (5) The transition of the central heliotail from a single lobe at ~1.1 keV to two lobes above ~2 keV is SC-dependent and directly reflects the IHS plasma properties, i.e., when ENA fluxes from fast SW from the polar coronal holes change over time. (6) The central lobe exhibits a substructure that is enhanced and offset from the downwind direction, possibly indicating an asymmetric ENA emission or an asymmetry in the parent plasma distribution. These results reveal the general stability of the heliotail structure over time and distinct variations in individual lobes' properties in relation to the SC phases. Furthermore, results show the effects of multiple SC phases in the tail, reflecting different ENA travel times and source histories.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Synthesis and Characterization of Carbon Microbeads

We report a microfluidic-based droplet generation platform for synthesizing micron-sized porous carbon microspheres. The setup employs carbon materials such as graphite, carbon nanotubes, graphene, fullerenes, and carbon black as starting materials. Custom composition, structure, and function are achieved through combinations of carbon materials, cross-linkers, and additives along with variations in process parameters. Carbon materials can be assembled into spheres with a mean diameter of units to hundreds of μm with relatively tight size distribution (<25% RSD). Pore structure and size (tens to hundreds of angstrom) can be modulated by incorporating porogen/coporogen dilutants during synthesis. The microbeads have excellent mechanical stability with an elastic modulus of hundreds of MPa. They can sustain high dynamic fluid flow pressures of up to 9000 psi. This work lays the foundation for synthesizing novel tailorable and customizable carbon microbeads. It opens avenues for applying these novel materials for composite and additive manufacturing, energy, life science, and biomedical applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Precise Modulation of CO 2 Sorption in Ti 8 Ce 2 –Oxo Clusters: Elucidating Lewis Acidity of the Ce Metal Sites and Structural Flexibility

Investigating the structure-property correlation in porous materials is a fundamental and consistent focus in various scientific domains, especially within sorption research. Metal oxide clusters with capping ligands, characterized by intrinsic cavities formed through specific solid-state packing, demonstrate significant potential as versatile platforms for sorption investigations due to their precisely tunable atomic structures and inherent long-range order. This study presents a series of Ti 8 Ce 2 -oxo clusters with subtle variations in coordinated linkers and explores their sorption behavior. Notably, Ti 8 Ce 2 -BA (BA denotes benzoic acid) manifests a distinctive two-step profile during CO 2 adsorption, accompanied by a hysteresis loop. This observation marks a new instance within the metal oxide cluster field. Of intrigue, the presence of unsaturated Ce(IV) sites was found to be correlated with the stepped sorption property. Moreover, the introduction of an electrophilic fluorine atom, positioned ortho or para to the benzoic acid, facilitated precise control over gate pressure and stepped sorption quantities. Advanced in-situ techniques systematically unraveled the underlying mechanism behind this unique sorption behavior. The findings elucidate that robust Lewis base-acid interactions are established between CO 2 molecules and Ce ions, consequently altering the conformation of coordinated linkers. Conversely, the F atoms primarily contribute to gate pressure variation by influencing the Lewis acidity of the Ce sites. This research advances the understanding in fabricating metal-oxo clusters with structural flexibility and provides profound insights into their host-guest interaction motifs. Furthermore, these insights hold substantial promise across diverse fields and offer valuable guidance for future adsorbent designs grounded in fundamental theories of structure-property relationships.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Multiple Wheat Genomes Reveal Novel Gli-2 Sublocus Location and Variation of Celiac Disease Epitopes in Duplicated α-Gliadin Genes

The seed protein α-gliadin is a major component of wheat flour and causes gluten-related diseases. However, due to the complexity of this multigene family with a genome structure composed of dozens of copies derived from tandem and genome duplications, little was known about the variation between accessions, and thus little effort has been made to explicitly target α-gliadin for bread wheat breeding. Here, we analyzed genomic variation in α-gliadins across 11 recently published chromosome-scale assemblies of hexaploid wheat, with validation using long-read data. We unexpectedly found that the Gli-B2 locus is not a single contiguous locus but is composed of two subloci, suggesting the possibility of recombination between the two during breeding. We confirmed that the number of immunogenic epitopes among 11 accessions varied. The D subgenome of a European spelt line also contained epitopes, in agreement with its hybridization history. Evolutionary analysis identified amino acid sites under diversifying selection, suggesting their functional importance. The analysis opens the way for improved grain quality and safety through wheat breeding.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Conditional diffusion machine-learning framework for mapping valence electron distribution from convergent beam electron diffraction

Quantitative convergent beam electron diffraction (CBED) enables determination of aspherical valence electron distributions through refinement of low-order structure factors, which are highly sensitive to chemical bonding and charge density variations. However, conventional quantitative CBED (QCBED) requires solving a highly nonlinear inverse problem with many coupled parameters, and computationally intensive dynamical diffraction calculations, making it time-consuming and difficult to apply to complex systems. More broadly, reconstructing charge density and orbital electron distribution from diffraction data has long been a central challenge in both x-ray and electron crystallography. Here, in this study, we introduce an artificial-intelligence (AI)-based framework that replaces traditional refinement with a data-driven inverse solver. Using a large synthetic CBED dataset generated by Bloch-wave simulations, we train a conditional diffusion model to directly infer crystal structural parameters and multipole density formalism parameters, and hence valence electron distributions, from CBED patterns alone. By learning from forward simulations across realistic parameter space, the model effectively solves the inverse problem. Compared with direct regression approaches, the diffusion-based framework provides posterior parameter distributions for rigorous uncertainty quantification while preserving quantitative fidelity and reducing analysis time by orders of magnitude. By eliminating the need for external single-crystal x-ray diffraction data and complex nonlinear refinement, this approach enables practical, high-throughput, and in situ quantitative CBED, enabling real-time mapping of valence electron distributions and their correlation with functional responses in quantum and energy materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Quantum algorithm for electronic band structures with local tight-binding orbitals

While the main thrust of quantum computing research in materials science is to accurately measure the classically intractable electron correlation effects due to Coulomb repulsion, designing optimal quantum algorithms for simpler problems with well-understood solutions is a useful tactic to advance our quantum “toolbox”. With this in mind, we consider the quantum calculation of a periodic system’s single-electron band structure over a path through reciprocal space. Previous efforts have used the Variational Quantum Eigensolver algorithm to solve the energy of each band, which involves numerically optimizing the parameters of a variational quantum circuit to minimize a cost function, constructed as the expectation value of a Hamiltonian operator. Traditionally, a unique Hamiltonian operator is constructed for each k-point, so that many cost functions, each with their own parameter space, must be optimized to generate a single band. Similarly, calculating higher bands than the first has traditionally involved modifying the cost function with additional overlap terms to ensure higher-energy eigenstates are orthogonal to those of lower bands. In this paper, we adopt a direct space approach, using a novel hybrid first/second-quantized qubit mapping which allows us to construct a single Hamiltonian, and a single cost-function, suitable for solving the entire electronic band structure. In contrast to previous approaches, the k-point and the band index are selected by additional parameters in our quantum circuit, rather than through modifications to the cost function. The result is a technically and conceptually simpler approach to band structure calculations on a quantum computer. Moreover, we expect that the tools developed herein will motivate new strategies for tackling highly-correlated materials beyond the grasp of classical computing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

The Milky Way tomography with APOGEE: intrinsic density distribution and structure of mono-abundance populations

ABSTRACT The spatial distribution of mono-abundance populations (MAPs, selected in [Fe/H] and [Mg/Fe]) reflect the chemical and structural evolution in a galaxy and impose strong constraints on galaxy formation models. In this paper, we use APOGEE data to derive the intrinsic density distribution of MAPs in the Milky Way, after carefully considering the survey selection function. We find that a single exponential profile is not a sufficient description of the Milky Way’s disc. Both the individual MAPs and the integrated disc exhibit a broken radial density distribution; densities are relatively constant with radius in the inner Galaxy and rapidly decrease beyond the break radius. We fit the intrinsic density distribution as a function of radius and vertical height with a 2D density model that considers both a broken radial profile and radial variation of scale height (i.e. flaring). There is a large variety of structural parameters between different MAPs, indicative of strong structure evolution of the Milky Way. One surprising result is that high-α MAPs show the strongest flaring. The young, solar-abundance MAPs present the shortest scale height and least flaring, suggesting recent and ongoing star formation confined to the disc plane. Finally we derive the intrinsic density distribution and corresponding structural parameters of the chemically defined thin and thick discs. The chemical thick and thin discs have local surface mass densities of 5.62 ± 0.08 and 15.69 ± 0.32 M⊙pc−2, respectively, suggesting a massive thick disc with a local surface mass density ratio between thick to thin disc of 36 per cent.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Correlating energy density induced residual stress, porosity, and mechanical property variations in directed energy deposition using neutron diffraction and imaging techniques

In a directed energy deposition (DED) process, the energy density deposited into the materials drives the structure and properties. The energy density is influenced by several printing parameters, such as laser power, hatch spacing, scanning speed, etc. In this study, 316L stainless steels samples were fabricated with varying hatch spacing aiming to induce variations in sample microstructure and properties. A combination of techniques including engineering neutron diffraction, neutron computed tomography, and tensile testing, were employed to correlate the structure and macroscopic properties. Porosity was predominantly observed at the bottom of the deposited materials and was effectively reduced by decreased hatch spacing. Correlatively, compressive residual stress was observed at the bottom of the specimens, while internal stress is largely determined by the hatch spacing. The residual stress in general decreases as the hatch spacing decreases. The maximum ultimate tensile strength (UTS) was found to increase from 517.4 to 528.6 MPa with a strain increased from 0.847 to 0.897 as the hatch spacing decreased from 0.45 to 0.42 mm. However, further increase in energy density by reducing the hatch spacing resulted in a significant decrease in UTS (as low as 452 MPa). This work provides new insights from bulk non-destructive techniques into correlating energy density induced residual stress, porosity, and mechanical property variations in directed energy deposition.

316L stainless steel↗
DOE OSTI2022

Tuning Metal–Organic Framework (MOF) Topology by Regulating Ligand and Secondary Building Unit (SBU) Geometry: Structures Built on 8-Connected M 6 (M = Zr, Y) Clusters and a Flexible Tetracarboxylate for Propane-Selective Propane/Propylene Separation

Topology evolution originating from variations of linker and SBU (Secondary Building Unit) geometries could largely enrich the chemistry of metal–organic frameworks (MOFs). Here we report the synthesis and characterization of three MOF structures built on the same organic linker, N,N,N',N'-Tetrakis(4-carboxyphenyl)-1,4-phenylenediamine (tcppda) and similar 8-connected M 6 (M = Zr or Y) clusters. The three compounds, HIAM-402, HIAM-403, and HIAM-311, feature 4,8-connected sqc, scu, and flu topology, respectively. Detailed structural analysis revealed that different geometries of the inorganic M 6 SBUs as well as the organic linker have led to the formation of distinct MOF nets. In particular, HIAM-402 features exceptional framework stability and high porosity and acts as a propane-selective adsorbent for the discrimination of propane and propylene. Its balanced adsorption selectivity (S propane/propylene = 1.43) and capacity (Q propane = 133.3 cm 3 /g, 298 K and 1 bar) endow it with the capability of separating propane and propylene mixtures and one-step production of highly pure propylene (purity >99.9%), as validated by column breakthrough measurements, with the presence of moisture or propyne. As a result, ab initio calculations further confirm that the propane-selective behavior of HIAM-402 is a result of its higher binding energy toward propane compared to that of propylene.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

SCOMAP-XD : atomistic deuterium contrast matching for small-angle neutron scattering in biology

The contrast-variation method in small-angle neutron scattering (SANS) is a uniquely powerful technique for determining the structure of individual components in biomolecular systems containing regions of different neutron scattering length density ρ . By altering the ρ of the target solute and the solvent through judicious incorporation of deuterium, the scattering of desired solute features can be highlighted. Most contrast-variation methods focus on highlighting specific bulk solute elements, but not on how the scattering at specific scattering vectors q , which are associated with specific structural distances, changes with contrast. Indeed, many systems exhibit q -dependent contrast effects. Here, a method is presented for calculating both bulk contrast-match points and q -dependent contrast using 3D models with explicit solute and solvent atoms and SASSENA , an explicit-atom SANS calculator. The method calculates the bulk contrast-match points within 2.4% solvent D 2 O accuracy for test protein–nucleic acid and lipid nanodisc systems. The method incorporates a general model for the incorporation of deuterium at non-exchangeable sites that was derived by performing mass spectrometry on green fluorescent protein. The method also decomposes the scattering profile into its component parts and identifies structural features that change with contrast. The method is readily applicable to a variety of systems, will expand the understanding of q -dependent contrast matching and will aid in the optimization of next-generation neutron scattering experiments.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Evolutionary lineage explains trait variation among 75 coexisting grass species

Summary Evolutionary history plays a key role driving patterns of trait variation across plant species. For scaling and modeling purposes, grass species are typically organized into C 3 vs C 4 plant functional types (PFTs). Plant functional type groupings may obscure important functional differences among species. Rather, grouping grasses by evolutionary lineage may better represent grass functional diversity. We measured 11 structural and physiological traits in situ from 75 grass species within the North American tallgrass prairie. We tested whether traits differed significantly among photosynthetic pathways or lineages (tribe) in annual and perennial grass species. Critically, we found evidence that grass traits varied among lineages, including independent origins of C 4 photosynthesis. Using a rigorous model selection approach, tribe was included in the top models for five of nine traits for perennial species. Tribes were separable in a multivariate and phylogenetically controlled analysis of traits, owing to coordination of important structural and ecophysiological characteristics. Our findings suggest grouping grass species by photosynthetic pathway overlooks variation in several functional traits, particularly for C 4 species. These results indicate that further assessment of lineage‐based differences at other sites and across other grass species distributions may improve representation of C 4 species in trait comparison analyses and modeling investigations.

Donnelly, Ryan C.↗
DOE OSTI2023

Accelerating full-waveform inversion using source stacking: synthetic experiments at the global scale in a realistic 3-D earth model

SUMMARY The spectral element method is currently the method of choice for computing accurate synthetic seismic wavefields in realistic 3-D earth models at the global scale. However, it requires significantly more computational time, compared to normal mode-based approximate methods. Source stacking, whereby multiple earthquake sources are aligned on their origin time and simultaneously triggered, can reduce the computational costs by several orders of magnitude. We present the results of synthetic tests performed on a realistic radially anisotropic 3-D model, slightly modified from model SEMUCB-WM1 with three component synthetic waveform ‘data’ for a duration of 10 000 s, and filtered at periods longer than 60 s, for a set of 273 events and 515 stations. We consider two definitions of the misfit function, one based on the stacked records at individual stations and another based on station-pair cross-correlations of the stacked records. The inverse step is performed using a Gauss–Newton approach where the gradient and Hessian are computed using normal mode perturbation theory. We investigate the retrieval of radially anisotropic long wavelength structure in the upper mantle in the depth range 100–800 km, after fixing the crust and uppermost mantle structure constrained by fundamental mode Love and Rayleigh wave dispersion data. The results show good performance using both definitions of the misfit function, even in the presence of realistic noise, with degraded amplitudes of lateral variations in the anisotropic parameter ξ. Interestingly, we show that we can retrieve the long wavelength structure in the upper mantle, when considering one or the other of three portions of the cross-correlation time series, corresponding to where we expect the energy from surface wave overtone, fundamental mode or a mixture of the two to be dominant, respectively. We also considered the issue of missing data, by randomly removing a successively larger proportion of the available synthetic data. We replace the missing data by synthetics computed in the current 3-D model using normal mode perturbation theory. The inversion results degrade with the proportion of missing data, especially for ξ, and we find that a data availability of 45 per cent or more leads to acceptable results. We also present a strategy for grouping events and stations to minimize the number of missing data in each group. This leads to an increased number of computations but can be significantly more efficient than conventional single-event-at-a-time inversion. We apply the grouping strategy to a real picking scenario, and show promising resolution capability despite the use of fewer waveforms and uneven ray path distribution. Source stacking approach can be used to rapidly obtain a starting 3-D model for more conventional full-waveform inversion at higher resolution, and to investigate assumptions made in the inversion, such as trade-offs between isotropic, anisotropic or anelastic structure, different model parametrizations or how crustal structure is accounted for.

Geochemistry & Geophysics↗