Influence of a stellar cusp on the dynamics of young stellar discs and the origin of the S-stars in the Galactic Centre

Observations of the Galactic Centre show evidence of one or two disc-like structures of very young stars orbiting the central supermassive black hole within a distance of a few 0.1 pc. A number of analyses have been carried out to investigate the dynamical behaviour and consequences of these discs, including disc thickness and eccentricity growth as well as mutual interaction and warping. However, most of these studies have neglected the influence of the stellar cusp surrounding the black hole, which is believed to be one to two orders of magnitude more massive than the disc(s).
By means of N -body integrations using our bhint code, we study the impact of stellar cusps of different compositions. We find that although the presence of a cusp does have an important effect on the evolution of an otherwise isolated flat disc, its influence on the evolution of disc thickness and warping is rather mild in a two-disc configuration. However, we show that the creation of highly eccentric orbits strongly depends on the graininess of the cusp (i.e. the mean and maximum stellar masses). While Chang recently found that full cycles of Kozai resonance are prevented by the presence of an analytic cusp, we show that relaxation processes play an important role in such highly dense regions and support short-term resonances. We thus find that young disc stars on initially circular orbits can achieve high eccentricities by resonant effects also in the presence of a cusp of stellar remnants, yielding a mechanism to create S-stars and hypervelocity stars.
Furthermore, we discuss the underlying initial mass function (IMF) of the young stellar discs and find no definite evidence for a non-canonical IMF.     

Dark stars at the Galactic Centre – the main sequence

In regions of very high dark matter density such as the Galactic Centre, the capture and annihilation of WIMP dark matter by stars has the potential to significantly alter their evolution. We describe the dark stellar evolution code D ark S tars , and present a series of detailed grids of WIMP-influenced stellar models for main-sequence stars. We describe the changes in stellar structure and main-sequence evolution which occur as a function of the rate of energy injection by WIMPs, for masses of  0.3–2.0 M_⊙  and metallicities   Z = 0.0003–0.02  . We show what rates of energy injection can be obtained using realistic orbital parameters for stars at the Galactic Centre, including detailed consideration of the velocity and density profiles of dark matter. Capture and annihilation rates are strongly boosted when stars follow elliptical rather than circular orbits. If there is a spike of dark matter induced by the supermassive black hole at the Galactic Centre, single solar mass stars following orbits with periods as long as 50 yr and eccentricities as low as 0.9 could be significantly affected. Binary systems with similar periods about the Galactic Centre could be affected on even less eccentric orbits. The most striking observational effect of this scenario would be the existence of a binary consisting of a low-mass protostar and a higher mass evolved star. The observation of low-mass stars and/or binaries on such orbits would either provide a detection of WIMP dark matter, or place stringent limits on the combination of the WIMP mass, spin-dependent nuclear-scattering cross-section, halo density and velocity distribution near the Galactic Centre. In some cases, the derived limits on the WIMP mass and spin-dependent nuclear-scattering cross-section would be of comparable sensitivity to current direct-detection experiments.     

The effectiveness of the Kozai mechanism in the Galactic Centre

I examine the effectiveness of Kozai oscillations in the centres of galaxies and in particular the Galactic Centre (GC) using standard techniques from celestial mechanics. In particular, I study the effects of a stellar bulge potential and general relativity on Kozai oscillations, which are induced by stellar discs. Löckmann et al. recently suggested that Kozai oscillations induced by the two young massive stellar discs in the GC drive the orbits of the young stars to large eccentricity  ( e ≈ 1)  . If some of these young eccentric stars are in binaries, they would be disrupted near pericentre, leaving one star in a tight orbit around the central supermassive black hole and producing the S-star population. I find that the spherical stellar bulge suppresses Kozai oscillations, when its enclosed mass inside a test body is of the order of the mass in the stellar disc(s). Since the stellar bulge in the GC is much larger than the stellar discs, Kozai oscillations due to the stellar discs are likely suppressed. Whether Kozai oscillations are induced from other non-spherical components to the potential (e.g. a flattened stellar bulge) is yet to be determined.     

Simulations of the formation of stellar discs in the Galactic Centre via cloud–cloud collisions

