Search NASA⌕ Search

SEARCH · Search NASA

Results for “matched field processing”

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 775 records · Page 43

Assimilation of Wave Imaging Radar Observations for Real-time Wave-by-Wave Forecasting

This project addressed Topic 3: “Wave Measurement Instrumentation for Feed Forward Controls” under the FOA number DE-FOA-0000971. The overall goal of the program was to develop a phase-resolving wave forecasting technique for application to the active control of Wave Energy Conversion (WEC) devices. We have developed an approach that couples a wave imaging marine radar with a phase-resolving linear wave model for real-time wave field reconstruction and forward propagation of the wave field in space and time. The scope of the project was to develop and assess the performance of this novel forecasting system. Specific project goals were as follows: Develop and verify a fast, GPU-based (Graphical Processing Unit) wave propagation model suitable for phase-resolved computation of nearshore wave transformation over variable bathymetry; Compare the accuracy and speed of performance of the wave model against a deep water model in their ability to predict wave field transformation in the intermediate water depths (50 to 70 m) typical of planned WEC sites; Develop and implement a variational assimilation algorithm that can ingest wave imaging radar observations and estimate the time-varying wave conditions offshore of the domain of interest such that the observed wave field is best reconstructed throughout the domain and then use this to produce model forecasts for a given WEC location; Collect wave-resolving marine radar data, along with relevant in situ wave data, at a suitable wave energy test site, apply the algorithm to the field data, assess performance, and identify any necessary improvements; and Develop a production cost estimate that addresses the affordability of the wave forecasting technology and include in the Final Report. The developed forecasting algorithm (“Wavecast”) was evaluated for both speed and accuracy against a substantial synthetic dataset. Early in the project, performance tests definitively demonstrated that the system was capable of forecasting in real-time, as the GPU-based wave model backbone was very computationally efficient. The data assimilation algorithm was developed on a polar grid domain in order to match the sampling characteristics of the observation system (wave imaging marine radar). For verification purposes, a substantial set of synthetic wave data (i.e. forward runs of the wave model) were generated to be used as ground truth for comparison to the reconstructions and forecasts produced by Wavecast. For these synthetic cases, Wavecast demonstrated very good accuracy, for example, typical forecast correlation coefficients were between 0.84-0.95 when compared to the input data. Dependencies on shadowing, observational noise, and forecast horizon were also identified. During the second year of the project, a short field deployment was conducted in order to assess forecast accuracy under field conditions. For this, a radar was installed on a fishing vessel and observations were collected at the South Energy Test Site (SETS) off the coast of Newport, OR. At the SETS site, simultaneous in situ wave observations were also available owing to an ongoing field project funded separately. Unfortunately, the position and heading information that was available for the fishing vessel were not of sufficient accuracy in order to validate the forecast in a phase-resolving sense. Instead, a spectral comparison was made between the Wavecast forecast and the data from the in situ wave buoy. Although the wave and wind conditions during the field test were complex, the comparison showed a promising reconstruction of the wave spectral shape, where both peaks in the bimodal spectrum were represented. However, the total reconstructed spectral energy (across all directions and frequencies) was limited to 44% of the observed spectrum. Overall, wave-by-wave forecasting using a data assimilation approach based on wave imaging radar observations and a physics-based wave model shows promise for short-term phase-resolved predictions. Two recommendations for future work are as follows: first, we would recommend additional focused field campaigns for algorithm validation. The field campaign should be long enough to capture a range of wave conditions relevant to the target application and WEC site. In addition, it will be crucial to make sure the vessel of choice has high accuracy position and heading instrumentation (this instrumentation is commercially available but not standard on commercial fishing vessels). The second recommendation is to expand the model physics in the wave model backbone to include some nonlinear effects. Specifically, the third-order correction to the wave speed due to amplitude dispersion would be the next step in order to more accurately represent the phase speeds of large amplitude waves.

16 TIDAL AND WAVE POWER↗

Combined measurement and modeling of the hydrological impact of hydraulic redistribution using CLM4.5 at eight AmeriFlux sites

Effects of hydraulic redistribution (HR) on hydrological, biogeochemical, and ecological processes have been demonstrated in the field, but the current generation of standard earth system models does not include a representation of HR. Though recent studies have examined the effect of incorporating HR into land surface models, few (if any) have done cross-site comparisons for contrasting climate regimes and multiple vegetation types via the integration of measurement and modeling. Here, we incorporated the HR scheme of Ryel et al. (2002) into the NCAR Community Land Model Version 4.5 (CLM4.5), and examined the ability of the resulting hybrid model to capture the magnitude of HR flux and/or soil moisture dynamics from which HR can be directly inferred, to assess the impact of HR on land surface water and energy budgets, and to explore how the impact may depend on climate regimes and vegetation conditions. Eight AmeriFlux sites with contrasting climate regimes and multiple vegetation types were studied, including the Wind River Crane site in Washington State, the Santa Rita Mesquite savanna site in southern Arizona, and six sites along the Southern California Climate Gradient. HR flux, evapotranspiration (ET), and soil moisture were properly simulated in the present study, even in the face of various uncertainties. Our cross-ecosystem comparison showed that the timing, magnitude, and direction (upward or downward) of HR vary across ecosystems, and incorporation of HR into CLM4.5 improved the model-measurement matches of evapotranspiration, Bowen ratio, and soil moisture particularly during dry seasons. Lastly, our results also reveal that HR has important hydrological impact in ecosystems that have a pronounced dry season but are not overall so dry that sparse vegetation and very low soil moisture limit HR.

