Researcher profile

Mario Pasquato

Mario Pasquato contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
13works
0followers
5topics
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

13 published item(s)

preprint2022arXiv

3D Morphology of Open Clusters in the Solar Neighborhood with Gaia EDR3 II: Hierarchical Star Formation Revealed by Spatial and Kinematic Substructures

We identify members of 65 open clusters in the solar neighborhood using the machine-learning algorithm StarGO based on Gaia EDR3 data. After adding members of twenty clusters from previous studies (Pang et al. 2021a,b; Li et al. 2021) we obtain 85 clusters, and study their morphology and kinematics. We classify the substructures outside the tidal radius into four categories: filamentary (f1) and fractal (f2) for clusters $<100$ Myr, and halo (h) and tidal-tail (t) for clusters $>100$ Myr. The kinematical substructures of f1-type clusters are elongated; these resemble the disrupted cluster Group X. Kinematic tails are distinct in t-type clusters, especially Pleiades. We identify 29 hierarchical groups in four young regions (Alessi 20, IC 348, LP 2373, LP 2442); ten among these are new. The hierarchical groups form filament networks. Two regions (Alessi 20, LP 2373) exhibit global &#34;orthogonal&#34; expansion (stellar motion perpendicular to the filament), which might cause complete dispersal. Infalling-like flows (stellar motion along the filament) are found in UBC 31 and related hierarchical groups in the IC 348 region. Stellar groups in the LP 2442 region (LP 2442 gp 1-5) are spatially well-mixed but kinematically coherent. A merging process might be ongoing in the LP 2442 subgroups. For younger systems ($\lesssim30$ Myr), the mean axis ratio, cluster mass and half-mass radius tend to increase with age values. These correlations between structural parameters may imply two dynamical processes occurring in the hierarchical formation scenario in young stellar groups: (1) filament dissolution and (2) sub-group mergers.

preprint2022arXiv

Exploring X-ray variability with unsupervised machine learning I. Self-organizing maps applied to XMM-Newton data

XMM-Newton provides unprecedented insight into the X-ray Universe, recording variability information for hundreds of thousands of sources. Manually searching for interesting patterns in light curves is impractical, requiring an automated data-mining approach for the characterization of sources. Straightforward fitting of temporal models to light curves is not a sure way to identify them, especially with noisy data. We used unsupervised machine learning to distill a large data set of light-curve parameters, revealing its clustering structure in preparation for anomaly detection and subsequent searches for specific source behaviors (e.g., flares, eclipses). Self-organizing maps (SOMs) achieve dimensionality reduction and clustering within a single framework. They are a type of artificial neural network trained to approximate the data with a two-dimensional grid of discrete interconnected units, which can later be visualized on the plane. We trained our SOM on temporal-only parameters computed from more than 100,000 detections from the EXTraS catalog. The resulting map reveals that about 2500 most variable sources are clustered based on temporal characteristics. We find distinctive regions of the SOM map associated with flares, eclipses, dips, linear light curves, and others. Each group contains sources that appear similar by eye. We single out a handful of interesting sources for further study. The condensed view of our dataset provided by SOMs allowed us to identify groups of similar sources, speeding up manual characterization by orders of magnitude. Our method also highlights problems with fitting simple temporal models to light curves and can be used to mitigate them to an extent. This will be crucial for fully exploiting the high data volume expected from upcoming X-ray surveys, and may also help with interpreting supervised classification models.

preprint2022arXiv

NIHAO XXVIII: Collateral effects of AGN on dark matter concentration and stellar kinematics

Although active galactic nuclei (AGN) feedback is required in simulations of galaxies to regulate star formation, further downstream effects on the dark matter distribution of the halo and stellar kinematics of the central galaxy can be expected. We combine simulations of galaxies with and without AGN physics from the Numerical Investigation of a Hundred Astrophysical Objects (NIHAO) to investigate the effect of AGN on the dark matter profile and central stellar rotation of the host galaxies. Specifically, we study how the concentration-halo mass ($c-M$) relation and the stellar spin parameter ($λ_R$) are affected by AGN feedback. We find that AGN physics is crucial to reduce the central density of simulated massive ($\gtrsim 10^{12}$ M$_\odot$) galaxies and bring their concentration to agreement with results from the Spitzer Photometry & Accurate Rotation Curves (SPARC) sample. Similarly, AGN feedback has a key role in reproducing the dichotomy between slow and fast rotators as observed by the ATLAS$^{3\text{D}}$ survey. Without star formation suppression due to AGN feedback, the number of fast rotators strongly exceeds the observational constraints. Our study shows that there are several collateral effects that support the importance of AGN feedback in galaxy formation, and these effects can be used to constrain its implementation in numerical simulations.

