Showing 1–50 of 223 results for author: Athanassoula, E

The response of the inner dark matter halo to stellar bars

Authors: Daniel A. Marostica, Rubens E. G. Machado, E. Athanassoula, T. Manos

Abstract: Barred galaxies constitute about two thirds of observed disc galaxies. Bars affect not only the mass distribution of gas and stars, but also that of the dark matter. An elongation of the inner dark matter halo is known as the halo bar. We aim to characterise the structure of the halo bars, with the goal of correlating them with the properties of the stellar bars. We use a suite of simulated galaxi… ▽ More Barred galaxies constitute about two thirds of observed disc galaxies. Bars affect not only the mass distribution of gas and stars, but also that of the dark matter. An elongation of the inner dark matter halo is known as the halo bar. We aim to characterise the structure of the halo bars, with the goal of correlating them with the properties of the stellar bars. We use a suite of simulated galaxies with various bar strengths, including gas and star formation. We quantify strengths, shapes, and densities of these simulated stellar bars. We carry out numerical experiments with frozen and analytic potentials in order to understand the role played by a live responsive stellar bar. We find that the halo bar generally follows the trends of the disc bar. The strengths of the halo and stellar bars are tightly correlated. Stronger bars induce a slight increase of dark matter density within the inner halo. Numerical experiments show that a non-responsive frozen stellar bar would be capable of inducing a dark matter bar, but it would be weaker than the live case by a factor of roughly two. △ Less

Submitted 27 May, 2024; originally announced May 2024.

Comments: 18 pages, 9 figures, accepted for publication in Galaxies

MHONGOOSE -- A MeerKAT Nearby Galaxy HI Survey

Authors: W. J. G. de Blok, J. Healy, F. M. Maccagni, D. J. Pisano, A. Bosma, J. English, T. Jarrett, A. Marasco, G. R. Meurer, S. Veronese, F. Bigiel, L. Chemin, F. Fraternali, B. W. Holwerda, P. Kamphuis, H. R. Klöckner, D. Kleiner, A. K. Leroy, M. Mogotsi, K. A. Oman, E. Schinnerer, L. Verdes-Montenegro, T. Westmeier, O. I. Wong, N. Zabel , et al. (35 additional authors not shown)

Abstract: The MHONGOOSE (MeerKAT HI Observations of Nearby Galactic Objects: Observing Southern Emitters) survey maps the distribution and kinematics of the neutral atomic hydrogen (HI) gas in and around 30 nearby star-forming spiral and dwarf galaxies to extremely low HI column densities. The HI column density sensitivity (3 sigma over 16 km/s) ranges from ~ 5 x 10^{17} cm^{-2} at 90'' resolution to ~4 x 1… ▽ More The MHONGOOSE (MeerKAT HI Observations of Nearby Galactic Objects: Observing Southern Emitters) survey maps the distribution and kinematics of the neutral atomic hydrogen (HI) gas in and around 30 nearby star-forming spiral and dwarf galaxies to extremely low HI column densities. The HI column density sensitivity (3 sigma over 16 km/s) ranges from ~ 5 x 10^{17} cm^{-2} at 90'' resolution to ~4 x 10^{19} cm^{-2} at the highest resolution of 7''. The HI mass sensitivity (3 sigma over 50 km/s) is ~5.5 X 10^5 M_sun at a distance of 10 Mpc (the median distance of the sample galaxies). The velocity resolution of the data is 1.4 km/s. One of the main science goals of the survey is the detection of cold, accreting gas in the outskirts of the sample galaxies. The sample was selected to cover a range in HI masses, from 10^7 M_sun to almost 10^{11} M_sun, to optimally sample possible accretion scenarios and environments. The distance to the sample galaxies ranges from 3 to 23 Mpc. In this paper, we present the sample selection, survey design, and observation and reduction procedures. We compare the integrated HI fluxes based on the MeerKAT data with those derived from single-dish measurement and find good agreement, indicating that our MeerKAT observations are recovering all flux. We present HI moment maps of the entire sample based on the first ten percent of the survey data, and find that a comparison of the zeroth- and second-moment values shows a clear separation between the physical properties of the HI in areas with star formation and areas without, related to the formation of a cold neutral medium. Finally, we give an overview of the HI-detected companion and satellite galaxies in the 30 fields, five of which have not previously been catalogued. We find a clear relation between the number of companion galaxies and the mass of the main target galaxy. △ Less

Submitted 6 June, 2024; v1 submitted 2 April, 2024; originally announced April 2024.

Comments: Accepted for publication in Astronomy & Astrophysics

Bar Properties as a Function of Wavelength: A Local Baseline with S4G for High-Redshift Studies

Authors: Karín Menéndez-Delmestre, Thiago S. Gonçalves, Kartik Sheth, Tomás Düringer Jacques de Lima, Taehyun Kim, Dimitri A. Gadotti, Eva Schinnerer, E. Athanassoula, Albert Bosma, Debra Meloy Elmegreen, Johan H. Knapen, Rubens E. G. Machado, Heikki Salo

Abstract: The redshift evolution of bars is an important signpost of the dynamic maturity of disk galaxies. To characterize the intrinsic evolution safe from band-shifting effects, it is necessary to gauge how bar properties vary locally as a function of wavelength. We investigate bar properties in 16 nearby galaxies from the Spitzer Survey of Stellar Structure in Galaxies (S4G) at ultraviolet, optical and… ▽ More The redshift evolution of bars is an important signpost of the dynamic maturity of disk galaxies. To characterize the intrinsic evolution safe from band-shifting effects, it is necessary to gauge how bar properties vary locally as a function of wavelength. We investigate bar properties in 16 nearby galaxies from the Spitzer Survey of Stellar Structure in Galaxies (S4G) at ultraviolet, optical and mid-infrared wavebands. Based on the ellipticity and position angle profiles from fitting elliptical isophotes to the two-dimensional light distribution, we find that both bar length and ellipticity - the latter often used as a proxy for bar strength - increase at bluer wavebands. Bars are 9% longer in the B-band than at 3.6 um. Their ellipticity increases typically by 8% in the B-band, with a significant fraction (>40%) displaying an increase up to 35%. We attribute the increase in bar length to the presence of star forming knots at the end of bars: these regions are brighter in bluer bands, stretching the bar signature further out. The increase in bar ellipticity could be driven by the apparent bulge size: the bulge is less prominent at bluer bands, allowing for thinner ellipses within the bar region. Alternatively, it could be due to younger stellar populations associated to the bar. The resulting effect is that bars appear longer and thinner at bluer wavebands. This indicates that band-shifting effects are significant and need to be corrected for high-redshift studies to reliably gauge any intrinsic evolution of the bar properties with redshift. △ Less

Submitted 7 December, 2023; originally announced December 2023.

Comments: Accepted to MNRAS

A new tidal scenario for double bar formation

Authors: Marcin Semczuk, Ewa L. Łokas, Adriana de Lorenzo-Cáceres, E. Athanassoula

Abstract: Double bars make up a significant fraction of barred galaxies. We propose a new formation scenario for double bars that involves tidal interactions. We demonstrate the viability of this scenario using two examples of simulated galaxies from run TNG50-1 of the IllustrisTNG project. In the proposed scenario the inner bar forms first, either in isolation, via instabilities, or through previous tides.… ▽ More Double bars make up a significant fraction of barred galaxies. We propose a new formation scenario for double bars that involves tidal interactions. We demonstrate the viability of this scenario using two examples of simulated galaxies from run TNG50-1 of the IllustrisTNG project. In the proposed scenario the inner bar forms first, either in isolation, via instabilities, or through previous tides. The outer bar forms later from the material that is tidally distorted by a strong interaction. The inner and outer bars formed this way rotate with different pattern speeds and can be mistaken for a single bar when their phases align. The double-barred structure is stable and can last for at least 3 Gyr. The inner bars of the tidally induced double bars can also have big sizes, which can possibly explain the origin of sizable inner bars recently found in some galaxies. △ Less

Submitted 29 September, 2023; originally announced September 2023.

Comments: 5 pages, 6 figures, submitted to MNRAS Letters

A JWST investigation into the bar fraction at redshifts 1 < z < 3

Authors: Zoe A. Le Conte, Dimitri A. Gadotti, Leonardo Ferreira, Christopher J. Conselice, Camila de Sá-Freitas, Taehyun Kim, Justus Neumann, Francesca Fragkoudi, E. Athanassoula, Nathan J. Adams

Abstract: The presence of a stellar bar in a disc galaxy indicates that the galaxy hosts in its main part a dynamically settled disc and that bar-driven processes are taking place in shaping its evolution. Studying the cosmic evolution of the bar fraction in disc galaxies is therefore essential to understand galaxy evolution in general. Previous studies have found, using the Hubble Space Telescope (HST), th… ▽ More The presence of a stellar bar in a disc galaxy indicates that the galaxy hosts in its main part a dynamically settled disc and that bar-driven processes are taking place in shaping its evolution. Studying the cosmic evolution of the bar fraction in disc galaxies is therefore essential to understand galaxy evolution in general. Previous studies have found, using the Hubble Space Telescope (HST), that the bar fraction significantly declines from the local Universe to redshifts near one. Using the first four pointings from the James Webb Space Telescope (JWST) Cosmic Evolution Early Release Science Survey (CEERS) and the initial public observations for the Public Release Imaging for Extragalactic Research (PRIMER), we extend the studies of the bar fraction in disc galaxies to redshifts $1 \leq z \leq 3$, i.e., for the first time beyond redshift two. We only use galaxies that are also present in the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS) on the Extended Groth Strip (EGS) and Ultra Deep Survey (UDS) HST observations. An optimised sample of 368 close-to-face-on galaxies is visually classified to find the fraction of bars in disc galaxies in two redshift bins: $1 \leq z \leq 2$ and $2 < z \leq 3$. The bar fraction decreases from $\approx 17.8^{+ 5.1}_{- 4.8}$ per cent to $\approx 13.8^{+ 6.5}_{- 5.8}$ per cent (from the lower to the higher redshift bin), but is about twice the bar fraction found using bluer HST filters. Our results show that bar-driven evolution might commence at early cosmic times and that dynamically settled discs are already present at a lookback time of $\sim 11$ Gyrs. △ Less

Submitted 23 April, 2024; v1 submitted 18 September, 2023; originally announced September 2023.

Comments: Accepted for publication in MNRAS; 18 pages, 10 figures; final version includes new disc classification and mass cut, but results remain the same; Figs. 6 and 7 summarise the main results

Roman Early-Definition Astrophysics Survey Opportunity: Galactic Roman Infrared Plane Survey (GRIPS)

Authors: Roberta Paladini, Catherine Zucker, Robert Benjamin, David Nataf, Dante Minniti, Gail Zasowski, Joshua Peek, Sean Carey, Lori Allen, Javier Alonso-Garcia, Joao Alves, Friederich Anders, Evangelie Athanassoula, Timothy C. Beers, Jonathan Bird, Joss Bland-Hwathorn, Anthony Brown, Sven Buder, Luca Casagrande, Andrew Casey, Santi Cassisi, Marcio Catelan, Ranga-Ram Chary, Andre-Nicolas Chene, David Ciardi , et al. (45 additional authors not shown)

Abstract: A wide-field near-infrared survey of the Galactic disk and bulge/bar(s) is supported by a large representation of the community of Galactic astronomers. The combination of sensitivity, angular resolution and large field of view make Roman uniquely able to study the crowded and highly extincted lines of sight in the Galactic plane. A ~1000 deg2 survey of the bulge and inner Galactic disk would yiel… ▽ More A wide-field near-infrared survey of the Galactic disk and bulge/bar(s) is supported by a large representation of the community of Galactic astronomers. The combination of sensitivity, angular resolution and large field of view make Roman uniquely able to study the crowded and highly extincted lines of sight in the Galactic plane. A ~1000 deg2 survey of the bulge and inner Galactic disk would yield an impressive dataset of ~120 billion sources and map the structure of our Galaxy. The effort would foster subsequent expansions in numerous dimensions (spatial, depth, wavelengths, epochs). Importantly, the survey would benefit from early defintion by the community, namely because the Galactic disk is a complex environment, and different science goals will require trade offs. △ Less

