Superbubbles in magnetized interstellar media: blowout or confinement?

The importance of the interstellar magnetic field is studied in relation to the evolution of superbubbles with a three-dimensional (3D) numerical magnetohydrodynamical (MHD) simulation. A superbubble is a large supernova remnant driven by sequential supernova explosions in an OB association. Its evolution is affected by the density stratification in the galactic disc. After the superbubble size reaches 2–3 times the density scaleheight, it expands preferentially in the z -direction, until finally it can punch out a hole in the gas disc (blowout). On the other hand, the magnetic field running parallel to the galactic disc has the effect of preventing it from expanding in the direction perpendicular to the field. The density stratification and the magnetic fields have completely opposite effects on the evolution of the superbubble. We present results of a 3D MHD simulation in which both effects are included. As a result, it is concluded that when the magnetic field has a much larger scaleheight than the density, even for a model in which the bubble would blow out from the disc if the magnetic field were absent, a magnetic field with a strength of 5 μG can confine the bubble in | z |≲300 pc for ≃ 20 Myr (confinement). In a model in which the field strength decreases in the halo as B  ∝ ρ^1/2, the superbubble eventually blows out like a model with B  = 0 even if the magnetic field in the mid-plane is as strong as B  = 5 μG.     

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.     

Large-scale magnetic fields: galaxy two-point correlation function

We study the effect of large-scale tangled magnetic fields on the galaxy two-point correlation function in the redshift space. We show that (i) the magnetic field effects can be comparable to the gravity-induced clustering for present magnetic field strength   B ₀≃ 5 × 10⁻⁸ G  , (ii) the absence of this signal from the present data gives an upper bound   B ₀≲ 3 × 10⁻⁸ G  and (iii) the future data can probe the magnetic fields of  ≃10⁻⁸ G  . A comparison with other constraints on the present magnetic field shows that they are marginally compatible. However, if the magnetic fields corresponding to   B ₀≃ 10⁻⁸ G  existed at the last scattering surface, they will cause unacceptably large cosmic microwave background radiation anisotropies.     

Formation of interstellar clouds: Parker instability with phase transitions

We numerically follow the nonlinear evolution of the Parker instability in the presence of phase transitions from a warm to a cold H  i interstellar medium in two spatial dimensions. The nonlinear evolution of the system favours modes that allow the magnetic field lines to cross the galactic plane. Cold H  i clouds form with typical masses  ≃10⁵ M_⊙  , mean densities  ≃20 cm⁻³  , mean magnetic-field strengths  ≃4.3 μG  (rms field strengths  ≃6.4 μG  ), mass-to-flux ratios  ≃0.1–0.3  relative to critical, temperatures  ≃50 K  , (two-dimensional) turbulent velocity dispersions  ≃1.6 km s⁻¹  and separations  ≃500 pc  , in agreement with observations. The maximum density and magnetic-field strength are  ≃10³ cm⁻³  and  ≃20 μG  , respectively. Approximately 60 per cent of all H  i mass is in the warm neutral medium. The cold neutral medium is arranged into sheet-like structures both perpendicular and parallel to the galactic plane, but it is also found almost everywhere in the galactic plane, with the density being highest in valleys of the magnetic field lines. 'Cloudlets' also form whose physical properties are in quantitative agreement with those observed for such objects by Heiles. The nonlinear phase of the evolution takes ≲30 Myr, so that, if the instability is triggered by a nonlinear perturbation such as a spiral density shock wave, interstellar clouds can form within a time suggested by observations.     

The galaxy luminosity–size relation and selection biases in the Hubble Ultra Deep Field

We use the Hubble Ultra Deep Field to study the galaxy luminosity–size  ( M – R _e )  distribution. With a careful analysis of selection effects due to both detection completeness and measurement reliability, we identify bias-free regions in the   M – R _e   plane for a series of volume-limited samples. By comparison to a nearby survey also having well-defined selection limits, namely the Millennium Galaxy Catalogue, we present clear evidence for evolution in surface brightness since   z ∼ 0.7  . Specifically, we demonstrate that the mean, rest-frame B -band  〈μ〉 _e   for galaxies in a sample spanning 8 mag in luminosity between   M _B =−22  and −14 mag increases by ∼1.0 mag arcsec⁻² from   z ∼ 0.1  to 0.7. We also highlight the importance of considering surface brightness-dependent measurement biases in addition to incompleteness biases. In particular, the increasing, systematic underestimation of Kron fluxes towards low surface brightnesses may cause diffuse, yet luminous, systems to be mistaken for faint, compact objects.     

