Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “transient detection”

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,405 records · Page 78

The rebrightening of a ROSAT -selected tidal disruption event: repeated weak partial disruption flares from a quiescent galaxy?

The ROSAT-selected tidal disruption event (TDE) candidate RX J133157.6-324319.7 (J1331) was detected in 1993 as a bright [0.2–2 keV flux of (1.0 ± 0.1) × 10 -12 erg s -1 cm -2 ], ultra-soft (kT = 0.11 ± 0.03 keV) X-ray flare from a quiescent galaxy (z = 0.051 89). During its fifth all-sky survey (eRASS5) in 2022, Spectrum-Roentgen-Gamma (SRG)/ eROSITA detected the repeated flaring of J1331, where it had rebrightened to an observed 0.2–2 keV flux of (6.0 ± 0.7) × 10 -13 erg s -1 cm -2 , with spectral properties (kT = 0.115 ± 0.007 keV) consistent with the ROSAT-observed flare ~30 yr earlier. In this work, we report on X-ray, ultraviolet, optical, and radio observations of this system. During a pointed XMM observation ~17 d after the eRASS5 detection, J1331 was not detected in the 0.2–2 keV band, constraining the 0.2–2 keV flux to have decayed by a factor of ≳40 over this period. Given the extremely low probability (~5 × 10 -6 ) of observing two independent full TDEs from the same galaxy over a 30 yr period, we consider the variability seen in J1331 to be likely caused by two partial TDEs involving a star on an elliptical orbit around a black hole. J1331-like flares show faster rise and decay time-scales [$\mathcal {O}(\mathrm{d})$] compared to standard TDE candidates, with negligible ongoing accretion at late times post-disruption between outbursts.

79 ASTRONOMY AND ASTROPHYSICS↗

Robust Carbon Dioxide Plume Imaging Using Joint Tomographic Inversion of Seismic Onset Time and Distributed Pressure and Temperature Measurements (Final Report)