Young massive stars in the central parsec of our Galaxy are best explained by star formation within at least one, and possibly two, massive self-gravitating gaseous discs. With help of numerical simulations, we here consider whether the observed population of young stars could have originated from a large angle collision of two massive gaseous clouds at   R ≃ 1 pc  from Sgr A*. In all the simulations performed, the post-collision gas flow forms an inner, nearly circular gaseous disc and one or two eccentric outer filaments, consistent with the observations. Furthermore, the radial stellar mass distribution is always very steep,  Σ_*∝ R ⁻²  , again consistent with the observations. All of our simulations produce discs that are warped by between 30° and 60°, in accordance with the most recent observations. The three-dimensional velocity structure of the stellar distribution is sensitive to initial conditions (e.g. the impact parameter of the clouds) and gas cooling details. For example, the runs in which the inner disc is fed intermittently with material possessing fluctuating angular momentum result in multiple stellar discs with different orbital orientations, contradicting the observed data. In all the cases the amount of gas accreted by our inner boundary condition is large, enough to allow Sgr A* to radiate near its Eddington limit over ∼10⁵ yr. This suggests that a refined model would have physically larger clouds (or a cloud and a disc such as the circumnuclear disc) colliding at a distance of a few parsecs rather than 1 pc as in our simulations.     

Numerical simulations of protostellar encounters — I. Star–disc encounters

It appears that most stars are born in clusters, and that at birth most stars have circumstellar discs which are comparable in size to the separations between the stars. Interactions between neighbouring stars and discs are therefore likely to play a key role in determining disc lifetimes, stellar masses, and the separations and eccentricities of binary orbits. Such interactions may also cause fragmentation of the discs, thereby triggering the formation of additional stars.   We have carried out a series of simulations of star–disc interactions using an SPH code which treats self-gravity, hydrodynamic and viscous forces. We find that interactions between discs and stars provide a mechanism for removing energy from, or adding energy to, the orbits of the stars, and for truncating the discs. However, capture during such encounters is unlikely to be an important binary formation mechanism.   A more significant consequence of such encounters is that they can trigger fragmentation of the disc, via tidally and compressionally induced gravitational instabilities, leading to the formation of additional stars and substellar objects. When the disc spins and stellar orbits are randomly oriented, encounters lead to the formation of new companions to the original star in 20 per cent of encounters. If most encounters are prograde and coplanar, as suggested by simulations of dynamically triggered star formation, then new companions are formed in approximately 50 per cent of encounters.     

Red giant stellar collisions in the Galactic Centre

We show that collisions with stellar-mass black holes can partially explain the absence of bright giant stars in the Galactic Centre, first noted by Genzel et al. We show that the missing objects are low-mass giants and asymptotic giant branch stars in the range  1–3 M_⊙  . Using detailed stellar evolution calculations, we find that to prevent these objects from evolving to become visible in the depleted K bands, we require that they suffer collisions on the red giant branch, and we calculate the fractional envelope mass losses required. Using a combination of smoothed particle hydrodynamic calculations, restricted three-body analysis and Monte Carlo simulations, we compute the expected collision rates between giants and black holes, and between giants and main-sequence stars in the Galactic Centre. We show that collisions can plausibly explain the missing giants in the  10.5 < K < 12  band. However, depleting the brighter  ( K < 10.5)  objects out to the required radius would require a large population of black hole impactors which would in turn deplete the  10.5 < K < 12  giants in a region much larger than is observed. We conclude that collisions with stellar-mass black holes cannot account for the depletion of the very brightest giants, and we use our results to place limits on the population of stellar-mass black holes in the Galactic Centre.     

High-velocity stars from decay of small stellar systems

In this paper we present numerical results on the decay of small stellar systems under different initial conditions (multiplicity 3 ≤  N  ≤ 10, and various mass spectra, initial velocities and initial configurations). The numerical treatment uses the CHAIN1 code (Mikkola &38; Aarseth). Particular attention is paid to the distribution of high-velocity escapers: we define these as stars with velocity above 30 km s⁻¹. These numerical experiments confirm that small N -body systems are dynamically unstable and produce cascades of escapers in the process of their decay. It is shown that the fraction of stars that escape from small dense stellar systems with an escape velocity greater than 30 km s⁻¹ is ∼1 per cent for all systems treated here. This relatively small fraction must be considered in relation to the rate of star formation in the Galaxy in small groups: this could explain some moderately high-velocity stars observed in the Galactic disc and possibly some young stars with relatively high metallicity in the thick disc.     

X-rays from cusps of compact remnants near galactic centres

Compact remnants – stellar mass black holes and neutron stars formed in the inner few parsec of galactic centres are predicted to sink into the central parsec due to dynamical friction on low-mass stars, forming a high concentration cusp. Same physical region may also contain very high-density molecular clouds and accretion discs that are needed to fuel supermassive black hole (SMBH) activity. Here we estimate gas capture rates on to the cusp of stellar remnants, and the resulting X-ray luminosity, as a function of the accretion disc mass. At low disc masses, most compact objects are too dim to be observable, whereas in the high disc case most of them are accreting at their Eddington rates. We find that for low accretion disc masses, compact remnant cusps may be more luminous than the central SMBHs. This 'diffuse' emission may be of importance for local moderately bright active galactic nuclei (AGNs), especially low-luminosity AGNs. We also briefly discuss how this expected emission can be used to put constraints on the black hole cusp near our Galactic Centre.     

