On the origin of warps and the role of the intergalactic medium

There is still no consensus as to what causes galactic discs to become warped. Successful models should account for the frequent occurrence of warps in quite isolated galaxies, their amplitude as well as the observed azimuthal and vertical distributions of the H  i layer. Intergalactic accretion flows and intergalactic magnetic fields may bend the outer parts of spiral galaxies. In this paper we consider the viability of these non-gravitational torques to take the gas off the plane. We show that magnetically generated warps are clearly flawed because they would wrap up into a spiral in less than two or three galactic rotations. The inclusion of any magnetic diffusivity to dilute the wrapping effect causes the amplitude of the warp to damp. We also consider the observational consequences of the accretion of an intergalactic plane-parallel flow at infinity. We have computed the amplitude and warp asymmetry in the accretion model, for a disc embedded in a flattened dark matter halo, including self-consistently the contribution of the modes with azimuthal wavenumbers   m = 0  and   m = 1  . Since the m = 0 component, giving a U-shaped profile, is not negligible compared to the m = 1 component, this model predicts quite asymmetric warps, maximum gas displacements on the two sides in the ratio 3 : 2 for the preferred Galactic parameters, and the presence of a fraction ∼3.5 per cent of U-shaped warps, at least. The azimuthal dependence of the moment transfer by the ram pressure would produce a strong asymmetry in the thickness of the H  i layer and asymmetric density distributions in z , in conflict with observational data for the warp in our Galaxy and in external galaxies. The amount of accretion that is required to explain the Galactic warp would give gas scaleheights in the far outer disc that are too small. We conclude that accretion of a flow with no net angular momentum cannot be the main and only cause of warps.     

Vertical distribution of Galactic disc stars and gas constrained by a molecular cloud complex

We investigate the dynamical effects of a molecular cloud complex with a mass ∼ 10⁷ M_⊙ and a size ∼ a few 100 pc on the vertical distribution of stars and atomic hydrogen gas in a spiral galactic disc. Such massive complexes have now been observed in a number of spiral galaxies. The extended mass distribution in a complex, with an average mass density 6 times higher than the Oort limit, is shown to dominate the local gravitational field. This results in a significant redistribution or clustering of the surrounding disc components towards the mid-plane, with a resulting decrease in their vertical scaleheights.
The modified, self-consistent stellar density distribution is obtained by solving the combined Poisson equation and the force equation along the z -direction for an isothermal stellar disc on which the complex is imposed. The effect of the complex is strongest at its centre, where the stellar mid-plane density increases by a factor of 2.6 and the vertical scaleheight decreases by a factor of 3.4 compared with the undisturbed stellar disc. A surprising result is the large radial distance of ∼ 500 pc from the complex centre over which the complex influences the disc; this is due to the extended mass distribution in a complex. The complex has a comparable effect on the vertical distribution of the atomic hydrogen gas in the galactic disc. This 'pinching' or constraining effect should be detectable in the nearby spiral galaxies, as for example has been done for NGC 2403 by Sicking. Thus the gravitational field of a complex results in local corrugations of the stellar and H  i vertical scaleheights, and the galactic disc potential is highly non-uniform on scales of the intercomplex separation of ∼ 1 kpc.     

The onset of warps in Spitzer observations of edge-on spiral galaxies

We analyse warps in the nearby edge-on spiral galaxies observed in the Spitzer /Infrared Array Camera (IRAC) 4.5-μm band. In our sample of 24 galaxies, we find evidence of warp in 14 galaxies. We estimate the observed onset radii for the warps in a subsample of 10 galaxies. The dark matter distribution in each of these galaxies are calculated using the mass distribution derived from the observed light distribution and the observed rotation curves. The theoretical predictions of the onset radii for the warps are then derived by applying a self-consistent linear response theory to the obtained mass models for six galaxies with rotation curves in the literature. By comparing the observed onset radii to the theoretical ones, we find that discs with constant thickness can not explain the observations; moderately flaring discs are needed. The required flaring is consistent with the observations. Our analysis shows that the onset of warp is not symmetric in our sample of galaxies. We define a new quantity called the onset-asymmetry index and study its dependence on galaxy properties. The onset asymmetries in warps tend to be larger in galaxies with smaller disc scalelengths. We also define and quantify the global asymmetry in the stellar light distribution, that we call the edge-on asymmetry in edge-on galaxies. It is shown that in most cases the onset asymmetry in warp is actually anticorrelated with the measured edge-on asymmetry in our sample of edge-on galaxies and this could plausibly indicate that the surrounding dark matter distribution is asymmetric.     