Continuous frequency spectrum of the global hydromagnetic oscillations of a magnetically confined mountain on an accreting neutron star

We compute the continuous part of the ideal-magnetohydrodynamic (ideal-MHD) frequency spectrum of a polar mountain produced by magnetic burial on an accreting neutron star. Applying the formalism developed by Hellsten & Spies, extended to include gravity, we solve the singular eigenvalue problem subject to line-tying boundary conditions. This spectrum divides into an Alfvén part and a cusp part. The eigenfunctions are chirped and anharmonic with an exponential envelope, and the eigenfrequencies cover the whole spectrum above a minimum ω_low. For equilibria with accreted mass  1.2 × 10⁻⁶≲ M _a/M_⊙≲ 1.7 × 10⁻⁴  and surface magnetic fields  10¹¹≲ B _*/G ≲ 10¹³, ω_low  is approximately independent of   B _*  , and increases with M _a. The results are consistent with the Alfvén spectrum excited in numerical simulations with the zeus-mp solver. The spectrum is modified substantially by the Coriolis force in neutron stars spinning faster than ∼100 Hz. The implications for gravitational-wave searches for low-mass X-ray binaries are considered briefly.     

The evolution of galactic discs

We use recent observations of high-redshift galaxies to study the evolution of galactic discs over the redshift range 0 <  z ≲1. The data are inconsistent with models in which discs were already assembled at z  = 1 and have evolved only in luminosity since that time. Assuming that disc properties change with redshift as powers of 1 +   z and analysing the observations assuming an Einstein–de Sitter universe, we find that for given rotation speed, disc scalelength decreases with z as ∼ (1 +  z )⁻¹, total B -band mass-to-light ratio decreases with z as ∼ (1 +  z )⁻¹, and disc luminosity (again in B ) depends only weakly on z . These scalings are consistent with current data on the evolution of disc galaxy abundance as a function of size and luminosity. Both the scalings and the abundance evolution are close to the predictions of hierarchical models for galaxy formation. If different cosmogonies are compared, the observed evolution in disc size and disc abundance favours a flat low-Ω₀ universe over an Einstein–de Sitter universe.     

Mass models of the Milky Way

A parametrized model of the mass distribution within the Milky Way is fitted to the available observational constraints. The most important single parameter is the ratio of the scalelength R _d* of the stellar disc to R ₀. The disc and bulge dominate v _c( R ) at R ≲ R ₀ only for R _d,*/ R ₀≲0.3. Since the only knowledge we have of the halo derives from studies like the present one, we allow it to contribute to the density at all radii. When allowed this freedom, however, the halo causes changes in assumptions relating to R  ≪  R ₀ to affect profoundly the structure of the best-fitting model at R  ≫  R ₀. For example, changing the disc slightly from an exponential surface-density profile significantly changes the form of v _c( R ) at R  ≫  R ₀, where the disc makes a negligible contribution to v _c. Moreover, minor changes in the constraints can cause the halo to develop a deep hole at its centre that is not physically plausible. These problems call into question the proposition that flat rotation curves arise because galaxies have physically distinct haloes rather than outwards-increasing mass-to-light ratios.   The mass distribution of the Galaxy and the relative importance of its various components will remain very uncertain until more observational data can be used to constrain mass models. Data that constrain the Galactic force field at z ≳ R and at R  >  R ₀ are especially important.     

Radial scalelengths of the galactic thin and thick disc with 2MASS data

