Showing 1–27 of 27 results for author: Di Cintio, P

Tidal mass loss in the Fornax dwarf spheroidal galaxy through N-body simulations with Gaia EDR3-based orbits

Authors: Pierfrancesco Di Cintio, Giuliano Iorio, Francesco Calura, Carlo Nipoti, Marcello Cantari

Abstract: The Fornax dwarf spheroidal galaxy (dSph) represents a challenge for some globular cluster (GC) formation models, because an exceptionally high fraction of its stellar mass is locked in its GC system. In order to shed light on our understanding of GC formation, we aim to constrain the amount of stellar mass that Fornax has lost via tidal interaction with the Milky Way (MW). Exploiting the flexibil… ▽ More The Fornax dwarf spheroidal galaxy (dSph) represents a challenge for some globular cluster (GC) formation models, because an exceptionally high fraction of its stellar mass is locked in its GC system. In order to shed light on our understanding of GC formation, we aim to constrain the amount of stellar mass that Fornax has lost via tidal interaction with the Milky Way (MW). Exploiting the flexibility of effective multi-component $N$-body simulations and relying on state-of-the-art estimates of Fornax's orbital parameters, we study the evolution of the mass distribution of the Fornax dSph in observationally justified orbits in the gravitational potential of the MW over 12 Gyr. We find that, though the dark-matter mass loss can be substantial, the fraction of stellar mass lost by Fornax to the MW is always negligible, even in the most eccentric orbit considered. We conclude that stellar-mass loss due to tidal stripping is not a plausible explanation for the unexpectedly high stellar mass of the GC system of the Fornax dSph and we discuss quantitatively the implications for GC formation scenarios. △ Less

Submitted 1 July, 2024; originally announced July 2024.

Comments: 13 pages, 11 figures, Submitted. Comments welcome

Dynamical friction in the quasi-linear formulation of modified Newtonian dynamics (QuMOND)

Authors: Pierfrancesco Di Cintio, Federico Re, Caterina Chiari

Abstract: Aims. We explore the dynamical friction on a test mass in gravitational systems in the Quasi linear formulation of Modified Newtonian Dynamics (QuMOND). Methods. Exploiting the quasi linearity of QuMOND we derive a simple expression for the dynamical friction in akin to its Newtonian counterpart in the standard Chandrasekhar derivation. Moreover, adopting a mean field approach based on the Liouvil… ▽ More Aims. We explore the dynamical friction on a test mass in gravitational systems in the Quasi linear formulation of Modified Newtonian Dynamics (QuMOND). Methods. Exploiting the quasi linearity of QuMOND we derive a simple expression for the dynamical friction in akin to its Newtonian counterpart in the standard Chandrasekhar derivation. Moreover, adopting a mean field approach based on the Liouville equation we obtain a more rigorous (though in integral form) dynamical friction formula that can be evaluated numerically for a given choice of the QuMOND interpolation function. Results. Consistently with previous work, we observe that dynamical friction is stronger in MOND with respect to a baryon only Newtonian system with the same mass distribution. This amounts to a correction of the Coulomb logarithmic factor via extra terms proportional to the MOND radius of the system. Moreover, with the aid of simple numerical experiments we confirm our theoretical predictions and those of previous work on MOND. △ Less

Submitted 26 June, 2024; v1 submitted 29 April, 2024; originally announced April 2024.

Comments: 12 pages, 3 Figures. Version accepted for publication in A&A

Sharpening the dark matter signature in gravitational waveforms II: Numerical simulations with the NbodyIMRI code

Authors: Bradley J. Kavanagh, Theophanes K. Karydas, Gianfranco Bertone, Pierfrancesco Di Cintio, Mario Pasquato

Abstract: Future gravitational wave observatories can probe dark matter by detecting the dephasing in the waveform of binary black hole mergers induced by dark matter overdensities. Such a detection hinges on the accurate modelling of the dynamical friction, induced by dark matter on the secondary compact object in intermediate and extreme mass ratio inspirals. In this paper, we introduce NbodyIMRI, a new p… ▽ More Future gravitational wave observatories can probe dark matter by detecting the dephasing in the waveform of binary black hole mergers induced by dark matter overdensities. Such a detection hinges on the accurate modelling of the dynamical friction, induced by dark matter on the secondary compact object in intermediate and extreme mass ratio inspirals. In this paper, we introduce NbodyIMRI, a new publicly available code designed for simulating binary systems within cold dark matter `spikes'. Leveraging higher particle counts and finer timesteps, we validate the applicability of the standard dynamical friction formalism and provide an accurate determination of the maximum impact parameter of particles which can effectively scatter with a compact object, across various mass ratios. We also show that in addition to feedback due to dynamical friction, the dark matter also evolves through a `stirring' effect driven by the time-dependent potential of the binary. We introduce a simple semi-analytical scheme to account for this effect and demonstrate that including stirring tends to slow the rate of dark matter depletion and therefore enhances the impact of dark matter on the dynamics of the binary. △ Less

Submitted 21 February, 2024; originally announced February 2024.

Comments: 11 pages, 9 figures + appendices. Companion paper to "Sharpening the dark matter signature in gravitational waveforms I: Accretion and eccentricity evolution" (arXiv:2402.13053). NbodyIMRI code available here: https://github.com/bradkav/NbodyIMRI