58 GEOSCIENCES↗

NGEE Arctic Meteorological Data from Instrumented PCDC Stations, Council Site, Seward Peninsula, Alaska, 2018 - 2023

This dataset provides meteorological data collected at two instrumented Power, Communication, and Data Collection (PCDC) Stations, located at Council Road Mile Marker 71 (CN_MM71) and Council Road Mile Marker 68 (Rocky Top). Measurements began in June of 2018 at CN_MM71 and in September 2018 at Rocky Top and ended in 2023. Raw data included. Meteorological data include air temperature, relative humidity, wind speed and direction, and snow depth. Measurements have been aggregated to both 10-minute and 30-minute reporting intervals provided in four comma-separated (*.csv) data files.• The 30-minute data files include only the primary meteorological variables, leaf wetness, and snow depth.. The 30-minute reporting interval matches other coincident site data collections.• The 10-minute files provide the primary meteorological variables and numerous system operational parameters such as battery charge and solar and wind turbine status. This 10-minute data product is valuable for analyzing solar and wind charging system performance in arctic landscapes.Daily photographic images (09:00, 12:00, and 15:00 AKST/AKDT) of the PCDC Stations and the surrounding landscape that are captured by an adjacent web camera, are also provided.The daily 12:00 images from CN_MM71 camera have been compiled into time-lapse movies:• October 2018 – June 12, 2019 = Council_winter_2018_2019.mov• January 2019 – December 2019 = Council_2019.mp4The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research.The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska.Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗

SU-E-T-323: The FLUKA Monte Carlo Code in Ion Beam Therapy

Purpose: Monte Carlo (MC) codes are increasingly used in the ion beam therapy community due to their detailed description of radiation transport and interaction with matter. The suitability of a MC code demands accurate and reliable physical models for the transport and the interaction of all components of the mixed radiation field. This contribution will address an overview of the recent developments in the FLUKA code oriented to its application in ion beam therapy. Methods: FLUKA is a general purpose MC code which allows the calculations of particle transport and interactions with matter, covering an extended range of applications. The user can manage the code through a graphic interface (FLAIR) developed using the Python programming language. Results: This contribution will present recent refinements in the description of the ionization processes and comparisons between FLUKA results and experimental data of ion beam therapy facilities. Moreover, several validations of the largely improved FLUKA nuclear models for imaging application to treatment monitoring will be shown. The complex calculation of prompt gamma ray emission compares favorably with experimental data and can be considered adequate for the intended applications. New features in the modeling of proton induced nuclear interactions also provide reliable cross section predictions for the production of radionuclides. Of great interest for the community are the developments introduced in FLAIR. The most recent efforts concern the capability of importing computed-tomography images in order to build automatically patient geometries and the implementation of different types of existing positron-emission-tomography scanner devices for imaging applications. Conclusion: The FLUA code has been already chosen as reference MC code in many ion beam therapy centers, and is being continuously improved in order to match the needs of ion beam therapy applications. Parts of this work have been supported by the European FP7 project ENVISION (grant agreement no. 241851)

07 ISOTOPE AND RADIATION SOURCES↗

Apolipoprotein E expression and behavioral toxicity of high charge, high energy (HZE) particle radiation

Apolipoprotein E (apoE) is a lipid binding protein that plays an important role in tissue repair following brain injury. In the present studies, we have investigated whether apoE affects the behavioral toxicity of high charge, high energy (HZE) particle radiation. METHODS: Sixteen male apoE knockout (KO) mice and sixteen genetically matched wild-type (WT) C57BL mice were used in this experiment. Half of the KO and half of the WT animals were irradiated with 600 MeV/amu iron particles (2 Gy whole body). The effect of irradiation on motor coordination and stamina (Rotarod test), exploratory behavior (open field test), and spatial working and reference memory (Morris water maze) was assessed. ROTAROD TEST: Performance was adversely affected by radiation exposure in both KO and WT groups at 30 d after irradiation. By 60 d after radiation, the radiation effect was lost in WT, but still apparent in irradiated KO mice. OPEN FIELD TEST: Radiation reduced open field exploratory activity 14, 28, 56, 84, and 168 d after irradiation of KO mice, but had no effect on WT mice. MORRIS WATER MAZE: Radiation adversely affected spatial working memory in the KO mice, but had no discernible effect in the WT mice as assessed 180 d after irradiation. In contrast, irradiated WT mice showed marked impairment of spatial reference memory in comparison to non-irradiated mice, while no effect of radiation was observed in KO mice. CONCLUSIONS: These studies show that apoE expression influences the behavioral toxicity of HZE particle radiation and suggest that apoE plays a role in the repair/recovery from radiation injury of the CNS. ApoE deficiency may exacerbate the previously reported effects of HZE particle radiation in accelerating the brain aging process.

