The spiral structure of the Galaxy revealed by CS sources and evidence for the 4:1 resonance

We present a map of the spiral structure of the Galaxy, as traced by molecular carbon monosulphide (CS) emission associated with IRAS sources which are believed to be compact H II regions. The CS line velocities are used to determine the kinematic distances of the sources in order to investigate their distribution in the galactic plane. This allows us to use 870 objects to trace the arms, a number larger than that of previous studies based on classical H II regions. The distance ambiguity of the kinematic distances, when it exists, is solved by different procedures, including the latitude distribution and an analysis of the longitude-velocity diagram. The study of the spiral structure is complemented with other tracers: open clusters, Cepheids, methanol masers and H II regions. The well-defined spiral arms are seen to be confined inside the corotation radius, as is often the case in spiral galaxies. We identify a square-shaped sub-structure in the CS map with that predicted by stellar orbits at the 4:1 resonance (four epicycle oscillations in one turn around the galactic centre). The sub-structure is found at the expected radius, based on the known pattern rotation speed and epicycle frequency curve. An inner arm presents an end with strong inwards curvature and intense star formation that we tentatively associate with the region where this arm surrounds the extremity of the bar, as seen in many barred galaxies. Finally, a new arm with concave curvature is found in the Sagitta to Cepheus region of the sky. The observed arms are interpreted in terms of perturbations similar to grooves in the gravitational potential of the disc, produced by crowding of stellar orbits.

Hydrodynamical simulations of Galactic fountains - I. Evolution of single fountains

The ejection of the gas out of the disc in late-type galaxies is related to star formation and is due mainly to Type II supernovae. In this paper, we studied in detail the development of the Galactic fountains in order to understand their dynamical evolution and their influence on the redistribution of the freshly delivered metals over the disc. To this aim, we performed a number of 3D hydrodynamical radiative cooling simulations of the gas in the Milky Way where the whole Galaxy structure, the Galactic differential rotation and the supernova explosions generated by a single OB association are considered. A typical fountain powered by 100 Type II supernovae may eject material up to similar to 2 kpc which than collapses back mostly in the form of dense, cold clouds and filaments. The majority of the gas lifted up by the fountains falls back on the disc remaining within a radial distance Delta R = 0.5 kpc from the place where the fountain originated. This localized circulation of disc gas does not influence the radial chemical gradients on large scale, as required by the chemical models of the Milky Way which reproduce the metallicity distribution without invoking large fluxes of metals. Simulations of multiple fountains fuelled by Type II supernovae of different OB associations will be presented in a companion paper.

The corotation gap in the Galactic H i distribution

We used the H i data from the LAB Survey to map the ring-shaped gap in H i density that lies slightly outside the solar circle. Adopting R(0) = 7.5 kpc, we find an average gap radius of 8.3 kpc and an average gap width of 0.8 kpc. The characteristics of the H i gap correspond closely to the expected ones, as predicted by theory and by numerical simulations of the gas flow near the corotation resonance.

Hydrodynamical simulations of Galactic fountains - II. Evolution of multiple fountains

The ejection of gas out of the disc in late-type galaxies is related to star formation and is mainly due to the explosion of Type II supernovae (SN II). In a previous paper, we considered the evolution of a single Galactic fountain, that is, a fountain powered by a single SN cluster. Using three-dimensional hydrodynamical simulations, we studied in detail the fountain flow and its dependence with several factors, such as the Galactic rotation, the distance to the Galactic centre and the presence of a hot gaseous halo. As a natural followup, this paper investigates the dynamical evolution of multiple generations of fountains generated by similar to 100 OB associations. We have considered the observed size-frequency distribution of young stellar clusters within the Galaxy in order to appropriately fuel the multiple fountains in our simulations. Most of the results of the previous paper have been confirmed, like for example the formation of intermediate velocity clouds above the disc by the multiple fountains. Also, this work confirms the localized nature of the fountain flows: the freshly ejected metals tend to fall back close to the same Galactocentric region where they are delivered. Therefore, the fountains do not change significantly the radial profile of the disc chemical abundance. The multiple fountain simulations also allowed us to consistently calculate the feedback of the star formation on the halo gas. We found that the hot gas gains about 10 per cent of all the SN II energy produced in the disc. Thus, the SN feedback more than compensate for the halo radiative losses and allow for a quasi steady-state disc-halo circulation to exist. Finally, we have also considered the possibility of mass infall from the intergalactic medium and its interaction with the clouds that are formed by the fountains. Though our simulations are not suitable to reproduce the slow rotational pattern that is typically observed in the haloes around the disc galaxies, they indicate that the presence of an external gas infall may help to slow down the rotation of the gas in the clouds and thus the amount of angular momentum that they transfer to the coronal gas, as previously suggested in the literature.

A mechanism for the formation of oxygen and iron bimodal radial distribution in the disc of our Galaxy