Temperature inversion in a confined plasma atmosphere: coarse-grained effect of temperature fluctuations at its base

Authors: Luca Barbieri, Emanuele Papini, Pierfrancesco Di Cintio, Simone Landi, Andrea Verdini, Lapo Casetti

Abstract: Prompted by the relevant problem of temperature inversion (i.e. gradient of density anti-correlated to the gradient of temperature) in astrophysics, we introduce a novel method to model a gravitationally confined multi-component collisionless plasma in contact with a fluctuating thermal boundary. We focus on systems with anti-correlated (inverted) density and temperature profiles, with application… ▽ More Prompted by the relevant problem of temperature inversion (i.e. gradient of density anti-correlated to the gradient of temperature) in astrophysics, we introduce a novel method to model a gravitationally confined multi-component collisionless plasma in contact with a fluctuating thermal boundary. We focus on systems with anti-correlated (inverted) density and temperature profiles, with applications to solar physics. The dynamics of the plasma is analytically described via the coupling of an appropriated coarse-grained distribution function and temporally coarse-grained Vlasov dynamics. We derive a stationary solution of the system and predict the inverted density and temperature profiles of the two-species for scenarios relevant for the corona. We validate our method by comparing the analytical results with kinetic numerical simulations of the plasma dynamics in the context of the two-species Hamiltonian mean-field model (HMF). Finally, we apply our theoretical framework to the problem of the temperature inversion in the solar corona obtaining density and temperature profiles in remarkably good agreement with the observations. △ Less

Submitted 18 June, 2024; v1 submitted 19 January, 2024; originally announced January 2024.

Comments: 27 pages, 12 Figs. Version accepted for publication in JPP

The search for the lost attractor

Authors: Mario Pasquato, Syphax Haddad, Pierfrancesco Di Cintio, Alexandre Adam, Pablo Lemos, Noé Dia, Mircea Petrache, Ugo Niccolò Di Carlo, Alessandro Alberto Trani, Laurence Perreault-Levasseur, Yashar Hezaveh

Abstract: N-body systems characterized by inverse square attractive forces may display a self similar collapse known as the gravo-thermal catastrophe. In star clusters, collapse is halted by binary stars, and a large fraction of Milky Way clusters may have already reached this phase. It has been speculated -- with guidance from simulations -- that macroscopic variables such as central density and velocity d… ▽ More N-body systems characterized by inverse square attractive forces may display a self similar collapse known as the gravo-thermal catastrophe. In star clusters, collapse is halted by binary stars, and a large fraction of Milky Way clusters may have already reached this phase. It has been speculated -- with guidance from simulations -- that macroscopic variables such as central density and velocity dispersion are governed post-collapse by an effective, low-dimensional system of ODEs. It is still hard to distinguish chaotic, low dimensional motion, from high dimensional stochastic noise. Here we apply three machine learning tools to state-of-the-art dynamical simulations to constrain the post collapse dynamics: topological data analysis (TDA) on a lag embedding of the relevant time series, Sparse Identification of Nonlinear Dynamics (SINDY), and Tests of Accuracy with Random Points (TARP). △ Less

Submitted 27 November, 2023; originally announced November 2023.

Comments: Accepted by ML4PS workshop at NeurIPS 2023

Dissipationless collapse and the dynamical mass-ellipticity relation of elliptical galaxies in Newtonian gravity and MOND

Authors: Pierfrancesco Di Cintio

Abstract: Context. Deur (2014) and Winters et al. (2023) proposed an empirical relation between the dark to total mass ratio and ellipticity in elliptical galaxies from their observed total dynamical mass-to-light ratio data M/L = (14.1 +/- 5.4)ε. In other words, the larger is the content of dark matter in the galaxy, the more the stellar component would be flattened. Such observational claim, if true, appe… ▽ More Context. Deur (2014) and Winters et al. (2023) proposed an empirical relation between the dark to total mass ratio and ellipticity in elliptical galaxies from their observed total dynamical mass-to-light ratio data M/L = (14.1 +/- 5.4)ε. In other words, the larger is the content of dark matter in the galaxy, the more the stellar component would be flattened. Such observational claim, if true, appears to be in stark contrast with the common intuition of the formation of galaxies inside dark halos with reasonably spherical symmetry. Aims. Comparing the processes of dissipationless galaxy formation in different theories of gravity, and emergence of the galaxy scaling relations therein is an important frame where, in principle one could discriminate them. Methods. By means of collisionless N-body simulations in modified Newtonian dynamics (MOND) and Newtonian gravity with and without active dark matter halos, with both spherical and clumpy initial structure, I study the trends of intrinsic and projected ellipticities, Sérsic index and anisotropy with the total dynamical to stellar mass ratio. Results. It is shown that, the end products of both cold spherical collapses and mergers of smaller clumps depart more and more from the spherical symmetry for increasing values of the total dynamical mass to stellar mass, at least in a range of halo masses. The equivalent Newtonian systems of the end products of MOND collapses show a similar behaviour. The M/L relation obtained from the numerical experiments in both gravities is however rather different from that reported by Deur and coauthors. △ Less

Submitted 6 February, 2024; v1 submitted 18 October, 2023; originally announced October 2023.

Comments: 11 pages, 8 figures, 1 table. Version accepted in A&A

