Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model fitting”

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 991 records · Page 55

3D structure and in situ arrangements of CatSper channel in the sperm flagellum

The sperm calcium channel CatSper plays a central role in successful fertilization as a primary Ca 2+ gateway. Here, we applied cryo-electron tomography to visualize the higher-order organization of the native CatSper complex in intact mammalian sperm. The repeating CatSper units form long zigzag-rows along mouse and human sperm flagella. Above each tetrameric channel pore, most of the extracellular domains form a canopy that interconnects to a zigzag-shaped roof. Murine CatSper contains an additional wing-structure connected to the tetrameric channel. The intracellular domains link two neighboring channels to a diagonal array, suggesting a dimer formation. Fitting of an atomic model of isolated monomeric CatSper to the in situ map reveals supramolecular interactions and assembly of the CatSper complex. Loss of EFCAB9-CATSPERζ alters the architecture and interactions of the channels, resulting in fragmentation and misalignment of the zigzag-rows and disruption of flagellar movement in Efcab9 -/- sperm. This work offers unique insights into the structural basis for understanding CatSper regulation of sperm motility.

59 BASIC BIOLOGICAL SCIENCES↗

Quantification of Dissolved O 2 in Bulk Aqueous Solutions and Porous Media Using NMR Relaxometry

Effects of dissolved paramagnetic oxygen (O 2 ) in water on 1 H nuclear magnetic resonance (NMR) Carr-Purcell-Meiboom-Gill (CPMG) experiments is evaluated at a 1 H Larmor frequency of 2 MHz. Dissolution of O 2 into water significantly reduces the 1 H transverse relaxation coefficient (T 2 ). For deoxygenated water, T 2 is 3388 ms, water at ambient atmospheric conditions (7.4 mg/L O2) exhibits a T2 of 2465 ms, and dissolution of 2710 mg/L O 2 further reduces T 2 to 36 ms. The results were fit with an empirical model to facilitate prediction of T 2 times for bulk water as a function of paramagnetic oxygen concentrations in solution. Dissolved O 2 also greatly influences 1 H NMR CPMG experiments of confined water in a model system composed of Berea sandstone. For this system, 90 mg/L O 2 in H 2 O enhances T 2 relaxation of bulk water such that the relaxation time is comparable to physically confined water in the sandstone pores. Given the sensitivity of NMR T 2 coefficients to paramagnetic oxygen, low-field NMR-based characterization of fluid and porous media structure requires control of dissolved oxygen, as geospatial variation in the partial pressure of O 2 alone is expected to perturb fluid and pore relaxation times by up to 60 and 36%, respectively.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Topological structure of complex predictions

Abstract Current complex prediction models are the result of fitting deep neural networks, graph convolutional networks or transducers to a set of training data. A key challenge with these models is that they are highly parameterized, which makes describing and interpreting the prediction strategies difficult. We use topological data analysis to transform these complex prediction models into a simplified topological view of the prediction landscape. The result is a map of the predictions that enables inspection of the model results with more specificity than dimensionality-reduction methods such as tSNE and UMAP. The methods scale up to large datasets across different domains. We present a case study of a transformer-based model previously designed to predict expression levels of a piece of DNA in thousands of genomic tracks. When the model is used to study mutations in the BRCA1 gene, our topological analysis shows that it is sensitive to the location of a mutation and the exon structure of BRCA1 in ways that cannot be found with tools based on dimensionality reduction. Moreover, the topological framework offers multiple ways to inspect results, including an error estimate that is more accurate than model uncertainty. Further studies show how these ideas produce useful results in graph-based learning and image classification.

Computer Science↗

High-temperature 2D ferromagnetism in conjugated microporous porphyrin-type polymers