Recently, it has been proposed that there are two type Ia supernova progenitors: short-lived and long-lived. On the basis of this idea, we develop a theory of a unified mechanism for the formation of the bimodal radial distribution of iron and oxygen in the Galactic disc. The underlying cause for the formation of the fine structure of the radial abundance pattern is the influence of the spiral arms, specifically the combined effect of the corotation resonance and turbulent diffusion. From our modelling, we conclude that in order to explain the bimodal radial distributions simultaneously for oxygen and iron and to obtain approximately equal total iron output from different types of supernovae, the mean ejected iron mass per supernova event should be the same as quoted in the literature if the maximum mass of stars, which eject heavy elements, is 50 M(circle dot). For the upper mass limit of 70 M(circle dot), the production of iron by a type II supernova explosion should increase by about 1.5 times.

Approximate Bayesian Computation for astronomical model analysis : a case study in galaxy demographics and morphological transformation at high redshift

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 (MCMC) simulation]. In this paper, 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 secondly, through an ABC-based comparison against the observed demographics of massive (Mgal > 1011 M⊙) galaxies (at 1.5 < z < 3) in the Cosmic Assembly Near-IR Deep Extragalatic Legacy Survey (CANDELS)/Extended Groth Strip (EGS) data set we derive posterior probability densities for the key parameters of this model. The ‘Sequential Monte Carlo’ 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 optimization 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 discs.

Dark matter dominance at all radii in the superthin galaxy UGC 7321

We model the shape and density profile of the dark matter halo of the low surface brightness, superthin galaxy UGC 7321, using the observed rotation curve and the H i scale height data as simultaneous constraints. We treat the galaxy as a gravitationally coupled system of stars and gas, responding to the gravitational potential of the dark matter halo. An isothermal halo of spherical shape with a core density in the range of View the MathML source and a core radius between 2.5 and 2.9 kpc, gives the best fit to the observations for a range of realistic gas parameters assumed. We find that the best-fit core radius is only slightly higher than the stellar disc scale length (2.1 kpc), unlike the case of the high surface brightness galaxies where the halo core radius is typically 3–4 times the disc scale length of the stars. Thus our model shows that the dark matter halo dominates the dynamics of the low surface brightness, superthin galaxy UGC 7321 at all radii, including the inner parts of the galaxy.

Why are some galaxy discs extremely thin?

Some low-surface-brightness galaxies are known to have extremely thin stellar discs with the vertical-to-planar axes ratio 0.1 or less, often referred to as superthin galaxies. Although their existence is now known for over three decades, the physical origin of the superthin discs is still not understood. We model the vertical thickness of the stellar disc using our model of a two-component (gravitationally coupled stars and gas) disc embedded in a dark matter halo, for a bulgeless, superthin galaxy UGC 7321 which has a dense, compact halo, and is compare with a typical dwarf irregular galaxy Holmberg II which has a low-density, non-compact halo. We show that while the presence of gas does constrain the stellar disc thickness and hence its axial ratio, it is the compact dark matter halo which plays the decisive role in determining the mean distribution of stars in the vertical direction in low-luminosity bulgeless galaxies like UGC 7321, and causes the stellar disc to be superthin. Thus, the compactness of the dark matter halo significantly affects the disc structure and this could be important for the early evolution of galaxies.

Narrow escape: how ionizing photons escape from disc galaxies

In this paper we calculate the escape fraction (f(esc)) of ionizing photons from starburst galaxies. Using 2D axisymmetric hydrodynamic simulations, we study superbubbles created by overlapping supernovae in OB associations. We calculate the escape fraction of ionizing photons from the centre of the disc along different angles through the superbubble and the gas disc. After convolving with the luminosity function of OB associations, we show that the ionizing photons escape within a cone of similar to 40 degrees, consistent with observations of nearby galaxies. The evolution of the escape fraction with time shows that it falls initially as cold gas is accumulated in a dense shell. After the shell crosses a few scaleheights and fragments, the escape fraction through the polar regions rises again. The angle-averaged escape fraction cannot exceed similar to1 - cos (1 rad)] = 0.5 from geometrical considerations (using the emission cone opening angle). We calculate the dependence of the time-and angle-averaged escape fraction on the mid-plane disc gas density (in the range n(0) = 0.15-50 cm(-3)) and the disc scaleheight (between z(0) = 10 and 600 pc). We find that the escape fraction is related to the disc parameters (the mid-plane disc density and scaleheight) roughly so that f(esc)(alpha)n(0)(2)z(0)(3) (with alpha approximate to 2.2) is a constant. For discs with a given warm neutral medium temperature, massive discs have lower escape fraction than low-mass galaxies. For Milky Way ISM parameters, we find f(esc) similar to 5 per cent, and it increases to approximate to 10 per cent for a galaxy 10 times less massive. We discuss the possible effects of clumpiness of the ISM on the estimate of the escape fraction and the implications of our results for the reionization of the Universe.

On the relativistic iron line and soft excess in the Seyfert 1 galaxy Markarian 335