preprint2021arXiv

Introducing a new multi-particle collision method for the evolution of dense stellar systems II. Core collapse

In a previous paper we introduced a new method for simulating collisional gravitational $N$-body systems with linear time scaling on $N$, based on the Multi-Particle Collision (MPC) approach. This allows us to simulate globular clusters with a realistic number of stellar particles in a matter of hours on a typical workstation. We evolve star clusters containing up to $10^6$ stars to core collapse and beyond. We quantify several aspects of core collapse over multiple realizations and different parameters, while always resolving the cluster core with a realistic number of particles. We run a large set of N-body simulations with our new code. The cluster mass function is a power-law with no stellar evolution, allowing us to clearly measure the effects of the mass spectrum. Leading up to core collapse, we find a power-law relation between the size of the core and the time left to core collapse. Our simulations thus confirm the theoretical self-similar contraction picture but with a dependence on the slope of the mass function. The time of core collapse has a non-monotonic dependence on the slope, which is well fit by a parabola. This holds also for the depth of core collapse and for the dynamical friction timescale of heavy particles. Cluster density profiles at core collapse show a broken power law structure, suggesting that central cusps are a genuine feature of collapsed cores. The core bounces back after collapse, and the inner density slope evolves to an asymptotic value. The presence of an intermediate-mass black hole inhibits core collapse. We confirm and expand on several predictions of star cluster evolution before, during, and after core collapse. Such predictions were based on theoretical calculations or small-size direct $N$-body simulations. Here we put them to the test on MPC simulations with a much larger number of particles, allowing us to resolve the collapsing core.

preprint2020arXiv

An astrophysically motivated ranking criterion for low-latency electromagnetic follow-up of gravitational wave events

We investigate the properties of the host galaxies of compact binary mergers across cosmic time. To this end, we combine population synthesis simulations together with galaxy catalogues from the hydrodynamical cosmological simulation EAGLE to derive the properties of the host galaxies of binary neutron star (BNS), black hole-neutron star (BHNS) and binary black hole (BBH) mergers. Within this framework, we derive the host galaxy probability, i.e., the probability that a galaxy hosts a compact binary coalescence as a function of its stellar mass, star formation rate, $K_s$ magnitude and $B$ magnitude. This quantity is particularly important for low-latency searches of gravitational wave (GW) sources as it provides a way to rank galaxies lying inside the credible region in the sky of a given GW detection, hence reducing the number of viable host candidates. Furthermore, even if no electromagnetic counterpart is detected, the proposed ranking criterion can still be used to classify the galaxies contained in the error box. Our results show that massive galaxies (or equivalently galaxies with a high luminosity in $K_s$ band) have a higher probability of hosting BNS, BHNS, and BBH mergers. We provide the probabilities in a suitable format to be implemented in future low-latency searches.

preprint2020arXiv

Binary black holes in young star clusters: the impact of metallicity

Young star clusters are the most common birth-place of massive stars and are dynamically active environments. Here, we study the formation of black holes (BHs) and binary black holes (BBHs) in young star clusters, by means of 6000 N-body simulations coupled with binary population synthesis. We probe three different stellar metallicities (Z=0.02, 0.002 and 0.0002) and two initial density regimes (density at the half-mass radius $ρ_{\rm h}\ge{}3.4\times10^4$ and $\ge{1.5\times10^2}$ M$_\odot$ pc$^{-3}$ in dense and loose star clusters, respectively). Metal-poor clusters tend to form more massive BHs than metal-rich ones. We find $\sim{}6$, $\sim{}2$, and $<1$% of BHs with mass $m_{\rm BH}>60$ M$_\odot$ at Z=0.0002, 0.002 and 0.02, respectively. In metal-poor clusters, we form intermediate-mass BHs with mass up to $\sim{}320$ M$_\odot$. BBH mergers born via dynamical exchanges (exchanged BBHs) can be more massive than BBH mergers formed from binary evolution: the former (latter) reach total mass up to $\sim{}140$ M$_\odot$ ($\sim{}80$ M$_\odot$). The most massive BBH merger in our simulations has primary mass $\sim{}88$ M$_\odot$, inside the pair-instability mass gap, and a mass ratio of $\sim{}0.5$. Only BBHs born in young star clusters from metal-poor progenitors can match the masses of GW170729, the most massive event in O1 and O2, and those of GW190412, the first unequal-mass merger. We estimate a local BBH merger rate density $\sim{}110$ and $\sim{}55$ Gpc$^{-3}$ yr$^{-1}$, if we assume that all stars form in loose and dense star clusters, respectively.

