Source author record

Octavio Valenzuela

Octavio Valenzuela appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

19works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

19 published item(s)

preprint2022arXiv

Boltzmann-Poisson-like approach to simulating the galactic halo response to satellite accretion Dependence on the halo density profile

Recent studies have reported the detection of the galactic stellar halo wake and dipole triggered by the Large Magellanic Cloud (LMC), mirroring the corresponding response from dark matter (DM). These studies open up the possibility of adding constraints on the global mass distribution of the Milky Way (MW), and even on the nature of DM itself, with current and upcoming stellar surveys reigniting the discussion on response modes in dynamical friction. However, the simulation of such features remains computationally challenging. We used a collisionless Boltzmann equation (CBE)+Poisson solver based on an existing method from the literature. We investigated the density and velocity response modes in simulations of Galactic-type DM halos accreting LMC-sized satellites, including the dependence on the halo density profile. We successfully captured both the local wake and the global over- and underdensity induced in the host halo. We also captured the velocity response. In line with previous studies, we find that the code can reproduce the core formation in the cuspy profile and the satellite core stalling. The angular power spectrum (APS) response is shown to be sensitive to each density profile. The cored Plummer density profile seems the most responsive, displaying a richness of modes. At the end of the simulation, the central halo acquires cylindrical rotation. The CBE description makes it tenable to capture the response modes with a better handling of noise in comparison to traditional N-body simulations. Hence, given a certain noise level, BPM has a lower computational cost than N-body simulations, making it feasible to explore large parameter sets. We anticipate that stellar spheroids in the MW or external galaxies could show central cylindrical rotation if they underwent a massive accretion event. The code can be adjusted to include a variety of DM physics.

preprint2022arXiv

Including relativistic and primordial Non-Gaussianity contributions in cosmological simulations by modifying the initial condition

We present a method to implement relativistic corrections to the evolution of dark matter structures in Newtonian simulations of a LCDM universe via the initial conditions. We take the nonlinear correspondence between the Lagrangian (Newtonian) evolution of dark matter inhomogeneities and the synchronous-comoving (relativistic) matter density description, and use it to promote the relativistic constraint as the initial condition for numerical simulations of structure formation. In this case, the incorporation of Primordial non-Gaussianity (PNG) contributions as initial conditions is straightforward. We implement the relativistic,fNL and gNLcontributions as initial conditions for the L-PICOLA code, and compute the power spectrum and bispectrum of the evolved matter field. We focus specifically on the case of largest values of non-Gaussianity allowed at 1-sigma by Planck observations(fNL=-4.2 and gNL=-7000). As a checkup, we show consistency with the one-loop perturbative prescription and with a fully relativistic simulation (GRAMSES) on the adequate scales. Our results confirm that both relativistic and PNG features are most prominent at very large scales and for squeezed triangulations. We discuss future prospects to probe these two contributions in the bispectrum of the matter density distribution

preprint2022arXiv

Phase spirals in cosmological simulations of Milky Way-size galaxies

We study the vertical perturbations in the galactic disc of the Milky Way-size high-resolution hydrodynamical cosmological simulation named GARROTXA. We detect phase spirals in the vertical projection $Z- V_{Z}$ of disc's stellar particles for the first time in this type of simulations. Qualitatively similar structures were detected in the recent Gaia data, and their origin is still under study. In our model the spiral-like structures in the phase space are present in a wide range of times and locations across the disc. By accounting for an evolving mix of stellar populations, we observe that, as seen in the data, the phase spirals are better observed in the range of younger-intermediate star particles. We measure the intensity of the spiral with a Fourier decomposition and find that these structures appear stronger near satellite pericenters. Current dynamical models of the phase spiral considering a single perturber required a mass at least of the order of 10$^{10}$ M$_\odot$, but all three of our satellites have masses of the order of $\sim$10$^8$ M$_\odot$. We suggest that there are other mechanisms at play which appear naturally in our model such as the physics of gas, collective effect of multiple perturbers, and a dynamically cold population that is continuously renovated by the star formation Complementing collisionless isolated N-body models with the use of fully-cosmological simulations with enough resolution can provide new insights into the nature/origin of the phase spiral.

preprint2022arXiv

Snowmass2021 Cosmic Frontier White Paper: Prospects for obtaining Dark Matter Constraints with DESI

