The habitability of the surface of any planet is determined by a complex evolution of its interior, surface, and atmosphere. The electromagnetic and particle radiation of stars drive thermal, chemical and physical alteration of planetary atmospheres, including escape. Many known extrasolar planets experience vastly different stellar environments than those in our Solar system: it is crucial to understand the broad range of processes that lead to atmospheric escape and evolution under a wide range of conditions if we are to assess the habitability of worlds around other stars. One problem encountered between the planetary and the astrophysics communities is a lack of common language for describing escape processes. Each community has customary approximations that may be questioned by the other, such as the hypothesis of H-dominated thermosphere for astrophysicists, or the Sun-like nature of the stars for planetary scientists. Since exoplanets are becoming one of the main targets for the detection of life, a common set of definitions and hypotheses are required. We review the different escape mechanisms proposed for the evolution of planetary and exoplanetary atmospheres. We propose a common definition for the different escape mechanisms, and we show the important parameters to take into account when evaluating the escape at a planet in time. We show that the paradigm of the magnetic field as an atmospheric shield should be changed and that recent work on the history of Xenon in Earth's atmosphere gives an elegant explanation to its enrichment in heavier isotopes: the so-called Xenon paradox.

We describe recent developments to the Castro astrophysics simulation code, focusing on new features that enable our simulations of X-ray bursts. Two highlights of Castro's ongoing development are the new integration technique to couple hydrodynamics and reactions to high order and GPU offloading. We discuss how these features will help offset some of the computational expense in X-ray burst models.

Simulations in stellar astrophysics involve the coupling of hydrodynamics and nuclear reactions under a wide variety of conditions, from simmering convective flows to explosive nucleosynthesis. Numerical techniques such as operator splitting (most notably Strang splitting) are usually employed to couple the physical processes, but this can affect the accuracy of the simulation, particularly when the burning is vigorous. Furthermore, Strang splitting does not have a straightforward extension to higher-order integration in time. We present a new temporal integration strategy based on spectral deferred corrections and describe the second- and fourth-order implementations in the open-source, finite-volume, compressible hydrodynamics code Castro. One notable advantage to these schemes is that they combine standard low-order discretizations for individual physical processes in a way that achieves an arbitrarily high order of accuracy. We demonstrate the improved accuracy of the new methods on several test problems of increasing complexity.