Non-NASA Center↗

Heavy-meson decay constants from QCD sum rules

We sketch a recent sum-rule extraction of the decay constants of the heavy pseudoscalar mesons D, D{sub s}, B, and B{sub s} from the two-point correlator of heavy-light pseudoscalar currents. Our main emphasis lies on the control over all the uncertainties in the decay constants, related both to the input QCD parameters and to the limited accuracy of the method of sum rules. Gaining this control has become possible by application of our new procedure of extracting hadron observables based on a dual threshold depending on the Borel parameter. For the charmed-meson decay constants, we find fD = (206.2{+-}7.3{sub (OPE)}{+-}5.1{sub (syst)}) MeV, fD{sub s} = (245.3{+-}15.7{sub (OPE)}{+-}4.5{sub (syst)}) MeV. For the beauty mesons, the decay constants turn out to be extremely sensitive to the precise value of the {ovr MS} mass of the b-quark, {bar m}{sub b}({bar m}{sub b}). By requiring our sum-rule estimate to match the average of the lattice determinations of f{sub B}, we extract the rather accurate value {bar m}{sub b}({bar m}{sub b}) = (4.245{+-}0.025) GeV. Feeding this parameter value into our sum-rule formalism leads to the beauty-meson decay constants fB = (193.4{+-}12.3{sub (OPE)}{+-}4.3{sub (syst)}) MeV, fB{sub s} = (232.5{+-}18.6{sub (OPE)}{+-}2.4{sub (syst)}) MeV.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Investigating Low-Altitude Constellations of Ad-Hoc Lunar PNT System for Distributed Spacecraft Autonomy

