Showing 1–50 of 87 results for author: Sinha, M

Exploring non-radial oscillation modes in dark matter admixed neutron stars

Authors: Pratik Thakur, Anil Kumar, Vivek Baruah Thapa, Vishal Parmar, Monika Sinha

Abstract: Because of their extreme densities and consequently, gravitational potential, compact objects such as neutron stars can prove to be excellent captors of dark matter particles. Considering purely gravitational interactions between dark and hadronic matter, we construct dark matter admixed stars composed of two-fluid matter subject to current astrophysical constraints of maximum mass and tidal defor… ▽ More Because of their extreme densities and consequently, gravitational potential, compact objects such as neutron stars can prove to be excellent captors of dark matter particles. Considering purely gravitational interactions between dark and hadronic matter, we construct dark matter admixed stars composed of two-fluid matter subject to current astrophysical constraints of maximum mass and tidal deformability. We choose a wide range of parameters to construct the dark matter equation of state, and the DDME2 parameterization for the hadronic equation of state. We then examine the effect of dark matter on the stellar structure, tidal deformability and non-radial modes considering the relativistic Cowling approximation. We find the effect on $p$-modes is substantial, with frequencies decreasing up to the typical $f-$mode frequency range for most stars with a dark matter halo. The effects on the $f-$mode frequency are less extreme. Finally, we find the most probable and $1σ$ values of the dark matter parameters used in this study. △ Less

Submitted 11 June, 2024; originally announced June 2024.

Comments: 23 pages and 5 figures

Identification and Quenching of Nugget Galaxies in the RESOLVE Survey at z = 0

Authors: Derrick S. Carr, Sheila J. Kannappan, Mark A. Norris, Manodeep Sinha, Michael L. Palumbo III, Kathleen D. Eckert, Amanda J. Moffett, Mugdha S. Polimera, Joel I. Bernstein, Zackary L. Hutchens

Abstract: We present a complete census of candidate nuggets, i.e., dense galaxies likely formed by compaction with intense gas influx, within the volume-limited $z \sim 0$ REsolved Spectroscopy Of a Local VolumE (RESOLVE) survey. These nuggets span all evolutionary stages and 3 orders of magnitude in stellar mass ($M_{*} \sim 10^{8} M_\odot$ to $10^{11} M_\odot$) from the dwarf to the giant regime. We devel… ▽ More We present a complete census of candidate nuggets, i.e., dense galaxies likely formed by compaction with intense gas influx, within the volume-limited $z \sim 0$ REsolved Spectroscopy Of a Local VolumE (RESOLVE) survey. These nuggets span all evolutionary stages and 3 orders of magnitude in stellar mass ($M_{*} \sim 10^{8} M_\odot$ to $10^{11} M_\odot$) from the dwarf to the giant regime. We develop selection criteria for our $z\sim0$ nugget candidates based on structure and introduce the use of environmental criteria to eliminate nugget-like objects with suspected non-compaction origins. The resulting $z\sim0$ nuggets follow expectations with respect to structure (i.e., density, size), population frequency, and likely origins. We show that the properties of our nugget census are consistent with permanent quenching above the gas-richness threshold scale (halo mass $M_{halo} \sim 10^{11.4} M_\odot$), cyclic temporary quenching below the threshold scale, and feedback from active galactic nuclei (AGN) assisting in permanent quenching. As predicted in simulations, most nuggets quench within the halo mass range $M_{halo} \sim 10^{11.45} M_\odot$ to $10^{11.9} M_\odot$. We find $\sim 0.29$ dex scatter around the star-forming main sequence for candidate blue nuggets below the threshold scale, which is consistent with temporary quenching as seen in simulations. A transitional population of green nuggets appears above the threshold scale. AGN also become more common in nuggets above this scale, and we see a likely AGN excess in nuggets vs. comparably selected non-nuggets. Our results provide the first observational confirmation of the mass-dependent, AGN-mediated shift from cyclic quenching to halo quenching in nuggets. △ Less

Submitted 6 June, 2024; originally announced June 2024.

Comments: 18 pages, 14 figures

Non-radial oscillations in newly born compact star considering effects of phase transition

Authors: Anil Kumar, Pratik Thakur, Monika Sinha

Abstract: The massive stars end their lives by supernova explosions leaving central compact objects that may evolve into neutron stars. Initially, after birth, the star remains hot and gradually cools down. We explore the matter and star properties during this initial stage of the compact stars considering the possibility of the appearance of deconfined quark matter in the core of the star. At the initial s… ▽ More The massive stars end their lives by supernova explosions leaving central compact objects that may evolve into neutron stars. Initially, after birth, the star remains hot and gradually cools down. We explore the matter and star properties during this initial stage of the compact stars considering the possibility of the appearance of deconfined quark matter in the core of the star. At the initial stage after the supernova explosion, the occurrence of non-radial oscillation in the newly born compact object is highly possible. Non-radial oscillations are an important source of GWs. There is a high chance for GWs from these oscillations, especially the nodeless fundamental (f-) mode to be detected by next-generation GW detectors. We study the evolution in frequencies of non-radial oscillation after birth considering phase transition and predicting the possible signature for different possibilities of theoretical compact star models. △ Less

Submitted 1 April, 2024; originally announced April 2024.

Comments: Accepted for publication in MNRAS

Dark Sage: Next-generation semi-analytic galaxy evolution with multidimensional structure and minimal free parameters

Authors: Adam R. H. Stevens, Manodeep Sinha, Alexander Rohl, Mawson W. Sammons, Boryana Hadzhiyska, César Hernández-Aguayo, Lars Hernquist

Abstract: After more than five years of development, we present a new version of Dark Sage, a semi-analytic model (SAM) of galaxy formation that breaks the mould for models of its kind. Included among the major changes is an overhauled treatment of stellar feedback that is derived from energy conservation, operates on local scales, affects gas gradually over time rather than instantaneously, and predicts a… ▽ More After more than five years of development, we present a new version of Dark Sage, a semi-analytic model (SAM) of galaxy formation that breaks the mould for models of its kind. Included among the major changes is an overhauled treatment of stellar feedback that is derived from energy conservation, operates on local scales, affects gas gradually over time rather than instantaneously, and predicts a mass-loading factor for every galaxy. Building on the model's resolved angular-momentum structure of galaxies, we now consider the heating of stellar discs, delivering predictions for disc structure both radially and vertically. We add a further dimension to stellar discs by tracking the distribution of stellar ages in each annulus. Each annulus--age bin has its own velocity dispersion and metallicity evolved in the model. This allows Dark Sage to make structural predictions for galaxies that previously only hydrodynamic simulations could. We present the model as run on the merger trees of the highest-resolution gravity-only simulation of the MillenniumTNG suite. Despite its additional complexity relative to other SAMs, Dark Sage only has three free parameters, the least of any SAM, which we calibrate exclusively against the cosmic star formation history and the $z=0$ stellar and HI mass functions using a particle-swarm optimisation method. The Dark Sage codebase, written in C and Python, is publicly available at https://github.com/arhstevens/DarkSage △ Less

Submitted 14 February, 2024; v1 submitted 7 December, 2023; originally announced December 2023.

Comments: Updates post referee's report. Get in touch for data access. 34 pages before references; 14 figures, 2 tables, 130 equations

Journal ref: PASA, 2024

Universal relations for compact stars with exotic degrees of freedom

Authors: Anil Kumar, Manoj Kumar Ghosh, Pratik Thakur, Vivek Baruah Thapa, Kamal Krishna Nath, Monika Sinha

Abstract: The nature of the highly dense matter inside the supernova remnant compact star is not constrained by terrestrial experiments and hence modeled phenomenologically to accommodate the astrophysical observations from compact stars. The observable properties of the compact stars are highly sensitive to the microscopic model of highly dense matter. However, some universal relations exist between some m… ▽ More The nature of the highly dense matter inside the supernova remnant compact star is not constrained by terrestrial experiments and hence modeled phenomenologically to accommodate the astrophysical observations from compact stars. The observable properties of the compact stars are highly sensitive to the microscopic model of highly dense matter. However, some universal relations exist between some macroscopic properties of compact stars independent of the matter model. We study the universal relation including the stars containing exotic degrees of freedom such as heavier strange and non-strange baryons, strange quark matter in normal and superconducting phases, etc. We examine the universal relations for quantities moment of inertia - tidal love number - quadrupole moment. We also study the correlation of non-radial f-mode and p-mode frequencies with stellar properties. We find the f-mode frequency observes the universal relation with dimensionless tidal deformability but the p-mode frequency does not show a good correlation with stellar properties. The p-mode frequency is sensitive to the composition of the matter. We find that universal relation is also applicable for stars with exotic matter in the core of the star with several models of exotic matter. △ Less

Submitted 27 March, 2024; v1 submitted 26 November, 2023; originally announced November 2023.

Comments: 15 pages, 5 figures

Fast neutron star cooling in light of the PREX-2 experiment

Authors: Trisha Sarkar, Vivek Baruah Thapa, Monika Sinha

Abstract: The nuclear symmetry energy and its behaviour with density has been recently evaluated with enhanced value by PREX-2 experiment. This new values enables direct Urca neutrino emission process to be functioning in the dense matter inside neutron stars. With this new outlook we study the cooling rate of canonical mass neutron stars and compare with available observational cooling data. We find most o… ▽ More The nuclear symmetry energy and its behaviour with density has been recently evaluated with enhanced value by PREX-2 experiment. This new values enables direct Urca neutrino emission process to be functioning in the dense matter inside neutron stars. With this new outlook we study the cooling rate of canonical mass neutron stars and compare with available observational cooling data. We find most of the isolated neutron star thermal profile is compatible with the cooling of canonical mass star including superfluidity suppression. △ Less

Submitted 31 August, 2023; originally announced August 2023.

Comments: 8 pages, 5 figures. Accepted for publication in Phys. Rev. C

Journal ref: Phys. Rev. C 108, 035801 (2023)

Thermal properties of the core of magnetar

Authors: Trisha Sarkar, Shalu Yadav, Monika Sinha

Abstract: During very early age of neutron stars, the core cools down faster compared to the crust creating a large thermal gradient in the interior of the star. During $10-100$ years, a cooling wave propagates from the core to the crust causing the interior of the star to thermalize. During this duration thermal properties of the core material is of great importance to understand the dynamics of the interi… ▽ More During very early age of neutron stars, the core cools down faster compared to the crust creating a large thermal gradient in the interior of the star. During $10-100$ years, a cooling wave propagates from the core to the crust causing the interior of the star to thermalize. During this duration thermal properties of the core material is of great importance to understand the dynamics of the interior of the star. The heat capacity and thermal conductivity of the core depends on the behaviour of matter inside the core. We investigate these two properties in case of magnetars. Due to presence of large magnetic field, the proton superconductivity is quenched partially inside the magnetars depending upon the comparative values of upper critical field and the strength of the magnetic field present. This produces non-uniformity in the behaviour of matter throughout the star. Moreover, such non-uniformity arises from the variation of nature of the pairing and values of the pairing gap energy. We find that the heat capacity is substantially reduced due to the presence of superfluidity. On the other hand, the thermal conductivity of neutron is enhanced due to proton superconductivity and gets reduced due to neutron superfluidity. Hence, the variation of the thermal properties due to superfluidity in presence of magnetic field is different at different radius inside the star. However, in all the cases the %minimum maximum variation is of the order one. This affects the thermal relaxation time of the star and eventually its the thermal evolution. △ Less

Submitted 14 June, 2023; originally announced June 2023.