Despite efforts over several decades, direct-detection experiments have not yet led to the discovery of the dark matter (DM) particle. This has led to increasing interest in alternatives to the Lambda CDM (LCDM) paradigm and alternative DM scenarios (including fuzzy DM, warm DM, self-interacting DM, etc.). In many of these scenarios, DM particles cannot be detected directly and constraints on their properties can ONLY be arrived at using astrophysical observations. The Dark Energy Spectroscopic Instrument (DESI) is currently one of the most powerful instruments for wide-field surveys. The synergy of DESI with ESA's Gaia satellite and future observing facilities will yield datasets of unprecedented size and coverage that will enable constraints on DM over a wide range of physical and mass scales and across redshifts. DESI will obtain spectra of the Lyman-alpha forest out to z~5 by detecting about 1 million QSO spectra that will put constraints on clustering of the low-density intergalactic gas and DM halos at high redshift. DESI will obtain radial velocities of 10 million stars in the Milky Way (MW) and Local Group satellites enabling us to constrain their global DM distributions, as well as the DM distribution on smaller scales. The paradigm of cosmological structure formation has been extensively tested with simulations. However, the majority of simulations to date have focused on collisionless CDM. Simulations with alternatives to CDM have recently been gaining ground but are still in their infancy. While there are numerous publicly available large-box and zoom-in simulations in the LCDM framework, there are no comparable publicly available WDM, SIDM, FDM simulations. DOE support for a public simulation suite will enable a more cohesive community effort to compare observations from DESI (and other surveys) with numerical predictions and will greatly impact DM science.

preprint2021arXiv

APOGEE Chemical Abundance Patterns of the Massive Milky Way Satellites

The SDSS-IV Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey has obtained high-resolution spectra for thousands of red giant stars distributed among the massive satellite galaxies of the Milky Way (MW): the Large and Small Magellanic Clouds (LMC/SMC), the Sagittarius Dwarf (Sgr), Fornax (Fnx), and the now fully disrupted \emph{Gaia} Sausage/Enceladus (GSE) system. We present and analyze the APOGEE chemical abundance patterns of each galaxy to draw robust conclusions about their star formation histories, by quantifying the relative abundance trends of multiple elements (C, N, O, Mg, Al, Si, Ca, Fe, Ni, and Ce), as well as by fitting chemical evolution models to the [$α$/Fe]-[Fe/H] abundance plane for each galaxy. Results show that the chemical signatures of the starburst in the MCs observed by Nidever et al. in the $α$-element abundances extend to C+N, Al, and Ni, with the major burst in the SMC occurring some 3-4 Gyr before the burst in the LMC. We find that Sgr and Fnx also exhibit chemical abundance patterns suggestive of secondary star formation epochs, but these events were weaker and earlier ($\sim$~5-7 Gyr ago) than those observed in the MCs. There is no chemical evidence of a second starburst in GSE, but this galaxy shows the strongest initial star formation as compared to the other four galaxies. All dwarf galaxies had greater relative contributions of AGB stars to their enrichment than the MW. Comparing and contrasting these chemical patterns highlight the importance of galaxy environment on its chemical evolution.

preprint2021arXiv

The DESI $N$-body Simulation Project I: Testing the Robustness of Simulations for the DESI Dark Time Survey

Analysis of large galaxy surveys requires confidence in the robustness of numerical simulation methods. The simulations are used to construct mock galaxy catalogs to validate data analysis pipelines and identify potential systematics. We compare three $N$-body simulation codes, ABACUS, GADGET, and SWIFT, to investigate the regimes in which their results agree. We run $N$-body simulations at three different mass resolutions, $6.25\times10^{8}$, $2.11\times10^{9}$, and $5.00\times10^{9}~h^{-1}$M$_{\odot}$, matching phases to reduce the noise within the comparisons. We find systematic errors in the halo clustering between different codes are smaller than the DESI statistical error for $s > 20\, h^{-1}$Mpc in the correlation function in redshift space. Through the resolution comparison we find that simulations run with a mass resolution of $2.1\times10^{9}~h^{-1}$M$_{\odot}$ are sufficiently converged for systematic effects in the halo clustering to be smaller than the DESI statistical error at scales larger than $20 \, h^{-1}$Mpc. These findings show that the simulations are robust for extracting cosmological information from large scales which is the key goal of the DESI survey. Comparing matter power spectra, we find the codes agree to within 1% for $k \leq 10~h$Mpc$^{-1}$. We also run a comparison of three initial condition generation codes and find good agreement. In addition, we include a quasi-$N$-body code, FastPM, since we plan use it for certain DESI analyses. The impact of the halo definition and galaxy-halo relation will be presented in a follow up study.