In this study, we examine a low-altitude Lunar Position, Navigation, and Timing (LPNT) constellations and the localization performance of Centralized Extended Kalman Filter (CEKF) and Decentralized Extended Kalman Filter (DEKF) algorithms. The primary investigation involves a 100-node swarm operating at a 100 km altitude, in contrast to previous studies that examined a 21-node asset in a frozen-orbit at 5,500 km. The autonomous operation of large-scale swarm is based on two-way Inter-Satellite Link (ISL) measurements, which involve pseudoranges and relative velocities among swarm nodes. We perform a numerical assessment of the two filtering approaches, utilizing ‘fully sampled’ measurements from all available assets as well as ‘two ISL’ measurements where each spacecraft is restricted to only two antennas. This research includes an analysis of CEKF under 2-ISL constraints and evaluates the performance of DEKF in a 100-node swarm, which has not been explored in previous studies. In addition, we examine the impact of increasing the sampling frequency for DEKF, showing that the update cycle can be shortened from a 10-minute interval. A novel approach for ‘2-ISL limited’ DEKF will also be introduced, using a matching formulation that exhaustively enumerates all potential matches. This study provides valuable insights into large-scale distributed swarm operations, considering various filter configurations, sampling frequencies, matching strategies, and scalability of CEKF and DEKF for low-altitude LPNT applications. The Lunar PNT technology plays a key role in providing reliable and robust navigation services on the Moon's surface and the South pole, where the primary Lunar missions are planned. To support upcoming Lunar missions, including small satellites from NASA's Commercial Lunar Payload Services program, the Lunar PNT system must be adaptable to smaller platforms like CubeSats. Driven by the growing involvement of public and private exploration partnerships, the traditional low Earth orbit missions are shifting to beyond geosynchronous orbit [1]. These upcoming missions aim to foster a sustainable and innovative exploration program, in collaboration with commercial and international partners, to facilitate human expansion throughout the solar system and return new knowledge and opportunities to Earth [2]. As part of this trend, there are increasing efforts to utilize science missions in Lunar orbit to develop a non-dedicated and ad-hoc PNT network system. Two traditional approaches, the Deep Space Network (DSN) and the weak signal Global Positioning System (GPS), are established deep-space navigation technologies for missions beyond the geosynchronous orbit. Beginning in 1958, the DSN was developed to communicate with the Explorer 1 spacecraft based on the use of radiometric tracking in spacecraft navigation [3]. The DSN is capable of providing nearly unfettered coverage to spacecraft beyond low-Earth orbit (LEO), however, increased space mission volume has created concerns about future expectations of DSN usage for spacecraft navigation [4]. For cislunar mission applications, the position accuracy using DSN achieves 100 m (3σ) with at least three geometrically diverse ground stations when using radiometric tracking alone [5]. The DSN's dependence on Earth-based ground stations restricts its operational capabilities to periods of Earth visibility. This limitation, coupled with its poor localization performance, renders the DSN unsuitable for future lunar missions that demand continuous tracking and precise positioning. To satisfy the increasing requirements of DSN in Lunar applications, spacecrafts are also required to improve their onboard antenna power and efficiency of the transmission. However, there is an important aggregate cost trade between adding capabilities to every spacecraft and adding to a capacity on the ground that serves multiple spacecraft [6]. A weak GPS system can provide PNT service while the user spacecraft is bound to the Moon, leveraging a single, steerable high gain antenna with the relatively narrow beam which includes all the sources in its field of view [7]. However, the higher the altitude the receiver is above the GPS constellations, the poorer and the weaker are the relative geometry and the received signal powers, respectively, leading to a significant navigation accuracy reduction [8]. The transmitted power becomes weaker with increasing distance from the Earth as well as signals tracked from one of the side lobes of the GPS antenna pattern. As a results, the number of visible satellites and relative geometric condition of the GPS satellites at very high altitude drops dramatically and reduces the navigation solution accuracy. Therefore, the weak GPS system is also not an ideal way to provide PNT service to upcoming Lunar missions when considering its limited geometric condition and the recued navigation accuracy. Another navigation approach on the Moon is being developed, similar to the Global Navigation Satellite System (GNSS) on Earth, aiming to offer navigation service with continuous 24/7 coverage across the entire Lunar surface. For example, lunar communications relay and navigation systems (LCRNS) by NASA and Lunar navigation satellite systems (LNSS) by JAXA are designed to serve as dedicated Position, Navigation, and Timing (PNT) systems for the Moon. However, designing a dedicated LNSS and PNT service involves additional challenges, which are unique to the lunar environment, including limited payload capacity for the CubeSat platform, i.e., the size, weight, and power (SWaP) of the onboard clock, limited lunar ground monitoring stations, and limited financial investment as compared to the legacy Earth-GPS [9]. NASA’s focus on utilizing CubeSat platforms on the Moon leads to an alternative Lunar navigation platform that leverages the existing Lunar science and exploration assets. The small satellites used in Lunar missions can be used to create a low-cost, autonomous, ad-hoc, and on-demand mission-centric Lunar PNT swarm capable of providing PNT services to these low-cost lunar missions [10]. As upcoming Lunar missions will often operate at low-altitude about 30 km to 100 km for scientific observations and mapping purposes, the low-altitude orbital constellations could be employed to create an ad-hoc Lunar PNT system. However, several issues must be addressed, such as the instability of these orbits, which often require maintenance or are only suitable for short-duration missions, operating for fewer than 90 days. Additionally, at an altitude of 100 km, the satellites have a limited period during which they are above the horizon and capable of providing PNT service to users. The implementation of a non-dedicated, ad-hoc Lunar navigation constellation facilitates on-demand PNT services. A preliminary study of ad-hoc Lunar PNT system was conducted using 21 spacecraft in 5,5000 km altitude frozen orbits to test its feasibility and a basic performance of orbital asset localization among ad-hoc Lunar constellations in small satellites format [10]. These swarm assets are designed for autonomous localization with minimal Earth interaction, reducing dependency on bandwidth and ground resources. The design in [10] demonstrated the feasibility of a decentralized PNT approach, specifically employing a DEKF approach for state estimation, which helps minimize onboard operating costs. The DEKF method distributes computation across individual satellites, which lightens the computational load while maintaining accuracy in orbit ephemeris and clock offsets, similar to centralized systems [11]. In a follow-on study [12], each spacecraft was limited to 2 communications antennae, forcing the selection of measurements and scheduling spacecraft activities to perform the measurements. A matching algorithm is implemented to select the best measurements and schedule position estimation updates. The decentralized localization performance is also investigated with increasing levels of network degradation for swarm assets considering the impact of intermittent and permanent communication failure, to demonstrate the robustness and fidelity of the decentralized Lunar PNT service [13]. This study confirmed that the ad-hoc PNT constellations in frozen orbit are highly robust and resilient to communication failures. However, unlike frozen orbit swarm assets, the low-altitude satellites have a limited ground view at an altitude of 100 km, where the ad-hoc Lunar constellation consists of 98 low-altitude satellites, evenly distributed across seven circular polar orbital planes, alongside two satellites in a frozen orbit at an altitude of 5,500 km (Figure 1). Therefore, the number of satellites visible to ground users is significantly limited in low-altitude orbit constellations. As each visibility of a spacecraft remains intact for only a few ticks before it moves out of the field of view, the ground user encounters challenges in maintaining continuous navigation service, resulting in sparse availability and provision of Lunar PNT system. Consequently, service availability is primarily restricted to the Lunar South Pole region (Figure 2). Given these limitations and concerns, the localization performance of low-altitude swarm assets will be assessed in this study. We focus on the investigation of the localization performance of low-altitude swarm assets and ground users near the Lunar South Pole. The overall flow of the Lunar PNT simulation incorporates the DEKF approach of asset localization and the weighted least-squares approach in user localization (Figure 3). The autonomous Lunar PNT simulation is primarily implemented in MATLAB, where the DEKF based on the matching scheduler is implemented with Google’s OR-tools as a model builder and Gurobi optimization tool as a backend solver. The General Mission Analysis Tool (GMAT) is utilized to generate ephemeris data for swarm assets, and accounts for satellite orbital details, mass, and perturbations like solar radiation pressure and drag coefficients. Each ephemeris dataset is produced in the Moon International Celestial Reference Frame (ICRF) inertial coordinate system. For state estimation, the distributed swarm assets rely on two-way Inter-Satellite Link (ISL) measurements, which involve tracking pseudoranges and relative velocities between visible satellites and anchor nodes during each observation. Numerical evaluations of the decentralized localization process are conducted to demonstrate the feasibility of the low-altitude PNT system in providing reliable navigation services. The main approach involves using DEKF and CEKF to localize 100 satellites in low-altitude constellations, where the CEKF is implemented to serve as a baseline for comparing the performance of distributed algorithms. In both cases, we evaluate ‘fully sampled’ measurements from all available assets, and ‘two ISL’ measurements when spacecraft are constrained to have only two antennas. We test four estimation techniques: CEKF fully sampled, CEKF two ISL, DEKF fully sampled, and DEKF two ISL filters. As the DEKF update cycle is comprised of network setup, communication, and computations, a global broadcast network and 2-way ISL network setup will take from 4 to 6 minutes as maximum [12]. In this simulation, the DEKF update cycle is set to 10 minutes, including a 4-minute latency for obtaining and computing the actual measurement updates. We experiment an increased update cycle to demonstrate the feasibility and evaluate the impact on localization performance using various tuning values for measurement noise covariances (Figures 4 and 5). By comparing centralized and decentralized approaches using a matching algorithm, we analyze the influence of cross-correlation factors in the covariance matrix, assuming 100% reliability of all assets and measurements. The increased frequency and the adjustments of tuning parameters reveal distinct error patterns between the two scenarios. The localization accuracy of the swarm assets and ground users is assessed by taking the median error across 100 assets and one ground user (84.9°S, 137.5°E) over 7-day simulation period (Table 1). Since the user localization accuracy is significantly affected by the performance of the swarm assets, it is crucial to maintain high localization accuracy within the swarm. This study will continue to explore decentralized filtering for autonomous LPNT operations, with further investigation of an 'iterative' matching approach which enumerates every valid matching pair, planned for the following month.