Global lopsided instability in a purely stellar galactic disc

It is shown that pure exponential discs in spiral galaxies are capable of supporting slowly varying discrete global lopsided modes, which can explain the observed features of lopsidedness in the stellar discs. Using linearized fluid dynamical equations with the softened self-gravity and pressure of the perturbation as the collective effect, we derive self-consistently a quadratic eigenvalue equation for the lopsided perturbation in the galactic disc. On solving this, we find that the ground-state mode shows the observed characteristics of the lopsidedness in a galactic disc, namely the fractional Fourier amplitude A ₁, increases smoothly with the radius. These lopsided patterns precess in the disc with a very slow pattern speed with no preferred sense of precession. We show that the lopsided modes in the stellar disc are long-lived because of a substantial reduction (approximately a factor of 10 compared to the local free precession rate) in the differential precession. The numerical solution of the equations shows that the ground-state lopsided modes are either very slowly precessing stationary normal mode oscillations of the disc or growing modes with a slow growth rate depending on the relative importance of the collective effect of the self-gravity. N -body simulations are performed to test the spontaneous growth of lopsidedness in a pure stellar disc. Both approaches are then compared and interpreted in terms of long-lived global   m = 1  instabilities, with almost zero pattern speed.     

Self-gravitating warped discs around supermassive black holes

We consider warped equilibrium configurations for stellar and gaseous discs in the Keplerian force field of a supermassive black hole, assuming that the self-gravity of the disc provides the only acting torques. Modelling the disc as a collection of concentric circular rings and computing the torques in the non-linear regime, we show that stable, strongly warped precessing equilibria are possible. These solutions exist for a wide range of disc-to-black-hole mass ratios   M _d/ M _bh  , can span large warp angles of up to  ±∼120°  , have inner and outer boundaries, and extend over a radial range of a factor of typically two to four. These equilibrium configurations obey a scaling relation such that in good approximation     where     is the (retrograde) precession frequency and Ω is a characteristic orbital frequency in the disc. Stability was determined using linear perturbation theory and, in a few cases, confirmed by numerical integration of the equations of motion. Most of the precessing equilibria are found to be stable, but some are unstable. The main result of this study is that highly warped discs near black holes can persist for long times without any persistent forcing other than by their self-gravity. The possible relevance of this to galactic nuclei is briefly discussed.     

Spatial and kinematical lopsidedness of atomic hydrogen in the Ursa Major group of galaxies

We have carried out the harmonic analysis of the atomic hydrogen (H  i ) surface density maps and the velocity fields for 11 galaxies belonging to the Ursa Major group, over a radial range of 4–6 disc scalelengths in each galaxy. This analysis gives the radial variation of spatial lopsidedness, quantified by the Fourier amplitude A ₁ of the   m = 1  component normalized to the average value. The kinematical analysis gives a value for the elongation of the potential to be ∼10 per cent. The mean amplitude of spatial lopsidedness is found to be ∼0.14 in the inner disc, similar to the field galaxies, and is smaller by a factor of ∼2 compared to the Eridanus group galaxies. It is also shown that the average value of A ₁ does not increase with the Hubble type, contrary to what is seen in field galaxies. We argue that the physical origin of lopsidedness in the Ursa Major group of galaxies is tidal interactions, albeit weaker and less frequent than in Eridanus. Thus systematic studies of lopsidedness in groups of galaxies can provide dynamical clues regarding the interactions and evolution of galaxies in a group environment.     

Submillimetre imaging of NGC 3079