The need for magnetic 2D materials that are stable to the enviroment and have high Curie temperatures is very important for various electronic and spintronic applications. We have found that two-dimensional porphyrin-type aza-conjugated microporous polymer crystals are such a material (Fe-aza-CMPs). Fe-aza-CMPs are stable to CO, CO 2 , and O 2 atmospheres and show unusual adsorption, electronic, and magnetic properties. Indeed, they are semiconductors with small energy band gaps ranging from 0.27 eV to 0.626 eV. CO, CO 2 , and O 2 molecules can be attached in three different ways where single, double, or triple molecules are bound to iron atoms in Fe-aza-CMPs. For different attachment configurations we find that for CO and CO 2 a uniform distribution of the molecules is most energetically favorable while for O 2 molecules aggregation is most energetically preferable. The magnetic moments decrease from 4 to 2 to 0 for singly, doubly, triply occupied configurations for all gasses respectively. The most interesting magnetic properties are found for O 2 molecules attached to the Fe-aza-CMP. For a single attachment configuration we find that an antiferromagnetic state is favorable. When two O 2 molecules are attached, the calculations show the highest exchange integral with a value of J = 1071 μeV. This value has been verified by two independent methods where in the first method J is calculated by the energy difference between ferromagnetic and anitferromagnetic configurations. The second method is based on the frozen magnon approach where the magnon dispersion curve has been fitted by the Ising model. Finally, for the second method J has been estimated at J = 1100 μeV in excellent agreement with the first method.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Sensitivity measurements for a 250 MHz quartz shear-horizontal surface acoustic wave biosensor under liquid viscous loading

Surface acoustic wave (SAW) devices have been used in biochemical assays due to their high sensitivity. The device sensitivity is a function of changes in the density and viscosity of the liquid. Here, we studied the effect of fluid viscosity using a 250 MHz quartz shear-horizontal (SH)-SAW biosensor by monitoring different concentrations of binary aqueous/glycerol solutions. In this study, the sensitivity of the biosensor was determined by fitting the data to models derived from perturbation theory. Measurements in water were used as the reference. For a 0% to 50% glycerol solution, an 87°–204° separation in the phase shift was observed. The slope of the plot of the phase shift vs (ηρ)0.5 was used to indicate the sensor’s sensitivity. The sensitivity for our 250 MHz quartz SH-SAW sensors was calculated to be 3.7×10−3m2sKg. The corresponding mass sensitivity was determined to be 9.25 × 105m2Kg. The limit of detection was calculated to be 36 picograms (pg), while the limit of quantification or LOQ was calculated to be 109 pg. Traditionally, liquid phase measurements have been challenging for SAW devices because liquids dampen the vibrating sensors severely. This problem has been largely solved using a transverse (shear) wave instead of the more popular longitudinal or Rayleigh waves. Liquid measurements are now possible using transverse waves, also known as shear waves, because transverse waves are only minimally attenuated by liquids. Shear-horizontal SAW sensors (SH-SAW) show great promise as label-free biosensors because of their ability to handle liquid samples. However, the viscosity of the liquid still induces loading effects and can be measured when the liquid is loaded onto the SH-SAW propagating surface (delay line). When the liquid above the delay line is perturbed by physical or chemical changes, such as binding to a receptor, it alters the propagating acoustic wave. The SH-SAW device can measure these changes in liquid properties as a change in the wave’s phase compared to the original wave. The device’s phase shift was recorded as a function of the changes in the density and viscosity of the binary glycerol solution and used to determine the sensitivity in the linear dynamic range of responses.

Materials Science↗

The Atacama Cosmology Telescope: DR4 maps and cosmological parameters

In this paper, we present new arcminute-resolution maps of the Cosmic Microwave Background temperature and polarization anisotropy from the Atacama Cosmology Telescope, using data taken from 2013–2016 at 98 and 150 GHz. The maps cover more than 17,000 deg 2 , the deepest 600 deg 2 with noise levels below 10µK-arcmin. We use the power spectrum derived from almost 6,000 deg 2 of these maps to constrain cosmology. The ACT data enable a measurement of the angular scale of features in both the divergence-like polarization and the temperature anisotropy, tracing both the velocity and density at last-scattering. From these one can derive the distance to the last-scattering surface and thus infer the local expansion rate, H 0 . By combining ACT data with large-scale information from WMAP we measure H 0 = 67.6±1.1 km/s/Mpc, at 68% confidence, in excellent agreement with the independently-measured Planck satellite estimate (from ACT alone we find H 0 = 67.9 ± 1.5 km/s/Mpc). The ΛCDM model provides a good fit to the ACT data, and we find no evidence for deviations: both the spatial curvature, and the departure from the standard lensing signal in the spectrum, are zero to within 1σ; the number of relativistic species, the primordial Helium fraction, and the running of the spectral index are consistent with ΛCDM predictions to within 1.5–2.2σ. We compare ACT, WMAP, and Planck at the parameter level and find good consistency; we investigate how the constraints on the correlated spectral index and baryon density parameters readjust when adding CMB large-scale information that ACT does not measure. The DR4 products presented here will be publicly released on the NASA Legacy Archive for Microwave Background Data Analysis.