preprint2016arXiv

GARROTXA Cosmological Simulations of Milky Way-sized Galaxies: General Properties, Hot Gas Distribution, and Missing Baryons

We introduce a new set of simulations of Milky Way-sized galaxies using the AMR code ART + hydrodynamics in a $Λ$CDM cosmogony. The simulation series is named GARROTXA and follow the formation of a halo/galaxy from z~$=$~60 to z~$=$~0. The final virial mass of the system is $\sim$7.4$\times$10$^{11}$M$_{\odot}$. Our results are as follows: (a) contrary to many previous studies, the circular velocity curve shows no central peak and overall agrees with recent MW observations. (b) Other quantities, such as M$\_{*}$(6$\times$10$^{10}$M$_\odot$) and R$_d$ (2.56 kpc), fall well inside the observational MW range. (c) We measure the disk-to-total ratio kinematically and find that D/T=0.42. (d) The cold gas fraction and star formation rate (SFR) at z=0, on the other hand, fall short from the values estimated for the Milky Way. As a first scientific exploitation of the simulation series, we study the spatial distribution of the hot X-ray luminous gas. We have found that most of this X-ray emitting gas is in a halo-like distribution accounting for an important fraction but not all of the missing baryons. An important amount of hot gas is also present in filaments. In all our models there is not a massive disk-like hot gas distribution dominating the column density. Our analysis of hot gas mock observations reveals that the homogeneity assumption leads to an overestimation of the total mass by factors 3 to 5 or to an underestimation by factors $0.7-0.1$, depending on the used observational method. Finally, we confirm a clear correlation between the total hot gas mass and the dark matter halo mass of galactic systems.

preprint2016arXiv

The DESI Experiment Part I: Science,Targeting, and Survey Design

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the Ly-$α$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.

preprint2016arXiv

The DESI Experiment Part II: Instrument Design

DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= λ/Δλ$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.

preprint2015arXiv

Detecting Triaxiality in the Galactic Dark Matter Halo through Stellar Kinematics II: Dependence on Dark Matter and Gravity Nature

Recent studies have presented evidence that the Milky Way global potential may be nonspherical. In this case, the assembling process of the Galaxy may have left long lasting stellar halo kinematic fossils due to the shape of the dark matter halo, potentially originated by orbital resonances. We further investigate such possibility, considering now potential models further away from $Λ$CDM halos, like scalar field dark matter halos, MOND, and including several other factors that may mimic the emergence and permanence of kinematic groups, such as, a spherical and triaxial halo with an embedded disk potential. We find that regardless of the density profile (DM nature), kinematic groups only appear in the presence of a triaxial halo potential. For the case of a MOND like gravity theory no kinematic structure is present. We conclude that the detection of these kinematic stellar groups could confirm the predicted triaxiality of dark halos in cosmological galaxy formation scenarios.

preprint2015arXiv

Detection of satellite remnants in the Galactic Halo with Gaia III. Detection limits for Ultra Faint Dwarf Galaxies

We present a method to identify Ultra Faint Dwarf Galaxy (UFDG) candidates in the halo of the Milky Way using the future Gaia catalogue and we explore its detection limits and completeness. The method is based on the Wavelet Transform and searches for over-densities in the combined space of sky coordinates and proper motions, using kinematics in the search for the first time. We test the method with a Gaia mock catalogue that has the Gaia Universe Model Snapshot (GUMS) as a background, and use a library of around 30 000 UFDGs simulated as Plummer spheres with a single stellar population. For the UFDGs we use a wide range of structural and orbital parameters that go beyond the range spanned by real systems, where some UFDGs may remain undetected. We characterize the detection limits as function of the number of observable stars by Gaia in the UFDGs with respect to that of the background and their apparent sizes in the sky and proper motion planes. We find that the addition of proper motions in the search improves considerably the detections compared to a photometric survey at the same magnitude limit. Our experiments suggest that Gaia will be able to detect UFDGs that are similar to some of the known UFDGs even if the limit of Gaia is around 2 magnitudes brighter than that of SDSS, with the advantage of having a full-sky catalogue. We also see that Gaia could even find some UFDGs that have lower surface brightness than the SDSS limit.

preprint2014arXiv

A novel method to bracket the corotation radius in galaxy disks: vertex deviation maps

