We present a full dynamical model of the Fornax dwarf spheroidal galaxy obtained with the spherically symmetric Schwarzschild orbit superposition method applied to the largest kinematic data set presently available. We modelled the total mass content of the dwarf with the mass-to-light ratio $\Upsilon$ varying with radius and found that Fornax is with high probability embedded in a massive and extended dark matter halo. We estimated the total mass contained within 1$\,$kpc to be $M(<1\,$kpc$)=1.25^{+0.06}_{-0.13} \times 10^8\,$M$_{\odot}$. The data are consistent with the constant mass-to-light ratio, i.e. the mass-follows-light model, only at 3$\sigma$ level, but still require a high value of $\Upsilon \approx 11.2\,$M$_{\odot}/$L$_{\odot}$. Our results are in general agreement with previous estimates of the dynamical mass profile of Fornax. As the Schwarzschild method does not require any assumptions on the orbital anisotropy of the stars, we obtained a profile of the anisotropy parameter $\beta$ as an output of our modelling. The derived anisotropy is close to zero in the centre of the galaxy and decreases with radius, but remains consistent with isotropic orbits at all radii at 1$\sigma$ confidence level.

In our previous work we confirmed the reliability of the spherically symmetric Schwarzschild orbit-superposition method to recover the mass and velocity anisotropy profiles of spherical dwarf galaxies. Here we investigate the effect of its application to intrinsically non-spherical objects. For this purpose we use a model of a dwarf spheroidal galaxy formed in a numerical simulation of a major merger of two disky dwarfs. The shape of the stellar component of the merger remnant is axisymmetric and prolate which allows us to identify and measure the bias caused by observing the spheroidal galaxy along different directions, especially the longest and shortest principal axis. The modelling is based on mock data generated from the remnant that are observationally available for dwarfs: projected positions and line-of-sight velocities of the stars. In order to obtain a reliable tool while keeping the number of parameters low we parametrize the total mass distribution as a radius-dependent mass-to-light ratio with just two free parameters we aim to constrain. Our study shows that if the total density profile is known, the true, radially increasing anisotropy profile can be well recovered for the observations along the longest axis whereas the data along the shortest axis lead to the inference an incorrect, isotropic model. On the other hand, if the density profile is derived from the method as well, the anisotropy is always underestimated but the total mass profile is well recovered for the data along the shortest axis whereas for the longest axis the mass content is overestimated.

We present a new study concerning the application of the Schwarzschild orbit superposition method to model spherical galaxies. The method aims to recover the mass and the orbit anisotropy parameter profiles of the objects using measurements of positions and line-of-sight velocities usually available for resolved stellar populations of dwarf galaxies in the Local Group. To test the reliability of the method, we used different sets of mock data extracted from four numerical realizations of dark matter haloes. The models shared the same density profile but differed in anisotropy profiles, covering a wide range of possibilities, from constant to increasing and decreasing with radius. The tests were done in two steps, first assuming that the mass profile of the dwarf is known and employing the method to retrieve the anisotropy only, and then varying also the mass distribution. We used two kinds of data samples: unrealistically large ones based on over 270 000 particles from the numerical realizations and small ones matching the amount of data available for the Fornax dwarf. For the large data samples we recover both the mass and the anisotropy profiles with very high accuracy. For the realistically small ones we also find a reasonably good agreement between the fitted and the input anisotropies, however the total density profiles can be significantly biased as a result of their oversensitivity to the available data. Our results therefore provide convincing evidence in favour of the applicability of the Schwarzschild method to break the mass-anisotropy degeneracy in dwarf galaxies.

We apply the Schwarzschild orbit superposition method to mock data in order to investigate the accuracy of recovering the profile of the orbit anisotropy. The mock data come from four numerical realizations of dark matter haloes with well defined anisotropy profiles. We show that when assuming a correct mass distribution we are able to determine the anisotropy with high precision and clearly distinguish between the models.

Using N-body simulations, we study the properties of a bar induced in a dwarf galaxy as a result of tidal interaction with the Milky Way. The dwarf is initially composed of a disk embedded in a dark matter halo and we follow its evolution on a typical orbit for 10 Gyr. It undergoes an evolution typical of tidally stirred dwarf galaxies: it loses mass, the stellar component transforms from a disk to a spheroid and the rotation of the stars is partially replaced by random motions. A tidally induced bar forms at the first pericentre passage and survives until the end of the evolution. Fourier decomposition of the distribution of stars reveals that only even modes are significant and preserve a hierarchy so that the bar mode is always the strongest. They show a characteristic profile with a maximum, similar to simulated bars forming in isolated galaxies and observed bars in real galaxies. We adopt the maximum of the bar mode as a measure of the bar strength and we estimate the bar length by comparing the density profiles along the bar and perpendicular to it. Both the bar strength and the bar length decrease with time, mainly at pericentres, as a result of tidal torques acting at those times and not to secular evolution. The pattern speed of the bar varies significantly on a time scale of 1 Gyr and seems to be fully controlled by the orientation of the tidal torque from the Milky Way. The bar is never tidally locked but we discover a hint of a 5/2 orbital resonance between the third and fourth pericentre passage. The speed of the bar expressed in terms of the ratio of the bar length to the corotation radius decreases in the long run so that the bar changes from initially rather fast to slow in the later stages. The boxy/peanut shape is clearly visible in the edge-on view between the second and third pericentre and its occurrence is preceded by a short period of buckling instability.