Comments: 19 pages, 9 figures, 8 tables, Accepted in Journal of High Energy Astrophysics

Journal ref: Journal of High Energy Astrophysics 39(2023)27-42

Using the Star Forming Main Sequence To Explore Quiescent Galaxies Across Cosmic Time

Authors: Tyler Houston, Darren J. Croton, Manodeep Sinha

Abstract: In this letter, we explore the quiescent lives of central galaxies using the SAGE galaxy model and Uchuu dark matter simulation. We ask three questions: (1) How much of a galaxy's life is spent in quiescence? (2) How often do galaxies transition off the main sequence? (3) What is the typical duration of a quiescent phase? We find low and high-mass galaxies spend the highest fraction of their lives… ▽ More In this letter, we explore the quiescent lives of central galaxies using the SAGE galaxy model and Uchuu dark matter simulation. We ask three questions: (1) How much of a galaxy's life is spent in quiescence? (2) How often do galaxies transition off the main sequence? (3) What is the typical duration of a quiescent phase? We find low and high-mass galaxies spend the highest fraction of their lives in quiescence: 45 \pm 19% for log10(Mstar) < 9.0 (3.68 \pm 1.80 Gyr) and 26 \pm 25% for log10(Mstar) > 11.5 (3.46 \pm 3.30 Gyr), falling to 7 \pm 13% for galaxies in-between (0.82 \pm 1.57 Gyr). Low mass galaxies move in and out of quiescence frequently, 2.8 \pm 1.3 times on average, though only for short periods, 1.49 \pm 1.04 Gyr. This can be traced to the influence of supernova feedback on their quite stochastic evolution. Galaxies of higher mass have fewer quiescent periods, ~0.7 \pm 0.9, and their length increases with mass, peaking at 1.97 \pm 2.27 Gyr. However, our high-mass population comprises star-forming and quiescent galaxies with diverging evolutionary paths, so the actual time may be even longer. These high-mass trends are driven by radio mode feedback from supermassive black holes, which, once active, tend to remain active for extended periods. △ Less

Submitted 13 March, 2023; originally announced March 2023.

Comments: 5 pages, 3 figures, accepted to MNRAS Letters, open access version see https://doi.org/10.1093/mnrasl/slad031

Hybrid stars are compatible with recent astrophysical observations

Authors: Anil Kumar, Vivek Baruah Thapa, Monika Sinha

Abstract: Compact stars (CS) are stellar remnants of massive stars. Inside CSs the density is so high that matter is in subatomic form composed of nucleons. With increase of density of matter towards the centre of the objects other degrees of freedom like hyperons, heavier non-strange baryons, meson condensates may appear. Not only that at higher densities, the nucleons may get decomposed into quarks and fo… ▽ More Compact stars (CS) are stellar remnants of massive stars. Inside CSs the density is so high that matter is in subatomic form composed of nucleons. With increase of density of matter towards the centre of the objects other degrees of freedom like hyperons, heavier non-strange baryons, meson condensates may appear. Not only that at higher densities, the nucleons may get decomposed into quarks and form deconfined strange quark matter (SQM). If it is so then CSs may contain SQM in the core surrounded by nucleonic matter forming hybrid stars (HSs). However, the nature and composition of matter inside CSs can only be inferred from the astrophysical observations of these CSs. Recent astrophysical observations in terms of CS mass-radius (M-R) relation and gravitational wave (GW) observation indicate that the matter should be soft in the intermediate density range and stiff enough at higher density range to attain the maximum possible mass above $2~M_\odot$ which is not compatible with pure hadronic equation of states (EOSs). Consequently, we study the HS properties with different models of SQM and find that within vector bag model considering density dependent bag parameter, the model goes well with the astrophysical observations so far. △ Less

Submitted 11 March, 2023; originally announced March 2023.

Comments: 12 pages, 9 figures; Accepted for publication in Phys. Rev. D

Journal ref: Physical Review D, 107, 063024 (2023)

Tension between implications from PREX-2 data and gravitational tidal response on dense matter equation of state

Authors: Vivek Baruah Thapa, Monika Sinha

Abstract: Recently an improved value of neutron skin thickness of $^{208}\text{Pb}$ was reported in Lead Radius EXperiment-2 (PREX-2) to be $R_{\text{skin}}=R_n - R_p=(0.283\pm 0.071)$ fm which corresponds to high estimations of nuclear symmetry energy ($E_{\text{sym}}$) and its slope ($L_{\text{sym}}$). The updated values of $E_{\text{sym}}$ and $L_{\text{sym}}$ commensurating to the neutron star observabl… ▽ More Recently an improved value of neutron skin thickness of $^{208}\text{Pb}$ was reported in Lead Radius EXperiment-2 (PREX-2) to be $R_{\text{skin}}=R_n - R_p=(0.283\pm 0.071)$ fm which corresponds to high estimations of nuclear symmetry energy ($E_{\text{sym}}$) and its slope ($L_{\text{sym}}$). The updated values of $E_{\text{sym}}$ and $L_{\text{sym}}$ commensurating to the neutron star observable estimations lie exterior to the astrophysical observed range. The higher values of $L_{\text{sym}}$ at $n_0$ deduced from recent PREX-2 data correlates to matter being easily deformable (yielding higher radius values) around intermediate matter densities leading to higher values of $\tildeΛ$ creating a tension between the terrestrial and astrophysical observations. In this study, we exploit this tension to constrain the $Δ$-scalar meson coupling parameter space. △ Less

Submitted 15 February, 2023; originally announced February 2023.

Comments: 6 pages, 3 figures. Conference proceeding of the "Nuclear Physics in Astrophysics-X 2022". To appear in EPJ Web of Conferences

Journal ref: EPJ Web of Conferences 279, 10003 (2023)

Dust contribution to the panchromatic galaxy emission

Authors: Dian P. Triani, Darren J. Croton, Manodeep Sinha, Edward N. Taylor, Camilla Pacifici, Eli Dwek

Abstract: We have developed a pipeline called \mentari to generate the far-ultraviolet to far-infrared spectral energy distribution (SED) of galaxies from the \dustysage semi-analytic galaxy formation model (SAM). \dustysage incorporates dust-related processes directly on top of the basic ingredients of galaxy formation like gas infall, cooling, star formation, feedback, and mergers. We derive a physically… ▽ More We have developed a pipeline called \mentari to generate the far-ultraviolet to far-infrared spectral energy distribution (SED) of galaxies from the \dustysage semi-analytic galaxy formation model (SAM). \dustysage incorporates dust-related processes directly on top of the basic ingredients of galaxy formation like gas infall, cooling, star formation, feedback, and mergers. We derive a physically motivated attenuation model from the computed dust properties in \dustysage, so each galaxy has a self-consistent set of attenuation parameters based on the complicated dust physics that occurred across the galaxy's assembly history. Then, we explore several dust emission templates to produce infrared spectra. Our results show that a physically-motivated attenuation model is better for obtaining a consistent multi-wavelength description of galaxy formation and evolution, compared to using a constant attenuation. We compare our predictions with a compilation of observations and find that the fiducial model is in reasonable agreement with: (i) the observed $z=0$ luminosity functions from the far-ultraviolet to far-infrared simultaneously, and hence (ii) the local cosmic SED in the same range, (iii) the rest-frame K-band luminosity function across $0 < z < 3$, and (iv) the rest-frame far-ultraviolet luminosity function across $0 < z < 1$. Our model underproduces the far-ultraviolet emission at $z=2$ and $z=3$, which can be improved by altering the AGN feedback and dust processes in \dustysage. However, this combination thus worses the agreement at $z=0$, which suggests that more detailed treatment of such processes is required. △ Less

Submitted 12 December, 2022; originally announced December 2022.

Comments: 20 pages, 14 figures, accepted for publication in MNRAS

Toward Accurate Modeling of Galaxy Clustering on Small Scales: Halo Model Extensions and Lingering Tension

Authors: Gillian D. Beltz-Mohrmann, Adam O. Szewciw, Andreas A. Berlind, Manodeep Sinha

Abstract: This paper represents an effort to provide robust constraints on the galaxy-halo connection and simultaneously test the Planck LCDM cosmology using a fully numerical model of small-scale galaxy clustering. We explore two extensions to the standard Halo Occupation Distribution model: assembly bias, whereby halo occupation depends on both halo mass and the larger environment, and velocity bias, wher… ▽ More This paper represents an effort to provide robust constraints on the galaxy-halo connection and simultaneously test the Planck LCDM cosmology using a fully numerical model of small-scale galaxy clustering. We explore two extensions to the standard Halo Occupation Distribution model: assembly bias, whereby halo occupation depends on both halo mass and the larger environment, and velocity bias, whereby galaxy velocities do not perfectly trace the velocity of the dark matter within the halo. Moreover, we incorporate halo mass corrections to account for the impact of baryonic physics on the halo population. We identify an optimal set of clustering measurements to constrain this "decorated" HOD model for both low- and high-luminosity galaxies in SDSS DR7. We find that, for low-luminosity galaxies, a model with both assembly bias and velocity bias provides the best fit to the clustering measurements, with no tension remaining in the fit. In this model we find evidence for both central and satellite galaxy assembly bias at the 99% and 95% confidence levels, respectively. In addition, we find evidence for satellite galaxy velocity bias at the 99.9% confidence level. For high luminosity galaxies, we find no evidence for either assembly bias or velocity bias, but our model exhibits significant tension with SDSS measurements. We find that all of these conclusions still stand when we include the effects of baryonic physics on the halo mass function, suggesting that the tension we find for high luminosity galaxies may be due to a problem with our assumed cosmological model. △ Less

Submitted 22 May, 2023; v1 submitted 29 November, 2022; originally announced November 2022.

Comments: 27 pages, 13 figures

Journal ref: The Astrophysical Journal, 948(2):100 (2023)

(Anti)kaon condensation in strongly magnetized dense matter

Authors: Debraj Kundu, Vivek Baruah Thapa, Monika Sinha

Abstract: Recent observations of several massive pulsars, with masses near and above $2~M_\odot$, point towards the existence of matter at very high densities, compared to normal matter that we are familiar with in our terrestrial world. This leads to the possibility of appearance of exotic degrees of freedom other than nucleons inside the core of the neutrons stars (NS). Another significant property of NSs… ▽ More Recent observations of several massive pulsars, with masses near and above $2~M_\odot$, point towards the existence of matter at very high densities, compared to normal matter that we are familiar with in our terrestrial world. This leads to the possibility of appearance of exotic degrees of freedom other than nucleons inside the core of the neutrons stars (NS). Another significant property of NSs is the presence of high surface magnetic field, with highest range of the order of $\sim~10^{16}$ G. We study the properties of highly dense matter with the possibility of appearance of heavier strange and non-strange baryons, and kaons in presence of strong magnetic field. We find that the presence of a strong magnetic field stiffens the matter at high density, delaying the kaon appearance and, hence, increasing the maximum attainable mass of NS family. △ Less

Submitted 22 March, 2023; v1 submitted 26 October, 2022; originally announced October 2022.

Comments: 16 pages, 10 figures; accepted for publication in Phys. Rev. C

Journal ref: Physical Review C, 107(3), 035807 (2023)

Uchuu-$ν^2$GC galaxies and AGN: Cosmic variance forecasts of high-redshift AGN for JWST, Euclid, and LSST

Authors: Taira Oogi, Tomoaki Ishiyama, Francisco Prada, Manodeep Sinha, Darren Croton, Sofía A. Cora, Eric Jullo, Anatoly A. Klypin, Masahiro Nagashima, J. López Cacheiro, José Ruedas, Masakazu A. R. Kobayashi, Ryu Makiya

