Search

To advance our understanding of the evolution of the interstellar medium (ISM) of our Galaxy, numerical models of Milky Way (MW) type galaxies are widely used. However, most models only vaguely resemble the MW (e.g. in total mass), and often use imposed analytic potentials (which cannot evolve dynamically). This poses a problem in asserting their applicability for the interpretation of observations of our own Galaxy. The goal of this work is to identify a numerical model that is not only a MW-type galaxy, but one that can mimic some of the main observed structures of our Galaxy, using dynamically evolving potentials, so that it can be used as a base model to study the ISM cycle in a galaxy like our own. This paper introduces a suite of 15 MW-type galaxy models developed using the {\sc arepo} numerical code, that are compared to Galactic observations of $^{12}$CO and \ion{H}{I} emission via longitude-velocity plots, from where we extract and compare the skeletons of major galactic features and the terminal gas velocities. We found that our best-fitting model to the overall structure, also reproduces some of the more specific observed features of the MW, including a bar with a pattern speed of $30.0 \pm 0.2$ km\,s$^{-1}$\,kpc$^{-1}$, a bar half-length of $3.2 \pm 0.8$\,kpc. Our model shows large streaming motions around spiral arms, and strong radial motions well beyond the inner bar. This model highlights the complex motions of a dynamic MW-type galaxy and has the potential to offer valuable insight into how our Galaxy regulates the ISM and star formation., Comment: Accepted for publication in MNRAS. 24 pages, 23 figures, 3 tables

Different mechanisms driving bar structure formation indicate that bar origins should be distinguishable in the stellar populations of galaxies. To study how these origins affect different bar morphologies and impact stellar orbits and migration, we analyse three simulated discs which are representative of bar formation under isolated evolution motivated by disc instability, and interaction driven tidal development. The first isolated disc and the tidally driven disc produce similar bar structure, while the second isolated disc, generated by the tidal initial condition without the companion, is visibly dissimilar. Changes to radial and vertical positions, angular momentum in the disc-plane, orbital eccentricity and the subsequent disc metallicities are assessed, as is the dependence on stellar age and formation radii. Bar origin is distinguishable, with the tidal disc displaying larger migration overall, higher metallicity difference between the inner and outer disc, as well as a population of inner disc stars displaced to large radii and below the disc-plane. The affect of closest approach on populations of stars formed before, after and during this period is evident. However, bar morphology is also found to be a significant factor in the evolution of disc stellar properties, with similar bars producing similar traits in migration tendency with radius, particularly in vertical stellar motion and in the evolution of central metallicity features., Comment: 18 pages, 12 figures, submitted to the Monthly Notices of the Royal Astronomical Society

Japan Astrometry Satellite Mission for INfrared Exploration (JASMINE) is a planned M-class science space mission by the Institute of Space and Astronautical Science, the Japan Aerospace Exploration Agency. JASMINE has two main science goals. One is the Galactic archaeology with Galactic Center Survey, which aims to reveal the Milky Way's central core structure and formation history from Gaia-level (~25 $\mu$as) astrometry in the Near-Infrared (NIR) Hw-band (1.0-1.6 $\mu$m). The other is the Exoplanet Survey, which aims to discover transiting Earth-like exoplanets in the habitable zone from NIR time-series photometry of M dwarfs when the Galactic center is not accessible. We introduce the mission, review many science objectives, and present the instrument concept. JASMINE will be the first dedicated NIR astrometry space mission and provide precise astrometric information of the stars in the Galactic center, taking advantage of the significantly lower extinction in the NIR. The precise astrometry is obtained by taking many short-exposure images. Hence, the JASMINE Galactic center survey data will be valuable for studies of exoplanet transits, asteroseismology, variable stars and microlensing studies, including discovery of (intermediate mass) black holes. We highlight a swath of such potential science, and also describe synergies with other missions., Comment: 50 pages, 22 figures, accepted for publication in PASJ

It remains unclear how galactic environment affects star formation and stellar cluster properties. This is difficult to address in Milky Way-mass galaxy simulations because of limited resolution and less accurate feedback compared to cloud-scale models. We carry out zoom-in simulations to re-simulate 100-300 pc regions of a Milky Way-like galaxy using smoothed particle hydrodynamics, including finer resolution (0.4 Msun per particle), cluster-sink particles, ray-traced photoionization from O stars, H$_2$/CO chemistry, and ISM heating/cooling. We select $10^6$ Msun cloud complexes from a galactic bar, inner spiral arm, outer arm, and inter-arm region (in order of galactocentric radius), retaining the original galactic potentials. The surface densities of star formation rate and neutral gas follow $\Sigma_{SFR} \propto \Sigma_{gas}^{1.3}$, with the bar lying higher up the relation than the other regions. However, the inter-arm region forms stars 2-3x less efficiently than the arm models at the same $\Sigma_{gas}$. The bar produces the most massive cluster, the inner arm the second, and the inter-arm the third. Almost all clusters in the bar and inner arm are small (radii < 5 pc), while 30-50 per cent of clusters in the outer arm and inter-arm have larger radii more like associations. Bar and inner arm clusters rotate at least twice as fast, on average, than clusters in the outer arm and inter-arm regions. The degree of spatial clustering also decreases from bar to inter-arm. Our results indicate that young massive clusters, potentially progenitors of globular clusters, may preferentially form near the bar/inner arm compared to outer arm/inter-arm regions., Comment: 14 pages, accepted for publication in MNRAS