We develop and demonstrate rapid and cost-effective methodologies for spatiotemporal tracking of CO2 plumes during geologic sequestration using joint inversion of seismic data and distributed pressure and temperature measurements. Key elements of our methodology are: (a) a computationally efficient approach to pressure and temperature propagation, (b) analysis of time lapse seismic data using a novel ‘seismic onset time’ approach to detect fluid front propagation, and (c) data assimilation and uncertainty assessment via joint inversion of pressure, temperature and time lapse seismic data, and (d) validating the numerical tomographic inversion using a CO2 injection demonstration projects, specifically data collected from the from the Petra Nova Parish Holdings CCUS project in the West Ranch Field, Texas and the Chester-16 reef CO2 injection site in Northern Michigan which is part of the DOE Midwestern Carbon Sequestration Project. The research team is led by Texas A&M University and includes Battelle as a subcontractor with support from Shell, Anadarko, Chevron and JX Nippon. A carbon dioxide (CO2) water-alternating-gas (WAG) pilot was conducted to gain insights into tertiary oil recovery potential via CO2 flood in the West Ranch Field as part of the Petra Nova project, the world’s largest post-combustion CO2 capture and utilization initiative. With a fluvial formation geology and large contrasts in permeability, this is a challenging and novel application of CO2 enhanced oil recovery (EOR). We build a predictive dynamic model of the subsurface that incorporates the multiphase and compositional data acquired during the pilot operation. The calibrated model is used for the carbon dioxide plume imaging. The study began with an initialization of the pilot sector model extracted from a calibrated full-field model. The pilot model calibration follows a two-step hierarchical workflow. First, we performed a large-scale update of the permeability distribution by integrating available bottomhole pressure and multiphase production data. In the second step, local permeability field is fine-tuned using a streamline-based method to match CO2 breakthrough times at the producers. The predictive capability of the calibrated model was verified through two blind validation tests: (1) the model showed good agreement with saturation logs acquired at two observation wells; and (2) the model reproduced the CO2 recovery as a fraction of the injected CO2. The use of seismic onset times has shown great promise for integrating near-continuous seismic surveys for updating geologic models. In this study, we analyze the impact of seismic survey frequency on the onset time approach aiming to extend the application of onset time to infrequent seismic surveys. In addition, we quantitatively examine the nonlinearity of the onset time method and compare it to the commonly used amplitude inversion method. We carry out a sensitivity analysis of seismic survey frequency based on the complete seismic survey data (over 175 surveys) of steam injection in a heavy oil reservoir (Peace River Unit) in Canada. Our results show that an adequate onset time map can be obtained from the infrequent seismic surveys by interpolation between seismic surveys as long as there is no change in the dominant underlying physics between the successive surveys. The study also shows that nonlinearity of the onset time method can be -smaller than that of the amplitude inversion method by several orders of magnitude. Application to the Brugge benchmark case shows that the onset time method obtains comparable permeability update as the traditional seismic amplitude inversion method with faster computation and improved convergence characteristics. We extend the streamline-based data integration approach to incorporate distributed temperature sensor (DTS) data using the concept of thermal tracer travel time. Then, a hierarchical workflow composed of evolutionary and streamline methods is employed to jointly history match the DTS and pressure data. Finally, CO2 saturation and streamline maps are used to visualize the CO2 plume movement during the sequestration process. The hierarchical workflow is applied to a carbon sequestration project in a carbonate reef reservoir within the Northern Niagaran Pinnacle Reef Trend in Michigan, USA. The monitoring data set consists of distributed temperature sensing (DTS) data acquired at the injection well and a monitoring well, flowing bottom-hole pressure data at the injection well, and time-lapse pressure measurements at several locations along the monitoring well. The history matching results indicate that the CO2 movement is mostly restricted to the intended zones of injection which is consistent with an independent warm-back analysis of the temperature data. In addition to employing simulation models and inverse methods for CO2 plume imaging, we also initialized a data-driven technology for detecting inter-well connectivity based on production and pressure data. Our machine-learning framework is built on the statistical recurrent unit (SRU) model and interprets well-based injection/production data into inter-well connectivity without relying on a geologic model. We test it on synthetic and field-scale CO2 EOR projects utilizing the water-alternating-gas (WAG) process. The validation of the proposed data-driven inter-well connectivity assessment is performed using synthetic data from simulation models where inter-well connectivity can be easily measured using the streamline-based flux allocation. The SRU model is shown to offer excellent prediction performance on the synthetic case. Despite significant measurement noise and frequent well shut-ins imposed in the field-scale case, the SRU model offers good prediction accuracy, the overall relative error of the phase production rates at most producers ranges from 10% to 30%. It is shown that the dominant connections identified by the data-driven method and streamline method are in close agreement. Texas A&M University, the lead organization in the project, was primarily responsible for the development of tomographic approaches for CO2 plume mapping in conjunction with distributed pressure, temperature and seismic onset time data. Battelle, as a subcontractor, was primarily responsible for the development of analytical and empirical methods for analyzing transient injection rate and pressure data from point/line sources such as injection and monitoring wells. An additional area of emphasis for Battelle was the use of machine learning for such tasks as inferring reservoir connectivity information from injection-production data, and identifying variable importance for machine learning-based proxy models developed from full-physics simulations. The two organizations also collaborated on the application of the tomographic inversion methodology for a field data set.

02 PETROLEUM↗

Core collapse supernova gravitational wave emission for progenitors of 9.6, 15, and 25 M ⊙

