Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “CrS”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

CERF: IM3 Projected Western US Power Plant Locations

Overview The Capacity Expansion Regional Feasibility (CERF) model is an open-source geospatial python package that provides new power plant locations at a 1km resolution. The model ingests U.S. state or regional-scale electricity system capacity expansion plans, such as those produced by the Global Change Analysis Model (GCAM-USA), and identifies feasible, site-specific locations for individual new power plants (renewable and non-renewable). CERF combines high-resolution geospatial suitability analyses with an economic algorithm that selects individual plant siting locations based on grid interconnection costs and the locational marginal value of new generation. The model incorporates a wide range of dynamic constraints and opportunities, such as protected lands, population density, existing infrastructure, and water availability. This dataset provides CERF power plant siting results for IM3 Phase 2 simulations across eight different scenarios for the Western US through 2055. The scenarios include combinations of two Shared Socioeconomic Pathways (SSP3 and SSP5) with four high-resolution climate projections specific to the United States (see, https://tgw-data.msdlive.org/). These climate projections include "hotter" and "cooler" variants for two Representative Concentration Pathways (RCP4.5 and RCP8.5). The resulting eight simulations are: rcp45cooler_ssp3 rcp45cooler_ssp5 rcp45hotter_ssp3 rcp45hotter_ssp5 rcp85cooler_ssp3 rcp85cooler_ssp5 rcp85hotter_ssp3 rcp85hotter_ssp5 CERF siting results in this dataset correspond to capacity expansion plans in the GCAM-USA IM3 Phase 2 simulation data and are available for each of the above scenarios. Data Details Temporal Range: 2015-2055 in 5-year timesteps. Note that 2015 is the experiment base year and 2020 and beyond represent model simulation years. Spatial Range: Plant locations are provided for the eleven states in the Western US including Arizona, California, Colorado, Idaho, Montana, New Mexico, Nevada, Oregon, Utah, Washington, and Wyoming. Spatial Resolution: 1 km-squared, provided in x and y coordinates Geospatial Projection: Albers Equal Area Conic (ESRI:102003) File Type: csv The dataset contains subdirectories for each of the eight scenarios described in the overview. Each scenario folder contains two subfolders with the following information: 1. Power Plant Data This directory contains a single .csv file of power plant locations for both pre-existing (non-CERF sited plants in operation in 2015) and new (CERF-sited) power plants across the temporal range along with additional CERF model output parameters for CERF-sited plants. Plant with a siting year earlier than 2020 correspond to facilities that are operational leading into the first timestep CERF simulation. For a more detailed description of CERF model output parameters, see the CERF model documentation. Note that the cerf_plant_id parameter is unique within each scenario file but not across scenario files. Parameter Descriptions scenario - Name of scenario cerf_plant_id - Unique siting identifier cerf_sited - If True, indicates that plant was sited by CERF model. If False, indicates pre-existing facility region_name - Name of region (state) tech_id - Technology ID tech_name - Full generation technology name inclusive of cooling type (if applicable) and additional characteristics tech_simple - Simplified generation technology type unit_size_mw - Power plant unit size (MW) xcoord - X coordinate in the default CRS (meters) ycoord - Y coordinate in the default CRS (meters) index - Index position in the flattend 2D array buffer_in_km - Exclusion buffer around site (km) sited_year - Year of siting retirement_year - Year of retirement lmp_zone - Locational marginal price (LMP) zone ID locational_marginal_price_usd_per_mwh - Locational marginal price ($/MWh) generation_mwh_per_year - Generation output (MWh/yr) operating_cost_usd_per_year - Cost of plant operations ($/yr) net_operational_value - Net operational value based on LMP and and operating costs ($/yr) interconnection_cost - Cost of interconnection for transmission & gas pipeline (if applicable) net_locational_cost -- Difference of interconnection cost and operating value ($/yr) capacity_factor_fraction - Capacity factor (fraction) carbon_capture_rate_fraction - Carbon capture rate (fraction) fuel_co2_content_tons_per_btu - Fuel CO2 content (tons/Btu) fuel_price_usd_per_mmbtu - Fuel price ($/MMBtu) fuel_price_esc_rate_fraction - Fuel price escalation rate (fraction) heat_rate_btu_per_kWh - Heat rate (Btu/kWh) lifetime_yrs - Technology lifetime for annuity (years) operational_life_yrs - Operational lifetime for retirement (years) variable_om_usd_per_mwh - Variable operation and maintenance costs of yearly capacity use ($/MWh) variable_om_esc_rate_fraction - Variable operation and maintenance costs escalation rate (fraction) carbon_tax_usd_per_ton - Carbon tax ($/ton) carbon_tax_esc_rate_fraction - Carbon tax escalation rate (fraction) 2. Storage Data This directory contains information on new and pre-existing energy storage facilities operational in each timestep along with various storage operational parameters. The 2015 timestep provides pre-existing energy storage data and corresponds with facilities that are operational leading into the first model simulation timestep. Note that coordinates in the storage files correspond to the interconnection point on the grid (substation location), not individual energy storage locations. Energy storage is added in a cumulative process at each given interconnection point. That is, each individual file provides the total operational storage capacity interconnected to the specified substation for the given timestep, inclusive of previously installed storage at that location and new storage installed in that timestep at that location. Parameters scenario - Name of scenario timestep - Simulation timestep name - Unique storage identifier s_typ - Type of energy storage technology (battery or pumped storage hydro) s_node - Node ID of interconnecting substation xcoord - X coordinate in the default CRS (meters) ycoord - Y coordinate in the default CRS (meters) charge_rate - Maximum charge rate (power capacity) of storage system (MW) discharge_rate - Maximum discharge rate (power capacity) of storage system (MW) duration - Duration of storage system (hours) max_SoC - Allowed maximum state of charge (energy capacity) of storage system (MWh) min_SoC -Allowed minimum state of charge (energy capacity) of storage system (MWh) charge_eff - Efficiency of charge (fraction between 0 and 1) discharge_eff - Efficiency of discharge (fraction between 0 and 1) Acknowledgment IM3 is a multi-institutional effort led by Pacific Northwest National Laboratory and supported by the U.S. Department of Energy's Office of Science as part of research in MultiSector Dynamics, Earth and Environmental Systems Modeling Program.

CERF↗
NASA NTRS1990

Particle astrophysics - The NASA cosmic ray program for the 1990s and beyond

Various papers on particle astrophysics are presented. Individual topics addressed include: the NASA cosmic ray (CR) program for the 1990s and beyond, SAMPEX Mission overview, the Advanced Composition Explorer, Positron Electron Magnet Spectrometer for the Eos Mission, Heavy Nucleus Collector for Space Station, the Astromag Facility, Large Isotope Spectrometer for Astromag, the Solar Probe Mission, the Mercury Dual Orbiter Mission, CRs in the heliosphere, origin of high-energy Galactic CRs, CR studies with the Gamma-Ray Observatory, gamma-ray astronomy at 1 TeV, experimental search for point sources above 1 TeV, the UMC Extensive Air Shower Array, status of the MACRO experiment. Also discussed are: CRs above 1 TeV/n and neutrino astronomy, abundance of ultraheavy nuclei in solar energetic particles, CR studies with an interstellar probe, isotopic composition of CR nuclei beyond the iron peak, experimental studies of CR isotopic composition up to Zr-40, use of accelerators in particle astrophysics, development of long-duration ballooning in Antarctica, Lunar-Based Heavy Nucleus Detector, neutrino astronomy on the moon, gamma rays at airplane altitudes, source composition of CRs.

Jones, W. V.↗
NASA NTRS1994

Removing cosmic-ray hits from multiorbit HST Wide Field Camera images

We present an optimized algorithm that removes cosmic rays ('CRs') from multiorbit Hubble Space Telescope (HST) Wide Field/Planetary Camera ('WF/PC') images. It computes the image noise in every iteration from the WF/PC CCD equation. This includes all known sources of random and systematic calibration errors. We test this algorithm on WF/PC stacks of 2-12 orbits as a function of the number of available orbits and the formal Poissonian sigma-clipping level. We find that the algorithm needs greater than or equal 4 WF/PC exposures to locate the minimal sky signal (which is noticeably affected by CRs), with an optimal clipping level at 2-2.5 x sigma(sub Poisson). We analyze the CR flux detected on multiorbit 'CR stacks,' which are constructed by subtracting the best CR filtered images from the unfiltered 8-12 orbit average. We use an automated object finder to determine the surface density of CRS as a function of the apparent magnitude (or ADU flux) they would have generated in the images had they not been removed. The power law slope of the CR 'counts' (gamma approximately = 0.6 for N(m) m(exp gamma)) is steeper than that of the faint galaxy counts down to V approximately = 28 mag. The CR counts show a drop off between 28 less than or approximately V less than or approximately 30 mag (the latter is our formal 2 sigma point source sensitivity without spherical aberration). This prevents the CR sky integral from diverging, and is likely due to a real cutoff in the CR energy distribution below approximately 11 ADU per orbit. The integral CR surface density is less than or approximately 10(exp 8)/sq. deg, and their sky signal is V approximately = 25.5-27.0 mag/sq. arcsec, or 3%-13% of our NEP sky background (V = 23.3 mag/sq. arcsec), and well above the EBL integral of the deepest galaxy counts (B(sub J) approximately = 28.0 mag/sq. arcsec). We conclude that faint CRs will always contribute to the sky signal in the deepest WF/PC images. Since WFPC2 has approximately 2.7x lower read noise and a thicker CCD, this will result in more CR detections than in WF/PC, potentially affecting approximately 10%-20% of the pixels in multiorbit WFPC2 data cubes.

Windhorst, Rogier A.↗
NASA NTRS2007

Validation and Determination of Ice Water Content - Radar Reflectivity Relationships during CRYSTAL-FACE: Flight Requirements for Future Comparisons

In order for clouds to be more accurately represented in global circulation models (GCM), there is need for improved understanding of the properties of ice such as the total water in ice clouds, called ice water content (IWC), ice particle sizes and their shapes. Improved representation of clouds in models will enable GCMs to better predict for example, how changes in emissions of pollutants affect cloud formation and evolution, upper tropospheric water vapor, and the radiative budget of the atmosphere that is crucial for climate change studies. An extensive cloud measurement campaign called CRYSTAL-FACE was conducted during Summer 2002 using instrumented aircraft and a variety of instruments to measure properties of ice clouds. This paper deals with the measurement of IWC using the Harvard water vapor and total water instruments on the NASA WB-57 high-altitude aircraft. The IWC is measured directly by these instruments at the altitude of the WB-57, and it is compared with remote measurements from the Goddard Cloud Radar System (CRS) on the NASA ER-2. CRS measures vertical profiles of radar reflectivity from which IWC can be estimated at the WB-57 altitude. The IWC measurements obtained from the Harvard instruments and CRS were found to be within 20-30% of each other. Part of this difference was attributed to errors associated with comparing two measurements that are not collocated in time an space since both aircraft were not in identical locations. This study provides some credibility to the Harvard and CRS-derived IWC measurements that are in general difficult to validate except through consistency checks using different measurement approaches.

Sayres, D. S.↗
NASA NTRS2016

Deep View of the Large Magellanic Cloud with Six Years of Fermi-LAT Observations

Context. The nearby Large Magellanic Cloud (LMC) provides a rare opportunity of a spatially resolved view of an external star-forming galaxy in gamma-rays. The LMC was detected at 0.1-100 GeV as an extended source with CGRO/EGRET and using early observations with the Fermi-LAT. The emission was found to correlate with massive star-forming regions and to be particularly bright towards 30 Doradus. Aims. Studies of the origin and transport of cosmic rays (CRs) in the Milky Way are frequently hampered by line-of-sight confusion and poor distance determination. The LMC offers a complementary way to address these questions by revealing whether and how the gamma-ray emission is connected to specific objects, populations of objects, and structures in the galaxy. Methods. We revisited the gamma-ray emission from the LMC using about 73 months of Fermi-LAT P7REP data in the 0.2-100 GeV range. We developed a complete spatial and spectral model of the LMC emission, for which we tested several approaches: a simple geometrical description, template-fitting, and a physically driven model for CR-induced interstellar emission. Results. In addition to identifying PSR J0540-6919 through its pulsations, we find two hard sources positionally coincident with plerion N 157B and supernova remnant N 132D, which were also detected at TeV energies with H.E.S.S. We detect an additional soft source that is currently unidentified. Extended emission dominates the total flux from the LMC. It consists of an extended component of about the size of the galaxy and additional emission from three to four regions with degree-scale sizes. If it is interpreted as CRs interacting with interstellar gas, the large-scale emission implies a large-scale population of approximately 1-100 GeV CRs with a density of approximately 30% of the local Galactic value. On top of that, the three to four small-scale emission regions would correspond to enhancements of the CR density by factors 2 to 6 or higher, possibly more energetic and younger populations of CRs compared to the large-scale population. An alternative explanation is that this is emission from an unresolved population of at least two dozen objects, such as pulsars and their nebulae or supernova remnants. This small-scale extended emission has a spatial distribution that does not clearly correlate with known components of the LMC, except for a possible relation to cavities and supergiant shells. Conclusions. The Fermi-LAT GeV observations allowed us to detect individual sources in the LMC. Three of the newly discovered sources are associated with rare and extreme objects. The 30 Doradus region is prominent in GeV gamma-rays because PSR J0540-6919 and N 157B are strong emitters. The extended emission from the galaxy has an unexpected spatial distribution, and observations at higher energies and in radio may help to clarify its origin.

gamma rays: galaxies – Magellanic Clouds – cos↗
NASA NTRS2021

Classifying planetary cloudiness with an updated set of MODIS Cloud Regimes

We present an updated Cloud Regime (CR) dataset based on Moderate resolution Imaging Spectroradiometer (MODIS) Collection 6.1 cloud products, specifically joint histograms that partition cloud fraction within distinct combinations of cloud top pressure and cloud optical thickness ranges. The paper focuses on an edition of the CR dataset derived from our own aggregation of MODIS pixel-level cloud retrievals on an equal area grid and pre-specified 3-hour UTC intervals that spatiotemporally match International Satellite Cloud Climatology Project (ISCCP) gridded cloud data. The other edition comes from the 1-degree daily aggregation provided by standard MODIS Level-3 data, as in previous versions of the MODIS CRs, for easier use with datasets mapped on equal angle grids. Both editions consist of 11 clusters whose centroids are nearly identical.We provide a physical interpretation of the new CRs and aspects of their climatology that have not been previously examined, such as seasonal and interannual variability of CR frequency of occurrence. We also examine the makeup and precipitation properties of the CRs assisted by independent datasets originating from active observations, and provide a first glimpse of how MODIS CRs relate to clouds as seen by ISCCP.

MODIS↗
DOE OSTI2021

Common risk segment mapping: Streamlining exploration for carbon storage sites, with application to coastal Texas and Louisiana

Large-scale deployment of Carbon Capture and Storage (CCS) will require a commensurately large number of sites. Efficient screening methods are needed to create investment assurance and focus efforts on the most promising sites. The problem is similar to petroleum exploration, for which there are well-developed (though seldom published) workflows, including Common Risk Segment (CRS) mapping. In brief, the process requires 1) defining the key play elements; 2) identifying candidate geologic intervals for each; 3) creating fact-based maps for those intervals; 4) determining minimum criteria for the success of each element; 5) reinterpreting the fact-based maps in terms of chance of success; and 6) combining the individual maps to form a composite, basin-scale view of prospectivity. We adapt the CRS process to screening for CO 2 storage sites. Critically, we redefine the process in terms of cost of characterization and development, rather than chance of success. For illustration, we apply the process to the example of the Lower Miocene on the Texas and Louisiana Gulf Coast. We show that the predictions are consistent with historic hydrocarbon production volumes and rates. The power of the CRS method is that it creates a systematic approach to geologic evaluation and translates complex, multidimensional analysis into clear, graphical and easily comprehended business inputs. The results highlight sweet spots and identifies critical risks, suggesting a focus for further data collection and analysis. Furthermore, the method developed here can be applied to both surface and subsurface factors anywhere that there is interest in geologic storage of CO 2 .

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Quantifying leaf symptoms of sorghum charcoal rot in images of field‐grown plants using deep neural networks

