Recent studies of Galactic evolution revealed that the dynamics of the stellar component might be one of the key factors when considering galactic habitability. We run an N-body simulation model of the Milky Way, which we evolve for 10 Gyr, to study the secular evolution of stellar orbits and the resulting galactic habitability related properties, i.e., the density of the stellar component and close stellar encounters. The results indicate that radial migrations are not negligible, even in a simple axisymmetric model with mild levels of dynamical heating, and that the net outward diffusion of the stellar component can populate galactic outskirts with habitable systems. Habitable environment is also likely even at sub-Solar galactocentric radii, because the rate of close encounters should not significantly degrade habitability. Stars that evolve from non-circular to stable nearly circular orbits typically migrate outwards, settling down in a broad Solar neighbourhood. The region between $R \approx 3$ kpc and $R \approx 12$ kpc represents the zone of radial mixing, which can blur the boundaries of the Galactic Habitable Zone (GHZ), as it has been conventionally understood. The present-day stable population of the stars in the Solar neighbourhood originates from this radial mixing zone, with most of the stars coming from the inner regions. The Solar system can be considered as a typical Milky Way habitable system because it migrated outwards from the metal-rich inner regions of the Disk and has a circular orbit in the present epoch. We conclude that the boundaries of the GHZ cannot be sharply confined for a given epoch because of the mixing caused by the stellar migrations and secular evolution of stellar orbits.

We present observations of the Mopra carbon monoxide (CO) survey of the Southern Galactic Plane, covering Galactic longitudes spanning $l = 250^{\circ}$ ($-110^{\circ}$) to $l = 355^{\circ}$ ($-5^{\circ}$), with a latitudinal coverage of at least $|b|<1^\circ$, totalling an area of $>$210 deg$^{2}$. These data have been taken at 0.6 arcmin spatial resolution and 0.1 km s$^{-1}$ spectral resolution, providing an unprecedented view of the molecular gas clouds of the Southern Galactic Plane in the 109–115 GHz $J = 1-0$ transitions of $^{12}$CO, $^{13}$CO, C$^{18}$O, and C$^{17}$O.

We present GECKOS (Generalising Edge-on galaxies and their Chemical bimodalities, Kinematics, and Outflows out to Solar environments), a new ESO VLT/MUSE large program. The main aim of GECKOS is to reveal the variation in key physical processes of disk formation by connecting Galactic Archaeology with integral field spectroscopic observations of nearby galaxies. Edge-on galaxies are ideal for this task: they allow us to disentangle the assembly history imprinted in thick disks and provide the greatest insights into outflows. The GECKOS sample of 35 nearby edge-on disk galaxies is designed to trace the assembly histories and properties of galaxies across a large range of star formation rates, bulge-to-total ratios, and boxy and non-boxy bulges. GECKOS will deliver spatially resolved measurements of stellar abundances, ages, and kinematics, as well as ionised gas metallicities, ionisation param- eters, pressure, and inflow and outflow kinematics; all key parameters for building a complete chemodynamical picture of disk galaxies. With these data, we aim to extend Galactic analysis methods to the wider galaxy population, reaping the benefits of detailed Milky Way studies, while probing the diverse mechanisms of galaxy evolution.

We have studied the validity of the historical Cygnus OB associations and have found that many do not show the kinematic coherence expected for true OB associations. We have revisited these groups by photometrically identifying thousands of OB stars across the region with an SED fitting process which combines photometry, astrometry, spectral and evolutionary models. We applied a flexible clustering method and identified seven kinematically-coherent new OB associations. We observe a distinct correlation between position and velocity for two sets of these associations that suggests an expansion pattern. Tracing the motion of the stars back in the past we find that the sets were at their closest around 7.9 and 8.5 Myr ago. We discuss whether this expansion is a natural by-product of the commonly observed size - velocity dispersion relation of molecular clouds, or requires feedback to initiate the dispersal.