preprint2020arXiv

Different Fates of Young Star Clusters After Gas Expulsion

We identify structures of the young star cluster NGC 2232 in the solar neighborhood (323.0 pc), and a newly discovered star cluster LP 2439 (289.1 pc). Member candidates are identified using the Gaia DR2 sky position, parallax and proper motion data, by an unsupervised machine learning method, \textsc{StarGO}. Member contamination from the Galactic disk is further removed using the color magnitude diagram. The four identified groups (NGC 2232, LP 2439 and two filamentary structures) of stars are coeval with an age of 25 Myr and were likely formed in the same giant molecular cloud. We correct the distance asymmetry from the parallax error with a Bayesian method. The 3D morphology shows the two spherical distributions of clusters NGC 2232 and LP 2439. Two filamentary structures are spatially and kinematically connected to NGC 2232. Both NGC 2232 and LP 2439 are expanding. The expansion is more significant in LP 2439, generating a loose spatial distribution with shallow volume number and mass density profiles. The expansion is suggested to be mainly driven by gas expulsion. NGC 2232, with 73~percent of the cluster mass bound, is currently experiencing a process of re-virialization, However, LP 2439, with 52 percent cluster mass being unbound, may fully dissolve in the near future. The different survivability traces different dynamical states of NGC 2232 and LP 2439 prior to the onset of gas expulsion. NGC 2232 may have been substructured and subvirial, while LP 2439 may either have been virial/supervirial, or it has experienced a much faster rate of gas removal.

preprint2020arXiv

Taking apart the dynamical clock. Fat-tailed dynamical kicks shape the blue-straggler star bimodality

In globular clusters, blue straggler stars are heavier than the average star, so dynamical friction strongly affects them. The radial distribution of BSS, normalized to a reference population, appears bimodal in a fraction of Galactic GCs, with a density peak in the core, a prominent zone of avoidance at intermediate radii, and again higher density in the outskirts. The zone of avoidance appears to be located at larger radii the more relaxed the host cluster, acting as a sort of dynamical clock. We use a new method to compute the evolution of the BSS radial distribution under dynamical friction and diffusion. We evolve our BSS in the mean cluster potential under dynamical friction plus a random fluctuating force, solving the Langevin equation with the Mannella quasi symplectic scheme. This amounts to a new simulation method which is much faster and simpler than direct N-body codes but retains their main feature: diffusion powered by strong, if infrequent, kicks. We compute the radial distribution of initially unsegregated BSS normalized to a reference population as a function of time. We trace the evolution of its minimum, corresponding to the zone of avoidance. We compare the evolution under kicks extracted from a Gaussian distribution to that obtained using a Holtsmark distribution. The latter is a fat tailed distribution which correctly models the effects of close gravitational encounters. We find that the zone of avoidance moves outwards over time, as expected based on observations, only when using the Holtsmark distribution. Thus the correct representation of near encounters is crucial to reproduce the dynamics of the system. We confirm and extend earlier results that showed how the dynamical clock indicator depends both on dynamical friction and effective diffusion powered by dynamical encounters.

preprint2019arXiv

Analytical solutions for the dynamical clock A+ indicator in a toy model of pure dynamical friction

Blue straggler stars are more massive than the average star in globular clusters, as they originate from the merger of two stars. Consequently, they experience dynamical friction, progressively sinking to the cluster center. Recently, several indicators of the degree of dynamical relaxation of a globular cluster have been proposed, based on the observed radial distribution of blue straggler stars. The most successful is the Alessandrini indicator, or A+ for short, which is the integral of the cumulative distribution of the blue straggler stars minus that of a lighter reference population. A+ correlates with the dynamical age of a cluster both in realistic simulations and in observations. Here I calculate the temporal dependence of the A+ indicator analytically in a simplified model of the evolution of the blue straggler star distribution under dynamical friction only.