Context. Filamentary structures in the interstellar medium are closely related to star formation. Dense gas mass fraction (DGMF) or clump formation efficiency in large-scale filaments possibly determine their hosting star formation activities. Aims. We aim to automatically identify large-scale filaments, characterize them, investigate their association with Galactic structures, and study their DGMFs. Methods. We use a modified minimum spanning tree (MST) algorithm to chain parsec-scale 13CO clumps previously extracted from the SEDIGISM (Structure, Excitation, and Dynamics of the Inner Galactic InterStellar Medium) survey. The MST connects nodes in a graph such that the sum of edge lengths is minimum. Modified MST also ensures velocity coherence between nodes, so the identified filaments are coherent in position-position-velocity (PPV) space. Results. We generate a catalog of 88 large-scale ($>10pc$) filaments in the inner Galactic plane (with $-60^\circ < l < 18^\circ and $|b| < 0.5^\circ$). These SEDIGISM filaments are larger and less dense than MST filaments previously identified from the BGPS and ATLASGAL surveys. We find that eight of the filaments run along spiral arms and can be regarded as "bones" of the Milky Way. We also find three bones associated with the Local Spur in PPV space. By compiling 168 large-scale filaments with available DGMF across the Galaxy, an order of magnitude more than previously investigated, we find that DGMFs do not correlate with Galactic location, but bones have higher DGMFs than other filaments., Comment: Accepted for publication in A&A. 30 pages, 8 figures, 3 tables

We present simulations of the formation and evolution of clusters in spiral arms. The simulations follow two different spiral arm regions, and the total gas mass is varied to produce a range of different mass clusters. We find that including photoionizing feedback produces the observed cluster mass radius relation, increasing the radii of clusters compared to without feedback. Supernovae have little impact on cluster properties. We find that in our high density, high gas mass simulations, star formation is less affected by feedback, as star formation occurs rapidly before feedback has much impact. In our lowest gas density simulation, the resulting clusters are completely different (e.g. the number of clusters and their masses) to the case with no feedback. The star formation rate is also significantly suppressed. The fraction of stars in clusters in this model decreases with time flattening at about 20\%. In our lowest gas simulation model, we see the formation of a star forming group with properties similar to an OB association, in particular similar to Orion Ia. We suggest that low densities, and stronger initial dynamics are conducive to forming associations rather than clusters. In all models cluster formation is complex with clusters merging and splitting. The most massive clusters which form have tended to undergo more mergers., Comment: 23 pages, 19 figures, accepted for publication in MNRAS, movies available at https://icybob.co.uk/movies.html (lowest 2 panels)

We study the effects of the magnetic fields on the formation of massive, self-gravitationally bound cores (MBCs) in high-speed cloud-cloud collisions (CCCs). Extending our previous work (Sakre et al. 2021), we perform magnetohydrodynamic simulations following the high-speed (20 - 40 km s$^{-1}$) collisions between two magnetized (4 $\mu$G initially), turbulent clouds of different sizes in the range of 7 - 20 pc. We show that a magnetic field effect hinders the core growth, particularly after a short-duration collision during which cores cannot get highly bound. In such a case, a shocked region created by the collision rapidly expands to the ambient medium owing to the enhanced magnetic pressure, resulting in the destruction of the highly unbound cores and suppression of gas accretion to massive cores. This negative effect on the MBC formation is a phenomenon not seen in the past hydrodynamic simulations of similar CCC models. Together with our previous work, we conclude that the magnetic fields provide the two competing effects on the MBC formation in CCC; while they promote the mass accumulation into cores during a collision, they operate to destroy cores or hinder the core growth after the collision. The duration of collision determines which effect prevails, providing the maximum collision speed for the MBC formation with given colliding clouds. Our results agree with the observed trend among CCC samples in the corresponding column density range; clouds with higher relative velocity require higher column density for the formation of massive stars (Enokiya et al. 2021)., Comment: Submitted to MNRAS journal

We present Atacama Large Millimeter/submillimeter Array observations of multiple CO(1-0), $^{13}$CO(1-0), and C$^{18}$O(1-0) lines and 2.9 mm and 1.3 mm continuum emission toward the nearby interacting luminous infrared galaxy NGC 3110, supplemented with similar spatial resolution H$\alpha$, 1.4GHz continuum, and $K$-band data. We estimate the typical CO-to-H$_2$ conversion factor of 1.7 $M_{\odot}$ (K km s$^{-1}$ pc$^2$)$^{-1}$ within the disk using LTE-based and dust-based H$_2$ column densities, and measure the 1-kpc scale surface densities of star formation rate ($\Sigma_{\rm SFR}$), super star clusters ($\Sigma_{\rm SSC}$), molecular gas mass, and star formation efficiency (SFE) toward the entire gas disk. These parameters show a peak at the southern part of the southern spiral arm (SFE $\sim$ 10$^{-8.2}$ yr$^{-1}$, $\Sigma_{\rm SFR}$ $\sim$ 10$^{-0.6}$ $M_{\odot}$ kpc$^{-2}$ yr$^{-1}$, $\Sigma_{\rm SSC}$ $\sim$ 6.0 kpc$^{-2}$), which is likely attributed to the on-going tidal interaction with the companion galaxy MCG-01-26-013, as well as toward the circumnuclear region. We also find that thermal free-free emission contributes to a significant fraction of the millimeter continuum emission at the southern peak position. Those measurements imply that the peak of the southern arm is an active and young star-forming region, whereas the central part of NGC 3110 is a site of long-continued star formation. We suggest that, during the early stage of the galaxy-galaxy interaction with large mass ratio that in NGC 3110, fragmentation along the main galaxy's arms is an important driver of merger-induced star formation and massive gas inflow results in dusty nuclear starbursts., Comment: 17 pages, 9 figures, 2 tables, accepted for publication in The Astrophysical Journal