Submillimetre mapping observations of the active edge-on spiral galaxy NGC 3079 are presented. These maps at 850 and 450 μm were made with the Submillimetre Common User Bolometer Array (SCUBA) at the James Clerk Maxwell Telescope (JCMT).
The source structure at these wavelengths consists of a central unresolved source embedded in diffuse disc emission, similar to that displayed at 1.2 mm. The disc emission is fitted with two optically thin, isothermal dust models which give temperatures of 12 and 31 K, similar to those derived previously by Braine et al. The core component is well described by a single-temperature fit (∼32 K). The combined dust mass from these observations, using the same mass absorption coefficient as Devereux & Young (1990) is 3.5×10⁸ M_⊙, of which ∼90 per cent resides in the cold component of the galactic disc. The effect of the cold dust component detected by SCUBA is thus to reduce the global gas-to-dust mass ratio from ∼1400 found in the above study to 85, very similar to the Galactic level. Calculations using the models of Draine & Lee and/or alternative molecular gas mass estimates yield gas-to-dust mass ratios in the range 60–190.
The data presented here, together with previously published 1.2-mm mapping observations and IRAS data, are inconsistent with detections made with the Infrared Space Observatory ( ISO ). In particular, the latter give an excess of flux at 200 and 180 μm relative to that predicted by our simple model fits (approximately a factor of 2–3).     

Structures in a class of magnetized scale-free discs

We construct analytically stationary global configurations for both aligned and logarithmic spiral coplanar magnetohydrodynamics (MHD) perturbations in an axisymmetric background MHD disc with a power-law surface mass density  Σ₀∝ r ^−α  , a coplanar azimuthal magnetic field   B ₀∝ r ^−γ  , a consistent self-gravity and a power-law rotation curve   v ₀∝ r ^−β  , where v ₀ is the linear azimuthal gas rotation speed. The barotropic equation of state  Π∝Σ ⁿ   is adopted for both MHD background equilibrium and coplanar MHD perturbations where Π is the vertically integrated pressure and n is the barotropic index. For a scale-free background MHD equilibrium, a relation exists among  α, β, γ  and n such that only one parameter (e.g. β) is independent. For a linear axisymmetric stability analysis, we provide global criteria in various parameter regimes. For non-axisymmetric aligned and logarithmic spiral cases, two branches of perturbation modes (i.e. fast and slow MHD density waves) can be derived once β is specified. To complement the magnetized singular isothermal disc analysis of Lou, we extend the analysis to a wider range of  −1/4 < β < 1/2  . As an illustrative example, we discuss specifically the  β= 1/4  case when the background magnetic field is force-free. Angular momentum conservation for coplanar MHD perturbations and other relevant aspects of our approach are discussed.     

An explanation for long flares from extragalactic globular cluster X-ray sources

We compare orbits in a thin axisymmetric disc potential in Modified Newtonian Dynamics (MOND) with those in a thin disc plus near-spherical dark matter halo predicted by a ΛCDM cosmology. Remarkably, the amount of orbital precession in MOND is nearly identical to that which occurs in a mildly oblate CDM Galactic halo (potential flattening   q = 0.9  ), consistent with recent constraints from the Sagittarius stream. Since very flattened mass distributions in MOND produce rounder potentials than in standard Newtonian mechanics, we show that it will be very difficult to use the tidal debris from streams to distinguish between a MOND galaxy and a standard CDM galaxy with a mildly oblate halo.
If a galaxy can be found with either a prolate halo or one that is more oblate than   q ∼ 0.9  this would rule out MOND as a viable theory. Improved data from the leading arm of the Sagittarius dwarf – which samples the Galactic potential at large radii – could rule out MOND if the orbital pole precession can be determined to an accuracy of the order of  ±1°  .     

Origin of disc lopsidedness in the Eridanus group of galaxies

