Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “waveforms”

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,261 records · Page 70
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Image-domain DAS 3D VSP elastic transmission tomography

SUMMARY Full-wavefield elastic imaging of active-source seismic data acquired by downhole receivers commonly offers higher-resolution subsurface images in the vicinity of a borehole compared to conventional surface seismic data sets, which can lack higher-frequency wavefield components due to longer travel paths and increased attenuation. An increasingly used approach for downhole acquisition is vertical seismic profiling (VSP), which has become more attractive when coupled with distributed acoustic sensing (DAS) using optical fibres installed in wells. The main difficulty for generating high-quality images with full-wavefield imaging tools for DAS VSP data, though, is the need for an accurate velocity model. To build plausible velocity models using active-source DAS VSP data, we adopt a 3-D image-domain elastic transmission tomography technique, originally developed for surface-recorded passive (microseismic) data, by exchanging the source and receiver positions (i.e. reciprocity) to mimic a passive-seismic surface monitoring scenario. The inversion approach exploits various images for each source constructed through time-reverse imaging (TRI) of downgoing P- and S-wave first-arrival waveforms. The TRI process uses the kinetic term of the (extended) PS energy imaging condition that exhibits sufficient sensitivity to velocity model errors. The method automatically updates the P- and S-wave velocity models to optimize image focusing via adjoint-state inversion. We illustrate the efficacy of the adopted elastic inversion technique using an active-source DAS 3-D VSP field data set acquired in the North Slope of Alaska. The numerical experiments demonstrate that the inverted elastic velocity models can be further used in full-wavefield acoustic/elastic imaging algorithms to obtain accurate subsurface images.

Oren, Can (ORCID:000000019534265X)↗
DOE OSTI2022

On the measurement of S diff splitting caused by lowermost mantle anisotropy

Seismic anisotropy has been detected at many depths of the Earth, including its upper layers, the lowermost mantle and the inner core. While upper mantle seismic anisotropy is relatively straightforward to resolve, lowermost mantle anisotropy has proven to be more complicated to measure. Due to their long, horizontal ray paths along the core–mantle boundary (CMB), S waves diffracted along the CMB (S diff ) are potentially strongly influenced by lowermost mantle anisotropy. S diff waves can be recorded over a large epicentral distance range and thus sample the lowermost mantle everywhere around the globe. S diff therefore represents a promising phase for studying lowermost mantle anisotropy; however, previous studies have pointed out some difficulties with the interpretation of differential SH diff –SV diff traveltimes in terms of seismic anisotropy. Here, we provide a new, comprehensive assessment of the usability of S diff waves to infer lowermost mantle anisotropy. Using both axisymmetric and fully 3-D global wavefield simulations, we show that there are cases in which S diff can reliably detect and characterize deep mantle anisotropy when measuring traditional splitting parameters (as opposed to differential traveltimes). First, we analyze isotropic effects on S diff polarizations, including the influence of realistic velocity structure (such as 3-D velocity heterogeneity and ultra-low velocity zones), the character of the lowermost mantle velocity gradient, mantle attenuation structure, and Earth’s Coriolis force. Secondly, we evaluate effects of seismic anisotropy in both the upper and the lowermost mantle on SH diff waves. In particular, we investigate how SH diff waves are split by seismic anisotropy in the upper mantle near the source and how this anisotropic signature propagates to the receiver for a variety of lowermost mantle models. We demonstrate that, in particular and predictable cases, anisotropy leads to S diff splitting that can be clearly distinguished from other waveform effects. Here, these results enable us to lay out a strategy for the analysis of S diff splitting due to anisotropy at the base of the mantle, which includes steps to help avoid potential pitfalls, with attention paid to the initial polarization of S diff and the influence of source-side anisotropy. We demonstrate our S diff splitting method using three earthquakes that occurred beneath the Celebes Sea, measured at many transportable array stations at a suitable epicentral distance. We resolve consistent and well-constrained S diff splitting parameters due to lowermost mantle anisotropy beneath the northeastern Pacific Ocean.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

