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.Comment: 8 pages, 6 figures. Version accepted in Astronomy & Astrophysic

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

We revisit the r\^{o}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.Comment: 13 pages 16 figures. Version accepted in MNRA

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

We continue to investigate the dynamics of collisionless systems of particles interacting via additive $r^{-\alpha}$ interparticle forces. Here we focus on the dependence of the radial-orbit instability on the force exponent $\alpha$. 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\alpha<3$. We determine, as a function of $\alpha$, the minimum value for stability of the anisotropy radius $r_{as}$ and of the maximum value of the associated stability indicator $\xi_s$. We find that, for decreasing $\alpha$, $r_{as}$ decreases and $\xi_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 $\alpha$. 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.Comment: 12 pages, 6 figure

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.Comment: 5 pages, 5 figure

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.Comment: 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 Materia