Showing 1–28 of 28 results for author: Smirnov, A A

Isochrone Fitting of Galactic Globular Clusters -- VI. High-latitude Clusters NGC5024 (M53), NGC5053, NGC5272 (M3), NGC5466, and NGC7099 (M30)

Authors: G. A. Gontcharov, S. S. Savchenko, A. A. Marchuk, C. J. Bonatto, O. S. Ryutina, M. Yu. Khovritchev, V. B. Il'in, A. V. Mosenkov, D. M. Poliakov, A. A. Smirnov

Abstract: We fit various colour-magnitude diagrams (CMDs) of the high-latitude Galactic globular clusters NGC\,5024 (M53), NGC\,5053, NGC\,5272 (M3), NGC\,5466, and NGC\,7099 (M30) by isochrones from the Dartmouth Stellar Evolution Database and Bag of Stellar Tracks and Isochrones for $α$-enrichment [$α$/Fe]$=+0.4$. For the CMDs, we use data sets from {\it Hubble Space Telescope}, {\it Gaia}, and other sour… ▽ More We fit various colour-magnitude diagrams (CMDs) of the high-latitude Galactic globular clusters NGC\,5024 (M53), NGC\,5053, NGC\,5272 (M3), NGC\,5466, and NGC\,7099 (M30) by isochrones from the Dartmouth Stellar Evolution Database and Bag of Stellar Tracks and Isochrones for $α$-enrichment [$α$/Fe]$=+0.4$. For the CMDs, we use data sets from {\it Hubble Space Telescope}, {\it Gaia}, and other sources utilizing, at least, 25 photometric filters for each cluster. We obtain the following characteristics with their statistic uncertainties for NGC\,5024, NGC\,5053, NGC\,5272, NGC\,5466, and NGC\,7099, respectively: metallicities [Fe/H]$=-1.93\pm0.02$, $-2.08\pm0.03$, $-1.60\pm0.02$, $-1.95\pm0.02$, and $-2.07\pm0.04$ dex with their systematic uncertainty 0.1 dex; ages $13.00\pm0.11$, $12.70\pm0.11$, $11.63\pm0.07$, $12.15\pm0.11$, and $12.80\pm0.17$ Gyr with their systematic uncertainty 0.8 Gyr; distances (systematic uncertainty added) $18.22\pm0.06\pm0.60$, $16.99\pm0.06\pm0.56$, $10.08\pm0.04\pm0.33$, $15.59\pm0.03\pm0.51$, and $8.29\pm0.03\pm0.27$ kpc; reddenings $E(B-V)=0.023\pm0.004$, $0.017\pm0.004$, $0.023\pm0.004$, $0.023\pm0.003$, and $0.045\pm0.002$ mag with their systematic uncertainty 0.01 mag; extinctions $A_\mathrm{V}=0.08\pm0.01$, $0.06\pm0.01$, $0.08\pm0.01$, $0.08\pm0.01$, and $0.16\pm0.01$ mag with their systematic uncertainty 0.03 mag, which suggest the total Galactic extinction $A_\mathrm{V}=0.08$ across the whole Galactic dust to extragalactic objects at the North Galactic pole. The horizontal branch morphology difference of these clusters is explained by their different metallicity, age, mass-loss efficiency, and loss of low-mass members in the evolution of the core-collapse cluster NGC\,7099 and loose clusters NGC\,5053 and NGC\,5466. △ Less

Submitted 23 April, 2024; originally announced April 2024.

Comments: 9 figures, published in Research in Astronomy and Astrophysics

Positive Lynden-Bell derivative as a ticket to the bar trap?

Authors: Viktor D. Zozulia, Anton A. Smirnov, Natalia Ya. Sotnikova

Abstract: We have translated the results of $N$-body simulations of one barred model into the language of action variables and frequencies. Using this language, we analysed the behaviour of all orbits in the model on a large time scale at the stage of a mature bar. We show that the orbits join the bar while preserving their adiabatic invariant, which takes into account the 3D structure of the orbits. This a… ▽ More We have translated the results of $N$-body simulations of one barred model into the language of action variables and frequencies. Using this language, we analysed the behaviour of all orbits in the model on a large time scale at the stage of a mature bar. We show that the orbits join the bar while preserving their adiabatic invariant, which takes into account the 3D structure of the orbits. This allows us to apply the concept of the Lynden-Bell derivative for each of these orbits and trace how the sign of the derivative changes, i.e. how asynchronous changes in angular momentum $L_z$ and orbital precession rate $Ω_\mathrm{pr}$ (normal orbital mode) change to synchronous (abnormal mode). The transition to the abnormal mode occurs when $Ω_\mathrm{pr}$ reaches the angular velocity of the pattern $Ω_\mathrm{p}$, after which the orbit becomes stuck in the bar trap. All this happens against the background of secular changes in actions ($L_z$ decreases, $J_\mathrm{R}$ and $J_z$ increase). At the same time, corotation particles near two stable Lagrange points are also subject to secular changes in their actions. They increase $L_z$ and drift to the periphery, shifting corotation outwards. We also show that a change in the orbital mode from normal to abnormal and the trapping of orbits in a bar is possible only when the bar speed decreases with time, regardless of what is causing the bar to slow down. Our findings clarify and expand the picture of bar formation and evolution in numerical models. △ Less

Submitted 13 March, 2024; originally announced March 2024.

Comments: Accepted in MNRAS

Galaxies decomposition with spiral arms -- II: A multiwavelength case study of M 51

Authors: Alexander A. Marchuk, Ilia V. Chugunov, George A. Gontcharov, Aleksandr V. Mosenkov, Vladimir B. Ilin, Sergey S. Savchenko, Anton A. Smirnov, Denis M. Poliakov, Jonah Seguine, Maxim I. Chazov

Abstract: Spiral structure can contribute significantly to a galaxy's luminosity. However, only rarely are proper photometric models of spiral arms used in decompositions. As we show in the previous work, including the spirals as a separate component in a photometric model of a galaxy would both allow to obtain their structural parameters, and reduce the systematic errors in estimating the parameters of oth… ▽ More Spiral structure can contribute significantly to a galaxy's luminosity. However, only rarely are proper photometric models of spiral arms used in decompositions. As we show in the previous work, including the spirals as a separate component in a photometric model of a galaxy would both allow to obtain their structural parameters, and reduce the systematic errors in estimating the parameters of other components. Doing so in different wavebands, one can explore how their properties vary with the wavelength. In this paper, second in this series, we perform decomposition of M 51 in 17 bands, from the far UV to far IR, using imaging from the DustPedia project. We use the same 2D photometric model of spiral structure where each arm is modelled independently. The complex and asymmetric spiral structure in M 51 is reproduced relatively well with our model. We analyze the differences between models with and without spiral arms, and investigate how the fit parameters change with wavelength. In particular, we find that the spiral arms demonstrate the largest width in the optical, whereas their contribution to the galaxy luminosity is most significant in the UV. The disk central intensity drops by a factor of 1.25--3 and its exponential scale changes by 5--10\% when spiral arms are included, depending on wavelength. Taking into account the full light distribution across the arms, we do not observe the signs of a long-lived density wave in the spiral pattern of M 51 as a whole. △ Less

Submitted 13 February, 2024; originally announced February 2024.

Comments: Accepted for publication in MNRAS

Isochrone fitting of Galactic globular clusters -- V. NGC6397 and NGC6809 (M55)

