Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distributed clustering methods”

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

A Hybrid DC Fault Primary Protection Algorithm for Multi-Terminal HVdc Systems

Protection against dc faults is one of the main technical hurdles faced when operating converter-based HVdc systems. Protection becomes even more challenging for multi-terminal dc (MTdc) systems with more than two terminals/converter stations. In this paper, a hybrid primary fault detection algorithm for MTdc systems is proposed to detect a broad range of failures. Sensor measurements, i.e., line currents and dc reactor voltages measured at local terminals, are first processed by a top-level context clustering algorithm. For each cluster, the best fault detector is selected among a detector pool according to a rule resulting from a learning algorithm. The detector pool consists of several existing detection algorithms, each performing differently across fault scenarios. The proposed hybrid primary detection algorithm: i) offers superior performance compared to an individual detector through a data-driven approach; ii) detects all major fault types including pole-to-pole (P2P), pole-to-ground (P2G), and external dc faults; iii) identifies faults with various fault locations and impedances; iv) is more robust to noisy sensor measurements compared to existing methods; v) does not require exhaustive simulation and sampling for training the model. Performance and effectiveness of the proposed algorithm are evaluated and verified based on time-domain simulations in the PSCAD/EMTDC software environment. The results confirm satisfactory operation, accuracy, and detection speed of the proposed algorithm under various fault scenarios.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Relationship between magnetoresistance behavior and magnetic states in intercalated compounds Fe x TiS 2

Here, temperature and field-dependent neutron powder diffraction (NPD) measurements have been performed to reveal the nature of the unusual evolution of the magnetoresistance behavior with increasing Fe content in the intercalated compounds Fe x TiS 2 (with x=0.25,0.33,0.50,0.55) synthesized by solid-phase reaction method with prolonged homogenization heat treatment. As derived from neutron diffraction measurements, both the Fe 0.25 TiS 2 and the Fe 0.50 TiS 2 compound exhibit an antiferromagnetic (AFM) order below their respective Néel temperatures T N ≈52 K and T N ≈140K, which results in the presence of a large magnetoresistance accompanying the field-induced phase transition from AFM to the ferromagnetic (FM) state. At low temperatures, this AFM-FM transition is irreversible, confirmed by the irreversibility of changes in the NPD patterns and the presence of remnant magnetoresistance. In contrast, Fe 0.33 TiS 2 shows short-range magnetic order at T f ≈44 K due to a triangular network of intercalated Fe atoms and frustrations of exchange interactions with a field-induced FM alignment of Fe magnetic moments in an applied magnetic field as revealed by NPD measurements. The field-induced transformations of the cluster glass magnetic state in this compound lead to a significant decrease in electrical resistivity. According to NPD data, a reduced impact of an external magnetic field on the electrical resistivity of the compound Fe 0.55 TiS 2 can be ascribed to the presence of ferromagnetic or ferrimagnetic order in compounds with the Fe concentrations above x=0.50. The results obtained indicate that the distribution of the Fe atoms, along with their concentration in Fe x TiS 2 layered compounds, plays a decisive role in the formation of the magnetic state and the behavior of the magnetoresistance.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Illustrating galaxy–halo connection in the DESI era with illustrisTNG

ABSTRACT We employ the hydrodynamical simulation illustrisTNG to inform the galaxy–halo connection of the Luminous Red Galaxy (LRG) and Emission Line Galaxy (ELG) samples of the Dark Energy Spectroscopic Instrument (DESI) survey at redshift z ∼ 0.8. Specifically, we model the galaxy colours of illustrisTNG and apply sliding DESI colour–magnitude cuts, matching the DESI target densities. We study the halo occupation distribution (HOD) model of the selected samples by matching them to their corresponding dark matter haloes in the illustrisTNG dark matter run. We find the HOD of both the LRG and ELG samples to be consistent with their respective baseline models, but also we find important deviations from common assumptions about the satellite distribution, velocity bias, and galaxy secondary biases. We identify strong evidence for concentration-based and environment-based occupational variance in both samples, an effect known as ‘galaxy assembly bias’. The central and satellite galaxies have distinct dependencies on secondary halo properties, showing that centrals and satellites have distinct evolutionary trajectories and should be modelled separately. These results serve to inform the necessary complexities in modelling galaxy–halo connection for DESI analyses and also prepare for building high-fidelity mock galaxies. Finally, we present a shuffling-based clustering analysis that reveals a 10–15 ${{\ \rm per\ cent}}$ excess in the LRG clustering of modest statistical significance due to secondary galaxy biases. We also find a similar excess signature for the ELGs, but with much lower statistical significance. When a larger hydrodynamical simulation volume becomes available, we expect our analysis pipeline to pinpoint the exact sources of such excess clustering signatures.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