Submitted 14 July, 2023; originally announced July 2023.

Comments: Submitted to the Roman Project on October 22 2021 in response to a call for white papers on early-definition Astrophysics opportunity

On the origin of the Galactic thin and thick discs, their abundance gradients and the diagnostic potential of their abundance ratios

Authors: Nikos Prantzos, Carlos Abia, Tianxiang Chen, Patrick de Laverny, Alejandra Recio-Blanco, E. Athanassoula, Lorenzo Roberti, Diego Vescovi, Marco Limongi, Alessandro Chieffi, Sergio Cristallo

Abstract: Using a semi-analytical model of the evolution of the Milky Way, we show how secular evolution can create distinct overdensities in the phase space of various properties (e.g. age vs metallicity or abundance ratios vs age) corresponding to the thin and thick discs. In particular, we show how key properties of the Solar vicinity can be obtained by secular evolution, with no need for external or spe… ▽ More Using a semi-analytical model of the evolution of the Milky Way, we show how secular evolution can create distinct overdensities in the phase space of various properties (e.g. age vs metallicity or abundance ratios vs age) corresponding to the thin and thick discs. In particular, we show how key properties of the Solar vicinity can be obtained by secular evolution, with no need for external or special events, like galaxy mergers or paucity in star formation. This concerns the long established double-branch behaviour of [alpha/Fe] vs metallicity and the recently found non-monotonic evolution of the stellar abundance gradient, evaluated at the birth radii of stars. We extend the discussion to other abundance ratios and we suggest a classification scheme, based on the nature of the corresponding yields (primary vs secondary or odd elements) and on the lifetimes of their sources (short-lived vs long-lived ones). The latter property is critical in determining the single- or double- branch behavior of an elementary abundance ratio in the Solar neighborhood. We underline the high diagnostic potential of this finding, which can help to separate clearly elements with sources evolving on different timescales and help determining the site of e.g. the r-process(es). We define the "abundance distance" between the thin and thick disc sequences as an important element for such a separation. We also show how the inside-out evolution of the Milky Way disc leads rather to a single-branch behavior in other disc regions. △ Less

Submitted 10 October, 2023; v1 submitted 22 May, 2023; originally announced May 2023.

Comments: 20 pages, 16 figures, typos corrected, matches published version in MNRAS

Spiral arms and the angular momentum gap in Milky Way Cepheids

Authors: Marcin Semczuk, Walter Dehnen, Ralph Schoenrich, E. Athanassoula

Abstract: The angular-momentum distribution of classical Cepheids in the outer Milky Way disc is bi-modal with a gap at $L_\mathrm{gap}=2950\,\mathrm{km}\,\mathrm{s}^{-1}\,\mathrm{kpc}$, corresponding to $R=13\,\mathrm{kpc}$, while no similar feature has been found in the general population of disc stars. We show that star formation in multiple spiral arm segments at the same azimuth leads to such multi-mod… ▽ More The angular-momentum distribution of classical Cepheids in the outer Milky Way disc is bi-modal with a gap at $L_\mathrm{gap}=2950\,\mathrm{km}\,\mathrm{s}^{-1}\,\mathrm{kpc}$, corresponding to $R=13\,\mathrm{kpc}$, while no similar feature has been found in the general population of disc stars. We show that star formation in multiple spiral arm segments at the same azimuth leads to such multi-modality which quickly dissolves and only shows in young stars. Unlike other explanations, such as a 1:1 orbital resonance with the Galactic bar, this also accounts for the observed steepening of the stellar warp at $L_\mathrm{gap}$, since the adjacent spiral arms represent different parts of the warped gas disc. In this scenario the gap is clearly present only in young stars, as observed, while most purely stellar dynamical origins would affect all disc populations, including older disc stars. △ Less

Submitted 21 October, 2022; originally announced October 2022.

Comments: 6 pages, 7 figures, submitted to MNRAS

Age-morphology dependence of the Milky Way boxy/peanut bulge seen in Mira variables

Authors: Marcin Semczuk, Walter Dehnen, Ralph Schoenrich, E. Athanassoula

Abstract: We analyse the distribution of Mira variable stars in the central region of the Milky Way. We find that with increasing period, i.e. decreasing age, the Miras shift towards negative Galactic longitudes $\ell$. Comparing to a cosmological zoom simulation of a barred galaxy, we find that this shift with age can be explained by an age-morphology dependence of the boxy peanut/bulge. Owing to a combina… ▽ More We analyse the distribution of Mira variable stars in the central region of the Milky Way. We find that with increasing period, i.e. decreasing age, the Miras shift towards negative Galactic longitudes $\ell$. Comparing to a cosmological zoom simulation of a barred galaxy, we find that this shift with age can be explained by an age-morphology dependence of the boxy peanut/bulge. Owing to a combination of projection effects and the limitation of the range of Galactic longitudes, the near hump at $\ell>0$ is more truncated for younger populations, and the far hump at $\ell<0$ dominates the observed distributions. △ Less

Submitted 9 June, 2022; originally announced June 2022.

Comments: 5 pages, 4 figures, submitted to MNRAS

Stellar masses, sizes, and radial profiles for 465 nearby early-type galaxies: an extension to the Spitzer Survey of Stellar Structure in Galaxies (S$^{4}$G)

Authors: A. E. Watkins, H. Salo, E. Laurikainen, S. Díaz-García, S. Comerón, J. Janz, A. H. Su, R. Buta, E. Athanassoula, A. Bosma, L. C. Ho, B. W. Holwerda, T. Kim, J. H. Knapen, S. Laine, K. Menéndez-Delmestre, R. F. Peletier, K. Sheth, D. Zaritsky

Abstract: The Spitzer Survey of Stellar Structure in Galaxies (S$^{4}$G) is a detailed study of over 2300 nearby galaxies in the near-infrared (NIR), which has been critical to our understanding of the detailed structures of nearby galaxies. Because the sample galaxies were selected only using radio-derived velocities, however, the survey favored late-type disk galaxies over lenticulars and ellipticals. A f… ▽ More The Spitzer Survey of Stellar Structure in Galaxies (S$^{4}$G) is a detailed study of over 2300 nearby galaxies in the near-infrared (NIR), which has been critical to our understanding of the detailed structures of nearby galaxies. Because the sample galaxies were selected only using radio-derived velocities, however, the survey favored late-type disk galaxies over lenticulars and ellipticals. A follow-up Spitzer survey was conducted to rectify this bias, adding 465 early-type galaxies (ETGs) to the original sample, to be analyzed in a manner consistent with the initial survey. We present the data release of this ETG extension, up to the third data processing pipeline (P3): surface photometry. We produce curves of growth and radial surface brightness profiles (with and without inclination corrections) using reduced and masked Spitzer IRAC 3.6$μ$m and 4.5$μ$m images produced through Pipelines 1 and 2, respectively. From these profiles, we derive the following integrated quantities: total magnitudes, stellar masses, concentration parameters, and galaxy size metrics. We showcase NIR scaling relations for ETGs among these quantities. We examine general trends across the whole S$^{4}$G and ETG extension among our derived parameters, highlighting differences between ETGs and late-type galaxies (LTGs). ETGs are, on average, more massive and more concentrated than LTGs, and also show subtle distinctions among ETG morphological sub-types. We also derive the following scaling relations and compare with previous results in visible light: mass--size (both half-light and isophotal), mass--concentration, mass--surface brightness (central, effective, and within 1 kpc), and mass--color. We find good agreement with previous works, though some relations (e.g., mass--central surface brightness) will require more careful multi-component decompositions to be fully understood. △ Less

Submitted 20 January, 2022; originally announced January 2022.

Comments: 25 pages, 17 figures, accepted for publication in A&A

Journal ref: A&A 660, A69 (2022)

A barred Milky Way surrogate from an N-body simulation

Authors: T. Tepper-Garcia, J. Bland-Hawthorn, E. Vasiliev, E. Athanassoula, O. Gerhard, A. Quillen, P. McMillan, K. Freeman, G. F. Lewis, R. Teyssier, S. Sharma, M. R. Hayden, S. Buder

Abstract: We present an N-body model for the barred Milky Way (MW) galaxy that reproduces many of its properties, including the overall mass distribution, the disc kinematics, and the properties of the central bar. Our high-resolution (N ~ 10^8 particles) simulation, performed with the Ramses code, starts from an axisymmetric non-equilibrium configuration constructed within the AGAMA framework. This is a se… ▽ More We present an N-body model for the barred Milky Way (MW) galaxy that reproduces many of its properties, including the overall mass distribution, the disc kinematics, and the properties of the central bar. Our high-resolution (N ~ 10^8 particles) simulation, performed with the Ramses code, starts from an axisymmetric non-equilibrium configuration constructed within the AGAMA framework. This is a self-consistent dynamical model of the MW defined by the best available parameters for the dark matter halo, the stellar disc and the bulge. For the known (stellar and gas) disc mass (4.5 x 10^10 Msun) and disc mass fraction at R ~ 2.2 R_d (f_d ~ 0.3 - 0.6), the low mass limit does not yield a bar in a Hubble time. The high mass limit adopted here produces a box/peanut bar within about 2 Gyr with the correct mass (~10^10 Msun), size (~5 kpc) and peak pattern speed (~ 40-45 km/s/kpc). In agreement with earlier work, the bar formation timescale scales inversely with f_d (i.e. log [T/Gyr] ~ 0.60/f_d - 0.83 for 1 < f_d < 0.3). The disc radial heating is strong, but, in contrast to earlier claims, we find that disc vertical heating outside of the box/peanut bulge structure is negligible. The synthetic barred MW exhibits long-term stability, except for the slow decline (roughly -2 km/s/kpc/Gyr) of the bar pattern speed, consistent with recent estimates. If our model is indicative of the Milky Way, we estimate that the bar first emerged 3-4 Gyr ago. △ Less

Submitted 9 November, 2021; originally announced November 2021.

Comments: We welcome comments and missing references

Cosmic Evolution of Barred Galaxies up to z ~ 0.84

Authors: Taehyun Kim, E. Athanassoula, Kartik Sheth, Albert Bosma, Myeong-Gu Park, Yun Hee Lee, Hong Bae Ann

Abstract: We explore the cosmic evolution of the bar length, strength, and light deficit around the bar for 379 barred galaxies at 0.2 < z $\leq$ 0.835 using F814W images from the COSMOS survey. Our sample covers galaxies with stellar mass 10.0 $\leq$ log(M*/Msun) $\leq$ 11.4 and various Hubble types. The bar length is strongly related to the galaxy mass, the disk scale length (h), R50, and R90, where the l… ▽ More We explore the cosmic evolution of the bar length, strength, and light deficit around the bar for 379 barred galaxies at 0.2 < z $\leq$ 0.835 using F814W images from the COSMOS survey. Our sample covers galaxies with stellar mass 10.0 $\leq$ log(M*/Msun) $\leq$ 11.4 and various Hubble types. The bar length is strongly related to the galaxy mass, the disk scale length (h), R50, and R90, where the last two are the radii containing 50 and 90% of total stellar mass, respectively. Bar length remains almost constant, suggesting little or no evolution in bar length over the last 7 Gyrs. The normalized bar lengths (Rbar/h, Rbar/R50, and Rbar/R90) do not show any clear cosmic evolution. Also, the bar strength (A2 and Qb) and the light deficit around the bar reveal little or no cosmic evolution. The constancy of the normalized bar lengths over cosmic time implies that the evolution of bars and of disks is strongly linked over all times. We discuss our results in the framework of predictions from numerical simulations. We conclude there is no strong disagreement between our results and up-to-date simulations. △ Less

Submitted 8 September, 2021; originally announced September 2021.

