
This whitepaper discusses the diversity of exoplanets that could be detected
by future observations, so that comparative exoplanetology can be performed in
the upcoming era of large spacebased flagship missions. The primary focus will
be on characterizing Earthlike worlds around Sunlike stars. However, we will
also be able to characterize companion planets in the system simultaneously.
This will not only provide a contextual picture with regards to our Solar
system, but also presents a unique opportunity to observe size dependent
planetary atmospheres at different orbital distances. We propose a preliminary
scheme based on chemical behavior of gases and condensates in a planet's
atmosphere that classifies them with respect to planetary radius and incident
stellar flux.

The energy levels of the first few lowlying states of carbon in intense
magnetic fields upwards of $\approx 10^7$ T are calculated in this study. We
extend our previously employed pseudospectral approach for calculating the
eigenstates of the carbon atom. We report data for the ground state and a
lowlying state that are in good agreement with findings elsewhere, as well as
new data for ten other states of the carbon atom that have not been
investigated until now. It is seen that these hitherto uninvestigated states
also become strongly bound with increasing magnetic field strengths. The data
presented in this study are relevant for astrophysical applications, such as
magnetized white dwarf and neutron star spectral analysis as well as opacity
calculations and absorption features, including in the context of material
accreting onto the surfaces of these compact objects.

A numerical method for solving the equations modeling acoustic scattering is
presented. The method is capable of handling several dozen scatterers, each of
which is several wavelengths long, on a personal work station. Even for
geometries involving cavities, solutions accurate to seven digits or better
were obtained. The method relies on a Boundary Integral Equation formulation of
the scattering problem, discretized using a highorder accurate Nystr\"om
method. A hybrid iterative/direct solver is used in which a local scattering
matrix for each body is computed, and then GMRES, accelerated by the Fast
Multipole Method, is used to handle reflections between the scatterers. The
main limitation of the method described is that it currently applies only to
scattering bodies that are rotationally symmetric.

An algorithm for the direct inversion of the linear systems arising from
Nystrom discretization of integral equations on onedimensional domains is
described. The method typically has O(N) complexity when applied to boundary
integral equations (BIEs) in the plane with nonoscillatory kernels such as
those associated with the Laplace and Stokes' equations. The scaling
coefficient suppressed by the "bigO" notation depends logarithmically on the
requested accuracy. The method can also be applied to BIEs with oscillatory
kernels such as those associated with the Helmholtz and Maxwell equations; it
is efficient at long and intermediate wavelengths, but will eventually become
prohibitively slow as the wavelength decreases. To achieve linear complexity,
rank deficiencies in the offdiagonal blocks of the coefficient matrix are
exploited. The technique is conceptually related to the H and H^2 matrix
arithmetic of Hackbusch and coworkers, and is closely related to previous work
on Hierarchically SemiSeparable matrices.

A survey that can cover the sky in optical bands over wide fields to faint
magnitudes with a fast cadence will enable many of the exciting science
opportunities of the next decade. The Large Synoptic Survey Telescope (LSST)
will have an effective aperture of 6.7 meters and an imaging camera with field
of view of 9.6 deg^2, and will be devoted to a tenyear imaging survey over
20,000 deg^2 south of +15 deg. Each pointing will be imaged 2000 times with
fifteen second exposures in six broad bands from 0.35 to 1.1 microns, to a
total pointsource depth of r~27.5. The LSST Science Book describes the basic
parameters of the LSST hardware, software, and observing plans. The book
discusses educational and outreach opportunities, then goes on to describe a
broad range of science that LSST will revolutionize: mapping the inner and
outer Solar System, stellar populations in the Milky Way and nearby galaxies,
the structure of the Milky Way disk and halo and other objects in the Local
Volume, transient and variable objects both at low and high redshift, and the
properties of normal and active galaxies at low and high redshift. It then
turns to farfield cosmological topics, exploring properties of supernovae to
z~1, strong and weak lensing, the largescale distribution of galaxies and
baryon oscillations, and how these different probes may be combined to
constrain cosmological models and the physics of dark energy.

The "collapsar" engine for gammaray bursts invokes as its energy source the
failure of a normal supernova and the formation of a black hole. Here we
present the results of the first threedimensional simulation of the collapse
of a massive star down to a black hole, including the subsequent accretion and
explosion. The explosion differs significantly from the axisymmetric scenario
obtained in twodimensional simulations; this has important consequences for
the nucleosynthetic yields. We compare the nucleosynthetic yields to those of
hypernovae. Calculating yields from threedimensional explosions requires new
strategies in postprocess nucleosynthesis; we discuss NuGrid's plan for
threedimensional yields.

We examine observed heavy element abundances in the Cassiopeia A supernova
remnant as a constraint on the nature of the Cas A supernova. We compare bulk
abundances from 1D and 3D explosion models and spatial distribution of elements
in 3D models with those derived from Xray observations. We also examine the
cospatial production of 26Al with other species. We find that the most reliable
indicator of the presence of 26Al in unmixed ejecta is a very low S/Si ratio
(~0.05). Production of N in O/S/Sirich regions is also indicative. The
biologically important element P is produced at its highest abundance in the
same regions. Proxies should be detectable in supernova ejecta with high
spatial resolution multiwavelength observations.

We obtained R and Iband CCD photometry of the soft Xray transient/neutron
star binary Aql X1 in 1998 June while it was at quiescence. We find that its
light curve is dominated by ellipsoidal variations, although the ellipsoidal
variations are severely distorted and have unequal maxima. After we correct for
the contaminating flux from a field star located only 0.46" away, the
peaktopeak amplitude of the modulation is ~0.25 mag in the R band, which
requires the orbital inclination to be greater than 36 degrees. The orbital
period we measure is consistent with the 18.95 h period measured by Chevalier &
Ilovaisky (1998). During its outbursts the light curve of Aql X1 becomes
single humped. The outburst light curve observed by Garcia et al. (1999) agrees
in phase with our quiescent light curve. We show that the single humped
variation is caused by a ``reflection effect,'' that is, by heating of the side
of the secondary star facing towards the neutron star.