Abstract: Measurements of the luminosity function of active galactic nuclei (AGN) at high redshift ($z\gtrsim 6$) are expected to suffer from field-to-field variance, including cosmic and Poisson variances. Future surveys, such as those from the Euclid telescope and James Webb Space Telescope (JWST), will also be affected by field variance. We use the Uchuu simulation, a state-of-the-art cosmological $N$-bo… ▽ More Measurements of the luminosity function of active galactic nuclei (AGN) at high redshift ($z\gtrsim 6$) are expected to suffer from field-to-field variance, including cosmic and Poisson variances. Future surveys, such as those from the Euclid telescope and James Webb Space Telescope (JWST), will also be affected by field variance. We use the Uchuu simulation, a state-of-the-art cosmological $N$-body simulation with 2.1 trillion particles in a volume of $25.7~\mathrm{Gpc}^3$, combined with a semi-analytic galaxy and AGN formation model, to generate the Uchuu-$ν^2$GC catalog, publicly available, that allows us to investigate the field-to-field variance of the luminosity function of AGN. With this Uchuu-$ν^2$GC model, we quantify the cosmic variance as a function of survey area, AGN luminosity, and redshift. In general, cosmic variance decreases with increasing survey area and decreasing redshift. We find that at $z\sim6-7$, the cosmic variance depends weakly on AGN luminosity. This is because the typical mass of dark matter haloes in which AGN reside does not significantly depend on luminosity. Due to the rarity of AGN, Poisson variance dominates the total field-to-field variance, especially for bright AGN. We also examine the effect of parameters related to galaxy formation physics on the field variance. We discuss uncertainties present in the estimation of the faint-end of the AGN luminosity function from recent observations, and extend this to make predictions for the expected number of AGN and their variance for upcoming observations with Euclid, JWST, and the Legacy Survey of Space and Time (LSST). △ Less

Submitted 5 August, 2023; v1 submitted 29 July, 2022; originally announced July 2022.

Comments: 18 pages, 16 figures, accepted for publication in MNRAS, Figs. 3, 4, 11, 13, and 14 revised, Figs. 14 and 16 in the original manuscript moved to the Appendix, All Uchuu-$ν^2$GC data used in this paper are publicly available on our Skies & Universes website https://www.skiesanduniverses.org/Simulations/Uchuu/GalaxyCatalogues/

The Astropy Project: Sustaining and Growing a Community-oriented Open-source Project and the Latest Major Release (v5.0) of the Core Package

Authors: The Astropy Collaboration, Adrian M. Price-Whelan, Pey Lian Lim, Nicholas Earl, Nathaniel Starkman, Larry Bradley, David L. Shupe, Aarya A. Patil, Lia Corrales, C. E. Brasseur, Maximilian Nöthe, Axel Donath, Erik Tollerud, Brett M. Morris, Adam Ginsburg, Eero Vaher, Benjamin A. Weaver, James Tocknell, William Jamieson, Marten H. van Kerkwijk, Thomas P. Robitaille, Bruce Merry, Matteo Bachetti, H. Moritz Günther, Thomas L. Aldcroft , et al. (111 additional authors not shown)

Abstract: The Astropy Project supports and fosters the development of open-source and openly-developed Python packages that provide commonly needed functionality to the astronomical community. A key element of the Astropy Project is the core package $\texttt{astropy}$, which serves as the foundation for more specialized projects and packages. In this article, we summarize key features in the core package as… ▽ More The Astropy Project supports and fosters the development of open-source and openly-developed Python packages that provide commonly needed functionality to the astronomical community. A key element of the Astropy Project is the core package $\texttt{astropy}$, which serves as the foundation for more specialized projects and packages. In this article, we summarize key features in the core package as of the recent major release, version 5.0, and provide major updates for the Project. We then discuss supporting a broader ecosystem of interoperable packages, including connections with several astronomical observatories and missions. We also revisit the future outlook of the Astropy Project and the current status of Learn Astropy. We conclude by raising and discussing the current and future challenges facing the Project. △ Less

Submitted 28 June, 2022; originally announced June 2022.

Comments: 43 pages, 5 figures. To appear in ApJ. The author list has two parts: the authors that made significant contributions to the writing and/or coordination of the paper, followed by maintainers of and contributors to the Astropy Project. The position in the author list does not correspond to contributions to the Astropy Project as a whole

Compact star merger events with stars composed of interacting strange quark matter

Authors: Anil Kumar, Vivek Baruah Thapa, Monika Sinha

Abstract: We investigate the properties of stars participating in double compact star merger events considering interacting model of stable strange quark matter. We model the matter making it compatible with the recent astrophysical observations of compact star mass-radius and gravitational wave events. In this context we consider modified MIT bag model and vector bag model with and without self interaction… ▽ More We investigate the properties of stars participating in double compact star merger events considering interacting model of stable strange quark matter. We model the matter making it compatible with the recent astrophysical observations of compact star mass-radius and gravitational wave events. In this context we consider modified MIT bag model and vector bag model with and without self interaction. We find new upper bound on tidal deformability of $1.4~M_\odot$ strange star corresponding to the upper bound of effective tidal deformability inferred from gravitational wave event. Range of compactness of $1.4~M_\odot$ strange star is obtained as ${0.175}\leq{C_{1.4}}\leq{0.199}$. Radius range of $1.5M_\odot$ primary star is deduced to be ${10.57}\leq{R_{1.5}}\leq{12.04}$ km, following stringent GW170817 constraints. GW190425 constraints provide with upper limit on radius of $1.7$ solar mass strange star that it should be less than $13.41$ $\text{km}$. △ Less

Submitted 15 May, 2022; v1 submitted 23 April, 2022; originally announced April 2022.

Comments: 10 pages, 9 figures, 4 tables. Corrected some typos, published version

Journal ref: Monthly Notices of the Royal Astronomical Society, 513(3): 3788 (2022)

Direct URCA process in light of PREX-2

Authors: Vivek Baruah Thapa, Monika Sinha

Abstract: We study the implications of the recent development in nuclear symmetry energy constraints from PREX-2 data on dense matter equation of state and its impact on dURCA threshold density. In this work, we construct the equation of state within the framework of covariant density functional theory implementing coupling schemes of non-linear and density-dependent models and exploring the coupling parame… ▽ More We study the implications of the recent development in nuclear symmetry energy constraints from PREX-2 data on dense matter equation of state and its impact on dURCA threshold density. In this work, we construct the equation of state within the framework of covariant density functional theory implementing coupling schemes of non-linear and density-dependent models and exploring the coupling parameter space of isovector-vector meson to baryons constrained by the isospin asymmetry parameter values deduced from recent PREX-2 data. The modified parameter sets are applied to evaluate the dense matter properties. We find that the updated data suggests the occurrence of dURCA process within neutron star even with mass as low as one solar mass. △ Less

Submitted 12 February, 2022; originally announced March 2022.

Comments: 8 pages, 7 figures, 3 tables

Influence of the nuclear symmetry energy slope on observables of compact stars with $Δ$-admixed hypernuclear matter

Authors: Vivek Baruah Thapa, Monika Sinha

Abstract: In this work, we study the effects of nuclear symmetry energy slope on neutron star dense matter equation of state and its impact on neutron star observables (mass-radius, tidal response). We construct the equation of state within the framework of covariant density functional theory implementing coupling schemes of non-linear and density-dependent models with viability of heavier non-nucleonic deg… ▽ More In this work, we study the effects of nuclear symmetry energy slope on neutron star dense matter equation of state and its impact on neutron star observables (mass-radius, tidal response). We construct the equation of state within the framework of covariant density functional theory implementing coupling schemes of non-linear and density-dependent models with viability of heavier non-nucleonic degrees of freedom. The slope of symmetry energy parameter ($L_{\text{sym}}$) is adjusted following density-dependence of isovector meson coupling to baryons. We find that smaller values of $L_{\text{sym}}$ at saturation favour early appearance of $Δ$-resonances in comparison to hyperons leading to latter's threshold at higher matter densities. We also investigate the dependence of $L_{\text{sym}}$ on tidal deformability and compactness parameter of a $1.4~M_\odot$ neutron star for different equation of states and observe similar converging behaviour for larger $L_{\text{sym}}$ values. △ Less

Submitted 5 January, 2022; v1 submitted 23 December, 2021; originally announced December 2021.

Comments: 14 pages, 10 figures, 8 tables; corrected some typos; matches with published version

Journal ref: Phys. Rev. C 105, 015802 (2022)

Toward Accurate Modeling of Galaxy Clustering on Small Scales: Constraining the Galaxy-Halo Connection with Optimal Statistics

Authors: Adam O. Szewciw, Gillian D. Beltz-Mohrmann, Andreas A. Berlind, Manodeep Sinha

Abstract: Applying halo models to analyze the small-scale clustering of galaxies is a proven method for characterizing the connection between galaxies and their host halos. Such works are often plagued by systematic errors or are limited to clustering statistics which can be predicted analytically. In this work, we employ a numerical mock-based modeling procedure to examine the clustering of SDSS DR7 galaxi… ▽ More Applying halo models to analyze the small-scale clustering of galaxies is a proven method for characterizing the connection between galaxies and their host halos. Such works are often plagued by systematic errors or are limited to clustering statistics which can be predicted analytically. In this work, we employ a numerical mock-based modeling procedure to examine the clustering of SDSS DR7 galaxies. We apply a standard halo occupation distribution (HOD) model to dark-matter-only simulations with a LCDM cosmology. To constrain the theoretical models, we utilize a combination of galaxy number density and selected scales of the projected correlation function, redshift-space correlation function, group multiplicity function, average group velocity dispersion, mark correlation function, and counts-in-cells statistics. We design an algorithm to choose an optimal combination of measurements that yields tight and accurate constraints on our model parameters. Compared to previous work using fewer clustering statistics, we find significant improvement in the constraints on all parameters of our halo model for two different luminosity-threshold galaxy samples. Most interestingly, we obtain unprecedented high-precision constraints on the scatter in the relationship between galaxy luminosity and halo mass. However, our best-fit model results in significant tension (>4 sigma) for both samples, indicating the need to add second-order features to the standard HOD model. To guarantee the robustness of these results, we perform an extensive analysis of the systematic and statistical errors in our modeling procedure, including a first-of-its-kind study of the sensitivity of our constraints to changes in the halo mass function due to baryonic physics. △ Less

Submitted 7 October, 2021; originally announced October 2021.

Comments: 36 pages, 17 figures, submitted to ApJ

Baryonic dense matter in view of gravitational-wave observations

Authors: Vivek Baruah Thapa, Anil Kumar, Monika Sinha

Abstract: The detection of gravitational waves from the merger of binary neutron stars events (GW170817, GW190425) and subsequent estimations of tidal deformability play a key role in constraining the behaviour of dense matter. In addition, massive neutron star candidates ($\sim 2 M_{\odot}$) along with NICER mass-radius measurements also, set sturdy constraints on the dense matter equation of state. Strict… ▽ More The detection of gravitational waves from the merger of binary neutron stars events (GW170817, GW190425) and subsequent estimations of tidal deformability play a key role in constraining the behaviour of dense matter. In addition, massive neutron star candidates ($\sim 2 M_{\odot}$) along with NICER mass-radius measurements also, set sturdy constraints on the dense matter equation of state. Strict bounds from gravitational waves and massive neutron stars observations constrain the theoretical models of nuclear matter comportment at large density regimes. On the other hand, model parameters providing the highly dense matter response are bounded by nuclear saturation properties. This work analyses coupling parametrizations from two classes based on covariant density functional models: Non-Linear and Density-Dependent schemes. Considering these constraints together, we study possible models and parametrization schemes with the feasibility of exotic degrees of freedom in dense matter which go well with the astrophysical observations as well as the terrestrial laboratory experiments. We show that most parametrizations with non-linear schemes do not support the observations and experiments while density-dependent scheme goes well with both. Astrophysical observations are well explained if the inclusion of heavier non-strange baryons is considered as one fraction of the dense matter particle spectrum. △ Less

