Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “cluster finding”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,531 records · Page 85
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

The behavior of methane–water mixtures under elevated pressures from simulations using many-body potentials

Non-polarizable empirical potentials have been proven to be incapable of capturing the mixing of methane–water mixtures at elevated pressures. Although density functional theory-based ab initio simulations may circumvent this discrepancy, they are limited in terms of the relevant time and length scales associated with mixing phenomena. Here, we show that the many-body MB-nrg potential, designed to reproduce methane–water interactions with coupled cluster accuracy, successfully captures this phenomenon up to 3 GPa and 500 K with varying methane concentrations. Two-phase simulations and long time scales that are required to fully capture the mixing, affordable due to the speed and accuracy of the MBX software, are assessed. Constructing the methane–water equation of state across the phase diagram shows that the stable mixtures are denser than the sum of their parts at a given pressure and temperature. As a result, we find that many-body polarization plays a central role, enhancing the induced dipole moments of methane by 0.20 D during mixing under pressure. Overall, the mixed system adopts a denser state, which involves a significant enthalpic driving force as elucidated by a systematic many-body energy decomposition analysis.

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↗
DOE OSTI2023

The Spatial Clustering of ROSAT All-Sky Survey Active Galactic Nuclei. V. The Evolution of Broad-line AGN Clustering Properties in the Last 6 Gyr

This is the fifth paper in a series of investigations of the clustering properties of luminous, broad-emission-line active galactic nuclei (AGNs) identified in the ROSAT All-Sky Survey (RASS) and Sloan Digital Sky Survey (SDSS). In this work we measure the cross-correlation function (CCF) between RASS/SDSS Data Release 14 AGNs with the SDSS CMASS galaxy sample at 0.44 < z < 0.64. We apply halo occupation distribution (HOD) modeling to the CCF along with the autocorrelation function of the CMASS galaxies. We find that X-ray-selected and optically selected AGNs at 0.44 < z < 0.64 reside in statistically identical halos with a typical dark matter halo (DMH) mass of $M$$^{typ,AGN}_{DMH}$ ~ 10 12.7 h –1 M ⊙ . The acceptable HOD parameter space for these two broad-line AGN samples have only statistically marginal differences caused by small deviations of the CCFs in the one-halo-dominated regime on small scales. In contrast to optically selected AGNs, the X-ray AGN sample may contain a larger population of satellites at MDMH ~ 1013h–1M⊙. We compare our measurements in this work with our earlier studies at lower independent redshift ranges, spanning a lookback time of 6 Gyr. The comparison over this wider redshift range of 0.07 < z < 0.64 reveals (i) no significant difference between the typical DMH masses of X-ray-selected and optically selected AGNs, (ii) weak positive clustering dependencies of $M$$^{typ,AGN}_{DMH}$ with L X and M BH , (iii) no significant dependence of $M$$^{typ,AGN}_{DMH}$ on Eddington ratio, and (iv) the same DMH masses host more-massive accreting black holes at high redshift than at low redshifts.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

The active CGCG 077-102 NED02 galaxy within the Abell 2063 galaxy cluster

Context.Within the framework of investigating the link between the central super massive black holes in the cores of galaxies and the galaxies themselves, we detected a variable X-ray source in the center of CGCG 077-102 NED02, which is a member of the CGCG 077-102 galaxy pair within the Abell 2063 cluster of galaxies. Aims.Our goal is to combine X-ray and optical data to demonstrate that this object harbors an active super massive black hole in its core, and to relate this to the dynamical status of the galaxy pair within the Abell 2063 cluster. Methods.We usedChandraandXMM-Newtonarchival data to derive the X-ray spectral shape and variability. We also obtained optical spectroscopy to detect the expected emission lines that are typically found in active galactic nuclei. Finally, we used public ZTF imaging data to investigate the optical variability. Results.There is no evidence of multiple X-ray sources or extended components within CGCG 077-102 NED02. Single X-ray spectral models fit the source well. We detect significant, nonrandom inter-observation 0.5–10 keV X-ray flux variabilities, for observations separated by ∼4 days for short-term variations and by up to ∼700 days for long-term variations. Optical spectroscopy points toward a passive galaxy for CGCG 077-102 NED01 and a Seyfert for CGCG 077-102 NED02. The classification of CGCG 077-102 NED02 is also consistent with its X-ray luminosity of over 10 42 erg s −1 . We do not detect short-term variability in the optical ZTF light curves. However, we find a significant long-term stochastic variability in theg-band that can be well described by the damped random walk model with a best-fit characteristic damping timescale ofτ DRW = 30 −12 +28 days. Finally, the CGCG 077-102 galaxy pair is deeply embedded within the Abell 2063 potential, with a long enough history within this massive structure to have been affected by the influence of this cluster for a long time. Conclusions.Our observations point toward a moderately massive black hole in the center of CGCG 077-102 NED02 of ∼10 6 M ⊙ . As compared to another similar pair in the literature, CGCG 077-102 NED02 is not heavily obscured, perhaps because of the surrounding intracluster medium ram-pressure stripping.