Structure of the equivalent Newtonian systems in MOND N-body simulations. Density profiles and the core-cusp problem

Authors: Federico Re, Pierfrancesco Di Cintio

Abstract: We investigate the core-cusp problem of the $Λ$ cold dark matter ($Λ$CDM) scenario in the context of Modified Newtonian Dynamics (MOND) paradigm exploiting the concept of equivalent Newtonian system (ENS). By means of particle-mesh $N-$body simulations in MOND we explore processes of galaxy formation via cold dissipationless collapse or merging of smaller substructures. From the end states of our… ▽ More We investigate the core-cusp problem of the $Λ$ cold dark matter ($Λ$CDM) scenario in the context of Modified Newtonian Dynamics (MOND) paradigm exploiting the concept of equivalent Newtonian system (ENS). By means of particle-mesh $N-$body simulations in MOND we explore processes of galaxy formation via cold dissipationless collapse or merging of smaller substructures. From the end states of our simulations we recover the associated ENS and study the properties of their dark matter halos. We compare the simulation results with simple analytical estimates with a family of $γ-$models. We find that the dark matter density of ENSs of most spherical cold collapses ha a markedly cored structure, in particular for the lowest values of the initial virial ratios. End states of some simulations with clumpy initial conditions have more complex profiles and some of their ENSs exhibit a moderate cusp, with logarithmic density slope always shallower than 1. These results seem to point towards the fact that the absence in most observed galaxies of a central DM cusp, at variance with what one would expect from theoretical and numerical arguments in $Λ$CDM, would be totally consistent in a MONDian description. △ Less

Submitted 28 August, 2023; v1 submitted 17 July, 2023; originally announced July 2023.

Comments: 10 pages, 8 figures 1 table. Matching the version accepted in A&A

Journal ref: A&A 678, A110 (2023)

Dynamics of intermediate mass black holes in globular clusters. Wander radius and anisotropy profiles

Authors: Pierfrancesco Di Cintio, Mario Pasquato, Luca Barbieri, Alessandro A. Trani, Ugo N. Di Carlo

Abstract: We recently introduced a new method for simulating collisional gravitational N-body systems with approximately linear time scaling with $N$, based on the Multi-Particle Collision (MPC) scheme, previously applied in Plasma Physics. We simulate globular clusters with a realistic number of stellar particles (at least up to several times $10^6$) on a standard workstation. We simulate clusters hosting… ▽ More We recently introduced a new method for simulating collisional gravitational N-body systems with approximately linear time scaling with $N$, based on the Multi-Particle Collision (MPC) scheme, previously applied in Plasma Physics. We simulate globular clusters with a realistic number of stellar particles (at least up to several times $10^6$) on a standard workstation. We simulate clusters hosting an intermediate mass black hole (IMBH), probing a broad range of BH-cluster and BH-average-star mass ratios, unrestricted by the computational constraints affecting direct N-body codes. We use either single mass models or models with a Salpeter mass function, with the IMBH initially sitting at the centre. The force exerted by and on the IMBH is evaluated with a direct scheme. We measure the evolution of the Lagrangian radii and core density and velocity dispersion over time. In addition, we study the evolution of the velocity anisotropy profiles. We find that models with an IMBH undergo core collapse at earlier times, the larger the IMBH mass the shallower, with an approximately constant central density at core collapse. The presence of an IMBH tends to lower the central velocity dispersion. These results hold independently of the mass function. For the models with Salpeter MF we observe that equipartition of kinetic energies is never achieved. Orbital anisotropy at large radii appears driven by energetic escapers on radial orbits. We measure the wander radius. Among the results we obtained, which mostly confirm or extend previously known trends that had been established over the range of parameters accessible to direct N-body simulations, we underline that the leptokurtic nature of the IMBH wander radius distribution might lead to IMBHs presenting as off-center more frequently than expected, with implications on observational IMBH detection. △ Less

Submitted 27 March, 2023; v1 submitted 10 February, 2023; originally announced February 2023.

Comments: 11 pages, 10 figures. Matching the version of the ms accepted in A&A

Journal ref: A&A 673, A8 (2023)

Relativistic dynamical friction in stellar systems

Authors: Caterina Chiari, Pierfrancesco Di Cintio

Abstract: We extend the classical formulation of the dynamical friction effect on a test star by Chandrasekhar to the case of relativistic velocities and velocity distributions also accounting for post-Newtonian corrections to the gravitational force. The original kinetic framework is revised and used to construct a special-relativistic dynamical friction formula where the relative velocities changes in sub… ▽ More We extend the classical formulation of the dynamical friction effect on a test star by Chandrasekhar to the case of relativistic velocities and velocity distributions also accounting for post-Newtonian corrections to the gravitational force. The original kinetic framework is revised and used to construct a special-relativistic dynamical friction formula where the relative velocities changes in subsequent encounters are added up with Lorentz transformation and the velocity distribution of the field stars accounts for relativistic velocities. Furthermore, a simple expression is obtained for systems where the post-Newtonian correction on the gravitational forces become relevant even at non-relativistic particle velocities. Finally, using a linearized Lagrangian we derive another expression for the dynamical friction expression in a more compact form than that of Lee (1969). Comparing our formulation with the classical one, we observe that a given test particle suffers a slightly stronger drag when moving through a distribution of field stars with relativistic velocity distribution. Vice versa, a purely classical treatment of a system where post-Newtonian (PN) corrections should be included, over estimates the effect of dynamical friction at low test particle velocity, regardless of the form of velocity distribution. Finally, a first order PN dynamical friction covariant formulation is less strong than its classical counterpart at small velocities but much higher for large velocities over a broad range of mass ratios △ Less