60 APPLIED LIFE SCIENCES↗

The microbial rare biosphere: current concepts, methods and ecological principles

ABSTRACT Our ability to describe the highly diverse pool of low abundance populations present in natural microbial communities is increasing at an unprecedented pace. Yet we currently lack an integrative view of the key taxa, functions and metabolic activity which make-up this communal pool, usually referred to as the ‘rare biosphere’, across the domains of life. In this context, this review examines the microbial rare biosphere in its broader sense, providing an historical perspective on representative studies which enabled to bridge the concept from macroecology to microbial ecology. It then addresses our current knowledge of the prokaryotic rare biosphere, and covers emerging insights into the ecology, taxonomy and evolution of low abundance microeukaryotic, viral and host-associated communities. We also review recent methodological advances and provide a synthetic overview on how the rare biosphere fits into different conceptual models used to explain microbial community assembly mechanisms, composition and function.

Pascoal, Francisco↗

A low-luminosity core-collapse supernova very similar to SN 2005cs

ABSTRACT We present observations and analysis of PSN J17292918+7542390, a low-luminosity Type II-P supernova (LL SN IIP). The observed sample of such events is still low, and their nature is still under debate. Such SNe are similar to SN 2005cs, a well-observed LL Type II-P event, having low expansion velocities, and small ejected 56Ni mass. We have developed a robust and relatively fast Monte Carlo code that fits semi-analytic models to light curves of core-collapse SNe. This allows the estimation of the most important physical parameters, like the radius of the progenitor star, the mass of the ejected envelope, the mass of the radioactive nickel synthesized during the explosion, among others. PSN J17292918+7542390 has $R_0 = 91_{-70}^{+119} \times 10^{11} \, \text{cm}$, $M_\text{ej} = 9.89_{-1.00}^{+2.10} \, \mathrm{ M}_{\odot }$, $E_{\mbox{kin}} = 0.65_{-0.18}^{+0.19} \, \text{foe}$, and $v_{\mbox{exp}} = 3332_{-347}^{+216}$ km s−1, for its progenitor radius, ejecta mass, kinetic energy, and expansion velocity, respectively. The initial nickel mass of the PSN J17292918+7542390 turned out to be $1.55_{-0.70}^{+0.75} \times 10^{-3} \, \mathrm{M}_{\odot }$. The measured photospheric velocity at the earliest observed phase is 7000 km s−1. As far as we can tell based on the small population of observed LL SNe IIP, the determined values are typical for these events.

Jäger, Jr, Zoltán↗

Discovery of a Candidate Binary Supermassive Black Hole in a Periodic Quasar from Circumbinary Accretion Variability

Binary supermassive black holes (BSBHs) are expected to be a generic byproduct from hierarchical galaxy formation. The final coalescence of BSBHs is thought to be the loudest gravitational wave (GW) siren, yet no confirmed BSBH is known in the GW-dominated regime. While periodic quasars have been proposed as BSBH candidates, the physical origin of the periodicity has been largely uncertain. Here we report discovery of a periodicity (P=1607±7 days) at 99.95% significance (with a global p-value of ~10 –3 accounting for the look elsewhere effect) in the optical light curves of a redshift 1.53 quasar, SDSS J025214.67–002813.7. Combining archival Sloan Digital Sky Survey data with new, sensitive imaging from the Dark Energy Survey, the total ~20-yr time baseline spans ~4.6 cycles of the observed 4.4-yr (restframe 1.7-yr) periodicity. The light curves are best fit by a bursty model predicted by hydrodynamic simulations of circumbinary accretion disks. The periodicity is likely caused by accretion rate modulation by a milli-parsec BSBH emitting GWs, dynamically coupled to the circumbinary accretion disk. A bursty hydrodynamic variability model is statistically preferred over a smooth, sinusoidal model expected from relativistic Doppler boost, a kinematic effect proposed for PG1302–102. Furthermore, the frequency dependence of the variability amplitudes disfavors Doppler boost, lending independent support to the circumbinary accretion variability hypothesis. Given our detection rate of one BSBH candidate from circumbinary accretion variability out of 625 quasars, it suggests that future large, sensitive synoptic surveys such as the Vera C. Rubin Observatory Legacy Survey of Space and Time may be able to detect hundreds to thousands of candidate BSBHs from circumbinary accretion with direct implications for Laser Interferometer Space Antenna.