Galactic bars are important drivers of galactic evolution, and yet how they impact the interstellar medium and correspondingly star formation, remains unclear. We present simulation results for two barred galaxies with different formation mechanisms, bars formed in isolation or via a tidal interaction, to consider the spatially and temporally varying trends of star formation. We focus on the early (< 1Gyr) epoch of bar formation so that the interaction is clearly identifiable. The nearby NGC 4303 (isolated) and NGC 3627 (interaction history) are selected as observational analogues to tailor these simulations. Regardless of formation mechanism, both models show similar internal dynamical features, although the interaction appears to promote bar-arm disconnection in the outer disc velocity structure. Both bars trigger similar boosts in star formation (79%; 66%), while the interaction also triggers an earlier 31% burst. Significant morphological dependence is observed in the relation between surface gas and star formation rate. In both cases, the bar component is notably steepest; the arm is similar to the overall disc average; and the inter-arm clearly the shallowest. A distinguishable feature of the tidal disc is the presence of moderately dense, inefficiently star forming gas mostly confined to tidal debris outside the optical disc. The tidal disc also exhibits a unique trend of radially increasing star formation efficiency and a clear dearth of star formation which persists along the bar between the centre and bar ends. These are potential signatures for identifying a barred system post-interaction., Comment: 19 pages, 15 figures, accepted for publication in MNRAS

We have performed simulations of cluster formation along two regions of a spiral arm taken from a global Milky Way simulation, including photoionising feedback. One region is characterised by strongly converging flows, the other represents a more typical spiral arm region. We find that more massive clusters are able to form on shorter timescales for the region with strongly converging flows. Mergers between clusters are frequent in the case of the strongly converging flows and enable the formation of massive clusters. We compare equivalent clusters formed in simulations with and without ionisation. Photoionisation does not prevent massive cluster formation, but can be seen to limit the masses of the clusters. On average the mass is reduced by around 20%, but we see a large spread from ionisation having minimal difference to leading to a 50% reduction in mass. Photoionisation is also able to clear out the gas in the vicinity of the clusters on Myr timescales, which can produce clusters with larger radii that are surrounded by more massive stellar halos. We find that the ionising feedback has more impact in our second region which is less dense and has less strongly converging flows., Comment: 20 pages, 18 figures, accepted for publication in MNRAS

The morphology of the Milky Way is still a matter of debate. In order to shed light on uncertainties surrounding the structure of the Galaxy, in this paper, we study the imprint of spiral arms on the distribution and properties of its molecular gas. To do so, we take full advantage of the SEDIGISM survey that observed a large area of the inner Galaxy in the $^{13}$CO(2-1) line at an angular resolution of 28". We analyse the influences of the spiral arms by considering the features of the molecular gas emission as a whole across the longitude-velocity map built from the full survey. Additionally, we examine the properties of the molecular clouds in the spiral arms compared to the properties of their counterparts in the inter-arm regions. Through flux and luminosity probability distribution functions, we find that the molecular gas emission associated with the spiral arms does not differ significantly from the emission between the arms. On average, spiral arms show masses per unit length of $\sim10^5-10^6$ M$_{\odot} $kpc$^{-1}$. This is similar to values inferred from data sets in which emission distributions were segmented into molecular clouds. By examining the cloud distribution across the Galactic plane, we infer that the molecular mass in the spiral arms is a factor of 1.5 higher than that of the inter-arm medium, similar to what is found for other spiral galaxies in the local Universe. We observe that only the distributions of cloud mass surface densities and aspect ratio in the spiral arms show significant differences compared to those of the inter-arm medium; other observed differences appear instead to be driven by a distance bias. By comparing our results with simulations and observations of nearby galaxies, we conclude that the measured quantities would classify the Milky Way as a flocculent spiral galaxy, rather than as a grand-design one., Comment: Accepted for publication in A&A. 38 pages (17 of Appendices), 26 figures, 7 tables. The updated SEDIGISM cloud catalogue, containing spiral arm association information, will be available as part of the SEDIGISM database ( https://sedigism.mpifr-bonn.mpg.de/index.html )

The SEDIGISM (Structure, Excitation and Dynamics of the Inner Galactic Interstellar Medium) survey used the APEX telescope to map 84 deg^2 of the Galactic plane between l = -60 deg and l = +31 deg in several molecular transitions, including 13CO(2-1) and C18O(2-1), thus probing the moderately dense (~10^3 cm^-3) component of the interstellar medium. With an angular resolution of 30'' and a typical 1-sigma sensitivity of 0.8-1.0 K at 0.25 km/s velocity resolution, it gives access to a wide range of structures, from individual star-forming clumps to giant molecular clouds and complexes. The coverage includes a good fraction of the first and fourth Galactic quadrants, allowing us to constrain the large scale distribution of cold molecular gas in the inner Galaxy. In this paper we provide an updated overview of the full survey and the data reduction procedures used. We also assess the quality of these data and describe the data products that are being made publicly available as part of this first data release (DR1). We present integrated maps and position-velocity maps of the molecular gas and use these to investigate the correlation between the molecular gas and the large scale structural features of the Milky Way such as the spiral arms, Galactic bar and Galactic centre. We find that approximately 60 per cent of the molecular gas is associated with the spiral arms and these appear as strong intensity peaks in the derived Galactocentric distribution. We also find strong peaks in intensity at specific longitudes that correspond to the Galactic centre and well known star forming complexes, revealing that the 13CO emission is concentrated in a small number of complexes rather than evenly distributed along spiral arms., Comment: MNRAS, in press