High-resolution source imaging and moment tensor estimation of acoustic emissions during brittle creep of basalt undergoing carbonation

SUMMARY As the high-frequency analogue to field-scale earthquakes, acoustic emissions (AEs) provide a valuable complement to study rock deformation mechanisms. During the load-stepping creep experiments with CO2-saturated water injection into a basaltic sample from Carbfix site in Iceland, 8791 AE events are detected by at least one of the seven piezoelectric sensors. Here, we apply a cross-correlation-based source imaging method, called geometric-mean reverse-time migration (GmRTM) to locate those AE events. Besides the attractive picking-free feature shared with other waveform-based methods (e.g. time-reversal imaging), GmRTM is advantageous in generating high-resolution source images with reduced imaging artefacts, especially for experiments with relatively sparse receivers. In general, the imaged AE locations are found to be scattered across the sample, suggesting a complicated fracture network rather than a well-defined major shear fracture plane, in agreement with X-ray computed tomography imaging results after retrieval of samples from the deformation apparatus. Clustering the events in space and time using the nearest-neighbour approach revealed a group of ‘repeaters’, which are spatially co-located over an elongated period of time and likely indicate crack, or shear band growth. Furthermore, we select 2196 AE events with high signal-to-noise-ratio (SNR) and conduct moment tensor estimation using the adjoint (backpropagated) strain tensor fields at the locations of AE sources. The resulting AE locations and focal mechanisms support our previously assertion that creep of basalt at the experimental conditions is accommodated dominantly by distributed microcracking.

58 GEOSCIENCES↗
DOE OSTI2024

Numerical modelling of impact seismic sources using the stress glut theory

SUMMARY Meteorite impacts have proved to be a significant source of seismic signal on the Moon, and have now been recorded on Mars by InSight seismometers. Understanding how impacts produce seismic signal is key to the interpretation of this unique data, and to improve their identification in continuous seismic records. Here, we use the seismic Representation Theorem, and particularly the stress glut theory, to model the seismic motion resulting from impact cratering. The source is described by equivalent forces, some resulting from the impactor momentum transfer, and others from the stress glut, which represents the mechanical effect of plasticity and non linear processes in the source region. We condense these equivalent forces into a point-source with a time-varying single force and nine-component moment tensor. This analytical representation bridges the gap between the complex dynamics of crater formation, and the linear point-source representation classically used in seismology. Using the multiphysics modelling software HOSS, we develop a method to compute the stress glut of an impact, and the associated point-source from hypervelocity impact simulations. For a vertical and an oblique impact at 1000 m s−1, we show that the moment tensor presents a significant deviatoric component. Hence, the source is not an ideal isotropic explosion contrary to previous assumptions, and draws closer to a double couple for the oblique impact. The contribution of the point force to the seismic signal appears negligible. We verify this model by comparing two signals: (1) HOSS is coupled to SPECFEM3D to propagate the near-source signal elastically to remote seismic stations; (2) the point-source model derived from the stress-glut theory is used to generate displacements at the same distance. The comparison shows that the point-source model is accurately simulating the low-frequency impact seismic waveform, and its seismic moment is in trend with Lunar and Martian impact data. High-frequencies discrepancies exist, which are partly related to finite-source effects, but might be further explained by the difference in mathematical framework between classical seismology and HOSS’ numerical modelling.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Comprehensive earthquake catalogue update and spatiotemporal distribution analysis for Iraq and surrounding regions, northeastern Arabian Plate