We use the AREPO numerical code to model the structure of a Milky Way like galaxy (MW) via a suite of simulations composed of a stellar disc and bulge, a dark matter halo, and a gaseous disc under isothermal conditions. For each model, we produce longitude velocity (l-v) maps of the gas surface densities to extract the skeletons of the main features (arms, bar), and the contours defining the terminal velocities of the gas. We compare these with observations via a number of diagnostic tools, and select the model that best reproduces the main observed features of the Milky Way.

The velocity-space distribution of the solar neighborhood stars shows complex substructures (moving groups) including the well-known Hercules stream. Recently, the Gaia observation revealed their detailed structures, but their origins are still in debate. We analyzed a high-resolution N-body simulation of a Milky Way (MW)-like galaxy. To find velocity-space distributions similar to that of the solar neighborhood stars, we used Kullback-Leibler divergence (KLD), which is a metric to measure similarities between probability distributions. The KLD analysis shows the time evolution and the spatial variation of the velocity-space distribution. Velocity-space distributions with small KLDs (i.e. high similarities) are frequently but not always detected around in the simulated MW. In the velocity-map with smallest KLD, the velocity-space substructures are made from bar resonances.

In the present work, we have developed a two-dimensional gravitational model of barred galaxies to analyse the fate of escaping stars from the central barred region. For that, the model has been analysed for two different bar profiles viz. strong and weak. Here the phenomena of stellar escape from the central barred region have been studied from the perspective of an open Hamiltonian dynamical system. We observed that the escape routes correspond to the escape basins of the two index-1 saddle points. Our results show that the formation of spiral arms is encouraged for the strong bars. Also, the formation of grand design spirals is more likely for strong bars if they host central super massive black holes (SMBHs). In the absence of central SMBHs, the formation of less-prominent spiral arms is more likely. Again, for weak bars, the formation of inner disc rings is more probable.

We used the spectroscopic and astrometric data provided from the GALactic Archaeology with HERMES (GALAH) Data Release (DR2) and Gaia DR2, respectively, for a large sample of stars to investigate the behaviour of the [ $\alpha$ /Fe] abundances via two procedures, that is, kinematically and spectroscopically. With the kinematical procedure, we investigated the distribution of the [ $\alpha$ /Fe] abundances into the high-/low-probability thin disc, and high-/low-probability thick-disc populations in terms of total space velocity, [Fe/H] abundance, and age. The high-probability thin-disc stars dominate in all sub-intervals of [ $\alpha$ /Fe], including the rich ones: [ $\alpha$ /Fe] $\,>\,0.3$ dex, where the high-probability thick-disc stars are expected to dominate. This result can be explained by the limiting apparent magnitude of the GALAH DR2 ( $V \lt 14$ mag) and intermediate galactic latitude of the star sample. Stars in the four populations share equivalent [ $\alpha$ /Fe] and [Fe/H] abundances, total space velocities, and ages. Hence, none of these parameters can be used alone for separation of a sample of stars into different populations. High-probability thin-disc stars with abundance $-1.3 \lt {\rm[Fe/H]}\leq -0.5$ dex and age $9 \lt \tau\leq13$ Gyr are assumed to have different birth places relative to the metal-rich and younger ones. With the spectroscopic procedure, we separated the sample stars into $\alpha$ -rich and $\alpha$ -poor categories by means of their ages as well as their [ $\alpha$ /Fe] and [Fe/H] abundances. Stars older than 8 Gyr are richer in [ $\alpha$ /Fe] than the younger ones. We could estimate the abundance [ $\alpha$ /Fe] = 0.14 dex as the boundary separating the $\alpha$ -rich and $\alpha$ -poor sub-samples in the [ $\alpha$ /Fe] $\,\times\,$ [Fe/H] plane.