We use the 13CO(2-1) emission from the SEDIGISM high-resolution spectral-line survey of the inner Galaxy, to extract the molecular cloud population with a large dynamic range in spatial scales, using the SCIMES algorithm. This work compiles a cloud catalogue with a total of 10663 molecular clouds, 10300 of which we were able to assign distances and compute physical properties. We study some of the global properties of clouds using a science sample, consisting of 6664 well resolved sources and for which the distance estimates are reliable. In particular, we compare the scaling relations retrieved from SEDIGISM to those of other surveys, and we explore the properties of clouds with and without high-mass star formation. Our results suggest that there is no single global property of a cloud that determines its ability to form massive stars, although we find combined trends of increasing mass, size, surface density and velocity dispersion for the sub-sample of clouds with ongoing high-mass star formation. We then isolate the most extreme clouds in the SEDIGISM sample (i.e. clouds in the tails of the distributions) to look at their overall Galactic distribution, in search for hints of environmental effects. We find that, for most properties, the Galactic distribution of the most extreme clouds is only marginally different to that of the global cloud population. The Galactic distribution of the largest clouds, the turbulent clouds and the high-mass star-forming clouds are those that deviate most significantly from the global cloud population. We also find that the least dynamically active clouds (with low velocity dispersion or low virial parameter) are situated further afield, mostly in the least populated areas. However, we suspect that part of these trends may be affected by some observational biases, and thus require further follow up work in order to be confirmed., Comment: 25 pages (+ appendices, 15 pages), 26 figures, MNRAS

By combining two surveys covering a large fraction of the molecular material in the Galactic disk we investigate the role the spiral arms play in the star formation process. We have matched clumps identified by ATLASGAL with their parental GMCs as identified by SEDIGISM, and use these giant molecular cloud (GMC) masses, the bolometric luminosities, and integrated clump masses obtained in a concurrent paper to estimate the dense gas fractions (DGF$_{\rm gmc}=\sum M_{\rm clump}/M_{\rm gmc}$) and the instantaneous star forming efficiencies (i.e., SFE$_{\rm gmc} = \sum L_{\rm clump}/M_{\rm gmc}$). We find that the molecular material associated with ATLASGAL clumps is concentrated in the spiral arms ($\sim$60% found within $\pm$10 km s$^{-1}$ of an arm). We have searched for variations in the values of these physical parameters with respect to their proximity to the spiral arms, but find no evidence for any enhancement that might be attributable to the spiral arms. The combined results from a number of similar studies based on different surveys indicate that, while spiral-arm location plays a role in cloud formation and HI to H$_2$ conversion, the subsequent star formation processes appear to depend more on local environment effects. This leads us to conclude that the enhanced star formation activity seen towards the spiral arms is the result of source crowding rather than the consequence of a any physical process., Comment: 15 pages, 12 figures. Full list of affiliations can be found at the end of the paper

The nature of galactic spiral arms in disc galaxies remains elusive. Regardless of the spiral model, arms are expected to play a role in sculpting the star-forming interstellar medium. As such, different arm models may result in differences in the structure of the interstellar medium and molecular cloud properties. In this study we present simulations of galactic discs subject to spiral arm perturbations of different natures. We find very little difference in how the cloud population or gas kinematics vary between the different grand-design spirals, indicting that the interstellar medium on cloud scales cares little about where spiral arms come from. We do, however, see a difference in the interarm/arm mass spectra, {and minor differences in tails of the distributions of cloud properties} (as well as radial variations in the stellar/gaseous velocity dispersions). These features can be attributed to differences in the radial dependence of the pattern speeds between the different spiral models, and could act as a metric of the nature of spiral structure in observational studies., Comment: 18 pages, 17 figures. Accepted for publication in MNRAS

$Gaia$ DR2 has provided an unprecedented wealth of information about the positions and motions of stars in our Galaxy, and has highlighted the degree of disequilibria in the disc. As we collect data over a wider area of the disc it becomes increasingly appealing to start analysing stellar actions and angles, which specifically label orbit space, instead of their current phase space location. Conceptually, while $\bar{x}$ and $\bar{v}$ tell us about the potential and local interactions, grouping in action puts together stars that have similar frequencies and hence similar responses to dynamical effects occurring over several orbits. Grouping in actions and angles refines this further to isolate stars which are travelling together through space and hence have shared histories. Mixing these coordinate systems can confuse the interpretation. For example, it has been suggested that by moving stars to their guiding radius, the Milky Way spiral structure is visible as ridge-like overdensities in the $Gaia$ data \citep{Khoperskov+19b}. However, in this work, we show that these features are in fact the known kinematic moving groups, both in the $L_z-\phi$ and the $v_{\mathrm{R}}-v_{\phi}$ planes. Using simulations we show how this distinction will become even more important as we move to a global view of the Milky Way. As an example, we show that the radial velocity wave seen in the Galactic disc in $Gaia$ and APOGEE should become stronger in the action-angle frame, and that it can be reproduced by transient spiral structure., Comment: 12 pages, 10 Figures