Comments: 16 pages, 8 figures, Accepted for publication in ApJ

The small boxy/peanut structure of the Milky Way traced by old stars

Authors: Marcin Semczuk, Walter Dehnen, Ralph Schoenrich, E. Athanassoula

Abstract: We analyse the positions of RR Lyrae stars in the central region of the Milky Way. In addition to the overall bar shape detected previously, we find evidence for a peanut shaped structure, in form of overdensities near $\ell=-2$ deg and $1$ deg at $b\sim-3.5$ deg. The corresponding physical distance between the two peaks of the peanut is $\sim0.7$ kpc, significantly shorter than that found from ne… ▽ More We analyse the positions of RR Lyrae stars in the central region of the Milky Way. In addition to the overall bar shape detected previously, we find evidence for a peanut shaped structure, in form of overdensities near $\ell=-2$ deg and $1$ deg at $b\sim-3.5$ deg. The corresponding physical distance between the two peaks of the peanut is $\sim0.7$ kpc, significantly shorter than that found from near-IR images (3.3 kpc) and red-clump stars. Qualitatively this is expected from `fractionation' of bar orbits, which we demonstrate to be matched in a simulation of an inside-out growing disc subsequently forming a bar. △ Less

Submitted 11 November, 2021; v1 submitted 22 July, 2021; originally announced July 2021.

Comments: 6 pages, 5 figures, accepted to MNRAS

Buckling bars in nearly face-on galaxies observed with MaNGA

Authors: K. M. Xiang, D. M. Nataf, E. Athanassoula, N. L. Zakamska, K. Rowlands, K. Masters, A. Fraser-McKelvie, N. Drory, K. Kraljic

Abstract: Over half of disk galaxies are barred, yet the mechanisms for bar formation and the life-time of bar buckling remain poorly understood. In simulations, a thin bar undergoes a rapid (<1 Gyr) event called "buckling," during which the inner part of the bar is asymmetrically bent out of the galaxy plane and eventually thickens, developing a peanut/X-shaped profile when viewed side-on. Through analyzin… ▽ More Over half of disk galaxies are barred, yet the mechanisms for bar formation and the life-time of bar buckling remain poorly understood. In simulations, a thin bar undergoes a rapid (<1 Gyr) event called "buckling," during which the inner part of the bar is asymmetrically bent out of the galaxy plane and eventually thickens, developing a peanut/X-shaped profile when viewed side-on. Through analyzing stellar kinematics of N-body model snapshots of a galaxy before, during, and after the buckling phase, we confirm a distinct quadrupolar pattern of out-of-plane stellar velocities in nearly face-on galaxies. This kinematic signature of buckling allows us to identify five candidates of currently buckling bars among 434 barred galaxies in the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) Survey, an integral field unit (IFU) spectroscopic survey that measures the composition and kinematic structure of nearby galaxies. The frequency of buckling events detected is consistent with the 0.5-1 Gyr timescale predicted by simulations. The five candidates we present more than double the total number of candidate buckling bars, and are the only ones found using the kinematic signature. △ Less

Submitted 11 January, 2021; originally announced January 2021.

Comments: 22 pages, 12 figures. Accepted for publication in ApJ

MeerKAT HI commissioning observations of MHONGOOSE galaxy ESO 302-G014

Authors: W. J. G. de Blok, E. Athanassoula, A. Bosma, F. Combes, J. English, G. H. Heald, P. Kamphuis, B. S. Koribalski, G. R. Meurer, J. Román, A. Sardone, L. Verdes-Montenegro, F. Bigiel, E. Brinks, L. Chemin, F. Fraternali, T. Jarrett, D. Kleiner, F. M. Maccagni, D. J. Pisano, P. Serra, K. Spekkens, P. Amram, C. Carignan, R-J. Dettmar , et al. (21 additional authors not shown)

Abstract: We present the results of three commissioning HI observations obtained with the MeerKAT radio telescope. These observations make up part of the preparation for the forthcoming MHONGOOSE nearby galaxy survey, which is a MeerKAT large survey project that will study the accretion of gas in galaxies and the link between gas and star formation. We used the available HI data sets, along with ancillary d… ▽ More We present the results of three commissioning HI observations obtained with the MeerKAT radio telescope. These observations make up part of the preparation for the forthcoming MHONGOOSE nearby galaxy survey, which is a MeerKAT large survey project that will study the accretion of gas in galaxies and the link between gas and star formation. We used the available HI data sets, along with ancillary data at other wavelengths, to study the morphology of the MHONGOOSE sample galaxy, ESO 302-G014, which is a nearby gas-rich dwarf galaxy. We find that ESO 302-G014 has a lopsided, asymmetric outer disc with a low column density. In addition, we find a tail or filament of HI clouds extending away from the galaxy, as well as an isolated HI cloud some 20 kpc to the south of the galaxy. We suggest that these features indicate a minor interaction with a low-mass galaxy. Optical imaging shows a possible dwarf galaxy near the tail, but based on the current data, we cannot confirm any association with ESO 302-G014. Nonetheless, an interaction scenario with some kind of low-mass companion is still supported by the presence of a significant amount of molecular gas, which is almost equal to the stellar mass, and a number of prominent stellar clusters, which suggest recently triggered star formation. These data show that MeerKAT produces exquisite imaging data. The forthcoming full-depth survey observations of ESO 302-G014 and other sample galaxies will, therefore, offer insights into the fate of neutral gas as it moves from the intergalactic medium onto galaxies. △ Less

Submitted 21 September, 2020; originally announced September 2020.

Comments: Accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 643, A147 (2020)

Kinematics of RR Lyrae stars in the Galactic bulge with OGLE-IV and Gaia DR2

Authors: Hangci Du, Shude Mao, E. Athanassoula, Juntai Shen, Pawel Pietrukowicz

Abstract: We analyze the kinematics and spatial distribution of 15,599 fundamental-mode RR Lyrae (RRL) stars in the Milky Way bulge by combining OGLE-IV photometric data and Gaia DR2 proper motions. We show that the longitudinal proper motions and the line-of-sight velocities can give similar results for the rotation in the Galactic central regions. The angular velocity of bulge RRLs is found to be around… ▽ More We analyze the kinematics and spatial distribution of 15,599 fundamental-mode RR Lyrae (RRL) stars in the Milky Way bulge by combining OGLE-IV photometric data and Gaia DR2 proper motions. We show that the longitudinal proper motions and the line-of-sight velocities can give similar results for the rotation in the Galactic central regions. The angular velocity of bulge RRLs is found to be around $35$ km s$^{-1}$ kpc$^{-1}$, significantly smaller than that for the majority of bulge stars ($50-60$ km s$^{-1}$ kpc$^{-1}$); bulge RRLs have larger velocity dispersion (120$-$140 km s$^{-1}$) than younger stars. The dependence of the kinematics of the bulge RRLs on their metallicities is shown by their rotation curves and spatial distributions. Metal-poor RRLs ([Fe/H]<$-1$) show a smaller bar angle than metal-rich ones. We also find clues suggesting that RRLs in the bulge are not dominated by halo stars. These results might explain some previous conflicting results over bulge RRLs and help understand the chemodynamical evolution of the Galactic bulge. △ Less

Submitted 13 October, 2020; v1 submitted 2 July, 2020; originally announced July 2020.

Comments: 15 pages, 15 figures. Published by MNRAS

Basis function expansions for galactic dynamics: Spherical versus cylindrical coordinates

Authors: Yougang Wang, E. Athanassoula, Shude Mao

Abstract: Aims. The orbital structure of galaxies is strongly influenced by the accuracy of the force calculation during orbit integration. We explore the accuracy of force calculations for two expansion methods and determine which one is preferable for orbit integration. Methods. We specifically compare two methods, one was introduced by Hernquist \& Ostriker (HO), which uses a spherical coordinate syste… ▽ More Aims. The orbital structure of galaxies is strongly influenced by the accuracy of the force calculation during orbit integration. We explore the accuracy of force calculations for two expansion methods and determine which one is preferable for orbit integration. Methods. We specifically compare two methods, one was introduced by Hernquist \& Ostriker (HO), which uses a spherical coordinate system and was built specifically for the Hernquist model, and the other by Vasiliev \& Athanassoula (CylSP) has a cylindrical coordinate system. Our comparisons include the Dehnen profile, its triaxial extension (of which the Hernquist profile is a special case) and a multicomponent system including a bar and disk density distributions for both analytical models and $N$-body realizations.We specifically compare two methods, one was introduced by Hernquist \& Ostriker (HO), which uses a spherical coordinate system and was built specifically for the Hernquist model, and the other by Vasiliev \& Athanassoula (CylSP) has a cylindrical coordinate system. Our comparisons include the Dehnen profile, its triaxial extension (of which the Hernquist profile is a special case) and a multicomponent system including a bar and disk density distributions for both analytical models and $N$-body realizations. Results. For the generalized Dehnen density, the CylSP method is more accurate than the HO method for nearly all inner power-law indices and shapes at all radii. For $N$-body realizations of the Dehnen model, or snapshots of an $N$-body simulation, the CylSP method is more accurate than the HO method in the central region for the oblate, prolate, and triaxial Hernquist profiles if the particle number is more than $5\times 10^5$. For snapshots of the Hernquist models with spherical shape, the HO method is preferred. (Abridged) △ Less

Submitted 25 May, 2020; v1 submitted 10 May, 2020; originally announced May 2020.

Comments: 14 pages, 17 figures, Accepted for publication in A&A

Journal ref: A&A 639, A38 (2020)

Suppressed or enhanced central star formation rates in late-type barred galaxies

Authors: Jing Wang, E. Athanassoula, Si-Yue Yu, Christian Wolf, Li Shao, Hua Gao, T. H. Randriamampandry

Abstract: Bars in disc-dominated galaxies are able to drive gas inflow inside the corotation radius, thus enhancing the central star formation rate (SFR). Previous work, however, has found that disc-dominated galaxies with centrally suppressed SFR frequently host a bar. Here we investigate possible causes for the suppression of central SFR in such cases. We compare physical properties of a sample of disc-do… ▽ More Bars in disc-dominated galaxies are able to drive gas inflow inside the corotation radius, thus enhancing the central star formation rate (SFR). Previous work, however, has found that disc-dominated galaxies with centrally suppressed SFR frequently host a bar. Here we investigate possible causes for the suppression of central SFR in such cases. We compare physical properties of a sample of disc-dominated barred galaxies with high central SFR (HC galaxies) with those of a sample of disc-dominated barred galaxies with low central SFR (LC galaxies). We find that the two samples have on average similar HI content and bars of similar strength. But we also find that the HCs have bluer colors than LCs, and that outside the bar region they host stronger spiral arms than the LCs where closed rings are more often seen. We discuss and evaluate the possible causes for the suppression of the central SFR in the LC galaxies as opposed to its enhancement in the HC galaxies. △ Less

Submitted 11 March, 2020; originally announced March 2020.

Comments: 22 pages, 11 figures, accepted for publication in ApJ

Tidally induced warps of spiral galaxies in IllustrisTNG

Authors: Marcin Semczuk, Ewa L. Lokas, Elena D'Onghia, E. Athanassoula, Victor P. Debattista, Lars Hernquist

Abstract: Warps are common features in both stellar and gaseous disks of nearby spiral galaxies with the latter usually easier to detect. Several theories have been proposed in the literature to explain their formation and prevalence, including tidal interactions with external galaxies. Observational correlations also suggest the importance of tides for warp formation. Here, we use the TNG100 run from the m… ▽ More Warps are common features in both stellar and gaseous disks of nearby spiral galaxies with the latter usually easier to detect. Several theories have been proposed in the literature to explain their formation and prevalence, including tidal interactions with external galaxies. Observational correlations also suggest the importance of tides for warp formation. Here, we use the TNG100 run from the magnetohydrodynamical cosmological simulation suite IllustrisTNG to investigate the connection between interactions and the formation of gas warps. We find that in the sample of well-resolved gas-rich spiral galaxies ($10^{10}\lesssim\mathrm{M_{*}/M_{\odot}}\lesssim10^{11}$ at $z=0$) from the simulation TNG100-1, about $16\%$ possess the characteristic S-shaped warp. Around one third of these objects have their vertical morphology induced by interactions with other galaxies. Half of these interactions end with the perturber absorbed by the host. Warps formed in interactions are more asymmetrical than the remaining sample, however after the interaction the asymmetry decreases with time. We find that warps induced by interactions survive on average for $<1$ Gyr. The angle between the orbital angular momentum of the perturber and the angular momentum of the host's disk that most likely leads to warp formation is around 45 degrees. While our main goal is to investigate tidally induced warps, we find that during interactions in addition to tides, new gas that is accreted from infalling satellites also can contribute to warp formation. △ Less

Submitted 17 February, 2020; originally announced February 2020.

Comments: 13 pages, 11 figures, submitted to MNRAS

The orbital content of bars: The origin of 'non-x1-tree', bar-supporting orbits

Authors: P. A. Patsis, E. Athanassoula

Abstract: Recently, many orbital studies in barred galaxy potentials have revealed the existence of orbits which are not trapped around x1 tree orbits, but could be potentially appropriate building blocks for bars. These findings question the uniqueness of the x1 family as the standard paradigm of orbital motion in galactic bars. The main goal of this paper is to investigate the role that such orbits could… ▽ More Recently, many orbital studies in barred galaxy potentials have revealed the existence of orbits which are not trapped around x1 tree orbits, but could be potentially appropriate building blocks for bars. These findings question the uniqueness of the x1 family as the standard paradigm of orbital motion in galactic bars. The main goal of this paper is to investigate the role that such orbits could play in shaping the morphology of bars. We trace the morphological patterns appearing in the face-on and edge-on views of the non-periodic orbits presented in these studies and we show that they are introduced in the system by second type bifurcations of x1. For this purpose, we use a typical 3D Ferrers bar model and follow the radial and vertical bifurcations of the x1 family considered as being multi periodic, with multi = 2, 3, 5. The variation of the stability indices of x1 in the multi = 2, 3 cases give us also the four- and six- periodic orbits, respectively. We tabulate these orbits including all information necessary to assess their role as appropriate building blocks. We discuss their stability and their extent, as well as their size and morphological evolution, as a function of energy. We conclude that even the most important of the multi-periodic orbits presented in Tables 2 to 5 are less appropriate building blocks for bars than the families of the x1 tree at the same energy. △ Less

Submitted 14 February, 2020; originally announced February 2020.

Comments: 20 pages, 16 figures, 4 tables with small figures, published in Monthly Notices of the Royal Astronomical Society

Journal ref: MNRAS 490, pp.2740-2759, 2020

The sequence of spiral arm classes: Observational signatures of persistent spiral density waves in grand-design galaxies

Authors: A. Bittner, D. A. Gadotti, B. G. Elmegreen, E. Athanassoula, D. M. Elmegreen, A. Bosma, J. Munoz-Mateos

Abstract: We investigate how the properties of spiral arms relate to other fundamental galaxy properties. To this end, we use previously published measurements of those properties, and our own measurements of arm-interarm luminosity contrasts for a large sample of galaxies, using 3.6$μ$m images from the Spitzer Survey of Stellar Structure in Galaxies. Flocculent galaxies are clearly distinguished from other… ▽ More We investigate how the properties of spiral arms relate to other fundamental galaxy properties. To this end, we use previously published measurements of those properties, and our own measurements of arm-interarm luminosity contrasts for a large sample of galaxies, using 3.6$μ$m images from the Spitzer Survey of Stellar Structure in Galaxies. Flocculent galaxies are clearly distinguished from other spiral arm classes, especially by their lower stellar mass and surface density. Multi-armed and grand-design galaxies are similar in most of their fundamental parameters, excluding some bar properties and the bulge-to-total luminosity ratio. Based on these results, we discuss dense, classical bulges as a necessary condition for standing spiral wave modes in grand-design galaxies. We further find a strong correlation between bulge-to-total ratio and bar contrast, and a weaker correlation between arm and bar contrasts. △ Less

Submitted 2 October, 2019; originally announced October 2019.

Comments: 4 pages, 3 figures, to appear in Galactic Dynamics in the Era of Large Surveys, M. Valluri & J. A. Sellwood (eds.)

Disc galaxies formed from major mergers in Illustris

Authors: Nicolas Peschken, Ewa L. Łokas, E. Athanassoula

Abstract: We show how wet major mergers can create disc galaxies in a cosmological context, using the Illustris simulation. We select a sample of 38 disc galaxies having experienced a major merger in their history with no subsequent significant minor merger, and study how they transform into discs after the merger. In agreement with what was previously found in controlled simulations of such mergers, we fin… ▽ More We show how wet major mergers can create disc galaxies in a cosmological context, using the Illustris simulation. We select a sample of 38 disc galaxies having experienced a major merger in their history with no subsequent significant minor merger, and study how they transform into discs after the merger. In agreement with what was previously found in controlled simulations of such mergers, we find that their disc is built gradually from young stars formed after the merger in the disc region, while the old stars born before the merger form an ellipsoidal component. Focusing on one fiducial case from our sample, we show how the gas was initially dispersed in the halo region right after the merger, but is then accreted onto a disc to form stars, and builds the disc component. We then select a sample of major mergers creating elliptical galaxies, to show that those cases correspond mainly to dry mergers, where the lack of star formation prevents the formation of a disc in the remnant galaxy. The amount of gas in the remnant galaxy after the merger is therefore essential to determine the final outcome of a major merger. △ Less

Submitted 3 March, 2020; v1 submitted 3 September, 2019; originally announced September 2019.

Comments: 14 pages, 14 figures, 1 table, accepted for publication in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 493, Issue 1, March 2020, Pages 1375-1387

Linking bar- and interaction-driven molecular gas concentration with centrally-enhanced star formation in EDGE-CALIFA galaxies

Authors: Ryan Chown, Cheng Li, E. Athanassoula, Niu Li, Christine D. Wilson, Lin Lin, Houjun Mo, Laura C. Parker, Ting Xiao

Abstract: We study the spatially resolved star formation history and molecular gas distribution of 58 nearby galaxies, using integral field spectroscopy from the CALIFA survey and CO $J=1\rightarrow 0$ intensity mapping from the CARMA EDGE survey. We use the 4000 Å break (D$_n$4000), the equivalent width of the H$δ$ absorption line (EW H$δ_A$), and the equivalent width of the H$α$ emission line (EW H$α$) to… ▽ More We study the spatially resolved star formation history and molecular gas distribution of 58 nearby galaxies, using integral field spectroscopy from the CALIFA survey and CO $J=1\rightarrow 0$ intensity mapping from the CARMA EDGE survey. We use the 4000 Å break (D$_n$4000), the equivalent width of the H$δ$ absorption line (EW H$δ_A$), and the equivalent width of the H$α$ emission line (EW H$α$) to measure the recent star formation history (SFH) of these galaxies. We measure radial profiles of the three SFH indicators and molecular gas mass surface density, from which we measure the level of centrally enhanced star formation and the molecular gas concentration. When we separate our galaxies into categories of barred (17 galaxies), unbarred (24 galaxies), and merging/paired (17 galaxies) we find that the galaxies which have centrally-enhanced star formation (19/58) are either barred (13/19) or in mergers/pairs (6/19) with relatively high molecular gas concentrations. A comparison between our barred galaxies and a snapshot of a hydrodynamic $N$-body simulation of a barred galaxy shows that the current theory of bar formation and evolution can qualitatively reproduce the main features of the observed galaxies in our sample, including both the sharp decrease of stellar age in the galactic center and the gradual decrease of age with increasing distance from center. These findings provide substantial evidence for a picture in which cold gas is transported inward by a bar or tidal interaction, which leads to the growth and rejuvenation of star formation in the central region. △ Less

Submitted 4 February, 2019; v1 submitted 19 October, 2018; originally announced October 2018.

Comments: Accepted by MNRAS

SDSS-IV MaNGA: pattern speeds of barred galaxies

Authors: Rui Guo, Shude Mao, E. Athanassoula, Hongyu Li, Junqiang Ge, R. J. Long, Michael Merrifield, Karen Masters

Abstract: The MaNGA project has obtained IFU data for several thousand nearby galaxies, including barred galaxies. With the two dimensional spectral and kinematic information provided by IFUs, we can measure the pattern speed of a barred galaxy, which determines the bar dynamics. We apply the non-parametric method proposed by Tremaine \& Weinberg to estimate the bar pattern speed for 53 barred galaxies, mak… ▽ More The MaNGA project has obtained IFU data for several thousand nearby galaxies, including barred galaxies. With the two dimensional spectral and kinematic information provided by IFUs, we can measure the pattern speed of a barred galaxy, which determines the bar dynamics. We apply the non-parametric method proposed by Tremaine \& Weinberg to estimate the bar pattern speed for 53 barred galaxies, making this the largest sample studied so far in this way. Our sample is selected from the MaNGA first public data release as part of SDSS Data Release 13 according mainly to the axis ratio and position angle difference between the bar and disc, while kinematic data is from the later SDSS Data Release 14. We have used both the photometric position angle from the photometric image and the kinematic position angle from the stellar velocity map to derive the pattern speed. Combining three independent bar length measurements and the circular velocity from Jeans Anisotropic modelling (JAM), we also determine the dimensionless ratio $\cal{R}$ of the corotation radius to the bar length. We find that the galaxy's position angle is the main uncertainty in determining the bar pattern speed. The kinematic position angle leads to fewer ultrafast bars than the photometric position angle, and this could be due to the method of measuring the kinematic position angle. We study the dependence of $\cal{R}$ values on galaxy properties such as the dark matter fraction from JAM modelling and the stellar age and metallicity from stellar population synthesis (SPS). A positive correlation between the bar length and bar strength is found: the longer the bar, the stronger the bar. However, no other significant correlations are found. This may result from errors in deriving the $\cal{R}$ values or from the complex formation and slowdown processes of galactic bars. △ Less

Submitted 7 October, 2018; originally announced October 2018.

Comments: 27 pages, 21 figures. Accepted for publication in MNRAS

Fundamental Physics with the Square Kilometre Array

Authors: A. Weltman, P. Bull, S. Camera, K. Kelley, H. Padmanabhan, J. Pritchard, A. Raccanelli, S. Riemer-Sørensen, L. Shao, S. Andrianomena, E. Athanassoula, D. Bacon, R. Barkana, G. Bertone, C. Bonvin, A. Bosma, M. Brüggen, C. Burigana, C. Bœhm, F. Calore, J. A. R. Cembranos, C. Clarkson, R. M. T. Connors, Á. de la Cruz-Dombriz, P. K. S. Dunsby , et al. (28 additional authors not shown)

Abstract: The Square Kilometre Array (SKA) is a planned large radio interferometer designed to operate over a wide range of frequencies, and with an order of magnitude greater sensitivity and survey speed than any current radio telescope. The SKA will address many important topics in astronomy, ranging from planet formation to distant galaxies. However, in this work, we consider the perspective of the SKA a… ▽ More The Square Kilometre Array (SKA) is a planned large radio interferometer designed to operate over a wide range of frequencies, and with an order of magnitude greater sensitivity and survey speed than any current radio telescope. The SKA will address many important topics in astronomy, ranging from planet formation to distant galaxies. However, in this work, we consider the perspective of the SKA as a facility for studying physics. We review four areas in which the SKA is expected to make major contributions to our understanding of fundamental physics: cosmic dawn and reionisation; gravity and gravitational radiation; cosmology and dark energy; and dark matter and astroparticle physics. These discussions demonstrate that the SKA will be a spectacular physics machine, which will provide many new breakthroughs and novel insights on matter, energy and spacetime. △ Less

Submitted 22 October, 2019; v1 submitted 5 October, 2018; originally announced October 2018.

Comments: 69 pages, PASA, in press

Journal ref: Publ. Astron. Soc. Aust. 37 (2020) e002

Unveiling the environment and faint features of the isolated galaxy CIG 96 with deep optical and HI observations

Authors: P. Ramirez-Moreta, L. Verdes-Montenegro, J. Blasco-Herrera, S. Leon, A. Venhola, M. Yun, V. Peris, R. Peletier, G. Verdoes Kleijn, E. Unda-Sanzana, D. Espada, A. Bosma, E. Athanassoula, M. Argudo-Fernandez, J. Sabater, J. C. Munoz-Mateos, M. G. Jones, W. Huchtmeier, J. E. Ruiz, J. Iglesias-Paramo, M. Fernandez-Lorenzo, J. Beckman, S. Sanchez-Exposito, J. Garrido

Abstract: Asymmetries in HI in galaxies are often caused by the interaction with close companions, making isolated galaxies an ideal framework to study secular evolution. The AMIGA project has demonstrated that isolated galaxies show the lowest level of asymmetry in their HI integrated profiles, yet some present significant asymmetries. CIG 96 (NGC 864) is a representative case reaching a 16% level. Our aim… ▽ More Asymmetries in HI in galaxies are often caused by the interaction with close companions, making isolated galaxies an ideal framework to study secular evolution. The AMIGA project has demonstrated that isolated galaxies show the lowest level of asymmetry in their HI integrated profiles, yet some present significant asymmetries. CIG 96 (NGC 864) is a representative case reaching a 16% level. Our aim is to investigate the HI asymmetries of this spiral galaxy and what processes have triggered the star-forming regions observed in the XUV pseudoring. We performed deep optical observations at CAHA 1.23m, 2.2m and VST telescopes. We reach surface brightness (SB) limits of mu_2.2m = 27.5 mag arcsec-2 (Cous R) and mu_VST = 28.7mag arcsec-2 (r) that show the XUV pseudoring of the galaxy in detail. Additionally, a wavelet filtering of the HI data cube from our deep observations with E/VLA telescope allowed us to reach a column density of N_HI = 8.9x10^18 cm -2 (5sigma) (28"x28" beam), lower than in any isolated galaxy. We confirm that the HI extends farther than 4xr_25 in all directions. Furthermore, we detect for the first time two gaseous structures (10^6 Msol) in the outskirts. The g-r colour index image from 1.23m shows extremely blue colours in certain regions of the pseudoring where N_HI>8.5x10^20 cm-2 , whereas the rest show red colours. Galactic cirrus contaminate the field, setting an unavoidable detection limit at 28.5mag arcsec-2 (r). We detect no stellar link within 1degx1deg or gaseous link within 40'x40' between CIG 96 and any companion. The isolation criteria rule out interactions with other similar-sized galaxies for at least 2.7Gyr. Using existing stellar evolution models, the age of the pseudoring is estimated at 1Gyr or older. Undetected previously accreted companions and cold gas accretion remain as the main hypothesis to explain the optical pseudoring and HI features of CIG 96. △ Less

Submitted 26 September, 2018; originally announced September 2018.

Comments: 23 pages, 18 figures, 4 tables

Journal ref: A&A 619, A163 (2018)

Tidally induced bars in gas-rich dwarf galaxies orbiting the Milky Way

Authors: Grzegorz Gajda, Ewa L. Łokas, E. Athanassoula

Abstract: Bars in galaxies may form not only through instability, but also as a result of an interaction with another galaxy. In particular, they may appear in disky dwarf galaxies interacting with Milky Way-like galaxies. Here we report the results of N-body/SPH simulations of such dwarfs orbiting in the static potential of a larger galaxy. We used several models of the dwarf galaxy, all of the same mass,… ▽ More Bars in galaxies may form not only through instability, but also as a result of an interaction with another galaxy. In particular, they may appear in disky dwarf galaxies interacting with Milky Way-like galaxies. Here we report the results of N-body/SPH simulations of such dwarfs orbiting in the static potential of a larger galaxy. We used several models of the dwarf galaxy, all of the same mass, but covering a large range of gas fractions: 0, 30 and 70%. We also tested the impact of subgrid star formation processes. In all cases bars of similar length formed in the stellar disk of the dwarfs at the first pericenter passage. However, unexpectedly, the gaseous component remained approximately axisymmetric and unaffected by the bar potential. The bar properties did not change significantly between two consecutive pericenters. The impact of the later encounters with the host depends strongly on the exact orientation of the bar at the pericenter. When the bar is spun up by the tidal force torque, it is also shortened. Conversely, if it is slowed down, it gets longer. In the models with a low gas fraction the bars were more pronounced and survived until the end of the simulations, while in the dwarfs with a high gas fraction the bars were destroyed after the second or third pericenter passage. In terms of the ratio of the corotation radius to the bar length, the bars are slow, and remain so independently of the encounters with the host. △ Less

Submitted 23 October, 2018; v1 submitted 2 July, 2018; originally announced July 2018.

Comments: 19 pages, 18 figures, accepted to ApJ

Journal ref: ApJ (2018), 868, 100

Tidally induced morphology of M33 in hydrodynamical simulations of its recent interaction with M31

Authors: Marcin Semczuk, Ewa L. Lokas, Jean-Baptiste Salomon, E. Athanassoula, Elena D'Onghia

Abstract: We present a hydrodynamical model of M33 and its recent interaction with M31. This scenario was previously proposed in the literature in order to explain the distorted gaseous and stellar disks of M33, as well as the increased star formation rate in both objects around 2 Gyr ago. We used an orbit integration scheme to find which estimate of the transverse velocity of M31 more favors the interactio… ▽ More We present a hydrodynamical model of M33 and its recent interaction with M31. This scenario was previously proposed in the literature in order to explain the distorted gaseous and stellar disks of M33, as well as the increased star formation rate in both objects around 2 Gyr ago. We used an orbit integration scheme to find which estimate of the transverse velocity of M31 more favors the interaction scenario and then tried to reproduce it in our simulations. M33 was modeled as a stellar and gaseous disk embedded in a live dark-matter halo, while M31 was approximated only with a live dark halo. In the simulations the two galaxies passed each other with a pericenter distance of 37 kpc. Tides excited a two-armed spiral structure in the M33 disk, which is found to be the predominant spiral signal in the observed galaxy and has been long known as a feature easily induced by tidal interactions. We found that the gaseous warp produced by the interaction did not resemble enough the observed one and we performed an additional simulation including the hot gas halo of M31 to show that this feature can be properly reproduced by tidal forces and ram pressure stripping acting simultaneously on the gaseous disk. In addition to the spiral arms, tidal forces produced the stellar stream similar to the observed one and triggered a star formation burst at similar radii as it is observed. △ Less

Submitted 18 July, 2018; v1 submitted 12 April, 2018; originally announced April 2018.

Comments: 22 pages, 17 figures, accepted for publication in ApJ

Using N-body simulations to understand the chemo-dynamical evolution of the inner Milky Way

Authors: E. Athanassoula

Abstract: I present examples of how chemo-dynamical N-body simulations can help understanding the structure and evolution of the inner Galaxy. Such simulations reproduce the observed links between kinematics, morphology and chemistry in the bar/bulge region and explain them by the self-consistent cohabitation of a number of components. Galactic archaeology, applied to simulation snapshots, explains the sequ… ▽ More I present examples of how chemo-dynamical N-body simulations can help understanding the structure and evolution of the inner Galaxy. Such simulations reproduce the observed links between kinematics, morphology and chemistry in the bar/bulge region and explain them by the self-consistent cohabitation of a number of components. Galactic archaeology, applied to simulation snapshots, explains the sequence in which the stars of the various components were formed. The thick disc stars form earlier than those of the thin disc and in a much shorter time scale. The bar in the thick disc is horizontally thicker than that of the thin disc and has a different vertical morphology. The Galaxy's inner disc scalelength is much smaller than what is expected from nearby galaxies of similar stellar mass. △ Less

Submitted 23 January, 2018; originally announced January 2018.

Comments: 8 pages, 3 figures, to be published in the Proceedings of IAU Symposium 334,"Rediscovering our Galaxy", Eds. C. Chiappini, I. Minchev, E. Starkenburg & M. Valentini, Cambridge University Press, Cambridge

French SKA White Book - The French Community towards the Square Kilometre Array

Authors: F. Acero, J. -T. Acquaviva, R. Adam, N. Aghanim, M. Allen, M. Alves, R. Ammanouil, R. Ansari, A. Araudo, E. Armengaud, B. Ascaso, E. Athanassoula, D. Aubert, S. Babak, A. Bacmann, A. Banday, K. Barriere, F. Bellossi, J. -P. Bernard, M. G. Bernardini, M. Béthermin, E. Blanc, L. Blanchet, J. Bobin, S. Boissier , et al. (153 additional authors not shown)

Abstract: The "Square Kilometre Array" (SKA) is a large international radio telescope project characterised, as suggested by its name, by a total collecting area of approximately one square kilometre, and consisting of several interferometric arrays to observe at metric and centimetric wavelengths. The deployment of the SKA will take place in two sites, in South Africa and Australia, and in two successive p… ▽ More The "Square Kilometre Array" (SKA) is a large international radio telescope project characterised, as suggested by its name, by a total collecting area of approximately one square kilometre, and consisting of several interferometric arrays to observe at metric and centimetric wavelengths. The deployment of the SKA will take place in two sites, in South Africa and Australia, and in two successive phases. From its Phase 1, the SKA will be one of the most formidable scientific machines ever deployed by mankind, and by far the most impressive in terms of data throughput and required computing power. With the participation of almost 200 authors from forty research institutes and six private companies, the publication of this French SKA white paper illustrates the strong involvement in the SKA project of the French astronomical community and of a rapidly growing number of major scientific and technological players in the fields of Big Data, high performance computing, energy production and storage, as well as system integration. △ Less

Submitted 28 March, 2018; v1 submitted 19 December, 2017; originally announced December 2017.

Comments: Editor in chief: C. Ferrari; Editors: M. Alves, S. Bosse, S. Corbel, A. Ferrari, K. Ferrière, S. Gauffre, E. Josselin, G. Lagache, S. Lambert, G. Marquette, J.-M. Martin, M.-A. Miville-Deschênes, L. Montier, B. Semelin, G. Theureau, S. Vergani, N. Vilmer, P. Zarka; Original file with high resolution figures at SKA-France link: https://ska-france.oca.eu/images/SKA-France-Media/FWB_051017.pdf

An Overview of the MHONGOOSE Survey: Observing Nearby Galaxies with MeerKAT

Authors: W. J. G. de Blok, E. A. K. Adams, P. Amram, E. Athanassoula, I. Bagetakos, C. Balkowski, M. A. Bershady, R. Beswick, F. Bigiel, S. -L. Blyth, A. Bosma, R. S. Booth, A. Bouchard, E. Brinks, C. Carignan, L. Chemin, F. Combes, J. Conway, E. C. Elson, J. English, B. Epinat, B. S. Frank, J. Fiege, F. Fraternali, J. S. Gallagher , et al. (43 additional authors not shown)

Abstract: MHONGOOSE is a deep survey of the neutral hydrogen distribution in a representative sample of 30 nearby disk and dwarf galaxies with HI masses from 10^6 to ~10^{11} M_sun, and luminosities from M_R ~ -12 to M_R ~ -22. The sample is selected to uniformly cover the available range in log(M_HI). Our extremely deep observations, down to HI column density limits of well below 10^{18} cm^{-2} - or a few… ▽ More MHONGOOSE is a deep survey of the neutral hydrogen distribution in a representative sample of 30 nearby disk and dwarf galaxies with HI masses from 10^6 to ~10^{11} M_sun, and luminosities from M_R ~ -12 to M_R ~ -22. The sample is selected to uniformly cover the available range in log(M_HI). Our extremely deep observations, down to HI column density limits of well below 10^{18} cm^{-2} - or a few hundred times fainter than the typical HI disks in galaxies - will directly detect the effects of cold accretion from the intergalactic medium and the links with the cosmic web. These observations will be the first ever to probe the very low-column density neutral gas in galaxies at these high resolutions. Combination with data at other wavelengths, most of it already available, will enable accurate modelling of the properties and evolution of the mass components in these galaxies and link these with the effects of environment, dark matter distribution, and other fundamental properties such as halo mass and angular momentum. MHONGOOSE can already start addressing some of the SKA-1 science goals and will provide a comprehensive inventory of the processes driving the transformation and evolution of galaxies in the nearby universe at high resolution and over 5 orders of magnitude in column density. It will be a Nearby Galaxies Legacy Survey that will be unsurpassed until the advent of the SKA, and can serve as a highly visible, lasting statement of MeerKAT's capabilities. △ Less

Submitted 25 September, 2017; originally announced September 2017.

Comments: To be published in "MeerKAT Science: On the Pathway to the SKA". Proceedings of Science. Workshop held 25-27 May, 2016 Stellenbosch, South Africa

Morphology and enhanced star formation in a Cartwheel-like ring galaxy

Authors: F. Renaud, E. Athanassoula, P. Amram, A. Bosma, F. Bournaud, P. -A. Duc, B. Epinat, J. Fensch, K. Kraljic, V. Perret, C. Struck

Abstract: We use hydrodynamical simulations of a Cartwheel-like ring galaxy, modelled as a nearly head-on collision of a small companion with a larger disc galaxy, to probe the evolution of the gaseous structures and flows, and to explore the physical conditions setting the star formation activity. Star formation is first quenched by tides as the companion approaches, before being enhanced shortly after the… ▽ More We use hydrodynamical simulations of a Cartwheel-like ring galaxy, modelled as a nearly head-on collision of a small companion with a larger disc galaxy, to probe the evolution of the gaseous structures and flows, and to explore the physical conditions setting the star formation activity. Star formation is first quenched by tides as the companion approaches, before being enhanced shortly after the collision. The ring ploughs the disc material as it radially extends, and almost simultaneously depletes its stellar and gaseous reservoir into the central region, through the spokes, and finally dissolve 200 Myr after the collision. Most of star formation first occurs in the ring before this activity is transferred to the spokes and then the nucleus. We thus propose that the location of star formation traces the dynamical stage of ring galaxies, and could help constrain their star formation histories. The ring hosts tidal compression associated with strong turbulence. This compression yields an azimuthal asymmetry, with maxima reached in the side furthest away from the nucleus, which matches the star formation activity distribution in our models and in observed ring systems. The interaction triggers the formation of star clusters significantly more massive than before the collision, but less numerous than in more classical galaxy interactions. The peculiar geometry of Cartwheel-like objects thus yields a star (cluster) formation activity comparable to other interacting objects, but with notable second order differences in the nature of turbulence, the enhancement of the star formation rate, and the number of massive clusters formed. △ Less

Submitted 8 September, 2017; originally announced September 2017.

Comments: Accepted for publication in MNRAS

How do spiral arm contrasts relate to bars, disk breaks and other fundamental galaxy properties?

Authors: Adrian Bittner, Dimitri A. Gadotti, Bruce G. Elmegreen, E. Athanassoula, Debra M. Elmegreen, Albert Bosma, Juan-Carlos Muñoz-Mateos

Abstract: We investigate how the properties of spiral arms relate to other fundamental galaxy properties, including bars and disc breaks. We use previously published measurements of those properties, and our own measurements of arm and bar contrasts for a large sample of galaxies, using $3.6 \mathrm{μm}$ images from the Spitzer Survey of Stellar Structure in Galaxies (S4G). Flocculent galaxies are clearly d… ▽ More We investigate how the properties of spiral arms relate to other fundamental galaxy properties, including bars and disc breaks. We use previously published measurements of those properties, and our own measurements of arm and bar contrasts for a large sample of galaxies, using $3.6 \mathrm{μm}$ images from the Spitzer Survey of Stellar Structure in Galaxies (S4G). Flocculent galaxies are clearly distinguished from other spiral arm classes, especially by their lower stellar mass and surface density. Multi-armed and grand-design galaxies are similar in most of their fundamental parameters, excluding some bar properties and the bulge-to-total ratio. Based on these results, we revisit the sequence of spiral arm classes, and discuss classical bulges as a necessary condition for standing spiral wave modes in grand-design galaxies. We find a strong correlation between bulge-to-total ratio and bar contrast, and a weaker correlation between arm and bar contrasts. Barred and unbarred galaxies exhibit similar arm contrasts, but the highest arm contrasts are found exclusively in barred galaxies. Interestingly, the bar contrast, and its increase from flocculent to grand-design galaxies, is systematically more significant than that of the arm contrast. We corroborate previous findings concerning a connection between bars and disc breaks. In particular, in grand-design galaxies the bar contrast correlates with the normalised disc break radius. This does not hold for other spiral arm classes or the arm contrast. Our measurements of arm and bar contrast and radial contrast profiles are publicly available. △ Less

Submitted 29 June, 2017; originally announced June 2017.

Comments: Accepted for publication in MNRAS; 20 pages, 17 figures, 4 tables

Journal ref: MNRAS 471, 1070-1087 (2017)

Tidally induced bars in dwarf galaxies on different orbits around a Milky Way-like host

Authors: Grzegorz Gajda, Ewa L. Lokas, E. Athanassoula

Abstract: Bars in galaxies may develop through a global instability or due to an interaction with another system. We study bar formation in disky dwarf galaxies orbiting a Milky Way-like galaxy. We employ $N$-body simulations to study the impact of initial orbital parameters: the size of the dwarf galaxy orbit and the inclination of its disc with respect to the orbital plane. In all cases a bar develops in… ▽ More Bars in galaxies may develop through a global instability or due to an interaction with another system. We study bar formation in disky dwarf galaxies orbiting a Milky Way-like galaxy. We employ $N$-body simulations to study the impact of initial orbital parameters: the size of the dwarf galaxy orbit and the inclination of its disc with respect to the orbital plane. In all cases a bar develops in the center of the dwarf during the first pericenter on its orbit around the host. Between subsequent pericenter passages the bars are stable, but at the pericenters they are usually weakened and shortened. The initial properties and details of the further evolution of the bars depend heavily on the orbital configuration. We find that for the exactly prograde orientation, the strongest bar is formed for the intermediate-size orbit. On the tighter orbit, the disc is too disturbed and stripped to form a strong bar. On the wider orbit, the tidal interaction is too weak. The dependence on the disc inclination is such that weaker bars form in more inclined discs. The bars experience either a very weak buckling or none at all. We do not observe any secular evolution, possibly because the dwarfs are perturbed at each pericenter passage. The rotation speed of the bars can be classified as slow ($R_\mathrm{CR}/l_\mathrm{bar}\sim2-3$). We attribute this to the loss of a significant fraction of the disc's rotation during the encounter with the host galaxy. △ Less

Submitted 22 May, 2017; v1 submitted 8 March, 2017; originally announced March 2017.

Comments: 17 pages, 14 figures, accepted to ApJ

Forming disc galaxies in major mergers: III. The effect of angular momentum on the radial density profiles of disc galaxies

Authors: N. Peschken, E. Athanassoula, S. A. Rodionov

Abstract: We study the effect of angular momentum on the surface density profiles of disc galaxies, using high resolution simulations of major mergers whose remnants have downbending radial density profiles (type II). As described in the previous papers of this series, in this scenario, most of the disc mass is acquired after the collision via accretion from a hot gaseous halo. We find that the inner and ou… ▽ More We study the effect of angular momentum on the surface density profiles of disc galaxies, using high resolution simulations of major mergers whose remnants have downbending radial density profiles (type II). As described in the previous papers of this series, in this scenario, most of the disc mass is acquired after the collision via accretion from a hot gaseous halo. We find that the inner and outer disc scalelengths, as well as the break radius, correlate with the total angular momentum of the initial merging system, and are larger for high angular momentum systems. We follow the angular momentum redistribution in our simulated galaxies, and find that, like the mass, the disc angular momentum is acquired via accretion, i.e. to the detriment of the gaseous halo. Furthermore, high angular momentum systems give more angular momentum to their discs, which affects directly their radial density profile. Adding simulations of isolated galaxies to our sample, we find that the correlations are valid also for disc galaxies evolved in isolation. We show that the outer part of the disc at the end of the simulation is populated mainly by inside-out stellar migration, and that in galaxies with higher angular momentum, stars travel radially further out. This, however, does not mean that outer disc stars (in type II discs) were mostly born in the inner disc. Indeed, generally the break radius increases over time, and not taking this into account leads to overestimating the number of stars born in the inner disc. △ Less

Submitted 2 March, 2017; originally announced March 2017.

Comments: 12 pages, 13 figures, accepted for publication in MNRAS

Forming disc galaxies in major mergers II. The central mass concentration problem and a comparison of GADGET3 with GIZMO

Authors: S. A. Rodionov, E. Athanassoula, N. Peschken

Abstract: Context: In a series of papers, we study the major merger of two disk galaxies in order to establish whether or not such a merger can produce a disc galaxy. Aims: Our aim here is to describe in detail the technical aspects of our numerical experiments. Methods: We discuss the initial conditions of our major merger, which consist of two protogalaxies on a collision orbit. We show that such merger s… ▽ More Context: In a series of papers, we study the major merger of two disk galaxies in order to establish whether or not such a merger can produce a disc galaxy. Aims: Our aim here is to describe in detail the technical aspects of our numerical experiments. Methods: We discuss the initial conditions of our major merger, which consist of two protogalaxies on a collision orbit. We show that such merger simulations can produce a non-realistic central mass concentration, and we propose simple, parametric, AGN-like feedback as a solution to this problem. Our AGN-like feedback algorithm is very simple: at each time-step we take all particles whose local volume density is above a given threshold value and increase their temperature to a preset value. We also compare the GADGET3 and GIZMO codes, by applying both of them to the same initial conditions. Results: We show that the evolution of isolated protogalaxies resembles the evolution of disk galaxies, thus arguing that our protogalaxies are well suited for our merger simulations. We demonstrate that the problem with the unphysical central mass concentration in our merger simulations is further aggravated when we increase the resolution. We show that our AGN-like feedback removes this non-physical central mass concentration, and thus allows the formation of realistic bars. Note that our AGN-like feedback mainly affects the central region of a model, without significantly modifying the rest of the galaxy. We demonstrate that, in the context of our kind of simulation, GADGET3 gives results which are very similar to those obtained with the PSPH (density independent SPH) flavor of GIZMO. Moreover, in the examples we tried, the differences between the results of the two flavors of GIZMO, namely PSPH, and MFM (mesh-less algorithm) are similar to and, in some comparisons, larger than the differences between the results of GADGET3 and PSPH. △ Less

Submitted 10 January, 2017; originally announced January 2017.

Comments: 13 pages, 11 figures, accepted for publication in A&A

Journal ref: A&A 600, A25 (2017)

Dynamical modelling of galactic disc outskirts

Authors: E. Athanassoula

Abstract: I review briefly some dynamical models of structures in the outer parts of disc galaxies, including models of polar rings, tidal tails and bridges. I then discuss the density distribution in the outer parts of discs. For this, I compare observations to results of a model in which the disc galaxy is in fact the remnant of a major merger, and find good agreement. This comparison includes radial prof… ▽ More I review briefly some dynamical models of structures in the outer parts of disc galaxies, including models of polar rings, tidal tails and bridges. I then discuss the density distribution in the outer parts of discs. For this, I compare observations to results of a model in which the disc galaxy is in fact the remnant of a major merger, and find good agreement. This comparison includes radial profiles of the projected surface density and of stellar age, as well as time evolution of the break radius and of the inner and outer disc scale lengths. I also compare the radial projected surface density profiles of dynamically motivated mono-age populations and find that, compared to older populations, younger ones have flatter density profiles in the inner region and steeper in the outer one. The break radius, however, does not vary with stellar age, again in good agreement with observations. △ Less

Submitted 19 December, 2016; originally announced December 2016.

Comments: 8 pages, 2 figures, invited review at IAU Symposium 321 "Formation and Evolution of Galaxy Outskirts", eds. A. Gil de Paz, J.C. Lee & J.H. Knapen

Metallicity-dependendent kinematics and morphology of the Milky Way bulge

Authors: E. Athanassoula, S. A. Rodionov, N. Prantzos

Abstract: We use N-body chemo-dynamic simulations to study the coupling between morphology, kinematics and metallicity of the bar/bulge region of our Galaxy. We make qualitative comparisons of our results with available observations and find very good agreement. We conclude that this region is complex, since it comprises several stellar components with different properties -- i.e. a boxy/peanut bulge, thin… ▽ More We use N-body chemo-dynamic simulations to study the coupling between morphology, kinematics and metallicity of the bar/bulge region of our Galaxy. We make qualitative comparisons of our results with available observations and find very good agreement. We conclude that this region is complex, since it comprises several stellar components with different properties -- i.e. a boxy/peanut bulge, thin and thick disc components, and, to lesser extents, a disky pseudobulge, a stellar halo and a small classical bulge -- all cohabiting in dynamical equilibrium. Our models show strong links between kinematics and metallicity, or morphology and metallicity, as already suggested by a number of recent observations. We discuss and explain these links. △ Less

Submitted 16 December, 2016; originally announced December 2016.

Comments: 5 pages, 4 figures, accepted for publication in MNRAS Letters

Constraining the dark matter content of NGC 1291 using hydrodynamic gas response simulations

Authors: F. Fragkoudi, E. Athanassoula, A. Bosma

Abstract: We present a pilot study on the nearby massive barred galaxy NGC 1291, in which we use dynamical modelling to constrain the disc mass-to-light ratio (M/L), thus breaking the degeneracy between the baryonic and dark matter in its central regions. We use the gas, specifically the morphology of the dust lanes on the leading side of the bar, as a tracer of the underlying gravitational potential. We ru… ▽ More We present a pilot study on the nearby massive barred galaxy NGC 1291, in which we use dynamical modelling to constrain the disc mass-to-light ratio (M/L), thus breaking the degeneracy between the baryonic and dark matter in its central regions. We use the gas, specifically the morphology of the dust lanes on the leading side of the bar, as a tracer of the underlying gravitational potential. We run a large number of hydrodynamic gas response simulations, in potentials obtained directly from near-infrared images of the galaxy, which have three free parameters: the M/L, the bar pattern speed and the height function. We explore the three-dimensional parameter space, by comparing the morphology of the shocks created in the gas response simulations with those of the observed dust lanes, and find the best-fitting models; these suggest that the M/L of NGC 1291 agrees with that predicted by stellar population synthesis models in the near-infrared ($\approx$0.6\,$M_{\odot}/L_{\odot}$), which leads to a borderline maximum disc for this galaxy. Furthermore, we find that the bar rotates fast, with a corotation radius which is $\leq$ 1.4 times the bar length. △ Less

Submitted 28 November, 2016; originally announced November 2016.

Comments: 16 pages including Appendix, 14 figures. Accepted for publication in MNRAS

Morpho-kinematics of z~1 galaxies probe the hierarchical scenario

Authors: M. Rodrigues, F. Hammer, H. Flores, M. Puech, E. Athanassoula

Abstract: We have studied a representative sample of intermediate-mass galaxies at z~1, observed by the kinematic survey KMOS3D. We have re-estimated the kinematical parameters from the published kinematic maps and analysed photometric data from HST to measure optical disk inclinations and PAs. We find that only half of the z~1 galaxies show kinematic properties consistent with rotating disks, using the sam… ▽ More We have studied a representative sample of intermediate-mass galaxies at z~1, observed by the kinematic survey KMOS3D. We have re-estimated the kinematical parameters from the published kinematic maps and analysed photometric data from HST to measure optical disk inclinations and PAs. We find that only half of the z~1 galaxies show kinematic properties consistent with rotating disks, using the same classification scheme than that adopted by the KMOS3D team. Because merger orbital motions can also brought rotation, we have also analysed galaxy morphologies from the available HST imagery. Combining these results to those from kinematics, it leads to a full morpho-kinematic classification. To test the robustness of the latter for disentangling isolated disks from mergers, we confronted the results with an analysis of pairs from the open-grism redshift survey 3D-HST. All galaxies found in pairs are affected by either kinematic and/or morphological perturbations. Conversely, all galaxies classified as virialized spirals are found to be isolated. A significant fraction (one fourth) of rotating disks classified from kinematics by the KMOS3D team are found in pairs, which further supports the need for a morpho-kinematic classification. It results that only one third of z~1 galaxies are isolated and virialized spirals, while 58% of them are likely involved in a merger sequence, from first approach to disk rebuilding. The later fraction is in good agreement with the results of semi-empirical ΛCDM models, supporting a merger-dominated hierarchical scenario as being the main driver of galaxy formation at least during the last 8 billion years. △ Less

Submitted 10 November, 2016; originally announced November 2016.

Orbital classification in an N-body bar

Authors: Yougang Wang, E. Athanassoula, Shude Mao

Abstract: The dynamics and evolution of any galactic structure are strongly influenced by the properties of the orbits that constitute it. In this paper, we compare two orbit classification schemes, one by Laskar (NAFF) , and the other by Carpintero and Aguilar (CA), by applying both of them to orbits obtained by following individual particles in a numerical simulation of a barred galaxy. We find that, at l… ▽ More The dynamics and evolution of any galactic structure are strongly influenced by the properties of the orbits that constitute it. In this paper, we compare two orbit classification schemes, one by Laskar (NAFF) , and the other by Carpintero and Aguilar (CA), by applying both of them to orbits obtained by following individual particles in a numerical simulation of a barred galaxy. We find that, at least for our case and some provisos, the main frequencies calculated by the two methods are in good agreement: for $80\%$ of the orbits the difference between the results of the two methods is less than $5\%$ for all three main frequencies. However, it is difficult to evaluate the amount of regular or chaotic bar orbits in a given system. The fraction of regular orbits obtained by the NAFF method strongly depends on the critical frequency drift parameter, while in the CA method the number of fundamental frequencies strongly depends on the frequency difference parameter $L_{\rm r}$ and the maximum integer used for searching the linear independence of the fundamental frequencies. We also find that, for a given particle, in general the projection of its motion along the bar minor axis is more regular than the other two projections, while the projection along the intermediate axis is the least regular. △ Less

Submitted 8 September, 2016; originally announced September 2016.

Comments: 14 pages, 21 figures; Accepted for publication in MNRAS

Evidence of Bar-induced Secular Evolution in the Inner Regions of Stellar Discs in Galaxies: What Shapes Disc Galaxies?

Authors: Taehyun Kim, Dimitri A. Gadotti, E. Athanassoula, Albert Bosma, Kartik Sheth, Myung Gyoon Lee

Abstract: We present evidence of bar-induced secular evolution in galactic discs using 3.6 ${μm}$ images of nearby galaxies from the Spitzer Survey of Stellar Structure in Galaxies (S4G). We find that among massive galaxies ($M_{star}/M_{sun}> 10^{10}$), longer bars tend to reside in inner discs having a flatter radial profile. Such galaxies show a light deficit in the disc surrounding the bar, within the b… ▽ More We present evidence of bar-induced secular evolution in galactic discs using 3.6 ${μm}$ images of nearby galaxies from the Spitzer Survey of Stellar Structure in Galaxies (S4G). We find that among massive galaxies ($M_{star}/M_{sun}> 10^{10}$), longer bars tend to reside in inner discs having a flatter radial profile. Such galaxies show a light deficit in the disc surrounding the bar, within the bar radius and often show a $Θ$-shaped morphology. We quantify this deficit and find that among all galaxies explored in this study (with $10^{9}<M_{star}/M_{sun}< 10^{11}$), galaxies with a stronger bar (i.e. longer and/or with a higher Bar/T) show a more pronounced deficit. We also examine simulation snapshots to confirm and extend results by Athanassoula and Misiriotis, showing that as bars evolve they become longer, while the light deficit in the disc becomes more pronounced. Theoretical studies have predicted that, as a barred galaxy evolves, the bar captures disc stars in its immediate neighbourhood so as to make the bar longer, stronger and thinner. Hence, we claim that the light deficit in the inner disc is produced by bars, which thus take part in shaping the mass distribution of their host galaxies. △ Less

Submitted 27 July, 2016; originally announced July 2016.

Comments: Accepted for publication in MNRAS, 12 pages, 8 figures

A close look at secular evolution: Boxy/peanut bulges reduce gas inflow to the central kiloparsec

Authors: F. Fragkoudi, E. Athanassoula, A. Bosma

Abstract: In this letter we investigate the effect of boxy/peanut (b/p) bulges on bar-induced gas inflow to the central kiloparsec, which plays a crucial role on the evolution of disc galaxies. We carry out hydrodynamic gas response simulations in realistic barred galaxy potentials, including or not the geometry of a b/p bulge, to investigate the amount of gas inflow induced in the different models. We find… ▽ More In this letter we investigate the effect of boxy/peanut (b/p) bulges on bar-induced gas inflow to the central kiloparsec, which plays a crucial role on the evolution of disc galaxies. We carry out hydrodynamic gas response simulations in realistic barred galaxy potentials, including or not the geometry of a b/p bulge, to investigate the amount of gas inflow induced in the different models. We find that b/p bulges can reduce the gas inflow rate to the central kiloparsec by more than an order of magnitude, which leads to a reduction in the amount of gas available in the central regions. We also investigate the effect of the dark matter halo concentration on these results, and find that for maximal discs, the effect of b/p bulges on gas inflow remains significant. The reduced amount of gas reaching the central regions due to the presence of b/p bulges could have significant repercussions on the formation of discy- (pseudo-) bulges, on the amount of nuclear star formation and feedback, on the fuel reservoir for AGN activity, and on the overall secular evolution of galaxies. △ Less

Submitted 14 June, 2016; originally announced June 2016.

Comments: Accepted for publication in Monthly Notices of the Royal Astronomical Society Letters. 5 pages, 6 figures

The orbital structure of a tidally induced bar

Authors: Grzegorz Gajda, Ewa L. Lokas, E. Athanassoula

Abstract: Orbits are the key building blocks of any density distribution and their study helps us understand the kinematical structure and the evolution of galaxies. Here we investigate orbits in a tidally induced bar of a dwarf galaxy, using an $N$-body simulation of an initially disky dwarf galaxy orbiting a Milky Way-like host. After the first pericenter passage, a tidally induced bar forms in the stella… ▽ More Orbits are the key building blocks of any density distribution and their study helps us understand the kinematical structure and the evolution of galaxies. Here we investigate orbits in a tidally induced bar of a dwarf galaxy, using an $N$-body simulation of an initially disky dwarf galaxy orbiting a Milky Way-like host. After the first pericenter passage, a tidally induced bar forms in the stellar component of the dwarf. The bar evolution is different than in isolated galaxies and our analysis focuses on the period before it buckles. We study the orbits in terms of their dominant frequencies, which we calculate in a Cartesian coordinate frame rotating with the bar. Apart from the well-known x$_1$ orbits we find many other types, mostly with boxy shapes of various degree of elongation. Some of them are also near-periodic, admitting frequency ratios of 4/3, 3/2 and 5/3. The box orbits have various degrees of vertical thickness but only a relatively small fraction of those have banana (i.e. smile/frown) or infinity-symbol shapes in the edge-on view. In the very center we also find orbits known from the potential of triaxial ellipsoids. The elongation of the orbits grows with distance from the center of the bar in agreement with the variation of the shape of the density distribution. Our classification of orbits leads to the conclusion that more than $80 \%$ of them have boxy shapes, while only $8 \%$ have shapes of classical x$_1$ orbits. △ Less

Submitted 17 August, 2016; v1 submitted 1 June, 2016; originally announced June 2016.

Comments: 15 pages, 15 figures, accepted for publication in ApJ

Journal ref: ApJ, 830, 108 (2016)

Forming first-ranked early-type galaxies through hierarchical dissipationless merging

Authors: IDILICO Collaboration, J. M. Solanes, J. D. Perea, L. Darriba, C. García-Gómez, A. Bosma, E. Athanassoula

Abstract: We have developed a computationally competitive N-body model of a previrialized aggregation of galaxies in a flat LambdaCDM universe to assess the role of the multiple mergers that take place during the formation stage of such systems in the configuration of the remnants assembled at their centres. An analysis of a suite of 48 simulations of low-mass forming groups (of about 1E13 solar masses) dem… ▽ More We have developed a computationally competitive N-body model of a previrialized aggregation of galaxies in a flat LambdaCDM universe to assess the role of the multiple mergers that take place during the formation stage of such systems in the configuration of the remnants assembled at their centres. An analysis of a suite of 48 simulations of low-mass forming groups (of about 1E13 solar masses) demonstrates that the gravitational dynamics involved in their hierarchical collapse is capable of creating realistic first-ranked galaxies without the aid of dissipative processes. Our simulations indicate that the brightest group galaxies (BGGs) constitute a distinct population from other group members, sketching a scenario in which the assembly path of these objects is dictated largely by the formation of their host system. We detect significant differences in the distribution of Sersic indices and total magnitudes, as well as a luminosity gap between BGGs and the next brightest galaxy that is positively correlated with the total luminosity of the parent group. Such gaps arise from both the grow of BGGs at the expense of lesser companions and the decrease in the relevance of second-ranked objects in equal measure. This results in a dearth of intermediate-mass galaxies which explains the characteristic central dip detected in their luminosity functions in dynamically young galaxy aggregations. The fact that the basic global properties of our BGGs define a thin mass fundamental plane strikingly similar to that followed giant early-type galaxies in the local universe reinforces confidence in the results obtained. △ Less

Submitted 26 May, 2016; originally announced May 2016.

Comments: 25 pages, 14 figures, 3 tables. Accepted to MNRAS

Breaking the disc-halo degeneracy in NGC 1291 using hydrodynamic simulations

Authors: F. Fragkoudi, E. Athanassoula, A. Bosma

Abstract: We present a pilot study on the nearby massive galaxy NGC 1291, in which we aim to constrain the dark matter in the inner regions, by obtaining a dynamical determination of the disc mass-to-light ratio (M/L). To this aim, we model the bar-induced dust lanes in the galaxy, using hydrodynamic gas response simulations. The models have three free parameters, the M/L of the disc, the bar pattern speed… ▽ More We present a pilot study on the nearby massive galaxy NGC 1291, in which we aim to constrain the dark matter in the inner regions, by obtaining a dynamical determination of the disc mass-to-light ratio (M/L). To this aim, we model the bar-induced dust lanes in the galaxy, using hydrodynamic gas response simulations. The models have three free parameters, the M/L of the disc, the bar pattern speed and the disc height function. We explore the parameter space to find the best fit models, i.e. those in which the morphology of the shocks in the gas simulations matches the observed dust lanes. The best-fit models suggest that the M/L of NGC 1291 agrees with that predicted by stellar population synthesis models in the near-infrared ($\approx$0.6\,$M_{\odot}/L_{\odot}$), which leads to a borderline maximum disc for this galaxy. The bar rotates fast, with corotation radius $\leq$ 1.4 times the bar length. Additionally, we find that the height function has a significant effect on the results, and can bias them towards lower or higher M/L. △ Less

Submitted 18 May, 2016; originally announced May 2016.

Comments: 3 pages. Proceedings of IAU Symposium 321, "Formation and evolution of galaxy outskirts", Eds. A. Gil de Paz, J. C. Lee & J. H. Knapen, Cambridge University Press, Cambridge

Tidally Induced Offset Disks in Magellanic Spiral Galaxies

Authors: Stephen A. Pardy, Elena D'Onghia, E. Athanassoula, Eric M. Wilcots, Kartik Sheth

Abstract: Magellanic spiral galaxies are a class of one-armed systems that often exhibit an offset stellar bar, and are rarely found around massive spiral galaxies. Using a set of N-body and hydrodynamic simulations we consider a dwarf-dwarf galaxy interaction as the driving mechanism for the formation of this peculiar class of systems. We investigate here the relation between the dynamical, stellar and gas… ▽ More Magellanic spiral galaxies are a class of one-armed systems that often exhibit an offset stellar bar, and are rarely found around massive spiral galaxies. Using a set of N-body and hydrodynamic simulations we consider a dwarf-dwarf galaxy interaction as the driving mechanism for the formation of this peculiar class of systems. We investigate here the relation between the dynamical, stellar and gaseous disk center and the bar. In all our simulations the bar center always coincides with the dynamical center, while the stellar disk becomes highly asymmetric during the encounter causing the photometric center of the Magellanic galaxy disk to become mismatched with both the bar and the dynamical center. The disk asymmetries persist for almost 2 Gyrs, the time that it takes for the disk to be re-centered with the bar, and well after the companion has passed. This explains the nature of the offset bar found in many Magellanic-type galaxies, including the Large Magellanic Cloud (LMC) and NGC 3906. In particular, these results, once applied to the LMC, suggest that the dynamical center should reside in the bar center instead of the HI center as previously assumed, pointing to a variation in the current estimate of the north component of the LMC proper motion. △ Less

Submitted 24 August, 2016; v1 submitted 24 February, 2016; originally announced February 2016.

Comments: 12 pages, 8 figures, corrected for referee comments

Journal ref: 2016 ApJ 827 149

Forming disk galaxies in wet major mergers. I. Three fiducial examples

Authors: E. Athanassoula, S. A. Rodionov, N. Peschken, J. C. Lambert

Abstract: Using three fiducial Nbody+SPH simulations, we follow the merging of two disk galaxies with a hot gaseous halo component each, and examine whether the merger remnant can be a spiral galaxy. The stellar progenitor disks are destroyed by violent relaxation during the merging and most of their stars form a classical bulge, while the remaining form a thick disk and its bar. A new stellar disk forms su… ▽ More Using three fiducial Nbody+SPH simulations, we follow the merging of two disk galaxies with a hot gaseous halo component each, and examine whether the merger remnant can be a spiral galaxy. The stellar progenitor disks are destroyed by violent relaxation during the merging and most of their stars form a classical bulge, while the remaining form a thick disk and its bar. A new stellar disk forms subsequently and gradually in the remnant from the gas accreted mainly from the halo. It is vertically thin and well extended in its equatorial plane. A bar starts forming before the disk is fully in place, contrary to what is assumed in idealised simulations of isolated bar-forming galaxies. It has morphological features such as ansae and boxy/peanut bulges. Stars of different ages populate different parts of the box/peanut. A disky pseudobulge forms also, so that by the end of the simulation, all three types of bulges coexist. The oldest stars are found in the classical bulge, followed by those of the thick disk, then by those in the thin disk. The youngest stars are in the spiral arms and the disky pseudobulge. The disk surface density profiles are of type II (exponential with downbending), and the circular velocity curves are flat and show that the disks are submaximum in these examples: two clearly so and one near-borderline between maximum and submaximum. On average, only roughly between 10 and 20% of the stellar mass is in the classical bulge of the final models, i.e. much less than in previous simulations. △ Less

Submitted 13 February, 2016; v1 submitted 9 February, 2016; originally announced February 2016.

Comments: 17 pages, 8 figures, accepted for publication in ApJ. V2: replaced Figure 4 with correct version

Halpha Kinematics of S4G Spiral Galaxies - III. Inner rotation curves

Authors: Santiago Erroz-Ferrer, Johan H. Knapen, Ryan Leaman, Simon Dıaz-Garcia, Heikki Salo, Eija Laurikainen, Miguel Querejeta, Juan Carlos Muñoz-Mateos, E. Athanassoula, Albert Bosma, Sebastien Comeron, Bruce G. Elmegreen, Inma Martınez-Valpuesta

Abstract: We present a detailed study of the shape of the innermost part of the rotation curves of a sample of 29 nearby spiral galaxies, based on high angular and spectral resolution kinematic Halpha Fabry-Perot observations. In particular, we quantify the steepness of the rotation curve by measuring its slope dRvc(0). We explore the relationship between the inner slope and several galaxy parameters, such… ▽ More We present a detailed study of the shape of the innermost part of the rotation curves of a sample of 29 nearby spiral galaxies, based on high angular and spectral resolution kinematic Halpha Fabry-Perot observations. In particular, we quantify the steepness of the rotation curve by measuring its slope dRvc(0). We explore the relationship between the inner slope and several galaxy parameters, such as stellar mass, maximum rotational velocity, central surface brightness (μ0), bar strength and bulge-to-total ratio. Even with our limited dynamical range, we find a trend for low-mass galaxies to exhibit shallower rotation curve inner slopes than high-mass galaxies, whereas steep inner slopes are found exclusively in high-mass galaxies. This trend may arise from the relationship between the total stellar mass and the mass of the bulge, which are correlated among them. We find a correlation between the inner slope of the rotation curve and the morphological T-type, complementary to the scaling relation between dRvc(0) and μ0 previously reported in the literature. Although we find that the inner slope increases with the Fourier amplitude A2 and decreases with the bar torque Qb, this may arise from the presence of the bulge implicit in both A2 and Qb. As previously noted in the literature, the more compact the mass in the central parts of a galaxy (more concretely, the presence of a bulge), the steeper the inner slopes. We conclude that the baryonic matter dominates the dynamics in the central parts of our sample galaxies. △ Less

Submitted 8 February, 2016; originally announced February 2016.

Comments: 11 pages, 11 figures, accepted for publication in MNRAS

Tidally induced bars of galaxies in clusters

Authors: Ewa L. Lokas, Ivana Ebrova, Andres del Pino, Agnieszka Sybilska, E. Athanassoula, Marcin Semczuk, Grzegorz Gajda, Sylvain Fouquet

Abstract: Using N-body simulations we study the formation and evolution of tidally induced bars in disky galaxies in clusters. Our progenitor is a massive, late-type galaxy similar to the Milky Way, composed of an exponential disk and an NFW dark matter halo. We place the galaxy on four different orbits in a Virgo-like cluster and evolve it for 10 Gyr. As a reference case we also evolve the same model in is… ▽ More Using N-body simulations we study the formation and evolution of tidally induced bars in disky galaxies in clusters. Our progenitor is a massive, late-type galaxy similar to the Milky Way, composed of an exponential disk and an NFW dark matter halo. We place the galaxy on four different orbits in a Virgo-like cluster and evolve it for 10 Gyr. As a reference case we also evolve the same model in isolation. Tidally induced bars form on all orbits soon after the first pericenter passage and survive until the end of the evolution. They appear earlier, are stronger, longer and have lower pattern speeds for tighter orbits. Only for the tightest orbit the properties of the bar are controlled by the orientation of the tidal torque from the cluster at pericenters. The mechanism behind the formation of the bars is the angular momentum transfer from the galaxy stellar component to its halo. All bars undergo extended periods of buckling instability that occur earlier and lead to more pronounced boxy/peanut shapes when the tidal forces are stronger. Using all simulation outputs of galaxies at different evolutionary stages we construct a toy model of the galaxy population in the cluster and measure the average bar strength and bar fraction as a function of clustercentric radius. Both are found to be mildly decreasing functions of radius. We conclude that tidal forces can trigger bar formation in cluster cores, but not in the outskirts, and thus cause larger concentrations of barred galaxies towards cluster center. △ Less

Submitted 30 May, 2016; v1 submitted 27 January, 2016; originally announced January 2016.

Comments: 13 pages, 13 figures, accepted for publication in ApJ

Journal ref: ApJ, 826, 227 (2016)