Submitted 4 September, 2021; v1 submitted 9 August, 2021; originally announced August 2021.

Comments: 13 pages, 13 figures, corrected some typos, matches the published version

Journal ref: Monthly Notices of the Royal Astronomical Society, 530, 1, 2024, 501-513

Exploring the relation between dust mass and galaxy properties using Dusty SAGE

Authors: Dian P. Triani, Manodeep Sinha, Darren J. Croton, Eli Dwek, Camilla Pacifici

Abstract: We explore the relation between dust and several fundamental properties of simulated galaxies using the Dusty SAGE semi-analytic model. In addition to tracing the standard galaxy properties, Dusty SAGE also tracks cold dust mass in the interstellar medium (ISM), hot dust mass in the halo and dust mass ejected by feedback activity. Based on their ISM dust content, we divide our galaxies into two ca… ▽ More We explore the relation between dust and several fundamental properties of simulated galaxies using the Dusty SAGE semi-analytic model. In addition to tracing the standard galaxy properties, Dusty SAGE also tracks cold dust mass in the interstellar medium (ISM), hot dust mass in the halo and dust mass ejected by feedback activity. Based on their ISM dust content, we divide our galaxies into two categories: ISM dust-poor and ISM dust-rich. We split the ISM dust-poor group into two subgroups: halo dust-rich and dust-poor (the latter contains galaxies that lack dust in both the ISM and halo). Halo dust-rich galaxies have high outflow rates of heated gas and dust and are more massive. We divide ISM dust-rich galaxies based on their specific star formation rate (sSFR) into star-forming and quenched subgroups. At redshift z=0, we find that ISM dust-rich galaxies have a relatively high sSFR, low bulge-to-total (BTT) mass ratio, and high gas metallicity. The high sSFR of ISM dust-rich galaxies allows them to produce dust in the stellar ejecta. Their metal-rich ISM enables dust growth via grain accretion. The opposite is seen in the ISM dust-poor group. Furthermore, ISM dust-rich galaxies are typically late-types, while ISM dust-poor galaxies resemble the early-type population, and we show how their ISM content evolves from being dust-rich to dust-poor. Finally, we investigate dust production from z=3 to z=0 and find that all groups evolve similarly, except for the quenched ISM dust-rich group. △ Less

Submitted 25 February, 2021; v1 submitted 24 February, 2021; originally announced February 2021.

Comments: 12 pages, 8 figure, accepted for publication in MNRAS

Massive $Δ$-resonance admixed hypernuclear stars with anti-kaon condensations

Authors: Vivek Baruah Thapa, Monika Sinha, Jia Jie Li, Armen Sedrakian

Abstract: In this work, we study the effect of (anti)kaon condensation on the properties of compact stars that develop hypernuclear cores with and without an admixture of $Δ$-resonances. We work within the covariant density functional theory with the parameters adjusted to $K$-atomic and kaon-nucleon scattering data in the kaonic sector. The density-dependent parameters in the hyperonic sector are adjusted… ▽ More In this work, we study the effect of (anti)kaon condensation on the properties of compact stars that develop hypernuclear cores with and without an admixture of $Δ$-resonances. We work within the covariant density functional theory with the parameters adjusted to $K$-atomic and kaon-nucleon scattering data in the kaonic sector. The density-dependent parameters in the hyperonic sector are adjusted to the data on $Λ$ and $Ξ^-$ hypernuclei data. The $Δ$-resonance couplings are tuned to the data obtained from their scattering off nuclei and heavy-ion collision experiments. We find that (anti)kaon condensate leads to a softening of the equation of state and lower maximum masses of compact stars than in the absence of the condensate. Both the $K^-$ and $\bar K^0$-condensations occur through a second-order phase transition, which implies no mixed-phase formation. For large values of (anti)kaon and $Δ$-resonance potentials in symmetric nuclear matter, we observe that condensation leads to an extinction of $ Ξ^{-,0}$ hyperons. We also investigate the influence of inclusion of additional hidden-strangeness $σ^{*}$ meson in the functional and find that it leads to a substantial softening of the equation of state and delay in the onset of (anti)kaons. △ Less

Submitted 17 February, 2021; originally announced February 2021.

Comments: 14 pages, 11 figures, accepted for publication in Phys. Rev. D

Journal ref: Phys. Rev. D 103, 063004 (2021)

Void Galaxies Follow a Distinct Evolutionary Path in the Environmental COntext Catalog

Authors: Jonathan Florez, Andreas A. Berlind, Sheila J. Kannappan, David V. Stark, Kathleen D. Eckert, Victor F. Calderon, Amanda J. Moffett, Duncan Campbell, Manodeep Sinha

Abstract: We measure the environmental dependence, where environment is defined by the distance to the third nearest neighbor, of multiple galaxy properties inside the Environmental COntext (ECO) catalog. We focus primarily on void galaxies, which we define as the $10 \%$ of galaxies having the lowest local density. We compare the properties of void and non-void galaxies: baryonic mass, color, fractional st… ▽ More We measure the environmental dependence, where environment is defined by the distance to the third nearest neighbor, of multiple galaxy properties inside the Environmental COntext (ECO) catalog. We focus primarily on void galaxies, which we define as the $10 \%$ of galaxies having the lowest local density. We compare the properties of void and non-void galaxies: baryonic mass, color, fractional stellar mass growth rate (FSMGR), morphology, and gas-to-stellar-mass ratio (estimated from a combination of HI data and photometric gas fractions calibrated with the RESOLVE survey). Our void galaxies typically have lower baryonic masses than galaxies in denser environments, and they display the properties expected of a lower mass population: they have more late-types, are bluer, have higher FSMGR, and are more gas rich. We control for baryonic mass and investigate the extent to which void galaxies are different at fixed mass. Void galaxies are bluer, more gas-rich, and more star forming at fixed mass than non-void galaxies, which is a possible signature of galaxy assembly bias. Furthermore, we show that these trends persist even at fixed mass and morphology, and we find that voids host a distinct population of early-types that are bluer and more star-forming than the typical red and quenched early-types. In addition to these empirical observational results, we also present theoretical results from mock catalogs with built-in galaxy assembly bias. We show that a simple matching of galaxy properties to (sub)halo properties, such as mass and age, can recover the observed environmental trends in ECO galaxies. △ Less

Submitted 16 November, 2020; originally announced November 2020.

Comments: 18 pages, 11 figures. Accepted for publication in ApJ

Dense matter equation of state of a massive neutron star with anti-kaon condensation

Authors: Vivek Baruah Thapa, Monika Sinha

Abstract: Recent measurements of neutron star mass from several candidates (PSR J$1614-2230$, PSR J$0348+0432$, MSP J$0740+6620$) set the lower bound on the maximum possible mass for this class of compact objects $\sim 2$ M$_\odot$. Existence of stars with high mass brings the possibility of existence of exotic matter (hyperons, meson condensates) at the core region of the objects. In this work, we investig… ▽ More Recent measurements of neutron star mass from several candidates (PSR J$1614-2230$, PSR J$0348+0432$, MSP J$0740+6620$) set the lower bound on the maximum possible mass for this class of compact objects $\sim 2$ M$_\odot$. Existence of stars with high mass brings the possibility of existence of exotic matter (hyperons, meson condensates) at the core region of the objects. In this work, we investigate the (anti)kaon ($K^-, \bar{K}^0$) condensation in $β-$equilibrated nuclear matter within the framework of covariant density functional theory. The functionals in the kaonic sector are constrained by the experimental studies on $K^-$ atomic, kaon-nucleon scattering data fits. We find that the equation of state softens with the inclusion of (anti)kaon condensates, which lowers the maximum mass of neutron star. In one of the density-independent coupling cases, the $K^-$ condensation is through a first-order phase transition type, which produces a $2$ M$_\odot$ neutron star. The first-order phase transition results in mixed phase region in the inner core of the stars. While $\bar{K}^0$ condensation appears via second-order phase transition for all the models we consider here. △ Less

Submitted 3 December, 2020; v1 submitted 12 November, 2020; originally announced November 2020.

Comments: 14 pages, 13 figures, corrected some typos, matches with the published version

Journal ref: Phys. Rev. D 102 (2020), 123007

Ambipolar diffusion velocity and magnetic field evolution in magnetar core: Generalised theoretical approach

Authors: Monika Sinha, Manoj Ghosh

Abstract: The magnetic field associated with neutron stars is generally believed to be threaded inside the star. In the presence of a magnetic field, the plasma present in the interior of the star goes through several processes that lead to magnetic field evolution. It is thought that magnetar activities are mainly due to field decay. The most important process of field decay inside the core of the star is… ▽ More The magnetic field associated with neutron stars is generally believed to be threaded inside the star. In the presence of a magnetic field, the plasma present in the interior of the star goes through several processes that lead to magnetic field evolution. It is thought that magnetar activities are mainly due to field decay. The most important process of field decay inside the core of the star is the ambipolar diffusion of the charged particles present in the interior plasma. The decay rate due to ambipolar diffusion is directly connected to the ambipolar velocity of the charged particles under the influence of the present magnetic field. The ambipolar velocity of the charged particles depends on the internal dynamics of the particles. We outline a general method to solve the particle dynamics in the presence of a magnetic field using a magnetohydrodynamic equation for ambipolar velocity. The equation is general and applies to all possible surrounding conditions \eg temperature, and matter states like normal or superfluid. △ Less

Submitted 29 September, 2023; v1 submitted 21 October, 2020; originally announced October 2020.

Comments: 15 pages, Calculation modified extensively, Appendix added

Equation of state of strongly magnetized matter with hyperons and $Δ$-resonances

Authors: Vivek Baruah Thapa, Monika Sinha, Jia-Jie Li, Armen Sedrakian

Abstract: We construct a new equation of state for the baryonic matter under an intense magnetic field within the framework of covariant density functional theory. The composition of matter includes hyperons as well as $ Δ$-resonances. The extension of the nucleonic functional to the hypernuclear sector is constrained by the experimental data on $Λ$ and $Ξ$-hypernuclei. We find that the equation of state st… ▽ More We construct a new equation of state for the baryonic matter under an intense magnetic field within the framework of covariant density functional theory. The composition of matter includes hyperons as well as $ Δ$-resonances. The extension of the nucleonic functional to the hypernuclear sector is constrained by the experimental data on $Λ$ and $Ξ$-hypernuclei. We find that the equation of state stiffens with the inclusion of the magnetic field, which increases the maximum mass of neutron star compared to the non-magnetic case. In addition, the strangeness fraction in the matter is enhanced. Several observables, like the Dirac effective mass, particle abundances, etc show typical oscillatory behavior as a function of the magnetic field and/or density which is traced back to the occupation pattern of Landau levels. △ Less

Submitted 16 October, 2020; v1 submitted 2 October, 2020; originally announced October 2020.

Comments: 16 pages, 8 figures, Accepted in journal "Particles", Corrected some typos. Matches with published version

Journal ref: Particles 3(4), 660-675 (2020)

Ambipolar decay of magnetic field in magnetars and the observed magnetar activities

