Semi-annual report on a space background study. Contract DA-30-069-ORD-2803. EK/ARD ED-442 [ED-575, ED-692 and ED-735] For the Army Rocket and Guided Missile Agency, U.S. Army Ordnance Missile Command, Redstone Arsenal, Alabama.

Includes bibliography.

1. Astronomical photometry at wavelengths of 8.5, 10.5, and 11.6 µm. 2. Infrared emission from asteroids at wavelengths of 8.5, 10.5, and 11.6 µm

The purpose of this thesis is to present new observations of thermal-infrared radiation from asteroids. Stellar photometry was performed to provide standards for comparison with the asteroid data. The details of the photometry and the data reduction are discussed in Part 1. A system of standard stars is derived for wavelengths of 8.5, 10.5 and 11.6 µm and a new calibration is adopted. Sources of error are evaluated and comparisons are made with the data of other observers.

The observations and analysis of the thermal-emission observations of asteroids are presented in Part 2. Thermal-emission lightcurve and phase effect data are considered. Special color diagrams are introduced to display the observational data. These diagrams are free of any model-dependent assumptions and show that asteroids differ in their surface properties.

On the basis of photometric models, (4) Vesta is thought to have a bolometric Bond albedo of about 0.1, an emissivity greater than 0.7 and a true radius that is close to the model value of 300^(+50)_(-30)km. Model albedos and model radii are given for asteroids 1, 2, 4, 5, 6, 7, 15, 19, 20, 27, 39, 44, 68, 80, 324 and 674. The asteroid (324) Bamberga is extremely dark with a model (~bolometric Bond) albedo in the 0.01 - 0.02 range, which is thought to be the lowest albedo yet measured for any solar-system body. The crucial question about such low-albedo asteroids is their number and the distribution of their orbits.

DEEP WASHINGTON PHOTOMETRY OF INCONSPICUOUS STAR CLUSTER CANDIDATES IN THE LARGE MAGELLANIC CLOUD

We present deep Washington photometry of 45 poorly populated star cluster candidates in the Large Magellanic Cloud (LMC). We have performed a systematic study to estimate the parameters of the cluster candidates by matching theoretical isochrones to the cleaned and dereddened cluster color-magnitude diagrams. We were able to estimate the basic parameters for 33 clusters, out of which 23 are identified as single clusters and 10 are found to be members of double clusters. The other 12 cluster candidates have been classified as possible clusters/asterisms. About 50% of the true clusters are in the 100-300 Myr age range, whereas some are older or younger. We have discussed the distribution of age, location, and reddening with respect to field, as well as the size of true clusters. The sizes and masses of the studied sample are found to be similar to that of open clusters in the Milky Way. Our study adds to the lower end of cluster mass distribution in the LMC, suggesting that the LMC, apart from hosting rich clusters, also has formed small, less massive open clusters in the 100-300 Myr age range.

Photometry of cometary nuclei: Rotation rates, colours and a comparison with Kuiper Belt Objects