Typical stars in the Milky Way galaxy have velocities of hundreds of kilometres per second and experience gravitational accelerations of $\sim\!10^{-10}~{\rm m\,s}^{-2}$, resulting in velocity changes of a few centimetres per second over a decade. Measurements of these accelerations would permit direct tests of the applicability of Newtonian dynamics on kiloparsec length scales and could reveal significant small-scale inhomogeneities within the galaxy, as well increasing the sensitivity of measurements of the overall mass distribution of the galaxy. Noting that a reasonable extrapolation of progress in exoplanet hunting spectrographs suggests that centimetre per second level precision will be attainable in the coming decade(s), we explore the possibilities such measurements would create. We consider possible confounding effects, including apparent accelerations induced by stellar motion and reflex velocities from planetary systems, along with possible strategies for their mitigation. If these issues can be satisfactorily addressed, it will be possible to use high-precision measurements of changing stellar velocities to perform a ‘blind search’ for dark matter, make direct tests of theories of non-Newtonian gravitational dynamics, detect local inhomogeneities in the dark matter density, and greatly improve measurements of the overall properties of the galaxy.

The Bulge Asymmetries and Dynamical Evolution (BAaDE) survey aims to explore the complex structure of the inner Galaxy and Galactic Bulge, by using the 43 GHz receivers at the Karl G. Jansky Very Large Array (VLA) and the 86 GHz receivers at the Atacama Large Millimeter/submillimeter Array (ALMA) to observe SiO maser lines in red giant stars. The goal is to construct a sample of stellar point-mass probes that can be used to test models of the gravitational potential, and the final sample is expected to provide at least 20,000 line-of-sight velocities and positions. A possible bias between the VLA and the ALMA SiO maser lines is explored, and the 86 GHz SiO line-peak velocities agree using either of the four sampled lines. Additionally, the SiO maser velocities agree with the OH maser derived velocities.

The unprecedented amount and accuracy of kinematic data from the second release of the Gaia mission have started revolutionizing our understanding of the dynamics of the Milky Way disk. The detailed stellar velocity field in the Galactic disk should allow us to constrain with unprecedented precision the parameters of the non-axisymmetric modes of the disk. We present here the status of our current modelling efforts in this area, and their implication on the dynamics of the Galactic bar in particular.

Although the stellar halo accounts for just ∼1% of the total stellar mass of the Milky Way, the kinematics of halo stars can tell us a lot about the origins and evolution of our Galaxy. It has been shown that the high transverse velocity stars in Gaia DR2 reveal a double sequence in the Hertzsprung-Russell (HR) diagram, indicating a duality in the local halo within 1 kpc. We fit these stars by updating the popular Besançon/Galaxia model, incorporating the latest observational results for the stellar halo. We are able to obtain a good match to the Gaia data and provide new constraints on the properties of the disc and halo. In particular, we show that the thick disc contribution to this high velocity tail is small, but not negligible, and likely has an influence on the red sequence of the HR diagram.

We have derived absolute proper motions of stars in the Galactic bulge region combining the VVV InfraRed Astrometric Catalogue (VIRAC) and Gaia. We use the proper motions to study the kinematic structure of the bulge both integrated along the line-of-sight and in magnitude intervals using red clump stars as standard candles. In parallel we compare to a made-to-measure barred dynamical model, folding in the VIRAC selection function, to understand and interpret the structures that we observe. The barred dynamical model, which contains a boxy/peanut bulge, and has a pattern speed of 37.5 kms−1 kpc−1, is able to reproduce all structures impressively well.

The presence of counter-rotating (CR) components in galaxies is not that rare but their origin is still unclear. Important clues to the formation and evolution of CR galaxies are provided by galaxy kinematics, such as the mass distribution and the shape of the gravitational potential. In order to better understand the origin and incidence of CR galaxies, we aim at modeling CR stellar disks, as they would be observed with Integral Field Units (IFU) instruments, and measuring the kinematics of these peculiar astrophysical objects to reveal the CR signatures. In the bi-dimensional maps of analysed models, the double sigma signature is the best diagnostic to spot the presence of a CR disk component.