The updated earthquake catalogue for Iraq covers the period from 1900 to the end of 2021 and includes over 37 000 recorded earthquakes. To create this comprehensive catalogue, five key steps were taken: compiling bulletins, calculating moment magnitudes, harmonizing magnitudes, establishing empirical conversion relations and evaluating the completeness of the catalogue. A notable enhancement in this update is the direct calculation of moment magnitudes for approximately 2800 earthquakes, achieved through the coda envelope technique and waveform data from the Mesopotamian Seismological Network (MPSN) in Iraq. This updated catalogue serves as a valuable resource for examining the spatiotemporal distribution of earthquakes, with respect to earthquake density, maximum moment magnitude and seismogenic depths. Additionally, the Gutenberg–Richter relationship was applied to calculate the a- and b-values specific to Iraq. The findings show that the Zagros Fold-Thrust Belt has a seismogenic layer (source) that ranges from 2 to 33 km deep and experiences high seismic activity. In contrast, the Mesopotamian Foredeep has a seismogenic layer ranging from 1 to 25 km deep and has lower seismic activity. The greatest seismic activity is concentrated around the Mandili-Badra-Teeb fault, which has experienced significant ruptures over time. The Outer Arabian Platform is identified as the main area of seismic activity, while additional activity occurs on the Inner Arabian Platform. Three major tectonic boundaries define the distribution of earthquakes in the northeastern Arabian Plate. These boundaries are defined by the Main Zagros Reverse Fault, the Zagros Foredeep Fault and the Anah Graben and Abu Jir-Euphrates Fault Zone. These boundaries highlight variations in seismicity levels and the spatial distribution of deformation in the region. The updated earthquake catalogue presented in this study is expected to play a vital role in regional seismicity assessments and seismic hazard analyses for Iraq and its surrounding areas.

58 GEOSCIENCES↗
DOE OSTI2015

Discovering new events beyond the catalogue—application of empirical matched field processing to Salton Sea geothermal field seismicity

Using empirical matched field processing (MFP), we compare 4 yr of continuous seismic data to a set of 195 master templates from within an active geothermal field and identify over 140 per cent more events than were identified using traditional detection and location techniques alone. In managed underground reservoirs, a substantial fraction of seismic events can be excluded from the official catalogue due to an inability to clearly identify seismic-phase onsets. Empirical MFP can improve the effectiveness of current seismic detection and location methodologies by using conventionally located events with higher signal-to-noise ratios as master events to define wavefield templates that could then be used to map normally discarded indistinct seismicity. Since MFP does not require picking, it can be carried out automatically and rapidly once suitable templates are defined. In this application, we extend MFP by constructing local-distance empirical master templates using Southern California Earthquake Data Center archived waveform data of events originating within the Salton Sea Geothermal Field. We compare the empirical templates to continuous seismic data collected between 1 January 2008 and 31 December 2011. The empirical MFP method successfully identifies 6249 additional events, while the original catalogue reported 4352 events. The majority of these new events are lower-magnitude events with magnitudes between M0.2–M0.8. Here, the increased spatial-temporal resolution of the microseismicity map within the geothermal field illustrates how empirical MFP, when combined with conventional methods, can significantly improve seismic network detection capabilities, which can aid in long-term sustainability and monitoring of managed underground reservoirs.

15 GEOTHERMAL ENERGY↗
DOE OSTI2016

Frequency domain analysis of errors in cross-correlations of ambient seismic noise

We analyse random errors (variances) in cross-correlations of ambient seismic noise in the frequency domain, which differ from previous time domain methods. Extending previous theoretical results on ensemble averaged cross-spectrum, we estimate confidence interval of stacked cross-spectrum of finite amount of data at each frequency using non-overlapping windows with fixed length. The extended theory also connects amplitude and phase variances with the variance of each complex spectrum value. Analysis of synthetic stationary ambient noise is used to estimate the confidence interval of stacked cross-spectrum obtained with different length of noise data corresponding to different number of evenly spaced windows of the same duration. This method allows estimating Signal/Noise Ratio (SNR) of noise cross-correlation in the frequency domain, without specifying filter bandwidth or signal/noise windows that are needed for time domain SNR estimations. Based on synthetic ambient noise data, we also compare the probability distributions, causal part amplitude and SNR of stacked cross-spectrum function using one-bit normalization or pre-whitening with those obtained without these pre-processing steps. Natural continuous noise records contain both ambient noise and small earthquakes that are inseparable from the noise with the existing pre-processing steps. Using probability distributions of random cross-spectrum values based on the theoretical results provides an effective way to exclude such small earthquakes, and additional data segments (outliers) contaminated by signals of different statistics (e.g. rain, cultural noise), from continuous noise waveforms. This technique is applied to constrain values and uncertainties of amplitude and phase velocity of stacked noise cross-spectrum at different frequencies, using data from southern California at both regional scale (~35 km) and dense linear array (~20 m) across the plate-boundary faults. A block bootstrap resampling method is used to account for temporal correlation of noise cross-spectrum at low frequencies (0.05–0.2 Hz) near the ocean microseismic peaks.