Astronomy & Astrophysics↗
DOE OSTI2024

Defect Production and Microstructural Feature Impact for Radiation Damage in Additively Manufactured 316 Stainless Steel

This milestone presents multi-scale modeling research results for additively manufactured 316 stainless steel. A combination of phase field, cluster dynamics, molecular dynamics, and density functional theory with machine learning is used, allowing for predictions of radiation-driven microstructural evolution in additively manufactured 316 stainless steel over a range of temperatures, damage rates, neutron spectra, and microstructures, and supporting the development of combined ion and neutron irradiation for material qualification. Informed by ion irradiation and neutron irradiation results across the Advanced Materials and Manufacturing Technologies program, we investigate the unique aspects of radiation-driven microstructure evolution in additively manufactured 316 stainless steel. In particular, we focus on understanding the impact of carbon concentration (varying, for example, between the 316L and 316H standards) on void formation; radiation-induced segregation at dislocation cells and grain boundaries; and the evolution of dislocation loops and network dislocation populations. We find that the unique characteristics of the additively manufactured microstructures must be accounted for in understanding the evolution of dislocation populations under thermal and irradiation conditions, such as the variation in sink strengths arising due to the variation in dislocation density. We also find that the radiation-induced segregation of Cr and Ni to grain boundaries and cell walls differs due to the differences in their defect sink biases. We also find that increasing the Ni content can slow vacancy diffusion, which may provide a mechanism for the observed reduction in transient swelling rate for austenitic Fe-Cr-Ni alloys with increasing Ni content. In addition, ion irradiations have shown that increasing carbon content in 316 SS results in a larger population of smaller voids, suggesting reduced vacancy diffusion. Our results show that the carbon content of additively manufactured 316 SS has a significant impact on the migration rate of defect clusters. The presence of carbon atoms results in carbon-vacancy trapping, significantly reducing the diffusion rate of vacancies. Carbon atoms may also be trapped near the surface of a void, which may reduce void growth by trapping vacancies that diffuse toward the void.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Fiducial-cosmology-dependent systematics for the DESI 2024 full-shape analysis

We assess the impact of the fiducial cosmology choice on cosmological inference from full-shape (FS) fits of the galaxy power spectrum in the DESI 2024 Data Release 1 (DR1). Using a suite of AbacusSummit DR1 mock catalogues based on the Planck 2018 best-fit cosmology, we quantify potential systematic shifts introduced by analysing the data under five secondary cosmologies — featuring variations in matter density, thawing dark energy, higher effective number of neutrino species, reduced clustering amplitude, and the DESI DR1 BAO best-fit w 0 w a CDM cosmology — relative to DESI's baseline Planck 2018 cosmology. We investigate two complementary FS analysis approaches: full-modelling (FM) and ShapeFit (SF), each with distinct sensitivities to the assumed fiducial model. Across all tracers, we find for FM that systematic shifts induced by fiducial cosmology mismatches remain well below the DESI DR1 statistical uncertainties, with maximum deviations of 0.22σ DR1 in ΛCDM scenarios and 0.12σ DR1+SN when including SN Ia mock data in extended w 0 w a CDM fits. For SF, the shifts in the compressed parameters remain below 0.45σ DR1 for all tracers and cosmologies.