A multi-dimensional parametric study of variability in multi-phase flow dynamics during geologic CO 2 sequestration accelerated with machine learning

Successful geologic CO 2 storage projects depend on numerical simulations to predict reservoir performance during site selection, injection verification, and post-injection monitoring phases of the project. These numerical simulations solve non-linear sets of coupled partial differential equations, while accounting for multi-phase fluid dynamics on the basis of constitutive equations that are embedded into the solution scheme. As a consequence, individual simulations often require tens to hundreds of hours to complete on high-performance computing clusters. Moreover, laboratory experiments reveal that parametric functions for capillary pressure and relative permeability exhibit substantial variability, even within the same rock type. This combination of computational expense and wide-ranging parametric variability means that there remains substantial uncertainty in the behavior of multi-phase CO 2 -water systems, particularly in the context of feedbacks between relative permeability and capillary pressure. To bridge this knowledge gap, here we develop a novel workflow that utilizes physics-based numerical simulation to train an artificial neural network (ANN) emulator for interrogating the multivariate parameter space that governs both capillary pressure and relative permeability. With this approach, the ANN is trained to emulate both fluid pressure distribution and CO 2 saturation, which are then interrogated quantitatively to generate parametric response surface mappings with high-fidelity resolution. Results from this study initially show that capillary entry pressure is the dominant control on both CO 2 plume geometry and fluid pressure propagation when considering the combined effects of capillary pressure and relative permeability, particularly when phase interference is low and residual CO 2 saturation is high. Moreover, the ANN emulator provides tremendous computational speed-up by computing 2691 individual simulations in several minutes; whereas, the same simulation ensemble would have required ~3 years of simulation time using only physics-based simulation methods (25,000 times speed up).

58 GEOSCIENCES↗
DOE OSTI2024

Impact of survey spatial variability on galaxy redshift distributions and the cosmological 3 × 2-point statistics for the Rubin Legacy Survey of Space and Time (LSST)

We investigate the impact of spatial survey non-uniformity on the galaxy redshift distributions for forthcoming data releases of the Rubin Observatory Legacy Survey of Space and Time (LSST). Specifically, we construct a mock photometry data set degraded by the Rubin OpSim observing conditions, and estimate photometric redshifts of the sample using a template-fitting photo-z estimator, BPZ, and a machine learning method, FlexZBoost. We select the Gold sample, defined as $i\lt 25.3$ for 10 yr LSST data, with an adjusted magnitude cut for each year and divide it into five tomographic redshift bins for the weak lensing lens and source samples. We quantify the change in the number of objects, mean redshift, and width of each tomographic bin as a function of the coadd i-band depth for 1-yr (Y1), 3-yr (Y3), and 5-yr (Y5) data. In particular, Y3 and Y5 have large non-uniformity due to the rolling cadence of LSST, hence provide a worst-case scenario of the impact from non-uniformity. We find that these quantities typically increase with depth, and the variation can be $10\!-\!40~{{\rm per\ cent}}$ at extreme depth values. Using Y3 as an example, we propagate the variable depth effect to the weak lensing $3\times 2$ pt analysis, and assess the impact on cosmological parameters via a Fisher forecast. We find that galaxy clustering is most susceptible to variable depth, and non-uniformity needs to be mitigated below 3 per cent to recover unbiased cosmological constraints. There is little impact on galaxy–shear and shear–shear power spectra, given the expected LSST Y3 noise.

cosmology↗
DOE OSTI2020

Differential Power Processing for Ultra-Efficient Data Storage