Authors: George A. Gontcharov, Charles J. Bonatto, Olga S. Ryutina, Sergey S. Savchenko, Aleksandr V. Mosenkov, Vladimir B. Il'in, Maxim Yu. Khovritchev, Alexander A. Marchuk, Denis M. Poliakov, Anton A. Smirnov, Jonah Seguine

Abstract: We fit various colour-magnitude diagrams (CMDs) of the Galactic globular clusters NGC\,6397 and NGC\,6809 (M55) by isochrones from the Dartmouth Stellar Evolution Database (DSED) and Bag of Stellar Tracks and Isochrones (BaSTI) for $α$-enhanced [$α$/Fe]$=+0.4$. For the CMDs, we use data sets from {\it HST}, {\it Gaia}, VISTA, and other sources utilizing 32 and 23 photometric filters for NGC\,6397… ▽ More We fit various colour-magnitude diagrams (CMDs) of the Galactic globular clusters NGC\,6397 and NGC\,6809 (M55) by isochrones from the Dartmouth Stellar Evolution Database (DSED) and Bag of Stellar Tracks and Isochrones (BaSTI) for $α$-enhanced [$α$/Fe]$=+0.4$. For the CMDs, we use data sets from {\it HST}, {\it Gaia}, VISTA, and other sources utilizing 32 and 23 photometric filters for NGC\,6397 and NGC\,6809, respectively, from the ultraviolet to mid-infrared. We obtain the following characteristics for NGC\,6397 and NGC\,6809, respectively: metallicities [Fe/H]$=-1.84\pm0.02\pm0.1$ and $-1.78\pm0.02\pm0.1$ (statistic and systematic uncertainties); distances $2.45\pm0.02\pm0.06$ and $5.24\pm0.02\pm0.18$ kpc; ages $12.9\pm0.1\pm0.8$ and $13.0\pm0.1\pm0.8$ Gyr; reddenings $E(B-V)=0.178\pm0.006\pm0.01$ and $0.118\pm0.004\pm0.01$ mag; extinctions $A_\mathrm{V}=0.59\pm0.01\pm0.02$ and $0.37\pm0.01\pm0.04$ mag; extinction-to-reddening ratio $R_\mathrm{V}=3.32^{+0.32}_{-0.28}$ and $3.16^{+0.66}_{-0.56}$. Our estimates agree with most estimates from the literature. BaSTI gives systematically higher [Fe/H] and lower reddenings than DSED. Despite nearly the same metallicity, age, and helium enrichment, these clusters show a considerable horizontal branch (HB) morphology difference, which must therefore be described by another parameter. This parameter must predominantly explain why the least massive HB stars (0.58-0.63 solar masses) are only found within NGC 6809. Probably they have been lost by the core-collapse cluster NGC\,6397 during its dynamical evolution and mass segregation. In contrast, NGC\,6809 has a very low central concentration and, hence, did not undergo this process. △ Less

Submitted 9 February, 2024; originally announced February 2024.

Comments: 19 figures. arXiv admin note: text overlap with arXiv:2211.12684

Journal ref: Monthly Notices of the Royal Astronomical Society, 2023, Volume 526, Issue 4, pp. 5628-5647

New Interstellar Extinction Maps Based on Gaia and Other Sky Surveys

Authors: G. A. Gontcharov, A. A. Marchuk, M. Yu. Khovrichev, A. V. Mosenkov, S. S. Savchenko, V. B. Il'in, D. M. Polyakov, A. A. Smirnov

Abstract: We present new three-dimensional (3D) interstellar extinction maps in the $V$ and Gaia $G$ filters within 2 kpc of the Sun, a 3D differential extinction (dust spatial distribution density) map along the lines of sight in the same space, a 3D map of variations in the ratio of the extinctions in the $V$ and Gaia $G$ filters within 800 pc of the Sun, and a 2D map of total Galactic extinction through… ▽ More We present new three-dimensional (3D) interstellar extinction maps in the $V$ and Gaia $G$ filters within 2 kpc of the Sun, a 3D differential extinction (dust spatial distribution density) map along the lines of sight in the same space, a 3D map of variations in the ratio of the extinctions in the $V$ and Gaia $G$ filters within 800 pc of the Sun, and a 2D map of total Galactic extinction through the entire dust half-layer from the Sun to extragalactic space for Galactic latitudes $|b|>13^{\circ}$. The 3D maps have a transverse resolution from 3.6 to 11.6 pc and a radial resolution of 50 pc. The 2D map has an angular resolution of 6.1 arcmin. We have produced these maps based on the Gaia DR3 parallaxes and Gaia, Pan-STARRS1, SkyMapper, 2MASS, and WISE photometry for nearly 100 million stars. We have paid special attention to the space within 200 pc of the Sun and high Galactic latitudes as regions where the extinction estimates have had a large relative uncertainty so far. Our maps estimate the extinction within the Galactic dust layer from the Sun to an extended object or through the entire dust half-layer from the Sun to extragalactic space with a precision $σ(A_\mathrm{V})=0.06$ mag. This gives a high relative precision of extinction estimates even at high Galactic latitudes, where, according to our estimates, the median total Galactic extinction through the entire dust half-layer from the Sun to extragalactic objects is $A_\mathrm{V}=0.12\pm0.06$ mag. We have shown that the presented maps are among the best ones in data amount, space size, resolution, precision, and other properties. △ Less

Submitted 9 February, 2024; originally announced February 2024.

Comments: 15 pages, 4 tables

Journal ref: Astronomy Letters, 2023, Vol. 49, No. 11, pp. 666-689

The Influence of the Bar on the Dynamics of Globular Clusters in the Central Region of the Milky Way. Frequency Analysis of Orbits According to Gaia EDR3 Data

Authors: A. T. Bajkova, A. A. Smirnov, V. V. Bobylev

Abstract: This work is devoted to studying the influence of the bar on the orbital dynamics of globular clusters. The orbits of 45 globular clusters in the central galactic region with a radius of 3.5 kpc were analyzed using spectral dynamics methods in order to identify objects captured by the bar. To form the 6D phase space required for orbit integration, the most accurate astrometric data to date from th… ▽ More This work is devoted to studying the influence of the bar on the orbital dynamics of globular clusters. The orbits of 45 globular clusters in the central galactic region with a radius of 3.5 kpc were analyzed using spectral dynamics methods in order to identify objects captured by the bar. To form the 6D phase space required for orbit integration, the most accurate astrometric data to date from the Gaia satellite (EDR3), as well as new refined average distances to globular clusters, were used. Since the parameters of the Milky Way bar are known with very great uncertainty, the orbits were constructed and their frequency analysis was carried out with varying the mass, length and angular velocity of rotation of the bar in a wide range of values with a fairly small step. The integration of orbits was carried out at 2.5 billion years ago. As a result, bar-supporting globular clusters were identified for each set of bar parameters. For the first time, an analytical expression has been obtained for the dependence of the dominant frequency $f_X$ on the angular velocity of rotation of the bar. In addition, the probabilities of capturing globular clusters by the bar were determined when the bar parameters were varied in certain ranges of values according to a random distribution law. A list of 14 globular clusters with the most significant capture probabilities is given, with five GCs - NGC6266, NGC6569, Terzan 5, NGC6522, NGC6540 - showing the probability capture by bar $\geq 0.2$. A conclusion is made about the regularity of the orbits of globular clusters based on the calculation of approximations of the maximum characteristic Lyapunov exponents. △ Less

Submitted 24 November, 2023; originally announced November 2023.