Here, we present gravitational wave emission predictions based on three core collapse supernova simulations corresponding to three different progenitor masses. The masses span a large range, between 9.6 and 25M ⊙ , are all initially nonrotating, and are of two metallicities: zero and solar. We compute both the temporal evolution of the gravitational wave strains for both the plus and the cross polarizations, as well as their spectral decomposition and characteristic strains. The temporal evolution of our zero metallicity 9.6M ⊙ progenitor model is distinct from the temporal evolution of our solar metallicity 15M ⊙ progenitor model and our zero metallicity 25M ⊙ progenitor model. In the former case, the high-frequency gravitational wave emission is largely confined to a brief time period ~75 ms after bounce, whereas in the latter two cases high-frequency emission does not commence until ~125 ms after bounce or later. The excitation mechanisms of the high-frequency emission in all three cases correspond to proto-neutron star convection and accretion onto the proto-neutron star from the convective gain layer above it, with the former playing the dominant role for most of the evolution. The low-frequency emission in all three models exhibits very similar behavior. At frequencies below ~250 Hz, gravitational waves are emitted by neutrino-driven convection and the standing accretion shock instability (SASI). This emission extends throughout the simulations when a gain region is present. In all three models, explosion is observed at ~125, ~500, and ~250 ms after bounce in the 9.6, 15, and 25M ⊙ progenitor models, respectively. At these times, the low-frequency gravitational wave emission is joined by very low-frequency emission, below ~10 Hz. These very low-frequency episodes are the result of explosion and begin at the above designated explosion times in each of our models. Our characteristic strains tell us that, in principle, all three gravitational wave signals would be detectable by current-generation detectors for a supernova at a distance of 10 kpc. However, our 9.6M ⊙ progenitor model is a significantly weaker source of gravitational waves, with strain amplitudes approximately 5–10 times less than in our other two models. The characteristic strain for this model tells us that such a supernova would be detectable only within a much more narrow frequency range around the maximum sensitivity of today’s detectors. Finally, in our 9.6M ⊙ progenitor model, we see very high-frequency gravitational radiation, extending up to ~2000 Hz. This feature results from the interaction of shock- and deleptonization-induced convection with perturbations introduced in the progenitor by nuclear burning during core collapse. While unique to the 9.6M ⊙ progenitor model analyzed here, this very high-frequency emission may, in fact, be a generic feature of the predictions for the gravitational wave emission from all core collapse supernova models when simulations are performed with three-dimensional progenitors.

79 ASTRONOMY AND ASTROPHYSICS↗

Active Galactic Nuclei Continuum Reverberation Mapping Based on Zwicky Transient Facility Light Curves

We perform a systematic survey of active galactic nuclei (AGNs) continuum lags using ~3 days cadence gri-band light curves from the Zwicky Transient Facility. We select a sample of 94 type 1 AGNs at z < 0.8 with significant and consistent inter-band lags based on the interpolated cross-correlation function method and the Bayesian method JAVELIN. Within the framework of the "lamp-post" reprocessing model, our findings are: (1) The continuum emission (CE) sizes inferred from the data are larger than the disk sizes predicted by the standard thin-disk model. (2) For a subset of the sample, the CE size exceeds the theoretical limit of the self-gravity radius (12 lt-days) for geometrically thin disks. (3) The CE size scales with continuum luminosity as R CE ∝ L 0.48±0.04 with a scatter of 0.2 dex, analogous to the well-known radius–luminosity relation of broad Hβ. These findings suggest a significant contribution of diffuse continuum emission from the broad-line region (BLR) to AGN continuum lags. We find that the R CE –L relation can be explained by a photoionization model that assumes ~23% of the total flux comes from the diffuse BLR emission. In addition, the ratio of the CE size and model-predicted disk size anticorrelates with the continuum luminosity, which is indicative of a potential nondisk BLR lag contribution evolving with the luminosity. Finally, a robust positive correlation between the CE size and black hole mass is detected.

79 ASTRONOMY AND ASTROPHYSICS↗

Plumelets: Dynamic Filamentary Structures in Solar Coronal Plumes

Solar coronal plumes long seemed to possess a simple geometry supporting spatially coherent, stable outflow without significant fine structure. Recent high-resolution observations have challenged this picture by revealing numerous transient, small-scale, collimated out-flows (“jetlets”) at the base of plumes. The dynamic filamentary structure of solar plumes above these outflows, and its relationship with the overall plume structure, have remained largely unexplored. We analyzed the statistics of continuously observed fine structure in-side a single representative bright plume within a mid-latitude coronal hole during 2016 July2-3. By applying advanced edge-enhancement and spatiotemporal analysis techniques to ex-tended series of high-resolution images from the Solar Dynamics Observatory’s Atmospheric Imaging Assembly, we determined that the plume was composed of numerous time-evolving filamentary substructures, referred to as “plumelets” in this paper, that accounted for most of the plume emission. The number of simultaneously identifiable plumelets was positively correlated with plume brightness, peaked in the fully formed plume, and remained saturated thereafter. The plumelets had transverse widths of 10 Mm and intermittently supported upwardly propagating periodic disturbances with phase speeds of 190-260 km s−1and longitudinal wavelengths of 55-65 Mm. The characteristic frequency (3.5 mHz) is commensurate with that of solar p-modes. Oscillations in neighboring plumelets are uncorrelated, indicating that the waves could be driven by p-mode flows at spatial scales smaller than the plumelet separation. Multiple independent sources of outflow within a single coronal plume should impart significant fine structure to the solar wind that may be detectable by Parker Solar Probe and Solar Orbiter.

