BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

David Izquierdo-Villalba

David Izquierdo-Villalba contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.GAastro-ph.COgr-qcastro-ph.IM

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
8works
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

8 published item(s)

preprint2022arXiv

Disc instability and bar formation: view from the IllustrisTNG simulations

We make use of z = 0 samples of strongly barred and unbarred disc galaxies from the TNG100 and TNG50 cosmological hydrodynamical simulations to assess the performance of the simple disc instability criterion proposed by Efstathiou, Lake & Negroponte (1982) (ELN-criterion). We find that strongly barred galaxies generally assemble earlier, are more star-dominated in their central regions, and have more massive and more compact discs than unbarred galaxies. The ELN-criterion successfully identifies ~75% and ~80% of the strongly barred and the unbarred galaxies, respectively. Strongly barred galaxies that the criterion fails to identify tend to have more extended discs, higher spin values and bars that assembled later than is typical for the bulk of the barred population. The bars in many of these cases appear to be produced by an interaction with a close neighbour (i.e. to be externally triggered) rather than to result from secular growth in the disc. On the other hand, we find that unbarred galaxies misclassified as barred by the ELN-criterion typically have stellar discs similar to those of barred galaxies, although more extended in the vertical direction and less star-dominated in their central regions, possibly reflecting later formation times. In addition, the bulge component of these galaxies is significantly more prominent at early times than in the strongly barred sample. Thus, the ELN-criterion robustly identifies secular bar instabilities in most simulated disc galaxies, but additional environmental criteria are needed to account for interaction-induced bar formation.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Extreme Mass Ratio Inspirals triggered by Massive Black Hole Binaries: from Relativistic Dynamics to Cosmological Rates

Extreme mass ratio inspirals (EMRIs) are compact binary systems characterized by a mass-ratio $q=m/M$ in the range $~10^{-9}-10^{-4}$ and represent primary gravitational wave (GW) sources for the forthcoming Laser Interferometer Space Antenna (LISA). While their standard formation channel involves relaxation processes deflecting compact objects on very low angular momentum orbits around the central massive black hole, a number of alternative formation channels has been proposed, including binary tidal break-up, migration in accretion disks and secular and chaotic dynamics around a massive black hole binary (MBHB). In this work, we take an extensive closer look at this latter scenario, investigating how EMRIs can be triggered by a MBHBs, formed in the aftermath of galaxy mergers. By employing a suite of relativistic three-body simulations, we evaluate the efficiency of EMRI formation for different parameters of the MBHB, assessing the importance of both secular and chaotic dynamics. By modelling the distribution of compact objects in galaxy nuclei, we estimate the resulting EMRI formation rate, finding that EMRI are produced in a sharp burst, with peak rates that are 10-100 times higher than the standard two-body relaxation channel, lasting for 10$^6$--10$^8$ years. By coupling our results with an estimate of the cosmic MBHB merger rate, we finally forecast that LISA could observe ${\cal O}(10)$ EMRIs per year formed by this channel.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Galaxy fields of LISA massive black hole mergers in a simulated Universe

LISA will extend the search for gravitational waves (GWs) at $0.1\,{-}\,100$ mHz where loud signals from coalescing binary black holes of $ 10^4 \,{-}\,10^7\,\rm M_{\odot}$ are expected. Depending on their mass and luminosity distance, the uncertainty in the LISA sky-localization decreases from hundreds of deg$^2$ during the inspiral phase to fractions of a deg$^2$ after the merger. By using the semi-analytical model L-Galaxies applied to the Millennium-I merger trees, we generate a simulated Universe to identify the hosts of $z\,{\leq}\,3$ coalescing binaries with total mass of $3\,{\times}\,10^{5}$, $3\,{\times}\,10^6$ and $3\,{\times}\,10^7\rm M_{\odot}$, and varying mass ratio. We find that, even at the time of merger, the number of galaxies around the LISA sources is too large (${\gtrsim}\,10^2$) to allow direct host identification. However, if an X-ray counterpart is associated to the GW sources at $z\,{<}\,1$, all LISA fields at merger are populated by ${\lesssim}\,10$ AGNs emitting above ${\sim}\, 10^{-17} \, \rm erg\,cm^{-2}\,s^{-1}$. For sources at higher redshifts, the poorer sky-localization causes this number to increase up to ${\sim}\, 10^3$. Archival data from eRosita will allow discarding ${\sim}\, 10\%$ of these AGNs, being too shallow to detect the dim X-ray luminosity of the GW sources. Inspiralling binaries in an active phase with masses ${\lesssim}\,10^6\rm M_{\odot}$ at $z\,{\leq}\,0.3$ can be detected, as early as $10$ hours before the merger, by future X-ray observatories in less than a few minutes. For these systems, ${\lesssim}\,10$ AGNs are within the LISA sky-localization area. Finally, the LISA-Taiji network would guarantee the identification of an X-ray counterpart $10$ hours before merger for all binaries at $z\,{\lesssim}\,1$.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Morphological decomposition of TNG50 galaxies: methodology and catalogue

