`Oumuamua (I1 2017) was the first macroscopic ($l\sim100\,{\rm m}$) body observed to traverse the inner solar system on an unbound hyperbolic orbit. Its light curve displayed strong periodic variation, and it showed no hint of a coma or emission from molecular outgassing. Astrometric measurements indicate that 'Oumuamua experienced non-gravitational acceleration on its outbound trajectory, but energy balance arguments indicate this acceleration is inconsistent with a water ice sublimation-driven jet of the type exhibited by solar system comets. We show that all of `Oumaumua's observed properties can be explained if it contained a significant fraction of molecular hydrogen (H$_{2}$) ice. H$_{2}$ sublimation at a rate proportional to the incident solar flux generates a surface-covering jet that reproduces the observed acceleration. Mass wasting from sublimation leads to monotonic increase in the body axis ratio, explaining `Oumuamua's shape. Back-tracing `Oumuamua's trajectory through the Solar System permits calculation of its mass and aspect ratio prior to encountering the Sun. We show that H$_{2}$-rich bodies plausibly form in the coldest dense cores of Giant Molecular Clouds, where number densities are of order $n\sim10^5$, and temperatures approach the $T=3\,{\rm K}$ background. Post-formation exposure to galactic cosmic rays implies a $\tau \sim 100$ Myr age, explaining the kinematics of `Oumuamua's inbound trajectory.

We report the clear detection of C$_2$ and of abundant NH$_2$ in the first prominently active interstellar comet, 2I/Borisov. We observed 2I on three nights in November 2019 at optical wavelengths 4800--9300 \AA with the Multi-Unit Spectroscopic Explorer (MUSE) integral-field spectrograph on the ESO/Very Large Telescope. These data, together with observations close in time from both 0.6-m TRAPPIST telescopes, provide constraints on the production rates of species of gas in 2I's coma. From the MUSE detection on all epochs of several bands of the optical emission of the C$_2$ Swan system, a rich emission spectrum of NH$_2$ with many highly visible bands, and the red (1-0) bandhead of CN, together with violet CN detections by TRAPPIST, we infer production rates of $Q$(C$_2$) = $1.1\times10^{24}$ mol s$^{-1}$, $Q$(NH$_2$) = $4.8\times10^{24}$ mol s$^{-1}$ and $Q$(CN) = $(1.8\pm0.2)\times 10^{24}$ mol s$^{-1}$. In late November at 2.03~au, 2I had a production ratio of C$_2$/CN$=0.61$, only barely carbon-chain depleted, in contrast to earlier reports measured further from the Sun of strong carbon-chain depletion. Thus, 2I has shown evolution in its C$_2$ production rate: a parent molecule reservoir has started sublimating. At $Q$(NH$_2$)/$Q$(CN) = 2.7, this second interstellar object is enriched in NH$_2$, relative to the known Solar System sample.

Recently Squire & Hopkins showed that charged dust grains moving through magnetized gas under the influence of any external force (e.g. radiation pressure, gravity) are subject to a spectrum of instabilities. Qualitatively distinct instability families are associated with different Alfvenic or magnetosonic waves and drift or gyro motion. We present a suite of simulations exploring these instabilities, for grains in a homogeneous medium subject to an external acceleration. We vary parameters such as the ratio of Lorentz-to-drag forces on dust, plasma $\beta$, size scale, and acceleration. All regimes studied drive turbulent motions and dust-to-gas fluctuations in the saturated state, can rapidly amplify magnetic fields into equipartition with velocity fluctuations, and produce instabilities that persist indefinitely (despite random grain motions). Different parameters produce diverse morphologies and qualitatively different features in dust, but the saturated gas state can be broadly characterized as anisotropic magnetosonic or Alfvenic turbulence. Quasi-linear theory can qualitatively predict the gas turbulent properties. Turbulence grows from small to large scales, and larger-scale modes usually drive more vigorous gas turbulence, but dust velocity and density fluctuations are more complicated. In many regimes, dust forms structures (clumps, filaments, sheets) that reach extreme over-densities (up to $\gg 10^{9}$ times mean), and exhibit substantial sub-structure even in nearly-incompressible gas. These can be even more prominent at lower dust-to-gas ratios. In other regimes, dust self-excites scattering via magnetic fluctuations that isotropize and amplify dust velocities, producing fast, diffusive dust motions.