This paper presents a global analysis of the 2MASS (Two Micron All Sky Survey) data as observed in seven fields at different galactic latitudes in our Galaxy. The data allow the preliminary determination of the scale parameters, which lead to strong constraints on the radial and vertical structure of the galactic thin and thick disc. The interpretation of star counts and colour distributions of stars in the near-infrared with the synthetic stellar population model gives strong evidence that the galactic thin disc density scalelength ( h _R ) is rather short (2.8±0.3 kpc). The galactic thick disc population is revisited in the light of new data. We find the thick disc to have a local density of 3.5±2.0 per cent of the thin disc, exponential scaleheight ( h _z ) of 860±200 pc and exponential scalelength ( h _R ) of 3.7±_0.5^0.8 kpc.     

Chemical evolution with radial mixing

Models of the chemical evolution of our Galaxy are extended to include radial migration of stars and flow of gas through the disc. The models track the production of both iron and α-elements. A model is chosen that provides an excellent fit to the metallicity distribution of stars in the Geneva–Copenhagen survey (GCS) of the solar neighbourhood and a good fit to the local Hess diagram. The model provides a good fit to the distribution of GCS stars in the age–metallicity plane, although this plane was not used in the fitting process. Although this model's star formation rate is monotonically declining, its disc naturally splits into an α-enhanced thick disc and a normal thin disc. In particular, the model's distribution of stars in the ([O/Fe], [Fe/H]) plane resembles that of Galactic stars in displaying a ridge line for each disc. The thin-disc's ridge line is entirely due to stellar migration, and there is the characteristic variation of stellar angular momentum along it that has been noted by Haywood in survey data. Radial mixing of stellar populations with high  σ _z   from inner regions of the disc to the solar neighbourhood provides a natural explanation of why measurements yield a steeper increase of  σ _z   with age than predicted by theory. The metallicity gradient in the interstellar medium is predicted to be steeper than in earlier models, but appears to be in good agreement with data for both our Galaxy and external galaxies. The models are inconsistent with a cut-off in the star formation rate at low gas surface densities. The absolute magnitude of the disc is given as a function of time in several photometric bands, and radial colour profiles are plotted for representative times.     

Primordial magnetic fields and formation of molecular hydrogen

We study the implications of primordial magnetic fields for the thermal and ionization history of the post-recombination era. In particular, we compute the effects of dissipation of primordial magnetic fields owing to ambipolar diffusion and decaying turbulence in the intergalactic medium (IGM) and the collapsing haloes, and compute the effects of the altered thermal and ionization history on the formation of molecular hydrogen. We show that, for magnetic field strengths in the range  2 × 10⁻¹⁰≲ B ₀≲ 2 × 10⁻⁹ G  , the molecular hydrogen fraction in IGM and collapsing halo can increase by a factor of 5 to 1000 over the case with no magnetic fields. We discuss the implication of the increased molecular hydrogen fraction on the radiative transfer of ultraviolet photons and the formation of first structures in the universe.     

Molecular hydrogen in damped Lyα systems: clues to interstellar physics at high redshift

In order to interpret H₂ quasar absorption-line observations of damped Lyα systems (DLAs) and subDLAs, we model their H₂ abundance as a function of dust-to-gas ratio, including H₂ self-shielding and dust extinction against dissociating photons. Then, we constrain the physical state of the gas by using H₂ data. Using H₂ excitation data for DLAs with H₂ detections, we derive a gas density  1.5 ≲ log n (cm⁻³) ≲ 2.5  , temperature  1.5 ≲ log T (K) ≲ 3  , and an internal ultraviolet (UV) radiation field (in units of the Galactic value)  0.5 ≲ log χ≲ 1.5  . We then find that the observed relation between the molecular fraction and the dust-to-gas ratio of the sample is naturally explained by the above conditions. However, it is still possible that H₂ deficient DLAs and subDLAs with H₂ fractions less than  ∼10⁻⁶  are in a more diffuse and warmer state. The efficient photodissociation by the internal UV radiation field explains the extremely small H₂ fraction  (≲10⁻⁶)  observed for  κ≲ 1/30  (κ is the dust-to-gas ratio in units of the Galactic value); H₂ self-shielding causes a rapid increase in, and large variations of, H₂ abundance for  κ≳ 1/30  . We finally propose an independent method to estimate the star formation rates of DLAs from H₂ abundances; such rates are then critically compared with those derived from other proposed methods. The implications for the contribution of DLAs to the cosmic star formation history are briefly discussed.     

