The Kuiper Belt is a distant region of the Solar System. On 1 January 2019, the New Horizons spacecraft flew close to (486958) 2014 MU69, a Cold Classical Kuiper Belt Object, a class of objects that have never been heated by the Sun and are therefore well preserved since their formation. Here we describe initial results from these encounter observations. MU69 is a bi-lobed contact binary with a flattened shape, discrete geological units, and noticeable albedo heterogeneity. However, there is little surface color and compositional heterogeneity. No evidence for satellites, ring or dust structures, gas coma, or solar wind interactions was detected. By origin MU69 appears consistent with pebble cloud collapse followed by a low velocity merger of its two lobes.

One of the principal bottlenecks to atmosphere characterisation in the era of all-sky surveys is the availability of fast, autonomous and robust atmospheric retrieval methods. We present a new approach using unsupervised machine learning to generate informed priors for retrieval of exoplanetary atmosphere parameters from transmission spectra. We use principal component analysis (PCA) to efficiently compress the information content of a library of transmission spectra forward models generated using the PLATON package. We then apply a $k$-means clustering algorithm in PCA space to segregate the library into discrete classes. We show that our classifier is almost always able to instantaneously place a previously unseen spectrum into the correct class, for low-to-moderate spectral resolutions, $R$, in the range $R~=~30-300$ and noise levels up to $10$~per~cent of the peak-to-trough spectrum amplitude. The distribution of physical parameters for all members of the class therefore provides an informed prior for standard retrieval methods such as nested sampling. We benchmark our informed-prior approach against a standard uniform-prior nested sampler, finding that our approach is up to a factor two faster, with negligible reduction in accuracy. We demonstrate the application of this method to existing and near-future observatories, and show that it is suitable for real-world application. Our general approach is not specific to transmission spectroscopy and should be more widely applicable to cases that involve repetitive fitting of trusted high-dimensional models to large data catalogues, including beyond exoplanetary science.

The large number of new planets expected from wide-area transit surveys means that follow-up transmission spectroscopy studies of their atmospheres will be limited by the availability of telescope assets. We argue that telescopes covering a broad range of apertures will be required, with even 1m-class instruments providing a potentially important contribution. Survey strategies that employ automated target selection will enable robust population studies. As part of such a strategy, we propose a decision metric to pair the best target to the most suitable telescope, and demonstrate its effectiveness even when only primary transit observables are available. Transmission spectroscopy target selection need not therefore be impeded by the bottle-neck of requiring prior follow-up observations to determine the planet mass. The decision metric can be easily deployed within a distributed heterogeneous network of telescopes equipped to undertake either broadband photometry or spectroscopy. We show how the metric can be used either to optimise the observing strategy for a given telescope (e.g. choice of filter) or to enable the selection of the best telescope to optimise the overall sample size. Our decision metric can also provide the basis for a selection function to help evaluate the statistical completeness of follow-up transmission spectroscopy datasets. Finally, we validate our metric by comparing its ranked set of targets against lists of planets that have had their atmospheres successfully probed, and against some existing prioritised exoplanet lists.

We consider a PDE system comprising compressible hydrodynamics, flux-limited diffusion radiation transport and chemical ionization kinetics in a cosmologically-expanding universe. Under an operator-split framework, the cosmological hydrodynamics equations are solved through the Piecewise Parabolic Method, as implemented in the Enzo community hydrodynamics code. The remainder of the model, including radiation transport, chemical ionization kinetics, and gas energy feedback, form a stiff coupled PDE system, which we solve using a fully-implicit inexact Newton approach, and which forms the crux of this paper. The inner linear Newton systems are solved using a Schur complement formulation, and employ a multigrid-preconditioned conjugate gradient solver for the inner Schur systems. We describe this approach and provide results on a suite of test problems, demonstrating its accuracy, robustness, and scalability to very large problems.

We use high resolution 2D hydrodynamic simulations to study the formation of spiral substructure in the gaseous disk of a galaxy. The obtained gaseous response is driven by a self-consistent non-axisymmetric potential obtained from an imposed spiral mass distribution. We highlight the importance of ultraharmonic resonances in generating these features. The temporal evolution of the system is followed with the parallel ZEUS-MP code, and we follow the steepening of perturbations induced by the spiral potential until large-scale shocks emerge. These shocks exhibit bifurcations that protrude from the gaseous arms and continue to steepen until new shocks are formed. When the contribution from the spiral potential relative to the axisymmetric background is increased from our default value, spurs protrude from the main arms after several revolutions of the gaseous disk. Such spurs overlap on top of the aforementioned shocks. These results support the hypothesis that a complicated gaseous response can coexist with an orderly spiral potential term, in the sense that the underlying background potential can be smooth yet drive a gaseous response that is far more spatially complex.