58 GEOSCIENCES↗
DOE OSTI2016

Fast computation of synthetic seismograms within a medium containing remote localized perturbations: a numerical solution to the scattering problem

We derive a fast discrete solution to the scattering problem. This solution allows us to compute accurate synthetic seismograms or waveforms for arbitrary locations of sources and receivers within a medium containing localized perturbations. The key to efficiency is that wave propagation modelling does not need to be carried out in the entire volume that encompasses the sources and the receivers but only within the subvolume containing the perturbations or scatterers. The proposed solution has important applications, for example, it permits the imaging of remote targets located in regions where no sources or receivers are present. Our solution relies on domain decomposition: within a small volume that contains the scatterers, wave propagation is modelled numerically, while in the surrounding volume, where the medium isn't perturbed, the response is obtained through wavefield extrapolation. The originality of this work is the derivation of discrete formulas for representation theorems and Kirchhof–Helmholtz integrals that naturally adapt to the numerical scheme employed for modelling wave propagation. Our solution applies, for example, to finite difference methods or finite/spectral elements methods. The synthetic seismograms obtained with our solution can be considered ‘exact’ as the total numerical error is comparable to that of the method employed for modelling wave propagation. We detail a basic implementation of our solution in the acoustic case using the finite difference method and present numerical examples that demonstrate the accuracy of the method. We show that ignoring some terms accounting for higher order scattering effects in our solution has a limited effect on the computed seismograms and significantly reduces the computational effort. Finally, we show that our solution can be used to compute localized sensitivity kernels and we discuss applications to target oriented imaging. Extension to the elastic case is straightforward and summarized in a dedicated section.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Source characteristics and geological implications of the January 2016 induced earthquake swarm near Crooked Lake, Alberta

On 2016 January 12, an intraplate earthquake with an initial reported local magnitude (M L ) of 4.8 shook the town of Fox Creek, Alberta. While there were no reported damages, this earthquake was widely felt by the local residents and suspected to be induced by the nearby hydraulic-fracturing (HF) operations. We determine the earthquake source parameters using moment tensor inversions, and then detect and locate the associated swarm using a waveform cross-correlation based method. The broad-band seismic recordings from regional arrays suggest a moment magnitude (M) 4.1 for this event, which is the largest in Alberta in the past decade. Similar to other recent M ~ 3 earthquakes near Fox Creek, the 2016 January 12 earthquake exhibits a dominant strike-slip (strike = 184°) mechanism with limited non-double-couple components (~22 per cent). This resolved focal mechanism, which is also supported by forward modelling and P-wave first motion analysis, indicates an NE–SW oriented compressional axis consistent with the maximum compressive horizontal stress orientations delineated from borehole breakouts. Further detection analysis on industry-contributed recordings unveils 1108 smaller events within 3 km radius of the epicentre of the main event, showing a close spatial-temporal relation to a nearby HF well. The majority of the detected events are located above the basement, comparable to the injection depth (3.5 km) on the Duvernay shale Formation. The spatial distribution of this earthquake cluster further suggests that (1) the source of the sequence is an N–S-striking fault system and (2) these earthquakes were induced by an HF well close to but different from the well that triggered a previous (January 2015) earthquake swarm. Reactivation of pre-existing, N-S oriented faults analogous to the Pine Creek fault zone, which was reported by earlier studies of active source seismic and aeromagnetic data, are likely responsible for the occurrence of the January 2016 earthquake swarm and other recent events in the Crooked Lake area.