We investigate the vorticity-preserving properties of the compressible, second-order residual-based scheme, "RBV2". The scheme has been extensively tested on hydrodynamical problems, and has been shown to exhibit remarkably accurate results on the propagation of inviscid compressible vortices, airfoil-vortex interactions on a curvilinear mesh, vortex mergers in an astrophysical accretion disk, and the establishment of a two-dimensional inverse cascade in high-resolution turbulent simulations. Here, we demonstrate that RBV2 sustains the analytic solution for a one-dimensional shear flow. We assess the fidelity by which the algorithm maintains a skewed shear flow, and present convergence tests to quantify the magnitude of the expected numerical dispersion. We propose an adjustment to the dissipation in the algorithm that retains the vorticity-preserving qualities, and accurately incorporates external body forces, and demonstrate that it indefinitely maintains a steady-state hydrostatic equilibrium between a generic acceleration and a density gradient. We present a novel numerical assessment of vorticity preservation for discrete-wavenumber vortical modes up to the Nyquist wavenumber. We apply this assessment to RBV2 in order to quantify the extent to which the scheme preserves vorticity for the full Euler equations. We find that RBV2 perfectly preserves vorticity for modes with symmetric wavenumbers, i.e., $k_x=k_y$, and that the error increases with asymmetry. We simulate the dynamical interaction of vortices in a protoplanetary disk to demonstrate the utility of the updated scheme for rendering astrophysical flows replete with vortices and turbulence. We conclude that RBV2 is a competitive treatment for evolving vorticity-dominated astrophysical flows with minimal dissipation.