Submitted 18 July, 2023; v1 submitted 12 July, 2022; originally announced July 2022.

Comments: 15 pages, 3 figures. Version accepted in A&A

Journal ref: A&A 677, A140 (2023)

Multiparticle collision simulations of dense stellar systems and plasmas

Authors: P. Di Cintio, M. Pasquato, L. Barbieri, H. Bufferand, L. Casetti, G. Ciraolo, U. N. Di Carlo, P. Ghendrih, J. P. Gunn, S. Gupta, H. Kim, S. Lepri, R. Livi, A. Simon-Petit, A. A. Trani, S. -J. Yoon

Abstract: We summarize a series of numerical experiments of collisional dynamics in dense stellar systems such as globular clusters (GCs) and in weakly collisional plasmas using a novel simulation technique, the so-called Multi-particle collision (MPC) method, alternative to Fokker-Planck and Monte Carlo approaches. MPC is related to particle-mesh approaches for the computation of self consistent long-range… ▽ More We summarize a series of numerical experiments of collisional dynamics in dense stellar systems such as globular clusters (GCs) and in weakly collisional plasmas using a novel simulation technique, the so-called Multi-particle collision (MPC) method, alternative to Fokker-Planck and Monte Carlo approaches. MPC is related to particle-mesh approaches for the computation of self consistent long-range fields, ensuring that simulation time scales with $N\log N$ in the number of particles, as opposed to $N^2$ for direct $N$-body. The collisional relaxation effects are modelled by computing particle interactions based on a collision operator approach that ensures rigorous conservation of energy and momenta and depends only on particles velocities and cell-based integrated quantities. △ Less

Submitted 11 February, 2022; v1 submitted 12 January, 2022; originally announced January 2022.

Comments: 7 pages, 2 figures. To appear in the proceedings of the 362 IAU symposium "Predictive Power of Computational Astrophysics as a Discovery Tool". Chamonix, France 8-12 Nov. 2021

Symplectic coarse graining approach to the dynamics of spherical self-gravitating systems

Authors: Luca Barbieri, Pierfrancesco Di Cintio, Guido Giachetti, Alicia Simon-Petit, Lapo Casetti

Abstract: We investigate the evolution of the phase-space distribution function around slightly perturbed stationary states and the process of violent relaxation in the context of the dissipationless collapse of an isolated spherical self-gravitating system. By means of the recently introduced symplectic coarse graining technique, we obtain an effective evolution equation that allows us to compute the scali… ▽ More We investigate the evolution of the phase-space distribution function around slightly perturbed stationary states and the process of violent relaxation in the context of the dissipationless collapse of an isolated spherical self-gravitating system. By means of the recently introduced symplectic coarse graining technique, we obtain an effective evolution equation that allows us to compute the scaling of the frequencies around a stationary state, as well as the damping times of Fourier modes of the distribution function, with the magnitude of the Fourier $k-$vectors themselves. We compare our analytical results with $N$-body simulations. △ Less

Submitted 17 February, 2022; v1 submitted 20 December, 2021; originally announced December 2021.

Comments: 17 pages 7 figures. Version accepted for publication in MNRAS

Noise, friction and the radial-orbit instability in anisotropic stellar systems: stochastic N-body simulations

Authors: Pierfrancesco Di Cintio, Lapo Casetti

Abstract: By means of numerical simulations we study the radial-orbit instability in anisotropic self-gravitating $N-$body systems under the effect of noise. We find that the presence of additive or multiplicative noise has a different effect on the onset of the instability, depending on the initial value of the orbital anisotropy. By means of numerical simulations we study the radial-orbit instability in anisotropic self-gravitating $N-$body systems under the effect of noise. We find that the presence of additive or multiplicative noise has a different effect on the onset of the instability, depending on the initial value of the orbital anisotropy. △ Less

Submitted 21 October, 2021; originally announced October 2021.

Comments: 6 pages, 3 figures. To appear in the proceedings of the 364 IAU symposium, "Multi-scale (time and mass) Dynamics of Space Objects" 18-22 October 2021, Iasi Roumania

Hierarchical generative models for star clusters from hydro-dynamical simulations

Authors: Stefano Torniamenti, Mario Pasquato, Pierfrancesco Di Cintio, Alessandro Ballone, Giuliano Iorio, M. Celeste Artale, Michela Mapelli

