No Arabic abstract
We present a significantly improved version of our numerical code JASMINE, that can now solve the Jeans equations for axisymmetric models of stellar systems, composed of an arbitrary number of stellar populations, a Dark Matter halo, and a central Black Hole. The stellar components can have different structural (density profile, flattening, mass, scale length), dynamical (rotational support, velocity dispersion anisotropy), and population (age, metallicity, Initial Mass Function, mass-to-light ratio) properties. These models, when combined with observations, will allow to investigate important issues, such as quantifying the systematic effects of IMF variations, of mass-to-light ratio gradients, and of different stellar kinematic components (e.g. counter rotating disks, kinematically decoupled cores) on luminosity-weighted properties. The developed analytical and numerical framework aims at modeling Early-Type Galaxies, but it can also be applied to dwarf Spheroidal galaxies and Globular Clusters.
Galaxies are self-gravitating structures composed by several components encompassing spherical, axial and triaxial symmetry. Although real systems feature heterogeneous components whose properties are intimately connected, semi-analytical approaches often exploit the linearity of the Poissons equation to represent the potential and mass distribution of a multi-component galaxy as the sum of the individual components. In this work, we expand the semi-analytical framework developed in Bonetti et al. (2020) by including both a detailed implementation of the gravitational potential of exponential disc (modelled with a ${rm sech}^2$ and an exponential vertical profile) and an accurate prescription for the dynamical friction experienced by massive perturbers in composite galaxy models featuring rotating disc structures. Such improvements allow us to evolve arbitrary orbits either within or outside the galactic disc plane. We validate the results obtained by our numerical model against public semi-analytical codes as well as full N-body simulations, finding that our model is in excellent agreement to the codes it is compared with. The ability to reproduce the relevant physical processes responsible for the evolution of massive perturber orbits and its computational efficiency make our framework perfectly suited for large parameter-space exploration studies.
We present new dynamical models of dwarf spheroidal galaxies (dSphs) in which both the stellar component and the dark halo are described by analytic distribution functions that depend on the action integrals. In their most general form these distribution functions can represent axisymmetric and possibly rotating stellar systems. Here, as a first application, we model the Fornax dSph, limiting ourselves, for simplicity, to the non rotating, spherical case. The models are compared with state-of-the-art spectroscopic and photometric observations of Fornax, exploiting the knowledge of the line-of-sight velocity distribution of the models and accounting for the foreground contamination from the Milky Way. The model that best fits the structural and kinematic properties of Fornax has a cored dark halo, with core size $r_{rm c}simeq1.03$ kpc. The dark-to-luminous mass ratio is $(M_{rm dm}/M_{star})|_{R_{rm eff}}simeq9.6$ within the effective radius $R_{rm eff} simeq 0.62,$kpc and $(M_{rm dm}/M_{star})|_{3 {rm kpc}} simeq 144$ within 3 kpc. The stellar velocity distribution is isotropic almost over the full radial range covered by the spectroscopic data and slightly radially anisotropic in the outskirts of the stellar distribution. The dark-matter annihilation $J$-factor and decay $D$-factor are, respectively, $log_{10}(J$ $[$GeV$^2$ cm$^{-5}])simeq18.34$ and $log_{10}(D$ $[$GeV cm$^{-2}])simeq18.55$, for integration angle $theta = 0.5^{circ}$. This cored halo model of Fornax is preferred, with high statistical significance, to both models with a Navarro, Frenk and White dark halo and simple mass-follows-light models.
Adaptive optics observations of the flattened nuclear star cluster in the nearby edge-on spiral galaxy NGC 4244 using the Gemini Near-Infrared Integral Field Spectrograph (NIFS) have revealed clear rotation. Using these kinematics plus 2MASS photometry we construct a series of axisymmetric two-component particle dynamical models with our improved version of NMAGIC, a flexible Chi^2-made-to-measure code. The models consist of a nuclear cluster disc embedded within a spheroidal particle population. We find a mass for the nuclear star cluster of M=1.6^+0.5_-0.2 x 10^7 M_sun within ~42.4 pc (2). We also explore the presence of an intermediate mass black hole and show that models with a black hole as massive as M_bh = 5.0 x 10^5 M_sun are consistent with the available data. Regardless of whether a black hole is present or not, the nuclear cluster is vertically anisotropic (beta_z < 0), as was found with earlier two-integral models. We then use the models as initial conditions for N-body simulations. These simulations show that the nuclear star cluster is stable against non-axisymmetric perturbations. We also explore the effect of the nuclear cluster accreting star clusters at various inclinations. Accretion of a star cluster with mass 13% that of the nuclear cluster is already enough to destroy the vertical anisotropy, regardless of orbital inclination.
We present a large sample of fully self-consistent hydrodynamical Nbody/Tree-SPH simulations of isolated dwarf spheroidal galaxies (dSphs). It has enabled us to identify the key physical parameters and mechanisms at the origin of the observed variety in the Local Group dSph properties. The initial total mass (gas + dark matter) of these galaxies is the main driver of their evolution. Star formation (SF) occurs in series of short bursts. In massive systems, the very short intervals between the SF peaks mimic a continuous star formation rate, while less massive systems exhibit well separated SF bursts, as identified observationally. The delay between the SF events is controlled by the gas cooling time dependence on galaxy mass. The observed global scaling relations, luminosity-mass and luminosity-metallicity, are reproduced with low scatter. We take advantage of the unprecedentedly large sample size and data homogeneity of the ESO Large Programme DART, and add to it a few independent studies, to constrain the star formation history of five Milky Way dSphs, Sextans, LeoII, Carina, Sculptor and Fornax. For the first time, [Mg/Fe] vs [Fe/H] diagrams derived from high-resolution spectroscopy of hundreds of individual stars are confronted with model predictions. We find that the diversity in dSph properties may well result from intrinsic evolution. We note, however, that the presence of gas in the final state of our simulations, of the order of what is observed in dwarf irregulars, calls for removal by external processes.
We use the Chandrasekhar formalism and direct N-body simulations to study the effect of dynamical friction on a test object only slightly more massive than the field stars, orbiting a spherically symmetric background of particles with a mass spectrum. The main goal is to verify whether the dynamical friction time (t_DF) develops a non-monotonic radial-dependence that could explain the bimodality of the Blue Straggler radial distributions observed in globular clusters. In these systems, in fact, relaxation effects lead to a mass and velocity radial segregation of the different mass components, so that mass-spectrum effects on t_DF are expected to be dependent on radius. We find that, in spite of the presence of different masses, t_DF is always a monotonic function of radius, at all evolutionary times and independently of the initial concentration of the simulated cluster. This because the radial dependence of t_DF is largely dominated by the total mass density profile of the background stars (which is monotonically decreasing with radius). Hence, a progressive temporal erosion of the BSS population at larger and larger distances from the cluster center remains the simplest and the most likely explanation of the shape of the observed BSS radial distributions, as suggested in previous works. We also confirm the theoretical expectation that approximating a multi-mass globular cluster as made of (averaged) equal-mass stars can lead to significant overestimates of t_DF within the half-mass radius.