V M Uritsky↗

Inferring small-sample-reactivity worth from oscillation experiment with auto-corrected local flux

To estimate reactivity changes as a function of time from the signal of detectors count rates, one usually uses point kinetic equations. However, when detectors are located inside the core, it is often necessary to compute a correction factor to account for the flux-shape deformation (i.e. the change of detector efficiency that may have been induced by neutron absorbers movements). Various simulation methods have been developed depending on transients observed. They range from the resolution of static equations up to full-time-dependent neutron transport in 3D. However, those methods require high computational-cost calculations of local reaction rates that rely on nuclear data evaluations. Therefore, a 'measured' reactivity is never determined through pure measurements. Taking advantage of the known oscillation movement of a small sample at the center of the MINERVE reactor, the innovative algorithm presented in this paper derives the correction factor needed to get no time drift on reactivity steps. Derived correction factors are in good agreement with those obtained independently with TRIPOLI-4 through sensitivity analysis. Low sensitivity to the sample composition and moderate sensitivity to the detectors position reinforce the reliability of those results. This application paves the way to almost pure nuclear data reactivity measurements. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Unsupervised learning-enabled pulsed infrared thermographic microscopy of subsurface defects in stainless steel

Metallic structures produced with laser powder bed fusion (LPBF) additive manufacturing method (AM) frequently contain microscopic porosity defects, with typical approximate size distribution from one to 100 microns. Presence of such defects could lead to premature failure of the structure. In principle, structural integrity assessment of LPBF metals can be accomplished with nondestructive evaluation (NDE). Pulsed infrared thermography (PIT) is a non-contact, one-sided NDE method that allows for imaging of internal defects in arbitrary size and shape metallic structures using heat transfer. PIT imaging is performed using compact instrumentation consisting of a flash lamp for deposition of a heat pulse, and a fast frame infrared (IR) camera for measuring surface temperature transients. However, limitations of imaging resolution with PIT include blurring due to heat diffusion, sensitivity limit of the IR camera. We demonstrate enhancement of PIT imaging capability with unsupervised learning (UL), which enables PIT microscopy of subsurface defects in high strength corrosion resistant stainless steel 316 alloy. PIT images were processed with UL spatial–temporal separation-based clustering segmentation (STSCS) algorithm, refined by morphology image processing methods to enhance visibility of defects. The STSCS algorithm starts with wavelet decomposition to spatially de-noise thermograms, followed by UL principal component analysis (PCA), fine-tuning optimization, and neural learning-based independent component analysis (ICA) algorithms to temporally compress de-noised thermograms. The compressed thermograms were further processed with UL-based graph thresholding K-means clustering algorithm for defects segmentation. The STSCS algorithm also includes online learning feature for efficient re-training of the model with new data. For this study, metallic specimens with calibrated microscopic flat bottom hole defects, with diameters in the range from 203 to 76 µm, were produced using electro discharge machining (EDM) drilling. While the raw thermograms do not show any material defects, using STSCS algorithm to process PIT images reveals defects as small as 101 µm in diameter. To the best of our knowledge, this is the smallest reported size of a sub-surface defect in a metal imaged with PIT, which demonstrates the PIT capability of detecting defects in the size range relevant to quality control requirements of LPBF-printed high-strength metals.

36 MATERIALS SCIENCE↗

Impedance-Based Prediction of SSR-Generated Harmonics in Doubly-Fed Induction Generators

Doubly-fed induction generators (DFIGs) are prone to subsynchronous resonance (SSR) with series-compensated transmission lines and experiencing SSR-generated harmonics of wide-ranging magnitudes. This paper presents a large-signal impedance-based approach for the prediction of the magnitude of SSR-generated harmonics in DFIGs. Note that the large-signal impedance of a network shows variation in its impedance response with the magnitude of perturbation at its terminals. It is discovered in this paper that the impedance of a series-compen- sated transmission line starts changing with the magnitude of SSR because of the metal-oxide varistor-based protection of series capacitors in the line. The changing line impedance halts the growth of the SSR beyond a certain point, where the DFIG forms a limit-cycle mode and emits harmonics at the SSR frequency. The discovered role of the protection of series capacitors in SSR not only enables the prediction of SSR-generated harmonics but also can be used for the timely detection of an SSR event to avoid system disruption and sudden loss of bulk generation. The large- signal impedance-based prediction of SSR-generated harmonics in DFIG-based wind and hydro generators is demonstrated using numerical simulations of a 2.5-MVA, 690-V DFIG supplying to a series-compensated transmission line.