Yeji Kim↗

INR-TEM: Robust cavity detection in multifocus TEM images via implicit neural representations

When characterizing materials using transmission electron microscopy (TEM) images, detecting and quantifying small features in microstructures, such as cavities, pose significant challenges. Off-the-shelf object detection models, including YOLOv8, show considerable performance degradation, particularly when images vary in resolution and the objects of interest possess a low percentage of the total image region of interest. In this study, we introduce a novel detection pipeline that incorporates an implicit neural representation (INR)-based detection method, INR-TEM, and two-modality imaging (e.g., under-focused and over-focused images typically acquired during materials characterization) to improve object detection performance. The INR-TEM method incorporates a pixel-wise prediction principle inspired by pixel-wise centerness weighting. INR-TEM demonstrates superior robustness to resolution variability, maintaining high detection accuracy even at low image resolutions compared to YOLOv8. To leverage INR-TEM effectively in real-world two-modality characterization applications, we further integrate a two-stage motion correction pipeline designed explicitly for aligning multifocus TEM images. The alignment process, comprising keypoint (based on scale-invariant feature transform, SIFT) and intensity matching, significantly mitigates the adverse effects of perceived motion-induced image degradation during through-focal TEM imaging, directly enhancing INR-TEM’s detection capability over conventional single-focus images. Our integrated INR-TEM cavity detection framework notably improves performance across various cavity sizes, outperforming off-the-shelf YOLOv8 detections that rely on a single image modality.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

General-purpose event generators for LHC physics