Abstract: Star formation in molecular clouds is clumpy, hierarchically subclustered. Fractal structure also emerges in hydro-dynamical simulations of star-forming clouds. Simulating the formation of realistic star clusters with hydro-dynamical simulations is a computational challenge, considering that only the statistically averaged results of large batches of simulations are reliable, due to the chaotic na… ▽ More Star formation in molecular clouds is clumpy, hierarchically subclustered. Fractal structure also emerges in hydro-dynamical simulations of star-forming clouds. Simulating the formation of realistic star clusters with hydro-dynamical simulations is a computational challenge, considering that only the statistically averaged results of large batches of simulations are reliable, due to the chaotic nature of the gravitational N-body problem. While large sets of initial conditions for N-body runs can be produced by hydro-dynamical simulations of star formation, this is prohibitively expensive in terms of computational time. Here we address this issue by introducing a new technique for generating many sets of new initial conditions from a given set of star masses, positions and velocities from a hydro-dynamical simulation. We use hierarchical clustering in phase space to learn a tree representation of the spatial and kinematic relations between stars. This constitutes the basis for the random generation of new sets of stars which share the same clustering structure of the original ones but have individually different masses, positions, and velocities. We apply this method to the output of a number of hydro-dynamical star-formation simulations, comparing the generated initial conditions to the original ones through a series of quantitative tests, including comparing mass and velocity distributions and fractal dimension. Finally, we evolve both the original and the generated star clusters using a direct N-body code, obtaining a qualitatively similar evolution. △ Less

Submitted 24 December, 2021; v1 submitted 1 June, 2021; originally announced June 2021.

Comments: Updated version of the manuscript: "Seeing the forest for the trees: hierarchical generative models for star clusters from hydro-dynamical simulations". 15 pages, 15 figures, 2 tables. Comments welcome

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

Authors: Pierfrancesco Di Cintio, Mario Pasquato, Alicia Simon-Petit, Suk-Jin Yoon

Abstract: 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… ▽ More 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. △ Less

Submitted 16 June, 2021; v1 submitted 3 March, 2021; originally announced March 2021.

Comments: 8 pages, 9 figs., 1 tab. Version accepted for publication in A&A

Introducing a new multi-particle collision method for the evolution of dense stellar systems. Crash-test N-body simulations

Authors: Pierfrancesco Di Cintio, Mario Pasquato, Hyunwoo Kim, Suk-Jin Yoon

Abstract: Stellar systems are broadly divided into collisional and non-collisional. The latter are large-N systems with long relaxation timescales and can be simulated disregarding two-body interactions, while either computationally expensive direct N-body simulations or approximate schemes are required to properly model the former. Large globular clusters and nuclear star clusters, with relaxation timescal… ▽ More Stellar systems are broadly divided into collisional and non-collisional. The latter are large-N systems with long relaxation timescales and can be simulated disregarding two-body interactions, while either computationally expensive direct N-body simulations or approximate schemes are required to properly model the former. Large globular clusters and nuclear star clusters, with relaxation timescales of the order of a Hubble time, are small enough to display some collisional behaviour and big enough to be impossible to simulate with direct $N$-body codes and current hardware. We introduce a new method to simulate collisional stellar systems, and validate it by comparison with direct $N$-body codes on small-$N$ simulations. The Multi-Particle collision for Dense stellar systems Code (MPCDSS) is a new code for evolving stellar systems with the Multi-Particle Collision method. Such method amounts to a stochastic collision rule that allows to conserve exactly the energy and momentum over a cluster of particles experiencing the collision. The code complexity scales with $N \log N$ in the number of particles. Unlike Monte-Carlo codes, MPCDSS can easily model asymmetric, non-homogeneous, unrelaxed and rotating systems, while allowing us to follow the orbits of individual stars. We evolve small ($N = 3.2 \times 10^4$) star clusters with MPCDSS and with the direct-summation code NBODY6, finding a similar evolution of key indicators. We then simulate different initial conditions in the $10^4 - 10^6$ star range. MPCDSS bridges the gap between small, collisional systems that can be simulated with direct $N$-body codes and large noncollisional systems. MPCDSS in principle allows us to simulate globular clusters such as Omega Cen and M54 and even the nuclear star cluster, beyond the limits of current direct N-body codes in terms of the number of particles. △ Less

Submitted 18 November, 2020; v1 submitted 29 June, 2020; originally announced June 2020.

Comments: 12 pages, 13 figures. Matching the version accepted for publication in A&A

Journal ref: A&A 649, A24 (2021)

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

Authors: Mario Pasquato, Pierfrancesco Di Cintio

Abstract: 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 t… ▽ More 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. △ Less

Submitted 23 August, 2020; v1 submitted 4 May, 2020; originally announced May 2020.

Comments: 8 pages, 6 figures. Version accepted in Astronomy & Astrophysics

Journal ref: A&A 640, A79 (2020)

Discreteness effects, $N-$body chaos and the onset of radial-orbit instability

Authors: Pierfrancesco Di Cintio, Lapo Casetti

Abstract: We study the stability of a family of spherical equilibrium models of self-gravitating systems, the so-called $γ-$models with Osipkov-Merritt velocity anisotropy, by means of $N-$body simulations. In particular, we analyze the effect of self-consistent $N-$body chaos on the onset of radial-orbit instability (ROI). We find that degree of chaoticity of the system associated to its largest Lyapunov e… ▽ More We study the stability of a family of spherical equilibrium models of self-gravitating systems, the so-called $γ-$models with Osipkov-Merritt velocity anisotropy, by means of $N-$body simulations. In particular, we analyze the effect of self-consistent $N-$body chaos on the onset of radial-orbit instability (ROI). We find that degree of chaoticity of the system associated to its largest Lyapunov exponent $Λ_{\rm max}$ has no appreciable relation with the stability of the model for fixed density profile and different values of radial velocity anisotropy. However, by studying the distribution of the Lyapunov exponents $λ_{\rm m}$ of the individual particles in the single-particle phase space, we find that more anisotropic systems have a larger fraction of orbits with larger $λ_{\rm m}$. △ Less