The influence of Hall drift on the magnetic fields of neutron stars

We consider the evolution of magnetic fields under the influence of Hall drift and Ohmic decay. The governing equation is solved numerically, in a spherical shell with   r _i / r _o = 0.75  . Starting with simple free-decay modes as initial conditions, we then consider the subsequent evolution. The Hall effect induces so-called helicoidal oscillations, in which energy is redistributed among the different modes. We find that the amplitude of these oscillations can be quite substantial, with some of the higher harmonics becoming comparable with the original field. Nevertheless, this transfer of energy to the higher harmonics is not sufficient to accelerate significantly the decay of the original field, at least not at the   R _B = O (100)  parameter values accessible to us, where this Hall parameter   R _B   measures the ratio of the Ohmic time-scale to the Hall time-scale. We do find clear evidence though of increasingly fine structures developing for increasingly large   R _B   , suggesting that perhaps this Hall-induced cascade to ever-shorter length-scales is eventually sufficiently vigorous to enhance the decay of the original field. Finally, the implications for the evolution of neutron star magnetic fields are discussed.     

High-energy emission from the starburst galaxy NGC 253

Measurement sensitivity in the energetic γ-ray region has improved considerably and is about to increase further in the near future, motivating a detailed calculation of high-energy (HE; ≥100 MeV) and very high-energy (VHE; ≥100 GeV) γ-ray emission from the nearby starburst galaxy NGC 253. Adopting the convection–diffusion model for energetic electron and proton propagation, and accounting for all the relevant hadronic and leptonic processes, we determine the steady-state energy distributions of these particles by a detailed numerical treatment. The electron distribution is directly normalized by the measured synchrotron radio emission from the central starburst region; a commonly expected theoretical relation is then used to normalize the proton spectrum in this region. Doing so fully specifies the electron spectrum throughout the galactic disc and, with an assumed spatial profile of the magnetic field, the predicted radio emission from the full disc matches well the observed spectrum, confirming the validity of our treatment. The resulting radiative yields of both particles are calculated; the integrated HE and VHE fluxes from the entire disc are predicted to be   f (≥100 MeV) ≃ (1.8^+1.5_−0.8) × 10⁻⁸ cm⁻² s⁻¹  and   f (≥100 GeV) ≃ (3.6^+3.4_−1.7) × 10⁻¹² cm⁻² s⁻¹  , with a central magnetic field value   B ₀≃ 190 ± 10 μ  G. We discuss the feasibility of measuring emission at these levels with the space-borne Fermi and ground-based Cherenkov telescopes.     

Seed magnetic Fields Generated by Primordial Supernova Explosions

The origin of the magnetic field in galaxies is an open question in astrophysics. Several mechanisms have been proposed related, in general, to the generation of small seed fields amplified by a dynamo mechanism. In general, these mechanisms have difficulty in satisfying both the requirements of a sufficiently high strength for the magnetic field and the necessary large coherent scales. We show that the formation of dense and turbulent shells of matter, in the multiple explosion scenario of Miranda &38; Opher for the formation of the large-scale structures of the Universe, can naturally act as a seed for the generation of a magnetic field. During the collapse and explosion of Population III objects, a temperature gradient not parallel to a density gradient can naturally be established, producing a seed magnetic field through the Biermann battery mechanism. We show that seed magnetic fields ∼ 10⁻¹²–10⁻¹⁴ G can be produced in this multiple explosion scenario on scales of the order of clusters of galaxies (with coherence length L  ∼ 1.8 Mpc) and up to ∼ 4.5 × 10⁻¹⁰ G on scales of galaxies ( L  ∼ 100 kpc).     

Polarization of synchrotron emission from relativistic reconfinement shocks

