Researcher profile

Raul E. Angulo

Raul E. Angulo contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
11works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

11 published item(s)

preprint2025arXiv

J-PLUS: Turning Off the Bright Stars

Photometric surveys require precise point spread function (PSF) characterization, as it varies across filters and is crucial for accurate photometry and low surface brightness (LSB) studies. However, the small PSF size provided by default pipelines suits only barely resolved objects, making it difficult to analyze regions near bright stars (rendering those regions unusable). These components are then combined to generate a final PSF for each exposure and filter, spanning 15 mag arcsec-2 in surface brightness and 4 arcmin in radius in the broad bands. In narrow-band filters, the J-PLUS PSF exhibits two rings, whereas in broad-band filters, only one ring is observed. Additionally, the position of the ring shifts with filter wavelength: as the filters become redder, the ring radius increases. We find that there is no significant variation in the extended PSF observed as a function of time (within 2.5h) or position in the field of view. The radial profile of NGC 4212 (which is close to a star) is also studied before/after PSF-subtraction. We developed a novel method to determine the central coordinates of saturated stars, and classify stars without using Gaia magnitudes. Additionally, mirror reflections are automatically detected and masked. Furthermore, in combining different stars and various components of the PSF, we avoided the use of a fixed radius by introducing a new method that does not depend on radial measurements. Accurate characterization of the extended PSF and its subtraction improves sky subtraction, increases the effective area of the survey by about 10%, and enables the study of extended large LSB features in wide area surveys like J-PLUS. Our pipeline is published as free software (GNU GPLv3) an can be customized to other surveys such as J-PAS, where its impact will be even greater due to its depth. This paper is fully reproducible and produced from Commit 4860c70.

preprint2022arXiv

Accelerating Large-Scale-Structure data analyses by emulating Boltzmann solvers and Lagrangian Perturbation Theory

The linear matter power spectrum is an essential ingredient in all theoretical models for interpreting large-scale-structure observables. Although Boltzmann codes such as CLASS or CAMB are very efficient at computing the linear spectrum, the analysis of data usually requires $10^4$-$10^6$ evaluations, which means this task can be the most computationally expensive aspect of data analysis. Here, we address this problem by building a neural network emulator that provides the linear theory (total and cold) matter power spectrum in about one millisecond with 0.2% (0.5%) accuracy over redshifts $z \le 3$ ($z \le 9$), and scales $10^{-4} \le k \, [h {\rm Mpc^{-1}}] < 50$. We train this emulator with more than 200,000 measurements, spanning a broad cosmological parameter space that includes massive neutrinos and dynamical dark energy. We show that the parameter range and accuracy of our emulator is enough to get unbiased cosmological constraints in the analysis of a Euclid-like weak lensing survey. Complementing this emulator, we train 15 other emulators for the cross-spectra of various linear fields in Eulerian space, as predicted by 2nd-order Lagrangian Perturbation theory, which can be used to accelerate perturbative bias descriptions of galaxy clustering. Our emulators are specially designed to be used in combination with emulators for the nonlinear matter power spectrum and for baryonic effects, all of which are publicly available at http://www.dipc.org/bacco.

preprint2022arXiv

J-PLUS: A catalogue of globular cluster candidates around the M81/M82/NGC3077 triplet of galaxies

Globular clusters (GCs) are proxies of the formation assemblies of their host galaxies. However, few studies exist targeting GC systems of spiral galaxies up to several effective radii. Through 12-band Javalambre Photometric Local Universe Survey (J-PLUS) imaging, we study the point sources around the M81/M82/NGC3077 triplet in search of new GC candidates. We develop a tailored classification scheme to search for GC candidates based on their similarity to known GCs via a principal components analysis (PCA) projection. Our method accounts for missing data and photometric errors. We report 642 new GC candidates in a region of 3.5 deg$^2$ around the triplet, ranked according to their Gaia astrometric proper motions when available. We find tantalising evidence for an overdensity of GC candidate sources forming a bridge connecting M81 and M82. Finally, the spatial distribution of the GC candidates $(g-i)$ colours is consistent with halo/intra-cluster GCs, i.e. it gets bluer as they get further from the closest galaxy in the field. We further employ a regression-tree based model to estimate the metallicity distribution of the GC candidates based on their J-PLUS bands. The metallicity distribution of the sample candidates is broad and displays a bump towards the metal-rich end. Our list increases the population of GC candidates around the triplet by 3-fold, stresses the usefulness of multi-band surveys in finding these objects, and provides a testbed for further studies analysing their spatial distribution around nearby (spirals) galaxies.

preprint2022arXiv

Priors on Lagrangian bias parameters from galaxy formation modelling