The H  i surface density maps for a sample of 18 galaxies in the Eridanus group are Fourier analysed. This analysis gives the radial variation of the lopsidedness in the H  i spatial distribution. The lopsidedness is quantified by the Fourier amplitude A ₁ of the m = 1 component normalized to the average value. It is also shown that in the radial region where the stellar disc and H  i overlap, their A ₁ coefficients are comparable. All the galaxies studied show significant lopsidedness in H  i . The mean value of A ₁ in the inner regions of the galaxies (1.5–2.5 scalelengths) is ≥ 0.2. This value of A ₁ is twice the average value seen in the field galaxies. Also, the lopsidedness is found to be smaller for late-type galaxies; this is opposite to the trend seen in the field galaxies. These two results indicate a different physical origin for disc lopsidedness in galaxies in a group environment compared to the field galaxies. Further, a large fraction (∼30 per cent) shows a higher degree of lopsidedness ( A ₁≥ 0.3). It is also seen that the disc lopsidedness increases with the radius as demonstrated in earlier studies, but over a radial range that is two times larger than done in the previous studies. The average lopsidedness of the halo potential is estimated to be ∼10 per cent, assuming that the lopsidedness in H  i disc is due to its response to the halo asymmetry.     

A new superwind galaxy: XMM–Newton observations of NGC 6810

We present the first imaging X-ray observation of the highly inclined  ( i = 78°)  Sab Seyfert 2 galaxy NGC 6810 using XMM–Newton , which reveals soft X-ray emission that extends out to a projected height of ∼7 kpc away from the plane of the galaxy. The soft X-ray emission beyond the optical disc of the galaxy is most plausibly extraplanar, although it could instead come from large galactic radius. This extended X-ray emission is spatially associated with diffuse Hα emission, in particular with a prominent 5-kpc-long Hα filament on the north-west of the disc. A fraction ≲35 per cent of the total soft X-ray emission of the galaxy arises from projected heights  | z | ≥ 2 kpc  . Within the optical disc of the galaxy the soft X-ray emission is associated with the star-forming regions visible in ground-based Hα and XMM–Newton optical monitor near-UV imaging. The temperature, supersolar α-element-to-iron abundance ratio, soft X-ray/Hα correlation, and X-ray to far-infrared (FIR) flux ratio of NGC 6810 are all consistent with local starbursts with winds, although the large base radius of the outflow would make NGC 6810 one of the few 'disc-wide' superwinds currently known. Hard X-ray emission from NGC 6810 is weak, and the total   E = 2–10 keV  luminosity and spectral shape are consistent with the expected level of X-ray binary emission from the old and young stellar populations. The X-ray observations provide no evidence of any active galactic nucleus activity. We find that the optical, IR and radio properties of NGC 6810 are all consistent with a starburst galaxy, and that the old classification of this galaxy as a Seyfert 2 galaxy is probably incorrect.     

Inclination- and dust-corrected galaxy parameters: bulge-to-disc ratios and size–luminosity relations

While galactic bulges may contain no significant dust of their own, the dust within galaxy discs can strongly attenuate the light from their embedded bulges. Furthermore, such dust inhibits the ability of observationally determined inclination corrections to recover intrinsic (i.e. dust-free) galaxy parameters. Using the sophisticated 3D radiative transfer model of Popescu et al. and Tuffs et al., together with the recent determination of the average face-on opacity by Driver et al. in nearby disc galaxies, we provide simple equations to correct (observed) disc central surface brightness and scalelengths for the effects of both inclination and dust in the B , V , I , J and K passbands. We then collate and homogenize various literature data sets and determine the typical intrinsic scalelengths, central surface brightness and magnitudes of galaxy discs as a function of morphological type. All galaxies have been carefully modelled in their respective papers with a Sérsic   R ^{1/ n}   bulge plus an exponential disc. Using the bulge magnitude corrections from Driver et al., we additionally derive the average, dust-corrected, bulge-to-disc flux ratio as a function of galaxy type. With values typically less than 1/3, this places somewhat uncomfortable constraints on some current semi-analytic simulations. Typical bulge sizes, profile shapes, surface brightness and deprojected densities are provided. Finally, given the two-component nature of disc galaxies, we present luminosity–size and (surface brightness)–size diagrams for discs and bulges. We also show that the distribution of elliptical galaxies in the luminosity–size diagram is not linear but strongly curved.     

Do high-velocity clouds form by thermal instability?

