-
Far-infrared line emission from the outer Galaxy cluster Gy 3-7 with SOFIA/FIFI-LS: Physical conditions and UV fields
Authors:
N. Le,
A. Karska,
M. Figueira,
M. Sewiło,
A. Mirocha,
Ch. Fischer,
M. Kaźmierczak-Barthel,
R. Klein,
M. Gawroński,
M. Koprowski,
K. Kowalczyk,
W. J. Fischer,
K. M. Menten,
F. Wyrowski,
C. König,
L. E. Kristensen
Abstract:
(abridged) Far-infrared (FIR) line emission provides key information about the gas cooling and heating due to shocks and UV radiation associated with the early stages of star formation. Gas cooling via FIR lines might, however, depend on metallicity. We aim to quantify the FIR line emission and determine the spatial distribution of the CO rotational temperature, ultraviolet (UV) radiation field, a…
▽ More
(abridged) Far-infrared (FIR) line emission provides key information about the gas cooling and heating due to shocks and UV radiation associated with the early stages of star formation. Gas cooling via FIR lines might, however, depend on metallicity. We aim to quantify the FIR line emission and determine the spatial distribution of the CO rotational temperature, ultraviolet (UV) radiation field, and H2 number density toward the embedded cluster Gy 3-7 in the CMa-l224 star-forming region, whose metallicity is expected to be intermediate between that of the LMC and the Solar neighborhood. By comparing the total luminosities of CO and [O I] toward Gy 3-7 with values found for low- and high-mass protostars extending over a broad range of metallicities, we also aim to identify the possible effects of metallicity on the FIR line cooling within our Galaxy. We studied SOFIA/FIFI-LS spectra of Gy 3-7 covering several FIR lines. The spatial extent of CO high-J (J>14) emission resembles that of the elongated 160 um continuum emission detected with Herschel. The CO transitions from J=14-13 to J=16-15 are detected throughout the cluster and show a median rotational temperature of 170+/-30 K on Boltzmann diagrams. Comparisons to other protostars observed with Herschel show a good agreement with intermediate-mass sources in the inner Galaxy. Assuming an origin of the [O I] and high-J CO emission in UV-irradiated C-shocks, we obtained pre-shock H2 number densities of 10^4-5 cm-3 and UV radiation field strengths of 0.1-10 Habing fields. Far-IR line observations reveal ongoing star formation in Gy 3-7, dominated by intermediate-mass Class 0/I young stellar objects. The ratio of molecular-to-atomic far-IR line emission shows a decreasing trend with bolometric luminosities of the protostars. However, it does not indicate that the low-metallicity has an impact on the line cooling in Gy 3-7.
△ Less
Submitted 17 April, 2023;
originally announced April 2023.
-
Multiple stellar populations in Schwarzschild modeling and the application to the Fornax dwarf
Authors:
Klaudia Kowalczyk,
Ewa L. Lokas
Abstract:
Dwarf spheroidal (dSph) galaxies are believed to be strongly dark matter dominated and thus are considered perfect objects to study dark matter distribution and test theories of structure formation. They possess resolved, multiple stellar populations that offer new possibilities for modeling. A promising tool for the dynamical modeling of these objects is the Schwarzschild orbit superposition meth…
▽ More
Dwarf spheroidal (dSph) galaxies are believed to be strongly dark matter dominated and thus are considered perfect objects to study dark matter distribution and test theories of structure formation. They possess resolved, multiple stellar populations that offer new possibilities for modeling. A promising tool for the dynamical modeling of these objects is the Schwarzschild orbit superposition method. In this work we extend our previous implementation of the scheme to include more than one population of stars and a more general form of the mass-to-light ratio function. We tested the improved approach on a nearly spherical, gas-free galaxy formed in the cosmological context from the Illustris simulation. We modeled the binned velocity moments for stars split into two populations by metallicity and demonstrate that in spite of larger sampling errors the increased number of constraints leads to significantly tighter confidence regions on the recovered density and velocity anisotropy profiles. We then applied the method to the Fornax dSph galaxy with stars similarly divided into two populations. In comparison with our earlier work, we find the anisotropy parameter to be slightly increasing, rather than decreasing, with radius and more strongly constrained. We are also able to infer anisotropy for each stellar population separately and find them to be significantly different.
△ Less
Submitted 1 January, 2022;
originally announced January 2022.
-
Schwarzschild dynamical model of the Fornax dwarf spheroidal galaxy
Authors:
Klaudia Kowalczyk,
Andrés del Pino,
Ewa L. Lokas,
Monica Valluri
Abstract:
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 $Υ$ varying with radius and found that Fornax is with high probability embedded in a massive and extended d…
▽ More
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 $Υ$ 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$σ$ level, but still require a high value of $Υ\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 $β$ 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$σ$ confidence level.
△ Less
Submitted 13 November, 2018; v1 submitted 20 July, 2018;
originally announced July 2018.
-
The effect of non-sphericity on mass and anisotropy measurements in dSph galaxies with Schwarzschild method
Authors:
Klaudia Kowalczyk,
Ewa L. Lokas,
Monica Valluri
Abstract:
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 mer…
▽ More
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.
△ Less
Submitted 16 February, 2018; v1 submitted 30 August, 2017;
originally announced August 2017.
-
Recovering the mass profile and orbit anisotropy of mock dwarf galaxies with Schwarzschild modelling
Authors:
Klaudia Kowalczyk,
Ewa L. Lokas,
Monica Valluri
Abstract:
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 t…
▽ More
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.
△ Less
Submitted 19 June, 2017; v1 submitted 20 February, 2017;
originally announced February 2017.
-
Orbit anisotropy of dark matter haloes with Schwarzschild modelling
Authors:
Klaudia Kowalczyk,
Ewa L. Lokas,
Monica Valluri
Abstract:
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 disti…
▽ More
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.
△ Less
Submitted 22 December, 2015;
originally announced December 2015.
-
Estimating masses of dwarf spheroidal galaxies
Authors:
Klaudia Kowalczyk,
Ewa L. Lokas
Abstract:
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 initial…
▽ More
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.
△ Less
Submitted 22 August, 2014;
originally announced August 2014.
-
Adventures of a tidally induced bar
Authors:
E. L. Lokas,
E. Athanassoula,
V. P. Debattista,
M. Valluri,
A. del Pino,
M. Semczuk,
G. Gajda,
K. Kowalczyk
Abstract:
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 t…
▽ More
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.
△ Less
Submitted 23 September, 2014; v1 submitted 4 April, 2014;
originally announced April 2014.
-
Can we measure the slopes of density profiles in dwarf spheroidal galaxies?
Authors:
Klaudia Kowalczyk,
Ewa L. Lokas,
Stelios Kazantzidis,
Lucio Mayer
Abstract:
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 compone…
▽ More
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.
△ Less
Submitted 1 March, 2013; v1 submitted 14 December, 2012;
originally announced December 2012.
-
Stellar populations in tidally stirred dwarf galaxies
Authors:
Ewa L. Lokas,
Klaudia Kowalczyk,
Stelios Kazantzidis
Abstract:
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…
▽ More
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.
△ Less
Submitted 4 December, 2012;
originally announced December 2012.
-
The Optical Gravitational Lensing Experiment. Optical Counterparts to the X-ray Sources in the Galactic Bulge
Authors:
A. Udalski,
K. Kowalczyk,
I. Soszynski,
R. Poleski,
M. K. Szymanski,
M. Kubiak,
G. Pietrzynski,
S. Kozlowski,
P. Pietrukowicz,
K. Ulaczyk,
J. Skowron,
L. Wyrzykowski
Abstract:
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…
▽ More
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
△ Less
Submitted 1 August, 2012;
originally announced August 2012.