dark energy experiments↗
DOE OSTI2022

Discovering hidden geothermal signatures using non-negative matrix factorization with customized k-means clustering

Discovery of hidden geothermal resources is challenging. It requires the mining of large datasets with diverse data attributes representing subsurface hydrogeological and geothermal conditions. The commonly used play fairway analysis approach typically incorporates subject-matter expertise to analyze regional data to estimate geothermal characteristics and favorability. We demonstrate an alternative approach based on machine learning (ML) to process a geothermal dataset from southwest New Mexico (SWNM). The study region includes low- and medium-temperature hydrothermal systems. Several of these systems are not well characterized because of insufficient existing data and limited past explorative work. This study discovers hidden patterns and relations in the SWNM geothermal dataset to improve our understanding of the regional hydrothermal conditions and energy-production favorability. This understanding is obtained by applying an unsupervised ML algorithm based on non-negative matrix factorization coupled with customized k-means clustering (NMFk). NMFk can automatically identify (1) hidden signatures characterizing analyzed datasets, (2) the optimal number of these signatures, (3) the dominant data attributes associated with each signature, and (4) the spatial distribution of the extracted signatures. Here, in this study, NMFk is applied to analyze 18 geological, geophysical, hydrogeological, and geothermal attributes at 44 locations in SWNM. Using NMFk, we find data patterns and identify the spatial associations of hydrothermal signatures within two physiographic provinces (Colorado Plateau and Basin and Range) and two sub-regions of these provinces (the Mogollon-Datil volcanic field and the Rio Grande rift) in SWNM. The ML algorithm extracted five hydrothermal signatures in the SWNM datasets that differentiate between low (<90°C) and medium (90-150°C)-temperature hydrothermal systems. The algorithm also suggests that the Rio Grande rift and northern Mogollon-Datil volcanic field are the most favorable regions for future geothermal resource discovery. NMFk also identified critical attributes to identify medium-temperature hydrothermal systems in the study area. The resulting NMFk model can be applied to predict geothermal conditions and their uncertainties at new SWNM locations based on limited data from unexplored regions. The code to execute the performed analyses as well as the corresponding data can be found at https://github.com/SmartTensors/GeoThermalCloud.jl.

15 GEOTHERMAL ENERGY↗
DOE OSTI2026

Reconstruction and Dissolution of Copper Catalysts during Electrocatalytic Nitrate Reduction

Copper-based electrocatalysts are widely explored for electrochemical nitrate remediation, yet their stability under operating conditions remains poorly understood. While nanoscale and subnanoscale Cu motifs are known to restructure during the nitrate reduction reaction (NO 3 RR), how these transformations translate into irreversible material loss remains unclear. Here, we quantify Cu dissolution during NO 3 RR as a function of catalyst architecture and electrolyte chemistry using two model systems: single-atom Cu (Cu 1 ) and Cu nanoparticles (Cu NP ). Time-resolved leaching measurements, integrated with in situ X-ray absorption spectroscopy (XAS), reveal measurable Cu loss for both catalysts during NO 3 RR. Dissolution is concentrated at the initiation of electrolysis, coinciding with rapid restructuring, and depends strongly on morphology, with Cu NP consistently exhibiting greater Cu loss than Cu 1 . In situ XAS reveals that Cu 1 forms transient metallic clusters under reduction that largely redisperse upon returning to open-circuit voltage, whereas Cu NP undergoes reduction of an oxidized surface layer accompanied by sustained Cu loss during electrolysis. Notably, the presence of nitrate significantly intensifies restructuring and Cu loss, highlighting the critical role of electrolyte composition. Furthermore, these findings establish a direct link between electrochemical restructuring and Cu dissolution, unveiling electrolyte-dependent interactions as key determinants of catalyst durability in electrochemical nitrate conversion.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Third-Order Møller–Plesset Theory Made More Useful? The Role of Density Functional Theory Orbitals