Dynamical models will be key to exploitation of the incoming flood data for our Galaxy. Modelling techniques are reviewed with an emphasis on f(J) modelling.

The on-going phase mixing in the vertical direction of the Galactic disk has been discovered with the revolutionary Gaia DR2 data. It manifests itself as the snail shell in the Z–Vz phase space. To better understand the origin and properties of the phase mixing process, we study the phase-mixing signatures in moving groups (also known as the kinematic streams) with the Gaia DR2 data in the Galactic disk near the Solar circle. Interestingly, the phase space snail shell exists only in the main kinematic streams with |VR|≲ 50 km/s and |Vφ –VLSR|≲30 km/s, i.e., stars on dynamically “colder” orbits. Compared to the colder orbits, the hotter orbits may have phase-wrapped away already due to the much larger dynamical range in radial variation to facilitate faster phase mixing. These results help put tighter constraints on the vertical perturbation history of the Milky Way disk. To explain the lack of a well-defined snail shell in the hotter orbits, the disk should have been perturbed at least ∼400–500 Myr ago. Our results offer more support to the recent satellite-disk encounter scenario than the internal bar buckling perturbation scenario as the origin of the phase space mixing.

A major uncertainty in the determination of the mass profile of the Milky Way using stellar kinematics in the halo is the poorly determined anisotropy parameter, , where σr is the Galactocentric radial velocity dispersion, and σθ and σφ are the tangential components of the velocity dispersion. We have used a sample of over 24,000 Galactic halo K giant and blue horizontal branch stars from the LAMOST stellar spectroscopic survey and SDSS/SEGUE, combined with proper motions from Gaia Data Release 2, to measure β(rgc) over a wide range of Galactocentric distances rgc from 5 to 80 kpc. Kinematic substructures have been carefully removed to reveal the underlying diffuse stellar halo prior to measuring β. We find that orbits are generally radial (β > 0) and β is constant out to distances of about 40 kpc, with a dependence on metallicity of the stars, such that β declines with lower metallicity. Similar behavior is seen in both the K giant and BHB samples.

Recent advancements in astrometry and in cosmological models of dark matter halo growth have significantly changed our understanding of the dynamics of the Local Group. The most dramatic changes owe to a new picture of the structure and dynamics of the Milky Way’s most massive satellite galaxy, the Large Magellanic Cloud (LMC), which is most likely on its first passage about the Milky Way and ten times larger in mass than previously assumed. The LMC’s orbit through the Milky Way’s dark matter and stellar halo will leave characteristic signatures in both density and kinematics. Furthermore, the gravitational perturbations produced by both direct tidal forcing from the LMC and the response of the halo to its passage will together cause significant perturbations to the orbits of tracers of the Milky Way’s dark matter distribution. We advocate for the use of basis field expansion methods to fully capture and quantify these effects.

The Bulge Asymmetries and Dynamical Evolution (BAaDE) survey aims to use circumstellar SiO maser line-of-sight velocities as probes for the Galactic gravitational potential and dynamical structure. The SiO masers are detected at a high rate in specific color-selected MSX infrared sources. Furthermore, the SiO maser properties and line ratios, in combination with infrared spectral energy distributions and location in the Galaxy, will statistically yield detailed information on population and evolution of low- to intermediate-mass evolved stars in the Galaxy.

Since past few decades, observations have improved so strongly that when modelling Milky Way (MW) dynamics it is required to include small perturbations to the modelling process. It is difficult task that we try to solve by selecting regions to model so small that the perturbation can be considered to give nearly constant effect. We use Solar Neighbourhood (SN) as our test sample and assume that the bar effects show more or less constant contribution to SN. By extrapolating and smoothing observed stars on their orbits, and requiring that smoothed and observed phase space are consistent we were able to deduce acceleration vector. We conclude from non-radial acceleration component that the bar must cause about one third of total acceleration near SN.