Formation of ω Centauri from an ancient nucleated dwarf galaxy in the young Galactic disc

We first present a self-consistent dynamical model in which ω Cen is formed from an ancient nucleated dwarf galaxy merging with the first generation of the Galactic thin disc in a retrograde manner with respect to the Galactic rotation. Our numerical simulations demonstrate that during merging between the Galaxy and the ω Cen host dwarf with   M _B≃−14 mag  and its nucleus mass of  10⁷ M⊙  , the outer stellar envelope of the dwarf is nearly completely stripped, whereas the central nucleus can survive from the tidal stripping because of its compactness. The developed naked nucleus has a very bound retrograde orbit around the young Galactic disc, as observed for ω Cen, with apocentre and pericentre distances of ∼8 and ∼1 kpc, respectively. The Galactic tidal force can induce radial inflow of gas to the centre of the dwarf and consequently triggers moderately strong nuclear starbursts in a repetitive manner. This result implies that efficient nuclear chemical enrichment resulting from the later starbursts can be closely associated with the origin of the observed relatively young and metal-rich stars in ω Cen. Dynamical heating by the ω Cen host can transform the young thin disc into the thick disc during merging.     

An excess of very bright stars in the inner bulge

From an analysis of the stars remaining in central regions of the Galaxy after subtracting those belonging to the disc and the bulge, we deduce that the inner bulge must have an extra young population with respect to the rest of the bulge. It is shown that there is a higher ratio of very bright stars in the central bulge than there is in the outer bulge. This is interpreted as being an additional young component due to the presence of star formation regions near the Galactic Centre which is absent in the outer bulge.     

Metal-rich carbon stars in the Sagittarius dwarf spheroidal galaxy

We present spectroscopic observations from the Spitzer Space Telescope of six carbon-rich asymptotic giant branch (AGB) stars in the Sagittarius dwarf spheroidal galaxy (Sgr dSph) and two foreground Galactic carbon stars. The band strengths of the observed C₂H₂ and SiC features are very similar to those observed in Galactic AGB stars. The metallicities are estimated from an empirical relation between the acetylene optical depth and the strength of the SiC feature. The metallicities are higher than those of the Large Magellanic Cloud, and close to Galactic values. While the high metallicity could imply an age of around 1 Gyr, for the dusty AGB stars, the pulsation periods suggest ages in excess of 2 or 3 Gyr. We fit the spectra of the observed stars using the dusty radiative transfer model and determine their dust mass-loss rates to be in the range  1.0–3.3 × 10⁻⁸ M_⊙ yr⁻¹  . The two Galactic foreground carbon-rich AGB stars are located at the far side of the solar circle, beyond the Galactic Centre. One of these two stars shows the strongest SiC feature in our present Local Group sample.     

Stellar populations in the Canis Major overdensity

We performed a photometric multicolour survey of the core of the Canis Major overdensity at     , reaching   V ∼ 22  and covering  0.3 × 1.0  arcmin². The main aim is to unravel the complex mixture of stellar populations toward this Galactic direction, where in the recent past important signatures of an accretion event have been claimed to be detected. While our previous investigations were based on disjointed pointings aimed at revealing the large-scale structure of the third Galactic Quadrant, we now focus on a complete coverage of a smaller field centred on the Canis Major overdensity. A large wavelength baseline, in the UBVRI bands, allows us to build up a suite of colour–colour and colour–magnitude diagrams, providing a much better diagnostic tool to disentangle the stellar populations of the region. In fact, the simple use of one colour–magnitude diagram, widely employed in all the previous studies defending the existence of the Canis Major galaxy, does not allow one to separate the effects of the different parameters (reddening, age, metallicity and distance) involved in the interpretation of data, forcing to rely on heavy modelling. In agreement with our previous studies, in the same general region of the Milky Way, we recognize a young stellar population compatible with the expected structure and extension of the Local (Orion) and Outer (Norma–Cygnus) spiral arms in the Third Galactic Quadrant. Moreover, we interpret the conspicuous intermediate-age metal-poor population as belonging to the Galactic thick disc, distorted by the effect of strong disc warping at this latitude, and to the Galactic halo.     

Variable accretion and emission from the stellar winds in the Galactic Centre