58 GEOSCIENCES↗
DOE OSTI2018

Acoustic Event Location and Background Noise Characterization on a Free Flying Infrasound Sensor Network in the Stratosphere

We present that a variety of Earth surface and atmospheric sources generate low frequency sound waves that can travel great distances. Despite a rich history of ground-based sensor studies, very few experiments have investigated the prospects of free floating microphone arrays at high altitudes. However, recent initiatives have shown that such networks have very low background noise and may sample an acoustic wave field that is fundamentally different than that at Earth’s surface. The experiments have been limited to at most two stations at altitude, making acoustic event detection and localization difficult. We describe the deployment of four drifting microphone stations at altitudes between 21 and 24 km above sea level. The stations detected one of two regional ground-based chemical explosions as well as the ocean microbarom while traveling almost 500 km across the American Southwest. The explosion signal consisted of multiple arrivals; signal amplitudes did not correlate with sensor elevation or source range. The waveforms and propagation patterns suggest interactions with gravity waves in the 35-45 km altitude. A sparse network method that employed curved wave front corrections was able to determine the backazimuth from the free flying network to the acoustic source. Episodic signals similar to those seen on previous flights in the same region were noted, but their source remains unclear. Lastly, background noise levels were commensurate with those on infrasound stations in the International Monitoring System below 2 seconds.

58 GEOSCIENCES↗
DOE OSTI2018

Shallow three-dimensional structure of the San Jacinto fault zone revealed from ambient noise imaging with a dense seismic array

In this work, we use 1 month of continuous seismic waveforms from a very dense seismic network to image with unprecedented resolution the shallow damage structure of the San Jacinto fault zone across the Clark fault strand. After calculating noise correlations, high apparent velocity arrivals coming from below the array are removed using a frequency–wavenumber filter. This is followed by a double-beamforming analysis on multiple pairs of subarrays to extract phase and group velocity information across the study area. The phase and group velocity dispersion curves are regionalized into phase and group velocity maps at different frequencies, and these maps are inverted using a neighbourhood algorithm to build a 3-D shear wave velocity model around the Clark fault down to ~500 m depth. The model reveals strong lateral variations across the fault strike with pronounced low-velocity zones corresponding to a local sedimentary basin and a fault zone trapping structure. The results complement previous earthquake- and seismic noise-based imaging of the fault zone at greater depth and clarify properties of structural features near the surface.

58 GEOSCIENCES↗
DOE OSTI2021

Projecting the likely importance of weak-interaction-driven bulk viscosity in neutron star mergers

ABSTRACT In this work, we estimate how much bulk viscosity driven by Urca processes is likely to affect the gravitational wave signal of a neutron star coalescence. In the late inspiral, we show that bulk viscosity affects the binding energy at fourth post-Newtonian order. Even though this effect is enhanced by the square of the gravitational compactness, the coefficient of bulk viscosity is likely too small to lead to observable effects in the waveform during the late inspiral, when only considering the orbital motion itself. In the post-merger, however, the characteristic time-scales and spatial scales are different, potentially leading to the opposite conclusion. We post-process data from a state-of-the-art equal-mass binary neutron star merger simulation to estimate the effects of bulk viscosity (which was not included in the simulation itself). In that scenario, we find that bulk viscosity can reach high values in regions of the merger. We compute several estimates of how much it might directly affect the global dynamics of the considered merger scenario, and find that it could become significant. Even larger effects could arise in different merger scenarios or in simulations that include non-linear effects. This assessment is reinforced by a quantitative comparison with relativistic heavy-ion collisions where such effects have been explored extensively.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Optimizing infrasound observations for sample return capsule re-entry: Insights from OSIRIS-REx and Hayabusa2