Here this paper presents the hardware, software, and power codesign of an ultra-efficient data storage server with differential power processing (DPP). DPP can reduce the power conversion stress, improve the efficiency, and enhance the functionality of modular power electronics systems. The power inputs of a large number of hard disk drives (HDDs) were connected in series and supported by a multiport ac-coupled differential power processing (MAC-DPP) converter through a multiwinding transformer. Methods for controlling the multi-input multi-output power flow in the multiwinding transformer while avoiding core saturation were investigated. A ten-port MAC-DPP prototype with 700-W/in 3 power density was built to support a 450-W HDD storage system with ten series-stacked voltage domains. The prototype was tested on a 50-HDD server testbench, and the overall system loss is below 1 W (99.77% system efficiency). The server was able to maintain high-speed reading and writing operation of all 50 HDDs against the worst hot-swapping scenarios. A variety of hardware/software configurations and many cloud storage techniques were tested on the fully functioning server. Experimental results show that the energy efficiency of large-scale information systems (CPU/GPU clusters, memory banks, HDD arrays, etc.) can be greatly improved by software, hardware, and power codesign.

42 ENGINEERING↗
DOE OSTI2023

Quantifying uncertainties due to irreducible three-body forces in deuteron-nucleus reactions

Deuteron-induced nuclear reactions are an essential tool for probing the structure of nuclei as well as astrophysical information such as (n, γ) cross sections. The deuteron-nucleus system is typically described within a Faddeev three-body model consisting of a neutron (n), a proton (p), and the target nucleus (A) interacting through pairwise phenomenological potentials. While Faddeev techniques enable the exact description of the three-body dynamics, their predictive power is limited in part by the omission of irreducible neutron-proton-nucleus three-body force (n–p–A 3BF). Here, our goal is to quantify systematic uncertainties stemming from the reduction of deuteron-nucleus (d + A) dynamics to a picture of three pointlike nuclear clusters interacting via pairwise nucleon-nucleus forces, using as testing grounds d + α scattering and the 6 Li ground state. We particularly focus on quantifying uncertainties arising from the full antisymmetrization of the (A + 2)-body system with the target nucleus fixed in its ground state. We adopt the ab initio no-core shell model coupled with the resonating group method (NCSM/RGM) to compute microscopic n–α and p–α interactions, and use them in a three-body description of the d + α system by means of momentum-space Faddeev-type equations. Simultaneously, we also carry out ab initio calculations of d + α scattering and 6 Li ground state by means of six-body NCSM/RGM calculations to serve as a benchmark for the three-body model predictions given by the Faddeev calculations. By comparing the Faddeev and NCSM/RGM results, we show that the irreducible n–p–α 3BF has a non-negligible effect on bound state and scattering observables alike. Specifically, the Faddeev approach yields a 6 Li ground state that is approximately 600 keV shallower than the one obtained with the NCSM/RGM. Additionally, the Faddeev calculations for d + α scattering yield a 3 + resonance that is located approximately 400 keV higher in energy compared to the NCSM/RGM result. The shape of the d + α angular distributions computed using the two approaches also differ, owing to the discrepancy in the predictions of the 3 + resonance energy. The Faddeev three-body model predictions for d + α scattering and 6 Li using microscopic n–α and p–α potentials differ from those computed microscopically with the NCSM/RGM. These discrepancies are due to the n–p–α 3BF, which arises from two-nucleon exchange terms in the microscopic d–α interaction and are not accounted for in the three-body model Faddeev calculations. This study lays the foundation for future parametrizations of the 3BF due to Pauli exclusion principle effects in improved three-body calculations of deuteron-induced reactions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRSDate not supplied

Predicting the Seawater Chemistry of an Ocean World Using Machine Learning on Isotopic Measurements of Volatile CO2