We examine the proposal that the H  i 'high-velocity' clouds (HVCs) surrounding the Milky Way and other disc galaxies form by condensation of the hot galactic corona via thermal instability. Under the assumption that the galactic corona is well represented by a non-rotating, stratified atmosphere, we find that for this formation mechanism to work the corona must have an almost perfectly flat entropy profile. In all other cases, the growth of thermal perturbations is suppressed by a combination of buoyancy and thermal conduction. Even if the entropy profile were nearly flat, cold clouds with sizes smaller than  10 kpc  could form in the corona of the Milky Way only at radii larger than  100 kpc  , in contradiction with the determined distances of the largest HVC complexes. Clouds with sizes of a few kpc can form in the inner halo only in low-mass systems. We conclude that unless even slow rotation qualitatively changes the dynamics of a corona, thermal instability is unlikely to be a viable mechanism for formation of cold clouds around disc galaxies.     

The K-band Hubble diagram of submillimetre galaxies and hyperluminous galaxies

We present the K -band Hubble diagrams ( K – z relations) of submillimetre-selected galaxies and hyperluminous galaxies (HLIRGs). We report the discovery of a remarkably tight K – z relation of HLIRGs, indistinguishable from that of the most luminous radio galaxies. Like radio galaxies, the HLIRG K – z relation at   z ≲ 3  is consistent with a passively evolving ∼3 L _* instantaneous starburst starting from a redshift of   z ∼ 10  . In contrast, many submillimetre-selected galaxies are ≳2 mag fainter, and the population has a much larger dispersion. We argue that dust obscuration and/or a larger mass range may be responsible for this scatter. The galaxies so far proved to be hyperluminous may have been biased towards higher AGN bolometric contributions than submillimetre-selected galaxies due to the 60-μm selection of some, so the location on the K – z relation may be related to the presence of the most massive active galactic nucleus. Alternatively, a particular host galaxy mass range may be responsible for both extreme star formation and the most massive active nuclei.     

Stability of power-law discs — II. The global spiral modes

This paper reports on the in-plane normal modes in the self-consistent and the cut-out power-law discs. Although the cut-out discs are remarkably stable to bisymmetric perturbations, they are very susceptible to one-armed modes. For this harmonic, there is no inner Lindblad resonance, thus removing a powerful stabilizing influence. A physical mechanism for the generation of the one-armed instabilities is put forward. Incoming trailing waves are reflected as leading waves at the inner cut-out, thus completing the feedback for the swing-amplifier. Growing three-armed and four-armed modes occur only at very low temperatures. However, neutral m  = 3 and m  = 4 modes are possible at higher temperatures for some discs. The rotation curve index β has a marked effect on stability. For all azimuthal wavenumbers, any unstable modes persist to higher temperatures and grow more vigorously if the rotation curve is rising (β < 0) than if the rotation curve is falling (β > 0). If the central regions or outer parts of the disc are carved out more abruptly, any instabilities become more virulent. The self-consistent power-law discs possess a number of unusual stability properties. There is no natural time-scale in the self-consistent disc. If a mode is admitted at some pattern speed and growth rate, then it must be present at all pattern speeds and growth rates. Our analysis — although falling short of a complete proof — suggests that such a two-dimensional continuum of non-axisymmetric modes does not occur and that the self-consistent power-law discs admit no global non-axisymmetric normal modes whatsoever. Without reflecting boundaries or cut-outs, there is no resonant cavity and no possibility of unstable growing modes. The self-consistent power-law discs certainly admit equi-angular spirals as neutral modes, together with a one-dimensional continuum of growing axisymmetric modes.     

On the nature of superoutbursts in dwarf novae

We determine a crucial feature of the dark halo density distribution from the fact that the luminous matter dominates the gravitational potential at about one disc scalelength R _d, but at the optical edge     the dark matter has already become the main component of the galaxy density. From the kinematics of 137 spirals we find that the dark matter halo density profiles are self-similar at least out to R _opt and show core radii much larger than the corresponding disc scalelengths. The luminous regions of spirals consist of stellar discs embedded in dark haloes with roughly constant density. This invariant dark matter profile is very difficult to reconcile with the fundamental properties of the density distribution of cold dark matter haloes. With respect to previous work, the present evidence is obtained by means of a robust method and for a large and complete sample of normal spirals.     

Time-dependent models of the structure and stability of self-gravitating protoplanetary discs