We present time-series data on Jupiter Family Comets (JFCs) 17P/Holmes, 47P/Ashbrook-Jackson and 137P/Shoemaker-Levy 2. In addition we also present results from `snap-shot' observations of comets 43P/Wolf-Harrington, 44P/Reinmuth 2, 103P/Hartley 2 and 104P/Kowal 2 taken during the same run. The comets were at heliocentric distances of between 3 and 7 au at this time. We present measurements of size and activity levels for the snap-shot targets. The time-series data allow us to constrain rotation periods and shapes, and thus bulk densities. We also measure colour indices (V - R) and (R - I) and reliable radii for these comets. We compare all of our findings to date with similar results for other comets and Kuiper Belt Objects (KBOs). We find that the rotational properties of nuclei and KBOs are very similar, that there is evidence for a cut-off in bulk densities at ~0.6 g cm-3 in both populations, and the colours of the two populations show similar correlations. For JFCs, there is no observational evidence for the optical colours being dependent on either position in the orbit or orbital parameters.

Time-resolved optical photometry of the ultracompact binary 4U 0614+091

We present a detailed optical study of the ultracompact X-ray binary 4U 0614+091. We have used 63 hr of time-resolved optical photometry taken with three different telescopes (IAC80, NOT, and SPM) to search for optical modulations. The power spectra of each data set reveals sinusoidal modulations with different periods, which are not always present. The strongest modulation has a period of 51.3 minutes, a semiamplitude of 4.6 mmag, and is present in the IAC80 data. The SPM and NOT data show periods of 42 minutes and 64 minutes, respectively, but with much weaker amplitudes, 2.6 mmag and 1.3 mmag, respectively. These modulations arise from either X-ray irradiation of the inner face of the secondary star and/or a superhump modulation from the accretion disk, or quasiperiodic modulations in the accretion disk. It is unclear whether these periods/quasi-periodic modulations are related to the orbital period; however, the strongest period of 51.3 minutes is close to earlier tentative orbital periods. Further observations taken over a long baseline are encouraged.

ULTRACAM photometry of the eclipsing cataclysmic variables GY Cnc, IR Com and HT Cas

We present high-speed, three-colour photometry of the eclipsing cataclysmic variables GY Cnc, IR Com and HT Cas. We find that the sharp eclipses in GY Cnc and IR Com are due to eclipses of the white dwarf. There is some evidence for a bright-spot on the edge of the accretion disc in GY Cnc, but not in IR Com. Eclipse mapping of HT Cas is presented which shows changes in the structure of the quiescent accretion disc. Observations in 2002 show the accretion disc to be invisible except for the presence of a bright-spot at the disc edge. 2003 observations, however, clearly show a bright inner disc and the bright-spot to be much fainter than in 2002. Although no outburst was associated with either set of quiescent observations, the system was similar to 0.6 mJy brighter in 2003, mainly due to the enhanced emission from the inner disc. We propose that these changes are due to variations in the mass-transfer rate from the secondary star and through the disc. The disc colours indicate that it is optically thin in both its inner and outer regions. We estimate the white dwarf temperature of HT Cas to be 15 000 +/- 1000 K in 2002 and 14 000 +/- 1000 K in 2003.

ULTRACAM photometry of the eclipsing cataclysmic variable OU Vir

We present high-speed, three-colour photometry of the faint eclipsing cataclysmic variable OU Vir. For the first time in OU Vir, separate eclipses of the white dwarf and the bright spot have been observed. We use timings of these eclipses to derive a purely photometric model of the system, obtaining a mass ratio of q=0.175+/-0.025, an inclination of i=79.degrees2+/-0.degrees7 and a disc radius of R-d/a=0.2315+/-0.0150. We separate the white dwarf eclipse from the light curve and, by fitting a blackbody spectrum to its flux in each passband, obtain a white dwarf temperature of T=13900+/-600 K and a distance of D=51+/-17 pc. Assuming that the primary obeys the Nauenberg mass-radius relation for white dwarfs and allowing for temperature effects, we also find a primary mass M-w/M-circle dot=0.89+/-0.20, a primary radius R-w/R-circle dot=0.0097+/-0.0031 and an orbital separation a/R-circle dot=0.74+/-0.05.

Optical and infrared photometry of the Type IInSN1998S: days 11-146

We present contemporaneous optical and infrared (IR) photometric observations of the Type IIn SN 1998S covering the period between 11 and 146 d after discovery. The IR data constitute the first ever IR light curves of a Type IIn supernova. We use blackbody and spline fits to the photometry to examine the luminosity evolution. During the first 2-3 months, the luminosity is dominated by the release of shock-deposited energy in the ejecta. After similar to 100 d the luminosity is powered mostly by the deposition of radioactive decay energy from 0.15 +/-0.05 M-. of Ni-56 which was produced in the explosion. We also report the discovery of an astonishingly high IR excess, K-L'=2.5, that was present at day 130. We interpret this as being due to thermal emission from dust grains in the vicinity of the supernova. We argue that to produce such a high IR luminosity so soon after the explosion, the dust must be pre-existing and so is located in the circumstellar medium of the progenitor. The dust could be heated either by the UV/optical flash (IR echo) or by the X-rays from the interaction of the ejecta with the circumstellar material.

Predicting spectral features in galaxy spectra from broad-band photometry

We explore the prospects of predicting emission-line features present in galaxy spectra given broad-band photometry alone. There is a general consent that colours, and spectral features, most notably the 4000 angstrom break, can predict many properties of galaxies, including star formation rates and hence they could infer some of the line properties. We argue that these techniques have great prospects in helping us understand line emission in extragalactic objects and might speed up future galaxy redshift surveys if they are to target emission-line objects only. We use two independent methods, Artificial Neural Networks (based on the ANNz code) and Locally Weighted Regression (LWR), to retrieve correlations present in the colour N-dimensional space and to predict the equivalent widths present in the corresponding spectra. We also investigate how well it is possible to separate galaxies with and without lines from broad-band photometry only. We find, unsurprisingly, that recombination lines can be well predicted by galaxy colours. However, among collisional lines some can and some cannot be predicted well from galaxy colours alone, without any further redshift information. We also use our techniques to estimate how much information contained in spectral diagnostic diagrams can be recovered from broad-band photometry alone. We find that it is possible to classify active galactic nuclei and star formation objects relatively well using colours only. We suggest that this technique could be used to considerably improve redshift surveys such as the upcoming Fibre Multi Object Spectrograph (FMOS) survey and the planned Wide Field Multi Object Spectrograph (WFMOS) survey.

Global properties of `ordinary` early-type galaxies: photometry and spectroscopy of stars and globular clusters in NGC 4494

We present a comprehensive analysis of the spatial, kinematic and chemical properties of stars and globular clusters (GCs) in the `ordinary` elliptical galaxy NGC 4494 using data from the Keck and Subaru telescopes. We derive galaxy surface brightness and colour profiles out to large galactocentric radii. We compare the latter to metallicities derived using the near-infrared Calcium Triplet. We obtain stellar kinematics out to similar to 3.5 effective radii. The latter appear flattened or elongated beyond similar to 1.8 effective radii in contrast to the relatively round photometric isophotes. In fact, NGC 4494 may be a flattened galaxy, possibly even an S0, seen at an inclination of similar to 45 degrees. We publish a catalogue of 431 GC candidates brighter than i(0) = 24 based on the photometry, of which 109 are confirmed spectroscopically and 54 have measured spectroscopic metallicities. We also report the discovery of three spectroscopically confirmed ultra-compact dwarfs around NGC 4494 with measured metallicities of -0.4 less than or similar to [Fe/H] less than or similar to -0.3. Based on their properties, we conclude that they are simply bright GCs. The metal-poor GCs are found to be rotating with similar amplitude as the galaxy stars, while the metal-rich GCs show marginal rotation. We supplement our analysis with available literature data and results. Using model predictions of galaxy formation, and a suite of merger simulations, we find that many of the observational properties of NGC 4494 may be explained by formation in a relatively recent gas-rich major merger. Complete studies of individual galaxies incorporating a range of observational avenues and methods such as the one presented here will be an invaluable tool for constraining the fine details of galaxy formation models, especially at large galactocentric radii.

TRIPPy: Trailed Image Photometry in Python

Photometry of moving sources typically suffers from a reduced signal-to-noise ratio (S/N) or flux measurements biased to incorrect low values through the use of circular apertures. To address this issue, we present the software package, TRIPPy: TRailed Image Photometry in Python. TRIPPy introduces the pill aperture, which is the natural extension of the circular aperture appropriate for linearly trailed sources. The pill shape is a rectangle with two semicircular end-caps and is described by three parameters, the trail length and angle, and the radius. The TRIPPy software package also includes a new technique to generate accurate model point-spread functions (PSFs) and trailed PSFs (TSFs) from stationary background sources in sidereally tracked images. The TSF is merely the convolution of the model PSF, which consists of a moffat profile, and super-sampled lookup table. From the TSF, accurate pill aperture corrections can be estimated as a function of pill radius with an accuracy of 10 mmag for highly trailed sources. Analogous to the use of small circular apertures and associated aperture corrections, small radius pill apertures can be used to preserve S/Ns of low flux sources, with appropriate aperture correction applied to provide an accurate, unbiased flux measurement at all S/Ns.

Characterising transiting exoplanets by way of differential photometry

This paper describes a simple activity for plotting and characterising the light curve from an exoplanet transit event by way of differential photometry analysis. Using free digital imaging software, participants analyse a series of telescope images with the goal of calculating various exoplanet parameters, including its size, orbital radius and habitability. The activity has been designed for a high-school or undergraduate university level and introduces fundamental concepts in astrophysics and an understanding of the basis for exoplanetary science, the transit method and digital photometry.

Approximate Bayesian Computation for astronomical model analysis : a case study in galaxy demographics and morphological transformation at high redshift

Approximate Bayesian Computation’ (ABC) represents a powerful methodology for the analysis of complex stochastic systems for which the likelihood of the observed data under an arbitrary set of input parameters may be entirely intractable – the latter condition rendering useless the standard machinery of tractable likelihood-based, Bayesian statistical inference [e.g. conventional Markov chain Monte Carlo (MCMC) simulation]. In this paper, we demonstrate the potential of ABC for astronomical model analysis by application to a case study in the morphological transformation of high-redshift galaxies. To this end, we develop, first, a stochastic model for the competing processes of merging and secular evolution in the early Universe, and secondly, through an ABC-based comparison against the observed demographics of massive (Mgal > 1011 M⊙) galaxies (at 1.5 < z < 3) in the Cosmic Assembly Near-IR Deep Extragalatic Legacy Survey (CANDELS)/Extended Groth Strip (EGS) data set we derive posterior probability densities for the key parameters of this model. The ‘Sequential Monte Carlo’ implementation of ABC exhibited herein, featuring both a self-generating target sequence and self-refining MCMC kernel, is amongst the most efficient of contemporary approaches to this important statistical algorithm. We highlight as well through our chosen case study the value of careful summary statistic selection, and demonstrate two modern strategies for assessment and optimization in this regard. Ultimately, our ABC analysis of the high-redshift morphological mix returns tight constraints on the evolving merger rate in the early Universe and favours major merging (with disc survival or rapid reformation) over secular evolution as the mechanism most responsible for building up the first generation of bulges in early-type discs.

Light scattering in random media with wavelength-scale structures : astronomical and industrial applications

Light scattering, or scattering and absorption of electromagnetic waves, is an important tool in all remote-sensing observations. In astronomy, the light scattered or absorbed by a distant object can be the only source of information. In Solar-system studies, the light-scattering methods are employed when interpreting observations of atmosphereless bodies such as asteroids, atmospheres of planets, and cometary or interplanetary dust. Our Earth is constantly monitored from artificial satellites at different wavelengths. With remote sensing of Earth the light-scattering methods are not the only source of information: there is always the possibility to make in situ measurements. The satellite-based remote sensing is, however, superior in the sense of speed and coverage if only the scattered signal can be reliably interpreted. The optical properties of many industrial products play a key role in their quality. Especially for products such as paint and paper, the ability to obscure the background and to reflect light is of utmost importance. High-grade papers are evaluated based on their brightness, opacity, color, and gloss. In product development, there is a need for computer-based simulation methods that could predict the optical properties and, therefore, could be used in optimizing the quality while reducing the material costs. With paper, for instance, pilot experiments with an actual paper machine can be very time- and resource-consuming. The light-scattering methods presented in this thesis solve rigorously the interaction of light and material with wavelength-scale structures. These methods are computationally demanding, thus the speed and accuracy of the methods play a key role. Different implementations of the discrete-dipole approximation are compared in the thesis and the results provide practical guidelines in choosing a suitable code. In addition, a novel method is presented for the numerical computations of orientation-averaged light-scattering properties of a particle, and the method is compared against existing techniques. Simulation of light scattering for various targets and the possible problems arising from the finite size of the model target are discussed in the thesis. Scattering by single particles and small clusters is considered, as well as scattering in particulate media, and scattering in continuous media with porosity or surface roughness. Various techniques for modeling the scattering media are presented and the results are applied to optimizing the structure of paper. However, the same methods can be applied in light-scattering studies of Solar-system regoliths or cometary dust, or in any remote-sensing problem involving light scattering in random media with wavelength-scale structures.