Abstract Charcoal rot of sorghum (CRS) is a significant disease affecting sorghum crops, with limited genetic resistance available. The causative agent, Macrophomina phaseolina (Tassi) Goid, is a highly destructive fungal pathogen that targets over 500 plant species globally, including essential staple crops. Utilizing field image data for precise detection and quantification of CRS could greatly assist in the prompt identification and management of affected fields and thereby reduce yield losses. The objective of this work was to implement various machine learning algorithms to evaluate their ability to accurately detect and quantify CRS in red‐green‐blue images of sorghum plants exhibiting symptoms of infection. EfficientNet‐B3 and a fully convolutional network emerged as the top‐performing models for image classification and segmentation tasks, respectively. Among the classification models evaluated, EfficientNet‐B3 demonstrated superior performance, achieving an accuracy of 86.97%, a recall rate of 0.71, and an F1 score of 0.73. Of the segmentation models tested, FCN proved to be the most effective, exhibiting a validation accuracy of 97.76%, a recall rate of 0.68, and an F1 score of 0.66. As the size of the image patches increased, both models’ validation scores increased linearly, and their inference time decreased exponentially. This trend could be attributed to larger patches containing more information, improving model performance, and fewer patches reducing the computational load, thus decreasing inference time. The models, in addition to being immediately useful for breeders and growers of sorghum, advance the domain of automated plant phenotyping and may serve as a foundation for drone‐based or other automated field phenotyping efforts. Additionally, the models presented herein can be accessed through a web‐based application where users can easily analyze their own images.