We review the physics basis, main features and use of general-purpose Monte Carlo event generators for the simulation of proton-proton collisions at the Large Hadron Collider. Topics included are: the generation of hard-scattering matrix elements for processes of interest, at both leading and next-to-leading QCD perturbative order; their matching to approximate treatments of higher orders based on the showering approximation; the parton and dipole shower formulations; parton distribution functions for event generators; non-perturbative aspects such as soft QCD collisions, the underlying event and diffractive processes; the string and cluster models for hadron formation; the treatment of hadron and tau decays; the inclusion of QED radiation and beyond-Standard-Model processes. We describe the principal features of the Ariadne, Herwig++, Pythia 8 and Sherpa generators, together with the Rivet and Professor validation and tuning tools, and discuss the physics philosophy behind the proper use of these generators and tools. This review is aimed at phenomenologists wishing to understand better how parton-level predictions are translated into hadron-level events as well as experimentalists wanting a deeper insight into the tools available for signal and background simulation at the LHC.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Heavy states in 3d gravity and 2d CFT

We discuss correlators of light fields in heavy states in AdS 3 gravity and holographic 2d CFTs. In the bulk, the propagator of free fields in AdS backgrounds containing a conical defect or a BTZ black hole can be obtained by solving a wave equation, as well as by the method of images. On the boundary, these geometries are sourced by heavy operator insertions, and the propagator is dual to a heavy-light (HHLL) correlator. By matching its expansion in Virasoro blocks to our bulk results, we determine the OPE coefficients of all contributing states in both the s and t channels. In the s channel, these states are excitations of the light field on top of the heavy state, and their OPE coefficients are the amplitudes to create them. The t-channel OPE is dominated by the Virasoro vacuum block, but there is also an infinite family of light two-particle states that contribute to the correlator. The OPE coefficients that couple these states to heavy operators represent their expectation values in heavy backgrounds. We determine them exactly, derive their asymptotic form at large twist, and discuss their behavior near and above the BTZ threshold, where they become thermal one-point functions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

LHC hadronic jet generation using convolutional variational autoencoders with normalizing flows

Abstract In high energy physics, one of the most important processes for collider data analysis is the comparison of collected and simulated data. Nowadays the state-of-the-art for data generation is in the form of Monte Carlo (MC) generators. However, because of the upcoming high-luminosity upgrade of the Large Hadron Collider (LHC), there will not be enough computational power or time to match the amount of needed simulated data using MC methods. An alternative approach under study is the usage of machine learning generative methods to fulfill that task. Since the most common final-state objects of high-energy proton collisions are hadronic jets, which are collections of particles collimated in a given region of space, this work aims to develop a convolutional variational autoencoder (ConVAE) for the generation of particle-based LHC hadronic jets. Given the ConVAE’s limitations, a normalizing flow (NF) network is coupled to it in a two-step training process, which shows improvements on the results for the generated jets. The ConVAE+NF network is capable of generating a jet in 18.30 ± 0.04 μ s , making it one of the fastest methods for this task up to now.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Carbon Storage and Enhanced Oil Recovery in Pennsylvanian Morrow Formation Clastic Reservoirs: Controls on Oil–Brine and Oil–CO2 Relative Permeability from Diagenetic Heterogeneity and Evolving Wettability

The efficiency of carbon utilization and storage within the Pennsylvanian Morrow B sandstone, Farnsworth Unit, Texas, is dependent on three-phase oil, brine, and CO2 flow behavior, as well as spatial distributions of reservoir properties and wettability. We show that end member two-phase flow properties, with binary pairs of oil–brine and oil–CO2, are directly dependent on heterogeneity derived from diagenetic processes, and evolve progressively with exposure to CO2 and changing wettability. Morrow B sandstone lithofacies exhibit a range of diagenetic processes, which produce variations in pore types and structures, quantified at the core plug scale using X-ray micro computed tomography imaging and optical petrography. Permeability and porosity relationships in the reservoir permit the classification of sedimentologic and diagenetic heterogeneity into five distinct hydraulic flow units, with characteristic pore types including: macroporosity with little to no clay filling intergranular pores; microporous authigenic clay-dominated regions in which intergranular porosity is filled with clay; and carbonate–cement dominated regions with little intergranular porosity. Steady-state oil–brine and oil–CO2 co-injection experiments using reservoir-extracted oil and brine show that differences in relative permeability persist between flow unit core plugs with near-constant porosity, attributable to contrasts in and the spatial arrangement of diagenetic pore types. Core plugs “aged” by exposure to reservoir oil over time exhibit wettability closer to suspected in situ reservoir conditions, compared to “cleaned” core plugs. Together with contact angle measurements, these results suggest that reservoir wettability is transient and modified quickly by oil recovery and carbon storage operations. Reservoir simulation results for enhanced oil recovery, using a five-spot pattern and water-alternating-with-gas injection history at Farnsworth, compare models for cumulative oil and water production using both a single relative permeability determined from history matching, and flow unit-dependent relative permeability determined from experiments herein. Both match cumulative oil production of the field to a satisfactory degree but underestimate historical cumulative water production. Differences in modeled versus observed water production are interpreted in terms of evolving wettability, which we argue is due to the increasing presence of fast paths (flow pathways with connected higher permeability) as the reservoir becomes increasingly water-wet. The control of such fast-paths is thus critical for efficient carbon storage and sweep efficiency for CO2-enhanced oil recovery in heterogeneous reservoirs.