With backwards orbit integration we estimate birth locations of young stellar associations and moving groups identified in the solar neighborhood that are younger than 70 Myr. The birth locations of most of these stellar associations are at smaller galactocentric radius than the Sun, implying that their stars moved radially outwards after birth. Exceptions to this rule include are the Argus and Octans associations which formed outside the Sun's Galactocentric radius. Variations in birth heights of the stellar associations suggest that they were born in a corrugated disk of molecular clouds, similar to that inferred from the current filamentary molecular cloud distribution and dust extinction maps. Multiple spiral arm features with different but near corotation pattern speeds and at different heights could account for the stellar association birth sites. We find that the young stellar associations are located in between peaks in the UV stellar velocity distribution for stars in the solar neighborhood. This would be expected if they were born in a spiral arm which perturbs stellar orbits that cross it. In contrast, stellar associations seem to be located near peaks in the vertical phase space distribution, suggesting that the gas in which stellar associations are born moves vertically together with the low velocity dispersion disk stars., Comment: The animated movie is available here: http://astro.pas.rochester.edu/~aquillen/figures/movie3D.gif

We study effect of magnetic field on massive dense core formation in colliding unequal molecular clouds by performing magnetohydrodynamic simulations with sub-parsec resolution (0.015 pc) that can resolve the molecular cores. Initial clouds with the typical gas density of the molecular clouds are immersed in various uniform magnetic fields. The turbulent magnetic fields in the clouds consistent with the observation by Crutcher et al. (2010) are generated by the internal turbulent gas motion before the collision, if the uniform magnetic field strength is 4.0 $\mu$G. The collision speed of 10 km s$^{-1}$ is adopted, which is much larger than the sound speeds and the Alfv\'{e}n speeds of the clouds. We identify gas clumps with gas densities greater than 5 $\times$ 10$^{-20}$ g cm$^{-3}$ as the dense cores and trace them throughout the simulations to investigate their mass evolution and gravitational boundness. We show that a greater number of massive, gravitationally bound cores are formed in the strong magnetic field (4.0 $\mu$G) models than the weak magnetic field (0.1 $\mu$G) models. This is partly because the strong magnetic field suppresses the spatial shifts of the shocked layer that should be caused by the nonlinear thin shell instability. The spatial shifts promote formation of low-mass dense cores in the weak magnetic field models. The strong magnetic fields also support low-mass dense cores against gravitational collapse. We show that the numbers of massive, gravitationally bound cores formed in the strong magnetic field models are much larger than the isolated, non-colliding cloud models, which are simulated for comparison. We discuss the implications of our numerical results on massive star formation., Comment: Published in PASJ (for special issue on cloud-cloud collision). This arXiv version is similar to the published paper ( https://academic.oup.com/pasj/article/doi/10.1093/pasj/psaa059/5867032 )

Far-Ultraviolet (FUV) radiation greatly exceeds ultraviolet, supernovae and winds in the energy budget of young star clusters but is poorly modeled in galaxy simulations. We present results of the first isolated galaxy disk simulations to include photo-electric heating of gas via dust grains from FUV radiation self-consistently, using a ray-tracing approach that calculates optical depths along the source-receiver sight-line. This is the first science application of the TREVR radiative transfer algorithm. We find that FUV radiation alone cannot regulate star formation. However, FUV radiation produces warm neutral gas and is able to produce regulated galaxies with realistic scale heights. FUV is also a long-range feedback and is more important in the outer disks of galaxies. We also use the super-bubble feedback model, which depends only on the supernova energy per stellar mass, is more physically realistic than common, parameter-driven alternatives and thus better constrains supernova feedback impacts. FUV and supernovae together can regulate star formation without producing too much hot ionized medium and with less disruption to the ISM compared to supernovae alone., Comment: accepted by MNRAS

Identifying the structure of our Galaxy has always been fraught with difficulties, and while modern surveys continue to make progress building a map of the Milky Way, there is still much to understand. The arm and bar features are important drivers in shaping the interstellar medium, but their exact nature and influence still require attention. We present results of smoothed particle hydrodynamic simulations of gas in the Milky Way including star formation, stellar feedback, and ISM cooling, when exposed to different arm and bar features, with the aim of better understanding how well newly formed stars trace out the underlying structure of the Galaxy. The bar is given a faster pattern speed than the arms, resulting in a complex, time-dependent morphology and star formation. Inter-arm branches and spurs are easily influenced by the bar, especially in the two-armed spiral models where there is a wide region of resonance overlap in the disc. As the bar over-takes the spiral arms it induces small boosts in star formation and enhances spiral features, which occur at regularly spaced beat-like intervals. The locations of star formation events are similar to those seen in observational data, and do not show a perfect 1:1 correspondence with the underlying spiral potential, though arm tangencies are generally well traced by young stars. Stellar velocity fields from the newly formed stars are compared to data from Gaia DR2, showing that the spiral and bar features can reproduce many of the non-axisymmetric features seen in the data. A simple analytical model is used to show many of these feature are a natural response of gas to rigidly rotating spiral and bar potentials., Comment: 21 pages, 26 figures, accepted for publication in MNRAS. Movies of simulations included in this work can be found at: https://zenodo.org/record/3239994