Gonzalez, Emmanuel M.↗
DOE OSTI2023

Geothermal Play Fairway Analysis, Part 2: GIS methodology

Play Fairway Analysis (PFA) in geothermal exploration originates from a systematic methodology developed within the petroleum industry and is based on a geologic, geophysical, and hydrologic framework of identified geothermal systems. We tailored this methodology to study the geothermal resource potential of the Snake River Plain and surrounding region, but it can be adapted to other geothermal resource settings. We adapted the PFA approach to geothermal resource exploration by cataloging the critical elements controlling exploitable hydrothermal systems, establishing risk matrices that evaluate these elements in terms of both probability of success and level of knowledge, and building a code-based ‘processing model’ to process results. A geographic information system was used to compile a range of different data types, which we refer to as elements (e.g., faults, vents, heat flow, etc.), with distinct characteristics and measures of confidence. Discontinuous discrete data (points, lines, or polygons) for each element were transformed into continuous interpretive 2D grid surfaces called evidence layers. Because different data types have varying uncertainties, most evidence layers have an accompanying confidence layer which reflects spatial variations in these uncertainties. Risk layers, as defined here, are the product of evidence and confidence layers, and are the building blocks used to construct Common Risk Segment (CRS) maps for heat, permeability, and seal, using a weighted sum for permeability and heat, but a different approach with seal. CRS maps quantify the variable risk associated with each of these critical components. In a final step, the three CRS maps were combined into a Composite Common Risk Segment (CCRS) map, using a modified weighted sum, for results that reveal favorable areas for geothermal exploration. Additional maps are also presented that do not mix contributions from evidence and confidence (to allow an isolated view of evidence and confidence), as well as maps that calculate favorability using the product of components instead of a weighted sum (to highlight where all components are present). Our approach helped to identify areas of high geothermal favorability in the western and central Snake River Plain during the first phase of study and helped identify more precise local drilling targets during the second phase of work. By identifying favorable areas, this methodology can help to reduce uncertainty in geothermal energy exploration and development.

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