Introduction: Given the long time intervals required for data transmission to and from ocean worlds targets, low bandwidth for data transmission, time required for data processing and analysis, and potentially extreme radiation environments (e.g., Europa), it is clear that ocean worlds missions will need more autonomous flight instruments and software in order to achieve established science goals. Protracted time intervals for data analysis (e.g., Europa Lander) strongly motivates the development of rapid, consistent and streamlined methods for interpreting data from flight mass spectrometers to e.g., determine how mass spectra from a plume or surface liquid/ice relates to the surface/subsurface. Since mass spectrometry also has the potential to correctly identify biosignatures[1], it is imperative that such methods for interpreting data are consistent and accurate. We used 848 isotope ratio mass spectra from laboratory analyses of CO2 that interacted with ocean worlds-relevant seawaters as a ‘training’ dataset for ‘unsupervised’ machine learning. In unsupervised learning, characteristics of the data are not labeled or linked, and any similarities found only result from the neural network. CO2 isotopologues analyzed for this dataset mimic the remote measurements of CO2 by a flight mass spectrometer, and are detailed in Theiling [2]. From this dataset, we used measured features of the spectra, such as retention time, intensity, and (isotopologue) mass ratios as inputs for our autoencoder neural network. Our neural network was trained to find similarities in these and other spectral features for seawaters of a particular composition and amount of initial CO2. Successful training then created an output of these similarities for various seawaters, which included MgSO4, Na2SO4, NaCl, MgCl2, KCl, and NaHCO3, and combinations of these salts. We then applied dimensionality reduction techniques such as Principal Component Analysis (PCA), T-Distributed Stochastic Neighbor Embedding (TSNE), and Uniform Manifold Approximation and Projection (UMAP) to demonstrate latent data features as a two-dimensional projection in a unitless, high-dimensional space. In this projection, a data point represents the combined effect of spectral features such as intensity, retention time, and isotope ratio. Our initial UMAP demonstrates data clustering (organization of the data by the neural network) based on the amount of CO2 that had initially interacted with each seawater. Further training using more ‘supervised’ learning techniques demonstrate strong clustering of preliminary data based on initial CO2 concentration, seawater chemical composition, and ionic strength (salinity). Our preliminary work therefore suggests that machine learning has the potential to identify compositional variants of an ocean world seawater based on mass spectra from volatile CO2 measurements. Acknowledgments: This work was funded through a Strategic Task Group at NASA Goddard Space Flight Center. The training dataset was collected through funding from the Oklahoma Space Grant Consortium. References: [1] Pappalardo, R. et al. (2013) Astrobiology, 13, 740–773. [2] Theiling (2020) Icarus, 114216.

Europa↗
DOE OSTI2020

Development and Clustering of Rate-Oriented Load Metrics for Customer Price-Plan Analysis

One of the few methods electric utilities can use to motivate and change customer energy consumption is through retail rate structures. Utilities are increasingly moving toward more dynamic rate plans to encourage energy conservation, utilization of onsite renewable generation, peak demand reduction and flattening of demand profiles. This paper creates a set of rate-oriented load metrics that are the determinants of customers' bills under four unique rate plans. These metrics are not only indicative of which rate structure can provide customer bill reductions based on their load profile characteristics, but also convey useful information about load consumption behavior. With these metrics, utilities can analyze their customers and identify classes that are rewarded under each rate plan. This can help inform utilities whether the customers rewarded under each rate plans are meeting their original objectives. To develop these customer classes, we calculate these rate-oriented load metrics for each customer and perform k-means clustering. The analysis is conducted on a set of 300 customer profiles, examining four different rate plans, different numbers of clusters, customer bills and cluster load profile characteristics.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Symphony: Cosmological Zoom-in Simulation Suites over Four Decades of Host Halo Mass

Abstract We present Symphony, a compilation of 262 cosmological, cold-dark-matter-only zoom-in simulations spanning four decades of host halo mass, from 10 11 –10 15 M ⊙ . This compilation includes three existing simulation suites at the cluster and Milky Way–mass scales, and two new suites: 39 Large Magellanic Cloud-mass (10 11 M ⊙ ) and 49 strong-lens-analog (10 13 M ⊙ ) group-mass hosts. Across the entire host halo mass range, the highest-resolution regions in these simulations are resolved with a dark matter particle mass of ≈3 × 10 −7 times the host virial mass and a Plummer-equivalent gravitational softening length of ≈9 × 10 −4 times the host virial radius, on average. We measure correlations between subhalo abundance and host concentration, formation time, and maximum subhalo mass, all of which peak at the Milky Way host halo mass scale. Subhalo abundances are ≈50% higher in clusters than in lower-mass hosts at fixed sub-to-host halo mass ratios. Subhalo radial distributions are approximately self-similar as a function of host mass and are less concentrated than hosts’ underlying dark matter distributions. We compare our results to the semianalytic model Galacticus , which predicts subhalo mass functions with a higher normalization at the low-mass end and radial distributions that are slightly more concentrated than Symphony. We use UniverseMachine to model halo and subhalo star formation histories in Symphony, and we demonstrate that these predictions resolve the formation histories of the halos that host nearly all currently observable satellite galaxies in the universe. To promote open use of Symphony, data products are publicly available at http://web.stanford.edu/group/gfc/symphony .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