24 POWER TRANSMISSION AND DISTRIBUTION↗

HOLISMOKES IX. Neural network inference of strong-lens parameters and uncertainties from ground-based images

Modeling of strong gravitational lenses is a necessity for further applications in astrophysics and cosmology. With the large number of detections in current and upcoming surveys, such as the Rubin Legacy Survey of Space and Time (LSST), it is pertinent to investigate automated and fast analysis techniques beyond the traditional and time-consuming Markov chain Monte Carlo sampling methods. Building upon our (simple) convolutional neural network (CNN), we present here another CNN, specifically a residual neural network (ResNet), that predicts the five mass parameters of a singular isothermal ellipsoid (SIE) profile (lens center x and y, ellipticity e x and e y , Einstein radius θ E ) and the external shear (γ ext, 1 , γ ext, 2 ) from ground-based imaging data. In contrast to our previous CNN, this ResNet further predicts the 1σ uncertainty for each parameter. To train our network, we use our improved pipeline to simulate lens images using real images of galaxies from the Hyper Suprime-Cam Survey (HSC) and from the Hubble Ultra Deep Field as lens galaxies and background sources, respectively. We find very good recoveries overall for the SIE parameters, especially for the lens center in comparison to our previous CNN, while significant differences remain in predicting the external shear. From our multiple tests, it appears that most likely the low ground-based image resolution is the limiting factor in predicting the external shear. Given the run time of milli-seconds per system, our network is perfectly suited to quickly predict the next appearing image and time delays of lensed transients. Therefore, we use the network-predicted mass model to estimate these quantities and compare to those values obtained from our simulations. Unfortunately, the achieved precision allows only a first-order estimate of time delays on real lens systems and requires further refinement through follow-up modeling. Nonetheless, our ResNet is able to predict the SIE and shear parameter values in fractions of a second on a single CPU, meaning that we are able to efficiently process the huge amount of galaxy-scale lenses expected in the near future.

79 ASTRONOMY AND ASTROPHYSICS↗

High-energy emission from a magnetar giant flare in the Sculptor galaxy

Magnetars are the most highly magnetized neutron stars in the cosmos (with magnetic field 10 13 –10 15 G). Giant flares from magnetars are rare, short-duration (about 0.1 s) bursts of hard X-rays and soft γ rays. Owing to the limited sensitivity and energy coverage of previous telescopes, no magnetar giant flare has been detected at gigaelectronvolt (GeV) energies. Here, we report the discovery of GeV emission from a magnetar giant flare on 15 April 2020 (refs. 3,4 and A. J. Castro-Tirado et al., manuscript in preparation). The Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope detected GeV γ rays from 19 s until 284 s after the initial detection of a signal in the megaelectronvolt (MeV) band. Our analysis shows that these γ rays are spatially associated with the nearby (3.5 megaparsecs) Sculptor galaxy and are unlikely to originate from a cosmological γ-ray burst. Thus, we infer that the γ rays originated with the magnetar giant flare in Sculptor. Furthermore, we suggest that the GeV signal is generated by an ultra-relativistic outflow that first radiates the prompt MeV-band photons, and then deposits its energy far from the stellar magnetosphere. After a propagation delay, the outflow interacts with environmental gas and produces shock waves that accelerate electrons to very high energies; these electrons then emit GeV γ rays as optically thin synchrotron radiation. This observation implies that a relativistic outflow is associated with the magnetar giant flare, and suggests the possibility that magnetars can power some short γ-ray bursts.

79 ASTRONOMY AND ASTROPHYSICS↗

Direct analysis of cotton swipes for uranium and plutonium isotopic determination by microextraction-ICP-MS