The Pluto Energetic Particle Spectrometer Science Investigation (PEPSSI) comprises the hardware and accompanying science investigation on the New Horizons spacecraft to measure pick-up ions from Pluto's outgassing atmosphere. To the extent that Pluto retains its characteristics similar to those of a "heavy comet" as detected in stellar occultations since the early 1980s, these measurements will characterize the neutral atmosphere of Pluto while providing a consistency check on the atmospheric escape rate at the encounter epoch with that deduced from the atmospheric structure at lower altitudes by the ALICE, REX, and SWAP experiments on New Horizons. In addition, PEPSSI will characterize any extended ionosphere and solar wind interaction while also characterizing the energetic particle environment of Pluto, Charon, and their associated system. First proposed for development for the Pluto Express mission in September 1993, what became the PEPSSI instrument went through a number of development stages to meet the requirements of such an instrument for a mission to Pluto while minimizing the required spacecraft resources. The PEPSSI instrument provides for measurements of ions (with compositional information) and electrons from 10s of keV to ~1 MeV in a 120 deg x 12 deg fan-shaped beam in six sectors for 1.5 kg and ~2.5 W.

The LOng-Range Reconnaissance Imager (LORRI) is the high resolution imaging instrument for the New Horizons mission to Pluto, its giant satellite Charon, its small moons Nix and Hydra, and the Kuiper Belt, which is the vast region of icy bodies extending roughly from Neptune's orbit out to 50 astronomical units (AU). New Horizons launched on January 19, 2006 as the inaugural mission in NASA's New Frontiers program. LORRI is a narrow angle (field of view=0.29 deg), high resolution (4.95 microrad pixels), Ritchey-Chretien telescope with a 20.8 cm diameter primary mirror, a focal length of 263 cm, and a three lens field-flattening assembly. A 1024 x 1024 pixel (optically active region), thinned, backside-illuminated charge-coupled device (CCD) detector is used in the focal plane unit and is operated in frame transfer mode. LORRI provides panchromatic imaging over a bandpass that extends approximately from 350 nm to 850 nm. LORRI operates in an extreme thermal environment, situated inside the warm spacecraft with a large, open aperture viewing cold space. LORRI has a silicon carbide optical system, designed to maintain focus over the operating temperature range without a focus adjustment mechanism. Moreover, the spacecraft is thruster-stabilized without reaction wheels, placing stringent limits on the available exposure time and the optical throughput needed to satisfy the measurement requirements.

The overwhelming evidence that the core collapse supernova mechanism is inherently multidimensional, the complexity of the physical processes involved, and the increasing evidence from simulations that the explosion is marginal presents great computational challenges for the realistic modeling of this event, particularly in 3 spatial dimensions. We have developed a code which is scalable to computations in 3 dimensions which couples PPM Lagrangian with remap hydrodynamics [1], multigroup, flux-limited diffusion neutrino transport [2], with many improvements), and a nuclear network [3]. The neutrino transport is performed in a ray-by-ray plus approximation wherein all the lateral effects of neutrinos are included (e.g., pressure, velocity corrections, advection) except the transport. A moving radial grid option permits the evolution to be carried out from initial core collapse with only modest demands on the number of radial zones. The inner part of the core is evolved after collapse along with the rest of the core and mantle by subcycling the lateral evolution near the center as demanded by the small Courant times. We present results of 2-D simulations of a symmetric and an asymmetric collapse of both a 15 and an 11 M progenitor. In each of these simulations we have discovered that once the oxygen rich material reaches the shock there is a synergistic interplay between the reduced ram pressure, the energy released by the burning of the shock heated oxygen rich material, and the neutrino energy deposition which leads to a revival of the shock and an explosion.

This paper describes ZEUS-MP, a multi-physics, massively parallel, message- passing implementation of the ZEUS code. ZEUS-MP differs significantly from the ZEUS-2D code, the ZEUS-3D code, and an early "version 1" of ZEUS-MP distributed publicly in 1999. ZEUS-MP offers an MHD algorithm better suited for multidimensional flows than the ZEUS-2D module by virtue of modifications to the Method of Characteristics scheme first suggested by Hawley and Stone (1995), and is shown to compare quite favorably to the TVD scheme described by Ryu et. al (1998). ZEUS-MP is the first publicly-available ZEUS code to allow the advection of multiple chemical (or nuclear) species. Radiation hydrodynamic simulations are enabled via an implicit flux-limited radiation diffusion (FLD) module. The hydrodynamic, MHD, and FLD modules may be used in one, two, or three space dimensions. Self gravity may be included either through the assumption of a GM/r potential or a solution of Poisson's equation using one of three linear solver packages (conjugate-gradient, multigrid, and FFT) provided for that purpose. Point-mass potentials are also supported. Because ZEUS-MP is designed for simulations on parallel computing platforms, considerable attention is paid to the parallel performance characteristics of each module. Strong-scaling tests involving pure hydrodynamics (with and without self-gravity), MHD, and RHD are performed in which large problems (256^3 zones) are distributed among as many as 1024 processors of an IBM SP3. Parallel efficiency is a strong function of the amount of communication required between processors in a given algorithm, but all modules are shown to scale well on up to 1024 processors for the chosen fixed problem size.