Rasmussen, Lindsey↗

Recent developments in high-temperature photonic crystals for energy conversion

After decades of intense studies focused on cryogenic and room temperature nanophotonics, scientific interest is furthermore growing in high-temperature nanophotonics aimed at solid-state energy conversion. These latest extensive research efforts are spurred by a renewed interest in high temperature thermal-to-electrical energy conversion schemes including thermophotovoltaics (TPV), solar–thermophotovoltaics, solar–thermal, and solar–thermochemical energy conversion systems. This field is profiting tremendously from the outstanding degree of control over the thermal emission properties that can be achieved with nanoscale photonic materials. The key to obtaining high efficiency in this class of high temperature energy conversion is the spectral and angular matching of the radiation properties of an emitter to those of an absorber. Together with the achievements in the field of high-performance narrow bandgap photovoltaic cells, the ability to tailor the radiation properties of thermal emitters and absorbers using nanophotonics facilitates a route to achieving the impressive efficiencies predicted by theoretical studies. In this review, we will discuss the possibilities of emission tailoring by nanophotonics in the light of high temperature thermal-to-electrical energy conversion applications, and give a brief introduction to the field of TPV. We will demonstrate how a class of large area 2D metallic photonic crystals can be designed and employed to efficiently control and tailor the spectral and angular emission properties, paving the way towards new and highly efficient thermophotovoltaic systems and enabling other energy conversion schemes based on high-performance high-temperature nanoscale photonic materials.

30 DIRECT ENERGY CONVERSION↗

The Role of Self-Organized Criticality in the Substorm Phenomenon and its Relation to Localized Reconnection in the Magnetospheric Plasma Sheet

Evidence is presented that suggests there is a significant self-organized criticality (SOC) component in the dynamics of substorms in the magnetosphere. Observations of BBFs, fast flows, localized dipolarizations, plasma turbulence, etc. are taken to show that multiple localized reconnection sites provide the basic avalanche phenomenon in the establishment of SOC in the plasma sheet. First results are presented from a continuing plasma physical study of this avalanche process. A one-dimensional resistive MHD model of a magnetic field reversal is discussed. Resistivity, in this model, is self-consistently generated in response to the excitation of an idealized current-driven instability. When forced by convection of magnetic flux into the field reversal region, the model yields rapid magnetic field annihilation through a dynamic behavior that is shown to exhibit many of the characteristics of SOC. Over a large range of forcing strengths, the annihilation rate is shown to self-adjust to balance the rate at which flux is convected into the reversal region. Several analogies to magnetotail dynamics are discussed: (1) It is shown that the presence of a localized criticality in the model produces a remarkable stability in the global configuration of the field reversal while simultaneously exciting extraordinarily dynamic internal evolution. (2) Under steady forcing, it is shown that a loading-unloading cycle may arise that, as a consequence of the global stability, is quasi-periodic and, therefore, predictable despite the presence of internal turbulence in the field distribution. Indeed, it is shown that the global loading-unloading cycle is a consequence of the internal turbulence. (3) It is shown that, under steady, strong forcing the loading-unloading cycle vanishes. Instead, a recovery from a single unloading persists indefinitely. The field reversal is globally very steady while internally it is very dynamic as field annihilation goes on at the rate necessary to match the strong forcing. From this result we speculate that steady magnetospheric convection events result when the plasma sheet has been driven close to criticality over an extended spatial domain. During these events, we would expect to find localized reconnection sites distributed over the spatial domain of near criticality and we would expect to find plasma sheet transport in that domain to be closely related to that of BBF and fast flow events.

Klimas, Alex J.↗

Origin, Bulk Chemical Composition and Physical Structure of the Galilean Satellites of Jupiter: A Post-Galileo Analysis