We study the relations among the parameters of the hybrid Lagrangian bias expansion model, fitting biased auto and cross power spectra up to $k_{\rm max} = 0.7 \, h \, \mathrm{Mpc}^{-1}$. We consider $\sim 8000$ halo and galaxy samples, with different halo masses, redshifts, galaxy number densities, and varying the parameters of the galaxy formation model. Galaxy samples are obtained through state-of-the-art extended subhalo abundance matching techniques and include both stellar-mass and star-formation-rate selected galaxies. All of these synthetic galaxies samples are publicly available at https://bacco.dipc.org/galpk.html. We find that the hybrid Lagrangian bias model provides accurate fits to all of our halo and galaxy samples. The coevolution relations between galaxy bias parameters, although roughly compatible with those obtained for haloes, show systematic shifts and larger scatter. We explore possible sources of this difference in terms of dependence on halo occupation and assembly bias of each sample. The bias parameter relations displayed in this work can be used as a prior for future Bayesian analyses employing the hybrid Lagrangian bias expansion model.

preprint2021arXiv

J-PLUS: Stellar Parameters, C, N, Mg, Ca and [α/Fe] Abundances for Two Million Stars from DR1

Context. The Javalambre Photometric Local Universe Survey (J-PLUS) has obtained precise photometry in twelve specially designed filters for large numbers of Galactic stars. Deriving their precise stellar atmospheric parameters and individual elemental abundances is crucial for studies of Galactic structure, and the assembly history and chemical evolution of our Galaxy. Aims. Our goal is to estimate not only stellar parameters (effective temperature, Teff, surface gravity, log g, and metallicity, [Fe/H]), but also [α/Fe] and four elemental abundances ([C/Fe], [N/Fe], [Mg/Fe], and [Ca/Fe]) using data from J-PLUS DR1. Methods. By combining recalibrated photometric data from J-PLUS DR1, Gaia DR2, and spectroscopic labels from LAMOST, we design and train a set of cost-sensitive neural networks, the CSNet, to learn the non-linear mapping from stellar colors to their labels. Results. We have achieved precisions of δTeff {\sim}55K, δlogg{\sim}0.15dex, and δ[Fe/H]{\sim}0.07dex, respectively, over a wide range of temperature, surface gravity, and metallicity. The uncertainties of the abundance estimates for [α/Fe] and the four individual elements are in the range 0.04-0.08 dex. We compare our parameter and abundance estimates with those from other spectroscopic catalogs such as APOGEE and GALAH, and find an overall good agreement. Conclusions. Our results demonstrate the potential of well-designed, high-quality photometric data for determinations of stellar parameters as well as individual elemental abundances. Applying the method to J-PLUS DR1, we have obtained the aforementioned parameters for about two million stars, providing an outstanding data set for chemo-dynamic analyses of the Milky Way. The catalog of the estimated parameters is publicly accessible.

preprint2021arXiv

Large-scale dark matter simulations

We review the field of collisionless numerical simulations for the large-scale structure of the Universe. We start by providing the main set of equations solved by these simulations and their connection with General Relativity. We then recap the relevant numerical approaches: discretization of the phase-space distribution (focusing on N-body but including alternatives, e.g., Lagrangian submanifold and Schrödinger-Poisson) and the respective techniques for their time evolution and force calculation (Direct summation, mesh techniques, and hierarchical tree methods). We pay attention to the creation of initial conditions and the connection with Lagrangian Perturbation Theory. We then discuss the possible alternatives in terms of the micro-physical properties of dark matter (e.g., neutralinos, warm dark matter, QCD axions, Bose-Einstein condensates, and primordial black holes), and extensions to account for multiple fluids (baryons and neutrinos), primordial non-Gaussianity and modified gravity. We continue by discussing challenges involved in achieving highly accurate predictions. A key aspect of cosmological simulations is the connection to cosmological observables, we discuss various techniques in this regard: structure finding, galaxy formation and baryonic modelling, the creation of emulators and light-cones, and the role of machine learning. We finalise with a recount of state-of-the-art large-scale simulations and conclude with an outlook for the next decade.

preprint2021arXiv

Modelling galaxy clustering in redshift space with a Lagrangian bias formalism and $N$-body simulations

Improving the theoretical description of galaxy clustering on small scales is an important challenge in cosmology, as it can considerably increase the scientific return of forthcoming galaxy surveys -- e.g. tightening the bounds on neutrino masses and deviations from general relativity. In this paper, we propose and test a new model for the clustering of galaxies that is able to accurately describe redshift-space distortions even down to small scales. This model corresponds to a second-order perturbative Lagrangian bias expansion which is advected to Eulerian space employing a displacement field extracted from $N$-body simulations. Eulerian coordinates are then transformed into redshift space by directly employing simulated velocity fields augmented with nuisance parameters capturing various possible satellite fractions and intra-halo small-scale velocities. We quantify the accuracy of our approach against samples of physically-motivated mock galaxies selected according to either Stellar Mass (SM) or Star Formation Rate (SFR) at multiple abundances and at $z=0$ and $1$. We find our model describes the monopole, quadrupole, and hexadecapole of the galaxy-power spectra down to scales of $k\approx 0.6 [h/$Mpc] within the accuracy of our simulations. This approach could pave the way to significantly increase the amount of cosmological information to be extracted from future galaxy surveys.

preprint2020arXiv

Modelling the large scale structure of the Universe as a function of cosmology and baryonic physics