79 ASTRONOMY AND ASTROPHYSICS↗

Blazar variability power spectra from radio up to TeV photon energies: Mrk 421 and PKS 2155−304

ABSTRACT We present the results of the power spectral density (PSD) analysis for the blazars Mrk 421 and PKS 2155−304, using good-quality, densely sampled light curves at multiple frequencies, covering 17 decades of the electromagnetic spectrum, and variability time-scales from weeks up to a decade. The data were collected from publicly available archives of observatories at radio from Owens Valley Radio Observatory, optical and infrared (B, V, R, I, J, H, and Kbands), X-rays from the Swift and the Rossi X-ray Timing Explorer, high and very high energy (VHE) γ-rays from the Fermi and Very Energetic Radiation Imaging Telescope Array System as well as the High Energy Stereoscopic System. Our results are: (1) the power-law form of the variability power spectra at radio, infrared, and optical frequencies have slopes ∼1.8, indicative of random-walk-type noise processes; (2) the power-law form of the variability power spectra at higher frequencies, from X-rays to VHE γ-rays, however, have slopes ∼1.2, suggesting a flicker noise-type process; and (3) there is significantly more variability power at X-rays, high and VHE γ-rays on time-scales ≲ 100 d, as compared to lower energies. Our results do not easily fit into a simple model, in which a single compact emission zone is dominating the radiative output of the blazars across all the time-scales probed in our analysis. Instead, we argue that the frequency-dependent shape of the variability power spectra points out a more complex picture, with highly inhomogeneous outflow producing non-thermal emission over an extended, stratified volume.

Goyal, Arti↗

Evidence for sub-Chandrasekhar Type Ia supernovae from the last major merger

ABSTRACT We investigate the contribution of sub-Chandrasekhar mass Type Ia supernovae to the chemical enrichment of the Gaia Sausage galaxy, the progenitor of a significant merger event in the early life of the Milky Way. Using a combination of data from Nissen & Schuster, the GALactic Archaeology with HERMES (GALAH) Data Release 3 [with 1D non-local thermal equilibrium (NLTE) abundance corrections], and the Apache Point Observatory Galactic Evolution Experiment (APOGEE) Data Release 16, we fit analytic chemical evolution models to a nine-dimensional chemical abundance space (Fe, Mg, Si, Ca, Cr, Mn, Ni, Cu, and Zn) in particular focusing on the iron-peak elements, Mn and Ni. We find that low [Mn/Fe] $\sim -0.15\, \mathrm{dex}$ and low [Ni/Fe] $\sim -0.3\, \mathrm{dex}$ Type Ia yields are required to explain the observed trends beyond the [α/Fe] knee of the Gaia Sausage (approximately at [Fe/H] $=-1.4\, \mathrm{dex}$). Comparison to theoretical yield calculations indicates a significant contribution from sub-Chandrasekhar mass Type Ia supernovae in this system (from ${\sim} 60\, \mathrm{per\, cent}$ to $100\, \mathrm{per\, cent}$ depending on the theoretical model with an additional ${\pm} 10\, \mathrm{per\, cent}$ systematic from NLTE corrections). We compare to results from other Local Group environments including dwarf spheroidal galaxies, the Magellanic Clouds, and the Milky Way’s bulge, finding the Type Ia [Mn/Fe] yield must be metallicity dependent. Our results suggest that sub-Chandrasekhar mass channels are a significant, perhaps even dominant, contribution to Type Ia supernovae in metal-poor systems, whilst more metal-rich systems could be explained by metallicity-dependent sub-Chandrasekhar mass yields, possibly with additional progenitor mass variation related to star formation history, or an increased contribution from Chandrasekhar mass channels at higher metallicity.

Sanders, Jason L.↗

The half-mass radius of MaNGA galaxies: effect of IMF gradients