Authors: Badal Bhalla, Monika Sinha

Abstract: Magnetars are comparatively young neutron stars with ultra-strong surface magnetic field in the range $10^{14-16}$ G. The old neutron stars have surface magnetic field some what less $\sim 10^8$ G which clearly indicates the decay of field with time. One possible way of magnetic field decay is by ambipolar diffusion. We describe the general procedure to solve for the ambipolar velocity inside the… ▽ More Magnetars are comparatively young neutron stars with ultra-strong surface magnetic field in the range $10^{14-16}$ G. The old neutron stars have surface magnetic field some what less $\sim 10^8$ G which clearly indicates the decay of field with time. One possible way of magnetic field decay is by ambipolar diffusion. We describe the general procedure to solve for the ambipolar velocity inside the star core without any approximation. With a realistic model of neutron star we determine the ambipolar velocity configuration inside the neutron star core and hence find the ambipolar decay rate and time scale which is consistent with the magnetar observations. △ Less

Submitted 21 June, 2021; v1 submitted 21 September, 2020; originally announced September 2020.

Comments: 12 pages, 4 figures, Accepted for publication in MPLA

Journal ref: Vol. 36, No. 20, 2150144 (2021)

Machine Learning the Fates of Dark Matter Subhalos: A Fuzzy Crystal Ball

Authors: Abigail Petulante, Andreas A. Berlind, J. Kelly Holley-Bockelmann, Manodeep Sinha

Abstract: The evolution of a dark matter halo in a dark matter only simulation is governed purely byNewtonian gravity, making a clean testbed to determine what halo properties drive its fate.Using machine learning, we predict the survival, mass loss, final position, and merging time of subhalos within a cosmological N-body simulation, focusing on what instantaneous initial features of the halo, interaction,… ▽ More The evolution of a dark matter halo in a dark matter only simulation is governed purely byNewtonian gravity, making a clean testbed to determine what halo properties drive its fate.Using machine learning, we predict the survival, mass loss, final position, and merging time of subhalos within a cosmological N-body simulation, focusing on what instantaneous initial features of the halo, interaction, and environment matter most. Survival is well predicted, with our model achieving 96.5% accuracy using only 3 model inputs from the initial interaction.However, the mass loss, final location, and merging times are much more stochastic processes, with significant margins of error between the true and predicted quantities for much of our sample. The redshift, impact angle, relative velocity, and the masses of the host and subhalo are the only relevant initial inputs for determining subhalo evolution. In general, subhalos that enter their hosts at a mid-range of redshifts (typically z = 0.67-0.43) are the most challenging to make predictions for, across all of our final outcomes. Subhalo orbits that come in more perpendicular to the host are also easier to predict, except for in the case of predicting disruption, where the opposite appears to be true. We conclude that the detailed evolution of individual subhalos within N-body simulations is quite difficult to predict, pointing to a stochasticity in the merging process. We discuss implications for both simulations and observations △ Less

Submitted 11 August, 2020; originally announced August 2020.

Comments: 19 pages, 11 figures

The Uchuu Simulations: Data Release 1 and Dark Matter Halo Concentrations

Authors: Tomoaki Ishiyama, Francisco Prada, Anatoly A. Klypin, Manodeep Sinha, R. Benton Metcalf, Eric Jullo, Bruno Altieri, Sofía A. Cora, Darren Croton, Sylvain de la Torre, David E. Millán-Calero, Taira Oogi, José Ruedas, Cristian A. Vega-Martínez

Abstract: We introduce the Uchuu suite of large high-resolution cosmological $N$-body simulations. The largest simulation, named Uchuu, consists of 2.1 trillion ($12800^3$) dark matter particles in a box of side-length 2.0 Gpc/h, with particle mass $3.27 \times 10^{8}$ Msun/h. The highest resolution simulation, Shin-Uchuu, consists of 262 billion ($6400^3$) particles in a box of side-length 140 Mpc/h, with… ▽ More We introduce the Uchuu suite of large high-resolution cosmological $N$-body simulations. The largest simulation, named Uchuu, consists of 2.1 trillion ($12800^3$) dark matter particles in a box of side-length 2.0 Gpc/h, with particle mass $3.27 \times 10^{8}$ Msun/h. The highest resolution simulation, Shin-Uchuu, consists of 262 billion ($6400^3$) particles in a box of side-length 140 Mpc/h, with particle mass $8.97 \times 10^{5}$ Msun/h. Combining these simulations we can follow the evolution of dark matter halos and subhalos spanning those hosting dwarf galaxies to massive galaxy clusters across an unprecedented volume. In this first paper, we present basic statistics, dark matter power spectra, and the halo and subhalo mass functions, which demonstrate the wide dynamic range and superb statistics of the Uchuu suite. From an analysis of the evolution of the power spectra we conclude that our simulations remain accurate from the Baryon Acoustic Oscillation scale down to the very small. We also provide parameters of a mass-concentration model, which describes the evolution of halo concentration and reproduces our simulation data to within 5 per cent for halos with masses spanning nearly eight orders of magnitude at redshift 0<z<14. There is an upturn in the mass-concentration relation for the population of all halos and of relaxed halos at z>0.5, whereas no upturn is detected at z<0.5. We make publicly available various $N$-body products as part of Uchuu Data Release 1 on the Skies & Universes site. Future releases will include gravitational lensing maps and mock galaxy, X-ray cluster, and active galactic nuclei catalogues. △ Less

Submitted 13 July, 2021; v1 submitted 29 July, 2020; originally announced July 2020.

Comments: 22 pages, 14 figures, accepted by MNRAS. We release various $N$-body products as data release 1 on http://skiesanduniverses.org/ such as subsets of simulation particles, matter power spectra, halo/subhalo catalogues, and their merger trees

Mentari: A pipeline to model the galaxy SED using semi analytic models

Authors: Dian Triani, Darren Croton, Manodeep Sinha

Abstract: We build a theoretical picture of how the light from galaxies evolves across cosmic time. In particular, we predict the evolution of the galaxy spectral energy distribution (SED) by carefully integrating the star formation and metal enrichment histories of semi-analytic model (SAM) galaxies and combining these with stellar population synthesis models which we call mentari. Our SAM combines prescri… ▽ More We build a theoretical picture of how the light from galaxies evolves across cosmic time. In particular, we predict the evolution of the galaxy spectral energy distribution (SED) by carefully integrating the star formation and metal enrichment histories of semi-analytic model (SAM) galaxies and combining these with stellar population synthesis models which we call mentari. Our SAM combines prescriptions to model the interplay between gas accretion, star formation, feedback process, and chemical enrichment in galaxy evolution. From this, the SED of any simulated galaxy at any point in its history can be constructed and compared with telescope data to reverse engineer the various physical processes that may have led to a particular set of observations. The synthetic SEDs of millions of simulated galaxies from mentari can cover wavelengths from the far UV to infrared, and thus can tell a near complete story of the history of galaxy evolution. \keywords{galaxies: evolution - galaxies: stellar content - galaxies.} △ Less

Submitted 16 June, 2020; originally announced June 2020.

Comments: 5 pages, 3 figures, Proceedings of IAU symposium #341: Challenges in Panchromatic Modelling with Next Generation Facilities

The origin of dust in galaxies across cosmic time

Authors: Dian P. Triani, Manodeep Sinha, Darren J. Croton, Camilla Pacifici, Eli Dwek

Abstract: We study the dust evolution in galaxies by implementing a detailed dust prescription in the SAGE semi-analytical model for galaxy formation. The new model, called Dusty SAGE, follows the condensation of dust in the ejecta of type II supernovae and asymptotic giant branch (AGB) stars, grain growth in the dense molecular clouds, destruction by supernovae shocks, and the removal of dust from the ISM… ▽ More We study the dust evolution in galaxies by implementing a detailed dust prescription in the SAGE semi-analytical model for galaxy formation. The new model, called Dusty SAGE, follows the condensation of dust in the ejecta of type II supernovae and asymptotic giant branch (AGB) stars, grain growth in the dense molecular clouds, destruction by supernovae shocks, and the removal of dust from the ISM by star formation, reheating, inflows and outflows. Our model successfully reproduces the observed dust mass function at redshift z = 0 and the observed scaling relations for dust across a wide range of redshifts. We find that the dust mass content in the present Universe is mainly produced via grain growth in the interstellar medium (ISM). By contrast, in the early Universe, the primary production mechanism for dust is the condensation in stellar ejecta. The shift of the significant production channel for dust characterises the scaling relations of dust-to-gas (DTG) and dust-to-metal (DTM) ratios. In galaxies where the grain growth dominates, we find positive correlations for DTG and DTM ratios with both metallicity and stellar mass. On the other hand, in galaxies where dust is produced primarily via condensation, we find negative or no correlation for DTM and DTG ratios with either metallicity or stellar mass. In agreement with observation showing that the circumgalactic medium (CGM) contains more dust than the ISM, our model also shows the same trend for z < 4. Our semi-analytic model is publicly available at https: //github.com/dptriani/dusty-sage. △ Less

Submitted 12 February, 2020; originally announced February 2020.

Comments: 19 pages, 14 figures, accepted for publication in MNRAS

Journal ref: 2020MNRAS.tmp..415T

Corrfunc: Blazing fast correlation functions with AVX512F SIMD Intrinsics

Authors: Manodeep Sinha, Lehman H. Garrison

Abstract: Correlation functions are widely used in extra-galactic astrophysics to extract insights into how galaxies occupy dark matter halos and in cosmology to place stringent constraints on cosmological parameters. A correlation function fundamentally requires computing pair-wise separations between two sets of points and then computing a histogram of the separations. Corrfunc is an existing open-source,… ▽ More Correlation functions are widely used in extra-galactic astrophysics to extract insights into how galaxies occupy dark matter halos and in cosmology to place stringent constraints on cosmological parameters. A correlation function fundamentally requires computing pair-wise separations between two sets of points and then computing a histogram of the separations. Corrfunc is an existing open-source, high-performance software package for efficiently computing a multitude of correlation functions. In this paper, we will discuss the SIMD AVX512F kernels within Corrfunc, capable of processing 16 floats or 8 doubles at a time. The latest manually implemented Corrfunc AVX512F kernels show a speedup of up to $\sim 4\times$ relative to compiler-generated code for double-precision calculations. The AVX512F kernels show $\sim 1.6\times$ speedup relative to the AVX kernels and compare favorably to a theoretical maximum of $2\times$. In addition, by pruning pairs with too large of a minimum possible separation, we achieve a $\sim 5-10\%$ speedup across all the SIMD kernels. Such speedups highlight the importance of programming explicitly with SIMD vector intrinsics for complex calculations that can not be efficiently vectorized by compilers. Corrfunc is publicly available at https://github.com/manodeep/Corrfunc/. △ Less

Submitted 15 November, 2019; originally announced November 2019.

Comments: Paper II for the Corrfunc software package, paper I is on arXiv here: arXiv:1911.03545. Appeared in the refereed proceedings for the "Second Workshop on Software Challenges to Exascale Computing"

emcee v3: A Python ensemble sampling toolkit for affine-invariant MCMC

Authors: Daniel Foreman-Mackey, Will M. Farr, Manodeep Sinha, Anne M. Archibald, David W. Hogg, Jeremy S. Sanders, Joe Zuntz, Peter K. G. Williams, Andrew R. J. Nelson, Miguel de Val-Borro, Tobias Erhardt, Ilya Pashchenko, Oriol Abril Pla