We present and test a framework that models the three-dimensional distribution of mass in the Universe as a function of cosmological and astrophysical parameters. Our approach combines two different techniques: a rescaling algorithm that modifies the cosmology of gravity-only N-body simulations, and a baryonification algorithm which mimics the effects of astrophysical processes induced by baryons, such as star formation and AGN feedback. We show how this approach can accurately reproduce the effects of baryons on the matter power spectrum of various state-of-the-art hydro-dynamical simulations (EAGLE, Illustris, Illustris-TNG, Horizon-AGN, and OWLS,Cosmo-OWLS and BAHAMAS), to percent level from very large down to small, highly nonlinear scales, k= 5 h/Mpc, and from z=0 up to z=2. We highlight that, thanks to the heavy optimisation of the algorithms, we can obtain these predictions for arbitrary baryonic models and cosmology (including massive neutrinos and dynamical dark energy models) with an almost negligible CPU cost. Therefore, this approach is efficient enough for cosmological data analyses. With these tools in hand we explore the degeneracies between cosmological and astrophysical parameters in the nonlinear mass power spectrum. Our findings suggest that after marginalising over baryonic physics, cosmological constraints inferred from weak gravitational lensing should be moderately degraded.

preprint2020arXiv

Simulating the Complexity of the Dark Matter Sheet I: Numerical Algorithms

At early times dark matter has a thermal velocity dispersion of unknown amplitude which, for warm dark matter models, can influence the formation of nonlinear structure on observable scales. We propose a new scheme to simulate cosmologies with a small-scale suppression of perturbations that combines two previous methods in a way that avoids the numerical artefacts which have so far prevented either from producing fully reliable results. At low densities and throughout most of the cosmological volume, we represent the dark matter phase-sheet directly using high-accuracy interpolation, thereby avoiding the artificial fragmentation which afflicts particle-based methods in this regime. Such phase-sheet methods are, however, unable to follow the rapidly increasing complexity of the denser regions of dark matter haloes, so for these we switch to an N-body scheme which uses the geodesic deviation equation to track phase-sheet properties local to each particle. In addition, we present a novel high-resolution force calculation scheme based on an oct-tree of cubic force resolution elements which is well suited to approximate the force-field of our combined sheet+particle distribution. Our hybrid simulation scheme enables the first reliable simulations of the internal structure of low-mass haloes in a warm dark matter cosmology.

preprint2020arXiv

The manifestation of secondary halo bias on the galaxy population from IllustrisTNG300

We use the improved IllustrisTNG300 magneto-hydrodynamical cosmological simulation to revisit the effect that secondary halo bias has on the clustering of the central galaxy population. With a side length of 205 $h^{-1}$Mpc and significant improvements on the sub-grid model with respect to the previous Illustris boxes, IllustrisTNG300 allows us to explore the dependencies of galaxy clustering over a large cosmological volume and wide halo-mass range. We show, at high statistical significance, that the halo assembly bias signal (i.e., the secondary dependence of halo bias on halo formation redshift) manifests itself on the clustering of the central galaxy population when this is split by stellar mass, colour, specific star formation rate, and surface density. A significant detection is also obtained for galaxy size: at fixed halo mass, larger central galaxies are more tightly clustered than smaller central galaxies in haloes of mass M$_{\rm vir} \lesssim 10^{12.5}$ $h^{-1}$M$_{\odot}$. This effect, however, seems to be uncorrelated with halo formation time, unlike the rest of the secondary dependencies analysed. We also explore the transmission of the halo spin bias signal, i.e., the secondary dependence of halo bias on halo spin. Although galaxy spin retains little information about the total spin of the halo, the correlation is enough to produce a significant galaxy spin bias signal. We discuss possible ways to probe the spin bias effects with observations.

preprint2009arXiv

Simulations of Quintessential Cold Dark Matter: beyond the cosmological constant

We study the nonlinear growth of cosmic structure in different dark energy models, using large volume N-body simulations. We consider a range of quintessence models which feature both rapidly and slowly varying dark energy equations of state, and compare the growth of structure to that in a universe with a cosmological constant. The adoption of a quintessence model changes the expansion history of the universe, the form of the linear theory power spectrum and can alter key observables, such as the horizon scale and the distance to last scattering. We incorporate these effects into our simulations in stages to isolate the impact of each on the growth of structure. The difference in structure formation can be explained to first order by the difference in growth factor at a given epoch; this scaling also accounts for the nonlinear growth at the 15% level. We find that quintessence models that are different from $Λ$CDM both today and at high redshifts $(z \sim 1000)$ and which feature late $(z<2)$, rapid transitions in the equation of state, can have identical baryonic acoustic oscillation (BAO) peak positions to those in $Λ$CDM. We find that these models have higher abundances of dark matter haloes at $z>0$ compared to $Λ$CDM and so measurements of the mass function should allow us to distinguish these quintessence models from a cosmological constant. However, we find that a second class of quintessence models, whose equation of state makes an early $(z>2)$ rapid transition to $w=-1$, cannot be distinguished from $Λ$CDM using measurements of the mass function or the BAO, even if these models have non-negligible amounts of dark energy at early times.