Comments: 17 pages, 13 figures, 4 tables, published in the Astrophysical Bulletin. Russian version: 2023, Astrofizicheskii Byulleten, 2023, vol 78, Issue 4, pp. 525-540

Galaxies decomposition with spiral arms -- II: 29 galaxies from S$^4$G

Authors: Ilia V. Chugunov, Aleksandr V. Mosenkov, Alexander A. Marchuk, Sergey S. Savchenko, Ekaterina V. Shishkina, Maxim I. Chazov, Aleksandra E. Nazarova, Maria N. Skryabina, Polina I. Smirnova, Anton A. Smirnov

Abstract: Spiral structure can occupy a significant part of the galaxy, but properly accounting for it in photometric decomposition is rarely done. This may lead to significant errors in the parameters determined. To estimate how exactly neglecting the presence of spiral arms affects the estimation of galaxy decomposition parameters, we perform fitting of 29 galaxies considering spiral arms as a separate co… ▽ More Spiral structure can occupy a significant part of the galaxy, but properly accounting for it in photometric decomposition is rarely done. This may lead to significant errors in the parameters determined. To estimate how exactly neglecting the presence of spiral arms affects the estimation of galaxy decomposition parameters, we perform fitting of 29 galaxies considering spiral arms as a separate component. In this study, we utilize 3.6$μ$m-band images from the S$^4$G survey and use a new 2D photometric model where each spiral arm is modeled independently. In our model, the light distribution both along and across the arm can be varied significantly, as well as its overall shape. We analyze the differences between models with and without spiral arms, and show that neglecting spiral arms in decomposition causes errors in estimating the parameters of the disk, the bulge, and the bar. We retrieve different parameters of the spiral arms themselves, including their pitch angles, widths, and spiral-to-total luminosity ratio, and examine various relations between them and other galaxy parameters. In particular, we find that the spiral-to-total ratio is higher for galaxies with more luminous discs and with higher bulge-to-total ratios. We report that the pitch angle of spiral arms decreases with increasing bulge or bar fraction. We measure the width of the spiral arms to be 53\% of the disc scale length, on average. We examine the contribution of the spiral arms to the azimuthally-averaged brightness profile and find that spiral arms produce a ``bump'' on this profile with a typical height of 0.3--0.7 mag. △ Less

Submitted 3 November, 2023; originally announced November 2023.

Comments: submitted to MNRAS

Globular clusters and bar: captured or not captured?

Authors: Anton A. Smirnov, Anisa T. Bajkova, Vadim V. Bobylev

Abstract: Studies of the dynamics of globular clusters assume different values of bar parameters (mass, velocity, size) and analyse the results of orbit classifications over the range of the chosen values. It is also a usual thing that a spherical bulge component is converted into the bar to obtain a non-axisymmetric potential from an axisymmetric one. The choice of bar parameters and the way the bar is con… ▽ More Studies of the dynamics of globular clusters assume different values of bar parameters (mass, velocity, size) and analyse the results of orbit classifications over the range of the chosen values. It is also a usual thing that a spherical bulge component is converted into the bar to obtain a non-axisymmetric potential from an axisymmetric one. The choice of bar parameters and the way the bar is converted from the bulge introduce systematics into the orbit classifications that we explore in the present study. We integrate orbits of 30 bulge globular clusters residing in the inner area of the Galaxy ($R \lesssim 5$ kpc) backwards in time for three different potentials, two of which are obtained by fitting the rotation curve, and one is taken from the surrogate $N$-body model representing our Galaxy. We analyse each orbit in terms of dominant frequencies obtained from its coordinate spectra. We find that the bar pattern speed is a key factor in orbital classification. With an increase of it, frequencies deviate more and more from the ``bar'' frequency ratio 2:1. The bar-to-bulge mass ratio (assuming the total mass of the bar plus the bulge is fixed) and size of the bar play a smaller role. We also find that, in the $N$-body potential, the fraction of orbits that follow the bar is higher than in those obtained from fitting the rotation curve. △ Less

Submitted 4 January, 2024; v1 submitted 27 October, 2023; originally announced October 2023.

Comments: Updated: now accepted in MNRAS

B/PS bulges and barlenses from a kinematic viewpoint. I

Authors: Daria Zakharova, Iliya S. Tikhonenko, Natalia Ya. Sotnikova, Anton A. Smirnov