On the Use of Field RR Lyrae As Galactic Probes: IV. New Insights Into and Around the Oosterhoff Dichotomy

We discuss the largest and most homogeneous spectroscopic data set of field RR Lyrae variables (RRLs) available to date. We estimated abundances using both high-resolution and low-resolution ($Δ$S method) spectra for fundamental (RRab) and first overtone (RRc) RRLs. The iron abundances for 7941 RRLs were supplemented with similar estimates that are available in the literature, ending up with 9015 RRLs (6150 RRab, 2865 RRc). The metallicity distribution shows a mean value of $\langle$[Fe/H]$\rangle$ = -1.51 ± 0.01, and $σ$(standard deviation) = 0.41 dex with a long metal-poor tail approaching [Fe/H] ≃ - 3 and a sharp metal-rich tail approaching solar iron abundance. The RRab variables are more metal-rich ( $\langle$[Fe/H]$\rangle$ ab = -1.48 ± 0.01, $σ$ = 0.41 dex) than RRc variables ($\langle$[Fe/H]$\rangle$ c = -1.58 ± 0.01, $σ$ = 0.40 dex). The relative fraction of RRab variables in the Bailey diagram (visual amplitude versus period) located along the short-period (more metal-rich) and the long-period (more metalpoor) sequences are 80% and 20%, while RRc variables display an opposite trend, namely 30% and 70%, respectively. We found that the pulsation period of both RRab and RRc variables steadily decreases when moving from the metal-poor to the metal-rich regime. Furthermore, the visual amplitude shows the same trend, but RRc amplitudes are almost two times more sensitive than RRab amplitudes to metallicity. We also investigated the dependence of the population ratio (N $c$ /N tot ) of field RRLs on the metallicity and we found that the distribution is more complex than in globular clusters. The population ratio steadily increases from ~0.25 to ~0.36 in the metal-poor regime, it decreases from ~0.36 to ~0.18 for -1.8 ≤ [Fe/H] ≤ -0.9 and it increases to a value of ~0.3 approaching solar iron abundance.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Zinc Speciation and Propane Dehydrogenation in Zn/H-ZSM-5 Catalysts

Zn/H-ZSM-5 catalysts have been frequently investigated for propane dehydrogenation (PDH); however, the active site remains unresolved due to the complexity of the system. We employed in situ FTIR spectroscopy and a kinetics method to correlate the Zn speciation and PDH activity in Zn/H-ZSM-5 with two Si/Al ratios (15 and 39) and a range of Zn/Al ratios (0–1.7). Incremental additions of zinc show that Zn 2+ sites are preferentially formed on H-ZSM-5 over a fraction of paired Al sites followed by [Zn-O-Zn] 2+ and [ZnOH] + sites and then ZnO x clusters. The [Zn-OH] + and [Zn-O-Zn] 2+ sites in H-ZSM-5 are more active and selective than isolated Zn 2+ for PDH. [Zn-OH] + species sublimate over time on stream, leading to catalyst deactivation, while [Zn-O-Zn] 2+ species are stable even after high-temperature reduction (750 °C for 60 min). Three distinct Zn sites ([Zn-O-Zn] 2+ , Zn 2+ , and [ZnOH] + ) show a similar propane reaction order (close to 1) and H 2 reaction order (close to 0). Combined with the lack of Zn hydride when propane flows over the catalyst at 550 °C, it is concluded that propane adsorption and dissociation is a rate-determining step and H 2 desorption is fast. This work indicates that the preparation of H-ZSM-5 with abundant Al pairs may be a strategy to form stable and selective Zn/H-ZSM-5 catalysts for propane dehydrogenation. Furthermore, it is also highlighted that examining the effect of both metal/Al ratios and Al distribution of the zeolite is crucial in identifying the metal cations in metal–zeolite systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2019