The atmospheric entry of meteoroids presents a rare and unpredictable phenomenon, posing challenges for systematic observation and detailed characterization. Such events are nonetheless critical for advancing understanding of acoustic wave propagation, atmospheric structure, and entry dynamics. In contrast, sample return capsules (SRCs) from space missions follow well constrained re-entry trajectories, enabling planned observations of shock wave generation and propagation under controlled conditions. This study compares two SRC atmospheric entries, Hayabusa2 in 2020 and OSIRIS-REx in 2023, to assess how different infrasound array configurations influence shock wave detection and trajectory validation. Hayabusa2’s re-entry was monitored using a distributed network of 28 portable infrasound sensors across seven arrays in Woomera, Australia, permitting three-dimensional reconstruction of the trajectory and analysis of wave characteristics. For OSIRIS-REx, a compact four-sensor array deployed near Eureka Airport provided trajectory confirmation through arrival-time differences and back-azimuth estimates. Spectral and waveform analyses revealed differences in signal properties associated with variations in entry angle and velocity. The results illustrate both the strengths and the limitations of the deployed array configurations. The dense and distributed arrays during Hayabusa2’s re-entry enabled detailed trajectory reconstruction, whereas the compact array at Eureka primarily provided confirmation of signal coherence and back-azimuth consistency. These case studies highlight that even relatively small arrays, if located close to the predicted ground track, can still capture useful information on arrival direction and timing. Such insights provide practical guidance for planning future observational campaigns of SRC returns and other controlled atmospheric entries, and may also inform approaches to opportunistic observations of meteoroid events and other atmospheric acoustic phenomena.

celestial mechanics: orbit determination↗
DOE OSTI2019

Tracking DO 2 with Compensatory Reserve during Whole Blood Resuscitation in Baboons

Background Hemorrhagic shock can be mitigated by timely and accurate resuscitation designed to restore adequate delivery of oxygen (DO 2 ) by increasing cardiac output (CO). However, standard care of using systolic blood pressure (SBP) as a guide for resuscitation has proven ineffective and can be associated with increased morbidity. We have developed a novel vital sign called the compensatory reserve (CRM) generated from analysis of arterial pulse waveform feature changes that has been validated in experimental and clinical models of hemorrhage. We tested the hypothesis that thresholds of DO 2 could be accurately defined by CRM while avoiding over resuscitation during whole blood resuscitation following a 25% hemorrhage in non‐human primates. Methods To accomplish this, adult male baboons (n=12) were exposed to a progressive controlled hemorrhage that resulted in an average (± SD) maximal reduction of 601 ± 18 ml of their estimated circulating blood volume of 1,262 ± 64 ml based on body weight. Values of DO 2 were calculated as the product of CO times O 2 carrying capacity of hemoglobin (i.e., Hb content × 1.34). Results Average (± SD) CRM increased from 6 ± 5 % at the end of hemorrhage to 70 ± 7 % at the end of resuscitation. DO 2 of 12.1, 10.3, and 7.0 ml O 2 ·kg −1 ·min −1 were calculated at 0% (baseline), 12.5% and 25% blood loss. By linear regression, CRM values of 6% (decompensation), 35%, and 68% corresponded to calculated DO 2 values of 6.0, 9.2, and 12.1 ml O 2 ·kg −1 ·min −1 during resuscitation. As such, return of CRM to baseline during resuscitation with only ~ 450 ml provided adequate DO 2 (12.5 ml O 2 ·kg −1 ·min −1 ) while total blood volume replacement (~600 ml) resulted in over resuscitation as indicated by an average increase in SBP (+17 mmHg) and CO (+0.1 L/min) over baseline values. Conclusion Consistent with our hypothesis, thresholds of DO 2 were associated with specific CRM values. A target resuscitation CRM value of 70% minimized the requirement for whole blood, while avoiding over resuscitation. Furthermore, CRM at 0% provided a noninvasive metric for determining critical DO 2 at approximately 5.4 ml O 2 ·kg −1 ·min −1 . Support or Funding Information This research was supported in part by appointment of Ms. Koons to the Postgraduate Research Fellowship Program at the US Army Institute of Surgical Research administered by the Oak Ridge Institute for Science and Education through an interagency agreement between the US Department of Energy and USAMRMC. Study funding was provided by a grant from the US Army Combat Casualty Care Research Program (D‐009‐2014‐USAISR). Linearregression of the relationship between DO 2 and compensatory reserve during whole blood resuscitation in baboons. Data are expressed as mean ± SE. image Linearregression of the relationship between DO 2 and compensatory reserve during whole blood resuscitation in baboons. Data are expressed as mean ± SE. This abstract is from the Experimental Biology 2019 Meeting. There is no full text article associated with this abstract published in The FASEB Journal .

