Cosmography with cluster strong lensing

By stacking an ensemble of strong lensing clusters, we demonstrate the feasibility of placing constraints on the dark energy equation of state. This is achieved by using multiple images of sources at two or more distinct redshift planes. The sample of smooth clusters in our simulations is based on observations of massive clusters and the distribution of background galaxies is constructed using the Hubble Deep Field . Our source distribution reproduces the observed redshift distribution of multiply imaged sources in Abell 1689. The cosmology recovery depends on the number of image families with known spectroscopic redshifts and the number of stacked clusters. Our simulations suggest that constraints comparable to those derived from other competing established techniques on a constant dark energy equation of state can be obtained using 10–40 clusters with five or more families of multiple images. We have also studied the observational errors in the image redshifts and positions. We find that spectroscopic redshifts and high-resolution Hubble Space Telescope ( HST ) images are required to eliminate confidence contour relaxation relative to the ideal case in our simulations. This suggests that the dark energy equation of state, and other cosmological parameters, can be constrained with existing HST images of lensing clusters coupled with dedicated ground-based arc spectroscopy.     

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.     

Annual cycles in the interstellar scintillation time-scales of PKS B1519–273 and PKS B1622–253

We have used the University of Tasmania's 30-m radio telescope at Ceduna in South Australia to regularly monitor the flux density of a number of southern blazars. We report the detection of an annual cycle in the variability time-scale of the centimetre radio emission of PKS B1622−253. Observations of PKS B1519−273 over a period of nearly 2 yr confirm the presence of an annual cycle in the variability time-scale in that source. These observations prove that interstellar scintillation is the principal cause of inter-day variability at radio wavelengths in these sources. The best-fitting annual cycle model for both sources implies a high degree of anisotropy in the scattering screen and that it has a large velocity offset with respect to the local standard of rest. This is consistent with a greater screen distance for these 'slow' intra-day variability (IDV) sources than for rapid scintillators such as PKS B0405−385 or J1819+3845.     

Optical and near-infrared colours as a discriminant of the age and metallicity of stellar populations

We present a comprehensive analysis of the ability of current stellar population models to reproduce the optical ( ugriz ) and near-infrared ( JHK ) colours of a small sample of well-studied nearby elliptical and S0 galaxies. We find broad agreement between the ages and metallicities derived using different population models, although different models show different systematic deviations from the measured broad-band fluxes. Although it is possible to constrain simple stellar population models to a well-defined area in age–metallicity space, there is a clear degeneracy between these parameters even with such a full range of precise colours. The precision to which age and metallicity can be determined independently, using only broad-band photometry with realistic errors, is  Δ[Fe/H]≃ 0.18  and  Δlog Age ≃ 0.25  . To constrain the populations and therefore the star formation history further, it will be necessary to combine broad-band optical–IR photometry with either spectral line indices, or else photometry at wavelengths outside this range.     

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.     

On the origin of the 511-keV emission in the Galactic Centre

Diffuse 511-keV line emission, from the annihilation of cold positrons, has been observed in the direction of the Galactic Centre for more than 30 yr. The latest high-resolution maps of this emission produced by the SPI instrument on INTEGRAL suggest at least one component of the emission is spatially coincident with the distribution of ∼70 luminous, low-mass X-ray binaries detected in the soft gamma-ray band. The X-ray band, however, is generally a more sensitive probe of X-ray binary populations. Recent X-ray surveys of the Galactic Centre have discovered a much larger population (>4000) of faint, hard X-ray point sources. We investigate the possibility that the positrons observed in the direction of the Galactic Centre originate in pair-dominated jets generated by this population of fainter accretion-powered X-ray binaries. We also consider briefly whether such sources could account for unexplained diffuse emission associated with the Galactic Centre in the microwave (the Wilkinson Microwave Anisotropy Probe 'haze') and at other wavelengths. Finally, we point out several unresolved problems in associating Galactic Centre 511-keV emission with the brightest X-ray binaries.     

The sensitivity of the next generation of lunar Cherenkov observations to UHE neutrinos and cosmic rays