Submitted 16 March, 2020; v1 submitted 16 December, 2019; originally announced December 2019.

Comments: 9 pages, 7 figures. Matching the version accepted for publication in MNRAS

Brownian motion of supermassive black holes in galaxy cores

Authors: Pierfrancesco Di Cintio, Luca Ciotti, Carlo Nipoti

Abstract: We investigate the dynamics of supermassive black holes (SMBHs) in galactic cores by means of a semi-analytic model based on the Langevin equation, including dynamical friction and stochastic noise accounting for the gravitational interactions with stars. The model is validated against direct $N$-body simulations of intermediate-mass black holes in stellar clusters where a realistic number of part… ▽ More We investigate the dynamics of supermassive black holes (SMBHs) in galactic cores by means of a semi-analytic model based on the Langevin equation, including dynamical friction and stochastic noise accounting for the gravitational interactions with stars. The model is validated against direct $N$-body simulations of intermediate-mass black holes in stellar clusters where a realistic number of particles is accessible. For the galactic case, we find that the SMBH experiences a Brownian-like motion with a typical displacement from the geometric center of the Galaxy of a few parsecs, for system parameters compatible with M87. \keywords{stellar dynamics, black hole physics, methods: n-body simulations, methods: statistical. △ Less

Submitted 12 August, 2019; originally announced August 2019.

Comments: 4 pages 2 figures. To appear in the Proceedings of the 351 IAU Symposium "Star Clusters: From the Milky Way to the Early Universe"

Journal ref: Proc. IAU 14 (2019) 93-96

$N$-body chaos, phase-space transport and relaxation in numerical simulations

Authors: Pierfrancesco Di Cintio, Lapo Casetti

Abstract: Using direct $N$-body simulations of self-gravitating systems we study the dependence of dynamical chaos on the system size $N$. We find that the $N$-body chaos quantified in terms of the largest Lyapunov exponent $Λ_{\rm max}$ decreases with $N$. The values of its inverse (the so-called Lyapunov time $t_λ$) are found to be smaller than the two-body collisional relaxation time but larger than the… ▽ More Using direct $N$-body simulations of self-gravitating systems we study the dependence of dynamical chaos on the system size $N$. We find that the $N$-body chaos quantified in terms of the largest Lyapunov exponent $Λ_{\rm max}$ decreases with $N$. The values of its inverse (the so-called Lyapunov time $t_λ$) are found to be smaller than the two-body collisional relaxation time but larger than the typical violent relaxation time, thus suggesting the existence of another collective time scale connected to many-body chaos. △ Less

Submitted 30 July, 2019; originally announced July 2019.

Comments: 4 pages, 2 figures. To appear in the proceedings of the 351 IAU Symposium "Star Clusters: from the Milky Way to the Early Universe"

$N$-body chaos and the continuum limit in numerical simulations of self-gravitating systems, revisited

Authors: Pierfrancesco Di Cintio, Lapo Casetti

Abstract: We revisit the rôle of discreteness and chaos in the dynamics of self-gravitating systems by means of $N$-body simulations with active and frozen potentials, starting from spherically symmetric stationary states and considering the orbits of single particles in a frozen $N$-body potential as well as the orbits of the system in the full $6N$-dimensional phase space. We also consider the intermediat… ▽ More We revisit the rôle of discreteness and chaos in the dynamics of self-gravitating systems by means of $N$-body simulations with active and frozen potentials, starting from spherically symmetric stationary states and considering the orbits of single particles in a frozen $N$-body potential as well as the orbits of the system in the full $6N$-dimensional phase space. We also consider the intermediate case where a test particle moves in the field generated by $N$ non-interacting particles, which in turn move in a static smooth potential. We investigate the dependence on $N$ and on the softening length of the largest Lyapunov exponent both of single particle orbits and of the full $N$-body system. For single orbits we also study the dependence on the angular momentum and on the energy. Our results confirm the expectation that orbital properties of single orbits in finite-$N$ systems approach those of orbits in smooth potentials in the continuum limit $N \to \infty$ and that the largest Lyapunov exponent of the full $N$-body system does decrease with $N$, for sufficiently large systems with finite softening length. However, single orbits in frozen models and active self-consistent models have different largest Lyapunov exponents and the $N$-dependence of the values in non-trivial, so that the use of frozen $N$-body potentials to gain information on large-$N$ systems or on the continuum limit may be misleading in certain cases. △ Less

Submitted 5 September, 2019; v1 submitted 25 January, 2019; originally announced January 2019.

Comments: 13 pages 16 figures. Version accepted in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society 489, 5876-5888 (2019)

Non-thermal states in models of filaments: a dynamical study

Authors: Pierfrancesco Di Cintio, Shamik Gupta, Lapo Casetti

Abstract: We study the origin of the non-thermal profiles observed in filamentary structures in galactic molecular clouds by means of numerical dynamical simulations. We find that such profiles are intrinsic features of the end products of dissipationless collapse in cylindrical symmetry. Moreover, for sufficiently cold initial conditions, we obtain end states characterized by markedly anticorrelated radial… ▽ More We study the origin of the non-thermal profiles observed in filamentary structures in galactic molecular clouds by means of numerical dynamical simulations. We find that such profiles are intrinsic features of the end products of dissipationless collapse in cylindrical symmetry. Moreover, for sufficiently cold initial conditions, we obtain end states characterized by markedly anticorrelated radial density and temperature profiles. Gravitational, dissipationless dynamics alone is thus sufficient to reproduce, at least qualitatively, many of the properties of the observed non-thermal structures. △ Less

Submitted 21 December, 2017; originally announced December 2017.

Comments: To appear in the proceedings of the international conference "Francesco's legacy - Star Formation in Space and Time" Firenze, Italy June 5-9 2017

Dynamical origin of non-thermal states in galactic filaments

Authors: Pierfrancesco Di Cintio, Shamik Gupta, Lapo Casetti

Abstract: Observations strongly suggest that filaments in galactic molecular clouds are in a non-thermal state. As a simple model of a filament we study a two-dimensional system of self-gravitating point particles by means of numerical simulations of the dynamics, with various methods: direct $N$-body integration of the equations of motion, particle-in-cell simulations and a recently developed numerical sch… ▽ More Observations strongly suggest that filaments in galactic molecular clouds are in a non-thermal state. As a simple model of a filament we study a two-dimensional system of self-gravitating point particles by means of numerical simulations of the dynamics, with various methods: direct $N$-body integration of the equations of motion, particle-in-cell simulations and a recently developed numerical scheme that includes multiparticle collisions in a particle-in-cell approach. Studying the collapse of Gaussian overdensities we find that after the damping of virial oscillations the system settles in a non-thermal steady state whose radial density profile is similar to the observed ones, thus suggesting a dynamical origin of the non-thermal states observed in real filaments. Moreover, for sufficiently cold collapses the density profiles are anticorrelated with the kinetic temperature, i.e., exhibit temperature inversion, again a feature that has been found in some observations of filaments. The same happens in the state reached after a strong perturbation of an initially isothermal cylinder. Finally, we discuss our results in the light of recent findings in other contexts (including non-astrophysical ones) and argue that the same kind of non-thermal states may be observed in any physical system with long-range interactions. △ Less

Submitted 28 December, 2017; v1 submitted 6 June, 2017; originally announced June 2017.

Comments: 11 pages, 8 figures, MNRAS LaTeX, accepted for publication in MNRAS

Journal ref: MNRAS 475, 1137 (2018)

Radially anisotropic systems with $r^{-α}$ forces. II: radial-orbit instability

Authors: Pierfrancesco Di Cintio, Luca Ciotti, Carlo Nipoti

Abstract: We continue to investigate the dynamics of collisionless systems of particles interacting via additive $r^{-α}$ interparticle forces. Here we focus on the dependence of the radial-orbit instability on the force exponent $α$. By means of direct $N$-body simulations we study the stability of equilibrium radially anisotropic Osipkov-Merritt spherical models with Hernquist density profile and with… ▽ More We continue to investigate the dynamics of collisionless systems of particles interacting via additive $r^{-α}$ interparticle forces. Here we focus on the dependence of the radial-orbit instability on the force exponent $α$. By means of direct $N$-body simulations we study the stability of equilibrium radially anisotropic Osipkov-Merritt spherical models with Hernquist density profile and with $1\leqα<3$. We determine, as a function of $α$, the minimum value for stability of the anisotropy radius $r_{as}$ and of the maximum value of the associated stability indicator $ξ_s$. We find that, for decreasing $α$, $r_{as}$ decreases and $ξ_s$ increases, i.e. longer-range forces are more robust against radial-orbit instability. The isotropic systems are found to be stable for all the explored values of $α$. The end products of unstable systems are all markedly triaxial with minor-to-major axial ratio $>0.3$, so they are never flatter than an E7 system. △ Less

Submitted 27 February, 2017; v1 submitted 12 December, 2016; originally announced December 2016.

Comments: 12 pages, 6 figures

Journal ref: Mon Not R Astron Soc 468, 2222 (2017)

Radially anisotropic systems with $r^{-α}$ forces: equilibrium states

Authors: Pierfrancesco Di Cintio, Luca Ciotti, Carlo Nipoti

Abstract: We continue the study of collisionless systems governed by additive $r^{-α}$ interparticle forces by focusing on the influence of the force exponent $α$ on radial orbital anisotropy. In this preparatory work we construct the radially anisotropic Osipkov-Merritt phase-space distribution functions for self-consistent spherical Hernquist models with $r^{-α}$ forces and $1\leqα<3$. The resulting syste… ▽ More We continue the study of collisionless systems governed by additive $r^{-α}$ interparticle forces by focusing on the influence of the force exponent $α$ on radial orbital anisotropy. In this preparatory work we construct the radially anisotropic Osipkov-Merritt phase-space distribution functions for self-consistent spherical Hernquist models with $r^{-α}$ forces and $1\leqα<3$. The resulting systems are isotropic at the center and increasingly dominated by radial orbits at radii larger than the anisotropy radius $r_a$. For radially anisotropic models we determine the minimum value of the anisotropy radius $r_{ac}$ as a function of $α$ for phase-space consistency (such that the phase-space distribution function is nowhere negative for $r_a\geq r_{ac}$). We find that $r_{ac}$ decreases for decreasing $α$, and that the amount of kinetic energy that can be stored in the radial direction relative to that stored in the tangential directions for marginally consistent models increases for decreasing $α$. In particular, we find that isotropic systems are consistent in the explored range of $α$. By means of direct $N$-body simulations we finally verify that the isotropic systems are also stable. △ Less

Submitted 21 May, 2015; originally announced May 2015.

Comments: 15 pages, 6 figures. Accepted for publication in the special issue of Journal of Plasma Physics "Complex Plasma Phenomena in the Laboratory and in the Universe"

Journal ref: J. Plasma Phys. 81 (2015) 495810504

Temperature inversion in long-range interacting systems

Authors: Tarcisio N. Teles, Shamik Gupta, Pierfrancesco Di Cintio, Lapo Casetti

Abstract: Temperature inversions occur in nature, e.g., in the solar corona and in interstellar molecular clouds: somewhat counterintuitively, denser parts of the system are colder than dilute ones. We propose a simple and appealing way to spontaneously generate temperature inversions in systems with long-range interactions, by preparing them in inhomogeneous thermal equilibrium states and then applying an… ▽ More Temperature inversions occur in nature, e.g., in the solar corona and in interstellar molecular clouds: somewhat counterintuitively, denser parts of the system are colder than dilute ones. We propose a simple and appealing way to spontaneously generate temperature inversions in systems with long-range interactions, by preparing them in inhomogeneous thermal equilibrium states and then applying an impulsive perturbation. In similar situations, short-range systems would typically relax to another thermal equilibrium, with uniform temperature profile. By contrast, in long-range systems, the interplay between wave-particle interaction and spatial inhomogeneity drives the system to nonequilibrium stationary states that generically exhibit temperature inversion. We demonstrate this mechanism in a simple mean-field model and in a two-dimensional self-gravitating system. Our work underlines the crucial role the range of interparticle interaction plays in determining the nature of steady states out of thermal equilibrium. △ Less

Submitted 22 July, 2015; v1 submitted 13 February, 2015; originally announced February 2015.

Comments: 5 pages + 6 pages of appendix, 5 figures, REVTeX 4-1. To appear in Physical Review E (Rapid Communications). Appendix will be published online-only as Supplemental Material

Journal ref: Physical Review E (Rapid Communications) 92, 020101(R) (2015)

Relaxation of N-body systems with additive r^-α interparticle forces

Authors: Pierfrancesco Di Cintio, Luca Ciotti, Carlo Nipoti

Abstract: In Newtonian gravity the final states of cold dissipationless collapses are characterized by several structural and dynamical properties remarkably similar to those of observed elliptical galaxies. Are these properties a peculiarity of the Newtonian force or a more general feature of long-range forces? We study this problem by means of $N-$body simulations of dissipationless collapse of systems of… ▽ More In Newtonian gravity the final states of cold dissipationless collapses are characterized by several structural and dynamical properties remarkably similar to those of observed elliptical galaxies. Are these properties a peculiarity of the Newtonian force or a more general feature of long-range forces? We study this problem by means of $N-$body simulations of dissipationless collapse of systems of particles interacting via additive $r^{-α}$ forces. We find that most of the results holding in Newtonian gravity are also valid for $α\neq2$. In particular the end products are triaxial and never flatter than an E7 system, their surface density profiles are well described by the Sérsic law, the global density slope-anisotropy inequality is obeyed, the differential energy distribution is an exponential over a large range of energies (for $α\geq1$), and the pseudo phase-space density is a power law of radius. In addition, we show that the process of virialization takes longer (in units of the system's dynamical time) for decreasing values of $α$, and becomes infinite for $α=-1$ (the harmonic oscillator). This is in agreement with the results of deep-MOND collapses (qualitatively corresponding to $α=1$) and it is due to the fact the force becomes more and more similar to the $α=-1$ case, where as well known no relaxation can happen and the system oscillates forever. △ Less

Submitted 1 March, 2013; v1 submitted 22 January, 2013; originally announced January 2013.

Comments: 13 pages, 8 figures. Final version submitted to MNRAS, new figure added

Relaxation of spherical systems with long-range interactions: a numerical investigation

Authors: P. Di Cintio, L. Ciotti

Abstract: The process of relaxation of a system of particles interacting with long-range forces is relevant to many areas of Physics. For obvious reasons, in Stellar Dynamics much attention has been paid to the case of 1/r^2 force law. However, recently the interest in alternative gravities emerged, and significant differences with respect to Newtonian gravity have been found in relaxation phenomena. Here w… ▽ More The process of relaxation of a system of particles interacting with long-range forces is relevant to many areas of Physics. For obvious reasons, in Stellar Dynamics much attention has been paid to the case of 1/r^2 force law. However, recently the interest in alternative gravities emerged, and significant differences with respect to Newtonian gravity have been found in relaxation phenomena. Here we begin to explore this matter further, by using a numerical model of spherical shells interacting with an 1/r^alpha force law obeying the superposition principle. We find that the virialization and phase-mixing times depend on the exponent alpha, with small values of alpha corresponding to longer relaxation times, similarly to what happens when comparing for N-body simulations in classical gravity and in Modified Newtonian Dynamics. △ Less

Submitted 28 March, 2011; originally announced March 2011.

Comments: 6 pages, 3 figures, accepted in the International Journal of Bifurcation and Chaos