The distributional impact of FEMA’s community rating system

Community Rating System (CRS) incentivizes investments in risk reduction above NFIP standards using discounts on insurance premiums. These discounts are cross-subsidized by increasing premiums in non-CRS communities. We examine the distribution of these subsidies and find that redistribution does occur, but the gains and losses are not economically large with 95% of households gaining or losing no more than 0.3% of household income. We also examine their relationship with other community characteristics and find that the strongest predictor of premium reductions is the underlying flood risk level within the community. Thus, CRS appears to reduce the cost of living in the riskier communities.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2020

Lightning Interferometry Uncertainty, Beam Steering Interferometry, and Evidence of Lightning Being Ignited by a Cosmic Ray Shower

Abstract We present an uncertainty analysis on correlation‐based time delay estimate, the basis for broadband lightning interferometry. A normal interferometry would yield much higher uncertainties than the theoretically predicted lower bound due to finite correlation window lengths. However, if the signals are aligned in time, the uncertainties approach the lower bound and can be used to indicate the interferometric uncertainties. Based on this, we introduce a beam steering interferometry technique. It first estimates a direction centroid for the lightning sources with a normal interferometry and computes the time delays among the signals. It then shifts the raw data with the time delays to align the signals roughly in time and reprocess the aligned signals. It finally pushes the reprocessed results to their correct positions, based on the time delays estimated in the first step. We apply this technique on a fast positive breakdown (FPB) process that started a normal intracloud lightning. The FPB process is shown to have a much more complex structure and development than a normal interferometry would provide. More importantly, from both interferometric and polarization analyses, the FPB appears to be ignited by a cosmic ray shower (CRS). We estimate the radio frequency strength and frequency content related to a presumed CRS in a thunderstorm electric field and find that they are in agreement with the observations. We examine the electric field effect of the CRS front and find that it could raise the field above the threshold for positive breakdown and is capable of igniting the FPB (and the lightning).