We compare the properties of clouds in simulated M33 galaxies to those observed in the real M33. We apply a friends of friends algorithm and CPROPS to identify clouds, as well as a pixel by pixel analysis. We obtain very good agreement between the number of clouds, and maximum mass of clouds. Both are lower than occurs for a Milky Way-type galaxy and thus are a function of the surface density, size and galactic potential of M33. We reproduce the observed dependence of molecular cloud properties on radius in the simulations, and find this is due to the variation in gas surface density with radius. The cloud spectra also show good agreement between the simulations and observations, but the exact slope and shape of the spectra depends on the algorithm used to find clouds, and the range of cloud masses included when fitting the slope. Properties such as cloud angular momentum, velocity dispersions and virial relation are also in good agreement between the simulations and observations, but do not necessarily distinguish between simulations of M33 and other galaxy simulations. Our results are not strongly dependent on the level of feedback used here (10 and 20%) although they suggest that 15% feedback efficiency may be optimal. Overall our results suggest that the molecular cloud properties are primarily dependent on the gas and mass surface density, and less dependent on the localised physics such as the details of stellar feedback, or the numerical code used., Comment: 15 pages, 12 figures, accepted by MNRAS

The $^{12}$CO $(J=1\rightarrow0)$ velocity fields of a sample of 20 nearby spiral galaxies, selected from the CO Multi-line Imaging of Nearby Galaxies (COMING) legacy project of Nobeyama Radio Observatory, have been analyzed by Fourier decomposition to determine their basic kinematic properties, such as circular and noncircular velocities. On average, the investigated barred (SAB and SB) galaxies exhibit a ratio of noncircular to circular velocities of molecular gas larger by a factor of 1.5-2 than non-barred (SA) spiral galaxies at radii within the bar semimajor axis $a_\mathrm{b}$ at 1 kpc resolution, with a maximum at a radius of $R/a_\mathrm{b}\sim0.3$. Residual velocity field images, created by subtracting model velocity fields from the data, reveal that this trend is caused by kpc-scale streaming motions of molecular gas in the bar region. Applying a new method based on radial velocity reversal, we estimated the corotation radius $R_\mathrm{CR}$ and bar pattern speed $\Omega_\mathrm{b}$ in seven SAB and SB systems. The ratio of the corotation to bar radius is found to be in a range of $\mathcal{R}\equiv R_\mathrm{CR}/a_\mathrm{b}\sim0.8\mathrm{-}1.6$, suggesting that intermediate (SBb-SBc), luminous barred spiral galaxies host fast and slow rotator bars. Tentative negative correlations are found for $\Omega_\mathrm{b}$ vs. $a_\mathrm{b}$ and $\Omega_\mathrm{b}$ vs. total stellar mass $M_\ast$, indicating that bars in massive disks are larger and rotate slower, possibly a consequence of angular momentum transfer. The kinematic properties of SAB and SB galaxies, derived from Fourier decomposition, are compared with recent numerical simulations that incorporate various rotation curve models and galaxy interactions., Comment: Accepted for publication in PASJ; 39 pages, 15 figures

Mergers of galaxies are an important mode for galaxy evolution because they serve as an efficient trigger of powerful starbursts. However, observational studies of the molecular gas properties during their early stages are scarce. We present interferometric CO(2-1) maps of two luminous infrared galaxies (LIRGs), NGC 3110 and NGC 232, obtained with the Submillimeter Array (SMA) with ~ 1 kpc resolution. While NGC 3110 is a spiral galaxy interacting with a minor (14:1 stellar mass) companion, NGC 232 is interacting with a similarly sized object. We find that such interactions have likely induced in these galaxies enhancements in the molecular gas content and central concentrations, partly at the expense of atomic gas. The obtained molecular gas surface densities in their circumnuclear regions are $\Sigma_{\rm mol}~\gtrsim10^{2.5}$ M$_\odot$ pc$^{-2}$, higher than in non-interacting objects by an order of magnitude. Gas depletion times of ~ 0.5 - 1 Gyr are found for the different regions, lying in between non-interacting disk galaxies and the starburst sequence. In the case of NGC 3110, the spiral arms show on average 0.5 dex shorter depletion times than in the circumnuclear regions if we assume a similar H$_2$-CO conversion factor. We show that even in the early stages of the interaction with a minor companion, a starburst is formed along the circumnuclear region and spiral arms, where a large population of SSCs is found (~350), and at the same time a large central gas concentration is building up which might be the fuel for an active galactic nucleus. The main morphological properties of the NGC 3110 system are reproduced by our numerical simulations and allow us to estimate that the current epoch of the interaction is at ~ 150 Myrs after closest approach., Comment: 30 pages. Accepted for publication in ApJ

With the advent of modern observational efforts providing extensive giant molecular cloud catalogues, understanding the evolution of such clouds in a galactic context is of prime importance. While numerous previous numerical and theoretical works have focused on the cloud properties in isolated discs, few have looked into the cloud population in an interacting disc system. We present results of the first study investigating the evolution of the cloud population in galaxy experiencing an M51-like tidal fly-by using numerical simulations including star formation, interstellar medium cooling and stellar feedback. We see the cloud population shift to large unbound clouds in the wake of the companion passage, with the largest clouds appearing as fleeting short-lived agglomerations of smaller clouds within the tidal spiral arms, brought together by large scale streaming motions. These are then sheared apart as they leave the protection of the spiral arms. Clouds appear to lead diverse lives, even within similar environments, with some being born from gas shocked by filaments streaming into the spiral arms, and others from effectively isolated smaller colliding pairs. Overall this cloud population produces a shallower mass function than the disc in isolation, especially in the arms compared to the inter-arm regions. Direct comparisons to M51 observations show similarities between cloud populations, though models tailored to the mass and orbital models of M51 appear necessary to precisely reproduce the cloud population., Comment: 21 pages, 25 figures, accepted for publication in MNRAS July 25 2018