Abstract: A significant part of barred disc galaxies exhibits boxy/peanut-shaped structures (B/PS bulges) at high inclinations. Another structure also associated with the bar is a barlens, often observed in galaxies in a position close to face-on. At this viewing angle, special kinematic tests are required to detect a 3D extension of the bars in the vertical direction (B/PS bulges). We use four pure $N$-bod… ▽ More A significant part of barred disc galaxies exhibits boxy/peanut-shaped structures (B/PS bulges) at high inclinations. Another structure also associated with the bar is a barlens, often observed in galaxies in a position close to face-on. At this viewing angle, special kinematic tests are required to detect a 3D extension of the bars in the vertical direction (B/PS bulges). We use four pure $N$-body models of galaxies with B/PS bulges, which have different bar morphology from bars with barlenses to the so-called face-on peanut bars. We analyse the kinematics of our models to establish how the structural features of B/PS bulges manifest themselves in the kinematics for galaxies at intermediate inclinations and whether these features are related to the barlenses. We apply the dissection of the bar into different orbital groups to determine which of them are responsible for the features of the LOSVD (line-of-sight velocity distribution), i.e., for the deep minima of the $h_4$ parameter along the major axis of the bar. As a result, we claim that for our models at the face-on position, the kinematic signatures of a `peanut' indeed track the vertical density distribution features. We conclude that orbits responsible for such kinematic signatures differ from model to model. We pay special attention to the barlens model. We show that orbits assembled into barlens are not responsible for the kinematic signatures of B/PS bulges. The results presented in this work are applicable to the interpretation of IFU observations of real galaxies. △ Less

Submitted 5 September, 2023; originally announced September 2023.

Comments: 18 pages

Report number: stad2562

Barless flocculent galaxies: a dynamic puzzle

Authors: Daria Zakharova, Natalia Ya. Sotnikova, Anton A. Smirnov, Sergey S. Savchenko

Abstract: We draw attention to the bright galaxies that do not show a bar in their structure but have a flocculent spiral structure. Using the THINGS' and HERACLES' kinematic data for four barless galaxies (NGC~2841, NGC~3512, NGC~5055, NGC~7331) we built their mass models including dark halos. We concluded that the fraction of the dark matter does not exceed 50\% within the optical radii of the galaxies. T… ▽ More We draw attention to the bright galaxies that do not show a bar in their structure but have a flocculent spiral structure. Using the THINGS' and HERACLES' kinematic data for four barless galaxies (NGC~2841, NGC~3512, NGC~5055, NGC~7331) we built their mass models including dark halos. We concluded that the fraction of the dark matter does not exceed 50\% within the optical radii of the galaxies. This is too little to explain the lack of a bar in these galaxies. In an attempt to understand the featureless structure of these galaxies we constructed several $N$-body models with an initially reduced content of dark matter. We concluded that, in addition to the low mass of the dark halo, the decisive factor that leads to a barless disc is the start from an initially unstable state. An isolated dynamically cold disc (with the Toomre parameter $Q < 0.5$) settled into rotational equilibrium passes trough the short stage of violent instability with fragmentation and formation of stellar clumps. After that, it evolves passively and ends up with a featureless structure. We assume that the barless flocculent galaxies studied in the present work may be descendants of galaxies at high redshifts with rotation curves which are consistent with the high mass fraction of baryons relative to the total dark matter halo. △ Less

Submitted 4 September, 2023; originally announced September 2023.

Comments: 18 pages

Globular clusters in the central region of the Milky Way galaxy I. Bar influence on the orbit parameters according to Gaia EDR3

Authors: A. T. Bajkova, A. A. Smirnov, V. V. Bobylev

Abstract: The work is devoted to the analysis of the influence of the galactic bar on the orbital motion of globular clusters in the central region of the Galaxy. For this task, 45 globular clusters were selected, 34 of which belong to the bulge/bar and 11 to the disk. The most accurate astrometric data from the Gaia satellite (Vasiliev and Baumgardt, 2021), as well as new refined average distances (Baumgar… ▽ More The work is devoted to the analysis of the influence of the galactic bar on the orbital motion of globular clusters in the central region of the Galaxy. For this task, 45 globular clusters were selected, 34 of which belong to the bulge/bar and 11 to the disk. The most accurate astrometric data from the Gaia satellite (Vasiliev and Baumgardt, 2021), as well as new refined average distances (Baumgardt and Vasiliev, 2021), were used to form the 6D-phase space required for orbit integration. The orbits of globular clusters are obtained both in an axisymmetric potential and in a potential including a bar. In this case, the mass, rotation velocity, shape and scale length of the bar were varied. A comparison is made of such orbital parameters as apocentric and pericentric distances, eccentricity and maximum distance from the galactic plane. It is shown that the mass of the bar exerts the greatest influence on the orbital motion, which is expressed mainly in an increase in both the apocentric and pericentric distances in the vast majority of globular clusters. The eccentricities of the orbits in the overwhelming majority also change significantly, and there is a change both upward and downward, especially in the range of values from 0.2 to 0.8. The greatest changes in parameters are observed in globular clusters with high radial velocities and small pericentric distances. The change in orbital parameters depending on the bar rotation velocity is less pronounced. The influence of the geometric parameters of the bar is insignificant in the accepted range of their changes. Several examples show that globular clusters in the bulge are more affected by the bar than those belonging to the disk. △ Less

Submitted 20 July, 2023; v1 submitted 8 May, 2023; originally announced May 2023.

Comments: 29 pages, 10 figures, 6 tables

Journal ref: Publications of the Pulkovo Observatory (ISSN 0367--7966), Issue 228, 2023

A Three-Dimensional Analytical Model of the Interstellar Extinction within the Nearest Kiloparsec

Authors: G. A. Gontcharov, A. V. Mosenkov, S. S. Savchenko, V. B. Il'in, A. A. Marchuk, A. A. Smirnov, P. A. Usachev, D. M. Polyakov, N. Hebdon

Abstract: We present a new version of our analytical model of the spatial interstellar extinction variations within the nearest kiloparsec. This model treats the 3D dust distribution as a superposition of three overlapping layers: (1) the layer along the Galactic midplane, (2) the layer in the Gould Belt, and (3) the layer passing through the Cepheus and Chamaeleon dust cloud complexes. In each layer the du… ▽ More We present a new version of our analytical model of the spatial interstellar extinction variations within the nearest kiloparsec. This model treats the 3D dust distribution as a superposition of three overlapping layers: (1) the layer along the Galactic midplane, (2) the layer in the Gould Belt, and (3) the layer passing through the Cepheus and Chamaeleon dust cloud complexes. In each layer the dust density decreases exponentially with increasing distance from the midplane of the layer. In addition, there are sinusoidal longitudinal extinction variations along the midplane of each layer. We have found the most probable values of 29 parameters of our model using four data sets: the 3D stellar reddening maps by Gontcharov (2017), Lallement et al. (2019), and Green et al. (2019) and the extinctions inferred by Anders et al. (2022) for 993291 giants from the Gaia EDR3. All of the data give similar estimates of the model parameters. The extinction for a star or a point in space is predicted by our model with an accuracy from $σ(A_\mathrm{V})=0.07$ to 0.37 for high and low Galactic latitudes, respectively. The natural fluctuations of the dust medium dominate in these values. When ignoring the fluctuations of the medium, the average extinction for an extended object (a galaxy, a star cluster, a dust cloud) or a small region of space is predicted by our model with an accuracy from $σ(A_\mathrm{V})=0.04$ to 0.15 for high and low Galactic latitudes, respectively. Green et al. (2019) and Anders et al. (2022) give in unison an estimate of $A_\mathrm{V}=0.12^m$ for the extinction at high latitudes across the whole Galactic dust half-layer above or below the Sun with the natural fluctuations of the medium $σ(A_\mathrm{V})=0.06^m$. Our model is a step to explain how a substantial amount of dust ended up far from the Galactic midplane. △ Less

Submitted 7 February, 2023; originally announced February 2023.

Comments: Published in Astronomy Letters, 2022, 48, 578-600

Journal ref: Astronomy Letters, 2022, 48, 578-600

Prospects for future studies using deep imaging: Analysis of individual Galactic cirrus filaments

Authors: Anton A. Smirnov, Sergey S. Savchenko, Denis M. Poliakov, Alexander A. Marchuk, Aleksandr V. Mosenkov, Vladimir B. Ilin, George A. Gontcharov, Javier Roman, Jonah Seguine

Abstract: The presence of Galactic cirrus is an obstacle for studying both faint objects in our Galaxy and low surface brightness extragalactic structures. With the aim of studying individual cirrus filaments in SDSS Stripe 82 data, we develop techniques based on machine learning and neural networks that allow one to isolate filaments from foreground and background sources in the entirety of Stripe 82 with… ▽ More The presence of Galactic cirrus is an obstacle for studying both faint objects in our Galaxy and low surface brightness extragalactic structures. With the aim of studying individual cirrus filaments in SDSS Stripe 82 data, we develop techniques based on machine learning and neural networks that allow one to isolate filaments from foreground and background sources in the entirety of Stripe 82 with a precision similar to that of the human expert. Our photometric study of individual filaments indicates that only those brighter than 26 mag arcsec$^{-2}$ in the SDSS $r$ band are likely to be identified in SDSS Stripe~82 data by their distinctive colours in the optical bands. We also show a significant impact of data processing (e.g. flat-fielding, masking of bright stars, and sky subtraction) on colour estimation. Analysing the distribution of filaments' colours with the help of mock simulations, we conclude that most filaments have colours in the following ranges: $0.55 \leq g-r \leq 0.73$ and $0.01 \leq r-i \leq 0.33$. Our work provides a useful framework for an analysis of all types of low surface brightness features (cirri, tidal tails, stellar streams, etc.) in existing and future deep optical surveys. For practical purposes, we provide the catalogue of dust filaments. △ Less

Submitted 29 January, 2023; originally announced January 2023.

Comments: Accepted for publication in MNRAS

Interstellar Extinction in Galactic Cirri in SDSS Stripe 82

Authors: G. A. Gontcharov, A. V. Mosenkov, S. S. Savchenko, V. B. Il'in, A. A. Marchuk, A. A. Smirnov, P. A. Usachev, D. M. Polyakov, Z. Shakespear

Abstract: We have applied the method of star counts with Wolf diagrams to determine the interstellar extinction in five Galactic cirri in Sloan Digital Sky Survey (SDSS) Stripe 82. For this purpose, we have used the photometry of stars in the GALEX NUV filter and the photometry of red dwarfs in five SDSS bands and four SkyMapper Southern Sky Survey DR2 bands. We have identified the cirri as sky regions with… ▽ More We have applied the method of star counts with Wolf diagrams to determine the interstellar extinction in five Galactic cirri in Sloan Digital Sky Survey (SDSS) Stripe 82. For this purpose, we have used the photometry of stars in the GALEX NUV filter and the photometry of red dwarfs in five SDSS bands and four SkyMapper Southern Sky Survey DR2 bands. We have identified the cirri as sky regions with an enhanced infrared emission from the Schlegel+1998 map. The extinction in them has been calculated relative to the nearby comparison regions with a reduced emission. The results for different filters agree well, giving the range of distances and the extinction law for each cirrus. The distances in the range 140--415 pc found are consistent with the 3D reddening maps. In the range between the $B$ and $V$ filters the extinctions found are consistent with the estimates from Schlegel+1998 for the Cardelli+1989 extinction law with $R_\mathrm{V}=3.1$. However, the extinctions found for all of the filters are best described not by the Cardelli+1989 extinction law with some $R_\mathrm{V}=3.1$, but by the inverse proportionality of the extinction and wavelength with its own coefficient for each cirrus. In one of the cirri our results suggest a very slight decrease in extinction with wavelength, i.e., a large contribution of gray extinction. In the remaining cirri a manifestation of gray extinction is not ruled out either. This is consistent with the previous measurements of the extinction law far from the Galactic midplane. △ Less

Submitted 25 January, 2023; originally announced January 2023.

Comments: 4 figures, published in Astronomy Letter

Journal ref: Astronomy Letters, 2022, 48, pp. 503-516

Isochrone fitting of Galactic globular clusters -- IV. NGC6362 and NGC6723

Authors: George A. Gontcharov, Maxim Yu. Khovritchev, Aleksandr V. Mosenkov, Vladimir B. Il'in, Alexander A. Marchuk, Denis M. Poliakov, Olga S. Ryutina, Sergey S. Savchenko, Anton A. Smirnov, Pavel A. Usachev, Jae-Woo Lee, Conner Camacho, Noah Hebdon

Abstract: We present new isochrone fits to the colour-magnitude diagrams of the Galactic globular clusters NGC\,6362 and NGC\,6723. We utilize 22 and 26 photometric filters for NGC\,6362 and NGC\,6723, respectively, from the ultraviolet to mid-infrared using data sets from {\it HST}, {\it Gaia}, unWISE, and other photometric sources. We use models and isochrones from the Dartmouth Stellar Evolution Database… ▽ More We present new isochrone fits to the colour-magnitude diagrams of the Galactic globular clusters NGC\,6362 and NGC\,6723. We utilize 22 and 26 photometric filters for NGC\,6362 and NGC\,6723, respectively, from the ultraviolet to mid-infrared using data sets from {\it HST}, {\it Gaia}, unWISE, and other photometric sources. We use models and isochrones from the Dartmouth Stellar Evolution Database (DSED) and Bag of Stellar Tracks and Isochrones (BaSTI) for $α$-enhanced [$α$/Fe]$=+0.4$ and different helium abundances. The metallicities [Fe/H]$=-1.04\pm0.07$ and $-1.09\pm0.06$ are derived from the red giant branch slopes in our fitting for NGC\,6362 and NGC\,6723, respectively. They agree with spectroscopic estimates from the literature. We find a differential reddening up to $ΔE(B-V)=0.13$ mag in the NGC\,6723 field due to the adjacent Corona Australis cloud complex. We derive the following for NGC\,6362 and NGC\,6723, respectively: distances $7.75\pm0.03\pm0.15$ (statistic and systematic error) and $8.15\pm0.04\pm0.15$ kpc; ages $12.0\pm0.1\pm0.8$ and $12.4\pm0.1\pm0.8$ Gyr; extinctions $A_\mathrm{V}=0.19\pm0.04\pm0.06$ and $0.24\pm0.03\pm0.06$ mag; reddenings $E(B-V)=0.056\pm0.01\pm0.02$ and $0.068\pm0.01\pm0.02$ mag. DSED provides systematically lower [Fe/H] and higher reddenings than BaSTI. However, the models agree in their relative estimates: NGC\,6723 is $0.44\pm0.04$ kpc further, $0.5\pm0.1$ Gyr older, $ΔE(B-V)=0.007\pm0.002$ more reddened, and with $0.05\pm0.01$ dex lower [Fe/H] than NGC\,6362. The lower metallicity and greater age of NGC\,6723 with respect to NGC\,6362 explain their horizontal branch morphology difference. This confirms age as the second parameter for these clusters. We provide lists of the cluster members from the {\it Gaia} Data Release 3. △ Less

Submitted 22 November, 2022; originally announced November 2022.

Comments: Accepted to be published in MNRAS, 20 pages, 3 online tables

The distribution of dust in edge-on galaxies: I. The global structure

Authors: Aleksandr V. Mosenkov, Pavel A. Usachev, Zacory Shakespear, Jacob Guerrette, Maarten Baes, Simone Bianchi, Emmanuel M. Xilouris, George A. Gontcharov, Vladimir B. Il'in, Alexander A. Marchuk, Sergey S. Savchenko, Anton A. Smirnov

Abstract: In this first paper in a series we present a study of the global dust emission distribution in nearby edge-on spiral galaxies. Our sample consists of 16 angularly large and 13 less spatially resolved galaxies selected from the DustPedia sample. To explore the dust emission distribution, we exploit the Herschel photometry in the range 100-500 $μ$m. We employ Sérsic and three-dimensional disc models… ▽ More In this first paper in a series we present a study of the global dust emission distribution in nearby edge-on spiral galaxies. Our sample consists of 16 angularly large and 13 less spatially resolved galaxies selected from the DustPedia sample. To explore the dust emission distribution, we exploit the Herschel photometry in the range 100-500 $μ$m. We employ Sérsic and three-dimensional disc models to fit the observed two-dimensional profiles of the galaxies. Both approaches give similar results. Our analysis unequivocally states the case for the presence of extraplanar dust in between 6 to 10 large galaxies. The results reveal that both the disc scale length and height increase as a function of wavelength between 100 and 500 $μ$m. The dust disc scale height positively correlates with the dust disc scale length, similar to what is observed for the stellar discs. We also find correlations between the scale lengths and scale heights in the near- and far-infrared which suggest that the stellar discs and their dust counterparts are tightly connected. Furthermore, the intrinsic flattening of the dust disc is inversely proportional to the maximum rotation velocity and the dust mass of the galaxy: more massive spiral galaxies host, on average, relatively thinner dust discs. Also, there is a tendency for the dust-to-stellar scale height ratio to decrease with the dust mass and rotation velocity. We conclude that low-mass spiral galaxies host a diffuse, puffed-up dust disc with a thickness similar to that of the stellar disc. △ Less

Submitted 23 July, 2022; originally announced July 2022.

Comments: 22 pages, 14 figures, accepted for publication in MNRAS

B/PS bulges in DESI Legacy edge-on galaxies I: Sample building

Authors: Alexander A. Marchuk, Anton A. Smirnov, Natalia Y. Sotnikova, Dmitriy A. Bunakalya, Vladimir P. Reshetnikov, Pavel A. Usachev, Iliya S. Tikhonenko, Viktor D. Zozulia, Daria A. Zakharova

Abstract: We present the biggest up-to-date sample of edge-on galaxies with B/PS bulges and X-structures. The sample was prepared using images from the DESI Legacy catalogue and contains about 2000 galaxies. To find suitable candidates in catalogue, we made the assumption that the residues (original images minus model) of galaxies with B/PS bulges should exhibit a characteristic X-shape. Galaxies with such… ▽ More We present the biggest up-to-date sample of edge-on galaxies with B/PS bulges and X-structures. The sample was prepared using images from the DESI Legacy catalogue and contains about 2000 galaxies. To find suitable candidates in catalogue, we made the assumption that the residues (original images minus model) of galaxies with B/PS bulges should exhibit a characteristic X-shape. Galaxies with such features were selected by eye and then used as input data for a neural network training, which was applied to a bigger sample of edge-on galaxies. Using the available data and the photometric models from the literature, we investigated the observational and statistical properties of the sample created. Comparing the $B/D$ ratios for galaxies with and without B/PS bulges, we found that the $B/D$ ratio for galaxies from our sample is statistically higher, with typical values in the range $\approx 0.2-0.5$ depending on the decomposition procedure. We studied how the opening angles $\varphi$ of the X-structure and the length of its rays are distributed in the formed sample and found them to be consistent with previous measurements and predictions from $N$-body models, e.g. $\varphi \gtrsim 25~°$, but measured here for a much larger number of galaxies. We found a sharp increase in the B/PS bulge fraction for stellar masses $\log M_{\star} \gtrsim 10.4$, but for edge-on galaxies, which complements the results of previous works. The sample can be used in future work to test various bar models and their relationship with B/PS bulges, as well as to study their stability and evolution. △ Less

Submitted 2 March, 2022; originally announced March 2022.

Fractal dimension of optical cirrus in Stripe82

Authors: Alexander A. Marchuk, Anton A. Smirnov, Aleksandr V. Mosenkov, Vladimir B. Ilín, George A. Gontcharov, Sergey S. Savchenko, Javier Román

Abstract: The geometric characteristics of dust clouds provide important information on the physical processes that structure such clouds. One of such characteristics is the $2D$ fractal dimension $D$ of a cloud projected onto the sky plane. In previous studies, which were mostly based on infrared (IR) data, the fractal dimension of individual clouds was found to be in a range from 1.1 to 1.7 with a preferr… ▽ More The geometric characteristics of dust clouds provide important information on the physical processes that structure such clouds. One of such characteristics is the $2D$ fractal dimension $D$ of a cloud projected onto the sky plane. In previous studies, which were mostly based on infrared (IR) data, the fractal dimension of individual clouds was found to be in a range from 1.1 to 1.7 with a preferred value of 1.2--1.4. In the present work, we use data from Stripe82 of the Sloan Digital Sky Survey to measure the fractal dimension of the cirrus clouds. This is done here for the first time for optical data with significantly better resolution as compared to IR data. To determine the fractal dimension, the perimeter-area method is employed. We also consider IR (IRAS and Herschel) counterparts of the corresponding optical fields to compare the results between the optical and IR. We find that the averaged fractal dimension across all clouds in the optical is $\langle D \rangle =1.69^{+0.05}_{-0.05}$ which is significantly larger then the fractal dimension of its IR counterparts $\langle D\rangle=1.38^{+0.07}_{-0.06}$. We examine several reasons for this discrepancy (choice of masking and minimal contour level, image and angular resolution, etc.) and find that for approximately half of our fields the different angular resolution (point spread function) of the optical and IR data can explain the difference between the corresponding fractal dimensions. For the other half of the fields, the fractal dimensions of the IR and visual data remain inconsistent, which can be associated with physical properties of the clouds, but further physical simulations are required to prove it. △ Less

Submitted 28 September, 2021; originally announced September 2021.

Isochrone fitting of Galactic globular clusters -- III. NGC\,288, NGC\,362, and NGC\,6218 (M12)

Authors: George A. Gontcharov, Maxim Yu. Khovritchev, Aleksandr V. Mosenkov, Vladimir B. Il'in, Alexander A. Marchuk, Sergey S. Savchenko, Anton A. Smirnov, Pavel A. Usachev, Denis M. Poliakov

Abstract: We present new isochrone fits to colour-magnitude diagrams of the Galactic globular clusters NGC\,288, NGC\,362, and NGC\,6218 (M12). We utilize a lot of photometric bands from the ultraviolet to mid-infrared by use of data from the {\it HST}, {\it Gaia}, unWISE, Pan-STARRS, and other photometric sources. In our isochrone fitting we use theoretical models and isochrones from the Dartmouth Stellar… ▽ More We present new isochrone fits to colour-magnitude diagrams of the Galactic globular clusters NGC\,288, NGC\,362, and NGC\,6218 (M12). We utilize a lot of photometric bands from the ultraviolet to mid-infrared by use of data from the {\it HST}, {\it Gaia}, unWISE, Pan-STARRS, and other photometric sources. In our isochrone fitting we use theoretical models and isochrones from the Dartmouth Stellar Evolution Program and Bag of Stellar Tracks and Isochrones for $α$-enhanced abundance [$α$/Fe]$=+0.40$, different helium abundances, and a metallicity of about [Fe/H]$=-1.3$ adopted from the literature. We derive the most probable distances $8.96\pm0.05$, $8.98\pm0.06$, and $5.04\pm0.05$ kpc, ages $13.5\pm1.1$, $11.0\pm0.6$, and $13.8\pm1.1$ Gyr, extinctions $A_\mathrm{V}=0.08\pm0.03$, $0.11\pm0.04$, and $0.63\pm0.03$ mag, and reddenings $E(B-V)=0.014\pm0.010$, $0.028\pm0.011$, and $0.189\pm0.010$ mag for NGC\,288, NGC\,362, and NGC\,6218, respectively. The distance estimates from the different models are consistent, while those of age, extinction, and reddening are not. The uncertainties of age, extinction, and reddening are dominated by some intrinsic systematic differences between the models. However, the models agree in their relative age estimates: NGC\,362 is $2.6\pm0.5$ Gyr younger than NGC\,288 and $2.8\pm0.5$ Gyr younger than NGC\,6218, confirming age as the second parameter for these clusters. We provide reliable lists of the cluster members and precise cluster proper motions from the {\it Gaia} Early Data Release 3. △ Less

Submitted 27 September, 2021; originally announced September 2021.

Comments: 19 pages, 16 figures, accepted to MNRAS

Study of the Influence of an Evolving Galactic Potential on the Orbital Properties of 152 Globular clusters with Data from the Gaia EDR3 Catalogue

Authors: A. T. Bajkova, A. A. Smirnov, V. V. Bobylev

Abstract: We have studied the influence of an evolving gravitational potential of the Milky Way Galaxy on the orbital motion of 152 globular clusters with proper motions from the Gaia EDR3 catalogue and mean distances from Baumgardt and Vasiliev (2021). To construct a semicosmological evolving model potential with changing masses and sizes of the Galactic components, we have used the algorithm described in… ▽ More We have studied the influence of an evolving gravitational potential of the Milky Way Galaxy on the orbital motion of 152 globular clusters with proper motions from the Gaia EDR3 catalogue and mean distances from Baumgardt and Vasiliev (2021). To construct a semicosmological evolving model potential with changing masses and sizes of the Galactic components, we have used the algorithm described in Haghi et al. (2015). The adopted axisymmetric three-component model potential of the Galaxy includes a spherical bulge, a flat Miyamoto--Nagai disk, and a spherical Navarro--Frenk--White dark matter halo.The orbits are integrated backward in time. We compare the orbital parameters of globular clusters derived in static and evolving potentials when integrating the orbits for 5 and 12 Gyr backward. For the first time we have studied the influence of separately a change in the masses and a change in the sizes of the Galactic components. The changes in the masses and sizes of the components are shown to act on the orbital parameters in the opposite way. At small Galactocentric distances this influence is maximally compensated for. The orbits of distant globular clusters and those with a large apocenter distance undergo the biggest changes. We show that on time scales up to $-5$~Gyr the orbits of globular clusters in the case of a potential with both changing masses and changing sizes of the components undergo, on average, minor changes compared to the case of a static potential. These changes fit into the limits of the statistical uncertainties caused by the errors in the data. So, on these time scales the Galactic potential may be deemed static. We provide tables with the orbital parameters of globular clusters derived in both static and evolving potentials. △ Less

Submitted 22 September, 2021; originally announced September 2021.

Comments: 21 pages, 10 figures, 3 tables

Journal ref: Astronomy Letters, 2021, Vol. 47, No 7, pp. 454--473

The structure of the Milky Way based on unWISE 3.4$μ$m integrated photometry

Authors: Aleksandr V. Mosenkov, Sergey S. Savchenko, Anton A. Smirnov, Peter Camps

Abstract: We present a detailed analysis of the Galaxy structure using an unWISE wide-field image at $3.4μ$m. We perform a 3D photometric decomposition of the Milky Way taking into account i) the projection of the Galaxy on the celestial sphere and ii) that the observer is located within the Galaxy at the solar radius. We consider a large set of photometric models starting with a pure disc model and ending… ▽ More We present a detailed analysis of the Galaxy structure using an unWISE wide-field image at $3.4μ$m. We perform a 3D photometric decomposition of the Milky Way taking into account i) the projection of the Galaxy on the celestial sphere and ii) that the observer is located within the Galaxy at the solar radius. We consider a large set of photometric models starting with a pure disc model and ending with a complex model which consists of thin and thick discs plus a boxy-peanut-shaped bulge. In our final model, we incorporate many observed features of the Milky Way, such as the disc flaring and warping, several over-densities in the plane, and the dust extinction. The model of the bulge with the corresponding X-shape structure is obtained from N-body simulations of a Milky Way-like galaxy. This allows us to retrieve the parameters of the aforementioned stellar components, estimate their contribution to the total Galaxy luminosity, and constrain the position angle of the bar. The mass of the thick disc in our models is estimated to be 0.4-1.3 of that for the thin disc. The results of our decomposition can be directly compared to those obtained for external galaxies via multicomponent photometric decomposition. △ Less

Submitted 23 August, 2021; originally announced August 2021.

Comments: 19 pages, 10 figures, Accepted for publication in MNRAS

First direct identification of the barlens vertical structure in galaxy models

Authors: Iliya S. Tikhonenko, Anton A. Smirnov, Natalia Ya. Sotnikova

Abstract: Applying spectral dynamics methods to one typical $N$-body model with a barlens, we dissect the modelled bar into separate components supported by completely different types of orbits. We identify at least four components: a narrow elongated bar, a boxy bar, and two components contributing to the barlens. We analyse the vertical structure of all components that make up the thick part of the bar, w… ▽ More Applying spectral dynamics methods to one typical $N$-body model with a barlens, we dissect the modelled bar into separate components supported by completely different types of orbits. We identify at least four components: a narrow elongated bar, a boxy bar, and two components contributing to the barlens. We analyse the vertical structure of all components that make up the thick part of the bar, which has a boxy/peanut shape (B/P bulge). We show that the `peanut' shape is mainly due to the orbits that assemble the boxy part of the face-on bar. We associate the X-shape with the narrow and elongated bar. The wider part of the barlens with square-like isophotes contributes to the boxy shape of the B/P bulge when we observe the galaxy edge-on. However, the part of the barlens with rounded isophotes in the face-on view is a rather flat structure in the vertical direction without any significant off-centre protrusions. Thus, for the first time, we demonstrate that the rounded face-on barlens cannot be entirely associated with the B/P bulge. △ Less