Shao, Xuan‐Min↗
DOE OSTI2025

3D Radio Frequency Mapping and Polarization Observations Show Lightning Flashes Were Ignited by Cosmic‐Ray Showers

Previous 2D radio frequency interferometric observations showed that lightning is commonly started with a positive fast discharge (+FD) before a normal negative leader continues from the origin of the +FD to the ensuing lightning flash. However, the inception and development of the +FD cannot be convincingly explained by existing discharge theories. With our new 3D broadband interferometric mapping and polarization system, we observed that the +FD was sometimes followed by an even faster and more extensive negative discharge (−FD) that propagated backward and overshoot the origin by a few hundred meters, as reported by an earlier 2D study. Surprisingly, the signal polarization, which measures the orientation of the discharge current, systematically slanted from the discharge propagation direction and rotated between the two opposite discharges, showing the +FD and −FD were driven by other storm-independent factors in addition to the storm electric field, or otherwise all would align in the same direction. Assuming a cosmic-ray shower (CRS) piercing through the cloud immediately before the +FD/−FD discharge, we found that their path is consistent with a pre-ionized path by the CRS, and the polarizations for the two opposite discharges are consistent with the respective deflected trajectories of high-energy positrons and electrons in the geomagnetic and an electric field. We further analyzed the more commonly observed +FD and showed that it is consistent with the CRS interpretation, suggesting these flashes in thunderstorms were ignited by cosmic-ray showers.