We perform simulations of isolated galaxies in order to investigate the likely origin of the spiral structure in M33. In our models, we find that gravitational instabilities in the stars and gas are able to reproduce the observed spiral pattern and velocity field of M33, as seen in HI, and no interaction is required. We also find that the optimum models have high levels of stellar feedback which create large holes similar to those observed in M33, whilst lower levels of feedback tend to produce a large amount of small scale structure, and undisturbed long filaments of high surface density gas, hardly detected in the M33 disc. The gas component appears to have a significant role in producing the structure, so if there is little feedback, both the gas and stars organise into clear spiral arms, likely due to a lower combined $Q$ (using gas and stars), and the ready ability of cold gas to undergo spiral shocks. By contrast models with higher feedback have weaker spiral structure, especially in the stellar component, compared to grand design galaxies. We did not see a large difference in the behaviour of $Q_{stars}$ with most of these models, however, because $Q_{stars}$ stayed relatively constant unless the disc was more strongly unstable. Our models suggest that although the stars produce some underlying spiral structure, this is relatively weak, and the gas physics has a considerable role in producing the large scale structure of the ISM in flocculent spirals., Comment: 17 pages, 17 figures, accepted for publication in MNRAS

We explore the effect of different galactic disc environments on the properties of star-forming clouds through variations in the background potential in a set of isolated galaxy simulations. Rising, falling and flat rotation curves expected in halo dominated, disc dominated and Milky Way-like galaxies were considered, with and without an additional two-arm spiral potential. The evolution of each disc displayed notable variations that are attributed to different regimes of stability, determined by shear and gravitational collapse. The properties of a typical cloud were largely unaffected by the changes in rotation curve, but the production of small and large cloud associations was strongly dependent on this environment. This suggests that while differing rotation curves can influence where clouds are initially formed, the average bulk properties are effectively independent of the global environment. The addition of a spiral perturbation made the greatest difference to cloud properties, successfully sweeping the gas into larger, seemingly unbound, extended structures and creating large arm-interarm contrasts., Comment: Accepted to MNRAS on 3rd December, 2017

We present simulations of the gaseous and stellar material in several different galaxy mass models under the influence of different tidal fly-bys to assess the changes in their bar and spiral morphology. Five different mass models are chosen to represent the variety of rotation curves seen in nature. We find a multitude of different spiral and bar structures can be created, with their properties dependent on the strength of the interaction. We calculate pattern speeds, spiral wind-up rates, bar lengths, and angular momentum exchange to quantify the changes in disc morphology in each scenario. The wind-up rates of the tidal spirals follow the 2:1 resonance very closely for the flat and dark matter dominated rotation curves, whereas the more baryon dominated curves tend to wind-up faster, influenced by their inner bars. Clear spurs are seen in most of the tidal spirals, most noticeable in the flat rotation curve models. Bars formed both in isolation and interactions agree well with those seen in real galaxies, with a mixture of "fast" and "slow" rotators. We find no strong correlation between bar length or pattern speed and the interaction strength. Bar formation is, however, accelerated/induced in four out of five of our models. We close by briefly comparing the morphology of our models to real galaxies, easily finding analogues for nearly all simulations presenter here, showing passages of small companions can easily reproduce an ensemble of observed morphologies., Comment: 30 pages, 29 colour figures, accepted for publication in MNRAS. Videos of simulations can be found at http://www.youtube.com/playlist?list=PLQKy--XcWrIVBc1sS2RNc-ekyfeBsGtDs

Tidal encounters are believed to be one of the key drivers of galactic spiral structure in the Universe. Such spirals are expected to produce different morphological and kinematic features compared to density wave and dynamic spiral arms. In this work we present high resolution simulations of a tidal encounter of a small mass companion with a disc galaxy. Included are the effects of gas cooling and heating, star formation and stellar feedback. The structure of the perturbed disc differs greatly from the isolated galaxy, showing clear spiral features that act as sites of new star formation, and displaying interarm spurs. The two arms of the galaxy, the bridge and tail, appear to behave differently; with different star formation histories and structure. Specific attention is focused on offsets between gas and stellar spiral features which can be directly compared to observations. We find some offsets do exist between different media, with gaseous arms appearing mostly on the convex side of the stellar arms, though the exact locations appear highly time dependent. These results further highlight the differences between tidal spirals and other theories of arm structure., Comment: 17 pages, 19 colour figures, accepted for publication in MNRAS

We study the evolution of galactic magnetic fields using 3D smoothed particle magnetohydrodynamics (SPMHD) simulations of galaxies with an imposed spiral potential. We consider the appearance of reversals of the field, and amplification of the field. We find magnetic field reversals occur when the velocity jump across the spiral shock is above $\approx$20km s$^{-1}$, occurring where the velocity change is highest, typically at the inner Lindblad resonance (ILR) in our models. Reversals also occur at corotation, where the direction of the velocity field reverses in the co-rotating frame of a spiral arm. They occur earlier with a stronger amplitude spiral potential, and later or not at all with weaker or no spiral arms. The presence of a reversal at a radii of around 4--6 kpc in our fiducial model is consistent with a reversal identified in the Milky Way, though we caution that alternative Galaxy models could give a similar reversal. We find that relatively high resolution, a few million particles in SPMHD, is required to produce consistent behaviour of the magnetic field. Amplification of the magnetic field occurs in the models, and while some may be genuinely attributable to differential rotation or spiral arms, some may be a numerical artefact. We check our results using Athena, finding reversals but less amplification of the field, suggesting that some of the amplification of the field with SPMHD is numerical., Comment: 15 pages, 13 figures, accepted for publication in MNRAS