The practical utility of Møller–Plesset (MP) perturbation theory is severely constrained by the use of Hartree–Fock (HF) orbitals. It has recently been shown that the use of regularized orbital-optimized MP2 orbitals and scaling of MP3 energy could lead to a significant reduction in MP3 error. In this work, we examine whether density functional theory (DFT)-optimized orbitals can be similarly employed to improve the performance of MP theory at both the MP2 and MP3 levels. We find that the use of DFT orbitals leads to significantly improved performance for prediction of thermochemistry, barrier heights, noncovalent interactions, and dipole moments relative to the standard HF-based MP theory. Indeed, MP3 (with or without scaling) with DFT orbitals is found to surpass the accuracy of coupled-cluster singles and doubles (CCSD) for several data sets. We also found that the results are not particularly functional sensitive in most cases (although range-separated hybrid functionals with low delocalization error perform the best). As such, MP3 based on DFT orbitals thus appears to be an efficient, noniterative O(N 6 ) scaling wave-function approach for single-reference electronic structure computations. Scaled MP2 with DFT orbitals is also found to be quite accurate in many cases, although modern double hybrid functionals are likely to be considerably more accurate.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Using Density-Corrected DFT to Understand Density-Driven and Functional-Dependent Errors in Ab Initio Simulations of the Hydrated Electron

The hydrated electron, an excess electron in liquid water, plays a crucial role in a plethora of chemical processes, motivating extensive research efforts to characterize its structure, dynamics, and reactivity in solution. Recent theoretical approaches to understanding this intriguing object have involved ab initio simulations based on density functional theory (DFT). Although DFT allows for the study of hydrated electron reactivity and quantum mechanical behavior, it is well-known that anionic systems can suffer from significant density-driven errors (DDEs). Density-corrected DFT (DC-DFT) provides a framework to mitigate such errors; the method reduces DDEs by replacing the self-consistent (SC) density associated with a given density functional with the Hartree–Fock (HF) density. Since HF densities tend to be more localized than DFT SC densities, the DC-DFT scheme significantly improves errors in calculations where the SC density is spuriously delocalized. Here, we investigate how the use of density correction affects the calculated properties of the DFT-simulated (PBEh) hydrated electron, a particularly challenging diffuse anionic system to simulate. First, we analyze charge delocalization in a system consisting of a model octahedral hydrated electron water cluster (the so-called Kevan structure) along with a spatially separated sulfur atom. We show that the use of density correction indeed reduces DDEs in comparison to a standard DFT global hybrid functional. We then propagate molecular dynamics trajectories of the hydrated electron using DC-DFT, where we find that DC further localizes electron density in the cavity region, a signature of reduced charge delocalization. Unfortunately, the decreased radius of gyration of the spin density and corresponding tightening of the local solvation structure from density correction causes predicted observables to deviate further from experimental measurements than when density correction is not employed. Here, we argue that DC’s worse agreement with experiment results from the removal of a fortuitous cancellation of errors that is intrinsic to the PBEh functional. This indicates that the difficulties with DFT to simulate hydrated electrons are primarily due to the inherent approximations in DFT rather than to density-driven errors.

Density functional theory↗
DOE OSTI2023

Mastering morphology of non-fullerene acceptors towards long-term stable organic solar cells

Despite the rapid progress of organic solar cells based on non-fullerene acceptors, simultaneously achieving high power conversion efficiency and long-term stability for commercialization requires sustainable research effort. Here, we demonstrate stable devices by integrating a wide bandgap electron-donating polymer (namely PTzBI-dF) and two acceptors (namely L8BO and Y6) that feature similar structures yet different thermal and morphological properties. The organic solar cell based on PTzBI-dF:L8BO:Y6 could achieve a promising efficiency of 18.26% in the conventional device structure. In the inverted structure, excellent long-term thermal stability over 1400 h under 85 °C continuous heating is obtained. The improved performance can be ascribed to suppressed charge recombination along with appropriate charge transport. We find that the morphological features in terms of crystalline coherence length of fresh and aged films can be gradually regulated by the weight ratio of L8BO:Y6. Additionally, the occurrence of melting point decrease and reduced enthalpy in PTzBI-dF:L8BO:Y6 films could prohibit the amorphous phase to cluster, and consequently overcome the energetic traps accumulation aroused by thermal stress, which is a critical issue in high efficiency non-fullerene acceptors-based devices. This work provides insight into understanding non-fullerene acceptors-based organic solar cells for improved efficiency and stability.