Abstract: emcee is a Python library implementing a class of affine-invariant ensemble samplers for Markov chain Monte Carlo (MCMC). This package has been widely applied to probabilistic modeling problems in astrophysics where it was originally published, with some applications in other fields. When it was first released in 2012, the interface implemented in emcee was fundamentally different from the MCMC li… ▽ More emcee is a Python library implementing a class of affine-invariant ensemble samplers for Markov chain Monte Carlo (MCMC). This package has been widely applied to probabilistic modeling problems in astrophysics where it was originally published, with some applications in other fields. When it was first released in 2012, the interface implemented in emcee was fundamentally different from the MCMC libraries that were popular at the time, such as PyMC, because it was specifically designed to work with "black box" models instead of structured graphical models. This has been a popular interface for applications in astrophysics because it is often non-trivial to implement realistic physics within the modeling frameworks required by other libraries. Since emcee's release, other libraries have been developed with similar interfaces, such as dynesty (Speagle 2019). The version 3.0 release of emcee is the first major release of the library in about 6 years and it includes a full re-write of the computational backend, several commonly requested features, and a set of new "move" implementations. △ Less

Submitted 18 November, 2019; originally announced November 2019.

Comments: Published in the Journal for Open Source Software

Journal ref: Journal of Open Source Software, 2 4(43), 1864 (2019)

Corrfunc --- A Suite of Blazing Fast Correlation Functions on the CPU

Authors: Manodeep Sinha, Lehman H. Garrison

Abstract: The two-point correlation function (2PCF) is the most widely used tool for quantifying the spatial distribution of galaxies. Since the distribution of galaxies is determined by galaxy formation physics as well as the underlying cosmology, fitting an observed correlation function yields valuable insights into both. The calculation for a 2PCF involves computing pair-wise separations and consequently… ▽ More The two-point correlation function (2PCF) is the most widely used tool for quantifying the spatial distribution of galaxies. Since the distribution of galaxies is determined by galaxy formation physics as well as the underlying cosmology, fitting an observed correlation function yields valuable insights into both. The calculation for a 2PCF involves computing pair-wise separations and consequently, the computing time scales quadratically with the number of galaxies. The next-generation galaxy surveys are slated to observe many millions of galaxies, and computing the 2PCF for such surveys would be prohibitively time-consuming. Additionally, modern modelling techniques require the 2PCF to be calculated thousands of times on simulated galaxy catalogues of {\em at least} equal size to the data and would be completely unfeasible for the next generation surveys. Thus, calculating the 2PCF forms a substantial bottleneck in improving our understanding of the fundamental physics of the universe, and we need high-performance software to compute the correlation function. In this paper, we present Corrfunc --- a suite of highly optimised, OpenMP parallel clustering codes. The improved performance of Corrfunc arises from both efficient algorithms as well as software design that suits the underlying hardware of modern CPUs. Corrfunc can compute a wide range of 2-D and 3-D correlation functions in either simulation (Cartesian) space or on-sky coordinates. Corrfunc runs efficiently in both single- and multi-threaded modes and can compute a typical 2-point projected correlation function ($w_p(r_p)$) for ~1 million galaxies within a few seconds on a single thread. Corrfunc is designed to be both user-friendly and fast and is publicly available at https://github.com/manodeep/Corrfunc. △ Less

Submitted 8 November, 2019; originally announced November 2019.

Comments: Accepted for publication to MNRAS

The Growing Importance of a Tech Savvy Astronomy and Astrophysics Workforce

Authors: Dara Norman, Kelle Cruz, Vandana Desai, Britt Lundgren, Eric Bellm, Frossie Economou, Arfon Smith, Amanda Bauer, Brian Nord, Chad Schafer, Gautham Narayan, Ting Li, Erik Tollerud, Brigitta Sipocz, Heloise Stevance, Timothy Pickering, Manodeep Sinha, Joseph Harrington, Jeyhan Kartaltepe, Dany Vohl, Adrian Price-Whelan, Brian Cherinka, Chi-kwan Chan, Benjamin Weiner, Maryam Modjaz , et al. (4 additional authors not shown)

Abstract: Fundamental coding and software development skills are increasingly necessary for success in nearly every aspect of astronomical and astrophysical research as large surveys and high resolution simulations become the norm. However, professional training in these skills is inaccessible or impractical for many members of our community. Students and professionals alike have been expected to acquire th… ▽ More Fundamental coding and software development skills are increasingly necessary for success in nearly every aspect of astronomical and astrophysical research as large surveys and high resolution simulations become the norm. However, professional training in these skills is inaccessible or impractical for many members of our community. Students and professionals alike have been expected to acquire these skills on their own, apart from formal classroom curriculum or on-the-job training. Despite the recognized importance of these skills, there is little opportunity to develop them - even for interested researchers. To ensure a workforce capable of taking advantage of the computational resources and the large volumes of data coming in the next decade, we must identify and support ways to make software development training widely accessible to community members, regardless of affiliation or career level. To develop and sustain a technology capable astronomical and astrophysical workforce, we recommend that agencies make funding and other resources available in order to encourage, support and, in some cases, require progress on necessary training, infrastructure and policies. In this white paper, we focus on recommendations for how funding agencies can lead in the promotion of activities to support the astronomy and astrophysical workforce in the 2020s. △ Less

Submitted 17 October, 2019; originally announced October 2019.

Comments: Submitted as a ASTRO2020 Decadal Survey APC position paper. arXiv admin note: substantial text overlap with arXiv:1905.05116

Astro2020 APC White Paper: Elevating the Role of Software as a Product of the Research Enterprise

Authors: Arfon M. Smith, Dara Norman, Kelle Cruz, Vandana Desai, Eric Bellm, Britt Lundgren, Frossie Economou, Brian D. Nord, Chad Schafer, Gautham Narayan, Joseph Harrington, Erik Tollerud, Brigitta Sipőcz, Timothy Pickering, Molly S. Peeples, Bruce Berriman, Peter Teuben, David Rodriguez, Andre Gradvohl, Lior Shamir, Alice Allen, Joel R. Brownstein, Adam Ginsburg, Manodeep Sinha, Cameron Hummels , et al. (20 additional authors not shown)

Abstract: Software is a critical part of modern research, and yet there are insufficient mechanisms in the scholarly ecosystem to acknowledge, cite, and measure the impact of research software. The majority of academic fields rely on a one-dimensional credit model whereby academic articles (and their associated citations) are the dominant factor in the success of a researcher's career. In the petabyte era o… ▽ More Software is a critical part of modern research, and yet there are insufficient mechanisms in the scholarly ecosystem to acknowledge, cite, and measure the impact of research software. The majority of academic fields rely on a one-dimensional credit model whereby academic articles (and their associated citations) are the dominant factor in the success of a researcher's career. In the petabyte era of astronomical science, citing software and measuring its impact enables academia to retain and reward researchers that make significant software contributions. These highly skilled researchers must be retained to maximize the scientific return from petabyte-scale datasets. Evolving beyond the one-dimensional credit model requires overcoming several key challenges, including the current scholarly ecosystem and scientific culture issues. This white paper will present these challenges and suggest practical solutions for elevating the role of software as a product of the research enterprise. △ Less

Submitted 14 July, 2019; originally announced July 2019.

Comments: arXiv admin note: substantial text overlap with arXiv:1905.05116

The 21cm bispectrum during reionization: a tracer of the ionization topology

Authors: Anne Hutter, Catherine A. Watkinson, Jacob Seiler, Pratika Dayal, Manodeep Sinha, Darren J. Croton

Abstract: We compute the bispectra of the 21cm signal during the Epoch of Reionization for three different reionization scenarios that are based on a dark matter N-body simulation combined with a self-consistent, semi-numerical model of galaxy evolution and reionization. Our reionization scenarios differ in their trends of ionizing escape fractions ($f_\mathrm{esc}$) with the underlying galaxy properties an… ▽ More We compute the bispectra of the 21cm signal during the Epoch of Reionization for three different reionization scenarios that are based on a dark matter N-body simulation combined with a self-consistent, semi-numerical model of galaxy evolution and reionization. Our reionization scenarios differ in their trends of ionizing escape fractions ($f_\mathrm{esc}$) with the underlying galaxy properties and cover the physically plausible range, i.e. $f_\mathrm{esc}$ effectively decreasing, being constant, or increasing with halo mass. We find the 21cm bispectrum to be sensitive to the resulting ionization topologies that significantly differ in their size distribution of ionized and neutral regions throughout reionization. From squeezed to stretched triangles, the 21cm bispectra features a change of sign from negative to positive values, with ionized and neutral regions representing below-average and above-average concentrations contributing negatively and positively, respectively. The position of the change of sign provides a tracer of the size distribution of the ionized and neutral regions, and allows us to identify three major regimes that the 21cm bispectrum undergoes during reionization. In particular the regime during the early stages of reionization, where the 21cm bispectrum tracks the peak of the size distribution of the ionized regions, provides exciting prospects for pinning down reionization with the forthcoming Square Kilometre Array. △ Less

Submitted 18 November, 2019; v1 submitted 9 July, 2019; originally announced July 2019.

Comments: 16 pages, 7 figures. Accepted for publication in MNRAS

The Escape Fraction of Ionizing Photons During the Epoch of Reionization: observability with the Square Kilometre Array

Authors: Jacob Seiler, Anne Hutter, Manodeep Sinha, Darren Croton

Abstract: One of the most important parameters in characterizing the Epoch of Reionization, the escape fraction of ionizing photons, $f_\mathrm{esc}$, remains unconstrained both observationally and theoretically. With recent work highlighting the impact of galaxy-scale feedback on the instantaneous value of $f_\mathrm{esc}$, it is important to develop a model in which reionization is self-consistently coupl… ▽ More One of the most important parameters in characterizing the Epoch of Reionization, the escape fraction of ionizing photons, $f_\mathrm{esc}$, remains unconstrained both observationally and theoretically. With recent work highlighting the impact of galaxy-scale feedback on the instantaneous value of $f_\mathrm{esc}$, it is important to develop a model in which reionization is self-consistently coupled to galaxy evolution. In this work, we present such a model and explore how physically motivated functional forms of $f_\mathrm{esc}$ affect the evolution of ionized hydrogen within the intergalactic medium. Using the $21$cm power spectrum evolution, we investigate the likelihood of observationally distinguishing between a constant $f_\mathrm{esc}$ and other models that depend upon different forms of galaxy feedback. We find that changing the underlying connection between $f_\mathrm{esc}$ and galaxy feedback drastically alters the large-scale $21$cm power. The upcoming Square Kilometre Array Low Frequency instrument possesses the sensitivity to differentiate between our models at a fixed optical depth, requiring only $200$ hours of integration time focused on redshifts $z = 7.5-8.5$. Generalizing these results to account for a varying optical depth will require multiple $800$ hour observations spanning redshifts $z = 7-10$. This presents an exciting opportunity to observationally constrain one of the most elusive parameters during the Epoch of Reionization. △ Less

Submitted 15 July, 2019; v1 submitted 5 February, 2019; originally announced February 2019.

Comments: 14 pages, 8 figures, 2 tables

Journal ref: MNRAS 487 (2019) 5739-5752

Model dispersion with PRISM; an alternative to MCMC for rapid analysis of models

Authors: Ellert van der Velden, Alan R. Duffy, Darren Croton, Simon J. Mutch, Manodeep Sinha