The determination of uranium and plutonium isotopic abundance on environmental samples collected by International Atomic Energy Agency (IAEA) inspectors is vital for the detection of undeclared nuclear activities and material under the international nuclear safeguards regime. Current analytical protocols require time-consuming sample preparation steps prior to subsequent measurement by inorganic mass spectrometry (MS). Recent efforts from this laboratory have focused on developing sample preparation methods for faster analysis, potentially allowing higher sample throughput[1]. Alternative methods including microextraction sampling in conjunction with inductively coupled plasma-mass spectrometry (ICP-MS) have been recently explored. This methodology, microextraction-ICP-MS, was developed such that uranium and plutonium could be extracted from the swipe surface and directed into the ICP-MS for an in-situ measurement, eliminating the need for swipe ashing and digestion. A commercial off-the-shelf microextraction system was customized with an automated movable XY stage that can be programmed to save sampling locations, allowing for automated rapid sampling of swipe surfaces. Additional efforts have focused on the utilization of collision cell technology to the microextraction ICP-MS method. This would eliminate the need for lengthy column chemistry procedures to purify separated uranium and plutonium fractions before analysis. Here, the extracted U/Pu analyte is measured by reacting the uranium ions with CO2 in the collision cell of an ICP-MS, shifting the uranium to UO+, which will not interfere with the plutonium isotopic determination. The developed method utilizing collision cell – ICP-MS technology has demonstrated the ability to measure plutonium isotope ratios in the presence of high uranium concentration on the transient signal from the microextraction system utilizing certified reference materials from JRC-Geel and the New Brunswick Laboratory Program Office.

Bradley, Veronica↗

Gravitational-Wave Physics and Astronomy in the 2020s and 2030s

The 100 years since the publication of Albert Einstein’s theory of General Relativity saw significant development of the understanding of the theory, the identification of potential astrophysical sources of sufficiently strong gravitational waves, and development of key technologies for gravitational-wave detectors. In 2015, the first gravitational wave signals were detected by the two US Advanced LIGO instruments. In 2017, Advanced LIGO and the European Advanced Virgo detectors pinpointed a binary neutron star coalescence that was also seen across the electromagnetic spectrum. The field of gravitational wave astronomy is just starting, and this Roadmap of future developments surveys the potential for growth in bandwidth and sensitivity of future gravitational-wave detectors and discusses the science results anticipated to come from upcoming instruments.

Compact astrophysical objects↗

Chromospheric evaporation and decimetric radio emission in solar flares

We have discovered decimetric signatures of the chromospheric evaporation process. Evidence for the radio detection of chromospheric evaporation is based on the radio-inferred values of (1) the electron density, (2) the propagation speed, and (3) the timing, which are found to be in good agreement with statistical values inferred from the blueshifted Ca XIX soft X-ray line. The physical basis of our model is that free-free absorption of plasma emission is strongly modified by the steep density gradient and the large temperature increase in the upflowing flare plasma. The steplike density increase at the chromospheric evaporation front causes a local discontinuity in the plasma frequency, manifested as almost infinite drift rate in decimetric type III bursts. The large temperature increase of the upflowing plasma considerably reduces the local free-free opacity (due to the T(exp -3/2) dependence) and thus enhances the brightness of radio bursts emitted at the local plasma frequency near the chromospheric evaporation front, while a high-frequency cutoff is expected in the high-density regions behind the front, which can be used to infer the velocity of the upflowing plasma. From model calculations we find strong evidence that decimetric bursts with a slowly drifting high-frequency cutoff are produced by fundamental plasma emission, contrary to the widespread belief that decimetric bursts are preferentially emitted at the harmonic plasma level. We analyze 21 flare episodes from 1991-1993 for which broadband (100-3000 MHz) radio dynamic spectra from Pheonix, hard X-ray data from (BATSE/CGRO) and soft X-ray data from Burst and Transient Source Experiment/Compton Gamma Ray Observatory (GOES) were available.

Aschwanden, Markus J.↗

Resolving femtosecond photoinduced energy flow: capture of nonadiabatic reaction pathway topography and wavepacket dynamics from photoexcitation through the conical intersection seam (Final Technical Report)