14 SOLAR ENERGY↗
DOE OSTI2022

Variational quantum eigensolver with reduced circuit complexity

Abstract The variational quantum eigensolver (VQE) is one of the most promising algorithms to find eigenstates of a given Hamiltonian on noisy intermediate-scale quantum devices (NISQ). The practical realization is limited by the complexity of quantum circuits. Here we present an approach to reduce quantum circuit complexity in VQE for electronic structure calculations. Our ClusterVQE algorithm splits the initial qubit space into clusters which are further distributed on individual (shallower) quantum circuits. The clusters are obtained based on mutual information reflecting maximal entanglement between qubits, whereas inter-cluster correlation is taken into account via a new “dressed” Hamiltonian. ClusterVQE therefore allows exact simulation of the problem by using fewer qubits and shallower circuit depth at the cost of additional classical resources, making it a potential leader for quantum chemistry simulations on NISQ devices. Proof-of-principle demonstrations are presented for several molecular systems based on quantum simulators as well as IBM quantum devices.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Computational multiphysics modeling of radioactive aerosol deposition in diverse human respiratory tract geometries

The evaluation of aerosol exposure relies on generic mathematical models that assume uniform particle deposition profiles over the human respiratory tract and do not account for subject-specific characteristics. Here we introduce a hybrid-automated computational workflow that generates personalized particle deposition profiles in 3D reconstructed human airways from computed tomography scans using Computational Fluid and Particle Dynamics simulations. This is the first large-scale study to consider realistic airways variability, where 380 lower and 40 upper human respiratory tract 3D geometries are reconstructed and parameterized. The data is clustered into nine groups using random forest regression. Computational fluid and particle dynamics simulations are conducted on these representative geometries using a realistic heavy-breathing respiratory cycle and radioactive iodine-131 as a source term. Monte Carlo radiation transport simulations are performed to obtain detailed energy deposition maps. Our findings emphasize the importance of personalized studies, as minor respiratory tract variations notably influence deposition patterns rather than global parameters of the lower airways, observing more than 30% variance in the mass deposition fraction.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2022

Accurate core excitation and ionization energies from a state-specific coupled-cluster singles and doubles approach

Here, we investigate the use of orbital-optimized references in conjunction with single-reference coupled-cluster theory with single and double substitutions (CCSD) for the study of core excitations and ionizations of 18 small organic molecules, without the use of response theory or equation-of-motion (EOM) formalisms. Three schemes are employed to successfully address the convergence difficulties associated with the coupled-cluster equations, and the spin contamination resulting from the use of a spin symmetry-broken reference, in the case of excitations. In order to gauge the inherent potential of the methods studied, an effort is made to provide reasonable basis set limit estimates for the transition energies. Overall, we find that the two best-performing schemes studied here for ΔCCSD are capable of predicting excitation and ionization energies with errors comparable to experimental accuracies. The proposed ΔCCSD schemes reduces statistical errors against experimental excitation energies by more than a factor of two when compared to the frozen-core core–valence separated (FC-CVS) EOM-CCSD approach – a successful variant of EOM-CCSD tailored towards core excitations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Primordial non-Gaussianities with weak lensing: information on non-linear scales in the Ulagam full-sky simulations

Abstract Primordial non-Gaussianities (PNGs) are signatures in the density field that encode particle physics processes from the inflationary epoch. Such signatures have been extensively studied using the Cosmic Microwave Background, through constraining their amplitudes,f X NL , with future improvements expected from large-scale structure surveys; specifically, the galaxy correlation functions. We show that weak lensing fields can be used to achieve competitive and complementary constraints. This is shown via theUlagamsuite of N-body simulations, a subset of which evolves primordial fields with four types of PNGs. We create full-sky lensing maps and estimate the Fisher information from three summary statistics measured on the maps: the moments, the cumulative distribution function, and the 3-point correlation function. We find that the year 10 sample from the Rubin Observatory Legacy Survey of Space and Time (LSST) can constrain PNGs toσ(f NL eq ) ≈ 110,σ(f NL or, lss ) ≈ 120,σ(f NL loc ) ≈ 40. For the former two, this is better than or comparable to expected galaxy clustering-based constraints from the Dark Energy Spectroscopic Instrument (DESI). The PNG information in lensing fields is on non-linear scales and at low redshifts (z≲ 1.25), with a clear origin in the evolution history of massive halos. The constraining power degrades by ∼60% under scale cuts of ≳ 20 Mpc, showing there is still significant information on scales mostly insensitive to small-scale systematic effects (e.g., baryons). We publicly release theUlagamsuite to enable more survey-focused analyses.