We map the kinematics of stars in simulated galaxy disks with spiral arms using the velocity ellipsoid vertex deviation (l$_v$). We use test particle simulations, and for the first time, fully self-consistent high resolution N-body models. We compare our maps with the Tight Winding Approximation model analytical predictions. We see that for all barred models spiral arms rotate closely to a rigid body manner and the vertex deviation values correlate with the density peaks position bounded by overdense and underdense regions. In such cases, vertex deviation sign changes from negative to positive when crossing the spiral arms in the direction of disk rotation, in regions where the spiral arms are in between corotation (CR) and the Outer Lindblad Resonance (OLR). By contrast, when the arm sections are inside the CR and outside the OLR, l$_v$ changes from negative to positive.We propose that measurements of the vertex deviations pattern can be used to trace the position of the main resonances of the spiral arms. We propose that this technique might exploit future data from Gaia and APOGEE surveys. For unbarred N-body simulations with spiral arms corotating with disk material at all radii, our analysis suggests that no clear correlation exists between l$_v$ and density structures.

preprint2014arXiv

Supermassive Black Hole Binaries: Environment and Galaxy Host Properties of PTA and eLISA sources

Supermassive black hole (BH) binaries would comprise the strongest sources of gravitational waves (GW) once they reach <<1 pc separations, for both pulsar timing arrays (PTAs) and space based (SB) detectors. While BH binaries coalescences constitute a natural outcome of the cosmological standard model and galaxy mergers, their dynamical evolution is still poorly understood and therefore their abundances at different stages. We use a dynamical model for the decay of BH binaries coupled with a cosmological simulation and semi-empirical approaches to the occupation of haloes by galaxies and BHs, in order to follow the evolution of the properties distribution of galaxies hosting BH binaries candidates to decay due to GWs emission. Our models allow us to relax simplifying hypothesis about the binaries occupation in galaxies and their mass, as well as redshift evolution. Following previously proposed electromagnetic (EM) signatures of binaries in the subpc regime, that include spectral features and variability, we model possible distributions of such signatures and also set upper limits to their lifespan. We found a bimodal distribution of hosts properties, corresponding to BH binaries suitable to be detected by PTA and the ones detectable only from space missions, as eLISA. Although it has been discussed that the peak of eLISA sources may happen at high z, we show that there must be a population of such sources in the nearby Universe that might show detectable EM signatures, representing an important laboratory for multimessenger astrophysics. We found a weak dependence of galaxy host properties on the binaries occupation, that can be traced back to the BH origin. The combination of the host correlations reported here with the expected EM signal, may be helpful to verify the presence of nearby GW candidates, and to distinguish them from 'regular' intrinsic AGN variability.

preprint2013arXiv

Constraining dark matter sub-structure with the dynamics of astrophysical systems

The accuracy of the measurements of some astrophysical dynamical systems allows to constrain the existence of incredibly small gravitational perturbations. In particular, the internal Solar System dynamics (planets, Earth-Moon) opens up the possibility, for the first time, to prove the abundance, mass and size, of dark sub-structures at the Earth vicinity. We find that adopting the standard dark matter density, its local distribution can be composed by sub-solar mass halos with no currently measurable dynamical consequences, regardless of the mini-halo fraction. On the other hand, it is possible to exclude the presence of dark streams with linear mass densities higher than $λ_{\rm st}> 10^{-10} \Msun/\AU$ (about the Earth mass spread along the diameter of the SS up to the Kuiper belt). In addition, we review the dynamics of wide binaries inside the dwarf spheroidal galaxies in the MW. The dynamics of such kind of binaries seem to be compatible with the presence of a huge fraction of dark sub-structure, thus their existence is not a sharp discriminant of the dark matter hypothesis as been claimed before. However, there are regimes where the constraints from different astrophysical systems may reveal the sub-structure mass function cut-off scale.

preprint2013arXiv

On the galaxy spiral arms' nature as revealed by rotation frequencies

High resolution N-body simulations using different codes and initial condition techniques reveal two different behaviours for the rotation frequency of transient spiral arms like structures. Whereas unbarred disks present spiral arms nearly corotatingwith disk particles, strong barred models (bulged or bulge-less) quickly develop a bar-spiral structure dominant in density, with a pattern speed almost constant in radius. As the bar strength decreases the arm departs from bar rigid rotation and behaves similar to the unbarred case. In strong barred models we detect in the frequency space other subdominant and slower modes at large radii, in agreement with previous studies, however we also detect them in the configuration space. We propose that the distinctive behaviour of the dominant spiral modes can be exploited in order to constraint the nature of Galactic spiral arms by the astrometric survey GAIA and by 2-D spectroscopic surveys like CALIFA and MANGA in external galaxies.