We present simulation results for the detection of ultra-high energy (UHE) cosmic ray (CR) and neutrino interactions in the Moon by radio-telescopes. We simulate the expected radio signal at Earth from such interactions, expanding on previous work to include interactions in the sub-regolith layer for single dish and multiple telescope systems. For previous experiments at Parkes, Goldstone (GLUE), and Kalyazin we recalculate the sensitivity to an isotropic flux of UHE neutrinos. We find the published sensitivity for the GLUE experiment to be too high (too optimistic) by an order of magnitude, and consequently the GLUE limit to be too low by an order of magnitude. Our predicted sensitivity for future experiments using the Australia Telescope Compact Array (ATCA) and the Australian SKA Pathfinder (ASKAP) indicate these instruments will be able to detect the more optimistic UHE neutrino flux predictions, while the square kilometre array (SKA) will also be sensitive to all bar one prediction of a diffuse ‘cosmogenic’, or ‘GZK’, neutrino flux.Outstanding theoretical uncertainties at both high-frequency and low-frequency limits currently prevent a reliable estimate of the sensitivity of the lunar Cherenkov technique for UHE cosmic ray (CR) astronomy. Here, we place limits on the effects of large-scale surface roughness on UHE CR detection, and find that when near-surface ‘formation-zone’ effects are ignored, the proposed SKA low-frequency aperture array could detect CR events above 56 EeV at a rate between 15 and 40 times that of the current Pierre Auger Observatory. Should further work indicate that formation-zone effects have little impact on UHE CR sensitivity, observations of the Moon with the SKA would allow directional analysis of UHE cosmic rays, and investigation of correlations with putative cosmic ray source populations, to be conducted with very high statistics.     

Magnetically-regulated fragmentation induced by nonlinear flows and ambipolar diffusion

We present a parameter study of simulations of fragmentation regulated by gravity, magnetic fields, ambipolar diffusion, and nonlinear flows. The thin-sheet approximation is employed with periodic lateral boundary conditions, and the nonlinear flow field (“turbulence”) is allowed to freely decay. In agreement with previous results in the literature, our results show that the onset of runaway collapse (formation of the first star) in subcritical clouds is significantly accelerated by nonlinear flows in which a large-scale wave mode dominates the power spectrum. In addition, we find that a power spectrum with equal energy on all scales also accelerates collapse, but by a lesser amount. For a highly super-Alfvénic initial velocity field with most power on the largest scales, the collapse occurs promptly during the initial compression wave. However, for trans-Alfvénic perturbations, a subcritical magnetic field causes a rebound from the initial compression, and the system undergoes several oscillations before runaway collapse occurs. Models that undergo prompt collapse have highly supersonic infall motions at the core boundaries. Cores in magnetically subcritical models with trans-Alfvénic initial perturbations also pick up significant systematic speeds by inheriting motions associated with magnetically-driven oscillations. Core mass distributions are much broader than in models with small-amplitude initial perturbations, although the disturbed structure of cores that form due to nonlinear flows does not guarantee subsequent monolithic collapse. Our simulations also demonstrate that significant power (if present initially) can be maintained with negligible dissipation in large-scale compressive modes of a magnetic thin sheet, in the limit of perfect flux freezing.     

A simple unsplit Godunov method for multidimensional MHD

We describe a numerical algorithm based on Godunov methods for integrating the equations of compressible magnetohydrodynamics (MHD) in multidimensions. It combines a simple, dimensionally-unsplit integration method with the constrained transport (CT) discretization of the induction equation to enforce the divergence-free constraint. We present the results of a series of fully three-dimensional tests which indicate the method is second-order accurate for smooth solutions in all MHD wave families, and captures shocks, contact and rotational discontinuities well. However, it is also more diffusive than other more complex unsplit integrators combined with CT. Thus, the primary advantage of the method is its simplicity. It does not require a characteristic tracing step to construct interface values for the Riemann solver, it is straightforward to extend with additional physics, and it is suitable for use with nested and adaptive meshes. The method is implemented as one of two dimensionally unsplit MHD integrators in the Athena code, which is freely available for download from the web.     

A multiphysics and multiscale software environment for modeling astrophysical systems

We present MUSE, a software framework for combining existing computational tools for different astrophysical domains into a single multiphysics, multiscale application. MUSE facilitates the coupling of existing codes written in different languages by providing inter-language tools and by specifying an interface between each module and the framework that represents a balance between generality and computational efficiency. This approach allows scientists to use combinations of codes to solve highly coupled problems without the need to write new codes for other domains or significantly alter their existing codes. MUSE currently incorporates the domains of stellar dynamics, stellar evolution and stellar hydrodynamics for studying generalized stellar systems. We have now reached a “Noah’s Ark” milestone, with (at least) two available numerical solvers for each domain. MUSE can treat multiscale and multiphysics systems in which the time- and size-scales are well separated, like simulating the evolution of planetary systems, small stellar associations, dense stellar clusters, galaxies and galactic nuclei. In this paper we describe three examples calculated using MUSE: the merger of two galaxies, the merger of two evolving stars, and a hybrid N-body simulation. In addition, we demonstrate an implementation of MUSE on a distributed computer which may also include special-purpose hardware, such as GRAPEs or GPUs, to accelerate computations. The current MUSE code base is publicly available as open source at http://muse.li.