An ALMA+ACA measurement of the shock in the Bullet Cluster

Context: The thermal Sunyaev-Zeldovich (SZ) effect presents a relatively new tool for characterizing galaxy cluster merger shocks, traditionally studied through X-ray observations. Widely regarded as the “textbook example” of a cluster merger bow shock, the western, most-prominent shock front in the Bullet Cluster (1E0657-56) represents the ideal test case for such an SZ study. Aims: We aim to characterize the shock properties using deep, high-resolution interferometric SZ effect observations in combination with priors from an independent X-ray analysis. Methods: Our analysis technique relies on the reconstruction of a parametric model for the SZ signal by directly and jointly fitting data from the Atacama Large Millimeter/submillimeter Array (ALMA) and Atacama Compact Array (ACA) in Fourier space. Results: The ALMA+ACA data are primarily sensitive to the electron pressure difference across the shock front. To estimate the shock Mach number M, this difference can be combined with the value for the upstream electron pressure derived from an independent Chandra X-ray analysis. In the case of instantaneous electron-ion temperature equilibration, we find M = 2.08 (+0.12,−0.12), in ≈ 2.4σ tension with the independent constraint from Chandra, M(X)= 2.74 ± 0.25. The assumption of purely adiabatic electron temperature change across the shock leads to M = 2.53 (+0.33,−0.25), in better agreement with the X-ray estimate M(X) = 2.57 ± 0.23 derived for the same heating scenario. Conclusion: We have demonstrated that interferometric observations of the thermal SZ effect provide constraints on the properties of the shock in the Bullet Cluster that are highly complementary to X-ray observations. The combination of X-ray and SZ data yields a powerful probe of the shock properties, capable of measuring M and addressing the question of electron-ion equilibration in cluster shocks. Our analysis is however limited by systematics related to the overall cluster geometry and the complexity of the post-shock gas distribution. To overcome these limitations, a simultaneous, joint-likelihood analysis of SZ and X-ray data is needed.

Luca Di Mascolo↗
NASA NTRS2021

Atacama Cosmology Telescope measurements of a large sample of candidates from the Massive and Distant Clusters of WISE Survey: Sunyaev-Zeldovich effect confirmation of MaDCoWS candidates using ACT