Submitted 3 April, 2021; v1 submitted 4 March, 2021; originally announced March 2021.

Comments: 5 pages, 5 figures, accepted to Astronomy and Astrophysics on March 26, 2021

Journal ref: A&A 648, L4 (2021)

New X-shaped bulge photometric model as a tool for measuring B/PS bulges and their X-structures in photometric studies

Authors: Anton A. Smirnov, Sergey S. Savchenko

Abstract: Recent orbital studies of 3D bar structure in various numerical and analytical models show that X-structures that reside in boxy/peanut-shaped (B/PS) bulges are not delineated by some specific type of orbits, but are natural parts of them and formed by the same orbits that constitute such bulges. This implies that to accurately account for B/PS bulges and their X-structures in photometric studies,… ▽ More Recent orbital studies of 3D bar structure in various numerical and analytical models show that X-structures that reside in boxy/peanut-shaped (B/PS) bulges are not delineated by some specific type of orbits, but are natural parts of them and formed by the same orbits that constitute such bulges. This implies that to accurately account for B/PS bulges and their X-structures in photometric studies, one needs the photometric model of B/PS bulge that includes an X-structure as its natural part. To find such a model, we considered a self-consistent numerical galaxy model where a typical B/PS bulge arises. Using spectral characteristics of particle-"stars", we decomposed the galaxy model onto the bar and non-bar components. We used the extracted 3D bar component to find an appropriate B/PS bulge photometric model, which can account for X-structures residing in such bulges. The resulted B/PS bulge photometric model has a truncated 2D Sersic profile with truncations introduced above (in the upper half-plane) and below (in the bottom half-plane) the rays of X-structures. We applied this model to represent B/PS bulges of various numerical models and some real galaxies. The comparison with previous works revealed that there are systematic shifts between the X-structure parameters of the same galaxies measured within the different approaches. We found that the geometric parameters of X-structures of real and modelled galaxies are consistent with each other if we measure them using our new model. △ Less

