Source author record

Pierfrancesco Di Cintio

Pierfrancesco Di Cintio 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

16works
9topics
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

16 published item(s)

preprint2023arXiv

Modified Thirring model beyond the excluded-volume approximation

Long-range interacting systems may exhibit ensemble inequivalence and can possibly attain equilibrium states under completely open conditions, for which energy, volume and number of particles simultaneously fluctuate. Here we consider a modified version of the Thirring model for self-gravitating systems with attractive and repulsive long-range interactions in which particles are treated as hard spheres in dimension d=1,2,3. Equilibrium states of the model are studied under completely open conditions, in the unconstrained ensemble, by means of both Monte Carlo simulations and analytical methods and are compared with the corresponding states at fixed number of particles, in the isothermal-isobaric ensemble. Our theoretical description is performed for an arbitrary local equation of state, which allows us to examine the system beyond the excluded-volume approximation. The simulations confirm the theoretical prediction of the possible occurrence of first-order phase transitions in the unconstrained ensemble. This work contributes to the understanding of long-range interacting systems exchanging heat, work and matter with the environment.

preprint2022arXiv

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

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.

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.

preprint2021arXiv

Kinetic and hydrodynamic regimes in multi-particle-collision dynamics of a one-dimensional fluid with thermal walls

We study the non-equilibrium steady-states of a one-dimensional ($1D1V$) fluid in a finite space region of length $L$. Particles interact among themselves by multi-particle collisions and are in contact with two thermal-wall heat reservoirs, located at the boundaries of the region. After an initial ballistic regime, we find a crossover from a normal (kinetic) transport regime to an anomalous (hydrodynamic) one, above a characteristic size $L_*$. We argue that $L_*$ is proportional to the cube of the collision time among particles. Motivated by the physics of emissive divertors in fusion plasma, we consider the more general case of thermal walls injecting particles with given average (non-thermal) velocity. For fast and relatively cold particles, short systems fail to establish local equilibrium and display non-Maxwellian distributions of velocities.

preprint2020arXiv

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

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}$.

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

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

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.

preprint2019arXiv

Brownian motion of supermassive black holes in galaxy cores

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.

preprint2016arXiv

Pfaffian-like ground states for bosonic atoms and molecules in one-dimensional optical lattices

We study ground states and elementary excitations of a system of bosonic atoms and diatomic Feshbach molecules trapped in a one-dimensional optical lattice using exact diagonalization and variational Monte Carlo methods. We primarily study the case of an average filling of one boson per site. In agreement with bosonization theory, we show that the ground state of the system in the thermodynamic limit corresponds to the Pfaffian-like state when the system is tuned towards the superfluid-to-Mott insulator quantum phase transition. Our study clarifies the possibility of the creation of exotic Pfaffian-like states in realistic one-dimensional systems. We also present preliminary evidence that such states support non-Abelian anyonic excitations that have potential application for fault-tolerant topological quantum computation.

preprint2015arXiv

Anomalous dynamical scaling in anharmonic chains and plasma models with multiparticle collisions

We study the anomalous dynamical scaling of equilibrium correlations in one dimensional systems. Two different models are compared: the Fermi-Pasta-Ulam chain with cubic and quartic nonlinearity and a gas of point particles interacting stochastically through the multiparticle collision dynamics. For both models -that admit three conservation laws- by means of detailed numerical simulations we verify the predictions of nonlinear fluctuating hydrodynamics for the structure factors of density and energy fluctuations at equilibrium. Despite this, violations of the expected scaling in the currents correlation are found in some regimes, hindering the observation of the asymptotic scaling predicted by the theory. In the case of the gas model this crossover is clearly demonstrated upon changing the coupling constant.

preprint2015arXiv

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

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.

preprint2015arXiv

Temperature inversion in long-range interacting systems

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.

preprint2013arXiv

Proton ejection from molecular hydride clusters exposed to strong X-ray pulses

Clusters consisting of small molecules containing hydrogen do eject fast protons when illuminated by short X-ray pulses. A suitable overall charging of the cluster controlled by the X-ray intensity induces electron migration from the surface to the bulk leading to efficient segregation of the protons and to a globally hindered explosion of the heavy atoms even outside the screened volume. We investigate this peculiar effect systematically along the iso-electronic sequence of methane over ammonia and water to the atomic limit of neon as a reference. In contrast to core-shell systems where the outer shell is sacrificed to reduce radiation damage, the intricate proton dynamics of hydride clusters allows one to keep the entire backbone of heavy atoms intact.

preprint2013arXiv

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

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.