Abstract: We have built PRISM, a "Probabilistic Regression Instrument for Simulating Models". PRISM uses the Bayes linear approach and history matching to construct an approximation ('emulator') of any given model, by combining limited model evaluations with advanced regression techniques, covariances and probability calculations. It is designed to easily facilitate and enhance existing Markov chain Monte C… ▽ More We have built PRISM, a "Probabilistic Regression Instrument for Simulating Models". PRISM uses the Bayes linear approach and history matching to construct an approximation ('emulator') of any given model, by combining limited model evaluations with advanced regression techniques, covariances and probability calculations. It is designed to easily facilitate and enhance existing Markov chain Monte Carlo (MCMC) methods by restricting plausible regions and exploring parameter space efficiently. However, PRISM can additionally be used as a standalone alternative to MCMC for model analysis, providing insight into the behavior of complex scientific models. With PRISM, the time spent on evaluating a model is minimized, providing developers with an advanced model analysis for a fraction of the time required by more traditional methods. This paper provides an overview of the different techniques and algorithms that are used within PRISM. We demonstrate the advantage of using the Bayes linear approach over a full Bayesian analysis when analyzing complex models. Our results show how much information can be captured by PRISM and how one can combine it with MCMC methods to significantly speed up calibration processes (>15 times faster). PRISM is an open-source Python package that is available under the BSD 3-Clause License (BSD-3) at https://github.com/1313e/PRISM and hosted at https://prism-tool.readthedocs.io. PRISM has also been reviewed by "The Journal of Open Source Software" (https://doi.org/10.21105/joss.01229). △ Less

Submitted 11 June, 2019; v1 submitted 24 January, 2019; originally announced January 2019.

Comments: 28 pages, 13 figures, 1 table. Updated to reflect some changes made to published version

Journal ref: ApJS 242 22 (2019)

The Secondary Spin Bias of Dark Matter Haloes

Authors: James W. Johnson, Ariyeh H. Maller, Andreas A. Berlind, Manodeep Sinha, J. Kelly Holley-Bockelmann

Abstract: We investigate the role of angular momentum in the clustering of dark matter haloes. We make use of data from two high-resolution N-body simulations spanning over four orders of magnitude in halo mass, from $10^{9.8}$ to $10^{14}\ h^{-1}\ \text{M}_\odot$. We explore the hypothesis that mass accretion in filamentary environments alters the angular momentum of a halo, thereby driving a correlation b… ▽ More We investigate the role of angular momentum in the clustering of dark matter haloes. We make use of data from two high-resolution N-body simulations spanning over four orders of magnitude in halo mass, from $10^{9.8}$ to $10^{14}\ h^{-1}\ \text{M}_\odot$. We explore the hypothesis that mass accretion in filamentary environments alters the angular momentum of a halo, thereby driving a correlation between the spin parameter $λ$ and the strength of clustering. However, we do not find evidence that the distribution of matter on large scales is related to the spin of haloes. We find that a halo's spin is correlated with its age, concentration, sphericity, and mass accretion rate. Removing these correlations strongly affects the strength of secondary spin bias at low halo masses. We also find that high spin haloes are slightly more likely to be found near another halo of comparable mass. These haloes that are found near a comparable mass neighbour - a \textit{twin} - are strongly spatially biased. We demonstrate that this \textit{twin bias}, along with the relationship between spin and mass accretion rates, statistically accounts for halo spin secondary bias. △ Less

Submitted 5 December, 2018; originally announced December 2018.

Comments: 11 pages, 6 figures; submitted to MNRAS, comments welcome

The Three Hundred project: a large catalogue of theoretically modelled galaxy clusters for cosmological and astrophysical applications

Authors: Weiguang Cui, Alexander Knebe, Gustavo Yepes, Frazer Pearce, Chris Power, Romeel Dave, Alexander Arth, Stefano Borgani, Klaus Dolag, Pascal Elahi, Robert Mostoghiu, Giuseppe Murante, Elena Rasia, Doris Stoppacher, Jesus Vega-Ferrero, Yang Wang, Xiaohu Yang, Andrew Benson, Sofía A. Cora, Darren J. Croton, Manodeep Sinha, Adam R. H. Stevens, Cristian A. Vega-Martínez, Jake Arthur, Anna S. Baldi , et al. (12 additional authors not shown)

Abstract: We introduce the THE THREE HUNDRED project, an endeavour to model 324 large galaxy clusters with full-physics hydrodynamical re-simulations. Here we present the data set and study the differences to observations for fundamental galaxy cluster properties and scaling relations. We find that the modelled galaxy clusters are generally in reasonable agreement with observations with respect to baryonic… ▽ More We introduce the THE THREE HUNDRED project, an endeavour to model 324 large galaxy clusters with full-physics hydrodynamical re-simulations. Here we present the data set and study the differences to observations for fundamental galaxy cluster properties and scaling relations. We find that the modelled galaxy clusters are generally in reasonable agreement with observations with respect to baryonic fractions and gas scaling relations at redshift z = 0. However, there are still some (model-dependent) differences, such as central galaxies being too massive, and galaxy colours (g - r) being bluer (about 0.2 dex lower at the peak position) than in observations. The agreement in gas scaling relations down to 10^{13} h^{-1} M_{\odot} between the simulations indicates that particulars of the sub-grid modelling of the baryonic physics only has a weak influence on these relations. We also include - where appropriate - a comparison to three semi-analytical galaxy formation models as applied to the same underlying dark-matter-only simulation. All simulations and derived data products are publicly available. △ Less

Submitted 12 September, 2018; originally announced September 2018.

Comments: 20 pages, 8 figures, 7 tables. MNRAS published version

Journal ref: Cui, W., Knebe, A., Yepes, G., et al.\ 2018, \mnras, 480, 2898

The Indirect Influence of Quasars on Reionization

Authors: Jacob Seiler, Anne Hutter, Manodeep Sinha, Darren Croton

Abstract: The exact role of quasars during the Epoch of Reionization remains uncertain. With consensus leaning towards quasars producing a negligible amount of ionizing photons, we pose an alternate question: Can quasars indirectly contribute to reionization by allowing ionizing photons from stars to escape more easily? Using the Semi-Analytic Galaxy Evolution model to evolve a galaxy population through cos… ▽ More The exact role of quasars during the Epoch of Reionization remains uncertain. With consensus leaning towards quasars producing a negligible amount of ionizing photons, we pose an alternate question: Can quasars indirectly contribute to reionization by allowing ionizing photons from stars to escape more easily? Using the Semi-Analytic Galaxy Evolution model to evolve a galaxy population through cosmic time, we construct an idealized scenario in which the escape fraction of stellar ionizing photons ($f_\mathrm{esc}$) is boosted following quasar wind events, potentially for several dynamical times. We find that under this scenario, the mean value of $f_\mathrm{esc}$ as a function of galaxy stellar mass peaks for intermediate mass galaxies. This mass dependence will have consequences for the 21cm power spectrum, enhancing power at small scales and suppressing it at large scales. This hints that whilst quasars may not directly contribute to the ionizing photon budget, they could influence reionization indirectly by altering the topology of ionized regions. △ Less

Submitted 29 June, 2018; originally announced June 2018.

Comments: 5 pages, 3 figures

Connecting and dissecting galaxies' angular momenta and neutral gas in a hierarchical universe: cue DARK SAGE

Authors: Adam R. H. Stevens, Claudia del P. Lagos, Danail Obreschkow, Manodeep Sinha

Abstract: We explore the connection between the atomic gas fraction, f_atm, and `global disc stability' parameter, q, of galaxies within a fully cosmological context by examining galaxies in the Dark Sage semi-analytic model. The q parameter is determined by the ratio of disc specific angular momentum to mass. Dark Sage is well suited to our study, as it includes the numerical evolution of one-dimensional d… ▽ More We explore the connection between the atomic gas fraction, f_atm, and `global disc stability' parameter, q, of galaxies within a fully cosmological context by examining galaxies in the Dark Sage semi-analytic model. The q parameter is determined by the ratio of disc specific angular momentum to mass. Dark Sage is well suited to our study, as it includes the numerical evolution of one-dimensional disc structure, making both j_disc and q predicted quantities. We show that Dark Sage produces a clear correlation between gas fraction and j_disc at fixed disc mass, in line with recent results from observations and hydrodynamic simulations. This translates to a tight q--f_atm sequence for star-forming central galaxies, which closely tracks the analytic prediction of Obreschkow et al. The scatter in this sequence is driven by the probability distribution function of mass as a function of j (PDF of j) within discs, specifically where it peaks. We find that halo mass is primarily responsible for the peak location of the PDF of j, at least for low values of q. Two main mechanisms of equal significance are then identified for disconnecting f_atm from q. Mergers in the model can trigger quasar winds, with the potential to blow out most of the gas disc, while leaving the stellar disc relatively unharmed. Ram-pressure stripping of satellite galaxies has a similar effect, where f_atm can drop drastically with only a minimal effect to q. We highlight challenges associated with following these predictions up with observations. △ Less

Submitted 17 August, 2018; v1 submitted 19 June, 2018; originally announced June 2018.

Comments: 13 pages, 8 figures (excluding references and appendices). Submitted to MNRAS. Revisions after referee's report

Likelihood Non-Gaussianity in Large-Scale Structure Analyses

Authors: ChangHoon Hahn, Florian Beutler, Manodeep Sinha, Andreas Berlind, Shirley Ho, David W. Hogg

