$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.

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.

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.

Bayesian inference involves two main computational challenges. First, in estimating the parameters of some model for the data, the posterior distribution may well be highly multi-modal: a regime in which the convergence to stationarity of traditional Markov Chain Monte Carlo (MCMC) techniques becomes incredibly slow. Second, in selecting between a set of competing models the necessary estimation of the Bayesian evidence for each is, by definition, a (possibly high-dimensional) integration over the entire parameter space; again this can be a daunting computational task, although new Monte Carlo (MC) integration algorithms offer solutions of ever increasing efficiency. Nested sampling (NS) is one such contemporary MC strategy targeted at calculation of the Bayesian evidence, but which also enables posterior inference as a by-product, thereby allowing simultaneous parameter estimation and model selection. The widely-used MultiNest algorithm presents a particularly efficient implementation of the NS technique for multi-modal posteriors. In this paper we discuss importance nested sampling (INS), an alternative summation of the MultiNest draws, which can calculate the Bayesian evidence at up to an order of magnitude higher accuracy than `vanilla' NS with no change in the way MultiNest explores the parameter space. This is accomplished by treating as a (pseudo-)importance sample the totality of points collected by MultiNest, including those previously discarded under the constrained likelihood sampling of the NS algorithm. We apply this technique to several challenging test problems and compare the accuracy of Bayesian evidences obtained with INS against those from vanilla NS.

Far-Ultraviolet (FUV) radiation greatly exceeds ultraviolet, supernovae and winds in the energy budget of young star clusters but is poorly modelled in galaxy simulations. We present results of the first full isolated galaxy disk simulations to include FUV radiation self-consistently. 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 is crucial for producing warm neutral gas. 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, reproducing the tight Kennicutt-Schmidt relation observed for an individual galaxy, without producing too much hot ionized medium and with less disruption to the ISM compared to supernovae alone.

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.

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.

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.

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.

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.

We present Phantom, a fast, parallel, modular and low-memory smoothed particle hydrodynamics and magnetohydrodynamics code developed over the last decade for astrophysical applications in three dimensions. The code has been developed with a focus on stellar, galactic, planetary and high energy astrophysics and has already been used widely for studies of accretion discs and turbulence, from the birth of planets to how black holes accrete. Here we describe and test the core algorithms as well as modules for magnetohydrodynamics, self-gravity, sink particles, H_2 chemistry, dust-gas mixtures, physical viscosity, external forces including numerous galactic potentials as well as implementations of Lense-Thirring precession, Poynting-Robertson drag and stochastic turbulent driving. Phantom is hereby made publicly available.

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.

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.

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.

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.

The origin and life-cycle of molecular clouds are still poorly constrained, despite their importance for understanding the evolution of the interstellar medium. We have carried out a systematic, homogeneous, spectroscopic survey of the inner Galactic plane, in order to complement the many continuum Galactic surveys available with crucial distance and gas-kinematic information. Our aim is to combine this data set with recent infrared to sub-millimetre surveys at similar angular resolutions. The SEDIGISM survey covers 78 deg^2 of the inner Galaxy (-60 deg < l < +18 deg, |b| < 0.5 deg) in the J=2-1 rotational transition of 13CO. This isotopologue of CO is less abundant than 12CO by factors up to 100. Therefore, its emission has low to moderate optical depths, and higher critical density, making it an ideal tracer of the cold, dense interstellar medium. The data have been observed with the SHFI single-pixel instrument at APEX. The observational setup covers the 13CO(2-1) and C18O(2-1) lines, plus several transitions from other molecules. The observations have been completed. Data reduction is in progress, and the final data products will be made available in the near future. Here we give a detailed description of the survey and the dedicated data reduction pipeline. Preliminary results based on a science demonstration field covering -20 deg < l < -18.5 deg are presented. Analysis of the 13CO(2-1) data in this field reveals compact clumps, diffuse clouds, and filamentary structures at a range of heliocentric distances. By combining our data with data in the (1-0) transition of CO isotopologues from the ThrUMMS survey, we are able to compute a 3D realization of the excitation temperature and optical depth in the interstellar medium. Ultimately, this survey will provide a detailed, global view of the inner Galactic interstellar medium at an unprecedented angular resolution of ~30".

We examine how three fundamentally different numerical hydrodynamics codes follow the evolution of an isothermal galactic disc with an external spiral potential. We compare an adaptive mesh refinement code (RAMSES), a smoothed particle hydrodynamics code (sphNG), and a volume-discretised meshless code (GIZMO). Using standard refinement criteria, we find that RAMSES produces a disc that is less vertically concentrated and does not reach such high densities as the sphNG or GIZMO runs. The gas surface density in the spiral arms increases at a lower rate for the RAMSES simulations compared to the other codes. There is also a greater degree of substructure in the sphNG and GIZMO runs and secondary spiral arms are more pronounced. By resolving the Jeans' length with a greater number of grid cells we achieve more similar results to the Lagrangian codes used in this study. Other alterations to the refinement scheme (adding extra levels of refinement and refining based on local density gradients) are less successful in reducing the disparity between RAMSES and sphNG/GIZMO. Although more similar, sphNG displays different density distributions and vertical mass profiles to all modes of GIZMO (including the smoothed particle hydrodynamics version). This suggests differences also arise which are not intrinsic to the particular method but rather due to its implementation. The discrepancies between codes (in particular, the densities reached in the spiral arms) could potentially result in differences in the locations and timescales for gravitational collapse, and therefore impact star formation activity in more complex galaxy disc simulations.

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.

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.

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.

We investigate the utility to computational Bayesian analyses of a particular family of recursive marginal likelihood estimators characterized by the (equivalent) algorithms known as "biased sampling" or "reverse logistic regression" in the statistics literature and "the density of states" in physics. Through a pair of numerical examples (including mixture modeling of the well-known galaxy data set) we highlight the remarkable diversity of sampling schemes amenable to such recursive normalization, as well as the notable efficiency of the resulting pseudo-mixture distributions for gauging prior sensitivity in the Bayesian model selection context. Our key theoretical contributions are to introduce a novel heuristic ("thermodynamic integration via importance sampling") for qualifying the role of the bridging sequence in this procedure and to reveal various connections between these recursive estimators and the nested sampling technique.

We present the results from simulations of GMC formation in spiral galaxies. First we discuss cloud formation by cloud-cloud collisions, and gravitational instabilities, arguing that the former is prevalent at lower galactic surface densities and the latter at higher. Cloud masses are also limited by stellar feedback, which can be effective before clouds reach their maximum mass. We show other properties of clouds in simulations with different levels of feedback. With a moderate level of feedback, properties such as cloud rotations and virial parameters agree with observations. Without feedback, an unrealistic population of overly bound clouds develops. Spiral arms are not found to trigger star formation, they merely gather gas into more massive GMCs. We discuss in more detail interactions of clouds in the ISM, and argue that these are more complex than early ideas of cloud-cloud collisions. Finally we show ongoing work to determine whether the Milky Way is a flocculent or grand design spiral.

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.

This work aims at improving the current understanding of the interaction between H ii regions and turbulent molecular clouds. We propose a new method to determine the age of a large sample of OB associations by investigating the development of their associated H ii regions in the surrounding turbulent medium. Using analytical solutions, one-dimensional (1D), and three-dimensional (3D) simulations, we constrained the expansion of the ionized bubble depending on the turbulent level of the parent molecular cloud. A grid of 1D simulations was then computed in order to build isochrone curves for H ii regions in a pressure-size diagram. This grid of models allowed to date large sample of OB associations and was used on the H ii Region Discovery Survey (HRDS). Analytical solutions and numerical simulations showed that the expansion of H ii regions is slowed down by the turbulence up to the point where the pressure of the ionized gas is in a quasi-equilibrium with the turbulent ram pressure. Based on this result, we built a grid of 1D models of the expansion of H ii regions in a profile based on Larson laws. The 3D turbulence is taken into account by an effective 1D temperature profile. The ages estimated by the isochrones of this grid agree well with literature values of well-known regions such as Rosette, RCW 36, RCW 79, and M16. We thus propose that this method can be used to give ages of young OB associations through the Galaxy such as the HRDS survey and also in nearby extra-galactic sources.

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.

We present the current standing of an investigation into the structure of the Milky Way. We use smoothed particle hydrodynamics (SPH) to simulate the ISM gas in the Milky Way under the effect of a number of different gravitational potentials representing the spiral arms and nuclear bars, both fixed and time-dependent. The gas is subject to ISM cooling and chemistry, enabling us to track the CO and HI density. We use a 3D grid-based radiative transfer code to simulate the emission from the SPH output, allowing for the construction of synthetic longitude-velocity maps as viewed from the Earth. By comparing these maps with the observed emission in CO and HI from the Milky Way (Dame et al. 2001, Kalberla et al. 2005), we can infer the arm/bar geometry that provides a best fit to our Galaxy. By doing so we aim to answer key questions concerning the morphology of the Milky Way such as the number of the spiral arms, the pattern speeds of the bar(s) and arms, the pitch angle of the arms and shape of the bar(s)

"Approximate Bayesian Computation" (ABC) represents a powerful methodology for the analysis of complex stochastic systems for which the likelihood of the observed data under an arbitrary set of input parameters may be entirely intractable-the latter condition rendering useless the standard machinery of tractable likelihood-based, Bayesian statistical inference (e.g. conventional Markov Chain Monte Carlo simulation; MCMC). In this article we demonstrate the potential of ABC for astronomical model analysis by application to a case study in the morphological transformation of high redshift galaxies. To this end we develop, first, a stochastic model for the competing processes of merging and secular evolution in the early Universe; and second, through an ABC-based comparison against the observed demographics of massive (M_gal > 10^11 M_sun) galaxies (at 1.5 < z < 3) in the CANDELS/EGS dataset we derive posterior probability densities for the key parameters of this model. The "Sequential Monte Carlo" (SMC) implementation of ABC exhibited herein, featuring both a self-generating target sequence and self-refining MCMC kernel, is amongst the most efficient of contemporary approaches to this important statistical algorithm. We highlight as well through our chosen case study the value of careful summary statistic selection, and demonstrate two modern strategies for assessment and optimisation in this regard. Ultimately, our ABC analysis of the high redshift morphological mix returns tight constraints on the evolving merger rate in the early Universe and favours major merging (with disc survival or rapid reformation) over secular evolution as the mechanism most responsible for building up the first generation of bulges in early-type disks.