We present numerical simulations of stellar wind dynamics in the central parsec of the Galactic Centre, studying in particular the accretion of gas on to Sgr A*, the supermassive black hole. Unlike our previous work, here we use state-of-the-art observational data on orbits and wind properties of individual wind-producing stars. Since wind velocities were revised upwards and non-zero eccentricities were considered, our new simulations show fewer clumps of cold gas and no conspicuous disc-like structure. The accretion rate is dominated by a few close 'slow-wind stars' ( v _w≤ 750 km s⁻¹), and is consistent with the Bondi estimate, but variable on time-scales of tens to hundreds of years. This variability is due to the stochastic infall of cold clumps of gas, as in earlier simulations, and to the eccentric orbits of stars. The present models fail to explain the high luminosity of Sgr A* a few hundred years ago implied by Integral observations, but we argue that the accretion of a cold clump with a small impact parameter could have caused it. Finally, we show the possibility of constraining the total mass-loss rate of the 'slow-wind stars' using near infrared observations of gas in the central few arcseconds.     

Kinematics of hypervelocity stars in the triaxial halo of the Milky Way

Hypervelocity stars (HVSs) ejected by the massive black hole at the Galactic Centre have unique kinematic properties compared to other halo stars. Their trajectories will deviate from being exactly radial because of the asymmetry of the Milky Way potential produced by the flattened disc and the triaxial dark matter halo, causing a change of angular momentum that can be much larger than the initial small value at injection. We study the kinematics of HVSs and propose an estimator of dark halo triaxiality that is determined only by instantaneous position and velocity vectors of HVSs at large Galactocentric distances ( r ≳ 50 kpc). We show that, in the case of a substantially triaxial halo, the distribution of deflection angles (the angle between the stellar position and velocity vector) for HVSs on bound orbits is spread uniformly over the range 10°–180°. Future astrometric and deep wide-field surveys should measure the positions and velocities of a significant number of HVSs, and provide useful constraints on the shape of the Galactic dark matter halo.     

The influence of multiple stars on the high-mass stellar initial mass function and age dating of young massive star clusters

The study of young stellar populations has revealed that most stars are in binary or higher order multiple systems. In this study, the influence on the stellar initial mass function (IMF) of large quantities of unresolved multiple massive stars is investigated by taking into account the stellar evolution and photometrically determined system masses. The models, where initial masses are derived from the luminosity and colour of unresolved multiple systems, show that even under extreme circumstances (100 per cent binaries or higher order multiples), the difference between the power-law index of the mass function (MF) of all stars and the observed MF is small (≲0.1). Thus, if the observed IMF has the Salpeter index  α= 2.35  , then the true stellar IMF has an index not flatter than  α= 2.25  . Additionally, unresolved multiple systems may hide between 15 and 60 per cent of the underlying true mass of a star cluster. While already a known result, it is important to point out that the presence of a large number of unresolved binaries amongst pre-main-sequence stars induces a significant spread in the measured ages of these stars even if there is none. Also, lower mass stars in a single-age binary-rich cluster appear older than the massive stars by about 0.6 Myr.     

Revisiting the relations: Galactic thin disc age–velocity dispersion relation

The velocity dispersion of stars in the solar neighbourhood thin disc increases with time after star formation. Nordström et al. performed the most recent observations to constrain the age–velocity dispersion relation. They fitted the age–velocity dispersion relations of each Galactic cardinal direction space velocity component, U (towards the Galactic Centre), V (in the direction of Galactic rotation) and W (towards the North Galactic Pole), with power laws and interpreted these as evidence for continuous heating of the disc in all directions throughout its lifetime. We revisit these relations with their data and use the results of Famaey et al. to show that structure in the local velocity distribution function distorts the in-plane ( U and V ) velocity distributions away from Gaussian so that a dispersion is not an adequate parametrization of their functions. The age–σ _W relation can however be constrained because the sample is well phase-mixed vertically. We do not find any local signature of the stellar warp in the Galactic disc. Vertical disc heating does not saturate at an early stage. Our new result is that a power law is not required by the data: disc heating models that saturate after ∼4.5 Gyr are equally consistent with observations.     

The destructive effects of binary encounters on red giants in the Galactic Centre