The dynamics that take place within just tens to hundreds of femtoseconds following the absorption of light by a molecule can play a critical role in how the absorbed energy is directed, allowing it to be used for a specific function or dissipated harmlessly. The form of chemical change that occurs rapidly in these molecules is called a “nonadiabatic electronic transition.” Such transitions are known to mediate energy flow in natural biological systems such as the ultraviolet photoprotection mechanism of DNA and the first step of the human vision response. Understanding how these mechanisms work precisely may help scientists achieve controlled manipulation of solar energy or optical control of a wide range of energy management functions in artificial systems. Experimental methods, however, have not yet allowed a precisely resolved and complete measurement of nonadiabatic electronic transitions. This constitutes a major obstacle to progress in the field. For progress to occur that would inform a wide body of research aiming to efficiently harness the energy of light for practical purposes, it is especially important to benchmark computational models of the molecules undergoing these rapid changes with experimental measurements, in order to learn which models are accurate. With Dept. of Energy funding, we have made strong progress towards establishing a new optical method for experimentally detecting the full nonadiabatic electronic transition. This requires having coordinated pulses of light covering the visible through the mid-infrared range of the electromagnetic spectrum that last only ten femtoseconds. We have developed a new, relatively simple approach for generating such pulses of laser light, and have incorporated them into a time-resolved spectrometer for measuring rapid changes in molecules. These tools can provide the greater precision and new types of data that are needed to benchmark computational models of molecular change and thus to make progress in the field. Our tools were tested on graphene, an excellent solid-state sample for verifying the capabilities and limitations of our instrumentation. The investment made in these tools by the Dept. of Energy Office of Science will allow new fundamental scientific understanding of energy dynamics in molecules in future studies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Modelling populations of kilonovae

Abstract The 2017 detection of a kilonova coincident with gravitational-wave emission has identified neutron star mergers as the major source of the heaviest elements and dramatically constrained alternative theories of gravity. Observing a population of such sources has the potential to transform cosmology, nuclear physics, and astrophysics. However, with only one confident multi-messenger detection currently available, modelling the diversity of signals expected from such a population requires improved theoretical understanding. In particular, models that are quick to evaluate and are calibrated with more detailed multi-physics simulations are needed to design observational strategies for kilonovae detection and to obtain rapid-response interpretations of new observations. We use grey-opacity models to construct populations of kilonovae, spanning ejecta parameters predicted by numerical simulations. Our modelling focuses on wavelengths relevant for upcoming optical surveys, such as the Rubin Observatory Legacy Survey of Space and Time (LSST). In these simulations, we implement heating rates that are based on nuclear reaction network calculations. We create a Gaussian-process emulator for kilonova grey opacities, calibrated with detailed radiative transfer simulations. Using recent fits to numerical relativity simulations, we predict how the ejecta parameters from binary neutron star (BNS) mergers shape the population of kilonovae, accounting for the viewing-angle dependence. Our simulated population of BNS mergers produce peak i-band absolute magnitudes of −20 ≤ Mi ≤ −11. A comparison with detailed radiative transfer calculations indicates that further improvements are needed to accurately reproduce spectral shapes over the full light curve evolution.

79 ASTRONOMY AND ASTROPHYSICS↗

Gamma-ray burst constraints on the galactic frequency of extrasolar Oort Clouds

With the strong Compton Gamma-Ray Observatory/Burst and Transient Source Experiment (CGRO/BATSE) evidence that most gamma-ray bursts do not come from galactic neutron stars, models involving the accretion of a comet onto a neutron star (NS) no longer appear to be strong contenders for explaining the majority of bursts. If this is the case, then it is worth asking whether the lack of an observed galactic gamma-ray burst population provides a useful constraint on the number of comets and comet clouds in the galaxy. Owing to the previously unrecognized structural weakness of cometary nuclei, we find the capture cross sections for comet-NS events to be much higher than previously published estimates, with tidal breakup at distances R(sub b) approx. equals 4 x 10(exp 10) cm from the NS. As a result, impacts of comets onto field NSs penetrating the Oort Clouds of other stars are found to dominate all other galactic NS-comet capture rates by a factor of 100. This in turn predicts that if comet clouds are common, there should be a significant population of repeater sources with (1) a galactic distribution, (2) space-correlated repetition, and (3) a wide range of peak luminosities and luminosity time histories. If all main sequence stars have Oort Clouds like our own, we predict approximately 4000 such repeater sources in the Milky Way at any time, each repeating on time scales of months to years. Based on estimates of the sensitivity of the CGRO/BATSE instrument and assuming isotropic gamma-ray beaming from such events, we estimate that a population of approximately 20-200 of these galactic NS-Oort Cloud gamma-ray repeater sources should be detectable by CGRO. In addition, if giant planet formation is common in the galaxy, we estimate that the accretion of isolated comets injected to the interstellar medium by giant planet formation should produce an additional source of galactic, nonrepeating, events. Comparing these estimates to the 3-4 soft gamma-ray repeater sources detected by BATSE, one is forced to conclude that (1) comet impacts on NSs are inefficient at producing gamma rays; or (2) the gamma rays from such events are highly beamed; or (3) the fraction of stars in the galaxy with Oort Clouds like our own is not higher than a few percent.