Astronomy & Astrophysics↗
DOE OSTI2021

Characterizing the target selection pipeline for the Dark Energy Spectroscopic Instrument Bright Galaxy Survey

ABSTRACT We present the steps taken to produce a reliable and complete input galaxy catalogue for the Dark Energy Spectroscopic Instrument (DESI) Bright Galaxy Survey (BGS) using the photometric Legacy Survey DR8 DECam. We analyse some of the main issues faced in the selection of targets for the DESI BGS, such as star–galaxy separation, contamination by fragmented stars and bright galaxies. Our pipeline utilizes a new way to select BGS galaxies using Gaia photometry and we implement geometrical and photometric masks that reduce the number of spurious objects. The resulting catalogue is cross-matched with the Galaxy And Mass Assembly (GAMA) survey to assess the completeness of the galaxy catalogue and the performance of the target selection. We also validate the clustering of the sources in our BGS catalogue by comparing with mock catalogues and the Sloan Digital Sky Survey (SDSS) data. Finally, the robustness of the BGS selection criteria is assessed by quantifying the dependence of the target galaxy density on imaging and other properties. The largest systematic correlation we find is a 7 per cent suppression of the target density in regions of high stellar density.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Gravothermal collapse and the diversity of galactic rotation curves

The rotation curves of spiral galaxies exhibit a great diversity that challenges our understanding of galaxy formation and the nature of dark matter. Previous studies showed that in self-interacting dark matter (SIDM) models with a cross section per unit mass of σ/m ≈ $\mathscr{O}$(1) cm 2 /g, the predicted dark matter central densities are a good match to the observed densities in galaxies. Here, in this work, we explore a regime with a larger cross section of σ/m ≈ 20−40 cm 2 /g in dwarf galactic halos. We will show that such strong dark matter self-interactions can further amplify the diversity of halo densities inherited from their assembly history. High concentration halos can enter the gravothermal collapse phase within 10 Gyr, resulting in a high density, while low concentration ones remain in the expansion phase and have a low density. We fit the rotation curves of 14 representative low surface brightness galaxies and demonstrate how the large range of observed central densities is naturally accommodated in the strong SIDM regime of σ/m ≈ 20 – 40 cm 2 /g. Galaxies that are outliers in the previous studies, due to their high halo central densities, are no longer outliers in this SIDM regime as their halos would be in the collapse phase. For galaxies with a low density, the SIDM fits are robust to the variation of the cross section. Our findings open up a new window for testing gravothermal collapse, the unique signature of strong dark matter self-interactions, and exploring a broader SIDM model space. As an example, we illustrate how the larger cross sections favored by our fits, together with upper limits from strong lensing observations in clusters, pick out the preferred SIDM model space for a dark matter particle coupled to a light gauge boson in the Born regime.

dark matter↗
DOE OSTI2020

Collective effects and pattern formation for directional locking of disks moving through obstacle arrays

We examine directional locking effects in an assembly of disks driven through a square array of obstacles as the angle of drive rotates from 0° to 90°. For increasing disk densities, the system exhibits a series of different dynamic patterns along certain locking directions, including one-dimensional or multiple-row chain phases and density-modulated phases. For nonlocking driving directions, the disks form disordered patterns or clusters. Here, when the obstacles are small or far apart, a large number of locking phases appear; however, as the number of disks increases, the number of possible locking phases drops due to the increasing frequency of collisions between the disks and obstacles. For dense arrays or large obstacles, we find an increased clogging effect in which immobile and moving disks coexist.

97 MATHEMATICS AND COMPUTING↗