Galaxies in Clusters: the Observational Characteristics of Bow-Shocks, Wakes and Tails

The dynamical signatures of the interaction between galaxies in clusters and the intracluster medium (ICM) can potentially yield significant information about the structure and dynamical history of clusters. To develop our understanding of this phenomenon we present results from numerical modelling of the galaxy/ICM interaction, as the galaxy moves through the cluster. The simulations have been performed for a broad range, of ICM temperatures (kT = 1,4 and 8 keV), representative of poor clusters or groups through to rich clusters. There are several dynamical features that can be identified in these simulations; for supersonic galaxy motion, a leading bow-shock is present, and also a weak gravitationally focussed wake or tail behind the galaxy (analogous to Bondi-Hoyle accretion). For galaxies with higher mass-replenishment rates and a denser interstellar medium (ISM), the dominant feature is a dense ram-pressure stripped tail. In line with other simulations, we find that the ICM/galaxy ISM interaction can result in complex time- dependent dynamics, with ram-pressure stripping occurring in an episodic manner. In order to facilitate this comparison between the observational consequences of dynamical studies and X-ray observations we have calculated synthetic X-ray flux and hardness maps from these simulations. These calculations predict that the ram-pressure stripped tail will usually be the most visible feature, though in nearby galaxies the bow-shock preceding the galaxy should also be apparent in deeper X-ray observations. We briefly discuss these results and compare with X-ray observations of galaxies where there is evidence of such interactions.Comment: 14 pages, 8 diagrams, MNRAS (in press

Simulations of the Effects of Stripping and Accretion on Galaxy Haloes in Clusters

We present results from a series of hydrodynamic simulations investigating ram pressure stripping of galactic haloes as the host galaxy falls radially into a cluster. We perform a parameter study comprising of variations in initial gas content, gas injection rate (via stellar mass loss processes), galaxy mass and amplitude of infall. From the simulation results we track variations in both physical quantities (e.g. gas mass) and directly observable quantities (e.g. X-ray luminosities). The luminosity of the galaxy's X-ray halo is found to compare favourably with the observationally determined correlation with optical blue band luminosity (L_X:L_B) relation. Factors affecting the X-ray luminosity are explored and it is found that the gas injection rate is a dominant factor in determining the integrated luminosity. Observational properties of the material stripped from the galaxy, which forms an X-ray wake, are investigated and it is found that wakes are most visible around galaxies with a substantial initial gas content, during their first passage though the cluster. We define a statistical skewness measure which may be used to determine the direction of motion of a galaxy using X-ray observations. Structures formed in these simulations are similar to the cold fronts seen in observation of cluster mergers where a sharp increase in surface brightness is accompanied by a transition to a cooler region.Comment: Accepted for publication in MNRAS. 19 pages, 21 figure

The morphology of the Milky Way - II. Reconstructing CO maps from disc galaxies with live stellar distributions

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 MNRA

Using the UM dynamical cores to reproduce idealised 3D flows

We demonstrate that both the current (New Dynamics), and next generation (ENDGame) dynamical cores of the UK Met Office global circulation model, the UM, reproduce consistently, the long-term, large-scale flows found in several published idealised tests. The cases presented are the Held-Suarez test, a simplified model of Earth (including a stratosphere), and a hypothetical tidally locked Earth. Furthermore, we show that using simplifications to the dynamical equations, which are expected to be justified for the physical domains and flow regimes we have studied, and which are supported by the ENDGame dynamical core, also produces matching long-term, large-scale flows. Finally, we present evidence for differences in the detail of the planetary flows and circulations resulting from improvements in the ENDGame formulation over New Dynamics.Comment: 34 Pages, 23 Figures. Accepted for publication in Geoscientific Model Development (pre-proof version

Results from a set of three-dimensional numerical experiments of a hot Jupiter atmosphere

We present highlights from a large set of simulations of a hot Jupiter atmosphere, nominally based on HD 209458b, aimed at exploring both the evolution of the deep atmosphere, and the acceleration of the zonal flow or jet. We find the occurrence of a super-rotating equatorial jet is robust to changes in various parameters, and over long timescales, even in the absence of strong inner or bottom boundary drag. This jet is diminished in one simulation only, where we strongly force the deep atmosphere equator-to-pole temperature gradient over long timescales. Finally, although the eddy momentum fluxes in our atmosphere show similarities with the proposed mechanism for accelerating jets on tidally-locked planets, the picture appears more complex. We present tentative evidence for a jet driven by a combination of eddy momentum transport and mean flow.Comment: 26 pages, 22 Figures. Accepted for publication in Astronomy and Astrophysic

Accuracy tests of radiation schemes used in hot Jupiter global circulation models

The treatment of radiation transport in global circulation models (GCMs) is crucial to correctly describe Earth and exoplanet atmospheric dynamics processes. The two-stream approximation and correlated-k method are currently state-of-the-art approximations applied in both Earth and hot Jupiter GCM radiation schemes to facilitate rapid calculation of fluxes and heating rates. Their accuracy have been tested extensively for Earth-like conditions, but verification of the methods' applicability to hot Jupiter-like conditions is lacking in the literature. We are adapting the UK Met Office GCM, the Unified Model (UM), for the study of hot Jupiters, and present in this work the adaptation of the Edwards-Slingo radiation scheme based on the two-stream approximation and the correlated-k method. We discuss the calculation of absorption coefficients from high temperature line lists and highlight the large uncertainty in the pressure-broadened line widths. We compare fluxes and heating rates obtained with our adapted scheme to more accurate discrete ordinate (DO) line-by-line (LbL) calculations ignoring scattering effects. We find that, in most cases, errors stay below 10 % for both heating rates and fluxes using ~ 10 k-coefficients in each band and a diffusivity factor D = 1.66. The two-stream approximation and the correlated-k method both contribute non-negligibly to the total error. We also find that using band-averaged absorption coefficients, which have previously been used in radiative-hydrodynamical simulations of a hot Jupiter, may yield errors of ~ 100 %, and should thus be used with caution.European Community’s Seventh Framework Programme FP7/2007-2013Science & Technology Facilities Council (STFC)Royal Societ

A Synthetic 21-cm Galactic Plane Survey of an SPH Galaxy Simulation

We have created synthetic neutral hydrogen (HI) Galactic Plane Survey data cubes covering 90 degrees < l < 180 degrees, using a model spiral galaxy from SPH simulations and the radiative transfer code TORUS. The density, temperature and other physical parameters are fed from the SPH simulation into TORUS, where the HI emissivity and opacity are calculated before the 21-cm line emission profile is determined. Our main focus is the observation of Outer Galaxy `Perseus Arm' HI, with a view to tracing atomic gas as it encounters shock motions as it enters a spiral arm interface, an early step in the formation of molecular clouds. The observation of HI self-absorption features at these shock sites (in both real observations and our synthetic data) allows us to investigate further the connection between cold atomic gas and the onset of molecular cloud formation.Comment: MNRAS accepted; 11 pages, 12 figure

The structure of HI in galactic disks: Simulations vs observations

We generate synthetic HI Galactic plane surveys from spiral galaxy simulations which include stellar feedback processes. Compared to a model without feedback we find an increased scale height of HI emission (in better agreement with observations) and more realistic spatial structure (including supernova blown bubbles). The synthetic data show HI self-absorption with a morphology similar to that seen in observations. The density and temperature of the material responsible for HI self-absorption is consistent with observationally determined values, and is found to be only weakly dependent on absorption strength and star formation efficiency.Comment: 12 pages, 7 figures. Accepted for publication in MNRA

Modelling circumstellar discs with 3D radiation hydrodynamics

We present results from combining a grid-based radiative transfer code with a Smoothed Particle Hydrodynamics code to produce a flexible system for modelling radiation hydrodynamics. We use a benchmark model of a circumstellar disc to determine a robust method for constructing a gridded density distribution from SPH particles. The benchmark disc is then used to determine the accuracy of the radiative transfer results. We find that the SED and the temperature distribution within the disc are sensitive to the representation of the disc inner edge, which depends critically on both the grid and SPH resolution. The code is then used to model a circumstellar disc around a T-Tauri star. As the disc adjusts towards equilibrium vertical motions in the disc are induced resulting in scale height enhancements which intercept radiation from the central star. Vertical transport of radiation enables these perturbations to influence the mid-plane temperature of the disc. The vertical motions decay over time and the disc ultimately reaches a state of simultaneous hydrostatic and radiative equilibrium.Comment: MNRAS accepted; 15 pages; 17 figures, 4 in colou