Gradients in the stellar populations (SP) of galaxies – e.g. in age, metallicity, stellar initial mass function (IMF) – can result in gradients in the stellar-mass-to-light ratio, M * / L . Such gradients imply that the distribution of the stellar mass and light is different. For old SPs, e.g. in early-type galaxies at z ∼ 0, the M * / L gradients are weak if driven by variations in age and metallicity, but significantly larger if driven by the IMF. A gradient which has larger M * / L in the centre increases the estimated total stellar mass ( M * ) and reduces the scale which contains half this mass ( R e,* ), compared to when the gradient is ignored. For the IMF gradients inferred from fitting MILES simple SP models to the H β, 〈Fe〉, [MgFe], and TiO 2SDSS absorption lines measured in spatially resolved spectra of early-type galaxies in the MaNGA survey, the fractional change in R e,* can be significantly larger than that in M * , especially when the light is more centrally concentrated. The R e,* – M * correlation which results from accounting for IMF gradients is offset to smaller sizes by 0.3 dex compared to when these gradients are ignored. Comparisons with ‘quiescent’ galaxies at higher z must account for evolution in SP gradients (especially age and IMF) and in the light profile before drawing conclusions about how R e,* and M * evolve. Here, the implied merging between higher z and the present is less contrived if R e,* / R e at z ∼ 0 is closer to our IMF-driven gradient calibration than to unity.

79 ASTRONOMY AND ASTROPHYSICS↗

Deciphering enhancer sequence using thermodynamics-based models and convolutional neural networks

Abstract Deciphering the sequence-function relationship encoded in enhancers holds the key to interpreting non-coding variants and understanding mechanisms of transcriptomic variation. Several quantitative models exist for predicting enhancer function and underlying mechanisms; however, there has been no systematic comparison of these models characterizing their relative strengths and shortcomings. Here, we interrogated a rich data set of neuroectodermal enhancers in Drosophila, representing cis- and trans- sources of expression variation, with a suite of biophysical and machine learning models. We performed rigorous comparisons of thermodynamics-based models implementing different mechanisms of activation, repression and cooperativity. Moreover, we developed a convolutional neural network (CNN) model, called CoNSEPT, that learns enhancer ‘grammar’ in an unbiased manner. CoNSEPT is the first general-purpose CNN tool for predicting enhancer function in varying conditions, such as different cell types and experimental conditions, and we show that such complex models can suggest interpretable mechanisms. We found model-based evidence for mechanisms previously established for the studied system, including cooperative activation and short-range repression. The data also favored one hypothesized activation mechanism over another and suggested an intriguing role for a direct, distance-independent repression mechanism. Our modeling shows that while fundamentally different models can yield similar fits to data, they vary in their utility for mechanistic inference. CoNSEPT is freely available at: https://github.com/PayamDiba/CoNSEPT.

59 BASIC BIOLOGICAL SCIENCES↗

Accessing the gluon momentum fraction of nucleons through the gradient flow

We calculate the gluon momentum fraction of the nucleon using lattice QCD, with a nonperturbative renormalization technique based on the gradient flow. The gluon momentum fraction is determined on a single Wilson-clover ensemble using 𝑁 𝑓 =2 +1 flavors with pion mass 358 MeV and lattice spacing 0.094 fm. We employ the variational method to reduce excited-state contamination and apply the distillation framework to ensure a large operator basis. To reduce systematic uncertainties, we apply Bayesian model averaging to all fit procedures. We apply matching coefficients to the flow-time dependent lattice results to recover the gluon momentum fraction in the $\overline{MS}$-scheme at 2 GeV. Our final result is ⟨𝑥⟩ 𝑔 ⁢(𝜇 =2 GeV) =0.482⁢(35), where we quote only statistical uncertainties.

Lattice QCD↗

Charge Transport and Space-Charge Formation in Cd 1 - x Zn x Te 1 - y Se y Radiation Detectors

The electron- and hole-transport properties in cadmium zinc telluride selenide (CZTS) crystals are studied using a laser-induced transient-current technique with pulsed and dc bias. The internal electric field profile and velocity of surface recombination are determined by Monte Carlo simulations of electron and hole transient currents combined with a numerical solution of the drift-diffusion equation coupled with Poisson’s equation. Electron and hole drift mobilities of μe = 830 cm 2 /Vs and μh = 40 cm 2 /Vs, respectively, are determined. We also develop a simple technique for evaluating surface recombination directly from measured current waveforms without the need for numerical simulation. The good quality of the prepared detector at pulsed bias, with electron- and hole-mobility-lifetime products of (μτ)e = 1.9 × 10 -3 cm 2 /V and (μτ)h = 1.4 × 10 -4 cm 2 /V, respectively, are observed. The formation of a positive space charge, originating from hole injection combined with a recombination level, is found. We observe a significant position dependence of the lifetime of electrons and holes in dc bias due to hole injection. The experiment is successfully fitted by a simple model dominated by a single deep recombination level with an energy of E t =E C -0.73eV; concentration of 7.3 × 10 11 cm -3 ; and electron- and hole-capture cross sections of 3.5 × 10 -14 cm 2 and 6.5 × 10 -14 cm 2 , respectively.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Tuning the magnetic properties of the CrMnFeCoNi Cantor alloy

