Wolf-Rayet stars are advanced evolutionary stages of massive stars. Despite their large mass-loss rates and high wind velocities, none of them display a bow shock, although a fraction of them are classified as runaway. Our 2.5-D numerical simulations of circumstellar matter around a 60Mo runaway star show that the fast Wolf-Rayet stellar wind is released into a wind-blown cavity filled with various shocks and discontinuities generated throughout the precedent evolutionary phases. The resulting fast-wind slow-wind interaction leads to the formation of spherical shells of swept-up dusty material similar to those observed in near infrared 24 micron with Spitzer, and which appear to be co-moving with the runaway massive stars, regardless of their proper motion and/or the properties of the local ambient medium. We interpret bright infrared rings around runaway Wolf-Rayet stars in the Galactic plane, like WR138a, as indication of their very high initial masses and a complex evolutionary history. Stellar-wind bow shocks become faint as stars run in diluted media, therefore, our results explain the absence of detected bow shocks around Galactic Wolf-Rayet stars such as the high-latitude, very fast-moving objects WR71, WR124 and WR148. Our results show that the absence of a bow shock is consistent with a runaway nature of some Wolf-Rayet stars. This questions the in-situ star formation scenario of high-latitude Wolf-Rayet stars in favor of dynamical ejection from birth sites in the Galactic plane.

A very small fraction of (runaway) massive stars have masses exceeding $60$-$70\, \rm M_{\odot}$ and are predicted to evolve as Luminous-Blue-Variable and Wolf-Rayet stars before ending their lives as core-collapse supernovae. Our 2D axisymmetric hydrodynamical simulations explore how a fast wind ($2000\, \rm km\, \rm s^{-1}$) and high mass-loss rate ($10^{-5}\, \rm M_{\odot}\, \rm yr^{-1}$) can impact the morphology of the circumstellar medium. It is shaped as 100 pc-scale wind nebula which can be pierced by the driving star when it supersonically moves with velocity $20$-$40\, \rm km\, \rm s^{-1}$ through the interstellar medium (ISM) in the Galactic plane. The motion of such runaway stars displaces the position of the supernova explosion out of their bow shock nebula, imposing asymmetries to the eventual shock wave expansion and engendering Cygnus-loop-like supernova remnants. We conclude that the size (up to more than $200\, \rm pc$) of the filamentary wind cavity in which the chemically enriched supernova ejecta expand, mixing efficiently the wind and ISM materials by at least $10\%$ in number density, can be used as a tracer of the runaway nature of the very massive progenitors of such $0.1\, \rm Myr$ old remnants. Our results motivate further observational campaigns devoted to the bow shock of the very massive stars BD+43 3654 and to the close surroundings of the synchrotron-emitting Wolf-Rayet shell G2.4+1.4.

How the Milky Way has accumulated its mass over the Hubble time, whether significant amounts of gas and stars were accreted from satellite galaxies, or whether the Milky Way has experienced an initial gas assembly and then evolved more-or-less in isolation is one of the burning questions in modern astronomy, because it has consequences for our understanding of galaxy formation in the cosmological context. Here we present the evolutionary model of a Milky Way-type satellite system zoomed into a cosmological large-scale simulation. Embedded into Dark Matter halos and allowing for baryonic processes these chemo-dynamical simulations aim at studying the gas and stellar loss from the satellites to feed the Milky Way halo and the stellar chemical abundances in the halo and the satellite galaxies.

We examine the fraction of massive asymptotic giant branch (AGB) stars remaining bound in their parent star clusters and the effect of irradiation of these stars by intracluster ultraviolet (UV) field. We employ a set of N-body models of dynamical evolution of star clusters rotating in a galactic potential at the solar galactocentric radius. The cluster models are combined with stellar evolution formulae, a library of stellar spectra, and simple models for SiO photodissociation in circumstellar environment (CSE). The initial stellar masses of clusters are varied from $50\rm M_\odot$ to $10^{5}\rm M_\odot$. Results derived for individual clusters are combined using a mass distribution function for young star clusters. We find that about 30% of massive AGB stars initially born in clusters become members of the field population, while the rest evolves in star clusters. They are irradiated by strong intracluster UV radiation resulting in the decrease of the photodissociation radius of SiO molecules, in many stars down to the dust formation zone. In absence of dust shielding, the UV photons penetrate in the CSE deeper than $10R_*$ in 64% and deeper than $2 R_*$ in 42% of all massive AGB stars. If this suppresses following dust formation, the current injection rate of silicate dust from AGB stars in the local Galaxy decreases from $2.2 \times 10^{-4}\rm M_\odot\,kpc^{-2}\,Gyr^{-1}$ to $1.8 \times 10^{-4}\rm M_\odot\,kpc^{-2}\,Gyr^{-1}$ at most. A lower revised value of 40% for the expected fraction of presolar silicate grains from massive AGB stars is still high to explain the non-detection of these grains in meteorites.