We report a 133-ks XMM-Newton observation of the Seyfert 1 galaxy Markarian 335. The 0.4-12 keV spectrum contains an underlying power-law continuum, a soft excess below 2 keV, and a double-peaked iron emission feature in the 6-7 keV range. We investigate the possibility that the double-peaked emission might represent the characteristic signature of the accretion disc. Detailed investigations show that a moderately broad accretion disc line is most likely present, but that the peaks may be due to narrower components from more distant material. The peaks at 6.4 and 7 keV can be identified, respectively, with the molecular torus in active galactic nucleus unification schemes, and very highly ionized, optically thin gas filling the torus. The X-ray variability spectra on both long (~100 ks) and short (~1 ks) time-scales do not support the recent suggestion that the soft excess is an artefact of variable, moderately ionized absorption. © 2007 The Authors. Journal compilation © 2007 RAS.

A simple disc wind model for broad absorption line quasars

Approximately 20 per cent of quasi-stellar objects (QSOs) exhibit broad, blue-shifted absorption lines in their ultraviolet spectra. Such features provide clear evidence for significant outflows from these systems, most likely in the form of accretion disc winds. These winds may represent the ‘quasar’ mode of feedback that is often invoked in galaxy formation/evolution models, and they are also key to unification scenarios for active galactic nuclei (AGN) and QSOs. To test these ideas, we construct a simple benchmark model of an equatorial, biconical accretion disc wind in a QSO and use a Monte Carlo ionization/radiative transfer code to calculate the ultraviolet spectra as a function of viewing angle. We find that for plausible outflow parameters, sightlines looking directly into the wind cone do produce broad, blue-shifted absorption features in the transitions typically seen in broad absorption line (BAL) QSOs. However, our benchmark model is intrinsically X-ray weak in order to prevent overionization of the outflow, and the wind does not yet produce collisionally excited line emission at the level observed in non-BAL QSOs. As a first step towards addressing these shortcomings, we discuss the sensitivity of our results to changes in the assumed X-ray luminosity and mass-loss rate, Ṁwind. In the context of our adopted geometry, Ṁwind ∼ Ṁacc is required in order to produce significant BAL features. The kinetic luminosity and momentum carried by such outflows would be sufficient to provide significant feedback.

The old metal-poor open cluster ESO 92-SC05: accreted from a dwarf galaxy ?

The study of old open clusters outside the solar circle can bring constraints on formation scenarios of the outer disc. In particular, accretion of dwarf galaxies has been proposed as a likely mechanism in the area. We use BVI photometry for determining fundamental parameters of the faint open cluster ESO 92-SC05. Colour-magnitude diagrams are compared with Padova isochrones, in order to derive age, reddening and distance. We derive a reddening E(B - V) = 0.17, and an old age of similar to 6.0 Gyr. It is one of the rare open clusters known to be older than 5 Gyr. A metallicity of Z similar to 0.004 or [M/H] similar to -0.7 is found. The rather low metallicity suggests that this cluster might be the result of an accretion episode of a dwarf galaxy.

Measurement of the epicycle frequency in the Galactic disc and initial velocities of open clusters

A new method to measure the epicycle frequency kappa in the Galactic disc is presented. We make use of the large data base on open clusters completed by our group to derive the observed velocity vector (amplitude and direction) of the clusters in the Galactic plane. In the epicycle approximation, this velocity is equal to the circular velocity given by the rotation curve, plus a residual or perturbation velocity, of which the direction rotates as a function of time with the frequency kappa. Due to the non-random direction of the perturbation velocity at the birth time of the clusters, a plot of the present-day direction angle of this velocity as a function of the age of the clusters reveals systematic trends from which the epicycle frequency can be obtained. Our analysis considers that the Galactic potential is mainly axis-symmetric, or in other words, that the effect of the spiral arms on the Galactic orbits is small; in this sense, our results do not depend on any specific model of the spiral structure. The values of kappa that we obtain provide constraints on the rotation velocity of the in particular, V(0) is found to be 230 +/- 15 km s(-1) even if the scale (R(0) = 7.5 kpc) of the Galaxy is adopted. The measured kappa at the solar radius is 43 +/- 5 km s(-1) kpc(-1). The distribution of initial velocities of open clusters is discussed.

NGC922 - a New Drop-through Ring Galaxy

We have found the peculiar galaxy NGC 922 to be a new drop-through ring galaxy using multiwavelength (ultraviolet-radio) imaging and spectroscopic observations. Its 'C'-shaped morphology and tidal plume indicate a recent strong interaction with its companion which was identified with these observations. Using numerical simulations we demonstrate that the main properties of the system can be generated by a high-speed off-axis drop-through collision of a small galaxy with a larger disc system, thus making NGC 922 one of the nearest known collisional ring galaxies. While these systems are rare in the local Universe, recent deep Hubble Space Telescope images suggest they were more common in the early Universe.