Magnetic properties of more than 20 Cantor alloy samples of varying composition were investigated over a temperature range of 5 K to 300 K and in fields of up to 70 kOe using magnetometry and muon spin relaxation. Two transitions are identified: a spin-glass-like transition that appears between 55K and 190K, depending on composition, and a ferrimagnetic transition that occurs at approximately 43K in multiple samples with widely varying compositions. The magnetic signatures at 43K are remarkably insensitive to chemical composition. A modified Curie-Weiss model was used to fit the susceptibility data and to extract the net effective magnetic moment for each sample. The resulting values for the net effective moment were either diminished with increasing Cr or Mn concentrations or enhanced with decreasing Fe, Co, or Ni concentrations. Beyond a sufficiently large effective moment, the magnetic ground state transitions from ferrimagnetism to ferromagnetism. The effective magnetic moments, together with the corresponding compositions, are used in a global linear regression analysis to extract element-specific effective magnetic moments, which are compared to the values obtained by ab initio based density functional theory calculations. Finally, these moments provide the information necessary to controllably tune the magnetic properties of Cantor alloy variants.

36 MATERIALS SCIENCE↗

Perturbation theory models for LSST-era galaxy clustering: Tests with subpercent mock catalog measurements in Fourier and configuration space

Here, we analyze the clustering of galaxies using the z = 1.006 snapshot of the CosmoDC2 simulation, a high-fidelity synthetic galaxy catalog designed to validate analysis methods for the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). We present sub-percent measurements of the galaxy auto-correlation and galaxy-dark matter cross correlations in Fourier space and configuration space for a magnitude-limited galaxy sample. At these levels of precision, the statistical errors of the measurement are comparable to the systematic effects present in the simulation and measurement procedure; nevertheless, using a hybrid-PT model, we are able to model non-linear galaxy bias with 0.5% precision up to scales of k max = 0.5 h/Mpc and r min = 4 Mpc/h. While the linear bias parameter is measured with 0.01% precision, other bias parameters are determined with considerably weaker constraints and sometimes bimodal posterior distributions. We compare our fiducial model with lower dimensional models where the higher-order bias parameters are fixed at their co-evolution values and find that leaving these parameters free provides significant improvements in our ability to model small scale information. We also compare bias parameters for galaxy samples defined using different magnitude bands and find agreement between samples containing equal numbers of galaxies. Finally, we compare bias parameters between Fourier space and configuration space and find moderate to significant tension between the two approaches. Although our model is often unable to fit the CosmoDC2 galaxy samples within the 0.1% precision of our measurements, our results suggest that the hybrid-PT model used in this analysis is capable of modeling non-linear galaxy bias within the percent level precision needed for upcoming galaxy surveys.

79 ASTRONOMY AND ASTROPHYSICS↗

Robustness of the Galactic Center excess morphology against masking

The Galactic Center excess (GCE) remains an enduring mystery, with leading explanations being annihilating dark matter or an unresolved population of millisecond pulsars. Analyzing the morphology of the GCE provides critical clues to identify its exact origin. We investigate the robustness of the inferred GCE morphology against the effects of masking, an important step in the analysis where the gamma-ray emission from point sources and the galactic disk are excluded. Using different masks constructed from Fermi point source catalogs and a wavelet method, we find that the GCE morphology, particularly its ellipticity and cuspiness, is relatively independent of the choice of mask for energies above 2–3 GeV. The GCE morphology systematically favors an approximately spherical shape, as expected for dark matter annihilation. Compared to various stellar bulge profiles, a spherical dark matter annihilation profile better fits the data across different masks and galactic diffuse emission backgrounds, except for the stellar bulge profile which provides a similar fit to the data. Modeling the GCE with two components, one from dark matter annihilation and one tracing the Coleman bulge, we find this two-component model outperforms any single component or combinations of dark matter annihilation and other stellar bulge profiles. Uncertainty remains about the exact fraction contributed by each component across different background models and masks. Furthermore, when the Coleman bulge dominates, its corresponding spectrum lacks characteristics typically associated with millisecond pulsars, suggesting that it mostly models the emission from other sources instead of the GCE that is still present and spherically symmetric.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