preprint2012arXiv

Detecting Triaxiality in the Galactic Dark Matter halo through Stellar Kinematics

Assuming the dark matter halo of the Milky Way as a non-spherical potential (i.e. triaxial, prolate, oblate), we show how the assembling process of the Milky Way halo, may have left long lasting stellar halo kinematic fossils only due to the shape of the dark matter halo. In contrast with tidal streams, associated with recent satellite accretion events, these stellar kinematic groups will typically show inhomogeneous chemical and stellar population properties. However, they may be dominated by a single accretion event for certain mass assembling histories. If the detection of these peculiar kinematic stellar groups is confirmed, they would be the smoking gun for the predicted triaxiality of dark halos in cosmological galaxy formation scenarios.

preprint2012arXiv

Hints on halo evolution in SFDM models with galaxy observations

A massive, self-interacting scalar field has been considered as a possible candidate for the dark matter in the universe. We present an observational constraint to the model arising from strong lensing observations in galaxies. The result points to a discrepancy in the properties of scalar field dark matter halos for dwarf and lens galaxies, mainly because halo parameters are directly related to physical quantities in the model. This is an important indication that it becomes necessary to have a better understanding of halo evolution in scalar field dark matter models, where the presence of baryons can play an important role.

preprint2011arXiv

Detection of Satellite Remnants in the Galactic Halo with Gaia - II. A modified Great Circle Cell Method

We propose an extension of the GC3 streamer finding method of Johnston et al. 1996 that can be applied to the future Gaia database. The original method looks for streamers along great circles in the sky, our extension adds the kinematical restriction that velocity vectors should also be constrained to lie along these great circles, as seen by a Galactocentric observer. We show how to use these combined criteria starting from heliocentric observables. We test it by using the mock Gaia catalogue of Brown et al. 2005, which includes a realistic Galactic background and observational errors, but with the addition of detailed star formation histories for the simulated satellites. We investigate its success rate as a function of initial satellite luminosity, star formation history and orbit. We find that the inclusion of the kinematical restriction vastly enhances the contrast between a streamer and the background, even in the presence of observational errors, provided we use only data with good astrometric quality (fractional errors of 30 per cent or better). The global nature of the method diminishes the erasing effect of phase mixing and permits the recovery of merger events of reasonable dynamical age. Satellites with a star formation history different to that of the Galactic background are also better isolated. We find that satellites in the range of 10^8-10^9 Lsun can be recovered even for events as old as ~10 Gyr. Even satellites with 4-5x10^7 Lsun can be recovered for certain combinations of dynamical ages and orbits.

preprint2010arXiv

Low-mass galaxy formation in cosmological AMR simulations: the effects of varying the sub-grid physics parameters

We present numerical simulations aimed at exploring the effects of varying the sub-grid physics parameters on the evolution and the properties of the galaxy formed in a low-mass dark matter halo (~7 times 10^10 Msun/h at redshift z=0). The simulations are run within a cosmological setting with a nominal resolution of 218 pc comoving and are stopped at z = 0.43. In all of our simulations, an extended old/intermediate-age stellar halo and a more compact younger stellar disk are formed. We found that a non negligible fraction of the halo stars are formed in situ in a spheroidal distribution. Changes in the sub-grid physics parameters affect significantly and in a complex way the evolution and properties of the galaxy: (i) Lower threshold densities nsf produce larger stellar effective radii Re, less peaked circular velocity curves V_c(R), and greater amounts of low-density and hot gas in the disk mid-plane; (ii) When stellar feedback is modeled by temporarily switching off radiative cooling in the star forming regions, Re increases (by a factor of ~ 2 in our particular model), the circular velocity curve becomes flatter, and a complex multi-phase gaseous disk structure develops; (iii) A more efficient local conversion of gas mass to stars, measured by a stellar particle mass distribution biased toward larger values, increases the strength of the feedback energy injection -driving outflows and inducing burstier SF histories; iv) If feedback is too strong, gas loss by galactic outflows -which are easier to produce in low-mass galaxies- interrupts SF, whose history becomes episodic. The simulations exhibit two important shortcomings: the baryon fractions are higher, and the specific SF rates are much smaller, than observationally inferred values for redshifts ~ 0.4-1.