The origin of Jupiter and the Galilean satellite system is examined in the light of the new data that has been obtained by the NASA Galileo Project. In particular, special attention is given to a theory of satellite origin which was put forward at the start of the Galileo Mission and on the basis of which several predictions have now been proven successful. These predictions concern the chemical composition of Jupiter's atmosphere and the physical structure of the satellites. According to the proposed theory of satellite origin, each of the Galilean satellites formed by chemical condensation and gravitational accumulation of solid grains within a concentric family of orbiting gas rings. These rings were cast off equatorially by the rotating proto-Jovian cloud (PJC which contracted gravitationally to form Jupiter some 4 1/2 billion years ago. The PJC formed from the gas and grains left over from the gas ring that had been shed at Jupiter's orbit by the contracting proto-solar cloud (PSC Supersonic turbulent convection provides the means for shedding discrete gas rings. The temperatures T (sub n) of the system of gas rings shed by the PSC and PJC vary with their respective mean orbital radii R (sub n) (n = 0, 1,2,...) according as T (sub n) proportional to R (sub n) (exp -0.9). If the planet Mercury condenses at 1640 K, so accounting for the high density of that planet via a process of chemical fractionation between iron and silicates, then T (sub n) at Jupiter's orbit is 158 K. Only 35% of the water vapour condenses out. Thus fractionation between rock and ice, together with an enhancement in the abundance of solids relative to gas which takes place through gravitational sedimentation of solids onto the mean orbit of the gas ring, ensures nearly equal proportions of rock and ice in each of Ganymede and Callisto. Io and Europa condense above the H20 ice point and consist solely of hydrated rock (h-rock). The Ganymedan condensate consists of h-rock and H20 ice. For Callisto, NH3 ice makes up -5% of the condensate mass next to h-rock (approximately 50%) and H2O ice (approximately 45%). Detailed thermal and structural models for each of Europa, Ganymede and Callisto are constructed on the basis of the above initial bulk chemical compositions. For Europa (E), a predicted 2-zone model consisting of a dehydrated rock core of mass 0.912 M (sub E) and a 150 km thick frozen mantle of salty H2O yields a moment-of-inertia coefficient which matches the Galileo Orbiter gravity measurement. For Ganymede (G), a 3-zone model possessing an inner core of solid FeS and mass approximately 0.116 M (sub G), and an outer H2O ice mantle of mass approximately 0.502 M (sub G) is needed to explain the gravity data. Ganymede's native magnetic field was formed by thermoremanent magnetization of Fe3O4. A new Callisto (C) model is proposed consisting of a core of mass 0.826 M (sub C) containing a uniform mixture of h-rock (60% by mass) and H2O and NH3 ices, and capped by a mantle of pure ice. This model may have the capacity to yield a thin layer of liquid NH3 (raised dot) 2H2O at the core boundary, in line with Galileo's discovery of an induced magnetic field.

Prentice, A. J. R.↗

An Application of Multi-band Forced Photometry to One Square Degree of SERVS: Accurate Photometric Redshifts and Implications for Future Science

We apply The Tractor image modeling code to improve upon existing multi-band photometry for the Spitzer Extragalactic Representative Volume Survey (SERVS). SERVS consists of post-cryogenic Spitzer observations at 3.6 and 4.5 μ m over five well-studied deep fields spanning 18 deg{sup 2}. In concert with data from ground-based near-infrared (NIR) and optical surveys, SERVS aims to provide a census of the properties of massive galaxies out to z ≈ 5. To accomplish this, we are using The Tractor to perform “forced photometry.” This technique employs prior measurements of source positions and surface brightness profiles from a high-resolution fiducial band from the VISTA Deep Extragalactic Observations survey to model and fit the fluxes at lower-resolution bands. We discuss our implementation of The Tractor over a square-degree test region within the XMM Large Scale Structure field with deep imaging in 12 NIR/optical bands. Our new multi-band source catalogs offer a number of advantages over traditional position-matched catalogs, including (1) consistent source cross-identification between bands, (2) de-blending of sources that are clearly resolved in the fiducial band but blended in the lower resolution SERVS data, (3) a higher source detection fraction in each band, (4) a larger number of candidate galaxies in the redshift range 5 < z < 6, and (5) a statistically significant improvement in the photometric redshift accuracy as evidenced by the significant decrease in the fraction of outliers compared to spectroscopic redshifts. Thus, forced photometry using The Tractor offers a means of improving the accuracy of multi-band extragalactic surveys designed for galaxy evolution studies. We will extend our application of this technique to the full SERVS footprint in the future.

79 ASTRONOMY AND ASTROPHYSICS↗

Renormalization and matching for the Collins-Soper kernel from lattice QCD

The Collins-Soper kernel, which governs the energy evolution of transverse- momentum dependent parton distribution functions (TMDPDFs), is required to accurately predict Drell-Yan like processes at small transverse momentum, and is a key ingredient for extracting TMDPDFs from experiment. Earlier we proposed a method to calculate this kernel from ratios of the so-called quasi-TMDPDFs determined with lattice QCD, which are defined as hadronic matrix elements of staple-shaped Euclidean Wilson line operators. Here we provide the one-loop renormalization of these operators in a regularization-independent momentum subtraction (RI'/MOM) scheme, as well as the conversion factor from the RI'/MOM-renormalized quasi-TMDPDF to the \( \overline{\mathrm{MS}} \) scheme. We also propose a procedure for calculating the Collins-Soper kernel directly from position space correlators, which simplifies the lattice determination.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