preprint2019arXiv

Clustering clusters: unsupervised machine learning on globular cluster structural parameters

Globular Clusters (GCs) have historically been subdivided in either two (disk/halo) or three (disk/inner-halo/outer-halo) groups based on their orbital, chemical and internal physical properties. The qualitative nature of this subdivision makes it impossible to determine whether the natural number of groups is actually two, three, or more. In this paper we use cluster analysis on the $(\log M, \log σ_0, \log R_e, [Fe/H], \log | Z |)$ space to show that the intrinsic number of GC groups is actually either $k=2$ or $k=3$, with the latter being favored albeit non-significantly. In the $k=2$ case, the Partitioning Around Medoids (PAM) clustering algorithm recovers a metal-poor halo GC group and a metal-rich disk GC group. With $k=3$ the three groups can be interpreted as disk/inner-halo/outer-halo families. For each group we obtain a medoid, i.e. a representative element (NGC $6352$, NGC $5986$, and NGC $5466$ for the disk, inner halo, and outer halo respectively), and a measure of how strongly each GC is associated to its group, the so-called silhouette width. Using the latter, we find a correlation with age for both disk and outer halo GCs where the stronger the association of a GC with the disk (outer halo) group, the younger (older) it is.

preprint2019arXiv

Further properties of the dynamical clock A+ indicator in a toy model of pure dynamical friction

The Alessandrini A$+$ indicator is a quantitative measure of star cluster dynamical evolution based on the mass-segregation of blue straggler stars. A$+$ is defined as the integral of the cumulative distribution of blue straggler stars with the radius measured in logarithmic scale, minus a term related to the reference population used. In a companion paper I introduced a simplified model of dynamical friction and calculated the A$+$ indicator analytically. Here I show further properties of the time evolution of A$+$, focusing in particular on the physical interpretation of its time derivative $d{A^+}/dt$. I find that $d{A^+}/dt$ is the mean of the reciprocal dyamical friction timescale, weighted by the density of blue straggler stars. I show that it is non-negative (as expected based on monotonicity) due to the density of blue-straggler stars being non-negative and that, for a dynamical friction timescale that is non-decreasing with radius, $d{A^+}/dt$ is also non-decreasing with time, making the A$+$ indicator a convex function.

preprint2019arXiv

Mass and star formation rate of the host galaxies of compact binary mergers across cosmic time

We investigate the properties of the host galaxies of compact binary mergers across cosmic time, by means of population-synthesis simulations combined with galaxy catalogues from the EAGLE suite. We analyze the merger rate per galaxy of binary neutron stars (BNSs), black hole--neutron star binaries (BHNSs) and binary black holes (BBHs) from redshift zero up to redshift six. The binary merger rate per galaxy strongly correlates with the stellar mass of the host galaxy at any redshift considered here. This correlation is significantly steeper for BNSs than for both BHNSs and BBHs. Moreover, we find that the merger rate per galaxy depends also on host galaxy&#39;s star formation rate and metallicity. We derive a robust fitting formula that relates the merger rate per galaxy with galaxy&#39;s star formation rate, stellar mass and metallicity at different redshifts. The typical masses of the host galaxies increase significantly as redshift decreases, as a consequence of the interplay between delay time distribution of compact binaries and cosmic assembly of galaxies. Finally, we study the evolution of the merger rate density with redshift. At low redshift ($z\leq{}0.1$) early-type galaxies give a larger contribution to the merger rate density than late-type galaxies. This trend reverts at $z\ge{}1$.

preprint2019arXiv

Stellar-Mass Black Holes in Globular Clusters: Dynamical Consequences and Observational Signatures

Sizeable number of stellar-mass black holes (BHs) in globular clusters (GCs) can strongly influence the dynamical evolution and observational properties of their host cluster. Using results from a large set of numerical simulations, we identify the key ingredients needed to sustain a sizeable population of BHs in GCs up to a Hubble time. We find that while BH natal kick prescriptions are essential in determining the initial retention fraction of BHs in GCs, the long-term survival of BHs is determined by the size, initial central density and half-mass relaxation time of the GC. Simulated GC models that contain many BHs are characterized by relatively low central surface brightness, large half-light and core radii values. We also discuss novel ways to compare simulated results with available observational data to identify GCs that are most likely to contain many BHs.