This paper presents a new code for performing multidimensional radiation hydrodynamic (RHD) simulations on parallel computers involving anisotropic radiation fields and nonequilibrium effects. The radiation evolution modules described here encapsulate the physics provided by the serial algorithm of Stone et. al (ApJSupp, vol 80, pp. 819-845), but add new functionality with regard to physics and numerics. Physics enhancments include the addition of time dependence to the computation of the variable tensor Eddington factor (VTEF) closure term, and a matter-radiation coupling scheme which is particularly robust for nonequilibrium problems. Numerical highlights include a discussion of how our code is implemented for parallel execution and a description of our scalable linear solver module. We present a suite of numerical tests from which the virtues and vices of our method may be gleaned; these include nonequilibrium Marshak waves, 2-D "shadow" tests showing the one-sided illumination of an opaque cloud, and full RHD+VTEF simulations of radiating shocks. We conclude that radiation moment solutions closed with variable tensor Eddington factors show a dramatic qualitative improvement over results obtained with flux-limited diffusion, and further that this approach has a bright future in the context of parallel RHD simulations in astrophysics.

The cluster RX J1347.5-1145, the most luminous cluster in the X-ray wavelengths, was imaged with the newly installed Space Telescope Imaging Spectrograph (STIS) on-board HST. Its relatively high redshift (0.451) and luminosity indicate that this is one of the most massive of all known clusters. The STIS images unambiguously show several arcs in the cluster. The largest two arcs (> 5 arcsec in length) are symmetrically situated on opposite sides of the cluster, at a distance of ~ 35 arcsec from the central galaxy. The STIS images also show approximately 100 faint galaxies within the radius of the arcs whose combined luminosity is ~ 4 x 10^11 Lsun. We also present ground-based spectroscopic observations of the northern arc which show one clear emission line at 6730 A, which is consistent with an identification as [OII] 3727 A, implying a redshift of 0.81 for this arc. The southern arc shows a faint continuum but no emission features. The surface mass within the radius of the arcs (240 kpc), as derived from the gravitational lensing, is 6.3 x 10^14 Msun. The resultant mass-to-light ratio of ~1200 is higher than what is seen in many clusters but smaller than the value recently derived for some `dark' X-ray clusters (Hattori et al. 1997). The total surface mass derived from the X-ray flux within the radius of the arcs is ~2.1 - 6.8 x 10^14 Msun, which implies that the ratio of the gravitational to the X-ray mass is ~1 to 3. The surface GAS mass within this radius is ~3.5 x 10^13 Msun, which implies that at least 6% of the total mass within this region is baryonic.

New data imply that the average velocity of radio pulsars is large \cite{hla93}. Under the assumption that these data imply that a pulsar is born with an ``intrinsic'' kick, we investigate whether such kicks can be a consequence of asymmetrical stellar collapse and explosion. We calculate the gravitational wave (GW) signature of such asymmetries due to anisotropic neutrino radiation and mass motions. We predict that any recoils imparted to the neutron star at birth will result in a gravitational wave strain, h$^{TT}_{zz}$, that does not go to zero with time. Hence, there may be ``memory'' \cite{bt87} in the gravitational waveform from a protoneutron star that is correlated with its recoil and neutrino emissions.

It is now believed that pulsars comprise the fastest population of stars in the galaxy. With inferred mean, root-mean-square, and maximum 3-D pulsar speeds of $\sim$300-500 km/s, $\sim$500 km/s, and $\sim$2000 km/s, respectively, the question of the origin of such singular proper motions becomes acute. What mechanism can account for speeds that range from zero to twice the galactic escape velocity? We speculate that a major vector component of a neutron star's proper motion comes from the hydrodynamic recoil of the nascent neutron star during the supernova explosion in which it is born. Recently, theorists have shown that asymmetries and instabilities are a natural aspect of supernova dynamics. In this paper, we highlight two phenomena: 1) the ``Brownian-like'' stochastic motion of the core in response to the convective ``boiling'' of the mantle of the protoneutron star during the post-bounce, pre-explosion accretion phase, and 2) the asymmetrical bounce and explosion of an aspherically collapsing Chandrasekhar core.In principle,either phenomenon can leave the young neutron star with a speed of hundreds of kilometers per second.

We investigate in this paper the core-collapse supernova explosion mechanism in both one and two dimensions. We verify the usefulness of neutrino-driven overturn (``convection'') between the shock and the neutrinosphere in igniting the supernova explosion. The 2-D simulation of the core of a 15$\mdot$ star that we present here indicates that the breaking of spherical symmetry may be central to the explosion itself and that a multitude of bent and broken fingers is a common feature of the ejecta. In 2-D, the steady-state shock radius is larger than its value in 1-D by 30\%--100\%. The critical condition for explosion is thereby relaxed. Since supernovae are driven by neutrino heating, they are coronal phenomena, akin to winds, though initially bounded by an accretion tamp. Neutrino energy is pumped into the supernova during the shock's propagation through the inner many thousands of kilometers and not instantaneously. Curiously, just after the explosion is triggered, the matter that will eventually be ejected is still bound. The high-speed fingers that emerge from the core seem a natural explanation for the nickel bullets seen in 87A and the shrapnel inferred in some supernova remnants. Within 100 ms of the explosion, a strong, neutrino-driven wind is blowing outward from the protoneutron star.