Context. Galaxy clusters are an important tool for cosmology, and their detection and characterization are key goals for current and future surveys. Using data from the Wide-field Infrared Survey Explorer (WISE), the Massive and Distant Clusters of WISE Survey (MaDCoWS) located 2839 significant galaxy overdensities at redshifts 0.7 . z . 1.5, which included extensive follow-up imaging from the Spitzer Space Telescope to determine cluster richnesses. Concurrently, the Atacama Cosmology Telescope (ACT) has produced large area millimeter-wave maps in three frequency bands along with a large catalog of Sunyaev-Zeldovich (SZ)-selected clusters as part of its Data Release 5 (DR5). Aims. We aim to verify and characterize MaDCoWS clusters using measurements of, or limits on, their thermal SZ effect signatures. We also use these detections to establish the scaling relation between SZ mass and the MaDCoWS-defined richness. Methods. Using the maps and cluster catalog from DR5, we explore the scaling between SZ mass and cluster richness. We do this by comparing cataloged detections and extracting individual and stacked SZ signals from the MaDCoWS cluster locations. We use complementary radio survey data from the Very Large Array, submillimeter data from Herschel, and ACT 224 GHz data to assess the impact of contaminating sources on the SZ signals from both ACT and MaDCoWS clusters. We use a hierarchical Bayesian model to fit the mass-richness scaling relation, allowing for clusters to be drawn from two populations: one, a Gaussian centered on the mass-richness relation, and the other, a Gaussian centered on zero SZ signal. Results. We find that MaDCoWS clusters have submillimeter contamination that is consistent with a gray-body spectrum, while the ACT clusters are consistent with no submillimeter emission on average. Additionally, the intrinsic radio intensities of ACT clusters are lower than those of MaDCoWS clusters, even when the ACT clusters are restricted to the same redshift range as the MaDCoWS clusters. We find the best-fit ACT SZ mass versus MaDCoWS richness scaling relation has a slope of p1 = 1.84+0.15 −0.14, where the slope is defined as M ∝ λ p1 15 and λ15 is the richness. We also find that the ACT SZ signals for a significant fraction (∼57%) of the MaDCoWS sample can statistically be described as being drawn from a noise-like distribution, indicating that the candidates are possibly dominated by low-mass and unvirialized systems that are below the mass limit of the ACT sample. Further, we note that a large portion of the optically confirmed ACT clusters located in the same volume of the sky as MaDCoWS are not selected by MaDCoWS, indicating that the MaDCoWS sample is not complete with respect to SZ selection. Finally, we find that the radio loud fraction of MaDCoWS clusters increases with richness, while we find no evidence that the submillimeter emission of the MaDCoWS clusters evolves with richness. Conclusions. We conclude that the original MaDCoWS selection function is not well defined and, as such, reiterate the MaDCoWS collaboration’s recommendation that the sample is suited for probing cluster and galaxy evolution, but not cosmological analyses. We find a best-fit mass-richness relation slope that agrees with the published MaDCoWS preliminary results. Additionally, we find that while the approximate level of infill of the ACT and MaDCoWS cluster SZ signals (1–2%) is subdominant to other sources of uncertainty for current generation experiments, characterizing and removing this bias will be critical for next-generation experiments hoping to constrain cluster masses at the sub-percent level.

large↗
NASA NTRSDate not supplied

Simulation of Radiation-Induced DNA Damage With the Code RITRACKS

INTRODUCTION DNA damage is one of the most physiologically important effects of ionizing radiation. Clustered DNA damage events, like double-strand breaks (DSBs), have the most notable biological consequences. DNA damage types depend on both the track structure of the radiation and the spatial organization of the DNA. High linear energy transfer (LET) charged nuclei, found in galactic cosmic rays (GCR), are known to produce large numbers of complex DNA damage events. The human genome is packaged into chromatin, which can take on locus-dependent and cell type-dependent spatial conformations that correspond to epigenetic states, such as more open, extended structures in transcriptionally active chromatin. These epigenetic differences can affect DNA break patterns in response to ionizing radiation, potentially creating distinct DNA repair and signaling outcomes across the genome in different cells. MATERIAL AND METHODS The code RITRACKS (Relativistic Ion Tracks), which simulates stochastic radiation track structures and radiation chemistry, was used to model damage on isolated and histone-bound DNA by various types of ions and photons. The changes made to the code to perform radiation-induced DNA damage, and simulation results on single nucleosomes are given in our recent paper. In this work, the DNA building capabilities of RITRACKS have been extended to simulate more complex DNA structures build on the coarse-grain simulation framework meso-WLCsim. This code can sample generic chromatin fiber conformation ensembles based on the geometry of nucleosomes and mechanical properties of DNA. Using RITRACKS, we simulated the fragment length distributions (FLD) of irradiated DNA structures built using the chromatin conformations of WLCsim and obtained results representative of those obtained with Radiation-Induced Correlated Cleavage with sequencing (RICC-Seq) experiments [6]. We have also performed Fe ion and photon irradiations of K562, IMR90, BJ and RPE-1 cells at NSRL to experimentally validate results. Sample processing and data analysis are in progress and any available preliminary results will be discussed. DISCUSSION The recent updates in the code RITRACKS allow the calculation of several quantities such as the DNA damage yield and the FLD. This approach can be used to model epigenetic state-specific chromatin structure parameters to leverage the epigenetic state data available for many human cell types to infer relative DNA damage sensitivity among genomic loci.

I Plante↗
NASA NTRS2001

Extrapolation of the dna fragment-size distribution after high-dose irradiation to predict effects at low doses