We consider the destructive effects of encounters between binaries and red giant stars in the Galactic Centre. Such encounters may explain the observed depletion of luminous red giants within the central 0.2 pc of the galaxy. We consider encounters involving 2- and 8-M⊙ red giants, and thus span the range of stellar masses contributing to the most luminous red giants observed in the Galactic Centre. To explore the phase space of encounters thoroughly, we simulate 18 × 10³ encounters using a modified four-body code in which the red giant core and components of the binary are treated as point masses, and where the envelope configuration is assumed to remain static throughout the encounter. We then rerun a small number of encounters with a smoothed particle hydrodynamics (SPH) code to confirm the reliability of conclusions drawn from the four-body runs. We see two possible pathways to red giant destruction. A large fraction of encounters lead to the formation of common-envelope systems, where two compact objects (drawn from the red giant core and the components of the original binary) form a binary within a common gaseous envelope, whilst the third body is ejected. The destruction of the red giant will then follow when the envelope is ejected as the binary hardens. In a smaller number of encounters, the intruding binary passes through the star and ejects the red giant core from the envelope. The red giant envelope will then disperse on short time-scales. We compute the time-scales for both of these processes to occur in the Galactic Centre for a variety of binary populations.     

Gravitational instability of discs with dissipative coronae around supermassive black holes

We study the dynamical structure of a self-gravitating disc with coronae around a supermassive black hole. Assuming that the magnetorotational instability responsible for generating the turbulent stresses inside the disc is also the source for a magnetically dominated corona, a fraction of the power released when the disc matter accretes is transported to and dissipated in the corona. This has a major effect on the structure of the disc and its gravitational (in)stability according to our analytical and self-consistent solutions. We determine the radius where the disc crosses the inner radius of gravitational instability and forms the first stars. Not only the location of this radius which may extend to very large distances from the central black hole, but also the mass of the first stars highly depends on the input parameters, notably the viscosity coefficient, the mass of the central object and the accretion rate. For accretion discs around quasi-stellar objects (QSOs) and the Galactic Centre, we determine the self-gravitating radius and the mass of the first clumps. Comparing the cases with a corona and without a corona for typical discs around QSOs or the Galactic Centre, when the viscosity coefficient is around 0.3, we show that the self-gravitating radius decreases by a factor of approximately 2, but the mass of the fragments increases with more or less the same factor. The existence of a corona implies a more gravitationally unstable disc according to our results. The effect of a corona on the instability of the disc is more effective when the viscosity coefficient increases.     

Something about the structure of the Galaxy

We analyse a sample of 507 evolved (OH/IR) stars in the region (10°>ℓ>−45°), (| b |<3°). We derive average ages for subsets of this sample, and use those sets as beacons for the evolution of the Galaxy. In the bulge, the oldest OH/IR stars in the plane are 7.5 Gyr (1.3 M_⊙), and in the disc 2.7 Gyr (2.3 M_⊙). The vertical distribution of almost all AGB stars in the disc is found to be nearly exponential, with scaleheight increasing from 100 pc for stars ≲1 Gyr old to 500 pc for stars ≳5 Gyr old. There may be a small, disjunct population of OH/IR stars. The radial distribution of AGB stars is dictated by the metallicity gradient. Unequivocal morphological evidence is presented for the existence of a central bar, but parameters can be constrained only for a given spatial-density model. Using a variety of indicators, we identify the radii of the inner ultraharmonic (2.5 kpc) and corotation resonance (3.5 kpc). We show that the 3-kpc arm is likely to be an inner ring, as observed in other barred galaxies, by identifying a group of evolved stars that is connected to the 3-kpc H  i filament. Also, using several observed features, we argue that an inner-Lindblad resonance exists, at ∼1–1.5 kpc. The compositions of OH/IR populations within 1 kpc of the Galactic Centre give insight into the bar-driven evolution of the inner regions. We suggest that the bar is ∼8 Gyr old, relatively weak (SAB), and may be in a final stage of its existence.     

A search for planets transiting the M-dwarf debris disc host, AU Microscopii

We present high-cadence, high-precision multiband photometry of the young, M1Ve, debris disc star, AU Microscopii. The data were obtained in three continuum filters spanning a wavelength range from 4500 to 6600 Å, plus Hα, over 28 nights in 2005. The light curves show intrinsic stellar variability due to star-spots with an amplitude in the blue band of 0.051 mag and a period of 4.847 d. In addition, three large flares were detected in the data which all occur near the minimum brightness of the star. We remove the intrinsic stellar variability and combine the light curves of all the filters in order to search for transits by possible planetary companions orbiting in the plane of the nearly edge-on debris disc. The combined final light curve has a sampling of 0.35 min and a standard deviation of 6.8 mmag. We performed Monte Carlo simulations by adding fake transits to the observed light curve and find with 95 per cent significance that there are no Jupiter mass planets orbiting in the plane of the debris disc on circular orbits with periods,   P ≤ 5  d. In addition, there are no young Neptune like planets (with radii 2.5 times smaller than the young Jupiter) on circular orbits with periods,   P ≤ 3  d.