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.     

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 Cosmology of Extra Dimensions and Varying Fundamental Constants

I briefly present the Organizing Committee's and my own motivation for organizing this workshop, and I suggest a few key questions for which we will try to find possible answers in the coming days. This revised version was published online in July 2006 with corrections to the Cover Date.     

Automated analysis of eclipsing binary light curves – II. Statistical analysis of OGLE LMC eclipsing binaries

In the first paper of this series, we presented EBAS – Eclipsing Binary Automated Solver, a new fully automated algorithm to analyse the light curves of eclipsing binaries, based on the ebop code. Here, we apply the new algorithm to the whole sample of 2580 binaries found in the Optical Gravitational Lensing Experiment (OGLE) Large Magellanic Cloud (LMC) photometric survey and derive the orbital elements for 1931 systems. To obtain the statistical properties of the short-period binaries of the LMC, we construct a well-defined subsample of 938 eclipsing binaries with main-sequence B-type primaries. Correcting for observational selection effects, we derive the distributions of the fractional radii of the two components and their sum, the brightness ratios and the periods of the short-period binaries. Somewhat surprisingly, the results are consistent with a flat distribution in log P between 2 and 10 d. We also estimate the total number of binaries in the LMC with the same characteristics, and not only the eclipsing binaries, to be about 5000. This figure leads us to suggest that  (0.7 ± 0.4)  per cent of the main-sequence B-type stars in the LMC are found in binaries with periods shorter than 10 d. This frequency is substantially smaller than the fraction of binaries found by small Galactic radial-velocity surveys of B stars. On the other hand, the binary frequency found by Hubble Space Telescope ( HST ) photometric searches within the late main-sequence stars of 47 Tuc is only slightly higher and still consistent with the frequency we deduced for the B stars in the LMC.     

LU Vel (GJ 375): A M3.5Ve Flare and Double-Lined Spectroscopic Binary

High resolution echelle spectroscopic observations taken with the FEROS spectrograph at the 2.2 m telescope ESO confirm the binary nature of the flare M3.5V star LU Vel (GJ 375, RE J0958-462) previously reported by Christian and Mathioudakis (2002). Emission of similar intensity from both components is detected in the Balmer, Na i D₁&D₂, He i D₃, Ca ii H&K, and Ca ii IRT lines. We have determined precise radial velocities by cross correlation with radial velocity standard stars, which have allowed us to obtain for the first time the orbital solution of the system. The binary consists of two near-equal M3.5V components with an orbital period shorter than 2 days. We have analyzed the behaviour of the chromospheric activity indicators (variability and possible flares). In addition, we have determined its rotational velocity and kinematics.     

A seasonal climate model of the atmospheres of the giant planets at the voyager encounter time: I. Saturn's stratosphere

A radiative seasonal model which incorporates a multilayer radiative transfer treatment at wave-lengths longward of 7 μm is presented and applied to Saturn's stratosphere. Opacities due to H₂-He, CH₄, C₂H₂, and C₂H₆ are included. Season-dependent insolation is shown to produce a strong hemispheric asymmetry decreasing with depth at the Voyager encounter times, and seasonal amplitudes of 30°K at the poles are predicted in the high stratosphere. The ring-modulated dependence of the insolation and the orbital eccentricity are shown to have a significant effect. Calculations agree closely with the Voyager 1 and 2 radio occultation ingress profiles recorded at 76°S and 36.5°S for CH₄/H₂ = 3.5 + 1.4/? 1.0 × 10^?3;the estimated errors include modeling systematic errors and uncertainties in the occultations profiles. The possible role of aerosols in the stratospheric heating is analyzed. The Voyager 2 egress profile recorded at 31°S cannot be reproduced by calculations. Some constraints on the C₂H₂ and C₂H₆ abundances are derived. The upper portion of the occultation profiles (p < 3mbar) can be matched for C₂H₂/H₂ = 1.0 + 1.3/?0.6 × 10^?7, C₂H₆/H₂ = 1.5 + 1.8/?0.9 × 10^?6 at 76°S and C₂H₂/H₂ = 4 + 6/?4 × 10^?8, C₂H₆/H₂ = 6 + 9/?6 × 10^?7 at 36.5°N. At the northern occultation latitude, the discrepancy with the concentrations derived from analysis of IRIS spectra by R. Courtin, D. Gautier, A. Marten, B. Bézard, and R. Hanel (1984, Astrophys. J.287) can be explained by a sharp variation of the mixing ratios of these gases with altitude in the upper stratosphere. Other interpretations are discussed.