The patterns of DSBs induced in the genome are different for sparsely and densely ionizing radiations: In the former case, the patterns are well described by a random-breakage model; in the latter, a more sophisticated tool is needed. We used a Monte Carlo algorithm with a random-walk geometry of chromatin, and a track structure defined by the radial distribution of energy deposition from an incident ion, to fit the PFGE data for fragment-size distribution after high-dose irradiation. These fits determined the unknown parameters of the model, enabling the extrapolation of data for high-dose irradiation to the low doses that are relevant for NASA space radiation research. The randomly-located-clusters formalism was used to speed the simulations. It was shown that only one adjustable parameter, Q, the track efficiency parameter, was necessary to predict DNA fragment sizes for wide ranges of doses. This parameter was determined for a variety of radiations and LETs and was used to predict the DSB patterns at the HPRT locus of the human X chromosome after low-dose irradiation. It was found that high-LET radiation would be more likely than low-LET radiation to induce additional DSBs within the HPRT gene if this gene already contained one DSB.

NASA Discipline Radiation Health↗
NASA NTRS2022

X-Ray Morphology of Cluster-Mass Haloes in Self-Interacting Dark Matter

We perform cosmological zoom-in simulations of 19 relaxed cluster-mass haloes with the inclusion of adiabatic gas in the cold dark matter (CDM) and self-interacting dark matter (SIDM) models. These clusters are selected as dynamically relaxed clusters from a parent simulation with M200(1-3) × 10 15 M. Both the dark matter and the intracluster gas distributions in SIDM appear more spherical than their CDM counterparts. Mock X-ray images are generated based on the simulations and are compared to the real X-ray images of 84 relaxed clusters selected from the Chandra and ROSAT archives. We perform ellipse fitting for the isophotes of mock and real X-ray images and obtain the ellipticities at cluster-centric radii of r0.1–0.2R 200 .The X-ray isophotes in SIDM models with increasing cross-sections are rounder than their CDM counterparts, which manifests as asystematic shift in the distribution function of ellipticities. Unexpectedly, the X-ray morphology of the observed non-cool-coreclusters agrees better with SIDM model swith cross-section (σ/m)=0.5–1cm2g−1 than CDM and SIDM with (σ/m)=0.1cm2g−1. Our statistical analysis indicates that the latter two models are disfavoured at the 68 percent confidence level (as conservative estimates). This conclusion is not altered by shifting the radial range of measurements or applying a temperature selection criterion. However, the primary uncertainty originates from the lack of baryonic physics in the adiabatic model, such as cooling, star formation and feedback effects, which still have the potential to reconcile CDM simulations with observations.

dark matter – cosmology↗
DOE OSTI2022

Tracer-field cross-correlations with k -nearest neighbour distributions

ABSTRACT In astronomy and cosmology significant effort is devoted to characterizing and understanding spatial cross-correlations between points – e.g galaxy positions, high energy neutrino arrival directions, X-ray and AGN sources, and continuous fields – e.g. weak lensing meiand Cosmic Microwave Background maps. Recently, we introduced the k-nearest neighbour (kNN) formalism to better characterize the clustering of discrete (point) data sets. Here, we extend it to the point – field cross-correlations analysis. It combines kNN measurements of the point data set with measurements of the field smoothed at many scales. The resulting statistics are sensitive to all orders in the joint clustering of the points and the field. We demonstrate that this approach, unlike the 2-pt cross-correlation, can measure the statistical dependence of two data sets even when there are no linear (Gaussian) correlations between them. We further demonstrate that this framework is far more effective than the two point function in detecting cross-correlations when the continuous field is contaminated by high levels of noise. For a particularly high level of noise, the cross-correlation between haloes and the underlying matter field in a cosmological simulation, between 10 h−1 Mpc and 30 h−1 Mpc, is detected at >5σ significance using the technique presented here, when the two-point cross-correlation significance is ∼1σ. Finally, we show that kNN cross-correlations of haloes and the matter field can be well modelled on quasi-linear scales using the Hybrid Effective Field Theory (HEFT) framework, with the same set of bias parameters as are used for 2-pt cross-correlations. The substantial improvement in the statistical power of detecting cross-correlations using this method makes it a promising tool for various cosmological applications.

79 ASTRONOMY AND ASTROPHYSICS↗