We present MORDOR (MORphological DecOmposeR, a new algorithm for structural decomposition of simulated galaxies based on stellar kinematics. The code measures the properties of up to five structural components (a thin/cold and a thick/warm disc, a classical and a secular bulge, and a spherical stellar halo), and determines the properties of a stellar bar (if present). A comparison with other algorithms presented in the literature yields overall good agreement, with MORDOR displaying a higher flexibility in correctly decomposing systems and identifying bars in crowded environments (e.g. with ongoing fly-bys, often observable in cosmological simulations). We use MORDOR to analyse galaxies in the TNG50 simulation and find the following: ($i$) the thick disc component undergoes the strongest evolution in the binding energy-circularity plane, as expected when disc galaxies decrease their turbulent-rotational support with cosmic time; ($ii$) smaller galaxies (with stellar mass, $10^{9} \lesssim M_{*} / {\rm M_{\odot}} \leq 5 \times 10^{9}$) undergo a major growth in their disc components after $z\sim 1$, whereas ($iii$) the most massive galaxies ($5 \times 10^{10} < M_{*} / {\rm M_{\odot}} \leq 5\times10^{11}$) evolve toward more spheroidal dominated objects down to $z=0$ due to frequent gravitational interactions with satellites; ($iv$) the fraction of barred galaxies grows rapidly at high redshift and stabilizes below $z\sim 2$, except for the most massive galaxies that show a decrease in the bar occupation fraction at low redshift; ($v$) galaxies with $M_{*} \sim 10^{11}~{\rm M_{\odot}}$ exhibit the highest relative occurrence of bars at $z=0$, in agreement with observational studies. We publicly release MORDOR and the morphological catalogue of TNG50 galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Multiband gravitational wave cosmology with stellar origin black hole binaries

Massive stellar origin black hole binaries (SBHBs), originating from stars above the pair-instability mass gap, are primary candidates for multiband gravitational wave (GW) observations. Here we study the possibility to use them as effective dark standard sirens to constrain cosmological parameters. The long lasting inspiral signal emitted by these systems is accessible by the future $Laser \; Interferometer \; Space \; Antenna$ (LISA), while the late inspiral and merger are eventually detected by third generation ground-based telescopes such as the $Einstein \; Telescope$ (ET). The direct measurement of the luminosity distance and the sky position to the source, together with the inhomogeneous redshift distribution of possible host galaxies, allow us to infer cosmological parameters by probabilistic means. The efficiency of this statistical method relies in high parameter estimation performances. We show that this multiband approach allows a precise determination of the Hubble constant H$_0$ with just ${\cal O}(10)$ detected sources. For selected SBHB population models, assuming $4$ ($10$) years of LISA observations, we find that H$_0$ is typically determined at $\sim 2\%$ ($\sim 1.5\%$), whereas $Ω_m$ is only mildly constrained with a typical precision of $30\%$ ($20\%$). We discuss the origin of some outliers in our final estimates and we comment on ways to reduce their presence.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The buildup of strongly barred galaxies in the TNG100 simulation

We analyse the properties of strongly barred disc galaxies using the TNG100 simulation, a cosmological hydrodynamical realisation of the IllustrisTNG suite. We identify 270 disc galaxies at $z=0$ in the stellar mass range $M_{*}=10^{10.4-11}M_{\odot}$, of which 40 per cent are barred. Of the detected bars, more than half are strong. We find that the fraction of barred galaxies increases with stellar mass, in agreement with observational results. Strongly barred galaxies exhibit, overall, lower gas-to-stellar mass ratios compared to unbarred galaxies. The majority of barred galaxies are quenched (sSFR $\sim10^{-11.7} $yr$^{-1}$), whereas unbarred galaxies continue to be active (sSFR $\sim10^{-10.3}$yr$^{-1}$) on the main sequence of star-forming galaxies. We explore the evolution of strongly barred and unbarred galaxies to investigate their formation and quenching histories. We find that strong bars form between $0.5< z< 1.5$, with more massive galaxies hosting older bars. Strong bars form in galaxies with an early-established prominent disc component, undergoing periods of enhanced star formation and black hole accretion, possibly assisted by cosmological inflows. Unbarred galaxies, on the other hand, assemble most of their mass and disc component at late times. The nuclear region of strongly barred galaxies quenches shortly after bar formation, while unbarred galaxies remain active across time. Our findings are indicative of bar quenching, possibly assisted by nuclear feedback processes. We conclude that the cosmological environment, together with small scale feedback processes, determine the chances of a galaxy to form a bar and to rapidly quench its central region.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The evolution of the barred galaxy population in the TNG50 simulation

We use the magnetic-hydrodynamical simulation TNG50 to study the evolution of barred massive disc galaxies. Massive spiral galaxies are already present as early as $z=4$, and bar formation takes place already at those early times. The bars grow longer and stronger as the host galaxies evolve, with the bar sizes increasing at a pace similar to that of the disc scale lengths. The bar fraction mildly evolves with redshift for galaxies with $M_{*}\geq10^{10}M\odot$, being greater than $\sim40\%$ at $0.5<z<3$ and $\sim30\%$ at $z=0$. When bars larger than a given physical size ($\geq 2\,\rm kpc$) or the angular resolution limit of twice the I-band angular PSF FWHM of the HST are considered, the bar fraction dramatically decreases with increasing redshift, reconciling the theoretical predictions with observational data. We find that barred galaxies have an older stellar population, lower gas fractions and star formation rates than unbarred galaxies. In most cases, the discs of barred galaxies assembled earlier and faster than the discs of unbarred galaxies. We also find that barred galaxies are typical in haloes with larger concentrations and smaller spin parameters than unbarred galaxies. Furthermore, the inner regions of barred galaxies are more baryon-dominated than those of unbarred galaxies but have comparable global stellar mass fractions. Our findings suggest that the bar population could be used as a potential tracer of the buildup of disc galaxies and their host haloes. With this paper, we release a catalogue of barred galaxies in TNG50 at $6$ redshifts between $z=4$ and $z=0$.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

From galactic nuclei to the halo outskirts: tracing supermassive black holes across cosmic history and environments

We study the mass assembly and spin evolution of supermassive black holes (BHs) across cosmic time as well as the impact of gravitational recoil on the population of nuclear and wandering black holes (wBHs) by using the semi-analytical model LGalaxies run on top of Millennium merger trees. We track spin changes that BHs experience during both coalescence events and gas accretion phases. For the latter, we assume that spin changes are coupled with the bulge assembly. This assumption leads to predictions for the median spin values of z=0 BHs that depend on whether they are hosted by pseudobulges, classical bulges or ellipticals, being a ~ 0.9, 0.7 and 0.4, respectively. The outcomes of the model display a good consistency with z<4 quasar luminosity functions and the z=0 BH mass function, spin values and black hole-bulge mass correlation. Regarding the wBHs, we assume that they can originate from both the disruption of satellite galaxies (orphan wBH) and ejections due to gravitational recoils (ejected wBH). The model points to a number density of wBHs that increases with decreasing redshift, although this population is always ~2 dex smaller than the one of nuclear black holes. At all redshifts, wBHs are typically hosted in M_{halo}>10^{13} M_{\odot} and M_{stellar}>10^{10} M_{\odot}, being orphan wBHs the dominant type. Besides, independently of redshift and halo mass, \textit{ejected} wBHs inhabit the central regions (<0.3R_{200}) of the host DM halo, while orphan wBHs linger at larger scales (>0.5R_{200}). Finally, we find that gravitational recoils cause a progressive depletion of nuclear BHs with decreasing redshift and stellar mass. Moreover, ejection events lead to changes in the predicted local black hole-bulge relation, in particular for BHs in pseudobulges, for which the relation is flattened at M_{bulge}>10^{10.2} M{\odot} and the scatter increase up to ~3dex.

0 citations0 reviewsNo code linkedOpen access