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.     

Asymmetry of 56Ni ejecta and polarization in the type-IIP supernova 2004dj

I study the question of whether the asymmetry of ⁵⁶Ni ejecta that results in the asymmetry of the Hα emission line at the nebular epoch of the type-IIP supernova SN 2004dj can account for the recently detected polarization of the supernova radiation. I have developed a model of the Hα profile and luminosity with nonthermal ionization and excitation in a spherically symmetric envelope for an asymmetric bipolar ⁵⁶Ni distribution. I have calculated the polarized radiation transfer against the background of the recovered electron density distribution. The observed polarization is shown to be reproduced at the nebular epoch around day 140 for the same parameters of the envelope, and the ⁵⁶Ni distribution for which the evolution of the Hα luminosity and profile is explained. Yet the model polarization decreases with time more slowly than is observed. The origin of the additional component responsible for the early polarization on day 107 is discussed.     

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.     

A synthesis map of the sky at 34.5 MHz

This paper describes a wide-field survey made at 34.5 MHz using GEETEE,1 the low frequency telescope at Gauribidanur (latitude 13°36′12′′N). This telescope was used in the transit mode and by per forming 1-D synthesis along the north-south direction the entire observable sky was mapped in a single day. This minimized the problems that hinder wide-field low-frequency mapping. This survey covers the declination range of-50° to + 70° (- 33° to +61° without aliasing) and the complete 24 hours of right ascension. The synthesized beam has a resolution of 26′ x 42′ sec (δ- 14°. 1). The sensitivity of the survey is 5 Jy/beam (1σ). Special care has been taken to ensure that the antenna responds to all angular scale structures and is suitable for studies of both point sources and extended objects This telescope is jointly operated by the Indian Institute of Astrophysics, Bangalore and the Roman Research Institute, Bangalore.