58 GEOSCIENCES↗
DOE OSTI2024

Ultraviolet H 2 luminescence in molecular clouds induced by cosmic rays

Context. Galactic cosmic rays (CRs) play a crucial role in ionisation, dissociation, and excitation processes within dense cloud regions where UV radiation is absorbed by dust grains and gas species. CRs regulate the abundance of ions and radicals, leading to the formation of more and more complex molecular species, and determine the charge distribution on dust grains. A quantitative analysis of these effects is essential for understanding the dynamical and chemical evolution of star-forming regions. Aims. The CR-induced photon flux has a significant impact on the evolution of the dense molecular medium in its gas and dust components. This study evaluates the flux of UV photons generated by CRs to calculate the photon-induced dissociation and ionisation rates of a vast number of atomic and molecular species, as well as the integrated UV photon flux. Methods. To achieve these goals, we took advantage of recent developments in the determination of the spectra of secondary electrons, in the calculation of state-resolved excitation cross sections of H 2 by electron impact, and of photodissociation and photoionisation cross sections. Results. We calculated the H 2 level population of each rovibrational level of the X, B, C, B′, D, B″, D′, and a states. We then computed the UV photon spectrum of H 2 in its line and continuum components between 72 and 700 nm, with unprecedented accuracy, as a function of the CR spectrum incident on a molecular cloud, the H 2 column density, the isomeric H 2 composition, and the dust properties. The resulting photodissociation and photoionisation rates are, on average, lower than previous determinations by a factor of about 2, with deviations of up to a factor of 5 for the photodissociation of species such as AlH, C 2 H 2 , C 2 H 3 , C 3 H 3 , LiH, N 2 , NaCl, NaH, O 2 + , S2, SiH, l-C 4 , and l-C 5 H. A special focus is given to the photoionisation rates of H 2 , HF, and N 2 , as well as to the photodissociation of H 2 , which we find to be orders of magnitude higher than previous estimates. We give parameterisations for both the photorates and the integrated UV photon flux as a function of the CR ionisation rate, which implicitly depends on the H 2 column density, as well as the dust properties.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Chemical Recommender System: Replacement Suggestions for Small Molecules