Angular momentum transport within young massive protoplanetary discs may be dominated by self-gravity at radii where the disc is too weakly ionized to allow the development of the magneto-rotational instability. We use time-dependent one-dimensional disc models, based on a local cooling time calculation of the efficiency of transport, to study the radial structure and stability (against fragmentation) of protoplanetary discs in which self-gravity is the sole transport mechanism. We find that self-gravitating discs rapidly attain a quasi-steady state in which the surface density in the inner disc is high and the strength of turbulence very low (  α∼ 10⁻³  or less inside 5 au). Temperatures high enough to form crystalline silicates may extend out to several astronomical units at early times within these discs. None of our discs spontaneously develop regions that would be unambiguously unstable to fragmentation into substellar objects, though the outer regions (beyond 20 au) of the most massive discs are close enough to the threshold that fragmentation cannot be ruled out. We discuss how the mass accretion rates through such discs may vary with disc mass and with mass of the central star, and note that a determination of the     relation for very young systems may allow a test of the model.     

Light curves of jetted gamma-ray burst afterglows in circumstellar clouds

We present a Chandra observation of the powerful radio galaxy 3C 294 showing clear evidence for a surrounding intracluster medium. At a redshift of 1.786 this is the most distant cluster of galaxies yet detected in X-rays. The radio core is detected as a point source, which has a spectrum consistent with a heavily absorbed power law, implying an intrinsic 2–10 keV luminosity of ∼10⁴⁵ erg s⁻¹. A small excess of emission is associated with the southern radio hotspots. The soft, diffuse emission from the intracluster medium is centred on the radio source. It has an hourglass shape in the north–south direction, extending to radii of at least 100 kpc, well beyond the radio source. The X-ray spectrum of this extended component is fitted by a thermal model with temperature ∼5 keV, or by gas cooling from above 7 keV at rates of ∼ 400–700 M_⊙ yr⁻¹. The rest-frame 0.3–10 keV luminosity of the cluster is ∼ 4.5×10⁴⁴ erg s⁻¹. The existence of such a cluster is consistent with a low-density universe.     

Time evolution of galactic warps in prolate haloes

A recent observation with the Hipparcos satellite and some numerical simulations imply that the interaction between an oblate halo and a disc is inappropriate for the persistence of galactic warps. Following on from this , we have compared the time evolution of galactic warps in a prolate halo with that in an oblate halo. The haloes were approximated as fixed potentials, while the discs were represented by N -body particles. We have found that the warping in the oblate halo continues to wind up, and finally disappears. On the other hand, for the prolate halo model, the precession rate of the outer disc increases when the precession of the outer disc recedes from that of the inner disc, and vice versa. Consequently, the warping in the prolate halo persisted to the end of the simulation by retaining the alignment of the line of nodes of the warped disc. Therefore, our results suggest that prolate haloes could sustain galactic warps. The physical mechanism of the persistence of warp is discussed on the basis of the torque between a halo and a disc and that between the inner and outer regions of the disc.     

Exploring the active galactic nucleus and starburst content of local ultraluminous infrared galaxies through 5–8 μm spectroscopy

We present a 5–8 μm analysis of the Spitzer spectra of 71 ultraluminous infrared galaxies (ULIRGs) with redshift   z < 0.15  , devoted to the study of the role of active galactic nuclei (AGN) and starbursts (SBs) as the power source of the extreme infrared emission. Around ∼5 μm, an AGN is much brighter (by a factor of ≈30) than an SB of equal bolometric luminosity. This allows us to detect the presence of even faint accretion-driven cores inside ULIRGs: signatures of AGN activity are found in ∼70 per cent of our sample (50/71 sources). Through a simple analytical model, we are also able to obtain a quantitative estimate of the AGN/SB contribution to the overall energy output of each source. Although the main fraction of ULIRG luminosity is confirmed to arise from star formation events, the AGN contribution is non-negligible (∼23 per cent) and is shown to increase with luminosity. The existence of a rather heterogeneous pattern in the composition and geometrical structure of the dust among ULIRGs is newly supported by the comparison between individual absorption features and continuum extinction.