Simulations of multimode bolometric interferometers

We describe a procedure for the numerical modelling of astronomical interferometers, with particular relevance to far-infrared and submillimetre wavelengths. The scheme is based on identifying a set of modes that carry power from the sky to the detector. The procedure is extremely general, and can be used to model scalar or vector fields, in any state of coherence and polarization, the only limitation being that the propagation of a coherent field through the system be described by an integral transform, a constraint that is in practise always met.
We present simulations of ideal, multimode two-dimensional interferometers, and show that the modal theory reproduces the correct behaviour of both Michelson and Fizeau interferometers. We calculate simulated visibility data for a multimode bolometric Michelson interferometer, with a synthesized source, and produce a dirty map, recovering the original source with the usual artefacts associated with interferometers.     

EUV and coronagraphic observations of coronal mass ejections

The Large Angle Spectrometric Coronagraph (LASCO) and Extreme-ultraviolet Imaging Telescope (EIT) onboard Solar and Heliospheric Observatory (SOHO) provide us with unprecedented multi-wavelength observations helping us to understand different dynamic phenomena on the Sun and in the corona. In this paper we discuss the association between post-eruptive arcades (PEAs) detected by EIT and white-light coronal mass ejections (CMEs) detected by LASCO/C2 telescope.     

Self-consistent response of a galactic disc to vertical perturbations

We study the self-consistent, linear response of a galactic disc to vertical perturbations, as induced, say, by a tidal interaction. We calculate the self-gravitational potential corresponding to a non-axisymmetric, self-consistent density response of the disc using the Green's function method. The response potential is shown to oppose the perturbation potential because the self-gravity of the disc resists the imposed potential, and this resistance is stronger in the inner parts of a galactic disc. For the   m = 1  azimuthal wavenumber, the disc response opposes the imposed perturbation up to a radius that spans a range of 4–6 disc scalelengths, so that the disc shows a net warp only beyond this region. This physically explains the well known but so far unexplained observation that warps typically set in beyond this range of radii. We show that the inclusion of a dark matter halo in the calculation only marginally changes (by ∼10 per cent) the radius for the onset of warps. For perturbations with higher azimuthal wavenumbers, the net signature of the vertical perturbations can only be seen at larger radii – for example, beyond 7 exponential disc scalelengths for   m = 10  . Also, for the high- m cases, the magnitude of the negative disc response due to the disc self-gravity is much smaller. This is shown to result in corrugations of the mid-plane density, which explains the puzzling scalloping with   m = 10  detected in H  i in the outermost regions ∼30 kpc in the Galaxy.     

Thick and thin models of the evolution of carbon dioxide on Mars

We present results from a new simulation code that accounts for the evolution of the reservoirs of carbon dioxide on Mars, from its early years to the present. We establish a baseline model parameter set that produces results compatible with the present (i.e., Patm?6.5 mbar with permanent CO₂ ice cap) for a wide range of initial inventories. We find that the initial inventory of CO₂ broadly determines the evolutionary course of the reservoirs of CO₂. The reservoirs include the atmosphere, ice cap, adsorbed CO₂ in the regolith, and carbonate rocks. We track the evolution of the free inventory: the atmosphere, ice cap and regolith. Simulations begin at 4.53 Gyr before present with a rapid loss of free inventory to space in the early Noachian. Models that assume a relatively small initial inventory (?5 bar) have pronounced minima in the free inventory of CO₂ toward the end of the Noachian. Under baseline parameters, initial inventories below ∼4.5 bar result in a catastrophic loss of the free inventory to space. The current free inventory would be then determined by the balance between outgassing, sputtering losses and chemical weathering following the end of the late bombardment. We call these “thin” models. They generically predict small current free inventories in line with expectations of a small present CO₂ ice cap. For “thick” models, with initial inventories ?5 bar, a surplus of 300-700 mbar of free CO₂ remains during the late-Noachian. The histories of free inventory in time for thick models tend to converge within the last 3.5 Gyr toward a present with an ice cap plus atmospheric inventory of about 100 mbar. For thick models, the convergence is largely due to the effects of chemical weathering, which draws down higher free inventories more rapidly than the low. Thus, thick models have ?450 mbar carbonate reservoirs, while thin models have ?200 mbar. Though both thick and thin scenarios can reproduce the current atmospheric pressure, the thick models imply a relatively large current CO₂ ice cap and thin models, little or none. While the sublimation of a massive cap at a high obliquity would create a climate swing of greenhouse warming for thick models, under the thin model, mean temperatures and pressures would be essentially unaffected by increases in obliquity.     

The Relation between the Amplitude and the Period of Solar Cycles

The maximum amplitudes of solar activity cycles are found to be well anti-correlated (r = -0.72) with the newly defined solar cycle lengths three cycles before (at lag -3) in 13-month running mean sunspot numbers during the past 190 years. This result could be used for predicting the maximum sunspot numbers. The amplitudes of Cycles 24 and 25 are estimated to be 149.5±27.6 and 144.3±27.6, respectively.     

Disk galaxy evolution: from the Milky Way to high-redshift disks

We develop a detailed model of the Milky Way (a `prototypical' disk galaxy) and extend it to other disks with the help of some simple scaling relations, obtained in the framework of Cold Dark Matter models. This phenomenological (`hybrid') approach to the study of disk galaxy evolution allows us to reproduce successfully a large number of observed properties of disk galaxies in the local Universe and up to redshift z ∼ 1. The important conclusion is that, on average, massive disks have formed the bulk of their stars earlier than their lower mass counterparts: the `star formation hierarchy' has been apparently opposite to the `dark matter assembly' hierarchy. It is not yet clear whether `feedback' (as used in semi-analytical models of galaxy evolution) can explain that discrepancy. This revised version was published online in August 2006 with corrections to the Cover Date.