The Chemical Recommender System (CRS) is an open-source, high-performance toolkit that enables real-time similarity searches across the complete PubChem database (over 50 million molecules) using commodity hardware. The CRS addresses critical limitations in existing chemical informatics platforms through a novel vector database infrastructure, extensible model integration capabilities, and complete algorithmic transparency. The system implements a vector database deployment with partitioned indexing that achieves a ~60x speedup over traditional approaches. A containerized model integration framework allows researchers to seamlessly incorporate custom predictive models into the full-scale search and scoring pipeline, while complete configurability of search parameters, filtering logic, and scoring functions provides capabilities not available in existing black-box solutions. Beyond structural similarity, the CRS integrates OPERA QSAR models for thermophysical and toxicity predictions, RDKit synthetic accessibility scoring, and user-defined models to compute weighted final replacement scores. The complete system is accessible through an interactive web application supporting real-time progress monitoring, post-processing score re-weighting, automated PDF reporting, and batch processing capabilities.

Nair, Parthiv Anand [Sandia National Laboratories ↗
DOE OSTI2024

Ultraviolet H 2 luminescence in molecular clouds induced by cosmic rays

Context. Galactic cosmic rays (CRs) play a crucial role in ionisation, dissociation, and excitation processes within dense cloud regions where UV radiation is absorbed by dust grains and gas species. CRs regulate the abundance of ions and radicals, leading to the formation of more and more complex molecular species, and determine the charge distribution on dust grains. A quantitative analysis of these effects is essential for understanding the dynamical and chemical evolution of star-forming regions. Aims. The CR-induced photon flux has a significant impact on the evolution of the dense molecular medium in its gas and dust components. This study evaluates the flux of UV photons generated by CRs to calculate the photon-induced dissociation and ionisation rates of a vast number of atomic and molecular species, as well as the integrated UV photon flux. Methods. To achieve these goals, we took advantage of recent developments in the determination of the spectra of secondary electrons, in the calculation of state-resolved excitation cross sections of H 2 by electron impact, and of photodissociation and photoionisation cross sections. Results. We calculated the H 2 level population of each rovibrational level of the X, B, C, B', D, B", D', and a states. We then computed the UV photon spectrum of H 2 in its line and continuum components between 72 and 700 nm, with unprecedented accuracy, as a function of the CR spectrum incident on a molecular cloud, the H 2 column density, the isomeric H 2 composition, and the dust properties. The resulting photodissociation and photoionisation rates are, on average, lower than previous determinations by a factor of about 2, with deviations of up to a factor of 5 for the photodissociation of species such as AlH, C 2 H 2 , C 2 H 3 , C 3 H 3 , LiH, N 2 , NaCl, NaH, O$^+_2$, S 2 , SiH, l-C 4 , and l-C 5 H. A special focus is given to the photoionisation rates of H 2 , HF, and N 2 , as well as to the photodissociation of H 2 , which we find to be orders of magnitude higher than previous estimates. We give parameterisations for both the photorates and the integrated UV photon flux as a function of the CR ionisation rate, which implicitly depends on the H 2 column density, as well as the dust properties.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Transport of Protostellar Cosmic Rays in Turbulent Dense Cores

We report recent studies have suggested that low-energy cosmic rays (CRs) may be accelerated inside molecular clouds by the shocks associated with star formation. We use a Monte Carlo transport code to model the propagation of CRs accelerated by protostellar accretion shocks through protostellar cores. We calculate the CR attenuation and energy losses and compute the resulting flux and ionization rate as a function of both radial distance from the protostar and angular position. We show that protostellar cores have nonuniform CR fluxes that produce a broad range of CR ionization rates, with the maximum value being up to two orders of magnitude higher than the radial average at a given distance. In particular, the CR flux is focused in the direction of the outflow cavity, creating a "flashlight" effect and allowing CRs to leak out of the core. The radially averaged ionization rates are less than the measured value for the Milky Way of ζ ≈ 10 -16 s -1 ; however, within r ≈ 0.03 pc from the protostar, the maximum ionization rates exceed this value. We show that variation in the protostellar parameters, particularly in the accretion rate, may produce ionization rates that are a couple of orders of magnitude higher or lower than our fiducial values. Finally, we use a statistical method to model unresolved subgrid magnetic turbulence in the core. We show that turbulence modifies the CR spectrum and increases the uniformity of the CR distribution but does not significantly affect the resulting ionization rates.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Reacceleration of Galactic Cosmic Rays beyond the Knee at the Termination Shock of a Cosmic-Ray-driven Galactic Wind

The origin of cosmic rays (CRs) above the knee in the spectrum is an unsolved problem. We present a wind model in which interstellar gas flows along a nonrotating, expanding flux tube with a changing speed and cross-sectional area. CRs from Galactic sources, such as supernova remnants, which are coupled to the plasma via Alfvén waves, provide the main pressure source for driving this outflow. These CRs are then subject to diffusive shock reacceleration at the Galactic wind termination shock, which is located at a distance ~200 kpc. Some of the highest-energy reaccelerated particles propagate upstream against the wind and can contribute to the petaelectronvolt to exaelectronvolt range of the spectrum. We analyze the conditions under which efficient reacceleration can occur and find that rigidities ~10–40 PV can be obtained and that the termination shock may account for between 10% and 50% of the proton spectrum measured in IceCube/IceTop experiment. In our model, the termination shock is unable to fully explain the CR spectrum in the petaelectronvolt to exaelectronvolt range. The highest-energy particles that escape downstream from our termination shock, and similar shocks surrounding most galaxies, can be further accelerated by intergalactic shock fronts.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

A Unified Model of Cosmic-Ray Propagation and Radio Extreme Scattering Events from Intermittent Interstellar Structures

Intermittent magnetic structures are a plausible candidate for explaining cosmic-ray (CR) diffusion rates derived from observed CR energy spectra. Independently, studies of extreme scattering events (ESEs) of radio quasars and pulsar scintillation have hinted that very straight, large aspect ratio magnetic current sheets may be responsible for the localized large scattering of radio waves. The required shortest axis of the typical structures producing ESEs is of the same scale (∼au) as the gyroradii of ∼GeV CRs. In this Letter, we propose that the same magnetic/density sheets can produce large scattering of both CRs and radio waves. We demonstrate that the geometry and volume-filling factor of the sheets derived from quasar ESEs can explain the observed mean free path of GeV CRs without introducing free parameters. The model places constraints on the sheet geometry, such as straightness and large aspect ratio, and assumes that the statistics of the sheets are similar throughout the Galactic volume. We therefore discuss observational tests of the sheet model, which includes observations of echoes in pulsars and fast radio bursts, gravitationally lensed quasars, the distribution of ESE durations, and spatial correlations between ESE events and rotation measure fluctuations. Such tests will be enabled by upcoming wide-field radio instruments, including the Canadian Hydrogen Observatory and Radio-transient Detector and Deep Synoptic Array 2000 Antennas.

Kempski, Philipp [Princeton Univ., NJ (United Stat↗