We show that the $P\sim8\,{\rm h}$ photometric period and the astrometrically measured $A_{\rm ng}\sim2.5\times10^{-4}\,{\rm cm\,s^{-2}}$ non-gravitational acceleration (at $r\sim1.4\,{\rm AU}$) of the interstellar object 1I/2017 (`Oumuamua) can be explained by a nozzle-like venting of volatiles whose activity migrated to track the sub-solar location on the object's surface. Adopting the assumption that `Oumuamua was an elongated $a\times b \times c$ ellipsoid, this model produces a pendulum-like rotation of the body and implies a long semi-axis $a\sim 5A_{\rm ng}P^2/4\pi^2 \sim 260\,{\rm m}$. This scale agrees with the independent estimates of `Oumuamua's size that stem from its measured brightness, assuming an albedo of $p\sim0.1$, appropriate to ices that have undergone long-duration exposure to the interstellar cosmic ray flux. Using ray-tracing, we generate light curves for ellipsoidal bodies that are subject to both physically consistent sub-solar torques and to the time-varying geometry of the Sun-Earth-`Oumuamua configuration. Our synthetic light curves display variations from chaotic tumbling and changing cross-sectional illumination that are consistent with the observations, while avoiding significant secular changes in the photometric periodicity. If our model is correct, `Oumuamua experienced mass loss that wasted $\sim 10\%$ of its total mass during the $\sim 100\,{\rm d}$ span of its encounter with the inner Solar System and had an icy composition with a very low $[{\rm C}/{\rm O}]\lesssim 0.003$. Our interpretation of `Oumuamua's behavior is consistent with the hypothesis that it was ejected from either the outer regions of a planetesimal disk after an encounter with an embedded $M_{\rm p} \sim M_{\rm Nep}$ planet or from an exo-Oort cloud.

We investigate, for the first time, the nonlinear evolution of the magnetized "resonant drag instabilities" (RDIs). We explore magnetohydrodynamic (MHD) simulations of gas mixed with (uniform) dust grains subject to Lorentz and drag forces, using the GIZMO code. The magnetized RDIs exhibit fundamentally different behaviour than the purely acoustic RDIs. The dust organizes into coherent structures and the system exhibits strong dust-gas separation. In the linear and early nonlinear regime, the growth rates agree with linear theory and the dust self-organizes into two-dimensional planes or "sheets." Eventually the gas develops fully nonlinear, saturated Alfv\'enic and compressible fast-mode turbulence, which fills the under-dense regions with a small amount of dust, and drives a dynamo which saturates at equipartition of kinetic and magnetic energy. The dust density fluctuations exhibit significant non-Gaussianity, and the power spectrum is strongly weighted towards the largest (box-scale) modes. The saturation level can be understood via quasi-linear theory, as the forcing and energy input via the instabilities becomes comparable to saturated tension forces and dissipation in turbulence. The magnetized simulation presented here is just one case; it is likely that the magnetic RDIs can take many forms in different parts of parameter space.

We report the first confirmation of a hot Jupiter discovered by the Transiting Exoplanet Survey Satellite (TESS) mission: HD 202772A b. The transit signal was detected in the data from TESS Sector 1, and was confirmed to be of planetary origin through radial-velocity measurements. HD 202772A b is orbiting a mildly evolved star with a period of 3.3 days. With an apparent magnitude of V = 8.3, the star is among the brightest known to host a hot Jupiter. Based on the 27days of TESS photometry, and radial velocity data from the CHIRON and HARPS spectrographs, the planet has a mass of 1.008+/-0.074 M_J and radius of 1.562+/-0.053 R_J , making it an inflated gas giant. HD 202772A b is a rare example of a transiting hot Jupiter around a quickly evolving star. It is also one of the most strongly irradiated hot Jupiters currently known.

A rapid accumulation of observations and interpretation have followed in the wake of 1I `Oumuamua's passage through the inner Solar System. We briefly outline the consequences that this first detection of an interstellar asteroid implies for the planet-forming process, and we assess the near-term prospects for detecting and observing (both remotely and in situ) future Solar System visitors of this type. Drawing on detailed heat-transfer calculations that take both `Oumuamua's unusual shape and its chaotic tumbling into account, we affirm that the lack of a detectable coma in deep images of the object very likely arises from the presence of a radiation-modified coating of high molecular weight material (rather than a refractory bulk composition). Assuming that `Oumuamua is a typical representative of a larger population with a kinematic distribution similar to Population I stars in the local galactic neighborhood, we calculate expected arrival rates, impact parameters and velocities of similar objects and assess their prospects for detection using operational and forthcoming facilities. Using `Oumuamua as a proof-of-concept, we assess the prospects for missions that intercept interstellar objects (ISOs) using conventional chemical propulsion. Using a "launch on detection" paradigm, we estimate wait times of order 10 years between favorable mission opportunities with the detection capabilities of the Large-Scale Synoptic Survey Telescope (LSST), a figure that will be refined as the population of interstellar asteroids becomes observationally better constrained.

Vortices, turbulence, and unsteady non-laminar flows are likely both prominent and dynamically important features of astrophysical disks. Such strongly nonlinear phenomena are often difficult, however, to simulate accurately, and are generally amenable to analytic treatment only in idealized form. In this paper, we explore the evolution of compressible two-dimensional flows using an implicit dual-time hydrodynamical scheme that strictly conserves vorticity (if applied to simulate inviscid flows for which Kelvin's Circulation Theorem is applicable). The algorithm is based on the work of Lerat, Falissard & Side (2007), who proposed it in the context of terrestrial applications such as the blade-vortex interactions generated by helicopter rotors. We present several tests of Lerat et al.'s vorticity-preserving approach, which we have implemented to second-order accuracy, providing side-by-side comparisons with other algorithms that are frequently used in protostellar disk simulations. The comparison codes include one based on explicit, second-order van-Leer advection, one based on spectral methods, and another that implements a higher-order Godunov solver. Our results suggest that Lerat et al's algorithm will be useful for simulations of astrophysical environments in which vortices play a dynamical role, and where strong shocks are not expected.

We compare the average photospheric current helicity $H_c$, photospheric twist parameter $\alpha$ (a well-known proxy for the full relative magnetic helicity), and subsurface kinetic helicity $H_k$ for 194 active regions observed between 2006-2013. We use 2440 Hinode photospheric vector magnetograms, and the corresponding subsurface fluid velocity data derived from GONG (2006-2012) and HMI (2010-2013) Dopplergrams. We find a significant hemispheric bias in all three parameters. The subsurface kinetic helicity is preferentially positive in the southern hemisphere and negative in the northern hemisphere. The photospheric current helicity and the $\alpha$ parameter have the same bias for strong fields ($|B|>1000$ G) and no significant bias for weak fields (100 G $<|B|<500$ G). We find no significant region-by-region correlation between the subsurface kinetic helicity and either the strong-field current helicity or $\alpha$. Subsurface fluid motions of a given handedness correspond to photospheric helicities of both signs in approximately equal numbers. However, common variations appear in annual averages of these quantities over all regions. Furthermore, in a subset of 77 regions we find significant correlations between the temporal profiles of the subsurface and photospheric helicities. In these cases, the sign of the linear correlation coefficient matches the sign relationship between the helicities, indicating that the photospheric magnetic field twist is sensitive to the twisting motions below the surface.