We study the polarization properties of relativistic reconfinement shocks with chaotic magnetic fields. Using our hydrodynamical model of their structure, we calculate synthetic polarization maps, longitudinal polarization profiles and discuss the spatially averaged polarization degree as a function of jet half-opening angle  Θ _j   , jet Lorentz factor  Γ _j   and observer inclination angle to the jet axis  θ_obs  . We find that for  θ_obs≲Θ _j   the wave electric vectors are parallel in the vicinity of the structure ends and perpendicular in between, while for  θ_obs > Θ _j   the polarization can only be perpendicular. The spatially averaged polarization degree does not exceed 30 per cent. Parallel average polarization, with polarization degrees lower than 10 per cent, has been found for  θ_obs < Θ _j   under the condition  Γ _j Θ _j > 1  . As earlier works predicted the parallel polarization from relativistic conical shocks, we explain our results by discussing conical shocks with divergent upstream flow.     

On the H i content, dust-to-gas ratio and nature of Mg ii absorbers

We estimate the mean dust-to-gas ratio of Mg  ii absorbers as a function of rest equivalent width W ₀ and redshift over the range  0.5 < z < 1.4  . Using the expanded Sloan Digital Sky Survey/ Hubble Space Telescope sample of low-redshift Lyman-α absorbers, we first show the existence of a 8σ correlation between the mean hydrogen column density  〈 N _{H  i} 〉  and W ₀, an indicator of gas velocity dispersion. By combining these results with recent dust-reddening measurements, we show that the mean dust-to-gas ratio of Mg  ii absorbers does not appreciably depend on rest equivalent width. Assuming that, on average, dust-to-gas ratio is proportional to metallicity, we find its redshift evolution to be consistent with that of   L ^★  galaxies from   z = 0.5  to 1.4, and we show that our constraints disfavour dwarf galaxies as the origin of such absorbers. We discuss other scenarios and favour galactic outflows from  ∼ L ^★  galaxies as the origin of the majority of strong Mg  ii absorbers. Finally, we show that, once evolutionary effects are taken into account, the Bohlin et al. relation between A _V and N _H is also satisfied by strong Mg  ii systems down to lower column densities than those probed in our Galaxy.     

RXTE observation of NGC 6240: a search for the obscured active nucleus

We present new data taken at 850 μm with SCUBA at the James Clerk Maxwell Telescope for a sample of 19 luminous infrared galaxies. Fourteen galaxies were detected. We have used these data, together with fluxes at 25, 60 and 100 μm from IRAS , to model the dust emission. We find that the emission from most galaxies can be described by an optically thin, single temperature dust model with an exponent of the dust extinction coefficient ( k _λ ∝ λ ^{− β} ) of β ≃1.4–2. A lower β ≃1 is required to model the dust emission from two of the galaxies, Arp 220 and NGC 4418. We discuss various possibilities for this difference and conclude that the most likely is a high dust opacity. In addition, we compare the molecular gas mass derived from the dust emission, M _850 μm, with the molecular gas mass derived from the CO emission, M _CO, and find that M _CO is on average a factor 2–3 higher than M _850 μm.     

Luminous hot accretion discs

We find a new two-temperature hot branch of equilibrium solutions for stationary accretion discs around black holes. In units of Eddington accretion rate defined as 10 L _Edd c ², the accretion rates to which these solutions correspond are within the range ̇ ₁≲ ̇ ≲1, where ̇ ₁ is the critical rate of advection-dominated accretion flow (ADAF). In these solutions, the energy loss rate of the ions by Coulomb energy transfer between the ions and electrons is larger than the viscously heating rate and it is the advective heating together with the viscous dissipation that balances the Coulomb cooling of ions. When ̇ ₁≲ ̇ ≲ ̇ ₂, where ̇ ₂∼5 ̇ ₁<1, the accretion flow remains hot throughout the disc. When ̇ ₂≲ ̇ ≲1, Coulomb interaction will cool the inner region of the disc within a certain radius ( r _tr∼several tens of Schwarzschild radii or larger depending on the accretion rate and the outer boundary condition) and the disc will collapse on to the equatorial plane and form an optically thick cold annulus. Compared with ADAF, these hot solutions are much more luminous because of the high accretion rate and efficiency; therefore, we call them luminous hot accretion discs.