Submitted 28 September, 2020; v1 submitted 23 July, 2020; originally announced July 2020.

Comments: Accepted in MNRAS

Face-on structure of barlenses and boxy bars: an insight from spectral dynamics

Authors: Anton A. Smirnov, Iliya Tikhonenko, Natalia Ya. Sotnikova

Abstract: Based on the spectral analysis of individual orbits of stars from different $N$-body models, we show that the face-on morphology of the so-called `face-on peanut' bars (boxy bars) and barlenses is supported by different types of orbits. For `face-on peanut' bars, the so-called boxy orbits come to the fore, and they are responsible for the unusual morphology of the bar in the central regions. In th… ▽ More Based on the spectral analysis of individual orbits of stars from different $N$-body models, we show that the face-on morphology of the so-called `face-on peanut' bars (boxy bars) and barlenses is supported by different types of orbits. For `face-on peanut' bars, the so-called boxy orbits come to the fore, and they are responsible for the unusual morphology of the bar in the central regions. In the models with compact bulges, the bars show a barlens morphology in their central parts. We found that the barlens supporting orbits come in two types, one of which gives a square-like shape and the other have a rosette-like shape in the frame co-rotating with the bar. Such a shape is typical for orbits around stable loop orbits in nearly axisymmetric potentials only slightly distorted by the bar. They were already known from some of the previous studies but their role in barlens shaping was barely investigated. Although quite simple, the rosette-like orbits are found to be the main building block of a barlens in our models. The detailed analysis of all bar orbits also allowed us to distinguish the x2 orbital family and isolate the structure supported by orbits trapped around this family. The x2 family is well-known, but, apparently, for the first time in $N$-body models we have revealed the structure it supports by means of spectral dynamics and highlight its contribution to the barlens. We found that the x2 family population increases with an increase in central matter concentration. △ Less

Submitted 9 March, 2021; v1 submitted 17 July, 2020; originally announced July 2020.

Comments: Published in MNRAS

Journal ref: MNRAS, Volume 502, Issue 4, April 2021, Pages 4689-4707

Tilted outer and inner structures in edge-on galaxies?

Authors: Aleksandr V. Mosenkov, Anton A. Smirnov, Olga K. Sil'chenko, R. Michael Rich, Vladimir P. Reshetnikov, John Kormendy

Abstract: Tilted and warped discs inside tilted dark matter haloes are predicted from numerical and semi-analytical studies. In this paper, we use deep imaging to demonstrate the likely existence of tilted outer structures in real galaxies. We consider two SB0 edge-on galaxies, NGC4469 and NGC4452, which exhibit apparent tilted outer discs with respect to the inner structure. In NGC4469, this structure has… ▽ More Tilted and warped discs inside tilted dark matter haloes are predicted from numerical and semi-analytical studies. In this paper, we use deep imaging to demonstrate the likely existence of tilted outer structures in real galaxies. We consider two SB0 edge-on galaxies, NGC4469 and NGC4452, which exhibit apparent tilted outer discs with respect to the inner structure. In NGC4469, this structure has a boxy shape, inclined by $Δ$PA$\approx$3$^{\circ}$ with respect to the inner disc, whereas NGC4452 harbours a discy outer structure with $Δ$PA$\approx$6$^{\circ}$. In spite of the different shapes, both structures have surface brightness profiles close to exponential and make a large contribution ($\sim30$%) to the total galaxy luminosity. In the case of NGC4452, we propose that its tilted disc likely originates from a former fast tidal encounter (probably with IC3381). For NGC4469, a plausible explanation may also be galaxy harassment, which resulted in a tilted or even a tumbling dark matter halo. A less likely possibility is accretion of gas-rich satellites several Gyr ago. New deep observations may potentially reveal more such galaxies with tilted outer structures, especially in clusters. We also consider galaxies, mentioned in the literature, where a central component (a bar or a bulge) is tilted with respect to the stellar disc. According to our numerical simulations, one of the plausible explanations of such observed "tilts" of the bulge/bar is a projection effect due to a not exactly edge-on orientation of the galaxy coupled with a skew angle of the triaxial bulge/bar. △ Less

Submitted 26 June, 2020; originally announced June 2020.

Comments: 20 pages, 19 figures, Accepted for publication in MNRAS

Orbital ingredients for cooking X-structures in edge-on galaxies

Authors: Hanna D. Parul, Anton A. Smirnov, Natalia Ya. Sotnikova

Abstract: X-structures are often observed in galaxies hosting the so-called B/PS (boxy/peanuts) bulges and are visible from the edge-on view. They are the most notable features of B/PS bulges and appear as four rays protruding from the disk of the host galaxy and distinguishable against the B/PS bulge background. In some works their origin is thought to be connected with the so-called banana-shaped orbits w… ▽ More X-structures are often observed in galaxies hosting the so-called B/PS (boxy/peanuts) bulges and are visible from the edge-on view. They are the most notable features of B/PS bulges and appear as four rays protruding from the disk of the host galaxy and distinguishable against the B/PS bulge background. In some works their origin is thought to be connected with the so-called banana-shaped orbits with a vertical resonance 2:1. A star in such an orbit performs two oscillations in the vertical direction per one revolution in the bar frame. Several recent studies that analyzed ensembles of orbits arising in different $N$-body models do not confirm the dominance of the resonant 2:1 orbits in X-structures. In our work we analyze two $N$-body models and show how the X-structure in our models is gradually assembled from the center to the periphery from orbits with less than 2:1 frequency ratio. The most number of such orbits is enclosed in a 'farfalle'-shape (Italian pasta) form and turns out to be non-periodic. We conclude that the X-structure is only the envelope of regions of high density caused by the crossing or folding of different types of orbits at their highest points, and does not have a "backbone" similar to that of the in-plane bar. Comparing the orbital structure of two different numerical models, we show that the dominance of one or another family of orbits with a certain ratio of the vertical oscillations frequency to the in-plane frequency depends on the parameters of the underlying galaxy and ultimately determines the morphology of the X-structure and the opening angle of its rays. △ Less

Submitted 16 February, 2020; originally announced February 2020.

Comments: Accepted for publication in ApJ

Is the late buckling stage inevitable in the bar life?

Authors: Anton A. Smirnov, Natalia Ya. Sotnikova

Abstract: By means of self-consistent numerical simulations we investigated the dynamical impact of classical bulges on the growth of the secondary buckling of a bar. Overall we considered 14 models with different disc and bulge parameters. We obtained that a bulge with a quite modest mass $B/D=0.1$ leads to completely symmetrical evolution of the bar almost independently of the initial stellar disc paramet… ▽ More By means of self-consistent numerical simulations we investigated the dynamical impact of classical bulges on the growth of the secondary buckling of a bar. Overall we considered 14 models with different disc and bulge parameters. We obtained that a bulge with a quite modest mass $B/D=0.1$ leads to completely symmetrical evolution of the bar almost independently of the initial stellar disc parameters and even can damp the first bending. At the same time, the bars in all our bulgeless models suffer from the short primary and prolonged secondary buckling. Given the smallness of the mass suppressing secondary buckling, we conclude that a classical bulge along with the gas central concentration may be the main culprits for the rarity of bars with ongoing buckling in the local Universe. △ Less

Submitted 22 February, 2019; originally announced February 2019.

Comments: Accepted for publication in MNRAS

What determines the flatness of X-shaped structures in edge-on galaxies

Authors: Anton A. Smirnov, Natalia Ya. Sotnikova

Abstract: Recent observational studies of X-shaped structures revealed that values of their opening angles lie in a narrow range: from 20 deg to 43 deg with smaller X-shaped structures appearing to follow a characteristic opening angle ~ 27 deg - 31 deg. We use self-consistent numerical simulations to uncover what parameters of host galaxies govern the opening angle spread. We constructed a series of equili… ▽ More Recent observational studies of X-shaped structures revealed that values of their opening angles lie in a narrow range: from 20 deg to 43 deg with smaller X-shaped structures appearing to follow a characteristic opening angle ~ 27 deg - 31 deg. We use self-consistent numerical simulations to uncover what parameters of host galaxies govern the opening angle spread. We constructed a series of equilibrium models of galaxies with high spatial resolution, varying the dark halo and bulge contribution in the overall gravitational potential, the initial disc thickness of models and the Toomre parameter Q and followed their evolution for almost 8 Gyrs. Each model demonstrated the formation of clear X-structures with different flatness. We have obtained that opening angles lie in the range from 25 deg to 42 deg throughout the entire evolution. These values are roughly consistent with observational data. The greatest variation in the opening angles is obtained by varying the mass of the dark halo. The initial thickness of the disc and the Toomre parameter are responsible for smaller variations of the angle and shorter X-structures. An increase of both parameters changes the morphology of in-plane bars and X-structures. In some cases we observed even double X-structures. The main effect of the bulge is to prevent buckling at intermediate and late stages of the disc evolution. Comparison of models with different halo masses indicates that the smallest observable values of opening angles can be associated with the presence of a heavy dark halo (more than 3 masses of the disc within the optical radius). △ Less

Submitted 17 September, 2018; originally announced September 2018.

Comments: 19 pages, 36 figures, accepted for publication in MNRAS