Precise measurements of mass in dark matter dominated dwarf spheroidal galaxies are of great importance for testing the theories of structure formation. We use $N$-body simulations of the tidal evolution of a dwarf galaxy orbiting the Milky Way to generate mock kinematical data sets and use them to test the reliability of a simple mass estimator proposed by Wolf et al. The evolution of the initially disky dwarf galaxy embedded in a dark matter halo was traced for 10 Gyr on a rather tight orbit. After about half of the time a dwarf spheroidal galaxy is formed that retains some remnant rotation and a non-spherical shape. Observing the triaxial galaxy along each of its principal axes we measure its half-light radius and the line-of-sight velocity dispersion and use them to estimate the mass. We find that the mass is significantly overestimated when the dwarf is seen along the longest axis of the stellar component and underestimated when observed along the shortest axis. We provide a formula that quantifies the systematic error in the estimated mass with respect to the true one as a function of the galaxy shape and line of sight.

Using collisionless N-body simulations of dwarf galaxies orbiting the Milky Way (MW) we construct realistic models of dwarf spheroidal (dSph) galaxies of the Local Group. The dwarfs are initially composed of stellar disks embedded in dark matter haloes with different inner density slopes and are placed on an eccentric orbit typical for MW subhaloes. After a few Gyr of evolution the stellar component is triaxial as a result of bar instability induced by tidal forces. Observing the simulated dwarfs along three principal axes of the stellar component we create mock data sets and determine the their half-light radii and line-of-sight velocity dispersions. Using the estimator proposed by Wolf et al. we calculate masses within half-light radii. The masses obtained this way are over(under)estimated by up to a factor of two when the line of sight is along the longest (shortest) axis of the stellar component. We then divide the initial stellar distribution into an inner and outer population and trace their evolution in time. The two populations, although affected by tidal forces, retain different density profiles even after a few Gyr. We measure the half-light radii and velocity dispersions of the stars in the two populations along different lines of sight and use them to estimate the slope of the mass distribution in the dwarfs following the method proposed by Walker & Penarrubia. The inferred slopes are systematically over- or underestimated, depending on the line of sight. In particular, when the dwarf is seen along the longest axis of the stellar component, a significantly shallower density profile is inferred than the real one measured from the simulations. Since most dSphs are non-spherical and their orientation with respect to our line of sight is unknown, the method can be reliably applied only to a large sample of dwarfs when these systematic errors are expected to be diminished.

Using N-body simulations we study the evolution of separate stellar populations in dwarf galaxies in the context of the tidal stirring scenario for the formation of dwarf spheroidal (dSph) galaxies in the Local Group. The dwarf galaxies, initially composed of a stellar disk and a dark matter halo, are placed on seven different orbits around the Milky Way. The stars are divided into two populations, within and outside the half-light radius, and their positions are followed for 10 Gyr. We find that the populations retain different density distributions even over such long timescales. Some of the stars of the outer population migrate to the central part of the dwarf forming an extended core while the stars of the inner population develop a tail in the outer parts. In addition, the outer population is more heavily stripped by tidal forces from the Milky Way and may become subdominant at all radii on tight enough orbits. We conclude that the tidal stirring model is fully compatible with the presence of multiple stellar populations in dSph galaxies.

We present a sample of 209 variable objects - very likely optical counterparts to the X-ray sources detected in the direction of the Galactic center by the Galactic Bulge Survey (GBS) carried out with the Chandra satellite. The variable sources were found in the databases of the OGLE long term survey monitoring regularly the Galactic bulge since 1992. The counterpart candidates were searched based on the X-ray source position in the radius of 3."9. Optical light curves of the candidates comprise a full variety of variability types: spotted stars, pulsating red giants (potentially secondary stars of symbiotic variables), cataclysmic variables, eclipsing binary systems, irregular non-periodic objects including an AGN (GRS 1734-292). Additionally, we find that positions of 19 non-variable stars brighter than 16.5 mag in the OGLE databases are so well aligned with the X-ray positions (<0."75) that these objects are also likely optical counterparts to the GBS X-ray sources. We provide the OGLE astrometric and photometric information for all selected objects and their preliminary classifications. Photometry of the candidates is available from the OGLE Internet archive, http://ogle.astrouw.edu.pl