Unexpected dynamic phenomena have surprised solar system observers in the past and have led to important discoveries about solar system workings. Observations at the initial stages of these events provide crucial information on the physical processes at work. We advocate for long-term/permanent programs on ground-based and space-based telescopes of all sizes - including Extremely Large Telescopes (ELTs) - to conduct observations of high-priority dynamic phenomena, based on a predefined set of triggering conditions. These programs will ensure that the best initial dataset of the triggering event are taken; separate additional observing programs will be required to study the temporal evolution of these phenomena. While not a comprehensive list, the following are notional examples of phenomena that are rare, that cannot be anticipated, and that provide high-impact advances to our understandings of planetary processes. Examples include: new cryovolcanic eruptions or plumes on ocean worlds; impacts on Jupiter, Saturn, Uranus, or Neptune; extreme eruptions on Io; convective superstorms on Saturn, Uranus, or Neptune; collisions within the asteroid belt or other small-body populations; discovery of an interstellar object passing through our solar system (e.g. 'Oumuamua); and responses of planetary atmospheres to major solar flares or coronal mass ejections.

We discuss the challenges of modeling X-ray bursts in multi-dimensions, review the different calculations done to date, and discuss our new set of ongoing simulations. We also describe algorithmic improvements that may help in the future to offset some of the expense of these simulations, and describe what may be possible with exascale computing.

We present a method to quantify the upper-limit of the energy transmitted from the intense stellar wind to the upper atmospheres of three of the Trappist-1 planets (e, f, and g). We use a formalism that treats the system as two electromagnetic regions, where the efficiency of the energy transmission between one region (the stellar wind at the planetary orbits) to the other (the planetary ionospheres) depends on the relation between the conductances and impedances of the two regions. Since the energy flux of the stellar wind is very high at these planetary orbits, we find that for the case of high transmission efficiency (when the conductances and impedances are close in magnitude), the energy dissipation in the upper planetary atmospheres is also very large. On average, the Ohmic energy can reach $0.5-1~W/m^2$, about 1\% of the stellar irradiance and 5-15 times the EUV irradiance. Here, using constant values for the ionospheric conductance, we demonstrate that the stellar wind energy could potentially drive large atmospheric heating in terrestrial planets, as well as in hot jupiters. More detailed calculations are needed to assess the ionospheric conductance and to determine more accurately the amount of heating the stellar wind can drive in close-orbit planets.

We describe the AMReX suite of astrophysics codes and their application to modeling problems in stellar astrophysics. Maestro is tuned to efficiently model subsonic convective flows while Castro models the highly compressible flows associated with stellar explosions. Both are built on the block-structured adaptive mesh refinement library AMReX. Together, these codes enable a thorough investigation of stellar phenomena, including Type Ia supernovae and X-ray bursts. We describe these science applications and the approach we are taking to make these codes performant on current and future many-core and GPU-based architectures.

We derive water vapor column abundances and aerosol properties from Mars Science Laboratory (MSL) ChemCam passive mode observations of scattered sky light. Each ChemCam passive sky observation acquires spectra at two different elevation angles. We fit these spectra with a discrete-ordinates multiple scattering radiative transfer model, using the correlated-k approximation for gas absorption bands. The retrieval proceeds by first fitting the continuum of the ratio of the two elevation angles to solve for aerosol properties, and then fitting the continuum-removed ratio to solve for gas abundances. The final step of the retrieval makes use of the observed CO2 absorptions and the known CO2 abundance to correct the retrieved water vapor abundance for the effects of the vertical distribution of scattering aerosols and to derive an aerosol scale height parameter. The ChemCam-retrieved water abundances show, with only a few exceptions, the same seasonal behavior and the same timing of seasonal minima and maxima as the TES, CRISM, and REMS-H data sets that we compare them to. However ChemCam-retrieved water abundances are generally lower than zonal and regional scale from-orbit water vapor data, while at the same time being significantly larger than pre-dawn REMS-H abundances. Pending further analysis of REMS-H volume mixing ratio uncertainties, the differences between ChemCam and REMS-H pre-dawn mixing ratios appear to be much too large to be explained by large scale circulations and thus they tend to support the hypothesis of substantial diurnal interactions of water vapor with the surface. Our preliminary aerosol results, meanwhile, show the expected seasonal pattern in dust particle size but also indicate a surprising inter-annual increase in water-ice cloud opacities.

Over the last decade block-structured adaptive mesh refinement (SAMR) has found increasing use in large, publicly available codes and frameworks. SAMR frameworks have evolved along different paths. Some have stayed focused on specific domain areas, others have pursued a more general functionality, providing the building blocks for a larger variety of applications. In this survey paper we examine a representative set of SAMR packages and SAMR-based codes that have been in existence for half a decade or more, have a reasonably sized and active user base outside of their home institutions, and are publicly available. The set consists of a mix of SAMR packages and application codes that cover a broad range of scientific domains. We look at their high-level frameworks, and their approach to dealing with the advent of radical changes in hardware architecture. The codes included in this survey are BoxLib, Cactus, Chombo, Enzo, FLASH, and Uintah.

The dynamics of helium shell convection driven by nuclear burning establish the conditions for runaway in the sub-Chandrasekhar mass, double detonation model for Type Ia supernovae, as well as for a variety of other explosive phenomena. We explore these convection dynamics for a range of white dwarf core and helium shell masses in three dimensions using the low Mach number hydrodynamics code Maestro. We present calculations of the bulk properties of this evolution, including time-series evolution of global diagnostics, lateral averages of the 3D state, and the global 3D state. We find a variety of outcomes including quasi-equilibrium, localized runaway, and nova-like runaway. Our results suggest the double detonation progenitor model is promising, that 3D, dynamic convection plays a key role, and that these systems warrant further study.

We study the interaction between the atmospheres of Venus-like, non-magnetized exoplanets orbiting an M-dwarf star, and the stellar wind using a multi-species Magnetohydrodynaic (MHD) model. We focus our investigation on the effect of enhanced stellar wind and enhanced EUV flux as the planetary distance from the star decreases. Our simulations reveal different topologies of the planetary space environment for sub- and super-Alfvenic stellar wind conditions, which could lead to dynamic energy deposition in to the atmosphere during the transition along the planetary orbit. We find that the stellar wind penetration for non-magnetized planets is very deep, up to a few hundreds of kilometers. We estimate a lower limit for the atmospheric mass-loss rate and find that it is insignificant over the lifetime of the planet. However, we predict that when accounting for atmospheric ion acceleration, a significant amount of the planetary atmosphere could be eroded over the course of a billion years.

We perform the first detailed three-dimensional simulation of low Mach number convection preceding thermonuclear ignition in a mixed H/He X-ray burst. Our simulations include a moderate-sized, approximate network that captures hydrogen and helium burning up through rp-process breakout. We look in detail at the difference between two- and three-dimensional convective fields, including the details of the turbulent convection.

We present new kinetics data on the radiolytic destruction of amino acids measured in situ with infrared spectroscopy. Samples were irradiated at 15, 100, and 140 K with 0.8-MeV protons, and amino-acid decay was followed at each temperature with and without H$_2$O present. Observed radiation products included CO$_2$ and amines, consistent with amino-acid decarboxylation. The half-lives of glycine, alanine, and phenylalanine were estimated for various extraterrestrial environments. Infrared spectral changes demonstrated the conversion from the non-zwitterion structure NH$_2$-CH$_2$(R)-COOH at 15 K to the zwitterion structure $^+$NH$_3$-CH$_2$(R)-COO$^-$ at 140 K for each amino acid studied.

We study the magnetospheric structure and the ionospheric Joule Heating of planets orbiting M-dwarf stars in the habitable zone using a set of magnetohydrodynamic (MHD) models. The stellar wind solution is used to drive a model for the planetary magnetosphere, which is coupled with a model for the planetary ionosphere. Our simulations reveal that the space environment around close-in habitable planets is extreme, and the stellar wind plasma conditions change from sub- to super-Alfvenic along the planetary orbit. As a result, the magnetospheric structure changes dramatically with a bow shock forming in the super-Alfvenic sectors, while no bow shock forms in the sub-Alfvenic sectors. The planets reside most of the time in the sub-Alfvenic sectors with poor atmospheric protection. A significant amount of Joule Heating is provided at the top of the atmosphere as a result of the planetary interaction with the stellar wind. For the steady-state solution, the heating is about 0.1-3\% of the total incoming stellar irradiation, and it is enhanced by 50\% for the time-dependent case. The significant Joule Heating obtained here should be considered in models for the atmospheres of habitable planets in terms of the thickness of the atmosphere, the top-side temperature and density, the boundary conditions for the atmospheric pressure, and particle radiation and transport.

Type I X-ray Bursts (XRBs) are thermonuclear explosions of accreted material on the surfaces of a neutron stars in low mass X-ray binaries. Prior to the ignition of a subsonic burning front, runaway burning at the base of the accreted layer drives convection that mixes fuel and heavy-element ashes. In this second paper in a series, we explore the behavior of this low Mach number convection in mixed hydrogen/helium layers on the surface of a neutron star using two-dimensional simulations with the Maestro code. Maestro takes advantage of the highly subsonic flow field by filtering dynamically unimportant sound waves while retaining local compressibility effects, such as those due to stratification and energy release from nuclear reactions. In these preliminary calculations, we find that the rp-process approximate network creates a convective region that is split into two layers. While this splitting appears artificial due to the approximations of the network regarding nuclear flow out of the breakout reaction 18Ne(a,p)21Na, these calculations hint at further simplifications and improvements of the burning treatment for use in subsequent calculations in three dimensions for a future paper.

Dust aerosol plays a fundamental role in the behavior and evolution of the Martian atmosphere. The first five Mars years of Mars Exploration Rover data provide an unprecedented record of the dust load at two sites. This record is useful for characterization of the atmosphere at the sites and as ground truth for orbital observations. Atmospheric extinction optical depths have been derived from solar images after calibration and correction for time-varying dust that has accumulated on the camera windows. The record includes local, regional, and globally extensive dust storms. Comparison with contemporaneous thermal infrared data suggests significant variation in the size of the dust aerosols, with a 1 {\mu}m effective radius during northern summer and a 2 {\mu}m effective radius at the onset of a dust lifting event. The solar longitude (LS) 20-136{\deg} period is also characterized by the presence of cirriform clouds at the Opportunity site, especially near LS=50 and 115{\deg}. In addition to water ice clouds, a water ice haze may also be present, and carbon dioxide clouds may be present early in the season. Variations in dust opacity are important to the energy balance of each site, and work with seasonal variations in insolation to control dust devil frequency at the Spirit site.

We investigate turbulence generated by cosmological structure formation by means of large eddy simulations using adaptive mesh refinement. In contrast to the widely used implicit large eddy simulations, which resolve a limited range of length scales and treat the effect of turbulent velocity fluctuations below the grid scale solely by numerical dissipation, we apply a subgrid-scale model for the numerically unresolved fraction of the turbulence energy. For simulations with adaptive mesh refinement, we utilize a new methodology that allows us to adjust the scale-dependent energy variables in such a way that the sum of resolved and unresolved energies is globally conserved. We test our approach in simulations of randomly forced turbulence, a gravitationally bound cloud in a wind, and the Santa Barbara cluster. To treat inhomogeneous turbulence, we introduce an adaptive Kalman filtering technique that separates turbulent velocity fluctuations on resolved length scales from the non-turbulent bulk flow. From the magnitude of the fluctuating component and the subgrid-scale turbulence energy, a total turbulent velocity dispersion of several 100 km/s is obtained for the Santa Barbara cluster, while the low-density gas outside the accretion shocks is nearly devoid of turbulence. The energy flux through the turbulent cascade and the dissipation rate predicted by the subgrid-scale model correspond to dynamical time scales around 5 Gyr, independent of numerical resolution.

We present high-resolution, full-star simulations of the post-ignition phase of Type Ia supernovae using the compressible hydrodynamics code Castro. Initial conditions, including the turbulent velocity field and ignition site, are imported directly from a simulation of the last few hours of presupernova convection using a low Mach number code, Maestro. Adaptive mesh refinement allows the initial burning front to be modeled with an effective resolution of 36,864^3 zones (~136 m/zone). The initial rise and expansion of the deflagration front are tracked until burning reaches the star's edge and the role of the background turbulence on the flame is investigated. The effect of artificially moving the ignition location closer to the star's center is explored. The degree to which turbulence affects the burning front decreases with increasing ignition radius since the buoyancy force is stronger at larger radii. Even central ignition --- in the presence of a background convective flow field --- is rapidly carried off-center as the flame is carried by the flow field. We compare our results to analytic models for burning thermals, and find that they reproduce the general trends of the bubble's size and mass, but underpredict the amount of buoyant acceleration due to simplifying assumptions of the bubble's properties. Overall, we find that the amount of mass that burns prior to flame break out is small, consistent with a "gravitationally confined detonation" occurring at a later epoch, but additional burning will occur following breakout that may modify this conclusion.

A leading model for Type Ia supernovae (SNe Ia) begins with a white dwarf near the Chandrasekhar mass that ignites a degenerate thermonuclear runaway close to its center and explodes. In a series of papers, we shall explore the consequences of ignition at several locations within such dwarfs. Here we assume central ignition, which has been explored before, however, the problem is worth revisiting, if only to validate those previous studies and to further elucidate the relevant physics for future work. A perturbed sphere of hot iron ash with a radius of ~100 km is initialized at the middle of the star. The subsequent explosion is followed in several simulations using a thickened flame model in which the flame speed is either fixed --- within the range expected from turbulent combustion --- or based on the local turbulent intensity. Global results, including the explosion energy and bulk nucleosynthesis (e.g. 56Ni of 0.48--0.56 $\Msun$) turn out to be insensitive to this speed. In all completed runs, the energy released by the nuclear burning is adequate to unbind the star, but not enough to give the energy and brightness of typical SNe Ia. As found previously, the chemical stratification observed in typical events is not reproduced. These models produce a large amount of unburned carbon and oxygen in central low velocity regions, which is inconsistent with spectroscopic observations, and the intermediate mass elements and iron group elements are strongly mixed during the explosion.

We present a new N-body and gas dynamics code, called Nyx, for large-scale cosmological simulations. Nyx follows the temporal evolution of a system of discrete dark matter particles gravitationally coupled to an inviscid ideal fluid in an expanding universe. The gas is advanced in an Eulerian framework with block-structured adaptive mesh refinement (AMR); a particle-mesh (PM) scheme using the same grid hierarchy is used to solve for self-gravity and advance the particles. Computational results demonstrating the validation of Nyx on standard cosmological test problems, and the scaling behavior of Nyx to 50,000 cores, are presented.

We assess the robustness of a low Mach number hydrodynamics algorithm for modeling helium shell convection on the surface of a white dwarf in the context of the sub-Chandrasekhar model for Type Ia supernovae. We use the low Mach number stellar hydrodynamics code, MAESTRO, to perform three-dimensional, spatially-adaptive simulations of convection leading up to the point of the ignition of a burning front. We show that the low Mach number hydrodynamics model provides a robust description of the system.

We present a formulation for multigroup radiation hydrodynamics that is correct to order $O(v/c)$ using the comoving-frame approach and the flux-limited diffusion approximation. We describe a numerical algorithm for solving the system, implemented in the compressible astrophysics code, CASTRO. CASTRO uses an Eulerian grid with block-structured adaptive mesh refinement based on a nested hierarchy of logically-rectangular variable-sized grids with simultaneous refinement in both space and time. In our multigroup radiation solver, the system is split into three parts, one part that couples the radiation and fluid in a hyperbolic subsystem, another part that advects the radiation in frequency space, and a parabolic part that evolves radiation diffusion and source-sink terms. The hyperbolic subsystem and the frequency space advection are solved explicitly with high-order Godunov schemes, whereas the parabolic part is solved implicitly with a first-order backward Euler method. Our multigroup radiation solver works for both neutrino and photon radiation.

We extend our previous three-dimensional, full-star simulations of the final hours of convection preceding ignition in Type Ia supernovae to higher resolution using the adaptive mesh refinement capability of our low Mach number code, MAESTRO. We report the statistics of the ignition of the first flame at an effective 4.34 km resolution, and general flow field properties at an effective 2.17 km resolution. We find that off-center ignition is likely, with radius of 50 km most favored and a likely range of 40 to 75 km. This is consistent with our previous coarser (8.68 km resolution) simulations, implying that we have achieved sufficient resolution in our determination of likely ignition radii. The dynamics of the last few hot spots preceding ignition suggest that a multiple ignition scenario is not likely. With improved resolution, we can more clearly see the general flow pattern in the convective region, characterized by a strong outward plume with a lower speed recirculation. We show that the convective core is turbulent with a Kolmogorov spectrum and has a lower turbulent intensity and larger integral length scale than previously thought (on the order of 16 km s$^{-1}$ and 200 km, respectively), and we discuss the potential consequences for the first flames.

We present our end-to-end capability for computing the convective phase through the explosion phase of Type Ia supernovae. We compute the convective phase up to the time of ignition using our low Mach number code, MAESTRO, and the subsequent explosion phase using our compressible code, CASTRO. Both codes share the same BoxLib software framework and use finite-volume, block-structured adaptive mesh refinement (AMR) to enable high-resolution, three-dimensional full-star simulations that scale to 100,000+ cores. We present preliminary results from the first-ever simulations of convection preceding ignition using MAESTRO with AMR. We also demonstrate our ability to initialize a compressible simulation of the explosion phase in CASTRO using data obtained directly from MAESTRO just before ignition. Some care must be taken during this initialization procedure when interpreting the size and distribution of hot spots.

We develop a one-dimensional theoretical model for thermals burning in Type Ia supernovae based on the entrainment assumption of Morton, Taylor and Turner. Extensions of the standard model are required to account for the burning and for the expansion of the thermal due to changes in the background stratification found in the full star. The model is compared with high-resolution three-dimensional numerical simulations, both in a uniform environment, and in a full-star setting. The simulations in a uniform environment present compelling agreement with the predicted power-laws and provide model constants for the full-star model, which then provides excellent agreement with the full-star simulation. The importance of the different components in the model are compared, and are all shown to be relevant. An examination of the effect of initial conditions was then conducted using the one-dimensional model, which would have been infeasible in three dimensions. More mass was burned when the ignition kernel was larger and closer to the center of the star. The turbulent flame speed was found to be important during the early-time evolution of the thermal, but played a diminished role at later times when the evolution is dominated by the large-scale hydrodynamics responsible for entrainment. However, a higher flame speed effectively gave a larger initial ignition kernel and so resulted in more mass burned. This suggests that future studies should focus on the early-time behavior of these thermals (in particular, the transition to turbulence), and that the choice of turbulent flame speed does not play a significant role in the dynamics once the thermal has become established.

We describe the development of a flux-limited gray radiation solver for the compressible astrophysics code, CASTRO. CASTRO uses an Eulerian grid with block-structured adaptive mesh refinement based on a nested hierarchy of logically-rectangular variable-sized grids with simultaneous refinement in both space and time. The gray radiation solver is based on a mixed-frame formulation of radiation hydrodynamics. In our approach, the system is split into two parts, one part that couples the radiation and fluid in a hyperbolic subsystem, and another parabolic part that evolves radiation diffusion and source-sink terms. The hyperbolic subsystem is solved explicitly with a high-order Godunov scheme, whereas the parabolic part is solved implicitly with a first-order backward Euler method.

In previous studies, we examined turbulence-flame interactions in carbon-burning thermonuclear flames in Type Ia supernovae. In this study, we consider turbulence-flame interactions in the trailing oxygen flames. The two aims of the paper are to examine the response of the inductive oxygen flame to intense levels of turbulence, and to explore the possibility of transition to detonation in the oxygen flame. Scaling arguments analogous to the carbon flames are presented and then compared against three-dimensional simulations for a range of Damk\"ohler numbers ($\Da_{16}$) at a fixed Karlovitz number. The simulations suggest that turbulence does not significantly affect the oxygen flame when $\Da_{16}<1$, and the flame burns inductively some distance behind the carbon flame. However, for $\Da_{16}>1$, turbulence enhances heat transfer and drives the propagation of a flame that is {\em narrower} than the corresponding inductive flame would be. Furthermore, burning under these conditions appears to occur as part of a combined carbon-oxygen turbulent flame with complex compound structure. The simulations do not appear to support the possibility of a transition to detonation in the oxygen flame, but do not preclude it either.

In the distributed burning regime, turbulence disrupts the internal structure of the flame, and so the idea of laminar burning propagated by conduction is no longer valid. The nature of the burning depends on the turbulent Damkohler number (Da), which steadily declines from much greater than one to less that one as the density decreases to a few 10^6 g/cc. Scaling arguments predict that the turbulent flame speed s, normalized by the turbulent intensity u, follows s/u=Da^1/2 for Da<1. The flame in this regime is a single turbulently-broadened structure that moves at a steady speed, and has a width larger than the integral scale of the turbulence. The scaling is predicted to break down at Da=1, and the flame burns as a turbulently-broadened effective unity Lewis number flame. We refer to this kind of flame as a lambda-flame. The burning becomes a collection of lambda-flames spread over a region approximately the size of the integral scale. While the total burning rate continues to have a well-defined average, s_{T} ~ u, the burning is unsteady. We present a theoretical framework, supported by both 1D and 3D numerical simulations, for the burning in these two regimes. Our results indicate that the average value of s can actually be roughly twice u for Da>1, and that localized excursions to as much as five times u can occur. The lambda-flame speed and width can be predicted based on the turbulence in the star and the turbulent nuclear burning time scale of the fuel. We propose a practical method for measuring these based on the scaling relations and small-scale computationally-inexpensive simulations. This suggests that a simple turbulent flame model can be easily constructed suitable for large-scale distributed supernovae flames.

We present multidimensional simulations of the early convective phase preceding ignition in a Type I X-ray burst using the low Mach number hydrodynamics code, MAESTRO. A low Mach number approach is necessary in order to perform long-time integration required to study such phenomena. Using MAESTRO, we are able to capture the expansion of the atmosphere due to large-scale heating while capturing local compressibility effects such as those due to reactions and thermal diffusion. We also discuss the preparation of one-dimensional initial models and the subsequent mapping into our multidimensional framework. Our method of initial model generation differs from that used in previous multidimensional studies, which evolved a system through multiple bursts in one dimension before mapping onto a multidimensional grid. In our multidimensional simulations, we find that the resolution necessary to properly resolve the burning layer is an order of magnitude greater than that used in the earlier studies mentioned above. We characterize the convective patterns that form and discuss their resulting influence on the state of the convective region, which is important in modeling the outburst itself.

Performing high-resolution, high-fidelity, three-dimensional simulations of Type Ia supernovae (SNe Ia) requires not only algorithms that accurately represent the correct physics, but also codes that effectively harness the resources of the most powerful supercomputers. We are developing a suite of codes that provide the capability to perform end-to-end simulations of SNe Ia, from the early convective phase leading up to ignition to the explosion phase in which deflagration/detonation waves explode the star to the computation of the light curves resulting from the explosion. In this paper we discuss these codes with an emphasis on the techniques needed to scale them to petascale architectures. We also demonstrate our ability to map data from a low Mach number formulation to a compressible solver.

We explore the dependence on spatial dimension of the viability of the neutrino heating mechanism of core-collapse supernova explosions. We find that the tendency to explode is a monotonically increasing function of dimension, with 3D requiring $\sim$40$-$50\% lower driving neutrino luminosity than 1D and $\sim$15$-$25\% lower driving neutrino luminosity than 2D. Moreover, we find that the delay to explosion for a given neutrino luminosity is always shorter in 3D than 2D, sometimes by many hundreds of milliseconds. The magnitude of this dimensional effect is much larger than the purported magnitude of a variety of other effects, such as nuclear burning, inelastic scattering, or general relativity, which are sometimes invoked to bridge the gap between the current ambiguous and uncertain theoretical situation and the fact of robust supernova explosions. Since real supernovae occur in three dimensions, our finding may be an important step towards unraveling one of the most problematic puzzles in stellar astrophysics. In addition, even though in 3D we do see pre-explosion instabilities and blast asymmetries, unlike the situation in 2D, we do not see an obvious axially-symmetric dipolar shock oscillation. Rather, the free energy available to power instabilites seems to be shared by more and more degrees of freedom as the dimension increases. Hence, the strong dipolar axisymmetry seen in 2D and previously identified as a fundamental characteristic of the shock hydrodynamics may not survive in 3D as a prominent feature.

We present a new code, CASTRO, that solves the multicomponent compressible hydrodynamic equations for astrophysical flows including self-gravity, nuclear reactions and radiation. CASTRO uses an Eulerian grid and incorporates adaptive mesh refinement (AMR). Our approach to AMR uses a nested hierarchy of logically-rectangular grids with simultaneous refinement in both space and time. The radiation component of CASTRO will be described in detail in the next paper, Part II, of this series.

Many astrophysical phenomena are highly subsonic, requiring specialized numerical methods suitable for long-time integration. In a series of earlier papers we described the development of MAESTRO, a low Mach number stellar hydrodynamics code that can be used to simulate long-time, low-speed flows that would be prohibitively expensive to model using traditional compressible codes. MAESTRO is based on an equation set derived using low Mach number asymptotics; this equation set does not explicitly track acoustic waves and thus allows a significant increase in the time step. MAESTRO is suitable for two- and three-dimensional local atmospheric flows as well as three-dimensional full-star flows. Here, we continue the development of MAESTRO by incorporating adaptive mesh refinement (AMR). The primary difference between MAESTRO and other structured grid AMR approaches for incompressible and low Mach number flows is the presence of the time-dependent base state, whose evolution is coupled to the evolution of the full solution. We also describe how to incorporate the expansion of the base state for full-star flows, which involves a novel mapping technique between the one-dimensional base state and the Cartesian grid, as well as a number of overall improvements to the algorithm. We examine the efficiency and accuracy of our adaptive code, and demonstrate that it is suitable for further study of our initial scientific application, the convective phase of Type Ia supernovae.

The inclusion of rotationally-induced mixing in stellar evolution can alter the structure and composition of presupernova stars. We survey the effects of progenitor rotation on nucleosynthetic yields in Population III and II supernovae using the new adaptive mesh refinement (AMR) code CASTRO. We examine spherical explosions in 15, 25 and 40 solar mass stars at Z = 0 and 10^-4 solar metallicity with three explosion energies and two rotation rates. Rotation in the Z = 0 models resulted in primary nitrogen production and a stronger hydrogen burning shell which led all models to die as red supergiants. On the other hand, the Z=10^-4 solar metallicity models that included rotation ended their lives as compact blue stars. Because of their extended structure, the hydrodynamics favors more mixing and less fallback in the metal free stars than the Z = 10^-4 models. As expected, higher energy explosions produce more enrichment and less fallback than do lower energy explosions, and less massive stars produce more enrichment and leave behind smaller remnants than do more massive stars. We compare our nucleosynthetic yields to the chemical abundances in the three most iron-poor stars yet found and reproduce the abundance pattern of one, HE 0557-4840, with a zero metallicity 15 solar mass, 2.4 x 10^51 erg supernova. A Salpeter IMF averaged integration of our yields for Z=0 models with explosion energies of 2.4x10^51 ergs or less is in good agreement with the abundances observed in larger samples of extremely metal-poor stars, provided 15 solar mass stars are included. Since the abundance patterns of extremely metal-poor stars likely arise from a representative sample of progenitors, our yields suggest that low-mass supernovae contributed the bulk of the metals to the early universe.

We present the first three-dimensional, full-star simulations of convection in a white dwarf preceding a Type Ia supernova, specifically the last few hours before ignition. For these long-time calculations we use our low Mach number hydrodynamics code, MAESTRO, which we have further developed to treat spherical stars centered in a three-dimensional Cartesian geometry. The main change required is a procedure to map the one-dimensional radial base state to and from the Cartesian grid. Our models recover the dipole structure of the flow seen in previous calculations, but our long-time integration shows that the orientation of the dipole changes with time. Furthermore, we show the development of gravity waves in the outer, stable portion of the star. Finally, we evolve several calculations to the point of ignition and discuss the range of ignition radii.

The large range of time and length scales involved in type Ia supernovae (SN Ia) requires the use of flame models. As a prelude to exploring various options for flame models, we consider, in this paper, high-resolution three-dimensional simulations of the small-scale dynamics of nuclear flames in the supernova environment in which the details of the flame structure are fully resolved. The range of densities examined, 1 to $8 \times 10^7$ g cm$^{-3}$, spans the transition from the laminar flamelet regime to the distributed burning regime where small scale turbulence disrupts the flame. The use of a low Mach number algorithm facilitates the accurate resolution of the thermal structure of the flame and the inviscid turbulent kinetic energy cascade, while implicitly incorporating kinetic energy dissipation at the grid-scale cutoff. For an assumed background of isotropic Kolmogorov turbulence with an energy characteristic of SN Ia, we find a transition density between 1 and $3 \times 10^7$ g cm$^{-3}$ where the nature of the burning changes qualitatively. By $1 \times 10^7$ g cm$^{-3}$, energy diffusion by conduction and radiation is exceeded, on the flame scale, by turbulent advection. As a result, the effective Lewis Number approaches unity. That is, the flame resembles a laminar flame, but is turbulently broadened with an effective diffusion coefficient, $D_T \sim u' l$, where $u'$ is the turbulent intensity and $l$ is the integral scale. For the larger integral scales characteristic of a real supernova, the flame structure is predicted to become complex and unsteady. Implications for a possible transition to detonation are discussed.

We report observations of the nova RS Ophiuchi (RS Oph) using the Keck Interferometer Nuller (KIN), approximately 3.8 days following the most recent outburst that occurred on 2006 February 12. These observations represent the first scientific results from the KIN, which operates in N-band from 8 to 12.5 microns in a nulling mode. By fitting the unique KIN data, we have obtained an angular size of the mid-infrared continuum of 6.2, 4.0, or 5.4 mas for a disk profile, gaussian profile (FWHM), and shell profile respectively. The data show evidence of enhanced neutral atomic hydrogen emission and atomic metals including silicon located in the inner spatial regime near the white dwarf (WD) relative to the outer regime. There are also nebular emission lines and evidence of hot silicate dust in the outer spatial region, centered at ! 17 AU from the WD, that are not found in the inner regime. Our evidence suggests that these features have been excited by the nova flash in the outer spatial regime before the blast wave reached these regions. These identifications support a model in which the dust appears to be present between outbursts and is not created during the outburst event. We further discuss the present results in terms of a unifying model of the system that includes an increase in density in the plane of the orbit of the two stars created by a spiral shock wave caused by the motion of the stars through the cool wind of the red giant star. These data show the power and potential of the nulling technique which has been developed for the detection of Earth-like planets around nearby stars for the Terrestrial Planet Finder Mission and Darwin missions.

Type Ia supernovae (SNe Ia) are the largest thermonuclear explosions in the Universe. Their light output can be seen across great distances and has led to the discovery that the expansion rate of the Universe is accelerating. Despite the significance of SNe Ia, there are still a large number of uncertainties in current theoretical models. Computational modeling offers the promise to help answer the outstanding questions. However, even with today's supercomputers, such calculations are extremely challenging because of the wide range of length and time scales. In this paper, we discuss several new algorithms for simulations of SNe Ia and demonstrate some of their successes.

The convective period leading up to a Type Ia supernova (SN Ia) explosion is characterized by very low Mach number flows, requiring hydrodynamical methods well-suited to long-time integration. We continue the development of the low Mach number equation set for stellar scale flows by incorporating the effects of heat release due to external sources. Low Mach number hydrodynamics equations with a time-dependent background state are derived, and a numerical method based on the approximate projection formalism is presented. We demonstrate through validation with a fully compressible hydrodynamics code that this low Mach number model accurately captures the expansion of the stellar atmosphere as well as the local dynamics due to external heat sources. This algorithm provides the basis for an efficient simulation tool for studying the ignition of SNe Ia.

We introduce a low Mach number equation set for the large-scale numerical simulation of carbon-oxygen white dwarfs experiencing a thermonuclear deflagration. Since most of the interesting physics in a Type Ia supernova transpires at Mach numbers from 0.01 to 0.1, such an approach enables both a considerable increase in accuracy and savings in computer time compared with frequently used compressible codes. Our equation set is derived from the fully compressible equations using low Mach number asymptotics, but without any restriction on the size of perturbations in density or temperature. Comparisons with simulations that use the fully compressible equations validate the low Mach number model in regimes where both are applicable. Comparisons to simulations based on the more traditional anelastic approximation also demonstrate the agreement of these models in the regime for which the anelastic approximation is valid. For low Mach number flows with potentially finite amplitude variations in density and temperature, the low Mach number model overcomes the limitations of each of the more traditional models and can serve as the basis for an accurate and efficient simulation tool.

We extend a low Mach number hydrodynamics method developed for terrestrial combustion, to the study of thermonuclear flames in Type Ia supernovae. We discuss the differences between 2-D and 3-D Rayleigh-Taylor unstable flame simulations, and give detailed diagnostics on the turbulence, showing that the kinetic energy power spectrum obeys Bolgiano-Obukhov statistics in 2-D, but Kolmogorov statistics in 3-D. Preliminary results from 3-D reacting bubble calculations are shown, and their implications for ignition are discussed.

Flame instabilities play a dominant role in accelerating the burning front to a large fraction of the speed of sound in a Type Ia supernova. We present a three-dimensional numerical simulation of a Rayleigh-Taylor unstable carbon flame, following its evolution through the transition to turbulence. A low Mach number hydrodynamics method is used, freeing us from the harsh timestep restrictions imposed by sound waves. We fully resolve the thermal structure of the flame and its reaction zone, eliminating the need for a flame model. A single density is considered, $1.5\times 10^7 \gcc$, and half carbon\slash half oxygen fuel--conditions under which the flame propagated in the flamelet regime in our related two-dimensional study. We compare to a corresponding two-dimensional simulation, and show that while fire-polishing keeps the small features suppressed in two dimensions, turbulence wrinkles the flame on far smaller scales in the three-dimensional case, suggesting that the transition to the distributed burning regime occurs at higher densities in three dimensions. Detailed turbulence diagnostics are provided. We show that the turbulence follows a Kolmogorov spectrum and is highly anisotropic on the large scales, with a much larger integral scale in the direction of gravity. Furthermore, we demonstrate that it becomes more isotropic as it cascades down to small scales. Based on the turbulent statistics and the flame properties of our simulation, we compute the Gibson scale. We show the progress of the turbulent flame through a classic combustion regime diagram, indicating that the flame just enters the distributed burning regime near the end of our simulation.

Planar flames are intrinsically unstable in open domains due to the thermal expansion across the burning front--the Landau-Darrieus instability. This instability leads to wrinkling and growth of the flame surface, and corresponding acceleration of the flame, until it is stabilized by cusp formation. We look at the Landau-Darrieus instability for C/O thermonuclear flames at conditions relevant to the late stages of a Type Ia supernova explosion. Two-dimensional direct numerical simulations of both single-mode and multi-mode perturbations using a low Mach number hydrodynamics code are presented. We show the effect of the instability on the flame speed as a function of both the density and domain size, demonstrate the existence of the small scale cutoff to the growth of the instability, and look for the proposed breakdown of the non-linear stabilization at low densities. The effects of curvature on the flame as quantified through measurements of the growth rate and computation of the corresponding Markstein number. While accelerations of a few percent are observed, they are too small to have any direct outcome on the supernova explosion.

A Type Ia supernova explosion likely begins as a nuclear runaway near the center of a carbon-oxygen white dwarf. The outward propagating flame is unstable to the Landau-Darrieus, Rayleigh-Taylor, and Kelvin-Helmholtz instabilities, which serve to accelerate it to a large fraction of the speed of sound. We investigate the Rayleigh-Taylor unstable flame at the transition from the flamelet regime to the distributed-burning regime, around densities of $10^7$ g/cc, through detailed, fully resolved simulations. A low Mach number, adaptive mesh hydrodynamics code is used to achieve the necessary resolution and long time scales. As the density is varied, we see a fundamental change in the character of the burning--at the low end of the density range the Rayleigh-Taylor instability dominates the burning, whereas at the high end the burning suppresses the instability. In all cases, significant acceleration of the flame is observed, limited only by the size of the domain we are able to study. We discuss the implications of these results on the potential for a deflagration to detonation transition.

The Parkes multibeam pulsar survey has unlocked vast areas of the Galactic plane which were previously invisible to earlier low-frequency and less-sensitive surveys. The survey has discovered more than 600 new pulsars so far, including many that are young and exotic. In this paper we report the discovery of 200 pulsars for which we present positional and spin-down parameters, dispersion measures, flux densities and pulse profiles. A large number of these new pulsars are young and energetic, and we review possible associations of $\gamma$-ray sources with the sample of about 1300 pulsars for which timing solutions are known. Based on a statistical analysis, we estimate that about $19\pm6$ associations are genuine. The survey has also discovered 12 pulsars with spin properties similar to those of the Vela pulsar, nearly doubling the known population of such neutron stars. Studying the properties of all known `Vela-like' pulsars, we find their radio luminosities to be similar to normal pulsars, implying that they are very inefficient radio sources. Finally, we review the use of the newly discovered pulsars as Galactic probes and discuss the implications of the new NE2001 Galactic electron density model for the determination of pulsar distances and luminosities.

We have analyzed observations of the Acidalia hemisphere of Mars taken by the Hubble Space Telescope's Near-Infrared Camera Multi-Object Spectrograph (HST/NICMOS) during July of 1997 (Ls = 152 degrees, northern Martian summer). The data consist of images at ~60 km/pixel resolution, using both narrow- and medium-band filters specifically selected to allow us to study the hydration state of the Martian surface. Calibration was performed by comparison to Phobos-2 ISM observations of overlapping regions, and atmospheric gas correction was performed by modeling the atmosphere for each pixel using a line-by-line radiative transfer code coupled with the MOLA altimetry data. Our results indicate the presence of at least three spectrally different large-scale (>1000 km diameter) terrains corresponding to the dark regions of northern Acidalia, the southern hemisphere classical dark terrain, and the classical intermediate terrain adjacent to southern Acidalia. We also identified two other spectrally unique terrains, corresponding to the northern polar ice cap, and to the southern winter polar hood. Comparisons with mineral spectra indicate the possibility of different H2O- or OH-bearing (i.e., hydroxides and/or hydrates) minerals existing both in northern Acidalia and in the nearby intermediate albedo terrain. Hydrated minerals do not appear to be spectrally important components of the southern hemisphere dark terrains imaged by HST in 1997.

The Parkes multibeam pulsar survey is a sensitive survey of a strip of the Galactic plane with $|b|<5\degr$ and $260\degr < l < 50\degr$ at 1374 MHz. Here we report the discovery of 120 new pulsars and subsequent timing observations, primarily using the 76-m Lovell radio telescope at Jodrell Bank. The main features of the sample of 370 published pulsars discovered during the multibeam survey are described. Furthermore, we highlight two pulsars: PSR J1734$-$3333, a young pulsar with the second highest surface magnetic field strength among the known radio pulsars, $B_s = 5.4\times10^{13}$ G, and PSR J1830$-$1135, the second slowest radio pulsar known, with a 6-s period.

The Parkes multibeam pulsar survey is covering a 10deg-wide strip of the southern Galactic plane from l=260deg to l=50deg. It utilizes a 13-beam receiver operating in the 20-cm band on the Parkes 64-m radio telescope and is much more sensitive than any previous large-scale survey. Most of the 608 pulsars discovered so far are relatively distant and many are young, with 37 having a characteristic age of less than 10^5 years. At least one of these is associated with a supernova remnant and four other probable associations are suggested. Several multibeam pulsars have high values of the parameter Edot/d^2 and are within the position error contours of unidentified EGRET gamma-ray sources. These possible associations will be tested with the advent of new gamma-ray telescopes.

We report the discovery of a young and energetic pulsar in the Parkes multibeam survey of the Galactic plane. PSR J1016-5857 has a rotation period of 107 ms and period derivative of 8e-14, implying a characteristic age of 21 kyr and spin-down luminosity of 2.6e36 erg/s. The pulsar is located just outside, and possibly interacting with, the shell supernova remnant G284.3-1.8. Archival X-ray data show a source near the pulsar position which is consistent with emission from a pulsar wind nebula. The pulsar is also located inside the error box of the unidentified EGRET source 3EG J1013-5915, for which it represents a plausible counterpart.

The Parkes multibeam pulsar survey is a sensitive survey of a strip along the Galactic plane with |b|<5 deg and l=260 deg to l=50 deg. It uses a 13-beam receiver on the 64-m Parkes radio telescope, receiving two polarisations per beam over a 288 MHz bandwidth centred on 1374 MHz. Receiver and data acquisition systems are described in some detail. For pulsar periods in the range 0.1 - 2 s and dispersion measures of less than 300 cm^{-3} pc, the nominal limiting flux density of the survey is about 0.2 mJy. At shorter or longer periods or higher dispersions, the sensitivity is reduced. Timing observations are carried out for pulsars discovered in the survey for 12 - 18 months after confirmation to obtain accurate positions, spin parameters, dispersion measures, pulse shapes and mean flux densities. The survey is proving to be extremely successful, with more than 600 pulsars discovered so far. We expect that, when complete, this one survey will come close to finding as many pulsars as all previous pulsar surveys put together. The newly discovered pulsars tend to be young, distant and of high radio luminosity. They will form a valuable sample for studies of pulsar emission properties, the Galactic distribution and evolution of pulsars, and as probes of interstellar medium properties. This paper reports the timing and pulse shape parameters for the first 100 pulsars timed at Parkes, including three pulsars with periods of less than 100 ms which are members of binary systems. These results are briefly compared with the parameters of the previously known population.

We report on the discovery of a binary pulsar, PSR J1740-3052, during the Parkes multibeam survey. Timing observations of the 570-ms pulsar at Jodrell Bank and Parkes show that it is young, with a characteristic age of 350 kyr, and is in a 231-day, highly eccentric orbit with a companion whose mass exceeds 11 M_sun. An accurate position for the pulsar was obtained using the Australia Telescope Compact Array. Near-infrared 2.2-um observations make with the telescopes at the Siding Spring observatory reveal a late-type star coincident with the pulsar position. However, we do not believe that this star is the pulsar's companion, because a typical star of this spectral type and required mass would extend beyond the pulsar's orbit. Furthermore, the measured advance of periastron of the pulsar suggests a more compact companion, for example, a main-sequence star with radius only a few times that of the sun. Such a companion is also more consistent with the small dispersion measure variations seen near periastron. Although we cannot conclusively rule out a black-hole companion, we believe the companion is probably an early B star, making the system similar to the binary PSR J0045-7319.

We report on five binary pulsars discovered in the Parkes multibeam Galactic plane survey. All of the pulsars are old, with characteristic ages 1-11 Gyr, and have relatively small inferred magnetic fields, 5-90e8 G. The orbital periods range from 1.3 to 15 days. As a group these objects differ from the usual low-mass binary pulsars (LMBPs): their spin periods of 9-88 ms are relatively long; their companion masses, 0.2-1.1 Msun, are, in at least some cases, suggestive of CO or more massive white dwarfs; and some of the orbital eccentricities, 1e-5 < e < 0.002, are unexpectedly large. We argue that these observed characteristics reflect binary evolution that is significantly different from that of LMBPs. We also note that intermediate-mass binary pulsars apparently have a smaller scale-height than LMBPs.

The growing level of radio frequency interference (RFI) is a recognized problem for research in radio astronomy. This paper describes an intuitive but powerful RFI cancellation technique that is suitable for radio spectroscopy where time-averages are recorded. An RFI ``reference signal,'' is constructed from the cross power spectrum of the signals from the two polarizations of a reference horn pointed at the source of the RFI signal. The RFI signal paths obey simple phase and amplitude closure relations, which allows computation of the RFI contamination in the astronomical data and the corrections to be applied to the astronomical spectra. Since the method is immune to the effects of multipath scattering in both the astronomy and reference signal channels, ``clean copies'' of the RFI signal are not required. The method could be generalized (1) to interferometer arrays, (2) to correct for scattered solar radiation that causes spectral ``standing waves'' in single-dish spectroscopy, and (3) to pulsar survey and timing applications where a digital correlator plays an important role in broadband pulse dedispersion. Future large radio telescopes, such as the proposed LOFAR and SKA arrays, will require a high degree of RFI suppression and could implement the technique proposed here with the benefit of faster electronics, greater digital precision and higher data rates.

We report on the discovery of PSR J1141-6545, a radio pulsar in an eccentric, relativistic 5-hr binary orbit. The pulsar shows no evidence for being recycled, having pulse period P = 394 ms, characteristic age tau_c = 1.4 x 10^6 yr, and inferred surface magnetic dipole field strength B = 1.3 x 10^12 G. From the mass function and measured rate of periastron advance, we determine the total mass in the system to be (2.300 +/- 0.012) solar masses, assuming that the periastron advance is purely relativistic. Under the same assumption, we constrain the pulsar's mass to be M_p < 1.348 solar masses and the companion's mass to be M_c > 0.968 solar masses (both 99% confidence). Given the total system mass and the distribution of measured neutron star masses, the companion is probably a massive white dwarf which formed prior to the birth of the pulsar. Optical observations can test this hypothesis.

The conference brought together expertise on a range of interference mitigation techniques from CSIRO, Australian and international industry and universities. Key goals were to: enhance the understanding of techniques and their inter-relationship, increase awareness of advanced technologies such as software radios and photonics, and foster a cooperative approach to the development of interference mitigation techniques. The foremost application in mind was the square kilometre array (SKA) and the need to find ways to develop an hierarchical scheme for removing unwanted signals from astronomical data. This paper gives an overview of the topics discussed at the conference and summarises some of the key ideas and results that were presented. Many of the interference challenges we face can be overcome with technical solutions. There will always be a place for regulatory and legal approaches to spectrum management, but too much reliance has been placed on these in the past. The strategy of avoidance has worked well in the past and will continue to be useful in the future, but is now limited by the growth of space based communications systems, which are visible everywhere on Earth. There is no technical solution which is likely to provide a silver bullet. Rather a hierarchical approach will be required, combining a range of techniques. We have just begun exploring technical solutions and their toxicity to astronomical data. We need to continue on this path, determining which techniques are best and how they interact in order to develop a flexible and powerful system for interference suppression.

We report the discovery of two young isolated radio pulsars with very high inferred magnetic fields. PSR J1119-6127 has period P = 0.407 s, and the largest period derivative known among radio pulsars, Pdot = 4.0e-12. Under standard assumptions these parameters imply a characteristic spin-down age of only tau = 1.6 kyr and a surface dipole magnetic field strength of B = 4.1e13 G. We have measured a stationary period-second-derivative for this pulsar, resulting in a braking index of n = 2.91+-0.05. We have also observed a glitch in the rotation of the pulsar, with fractional period change Delta_P/P = -4.4e-9. Archival radio imaging data suggest the presence of a previously uncataloged supernova remnant centered on the pulsar. The second pulsar, PSR J1814-1744, has P = 3.975 s and Pdot = 7.4e-13. These parameters imply tau = 85 kyr, and B = 5.5e13 G, the largest of any known radio pulsar. Both PSR J1119-6127 and PSR J1814-1744 show apparently normal radio emission in a regime of magnetic field strength where some models predict that no emission should occur. Also, PSR J1814-1744 has spin parameters similar to the anomalous X-ray pulsar (AXP) 1E 2259+586, but shows no discernible X-ray emission. If AXPs are isolated, high magnetic field neutron stars (``magnetars''), these results suggest that their unusual attributes are unlikely to be merely a consequence of their very high inferred magnetic fields.

We report on searches of the globular cluster Terzan 5 for low luminosity and accelerated radio pulsars using the 64-m Parkes radio telescope. One new millisecond pulsar, designated PSR J1748-2446C, was discovered, having a period of 8.44ms. Timing measurements using the 76-m Lovell radio telescope at Jodrell Bank show that it is a solitary pulsar and lies close to the core of the cluster. We also present the results of timing measurements which show that the longer-period pulsar PSR J1748-2444 (formerly known as PSR B1744-24B) lies 10 arcmin from the core of the cluster and is unlikely to be associated with the cluster. We conclude that there are further pulsars to be detected in the cluster.

Radio astronomy provides a unique window on the universe, allowing us to study: non-thermal processes (galactic nuclei, quasars, pulsars) at the highest angular resolution using VLBI, with low opacity. It is the most interesting wave band for SETI searches. To date it has yielded 3 Nobel prizes (microwave background, pulsars, gravitational radiation). There are both exciting possibilities and substantial challenges for radio astronomy to remain at the cutting edge over the next 3 decades. New instruments like ALMA and the SKA will open up new science if the challenge of dealing human generated interference can be met. We summarise some of the issues and technological developments that will be essential to the future success of radio astronomy.

Astronomers use the 1612 MHz OH spectral line emission as a unique window on properties of evolved stars, galactic dynamics, and putative proto-planetary disk systems around young stars. In recent years, experiments using this OH line have become more difficult because radio telescopes are very sensitive to transmissions from the GLONASS satellite system. The weak astronomical signals are often undetectable in the presence of these unwanted human generated signals. In this paper we demonstrate that GLONASS narrow band signals may be removed using digital signal processing in a manner that is robust and non-toxic to the weak astronomy signals, without using a reference antenna. We present results using real astronomy data and outline the steps required to implement useful systems on radio telescopes.

We describe the nature of the interference challenges facing radio astronomy in the next decade. These challenges will not be solved by regulation only, negotiation and mitigation will become vital. There is no silver bullet for mitigating against interference. A successful mitigation approach is most likely to be a hierarchical or progressive approach throughout the telescope and signal conditioning and processing systems. We summarise some of the approaches, including adaptive systems.

We have searched 92 unidentified sources from the FIRST and NVSS 1400 MHz radio survey catalogs for radio pulsations at 610 MHz. The selected radio sources are bright, have no identification with extragalactic objects, are point-like and are more than 5% linearly polarized. Our search was sensitive to sub-millisecond pulsations from pulsars with dispersion measures (DMs) less than 500 pc cm-3 in the absence of scattering. We have detected no pulsations from these sources and consider possible effects which might prevent detection. We conclude that as a population, these sources are unlikely to be pulsars.

We report on Keck and {\em Hubble Space Telescope} observations of the binary millisecond pulsar PSR B1855+09. We detect its white-dwarf companion and measure $\mv=25.90\pm0.12$ and $\mi=24.19\pm0.11$ (Vega system). From the reddening-corrected color, $(\mv-\mi)_0=1.06\pm0.21$, we infer a temperature $\Teff=4800\pm800 $K. The white-dwarf mass is known accurately from measurements of the Shapiro delay of the pulsar signal, $\Mcomp=0.258^{+0.028}_{-0.016} \Msun$. Hence, given a cooling model, one can use the measured temperature to determine the cooling age. The main uncertainty in the cooling models for such low-mass white dwarfs is the amount of residual nuclear burning, which is set by the thickness of the hydrogen layer surrounding the helium core. From the properties of similar systems, it has been inferred that helium white dwarfs form with thick hydrogen layers, with mass $\simgt3\times10^{-3} \Msun$, which leads to significant additional heating. This is consistent with expectations from simple evolutionary models of the preceding binary evolution. For PSR B1855+09, though, such models lead to a cooling age of $\sim10 $Gyr, which is twice the spin-down age of the pulsar. It could be that the spin-down age were incorrect, which would call the standard vacuum dipole braking model into question. For two other pulsar companions, however, ages well over 10 Gyr are inferred, indicating that the problem may lie with the cooling models. There is no age discrepancy for models in which the white dwarfs are formed with thinner hydrogen layers ($\simlt3\times10^{-4} \Msun$).

We report the discovery of two isolated radio pulsars having the largest inferred surface dipole magnetic fields yet seen in the population: 4.1e10^13 G and 5.5e10^13 G. These pulsars show apparently normal radio emission in a regime of magnetic field strength where some models predict no emission should occur. They have spin parameters and magnetic fields similar to those of some magnetar candidates, but exhibit very different radiative properties. This demonstrates that if the putative magnetars are indeed isolated neutron stars, their unusual attributes cannot be only a consequence of their large inferred magnetic fields.

A high-frequency survey of the Galactic plane for radio pulsars is in progress, using the multibeam receiver on the 64-m Parkes radiotelescope. We describe the survey motivations, the observing plan and the inital results. The survey is discovering many pulsars, more than 500 so far. Eight of the new pulsars are binary, one with a massive companion. At least eight are young, with characteristic ages of less than 100 kyr. Two of these (Kaspi et al, this Conference) have surface dipole magnetic field strengths greater than any other known radio pulsar.

The Square Kilometre Array (SKA) is a proposed next generation radio telescope. Between now and 2005 this project is in a technology development and prototyping phase, with construction likely to begin in $\sim 2010$. This paper describes what the SKA may be like, its key features, the motivation for building it and where you can access more details about it. Its is important to see any new facility in context, so other new facilities are also discussed. Avenues for future extensibility of the SKA other telescopes are covered, with some emphasis on multiple beam systems. Some suggestions for useful pulsar experiments and pulsar searches strategies are summarised. A conclusion is that the SKA may not be the most cost effective way to search for pulsars and that a 128 beam receiver system on an Arecibo like telescope working between 1 and 5 GHz may be a more cost effective approach.

The Parkes multibeam pulsar survey began in 1997 and is now about 50% complete. It has discovered more than 400 new pulsars so far, including a number of young, high magnetic field, and relativistic binary pulsars. Early results, descriptions of the survey and follow up timing programs can be found in papers by Lyne et al. (1999 MNRAS in press, astro-ph/9911313), Camilo et al. (astro-ph/9911185), and Manchester et al. (astro-ph/9911319). This paper describes the data release policy and how you can gain access to the raw data and details on the pulsars discovered.

Measurement of accurate positions, pulse periods and period derivatives is an essential follow-up to any pulsar survey. The procedures being used to obtain timing parameters for the pulsars discovered in the Parkes multibeam pulsar survey are described. Completed solutions have been obtained so far for about 80 pulsars. They show that the survey is preferentially finding pulsars with higher than average surface dipole magnetic fields. Eight pulsars have been shown to be members of binary systems and some of the more interesting results relating to these are presented.

We are undertaking a high-frequency survey of the Galactic plane for radio pulsars, using the 13-element multibeam receiver on the 64-m Parkes radio telescope. We describe briefly the survey system and some of the initial results. PSR J1811-1736, one of the first pulsars discovered with this system, has a rotation period of 104 ms. Subsequent timing observations using the 76-m radio telescope at Jodrell Bank show that it is in an 18.8-day, highly-eccentric binary orbit. We have measured the rate of advance of periastron which indicates a total system mass of 2.6 +- 0.9 Msun, and the minimum companion mass is about 0.7 Msun. This, the high orbital eccentricity and the recycled nature of the pulsar suggests that this system is composed of two neutron stars, only the fourth or fifth such system known in the disk of the Galaxy.

We have searched 92 unidentified sources from the FIRST and NVSS radio catalogs for pulsations at 610 MHz. The selected sources are bright, have no identifications, are unresolved and have significant linear polarization. Our search was sensitive to sub-millisecond pulsations from pulsars with a large range of dispersion measures. We have detected no radio pulsations from these sources and conclude that they are unlikely to be a population of previously undetected pulsars.

The Parkes multibeam pulsar survey uses a 13-element receiver operating at a wavelength of 20 cm to survey the inner Galactic plane with remarkable sensitivity. To date we have collected and analyzed data from 45% of the survey region (|b| < 5 deg.; 260 deg. < l < 50 deg.), and have discovered 440 pulsars, in addition to re-detecting 190 previously known ones. Most of the newly discovered pulsars are at great distances, as inferred from a median dispersion measure (DM) of 400 pc/cm^3.

The Parkes multibeam pulsar survey is a high-frequency, fast-sampled survey of the Galactic Plane, expected to discover at least 500 new pulsars. To date, over 200 pulsars have been found, including several young pulsars and at least one with a very high magnetic field. Seven of the new stars are in binary systems; this number includes one probable double-neutron-star system, and one pulsar with an extremely massive companion.

Young radio pulsars shed vast amounts of rotational energy, sometimes as high as 100,000 times the total energy loss rate from the sun. The wide range of phenomena resulting from this energy loss include: glitches, timing noise, jets, bow shocks, bullets and plerions and are be reviewed from an observational perspective. Past and proposed surveys for young radio pulsars are summarised along with pulsar birth velocities and associations with supernova remnants. There are now 4 radio pulsars with measured braking indices. The resulting constraints on the evolution of young radio pulsars are discussed in light of the presently observed population of pulsars. Observations at optical, X-ray and gamma ray energies which provide a unique opportunity to study the emission and magnetospheric processes are described briefly. The status of pulsar birth velocities and supernova remnant associations are summarised.

General relativity asserts that: energy and momentum conservation laws are valid, preferred frames do not exist, and the strong equivalence principle is obeyed. In this paper recent progress in testing these important principles using millisecond pulsars is summarised.

We present timing observations of four millisecond pulsars, using data obtained over three years at the ATNF Parkes and NRAL Jodrell Bank radio telescopes. Astrometric, spin, and binary parameters are updated, and substantially improved for three pulsars, PSRs J0613-0200, J1045-4509 and J1643-1224. We have measured the time variation of the projected semi-major axis of the PSR J0437-4715 orbit due to its proper motion, and use it to constrain the inclination of the orbit and the mass of the companion. Some evidence is found for changes in the dispersion measures of PSRs J1045-4509 and J1643-1224. Limits are placed on the existence of planetary mass companions, ruling out companions with masses and orbits similar to the terrestrial planets of the solar system for eight pulsars.

The observational evidence for jets and phenomena arising from rotation powered radio pulsars is reviewed, including many recent and exciting discoveries at X-ray wavelengths. The well studied jets of the Crab pulsar are summarised, including recent results from the HST. The evolutionary links between the known binary radio pulsars and jets sources in X-ray binaries are discussed.

The motivation for radio pulsar timing and its basic principles are reviewed. Present and future radio timing techniques and hardware are summarised and compared. The array of present timing programmes and their scientific goals are collated and described. Recent results and future prospects are discussed, with emphasis on multi-wavelength techniques where appropriate. Timing of radio pulsars at other wavelengths is summarised along with the provision of contemporary ephemerides for timing and searches at other wavelengths.

General relativity predicts that energy and momentum conservation laws hold and that preferred frames do not exist. The parametrised post-Newtonian formalism (PPN) phenomenologically quantifies possible deviations from general relativity. The PPN parameter alpha_3 (which identically vanishes in general relativity) plays a dual role in that it is associated both with a violation of the momentum conservation law, and with the existence of a preferred frame. By considering the effects of alpha_3 neq 0 in certain binary pulsar systems, it is shown that alpha_3 < 2.2 x 10^-20 (90% CL). This limit improves on previous results by several orders of magnitude, and shows that pulsar tests of alpha_3 rank (together with Hughes-Drever-type tests of local Lorentz invariance) among the most precise null experiments of physics.

Gravity being a long-range force, one might {\it a priori} expect the Universe's global matter distribution to select a preferred rest frame for local gravitational physics. The phenomenology of preferred-frame effects, in the strong-gravitational field context of binary pulsars, is described by two parameters $\hat{\alpha}_1$ and $\hat{\alpha}_2$. These parameters vanish identically in general relativity, and reduce, in the weak-field limit, to the two parametrized post-Newtonian (PPN) parameters ${\alpha}_1$ and ${\alpha}_2$. We derive a limit of $|\hat{\alpha}_1| < 1.7\times 10^{-4}$ (90\%~C.L.) using the very low eccentricity binary pulsar PSR J2317+1439, improving by a factor of 3 on previous limits.

We demonstrate how measuring orbital period derivatives can lead to more accurate distance estimates and transverse velocities for some nearby binary pulsars. In many cases this method will estimate distances more accurately than is possible by annual parallax, as the relative error decreases as t^-5/2. Unfortunately, distance uncertainties limit the degree to which nearby relativistic binary pulsars can be used for testing the general relativistic prediction of orbital period decay to a few percent. Nevertheless, the measured orbital period derivative of PSR B1534+12 agrees within the observational uncertainties with that predicted by general relativity if the proper-motion contribution is accounted for.

Some theories of gravity predict the existence of preferred-frame effects and violations of conservation of energy and momentum. General relativity predicts no such effects. In the parameterised post-Newtonian (PPN) formalism, the parameter, \al3 $\equiv 0$ if these effects do not exist. The period derivatives ($\dot{P}$) of millisecond pulsars (MSPs) are used to more tightly constrain these effects by showing that $|$\al3$|$ $< 5 \times 10^{-16}$.

We report on optical observations of the low-mass binary millisecond pulsar systems J0034-0534 and J2145-0750. A faint (I=23.5) object was found to be coincident with the timing position of PSR J2145-0750. While a galaxy or distant main-sequence star cannot be ruled out, its magnitude is consistent with an ancient white dwarf, as expected from evolutionary models. For PSR J0034-0534 no objects were detected to a limiting magnitude of R=25.0, suggesting that the white dwarf in this system is cold. Using white dwarf cooling models, the limit on the magnitude of the PSR J0034-0534 companion suggests that at birth the pulsar in this system may have rotated with a period as short as 0.6 ms. These observations provide further evidence that the magnetic fields of millisecond pulsars do not decay on time scales shorter than 1 Gyr.

Binary radio pulsars are superb tools for mapping binary orbits, because of the precision of the pulse timing method (Taylor and Weisberg 1989). To date, all orbital parameters for binary pulsars have been derived from observations of the pulsar alone. We present the first observations of the radial velocity variations due to the binary motion of a companion to a radio pulsar. Our results demonstrate that the companion to the Small Magellanic Cloud pulsar PSR J0045-7319 is the B1V star identified by Kaspi et al. (1994). The mass ratio of the system is 6.3 +/- 1.2, which, for a neutron star mass of 1.4 Mo, implies a mass of 8.8 +/- 1.8 Mo for the companion, consistent with the mass expected for a B1V star. The inclination angle for the binary system is therefore 44 +/- 5 degrees, and the projected rotational velocity of the companion is 113 +/- 10 km/s. The heliocentric radial velocity of the binary system is consistent with that of other stars and gas in the same region of the Small Magellanic Cloud.