Tidal interactions between disc galaxies and low mass companions are an established method for generating galactic spiral features. In this work we present a study of the structure and dynamics of spiral arms driven in interactions between disc galaxies and perturbing companions in 3-D N-body/smoothed hydrodynamical numerical simulations. Our specific aims are to characterize any differences between structures formed in the gas and stars from a purely hydrodynamical and gravitational perspective, and to find a limiting case for spiral structure generation. Through analysis of a number of different interacting cases, we find that there is very little difference between arm morphology, pitch angles and pattern speeds between the two media. The main differences are a minor offset between gas and stellar arms, clear spurring features in gaseous arms, and different radial migration of material in the stronger interacting cases. We investigate the minimum mass of a companion required to drive spiral structure in a galactic disc, finding the limiting spiral generation cases with companion masses of the order $1\times10^9M_\odot$, equivalent to only 4% of the stellar disc mass, or 0.5% of the total galactic mass of a Milky Way analogue., Comment: 20 pages, 23 figures, accepted for publication by MNRAS

The arm structure of the Milky Way remains somewhat of an unknown, with observational studies hindered by our location within the Galactic disc. In the work presented here we use smoothed particle hydrodynamics (SPH) and radiative transfer to create synthetic longitude-velocity observations. Our aim is to reverse-engineer a top down map of the Galaxy by comparing synthetic longitude-velocity maps to those observed. We set up a system of N-body particles to represent the disc and bulge, allowing for dynamic creation of spiral features. Interstellar gas, and the molecular content, is evolved alongside the stellar system. A 3D-radiative transfer code is then used to compare the models to observational data. The resulting models display arm features that are a good reproduction of many of the observed emission structures of the Milky Way. These arms however are dynamic and transient, allowing for a wide range of morphologies not possible with standard density wave theory. The best fitting models are a much better match than previous work using fixed potentials. They favour a 4-armed model with a pitch angle of approximately 20 degrees, though with a pattern speed that decreases with increasing Galactic radius. Inner bars are lacking however, which appear required to fully reproduce the central molecular zone., Comment: 16 pages, 15 figures, accepted by MNRAS

We present an investigation into the morphological features of the Milky Way. We use smoothed particle hydrodynamics (SPH) to simulate the interstellar medium (ISM) in the Milky Way under the effect of a number of different gravitational potentials representing spiral arms and bars, assuming the Milky Way is grand design in nature. The gas is subject to ISM cooling and chemistry, enabling us to track the evolution of molecular gas. We use a 3D radiative transfer code to simulate the emission from the SPH output, allowing for the construction of synthetic longitude-velocity (l-v) emission maps as viewed from the Earth. By comparing these maps with the observed emission in CO from the Milky Way, we infer the arm/bar geometry that provides a best fit to our Galaxy. We find that it is possible to reproduce nearly all features of the l-v diagram in CO emission. There is no model, however, that satisfactorily reproduces all of the features simultaneously. Models with 2 arms cannot reproduce all the observed arm features, while 4 armed models produce too bright local emission in the inner Galaxy. Our best fitting models favour a bar pattern speed within 50-60km/s/kpc and an arm pattern speed of approximately 20km/s/kpc, with a bar orientation of approximately 45 degrees and arm pitch angle between 10-15 degrees., Comment: 25 pages, 27 Figures. Accepted for publication in MNRAS. Animations of a 2-armed and 4-armed simulation can be found at http://youtube.com/watch?v=9dAZIvo9FH8 and http://youtube.com/watch?v=mmYOMgEoaE0

We present narrow-band photometric measurements of the exoplanet GJ 1214b using the 10.4 m Gran Telescopio Canarias (GTC) and the OSIRIS instrument. Using tuneable filters we observed a total of five transits, three of which were observed at two wavelengths nearly simultaneously, producing a total of eight individual light curves, six of these probed the possible existence of a methane absorption feature in the 8770 - 8850 {\AA} region at high resolution. We detect no increase in the planet-to-star radius ratio across the methane feature with a change in radius ratio of $\Delta$R = -0.0007 $\pm$ 0.0017 corresponding to a scale height (H) change of -0.5 $\pm$ 1.2 H across the methane feature, assuming a hydrogen dominated atmosphere. We find a variety of water and cloudy atmospheric models fit the data well, but find that cloud-free models provide poor fits. These observations support a flat transmission spectrum resulting from the presence of a high-altitude haze or a water-rich atmosphere, in agreement with previous studies. In this study the observations are predominantly limited by the photometric quality and the limited number of data points (resulting from a long observing cadence), which make the determination of the systematic noise challenging. With tuneable filters capable of high resolution measurements (R ~ 600 - 750) of narrow absorption features, the interpretation of our results are also limited by the absence of high resolution methane models below 1 $\mu$m., Comment: 12 pages, including 14 figures and 1 table; accepted for publication in MNRAS