Koons, Natalie Jane↗
DOE OSTI2025

Theory of terahertz pulse transmission through ferroelectric nanomembranes

An analytical model is developed to predict the temporal evolution of the lattice polarization in ferroelectric nanomembranes upon the excitation by a terahertz (THz) electromagnetic pulse of an arbitrary waveform and the concurrent transmission of the THz pulse in both linear and nonlinear regimes. It involves the use of the perturbation method to solve the equation of motion for the lattice polarization in both unclamped and strained ferroelectric nanomembranes within the framework of Landau-Ginzburg-Devonshire theory. The model is applicable to perovskite oxides such as BaTiO 3 and SrTiO 3 , wurtzite Al 1−𝑥 ⁢Sc 𝑥 ⁢N, and trigonal LiNbO 3 . Our analytical model provides a theoretical basis for determining the thermodynamic and kinetic parameters of ferroelectric materials through a THz transmission experiment. The calculation results also suggest an approach to reversing the chirality of a circularly polarized THz pulse by harnessing the resonant polarization-photon coupling in ferroelectrics. This capability of chirality reversal, along with the high tunability from a strain applied along any arbitrarily oriented in-plane axis, provides new opportunities for THz wave modulation without relying on complex metasurface designs.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Laser-Modulation-Driven X-Ray Pulse Shaping in Regenerative Amplifier Free-Electron Lasers

In this Letter, we present a robust method for generating custom-shaped, coherent hard x-ray pulses in regenerative amplifier free-electron lasers (RAFELs) using laser-induced energy modulation of the electron beam. A temporally shaped optical modulation imprints an optical-wavelength energy pattern on the beam, which interacts with the recirculated seed in the x-ray cavity to generate spectral sidebands. These sidebands lie outside the Bragg reflection bandwidth and are extracted as out-coupled radiation, while the seed is retained for continued amplification, enabling a novel mechanism for both output coupling and pulse shaping. Furthermore, this approach allows deterministic control over the spectral and temporal properties of the emitted x-rays, including multicolor generation, twin-pulse structures, pulse duration tuning, and complex temporal waveforms. Importantly, it does not require modification of the RAFEL cavity or optics and is fully compatible with high-repetition-rate operation, offering a practical pathway to deliver coherent, shaped, and on-demand x-ray pulses for advanced time-resolved and nonlinear x-ray science.

Free-electron lasers↗
DOE OSTI2025

Supertranslations from Scattering Amplitudes

On shell methods have found a new application to local observables such as asymptotic radiation fields and gravitational waveforms. While these observables are invariant under small gauge transformations, they are known to depend on a choice of asymptotic gauge; in gravity on asymptotically Minkowski spacetimes, this is a choice of Bondi-Van der Burg-Metzner-Sachs frame. In this Letter, we provide a method for capturing these supertranslations, to all orders in perturbation theory, using the on shell framework of scattering amplitudes.

Classical black holes↗