Abstract: Standard present day large-scale structure (LSS) analyses make a major assumption in their Bayesian parameter inference --- that the likelihood has a Gaussian form. For summary statistics currently used in LSS, this assumption, even if the underlying density field is Gaussian, cannot be correct in detail. We investigate the impact of this assumption on two recent LSS analyses: the Beutler et al. (… ▽ More Standard present day large-scale structure (LSS) analyses make a major assumption in their Bayesian parameter inference --- that the likelihood has a Gaussian form. For summary statistics currently used in LSS, this assumption, even if the underlying density field is Gaussian, cannot be correct in detail. We investigate the impact of this assumption on two recent LSS analyses: the Beutler et al. (2017) power spectrum multipole ($P_\ell$) analysis and the Sinha et al. (2017) group multiplicity function ($ζ$) analysis. Using non-parametric divergence estimators on mock catalogs originally constructed for covariance matrix estimation, we identify significant non-Gaussianity in both the $P_\ell$ and $ζ$ likelihoods. We then use Gaussian mixture density estimation and Independent Component Analysis on the same mocks to construct likelihood estimates that approximate the true likelihood better than the Gaussian $pseudo$-likelihood. Using these likelihood estimates, we accurately estimate the true posterior probability distribution of the Beutler et al. (2017) and Sinha et al. (2017) parameters. Likelihood non-Gaussianity shifts the $fσ_8$ constraint by $-0.44σ$, but otherwise, does not significantly impact the overall parameter constraints of Beutler et al. (2017). For the $ζ$ analysis, using the pseudo-likelihood significantly underestimates the uncertainties and biases the constraints of Sinha et al. (2017) halo occupation parameters. For $\log M_1$ and $α$, the posteriors are shifted by $+0.43σ$ and $-0.51σ$ and broadened by $42\%$ and $66\%$, respectively. The divergence and likelihood estimation methods we present provide a straightforward framework for quantifying the impact of likelihood non-Gaussianity and deriving more accurate parameter constraints. △ Less

Submitted 16 March, 2018; originally announced March 2018.

Comments: 33 pages, 7 figures

Small- and Large-Scale Galactic Conformity in SDSS DR7

Authors: Victor F. Calderon, Andreas A. Berlind, Manodeep Sinha

Abstract: Galactic conformity is the phenomenon whereby galaxy properties exhibit excess correlations across distance than that expected if these properties only depended on halo mass. We perform a comprehensive study of conformity at low redshift using a galaxy group catalogue from the SDSS DR7 spectroscopic sample. We study correlations both between central galaxies and their satellites (1-halo), and betw… ▽ More Galactic conformity is the phenomenon whereby galaxy properties exhibit excess correlations across distance than that expected if these properties only depended on halo mass. We perform a comprehensive study of conformity at low redshift using a galaxy group catalogue from the SDSS DR7 spectroscopic sample. We study correlations both between central galaxies and their satellites (1-halo), and between central galaxies in separate haloes (2-halo). We use the quenched fractions and the marked correlation function (MCF), to probe for conformity in three galaxy properties, $(g-r)$ colour, specific star formation rate (sSFR), and morphology. We assess the statistical significance of conformity signals with a suite of mock galaxy catalogues that have no built-in conformity, but contain the same group-finding and mass assignment errors as the real data. In the case of 1-halo conformity, quenched fractions show strong signals at all group masses. However, these signals are equally strong in mock catalogues, indicating that the conformity signal is spurious and likely entirely caused by group-finding systematics, calling into question previous claims of 1-halo conformity detection. The MCF reveals a significant detection of radial segregation within massive groups, but no evidence of conformity. In the case of 2-halo conformity, quenched fractions show no significant evidence of conformity in colour or sSFR once compared with mock catalogues, but a clear signal using morphology. In contrast, the MCF reveals a small, yet highly significant signal for all three properties in low mass groups and scales of $0.8-4\ h^{-1}\textrm{Mpc}$, possibly representing the first robust detection of 2-halo conformity. △ Less

Submitted 15 August, 2018; v1 submitted 7 December, 2017; originally announced December 2017.

Comments: 16 pages, 5 figures, accepted by MNRAS

MultiDark-Galaxies: data release and first results

Authors: Alexander Knebe, Doris Stoppacher, Francisco Prada, Christoph Behrens, Andrew Benson, Sofia A. Cora, Darren J. Croton, Nelson D. Padilla, Andrés N. Ruiz, Manodeep Sinha, Adam R. H. Stevens, Cristian A. Vega-Martínez, Peter Behroozi, Violeta Gonzalez-Perez, Stefan Gottlöber, Anatoly A. Klypin, Gustavo Yepes, Harry Enke, Noam I. Libeskind, Kristin Riebe, Matthias Steinmetz

Abstract: We present the public release of the MultiDark-Galaxies: three distinct galaxy catalogues derived from one of the Planck cosmology MultiDark simulations (i.e. MDPL2, with a volume of (1 Gpc/$h$)$^{3}$ and mass resolution of $1.5 \times 10^{9} M_{\odot}/h$) by applying the semi-analytic models GALACTICUS, SAG, and SAGE to it. We compare the three models and their conformity with observational data… ▽ More We present the public release of the MultiDark-Galaxies: three distinct galaxy catalogues derived from one of the Planck cosmology MultiDark simulations (i.e. MDPL2, with a volume of (1 Gpc/$h$)$^{3}$ and mass resolution of $1.5 \times 10^{9} M_{\odot}/h$) by applying the semi-analytic models GALACTICUS, SAG, and SAGE to it. We compare the three models and their conformity with observational data for a selection of fundamental properties of galaxies like stellar mass function, star formation rate, cold gas fractions, and metallicities - noting that they sometimes perform differently reflecting model designs and calibrations. We have further selected galaxy subsamples of the catalogues by number densities in stellar mass, cold gas mass, and star formation rate in order to study the clustering statistics of galaxies. We show that despite different treatment of orphan galaxies, i.e. galaxies that lost their dark-matter host halo due to the finite mass resolution of the N-body simulation or tidal stripping, the clustering signal is comparable, and reproduces the observations in all three models - in particular when selecting samples based upon stellar mass. Our catalogues provide a powerful tool to study galaxy formation within a volume comparable to those probed by on-going and future photometric and redshift surveys. All model data consisting of a range of galaxy properties - including broad-band SDSS magnitudes - are publicly available. △ Less

Submitted 23 October, 2017; originally announced October 2017.

Comments: 29 pages, 16 figures, 8 tables, accepted for publication in MNRAS. All data incl. the complete galaxy catalogues for all models are publicly available from the CosmoSim database (http://www.cosmosim.org); a selected set of galaxy properties is available via the Skies & Universes website (http://www.skiesanduniverses.org)

Spatial Clustering of Dark Matter Halos: Secondary Bias, Neighbor Bias, and the Influence of Massive Neighbors on Halo Properties

Authors: Andrés N. Salcedo, Ariyeh H. Maller, Andreas A. Berlind, Manodeep Sinha, Cameron K. McBride, Peter S. Behroozi, Risa H. Wechsler, David H. Weinberg

Abstract: We explore the phenomenon commonly known as halo assembly bias, whereby dark matter halos of the same mass are found to be more or less clustered when a second halo property is considered, for halos in the mass range $3.7 \times 10^{11} \; h^{-1} \mathrm{M_{\odot}} - 5.0 \times 10^{13} \; h^{-1} \mathrm{M_{\odot}}$. Using the Large Suite of Dark Matter Simulations (LasDamas) we consider nine commo… ▽ More We explore the phenomenon commonly known as halo assembly bias, whereby dark matter halos of the same mass are found to be more or less clustered when a second halo property is considered, for halos in the mass range $3.7 \times 10^{11} \; h^{-1} \mathrm{M_{\odot}} - 5.0 \times 10^{13} \; h^{-1} \mathrm{M_{\odot}}$. Using the Large Suite of Dark Matter Simulations (LasDamas) we consider nine commonly used halo properties and find that a clustering bias exists if halos are binned by mass or by any other halo property. This secondary bias implies that no single halo property encompasses all the spatial clustering information of the halo population. The mean values of some halo properties depend on their halo's distance to a more massive neighbor. Halo samples selected by having high values of one of these properties therefore inherit a neighbor bias such that they are much more likely to be close to a much more massive neighbor. This neighbor bias largely accounts for the secondary bias seen in halos binned by mass and split by concentration or age. However, halos binned by other mass-like properties still show a secondary bias even when the neighbor bias is removed. The secondary bias of halos selected by their spin behaves differently than that for other halo properties, suggesting that the origin of the spin bias is different than of other secondary biases. △ Less

Submitted 9 January, 2018; v1 submitted 28 August, 2017; originally announced August 2017.

Comments: 14 pages, LaTeX; minor revisions, and added references; results unchanged

Towards Accurate Modelling of Galaxy Clustering on Small Scales: Testing the Standard $Λ\mathrm{CDM}$ + Halo Model

Authors: Manodeep Sinha, Andreas A. Berlind, Cameron K. McBride, Roman Scoccimarro, Jennifer A. Piscionere, Benjamin D. Wibking

Abstract: Interpreting the small-scale clustering of galaxies with halo models can elucidate the connection between galaxies and dark matter halos. Unfortunately, the modelling is typically not sufficiently accurate for ruling out models statistically. It is thus difficult to use the information encoded in small scales to test cosmological models or probe subtle features of the galaxy-halo connection. In th… ▽ More Interpreting the small-scale clustering of galaxies with halo models can elucidate the connection between galaxies and dark matter halos. Unfortunately, the modelling is typically not sufficiently accurate for ruling out models statistically. It is thus difficult to use the information encoded in small scales to test cosmological models or probe subtle features of the galaxy-halo connection. In this paper, we attempt to push halo modelling into the "accurate" regime with a fully numerical mock-based methodology and careful treatment of statistical and systematic errors. With our forward-modelling approach, we can incorporate clustering statistics beyond the traditional two-point statistics. We use this modelling methodology to test the standard $Λ\mathrm{CDM}$ + halo model against the clustering of SDSS DR7 galaxies. Specifically, we use the projected correlation function, group multiplicity function and galaxy number density as constraints. We find that while the model fits each statistic separately, it struggles to fit them simultaneously. Adding group statistics leads to a more stringent test of the model and significantly tighter constraints on model parameters. We explore the impact of varying the adopted halo definition and cosmological model and find that changing the cosmology makes a significant difference. The most successful model we tried (Planck cosmology with Mvir halos) matches the clustering of low luminosity galaxies, but exhibits a 2.3$σ$ tension with the clustering of luminous galaxies, thus providing evidence that the "standard" halo model needs to be extended. This work opens the door to adding interesting freedom to the halo model and including additional clustering statistics as constraints. △ Less

Submitted 6 August, 2018; v1 submitted 16 August, 2017; originally announced August 2017.

Comments: Replaced to match the published version

Journal ref: http://adsabs.harvard.edu/abs/2018MNRAS.478.1042S

From microphysics to dynamics of magnetars

Authors: Armen Sedrakian, Xu-Guang Huang, Monika Sinha, John W. Clark

Abstract: MeV-scale magnetic fields in the interiors of magnetars suppress the pairing of neutrons and protons in the $S$-wave state. In the case of a neutron condensate the suppression is the consequence of the Pauli-paramagnetism of the neutron gas, i.e., the alignment of the neutron spins along the magnetic field. The proton $S$-wave pairing is suppressed because of the Landau diamagnetic currents of pro… ▽ More MeV-scale magnetic fields in the interiors of magnetars suppress the pairing of neutrons and protons in the $S$-wave state. In the case of a neutron condensate the suppression is the consequence of the Pauli-paramagnetism of the neutron gas, i.e., the alignment of the neutron spins along the magnetic field. The proton $S$-wave pairing is suppressed because of the Landau diamagnetic currents of protons induced by the field. The Ginzburg-Landau and BCS theories of the critical magnetic fields for unpairing are reviewed. The macrophysical implications of the suppression (unpairing) of the condensates are discussed for the rotational crust-core coupling in magnetars and the neutrino-dominated cooling era of their thermal evolution. △ Less

Submitted 3 January, 2017; originally announced January 2017.

Comments: 10 pages, 4 figures, Proceedings of "Compact Stars in the QCD phase diagram V", 23-27 May 2016 GSSI and LNGS, L'Aquila, Italy

Journal ref: Journal of Physics: Conf. Series 861 (2017) 012025

The Astropy Problem

Authors: Demitri Muna, Michael Alexander, Alice Allen, Richard Ashley, Daniel Asmus, Ruyman Azzollini, Michele Bannister, Rachael Beaton, Andrew Benson, G. Bruce Berriman, Maciej Bilicki, Peter Boyce, Joanna Bridge, Jan Cami, Eryn Cangi, Xian Chen, Nicholas Christiny, Christopher Clark, Michelle Collins, Johan Comparat, Neil Cook, Darren Croton, Isak Delberth Davids, Éric Depagne, John Donor , et al. (129 additional authors not shown)

Abstract: The Astropy Project (http://astropy.org) is, in its own words, "a community effort to develop a single core package for Astronomy in Python and foster interoperability between Python astronomy packages." For five years this project has been managed, written, and operated as a grassroots, self-organized, almost entirely volunteer effort while the software is used by the majority of the astronomical… ▽ More The Astropy Project (http://astropy.org) is, in its own words, "a community effort to develop a single core package for Astronomy in Python and foster interoperability between Python astronomy packages." For five years this project has been managed, written, and operated as a grassroots, self-organized, almost entirely volunteer effort while the software is used by the majority of the astronomical community. Despite this, the project has always been and remains to this day effectively unfunded. Further, contributors receive little or no formal recognition for creating and supporting what is now critical software. This paper explores the problem in detail, outlines possible solutions to correct this, and presents a few suggestions on how to address the sustainability of general purpose astronomical software. △ Less

Submitted 10 October, 2016; originally announced October 2016.