Within the Cold Dark Matter scenario the hierarchical merging paradigm is the natural result to form massive galactic halos by the minor mergers of sub-halos and, by this, inherently their stellar halo. Although this must be also invoked for the Milky Way, the context of chemical and kinematic coherence of halo stars and dwarf spheroidal galaxies is yet unsolved a focus of present-day research. To examine this issue we model the chemo-dynamical evolution of the system of satellites selected from the cosmological Via Lactea II simulations to be similar for the Milky Way environment but at an early epoch.

Context: A recent new approach to apply a simple dynamical mass estimate of tidally limited star clusters is based on the identification of the tidal radius in a King profile with the dynamical Jacobi radius. The application to an unbiased open cluster catalogue yields significantly higher cluster masses compared to the classical methods. Aims: We quantify the bias in the mass determination as function of projection direction and cluster age by analysing a simulated star cluster. Methods: We use direct $N$-body simulations of a star cluster including stellar evolution in an analytic Milky Way potential and apply a best fit to the projected number density of cluster stars. Results: We obtain significantly overestimated star cluster masses which depend strongly on the viewing direction. The overestimation is typically in the range of 10-50 percent and reaches a factor of 3.5 for young clusters. Mass segregation reduces the derived limiting radii systematically.

(abridged) In this paper we derive observed and modelled shape parameters (apparent ellipticity and orientation of the ellipse) of 650 Galactic open clusters identified in the ASCC-2.5 catalogue. We provide the observed shape parameters of Galactic open clusters, computed with the help of a multi-component analysis. For the vast majority of clusters these parameters are determined for the first time. High resolution ("star by star") N-body simulations are carried out with the specially developed $\phi$GRAPE code providing models of clusters of different initial masses, Galactocentric distances and rotation velocities. The comparison of models and observations of about 150 clusters reveals ellipticities of observed clusters which are too low (0.2 vs. 0.3), and offers the basis to find the main reason for this discrepancy. The models predict that after $\approx 50$ Myr clusters reach an oblate shape with an axes ratio of $1.65:1.35:1$, and with the major axis tilted by an angle of $q_{XY} \approx 30^\circ$ with respect to the Galactocentric radius due to differential rotation of the Galaxy. Unbiased estimates of cluster shape parameters require reliable membership determination in large cluster areas up to 2-3 tidal radii where the density of cluster stars is considerably lower than the background. Although dynamically bound stars outside the tidal radius contribute insignificantly to the cluster mass, their distribution is essential for a correct determination of cluster shape parameters. In contrast, a restricted mass range of cluster stars does not play such a dramatic role, though deep surveys allow to identify more cluster members and, therefore, to increase the accuracy of the observed shape parameters.

Tidal tails of star clusters are not homogeneous but show well defined clumps in observations as well as in numerical simulations. Recently an epicyclic theory for the formation of these clumps was presented. A quantitative analysis was still missing. We present a quantitative derivation of the angular momentum and energy distribution of escaping stars from a star cluster in the tidal field of the Milky Way and derive the connection to the position and width of the clumps. For the numerical realization we use star-by-star $N$-body simulations. We find a very good agreement of theory and models. We show that the radial offset of the tidal arms scales with the tidal radius, which is a function of cluster mass and the rotation curve at the cluster orbit. The mean radial offset is 2.77 times the tidal radius in the outer disc. Near the Galactic centre the circumstances are more complicated, but to lowest order the theory still applies. We have also measured the Jacobi energy distribution of bound stars and showed that there is a large fraction of stars (about 35%) above the critical Jacobi energy at all times, which can potentially leave the cluster. This is a hint that the mass loss is dominated by a self-regulating process of increasing Jacobi energy due to the weakening of the potential well of the star cluster, which is induced by the mass loss itself.

We have monitored two bright galaxies, M85 and NGC7331, on timescales as short as 0.01s. In the optical, we discovered an unusual burst coincident with the galaxy M85. We registered a sudden onset with a characteristic time of less than 10ms with subsequent quasi-exponential decay within approximately 1s and an amplitude of 2.5mag in the V-band. In the course of high-speed monitoring with two Crimean telescopes operated synchronously, in both independent instruments we have registered one coincident event occurring in NGC7331 with a duration of ~0.6s. The amplitudes range from ~3mag to ~0.3mag in the U-band and I-band, respectively. Merging of an intermediate-mass black hole with a small black hole or normal star seems to be the most plausible mechanism responsible for short bursts. Our observations support the hypothesis concerning the existence of intermediate-mass black holes in the centers of galaxies and in dense globular clusters.

In this paper we study the process of the subsequent (runaway) fragmentation of the rotating isothermal Giant Molecular Cloud (GMC) complex. Our own developed Smoothed Particle Hydrodynamics (SPH) gas-dynamical model successfully reproduce the observed Cloud Mass-distribution Function (CMF) in our Galaxy (even the differences between the inner and outer parts of our Galaxy). The steady state CMF is established during the collapse within a free-fall timescale of the GMC. We show that one of the key parameters, which defines the observed slope of the present day CMF, is the initial ratio of the rotational (turbulent) and gravitational energy inside the fragmented GMC.