Shull, J. Michael↗

Legacy Survey of Space and Time Data Preview 2: dia_object dataset type

We present Rubin Data Preview 2 (DP2), the second data preview from the NDF-DOE Vera C. Rubin Observatory. Data Preview 2 (DP2) comprises coadds, detection catalogs, and ancillary data products; and when fully released will also include single-epoch images and difference images. DP2 is derived from observations acquired by the LSST Science Camera (LSSTCam) on the Simonyi Survey Telescope at the Summit Facility on Cerro Pachón, Chile, primarily during the on-sky commissioning campaign between 2025-04-16 and 2025-09-21, supplemented by observations taken between 2025-10-25 and 2026-01-06 that overlap the commissioning footprint. The DP2 footprint comprises the Science Validation wide-area survey, five Deep Drilling Fields, and a number of targeted small-field regions, including Trifid and Lagoon, Prawn, M49, and New Horizons, all observed as part of the Rubin First Look campaign. Each field was imaged in up to six broad photometric bands, ugrizy, and coadded to produce deep imaging covering an estimated 3,000 deg2. The addition of single-visit-only areas expands the total DP2 footprint to an estimated 15,000 deg2, with coverage in at least one filter. The median per-visit PSF FWHM across the wide-area survey ranges from 1.17 arcsec in the z band to 1.26 arcsec in g and r bands. The deepest field, reaches estimated coadded 5σ depths of u=26 mag, g=26.8 mag, r=26.3 mag, i=26.1 mag, z=25.3 mag, y=23.9 mag. Based on a roughly five-month primary observing baseline and covering only part of the eventual LSST footprint, DP2's area, depth, and multiband coverage nonetheless support a broad range of early science investigations ahead of LSST Data Release This dataset is a subset of the full data release consisting of the dia_object dataset type. These are derived properties for transient and variable objects. This release contains 2,112 datasets of this type.

79 ASTRONOMY AND ASTROPHYSICS↗

The dark energy survey: detection of weak lensing magnification of supernovae and constraints on dark matter haloes

The residuals of the distance moduli of Type Ia supernovae (SNe Ia) relative to a Hubble diagram fit contain information about the inhomogeneity of the Universe, due to weak lensing magnification by foreground matter. By correlating the residuals of the Dark Energy Survey Year 5 SN Ia sample (DES-SN5YR) with extragalactic foregrounds from the DES Y3 Gold catalogue, we detect the presence of lensing at |$6.0 \sigma$| significance. This is the first detection with a significance level above |$5\sigma$|⁠. Constraints on the effective mass-to-light ratios and radial profiles of dark matter haloes surrounding individual galaxies are also obtained. We show that the scatter of SNe Ia around the Hubble diagram is reduced by modifying the standardization of the distance moduli to include an easily calculable de-lensing (i.e. environmental) term. We use the de-lensed distance moduli to recompute cosmological parameters derived from SN Ia, finding in Flat wcold dark matter a difference of |$\Delta \Omega _{\rm M} = +0.036$| and |$\Delta w = -0.056$| compared to the unmodified distance moduli, a change of |$\sim 0.3\sigma$|⁠. We argue that our modelling of SN Ia lensing will lower systematics on future surveys with higher statistical power. We use the observed dispersion of lensing in DES-SN5YR to constrain |$\sigma _8$|⁠, but caution that the fit is sensitive to uncertainties at small scales. Nevertheless, our detection of SN Ia lensing opens a new pathway to study matter inhomogeneity that complements galaxy–galaxy lensing surveys and has unrelated systematics.

79 ASTRONOMY AND ASTROPHYSICS↗