Resolving the phase structure of neutral hydrogen (HI) is crucial for understanding the life cycle of the interstellar medium (ISM). However, accurate measurements of HI temperature and density are limited by the availability of background continuum sources for measuring HI absorption. Here we test the use of deep learning for extracting HI properties over large areas without optical depth information. We train a 1D convolutional neural network using synthetic observations of 3D numerical simulations of the ISM to predict the fraction of cold neutral medium (f_CNM) and the correction to the optically-thin HI column density for optical depth (R_HI) from $21\rm\,cm$ emission alone. We restrict our analysis to high Galactic latitudes ($|b|>30 deg), where the complexity of spectral line profiles is minimized. We verify that the network accurately predicts f_CNM and R_HI by comparing the results with direct constraints from 21 cm absorption. By applying the network to the GALFA-HI survey, we generate large-area maps of f_CNM and R_HI. Although the overall contribution to the total HI column of cold neutral medium (CNM)-rich structures is small (~5%), we find that these structures are ubiquitous. Our results are consistent with the picture that small-scale structures observed in 21 cm emission aligned with the magnetic field are dominated by CNM. Finally, we demonstrate that the observed correlation between HI column density and dust reddening (E(B-V)) declines with increasing R_HI, indicating that future efforts to quantify foreground Galactic E(B-V) using HI, even at high latitudes, should increase fidelity by accounting for HI phase structure.

We report on the discovery of a transiting Earth-sized (0.95$R_\oplus$) planet around an M3.5 dwarf star at 57\,pc, EPIC~249631677. The planet has a period of $\sim$3.14 days, i.e. ${\sim}\pi$, with an instellation of 7.5\,S$_{\oplus}$. The detection was made using publicly available data from \textit{K2}'s Campaign 15. We observed three additional transits with SPECULOOS Southern and Northern Observatories, and a stellar spectrum from Keck/HIRES, which allowed us to validate the planetary nature of the signal. The confirmed planet is well suited for comparative terrestrial exoplanetology. While exoplanets transiting ultracool dwarfs present the best opportunity for atmospheric studies of terrestrial exoplanets with the \textit{James Webb Space Telescope}, those orbiting mid-M dwarfs within 100\,pc such as EPIC~249631677b will become increasingly accessible with the next generation of observatories (e.g., \textit{HabEx, LUVOIR, OST}).

Tidal effects in planetary systems are the main driver in the orbital migration of natural satellites. They result from physical processes occurring deep inside celestial bodies, whose effects are rarely observable from surface imaging. For giant planet systems, the tidal migration rate is determined by poorly understood dissipative processes in the planet, and standard theories suggest an orbital expansion rate inversely proportional to the power 11/2 in distance, implying little migration for outer moons such as Saturn's largest moon, Titan. Here, we use two independent measurements obtained with the Cassini spacecraft to measure Titan's orbital expansion rate. We find Titan migrates away from Saturn at 11.3 $\pm$ 2.0 cm/year, corresponding to a tidal quality factor of Saturn of Q $\simeq$ 100, and a migration timescale of roughly 10 Gyr. This rapid orbital expansion suggests Titan formed significantly closer to Saturn and has migrated outward to its current position. Our results for Titan and five other moons agree with the predictions of a resonance locking tidal theory, sustained by excitation of inertial waves inside the planet. The associated tidal expansion is only weakly sensitive to orbital distance, motivating a revision of the evolutionary history of Saturn's moon system. The resonance locking mechanism could operate in other systems such as stellar binaries and exoplanet systems, and it may allow for tidal dissipation to occur at larger orbital separations than previously believed.

We report the discovery of two transiting exoplanets from the WASP survey, WASP-150b and WASP-176b. WASP-150b is an eccentric ($e$ = 0.38) hot Jupiter on a 5.6 day orbit around a $V$ = 12.03, F8 main-sequence host. The host star has a mass and radius of 1.4 $\rm M_{\odot}$ and 1.7 $\rm R_{\odot}$ respectively. WASP-150b has a mass and radius of 8.5 $\rm M_J$ and 1.1 $\rm R_J$, leading to a large planetary bulk density of 6.4 $\rm \rho_J$. WASP-150b is found to be $\sim3$ Gyr old, well below its circularisation timescale, supporting the eccentric nature of the planet. WASP-176b is a hot Jupiter planet on a 3.9 day orbit around a $V$ = 12.01, F9 sub-giant host. The host star has a mass and radius of 1.3 $\rm M_{\odot}$ and 1.9 $\rm R_{\odot}$. WASP-176b has a mass and radius of 0.86 $\rm M_J$ and 1.5 $\rm R_J$ respectively, leading to a planetary bulk density of 0.23 $\rm \rho_J$.

SPECULOOS-South, an observatory composed of four independent 1m robotic telescopes, located at ESO Paranal, Chile, started scientific operation in January 2019. This Southern Hemisphere facility operates as part of SPECULOOS, an international network of 1m-class telescopes surveying for transiting terrestrial planets around the nearest and brightest ultra-cool dwarfs. To automatically and efficiently process the observations of SPECULOOS-South, and to deal with the specialised photometric requirements of ultra-cool dwarf targets, we present our automatic pipeline. This pipeline includes an algorithm for automated differential photometry and an extensive correction technique for the effects of telluric water vapour, using ground measurements of the precipitable water vapour. Observing very red targets in the near-infrared can result in photometric systematics in the differential lightcurves, related to the temporally-varying, wavelength-dependent opacity of the Earth's atmosphere. These systematics are sufficient to affect the daily quality of the lightcurves, the longer time-scale variability study of our targets and even mimic transit-like signals. Here we present the implementation and impact of our water vapour correction method. Using the 179 nights and 98 targets observed in the I+z' filter by SPECULOOS-South since January 2019, we show the impressive photometric performance of the facility (with a median precision of ~1.5 mmag for 30-min binning of the raw, non-detrended lightcurves) and assess its detection potential. We compare simultaneous observations with SPECULOOS-South and TESS, to show that we readily achieve high-precision, space-level photometry for bright, ultra-cool dwarfs, highlighting SPECULOOS-South as the first facility of its kind.

We have studied the small scale distribution of atomic hydrogen (HI) using 21-cm absorption spectra against multiple-component background radio continuum sources from the 21-SPONGE survey and the Millennium Arecibo Absorption Line Survey. We have found $>5\sigma$ optical depth variations at a level of $\sim0.03-0.5$ between 13 out of 14 adjacent sightlines separated by a few arcseconds to a few arcminutes, suggesting the presence of neutral structures on spatial scales from a few to thousands of AU (which we refer to as tiny scale atomic structure, TSAS). The optical depth variations are strongest in directions where the HI column density and the fraction of HI in the cold neutral medium (CNM) are highest, which tend to be at low Galactic latitudes. By measuring changes in the properties of Gaussian components fitted to the absorption spectra, we find that changes in both the peak optical depth and the linewidth of TSAS absorption features contribute to the observed optical depth variations, while changes in the central velocity do not appear to strongly impact the observed variations. Both thermal and turbulent motions contribute appreciably to the linewidths, but the turbulence does not appear strong enough to confine overpressured TSAS. In a majority of cases, the TSAS column densities are sufficiently high that these structures can radiatively cool fast enough to maintain thermal equilibrium with their surroundings, even if they are overpressured. We also find that a majority of TSAS is associated with the CNM. For TSAS in the direction of the Taurus molecular cloud and the local Leo cold cloud, we estimate densities over an order of magnitude higher than typical CNM densities.

Mass, radius, and age are three of the most fundamental parameters for celestial objects, enabling studies of the evolution and internal physics of stars, brown dwarfs, and planets. Brown dwarfs are hydrogen-rich objects that are unable to sustain core fusion reactions but are supported from collapse by electron degeneracy pressure. As they age, brown dwarfs cool, reducing their radius and luminosity. Young exoplanets follow a similar behaviour. Brown dwarf evolutionary models are relied upon to infer the masses, radii and ages of these objects. Similar models are used to infer the mass and radius of directly imaged exoplanets. Unfortunately, only sparse empirical mass, radius and age measurements are currently available, and the models remain mostly unvalidated. Double-line eclipsing binaries provide the most direct route for the absolute determination of the masses and radii of stars. Here, we report the SPECULOOS discovery of 2M1510A, a nearby, eclipsing, double-line brown dwarf binary, with a widely-separated tertiary brown dwarf companion. We also find that the system is a member of the $45\pm5$ Myr-old moving group, Argus. The system's age matches those of currently known directly-imaged exoplanets. 2M1510A provides an opportunity to benchmark evolutionary models of brown dwarfs and young planets. We find that widely-used evolutionary models do reproduce the mass, radius and age of the binary components remarkably well, but overestimate the luminosity by up to 0.65 magnitudes, which could result in underestimated photometric masses for directly-imaged exoplanets and young field brown dwarfs by 20 to 35%.

We investigate the kinematics of neutral gas in the Small Magellanic Cloud (SMC) and test the hypothesis that it is rotating in a disk. To trace the 3D motions of the neutral gas distribution, we identify a sample of young, massive stars embedded within it. These are stars with radial velocity measurements from spectroscopic surveys and proper motion measurements from Gaia, whose radial velocities match with dominant HI components. We compare the observed radial and tangential velocities of these stars with predictions from the state-of-the-art rotating disk model based on high-resolution 21 cm observations of the SMC from the Australian Square Kilometer Array Pathfinder telescope. We find that the observed kinematics of gas-tracing stars are inconsistent with disk rotation. We conclude that the kinematics of gas in the SMC are more complex than can be inferred from the integrated radial velocity field. As a result of violent tidal interactions with the LMC, non-rotational motions are prevalent throughout the SMC, and it is likely composed of distinct sub-structures overlapping along the line of sight.

We report the discovery and characterisation of WASP-180Ab, a hot Jupiter confirmed by the detection of its Doppler shadow and by measuring its mass using radial velocities. We find the 0.9 $\pm$ 0.1 $M_{\rm Jup}$, 1.24 $\pm$ 0.04 $R_{\rm Jup}$ planet to be in a misaligned, retrograde orbit around an F7 star with $T_{\rm eff}$ = 6500K and a moderate rotation speed of vsini = 19.9 km s$^{-1}$. The host star is the primary of a $V$ = 10.7 binary, where a secondary separated by 5$''$ ($\sim$1200 AU) contributes $\sim$30% of the light. WASP-180Ab therefore adds to a small sample of transiting hot Jupiters known in binary systems. A 4.6-day modulation seen in the WASP data is likely to be the rotational modulation of the companion star, WASP-180B.

The Probe of Inflation and Cosmic Origins (PICO) is a proposed probe-scale space mission consisting of an imaging polarimeter operating in frequency bands between 20 and 800 GHz. We describe the science achievable by PICO, which has sensitivity equivalent to more than 3300 Planck missions, the technical implementation, the schedule and cost.

Our understanding of the dynamics of the interstellar medium is informed by the study of the detailed velocity structure of emission line observations. One approach to study the velocity structure is to decompose the spectra into individual velocity components; this leads to a description of the dataset that is significantly reduced in complexity. However, this decomposition requires full automation lest it becomes prohibitive for large datasets, such as Galactic plane surveys. We developed GaussPy+, a fully automated Gaussian decomposition package that can be applied to emission line datasets, especially large surveys of HI and isotopologues of CO. We built our package upon the existing GaussPy algorithm and significantly improved its performance for noisy data. New functionalities of GaussPy+ include: i) automated preparatory steps, such as an accurate noise estimation, which can also be used as standalone applications; ii) an improved fitting routine; iii) an automated spatial refitting routine that can add spatial coherence to the decomposition results by refitting spectra based on neighbouring fit solutions. We thoroughly tested the performance of GaussPy+ on synthetic spectra and a test field from the Galactic Ring Survey. We found that GaussPy+ can deal with cases of complex emission and even low to moderate signal-to-noise values.

We report Arecibo 21 cm absorption-emission observations to characterise the physical properties of neutral hydrogen (HI) in the proximity of five giant molecular clouds (GMCs): Taurus, California, Rosette, Mon OB1, and NGC 2264. Strong HI absorption was detected toward all 79 background continuum sources in the ~60x20 square degree region. Gaussian decompositions were performed to estimate temperatures, optical depths and column densities of the cold and warm neutral medium (CNM, WNM). The properties of individual CNM components are similar to those previously observed along random Galactic sightlines and in the vicinity of GMCs, suggesting a universality of cold HI properties. The CNM spin temperature (Ts) histogram peaks at ~50K. The turbulent Mach numbers of CNM vary widely, with a typical value of ~4, indicating that their motions are supersonic. About 60% of the total HI gas is WNM, and nearly 40% of the WNM lies in thermally unstable regime 500-5000K. The observed CNM fraction is higher around GMCs than in diffuse regions, and increases with increasing column density (NHI) to a maximum of ~75%. On average, the optically thin approximation (N*(HI)) underestimates the total N(HI) by ~21%, but we find large regional differences in the relationship between N(HI) and the required correction factor, f=N(HI)/N*(HI). We examine two different methods (linear fit of f vs log10(N*(HI)) and uniform Ts) to correct for opacity effects using emission data from the GALFA-HI survey. We prefer the uniform Ts method, since the linear relationship does not produce convincing fits for all subregions.

We report the Transiting Exoplanet Survey Satellite (TESS) discovery of three terrestrial-sized planets transiting L 98-59 (TOI-175, TIC 307210830) -- a bright M dwarf at a distance of 10.6 pc. Using the Gaia-measured distance and broad-band photometry we find that the host star is an M3 dwarf. Combined with the TESS transits from three sectors, the corresponding stellar parameters yield planet radii ranging from 0.8REarth to 1.6REarth. All three planets have short orbital periods, ranging from 2.25 to 7.45 days with the outer pair just wide of a 2:1 period resonance. Diagnostic tests produced by the TESS Data Validation Report and the vetting package DAVE rule out common false positive sources. These analyses, along with dedicated follow-up and the multiplicity of the system, lend confidence that the observed signals are caused by planets transiting L 98-59 and are not associated with other sources in the field. The L 98-59 system is interesting for a number of reasons: the host star is bright (V = 11.7 mag, K = 7.1 mag) and the planets are prime targets for further follow-up observations including precision radial-velocity mass measurements and future transit spectroscopy with the James Webb Space Telescope; the near resonant configuration makes the system a laboratory to study planetary system dynamical evolution; and three planets of relatively similar size in the same system present an opportunity to study terrestrial planets where other variables (age, metallicity, etc.) can be held constant. L 98-59 will be observed in 4 more TESS sectors, which will provide a wealth of information on the three currently known planets and have the potential to reveal additional planets in the system.

We present an analytical model to study the dynamics of the outer edge of Saturn's A ring. The latter is influenced by 7:6 mean motion resonances with Janus and Epimetheus. Because of the horseshoe motion of the two co-orbital moons, the ring edge particles are alternately trapped in a corotation eccentricity resonance (CER) or a Lindblad eccentricity resonance (LER). However, the resonance oscillation periods are longer than the 4-year interval between the switches in the orbits of Janus and Epimetheus. Averaged equations of motion are used, and our model is numerically integrated to describe the effects of the periodic sweeping of the 7:6 CERs and LERs over the ring edge region. We show that four radial zones (ranges 136715-136723, 136738-136749, 136756-136768, 136783-136791 km) are chaotic on decadal timescales, within which particle semi-major axes have periodic changes due to partial libration motions around the CER fixed points. After a few decades, the maximum variation of semi-major axis is about 11 km (respectively 3 km) in the case of the CER with Janus (respectively Epimetheus). Similarly, particle eccentricities have partial oscillations forced by the LERs every 4 yr. For initially circular orbits, the maximum eccentricity reached is ~0.001. We apply our work to "Peggy", an object recently discovered at the ring edge, confirming that it is strongly perturbed by the Janus 7:6 LER. The CER has currently no effect on that body, nevertheless the fitted semi-major axes are just outside the chaotic zone of radial range 136756-136768 km.

During the thirteen years in orbit around Saturn before its final plunge, the Cassini spacecraft provided more than ten thousand astrometric measurements. Such large amounts of accurate data enable the search for extremely faint signals in the orbital motion of the saturnian moons. Among these, the detection of the dynamical feedback of the rotation of the inner moons of Saturn on their respective orbits becomes possible. Using all the currently available astrometric data associated with Atlas, Prometheus, Pandora, Janus and Epimetheus, we first provide a detailed analysis of the Cassini Imaging Science Subsystem (ISS) data, with special emphasis on their statistical behavior and sources of bias. Then, we give updated estimates of the moons' averaged densities and try to infer more details about their interior properties by estimating the physical librations for Prometheus, Pandora, Epimetheus and Janus from anomalies in their apsidal precession. Our results are compatible with a homogeneous interior for Janus and Epimetheus, within the uncertainty of the measurements. On the other hand, we found some inconsistency for Pandora and Prometheus, which might result from a dynamical mismodeling of Saturn's gravity field. Last, we show how the synergistic introduction of libration measurements directly derived from imaging should allow the moons' moments of inertia to be better constrained.

The Probe of Inflation and Cosmic Origins (PICO) is an imaging polarimeter that will scan the sky for 5 years in 21 frequency bands spread between 21 and 799 GHz. It will produce full-sky surveys of intensity and polarization with a final combined-map noise level of 0.87 $\mu$K arcmin for the required specifications, equivalent to 3300 Planck missions, and with our current best-estimate would have a noise level of 0.61 $\mu$K arcmin (6400 Planck missions). PICO will either determine the energy scale of inflation by detecting the tensor to scalar ratio at a level $r=5\times 10^{-4}~(5\sigma)$, or will rule out with more than $5\sigma$ all inflation models for which the characteristic scale in the potential is the Planck scale. With LSST's data it could rule out all models of slow-roll inflation. PICO will detect the sum of neutrino masses at $>4\sigma$, constrain the effective number of light particle species with $\Delta N_{\rm eff}<0.06~(2\sigma)$, and elucidate processes affecting the evolution of cosmic structures by measuring the optical depth to reionization with errors limited by cosmic variance and by constraining the evolution of the amplitude of linear fluctuations $\sigma_{8}(z)$ with sub-percent accuracy. Cross-correlating PICO's map of the thermal Sunyaev-Zeldovich effect with LSST's gold sample of galaxies will precisely trace the evolution of thermal pressure with $z$. PICO's maps of the Milky Way will be used to determine the make up of galactic dust and the role of magnetic fields in star formation efficiency. With 21 full sky legacy maps in intensity and polarization, which cannot be obtained in any other way, the mission will enrich many areas of astrophysics. PICO is the only single-platform instrument with the combination of sensitivity, angular resolution, frequency bands, and control of systematic effects that can deliver this compelling, timely, and broad science.

We report the discovery of WASP-190b, an exoplanet on a 5.37-day orbit around a mildly-evolved F6 IV-V star with V = 11.7, T_eff = 6400 $\pm$ 100 K, M$_{*}$ = 1.35 $\pm$ 0.05 M_sun and R$_{*}$ = 1.6 $\pm$ 0.1 R_sun. The planet has a radius of R_p = 1.15 $\pm$ 0.09 R_Jup and a mass of M_p = 1.0 $\pm$ 0.1 M_Jup, making it a mildly inflated hot Jupiter. It is the first hot Jupiter confirmed via Doppler tomography with an orbital period >5 days. The orbit is also marginally misaligned with respect to the stellar rotation, with $\lambda$ = 21 $\pm$ 6$^{\circ}$ measured using Doppler tomography.

We report the Transiting Exoplanet Survey Satellite ($TESS$) detection of a multi-planet system orbiting the $V=10.9$ K0 dwarf TOI 125. We find evidence for up to five planets, with varying confidence. Three high signal-to-noise transit signals correspond to sub-Neptune-sized planets ($2.76$, $2.79$, and $2.94\ R_{\oplus}$), and we statistically validate the planetary nature of the two inner planets ($P_b = 4.65$ days, $P_c = 9.15$ days). With only two transits observed, we report the outer object ($P_{.03} = 19.98$ days) as a high signal-to-noise ratio planet candidate. We also detect a candidate transiting super-Earth ($1.4\ R_{\oplus}$) with an orbital period of only $12.7$ hours and a candidate Neptune-sized planet ($4.2\ R_{\oplus}$) with a period of $13.28$ days, both at low signal-to-noise. This system is amenable to mass determination via radial velocities and transit timing variations, and provides an opportunity to study planets of similar size while controlling for age and environment. The ratio of orbital periods between TOI 125 b and c ($P_c/P_b = 1.97$) is slightly smaller than an exact 2:1 commensurability and is atypical of multiple planet systems from $Kepler$, which show a preference for period ratios just $wide$ of first-order period ratios. A dynamical analysis refines the allowed parameter space through stability arguments and suggests that, despite the nearly commensurate periods, the system is unlikely to be in resonance.

We present the discovery by the WASP-South transit survey of three new transiting hot Jupiters, WASP-161 b, WASP-163 b and WASP-170 b. Follow-up radial velocities obtained with the Euler/CORALIE spectrograph and high-precision transit light curves obtained with the TRAPPIST-North, TRAPPIST-South, SPECULOOS-South, NITES, and Euler telescopes have enabled us to determine the masses and radii for these transiting exoplanets. WASP-161\,b completes an orbit around its $V=11.1$ F6V-type host star in 5.406 days, and has a mass and radius of $2.5\pm 0.2$$M_{Jup}$ and $1.14\pm 0.06$ $R_{Jup}$ respectively. WASP-163\,b has an orbital period of 1.609 days, a mass of $1.9\pm0.2$ $M_{Jup}$, and a radius of $1.2\pm0.1$ $R_{Jup}$. Its host star is a $V=12.5$ G8-type dwarf. WASP-170\,b is on a 2.344 days orbit around a G1V-type star of magnitude $V=12.8$. It has a mass of $1.7\pm0.2$ $M_{Jup}$ and a radius of $1.14\pm0.09$ $R_{Jup}$. Given their irradiations ($\sim10^9$ erg.s$^{-1}$.cm$^{-2}$) and masses, the three new planets sizes are in good agreement with classical structure models of irradiated giant planets.

The TRAPPIST-1 planetary system represents an exceptional opportunity for the atmospheric characterization of temperate terrestrial exoplanets with the upcoming James Webb Space Telescope (JWST). Assessing the potential impact of stellar contamination on the planets' transit transmission spectra is an essential precursor step to this characterization. Planetary transits themselves can be used to scan the stellar photosphere and to constrain its heterogeneity through transit depth variations in time and wavelength. In this context, we present our analysis of 169 transits observed in the optical from space with K2 and from the ground with the SPECULOOS and Liverpool telescopes. Combining our measured transit depths with literature results gathered in the mid/near-IR with Spitzer/IRAC and HST/WFC3, we construct the broadband transmission spectra of the TRAPPIST-1 planets over the 0.8-4.5 $\mu$m spectral range. While planets b, d, and f spectra show some structures at the 200-300ppm level, the four others are globally flat. Even if we cannot discard their instrumental origins, two scenarios seem to be favored by the data: a stellar photosphere dominated by a few high-latitude giant (cold) spots, or, alternatively, by a few small and hot (3500-4000K) faculae. In both cases, the stellar contamination of the transit transmission spectra is expected to be less dramatic than predicted in recent papers. Nevertheless, based on our results, stellar contamination can still be of comparable or greater order than planetary atmospheric signals at certain wavelengths. Understanding and correcting the effects of stellar heterogeneity therefore appears essential to prepare the exploration of TRAPPIST-1's with JWST.

We report the discovery of NGTS-4b, a sub-Neptune-sized planet transiting a 13th magnitude K-dwarf in a 1.34d orbit. NGTS-4b has a mass M=$20.6\pm3.0$M_E and radius R=$3.18\pm0.26$R_E, which places it well within the so-called "Neptunian Desert". The mean density of the planet ($3.45\pm0.95$g/cm^3) is consistent with a composition of 100% H$_2$O or a rocky core with a volatile envelope. NGTS-4b is likely to suffer significant mass loss due to relatively strong EUV/X-ray irradiation. Its survival in the Neptunian desert may be due to an unusually high core mass, or it may have avoided the most intense X-ray irradiation by migrating after the initial activity of its host star had subsided. With a transit depth of $0.13\pm0.02$%, NGTS-4b represents the shallowest transiting system ever discovered from the ground, and is the smallest planet discovered in a wide-field ground-based photometric survey.

We present here SPECULOOS, a new exoplanet transit search based on a network of 1m-class robotic telescopes targeting the $\sim$1200 ultracool (spectral type M7 and later) dwarfs bright enough in the infrared ($K$-mag $\leq 12.5$) to possibly enable the atmospheric characterization of temperate terrestrial planets with next-generation facilities like the $\textit{James Webb Space Telescope}$. The ultimate goals of the project are to reveal the frequency of temperate terrestrial planets around the lowest-mass stars and brown dwarfs, to probe the diversity of their bulk compositions, atmospheres and surface conditions, and to assess their potential habitability.

We present 21-cm Spectral Line Observations of Neutral Gas with the VLA (21-SPONGE), a Karl G. Jansky Very Large Array (VLA) large project (~600 hours) for measuring the physical properties of Galactic neutral hydrogen (HI). 21-SPONGE is distinguished among previous Galactic HI studies as a result of: (1) exceptional optical depth sensitivity ($\sigma_{\tau} < 10^{-3}$ per $0.42\rm\,km\,s^{-1}$ channels over 57 lines of sight); (2) matching 21 cm emission spectra with highest-possible angular resolution (~4') from the Arecibo Observatory; (3) detailed comparisons with numerical simulations for assessing observational biases. We autonomously decompose 21 cm spectra and derive the physical properties (i.e., spin temperature, $T_s$, column density) of the cold neutral medium (CNM; $T_s<250\rm\,K$), thermally unstable medium (UNM; $250< T_s < 1000\rm\,K$) and warm neutral medium (WNM; $T_s > 1000\rm\,K$) simultaneously. We detect 50% of the total HI mass in absorption, the majority of which is CNM (56 +/- 10%, corresponding to 28% of the total HI mass). Although CNM is detected ubiquitously, the CNM fraction along most lines of sight is <50%. We find that 20% of the total HI mass is thermally unstable (41 +/- 10% of HI detected in absorption), with no significant variation with Galactic environment. Finally, although the WNM comprises 52% of the total HI mass, we detect little evidence for WNM absorption with $1000<T_s<4000\rm\,K$. Following spectral modeling, we detect a stacked residual absorption feature corresponding to WNM with $T_s\sim10^4\rm\,K$. We conclude that excitation in excess of collisions likely produces significantly higher WNM $T_s$ than predicted by steady-state models.

The local interstellar medium (ISM) is suffused with "dark" gas, identified by excess infrared and gamma ray emission, yet undetected by standard ISM tracers such as neutral hydrogen (HI) or carbon monoxide emission. Based on observed dust properties from Planck, recent studies have argued that HI mixed with dust is strongly saturated and that dark gas is dominated by optically-thick HI. We test this hypothesis by reproducing this model using data from Planck and new 21 cm emission maps from GALFA-HI -- the first large-area 21cm emission survey with comparable angular resolution to Planck. We compare the results with those from a large sample of HI column densities based on direct observations of HI optical depth, and find that the inferred column density corrections are significantly lower than those inferred by the Planck-based model. Further, we rule out the hypothesis that the pencil-beam HI absorption sight lines preferentially miss opaque "blobs" with small covering fraction, as these structures require densities and pressures which are incompatible with ISM conditions. Our results support the picture that excess dust emission in the local ISM is not dominated by optically-thick HI, but is rather a combination of intrinsic changes in dust grain emissivities and H2 missed by CO observations.

Observations of interstellar dust are often used as a proxy for total gas column density $N_\mathrm{H}$. By comparing $\textit{Planck}$ thermal dust data (Release 1.2) and new dust reddening maps from Pan-STARRS 1 and 2MASS (Green et al. 2018), with accurate (opacity-corrected) HI column densities and newly-published OH data from the Arecibo Millennium survey and 21-SPONGE, we confirm linear correlations between dust optical depth $\tau_{353}$, reddening $E(B{-}V)$ and the total proton column density $N_\mathrm{H}$ in the range (1$-$30)$\times$10$^{20}$cm$^{-2}$, along sightlines with no molecular gas detections in emission. We derive an $N_\mathrm{H}$/$E(B{-}V)$ ratio of (9.4$\pm$1.6)$\times$10$^{21}$cm$^{-2}$mag$^{-1}$ for purely atomic sightlines at $|b|$$>$5$^{\circ}$, which is 60$\%$ higher than the canonical value of Bohlin et al. (1978). We report a $\sim$40$\%$ increase in opacity $\sigma_{353}$=$\tau_{353}$/$N_\mathrm{H}$, when moving from the low column density ($N_\mathrm{H}$$<$5$\times$10$^{20}$cm$^{-2}$) to moderate column density ($N_\mathrm{H}$$>$5$\times$10$^{20}$cm$^{-2}$) regime, and suggest that this rise is due to the evolution of dust grains in the atomic ISM. Failure to account for HI opacity can cause an additional apparent rise in $\sigma_{353}$, of the order of a further $\sim$20$\%$. We estimate molecular hydrogen column densities $N_{\mathrm{H}_{2}}$ from our derived linear relations, and hence derive the OH/H$_2$ abundance ratio of $X_\mathrm{OH}$$\sim$1$\times$10$^{-7}$ for all molecular sightlines. Our results show no evidence of systematic trends in OH abundance with $N_{\mathrm{H}_{2}}$ in the range $N_{\mathrm{H}_{2}}$$\sim$(0.1$-$10)$\times$10$^{21}$cm$^{-2}$. This suggests that OH may be used as a reliable proxy for H$_2$ in this range, which includes sightlines with both CO-dark and CO-bright gas.

HI self absorption (HISA) clouds are clumps of cold neutral hydrogen (HI) visible in front of warm background gas, which makes them ideal places to study the properties of the cold atomic component of the interstellar medium (ISM). The Riegel-Crutcher (R-C) cloud is the most striking HISA feature in the Galaxy. It is one of the closest HISA clouds to us and is located in the direction of the Galactic Centre, which provides a bright background. High-resolution interferometric measurements have revealed the filamentary structure of this cloud, however it is difficult to accurately determine the temperature and the density of the gas without optical depth measurements. In this paper we present new HI absorption observations with the Australia Telescope Compact Array (ATCA) against 46 continuum sources behind the Riegel-Crutcher cloud to directly measure the optical depth of the cloud. We decompose the complex HI absorption spectra into Gaussian components using an automated machine learning algorithm. We find 300 Gaussian components, from which 67 are associated with the R-C cloud (0 < v_{LSR} < 10 kms, FWHM < 10 kms). Combining the new HI absorption data with HI emission data from previous surveys we calculate the spin temperature and find it to be between 20 and 80 K. Our measurements uncover a temperature gradient across the cloud with spin temperatures decreasing towards positive Galactic latitudes. We also find three new OH absorption lines associated with the cloud, which support the presence of molecular gas.

We present the discovery of NGTS-3Ab, a hot Jupiter found transiting the primary star of an unresolved binary system. We develop a joint analysis of multi-colour photometry, centroids, radial velocity (RV) cross-correlation function (CCF) profiles and their bisector inverse slopes (BIS) to disentangle this three-body system. Data from the Next Generation Transit Survey (NGTS), SPECULOOS and HARPS are analysed and modelled with our new blendfitter software. We find that the binary consists of NGTS-3A (G6V-dwarf) and NGTS-3B (K1V-dwarf) at <1 arcsec separation. NGTS-3Ab orbits every P = 1.675 days. The planet radius and mass are R_planet = 1.48+-0.37 R_J and M_planet = 2.38+-0.26 M_J, suggesting it is potentially inflated. We emphasise that only combining all the information from multi-colour photometry, centroids and RV CCF profiles can resolve systems like NGTS-3. Such systems cannot be disentangled from single-colour photometry and RV measurements alone. Importantly, the presence of a BIS correlation indicates a blend scenario, but is not sufficient to determine which star is orbited by the third body. Moreover, even if no BIS correlation is detected, a blend scenario cannot be ruled out without further information. The choice of methodology for calculating the BIS can influence the measured significance of its correlation. The presented findings are crucial to consider for wide-field transit surveys, which require wide CCD pixels (>5 arcsec) and are prone to contamination by blended objects. With TESS on the horizon, it is pivotal for the candidate vetting to incorporate all available follow-up information from multi-colour photometry and RV CCF profiles.

Hydroxyl (OH) is expected to be abundant in diffuse interstellar molecular gas as it forms along with $H_2$ under similar conditions and within a similar extinction range. We have analyzed absorption measurements of OH at 1665 MHz and 1667 MHz toward 44 extragalactic continuum sources, together with the J=1-0 transitions of $^{12}$CO, $^{13}$CO , and C$^{18}$O, and the J=2-1 of $^{12}$CO. The excitation temperature of OH were found to follow a modified log-normal distribution, $ f(T{\rm_{ex}}) \propto \frac{1}{ \sqrt{2\pi}\sigma } \rm{exp}\left[-\frac{[ln(\textit{T}_{ex})-ln(3.4\ K)]^2}{2\sigma^2}\right] $, the peak of which is close to the temperature of the Galactic emission background (CMB+synchron). In fact, 90% of the OH has excitation temperature within 2 K of the Galactic background at the same location, providing a plausible explanation for the apparent difficulty to map this abundant molecule in emission. The opacities of OH were found to be small and peak around 0.01. For gas at intermediate extinctions (A$\rm_V$ $\sim$ 0.05-2 mag), the detection rate of OH with detection limit $N(\mathrm{OH})\simeq 10^{12}$ cm$^{-2}$ is approximately independent of $A\rm_V$. We conclude that OH is abundant in the diffuse molecular gas and OH absorption is a good tracer of `dark molecular gas (DMG)'. The measured fraction of DMG depends on assumed detection threshold of the CO data set. The next generation of highly sensitive low frequency radio telescopes, FAST and SKA, will make feasible the systematic inventory of diffuse molecular gas, through decomposing in velocity the molecular (e.g. OH and CH) absorption profiles toward background continuum sources with numbers exceeding what is currently available by orders of magnitude.

We analyze synthetic neutral hydrogen (HI) absorption and emission spectral lines from a high- resolution, three-dimensional hydrodynamical simulation to quantify how well observational methods recover the physical properties of interstellar gas. We present a new method for uniformly decomposing HI spectral lines and estimating the properties of associated gas using the Autonomous Gaussian Decomposition (AGD) algorithm. We find that Hi spectral lines recover physical structures in the simulation with excellent completeness at high Galactic latitude, and this completeness declines with decreasing latitude due to strong velocity-blending of spectral lines. The temperature and column density inferred from our decomposition and radiative transfer method agree with the simulated values within a factor of < 2 for the majority of gas structures. We next compare synthetic spectra with observations from the 21-SPONGE survey at the Karl G. Jansky Very Large Array using AGD. We find more components per line of sight in 21-SPONGE than in synthetic spectra, which reflects insufficient simulated gas scale heights and the limitations of local box simulations. In addition, we find a significant population of low-optical depth, broad absorption components in the synthetic data which are not seen in 21-SPONGE. This population is not obvious in integrated or per-channel diagnostics, and reflects the benefit of studying velocity-resolved components. The discrepant components correspond to the highest spin temperatures (1000 < Ts < 4000 K), which are not seen in 21-SPONGE despite sufficient observational sensitivity. We demonstrate that our analysis method is a powerful tool for diagnosing neutral ISM conditions, and future work is needed to improve observational statistics and implementation of simulated physics.

There is evidence to suggest that clathrate hydrates have a significant effect on the surface geology of icy bodies in the Solar System. However the aqueous environments believed to be present on these bodies are likely to be saline rather than pure water. Laboratory work to underpin the properties of clathrates in such environments is lacking. We fill this gap by carrying out a laboratory investigation of the physical properties of CO$_2$ clathrates produced in weak aqueous solutions of MgSO$_4$. We use synchrotron X-ray powder diffraction to investigate clathrates formed at high CO$_2$ pressure in ice that has formed from aqueous solutions of MgSO$_4$. We measure the thermal expansion, density and dissociation properties of the clathrates under temperature conditions similar to those on icy Solar System bodies. We find that the sulphate solution inhibits the formation of clathrates by lowering their dissociation temperatures. Hysteresis is found in the thermal expansion coefficients as clathrates are cooled and heated; we attribute this to the presence of the salt in solution. The density derived from X-ray powder diffraction is temperature and pressure dependent. When comparing the density of CO$_2$ clathrates to that of the solution in which they formed, we conclude that they sink in the oceans in which they form. We also find that the polymorph of ice present at low temperatures is Ih rather than Ic, which we attribute to the presence of the MgSO$_4$. We 1) conclude that the clathrate density has implications for their behaviour in satellite oceans as their sinking and floating capabilities are temperature and pressure dependent, 2) conclude that the presence of MgSO$_4$ inhibits the formation of clathrates and in some cases may even affect their structure and 3) report the dominance of Ih throughout the experimental procedure despite Ic being the stable phase at low temperature.

We revisit the dynamics of Atlas. Using Cassini ISS astrometric observations spanning February 2004 to August 2013, Cooper et al. (2015) found evidence that Atlas is currently perturbed by both a 54:53 corotation eccentricity resonance (CER) and a 54:53 Lindblad eccentricity resonance (LER) with Prometheus. They demonstrated that the orbit of Atlas is chaotic, with a Lyapunov time of order 10 years, as a direct consequence of the coupled resonant interaction (CER/LER) with Prometheus. Here we investigate the interactions between the two resonances using the CoraLin analytical model (El Moutamid et al. 2014), showing that the chaotic zone fills almost all the corotation sites occupied by the satellite's orbit. Four 70:67 apse-type mean motion resonances with Pandora are also overlapping, but these resonances have a much weaker effect. Frequency analysis allows us to highlight the coupling between the 54:53 resonances, and confirms that a simplified system including the perturbations due to Prometheus and Saturn's oblateness only captures the essential features of the dynamics.

Using astrometric observations spanning more than a century and including a large set of Cassini data, we determine Saturn's tidal parameters through their current effects on the orbits of the eight main and four coorbital moons. We have used the latter to make the first determination of Saturn's Love number, $k_2=0.390 \pm 0.024$, a value larger than the commonly used theoretical value of 0.341 (Gavrilov & Zharkov, 1977), but compatible with more recent models (Helled & Guillot, 2013) for which $k_2$ ranges from 0.355 to 0.382. Depending on the assumed spin for Saturn's interior, the new constraint can lead to a reduction of up to 80% in the number of potential models, offering great opportunities to probe the planet's interior. In addition, significant tidal dissipation within Saturn is confirmed (Lainey et al., 2012) corresponding to a high present-day tidal ratio $k_2/Q=(1.59 \pm 0.74) \times 10^{-4}$ and implying fast orbital expansions of the moons. This high dissipation, with no obvious variations for tidal frequencies corresponding to those of Enceladus and Dione, may be explained by viscous friction in a solid core, implying a core viscosity typically ranging between $10^{14}$ and $10^{16}$ Pa.s (Remus et al., 2012). However, a dissipation increase by one order of magnitude at Rhea's frequency could suggest the existence of an additional, frequency-dependent, dissipation process, possibly from turbulent friction acting on tidal waves in the fluid envelope of Saturn (Ogilvie & Li, 2004). Alternatively, a few of Saturn's moons might themselves experience large tidal dissipation.

We present a study of the behavior of Saturn's A ring outer edge, using images and occultation data obtained by the Cassini spacecraft over a period of 8 years from 2006 to 2014. More than 5000 images and 170 occultations of the A ring outer edge are analyzed. Our fits confirm the expected response to the Janus 7:6 Inner Lindblad resonance (ILR) between 2006 and 2010, when Janus was on the inner leg of its regular orbit swap with Epimetheus. During this period, the edge exhibits a regular 7-lobed pattern with an amplitude of 12.8 km and one minimum aligned with the orbital longitude of Janus, as has been found by previous investigators. However, between 2010 and 2014, the Janus/Epimetheus orbit swap moves the Janus 7:6 LR away from the A ring outer edge, and the 7-lobed pattern disappears. In addition to several smaller-amplitudes modes, indeed, we found a variety of pattern speeds with different azimuthal wave numbers, and many of them may arise from resonant cavities between the ILR and the ring edge; also we found some other signatures consistent with tesseral resonances that could be associated with inhomogeneities in Saturn's gravity field. Moreover, these signatures do not have a fixed pattern speed. We present an analysis of these data and suggest a possible dynamical model for the behavior of the A ring's outer edge after 2010.

The shape of the probability distribution function (PDF) of molecular clouds is an important ingredient for modern theories of star formation and turbulence. Recently, several studies have pointed out observational difficulties with constraining the low column density (i.e. Av <1) PDF using dust tracers. In order to constrain the shape and properties of the low column density probability distribution function, we investigate the PDF of multiphase atomic gas in the Perseus molecular cloud using opacity-corrected GALFA-HI data and compare the PDF shape and properties to the total gas PDF and the N(H2) PDF. We find that the shape of the PDF in the atomic medium of Perseus is well described by a lognormal distribution, and not by a power-law or bimodal distribution. The peak of the atomic gas PDF in and around Perseus lies at the HI-H2 transition column density for this cloud, past which the N(H2) PDF takes on a powerlaw form. We find that the PDF of the atomic gas is narrow and at column densities larger than the HI-H2 transition the HI rapidly depletes, suggesting that the HI PDF may be used to find the HI-H2 transition column density. We also calculate the sonic Mach number of the atomic gas by using HI absorption line data, which yields a median value of Ms=4.0 for the CNM, while the HI emission PDF, which traces both the WNM and CNM, has a width more consistent with transonic turbulence.

We present the first detection of HCO$^+$ absorption in the Magellanic System. Using the Australia Telescope Compact Array (ATCA), we observed 9 extragalactic radio continuum sources behind the Magellanic System and detected HCO$^+$ absorption towards one source located behind the leading edge of the Magellanic Bridge. The detection is located at LSR velocity of $v=214.0 \pm 0.4\rm\,km\,s^{-1}$, with a full width at half maximum of $\Delta v=4.5\pm 1.0\rm\,km\,s^{-1}$ and optical depth of $\tau(\rm HCO^+)=0.10\pm 0.02$. Although there is abundant neutral hydrogen (HI) surrounding the sightline in position-velocity space, at the exact location of the absorber the HI column density is low, $<10^{20}\rm\,cm^{-2}$, and there is little evidence for dust or CO emission from Planck observations. While the origin and survival of molecules in such a diffuse environment remains unclear, dynamical events such as HI flows and cloud collisions in this interacting system likely play an important role.

We investigate the impact of high optical depth on the HI saturation observed in the Perseus molecular cloud by using Arecibo HI emission and absorption measurements toward 26 radio continuum sources. The spin temperature and optical depth of individual HI components are derived along each line-of-sight, enabling us to estimate the correction for high optical depth. We examine two different methods for the correction, Gaussian decomposition and isothermal methods, and find that they are consistent (maximum correction factor ~ 1.2) likely due to the relatively low optical depth and insignificant contribution from the diffuse radio continuum emission for Perseus. We apply the correction to the optically thin HI column density on a pixel-by-pixel basis, and find that the total HI mass increases by ~10%. Using the corrected HI column density image and far-infrared data from the IRIS Survey, we then derive the H2 column density on ~0.4 pc scales. For five dark and star-forming sub-regions, the HI surface density is uniform with Sigma_HI ~ 7-9 solar mass/pc2, in agreement with the minimum HI surface density required for shielding H2 against photodissociation. As a result, Sigma_H2/Sigma_HI and Sigma_HI+Sigma_H2 show a tight relation. Our results are consistent with predictions for H2 formation in steady state and chemical equilibrium, and suggest that H2 formation is mainly responsible for the Sigma_HI saturation in Perseus. We also compare the optically thick HI with the observed "CO-dark" gas, and find that the optically thick HI only accounts for ~20% of the "CO-dark" gas in Perseus.

We present methods and results from "21-cm Spectral Line Observations of Neutral Gas with the EVLA" (21-SPONGE), a large survey for Galactic neutral hydrogen (HI) absorption with the Karl G. Jansky Very Large Array (VLA). With the upgraded capabilities of the VLA, we reach median root-mean-square (RMS) noise in optical depth of $\sigma_{\tau}=9\times 10^{-4}$ per $0.42\rm\,km\,s^{-1}$ channel for the 31 sources presented here. Upon completion, 21-SPONGE will be the largest HI absorption survey with this high sensitivity. We discuss the observations and data reduction strategies, as well as line fitting techniques. We prove that the VLA bandpass is stable enough to detect broad, shallow lines associated with warm HI, and show that bandpass observations can be combined in time to reduce spectral noise. In combination with matching HI emission profiles from the Arecibo Observatory ($\sim3.5'$ angular resolution), we estimate excitation (or spin) temperatures ($\rm T_s$) and column densities for Gaussian components fitted to sightlines along which we detect HI absorption (30/31). We measure temperatures up to $\rm T_s\sim1500\rm\,K$ for individual lines, showing that we can probe the thermally unstable interstellar medium (ISM) directly. However, we detect fewer of these thermally unstable components than expected from previous observational studies. We probe a wide range in column density between $\sim10^{16}$ and $>10^{21}\rm\,cm^{-2}$ for individual HI clouds. In addition, we reproduce the trend between cold gas fraction and average $\rm T_s$ found by synthetic observations of a hydrodynamic ISM simulation by Kim et al. (2014). Finally, we investigate methods for estimating HI $\rm T_s$ and discuss their biases.

As we strive to understand how galaxies evolve it is crucial that we resolve physical processes and test emerging theories in nearby systems that we can observe in great detail. Our own Galaxy, the Milky Way, and the nearby Magellanic Clouds provide unique windows into the evolution of galaxies, each with its own metallicity and star formation rate. These laboratories allow us to study with more detail than anywhere else in the Universe how galaxies acquire fresh gas to fuel their continuing star formation, how they exchange gas with the surrounding intergalactic medium, and turn warm, diffuse gas into molecular clouds and ultimately stars. The $\lambda$21-cm line of atomic hydrogen (HI) is an excellent tracer of these physical processes. With the SKA we will finally have the combination of surface brightness sensitivity, point source sensitivity and angular resolution to transform our understanding of the evolution of gas in the Milky Way, all the way from the halo down to the formation of individual molecular clouds.

We present numerically-derived orbits and mass estimates for the inner Saturnian satellites, Atlas, Prometheus, Pandora, Janus and Epimetheus from a fit to 2580 new Cassini ISS astrometric observations spanning February 2004 to August 2013. The observations are provided in a supplementary table. We estimate GM_ Atlas=0.384+/-0.001 x 10^(-3)km^3s^(-2), a value 13% smaller than the previously published estimate but with an order of magnitude reduction in the uncertainty. We also find GM_ Prometheus=10.677+/-0.006x10(-3)km^3s^(-2), GM_Pandora=9.133+/-0.009x10^(-3)km^3s^(-2), GM_Janus=126.51+/-0.03x10^(-3)km^3s^(-2) and GM_Epimetheus=35.110+/-0.009x10^(-3)km^3s^(-2), consistent with previously published values, but also with significant reductions in uncertainties. We show that Atlas is currently librating in both the 54:53 co-rotation-eccentricity resonance (CER) and the 54:53 inner Lindblad (ILR) resonance with Prometheus, making it the latest example of a coupled CER-ILR system, in common with the Saturnian satellites Anthe, Aegaeon and Methone, and possibly Neptune's ring arcs. We further demonstrate that Atlas's orbit is chaotic, with a Lyapunov time of ~10 years, and show that its chaotic behaviour is a direct consequence of the coupled resonant interaction with Prometheus, rather than being an indirect effect of the known chaotic interaction between Prometheus and Pandora. We provide an updated analysis of the second-order resonant perturbations involving Prometheus, Pandora and Epimetheus based on the new observations, showing that these resonant arguments are librating only when Epimetheus is the innermost of the co-orbital pair, Janus and Epimetheus. We also find evidence that the known chaotic changes in the orbits of Prometheus and Pandora are not confined to times of apse anti-alignement.

We present astrometric observations of the Saturnian satellites Mimas, Enceladus, Tethys, Dione and Rhea from Cassini Imaging Science Subsystem (ISS) narrow-angle camera (NAC) images. Image sequences were designed to observe mutual occultations between these satellites. The positions of satellite centres were estimated by fitting ellipsoidal shape models to the measured limbs of the imaged satellites. Spacecraft pointing corrections were computed using the UCAC2 star catalogue. We provide a total of 2303 astrometric observations, resulting in 976 pairs, the remainder consisting of observations of a single satellite. Mean residuals for the individual satellite positions relative to the SAT360 ephemeris were 4.3 km in the line direction and -2.4 km in the sample direction, with standard deviations of 5.6 and 7.0 km respectively, an order of magnitude improvement in precision compared to published HST observations. By considering inter-satellite separations, uncertainties in camera pointing and spacecraft positioning along with possible biases in the individual positions of the satellites can be largely eliminated, resulting in an order-of-magnitude increase in accuracy compared to that achievable using the individual satellite positions. We show how factors relating to the viewing geometry cause small biases in the individual positions of order 0.28 pixel to become systematic across the dataset as a whole and discuss options for reducing their effects . The reduced astrometric data are provided in the form of individual positions for each satellite, together with the measured positions of reference stars, in order to allow more flexibility in the processing of the observations, taking into account possible future advances in limb-fitting techniques as well as the future availability of more accurate star catalogues, such as those from the GAIA mission.

We present a new algorithm, named Autonomous Gaussian Decomposition (AGD), for automatically decomposing spectra into Gaussian components. AGD uses derivative spectroscopy and machine learning to provide optimized guesses for the number of Gaussian components in the data, and also their locations, widths, and amplitudes. We test AGD and find that it produces results comparable to human-derived solutions on 21cm absorption spectra from the 21cm SPectral line Observations of Neutral Gas with the EVLA (21-SPONGE) survey. We use AGD with Monte Carlo methods to derive the HI line completeness as a function of peak optical depth and velocity width for the 21-SPONGE data, and also show that the results of AGD are stable against varying observational noise intensity. The autonomy and computational efficiency of the method over traditional manual Gaussian fits allow for truly unbiased comparisons between observations and simulations, and for the ability to scale up and interpret the very large data volumes from the upcoming Square Kilometer Array and pathfinder telescopes.

(Abridged) Using the Arecibo Observatory we have obtained neutral hydrogen (HI) absorption and emission spectral pairs in the direction of 26 background radio continuum sources in the vicinity of the Perseus molecular cloud. Strong absorption lines were detected in all cases allowing us to estimate spin temperature (T_s) and optical depth for 107 individual Gaussian components along these lines of sight. Basic properties of individual HI clouds (spin temperature, optical depth, and the column density of the cold and warm neutral medium, CNM and WNM) in and around Perseus are very similar to those found for random interstellar lines of sight sampled by the Millennium HI survey. This suggests that the neutral gas found in and around molecular clouds is not atypical. However, lines of sight in the vicinity of Perseus have on average a higher total HI column density and the CNM fraction, suggesting an enhanced amount of cold HI relative to an average interstellar field. Our estimated optical depth and spin temperature are in stark contrast with the recent attempt at using Planck data to estimate properties of the optically thick HI. Only ~15% of lines of sight in our study have a column density weighted average spin temperature lower than 50 K, in comparison with >85% of Planck's sky coverage. The observed CNM fraction is inversely proportional to the optical-depth weighted average spin temperature, in excellent agreement with the recent numerical simulations by Kim et al. While the CNM fraction is on average higher around Perseus relative to a random interstellar field, it is generally low, 10-50%. This suggests that extended WNM envelopes around molecular clouds, and/or significant mixing of CNM and WNM throughout molecular clouds, are present and should be considered in the models of molecule and star formation.

Remote sensing observations meet some limitations when used to study the bulk atmospheric composition of the giant planets of our solar system. A remarkable example of the superiority of in situ probe measurements is illustrated by the exploration of Jupiter, where key measurements such as the determination of the noble gases abundances and the precise measurement of the helium mixing ratio have only been made available through in situ measurements by the Galileo probe. This paper describes the main scientific goals to be addressed by the future in situ exploration of Saturn placing the Galileo probe exploration of Jupiter in a broader context and before the future probe exploration of the more remote ice giants. In situ exploration of Saturn's atmosphere addresses two broad themes that are discussed throughout this paper: first, the formation history of our solar system and second, the processes at play in planetary atmospheres. In this context, we detail the reasons why measurements of Saturn's bulk elemental and isotopic composition would place important constraints on the volatile reservoirs in the protosolar nebula. We also show that the in situ measurement of CO (or any other disequilibrium species that is depleted by reaction with water) in Saturn's upper troposphere would constrain its bulk O/H ratio. We highlight the key measurements required to distinguish competing theories to shed light on giant planet formation as a common process in planetary systems with potential applications to most extrasolar systems. In situ measurements of Saturn's stratospheric and tropospheric dynamics, chemistry and cloud-forming processes will provide access to phenomena unreachable to remote sensing studies. Different mission architectures are envisaged, which would benefit from strong international collaborations.

We use the Karl G. Jansky Very Large Array (VLA) to conduct a high-sensitivity survey of neutral hydrogen (HI) absorption in the Milky Way. In combination with corresponding HI emission spectra obtained mostly with the Arecibo Observatory, we detect a widespread warm neutral medium (WNM) component with excitation temperature <Ts>= 7200 (+1800,-1200) K (68% confidence). This temperature lies above theoretical predictions based on collisional excitation alone, implying that Ly-{\alpha} scattering, the most probable additional source of excitation, is more important in the interstellar medium (ISM) than previously assumed. Our results demonstrate that HI absorption can be used to constrain the Ly-{\alpha} radiation field, a critical quantity for studying the energy balance in the ISM and intergalactic medium yet notoriously difficult to model because of its complicated radiative transfer, in and around galaxies nearby and at high redshift.

We investigate turbulent properties of the non-star-forming, translucent molecular cloud, MBM16 by applying the statistical technique of a two-dimensional spatial power spectrum (SPS) on the neutral hydrogen (HI) observations obtained by the Galactic Arecibo L-Band Feed Array HI (GALFA-HI) survey. The SPS, calculated over the range of spatial scales from 0.1 to 17 pc, is well represented with a single power-law function, with a slope ranging from -3.3 to -3.7 and being consistent over the velocity range of MBM16 for a fixed velocity channel thickness. However, the slope varies significantly with the velocity slice thickness, suggesting that both velocity and density contribute to HI intensity fluctuations. By using this variation we estimate the slope of 3D density fluctuations in MBM16 to be -3.7\pm0.2. This is significantly steeper than what has been found for HI in the Milky Way plane, the Small Magellanic Cloud, or the Magellanic Bridge, suggesting that interstellar turbulence in MBM16 is driven on scales >17 pc and that the lack of stellar feedback could be responsible for the steep power spectrum.

Small (~50km scale), irregular features seen in Cassini images to be emanating from Saturn's F ring have been termed mini-jets by Attree et al. (2012). One particular mini-jet was tracked over half an orbital period, revealing its evolution with time and suggesting a collision with a local moonlet as its origin. In addition to these data we present here a much more detailed analysis of the full catalogue of over 800 F ring mini-jets, examining their distribution, morphology and lifetimes in order to place constraints on the underlying moonlet population. We find mini-jets randomly located in longitude around the ring, with little correlation to the moon Prometheus, and randomly distributed in time, over the full Cassini tour to date. They have a tendency to cluster together, forming complicated `multiple' structures, and have typical lifetimes of ~1d. Repeated observations of some features show significant evolution, including the creation of new mini-jets, implying repeated collisions by the same object. This suggests a population of <~1km radius objects with some internal strength and orbits spread over 100km in semi-major axis relative to the F ring but with the majority within 20km. These objects likely formed in the ring under, and were subsequently scattered onto differing orbits by, the perturbing action of Prometheus. This reinforces the idea of the F ring as a region with a complex balance between collisions, disruption and accretion.

Understanding the consequences of the gravitational interaction between a star and a planet is fundamental to the study of exoplanets. The solution of the two-body problem shows that the planet moves in an elliptical path around the star and that each body moves in an ellipse about the common center of mass. The basic properties of such a system are derived from first principles and described in the context of detecting exoplanets.

We present the first detection of a 6.7 GHz Class II methanol (CH3OH) maser in the Andromeda galaxy (M31). The CH3OH maser was found in a Very Large Array (VLA) survey during the fall of 2009. We have confirmed the methanol maser with the new Expanded VLA (EVLA), in operation since March 2010, but were unsuccessful in detecting a water maser at this location. A direct application for this methanol maser is the determination of the proper motion of M31, such as was obtained with water masers in M33 and IC10 previously. Unraveling the three-dimensional velocity of M31 would solve for the biggest unknown in the modeling of the dynamics and evolution of the Local Group of galaxies.

Aegaeon (Saturn LIII, S/2008 S1) is a small satellite of Saturn that orbits within a bright arc of material near the inner edge of Saturn's G ring. This object was observed in 21 images with Cassini's Narrow-Angle Camera between June 15 (DOY 166), 2007 and February 20 (DOY 51), 2009. If Aegaeon has similar surface scattering properties as other nearby small Saturnian satellites (Pallene, Methone and Anthe), then its diameter is approximately 500 m. Orbit models based on numerical integrations of the full equations of motion show that Aegaeon's orbital motion is strongly influenced by multiple resonances with Mimas. In particular, like the G-ring arc it inhabits, Aegaeon is trapped in the 7:6 corotation eccentricity resonance with Mimas. Aegaeon, Anthe and Methone therefore form a distinctive class of objects in the Saturn system: small moons in co-rotation eccentricity resonances with Mimas associated with arcs of debris. Comparisons among these different ring-arc systems reveal that Aegaeon's orbit is closer to the exact resonance than Anthe's and Methone's orbits are. This could indicate that Aegaeon has undergone significant orbital evolution via its interactions with the other objects in its arc, which would be consistent with the evidence that Aegaeon's mass is much smaller relative to the total mass in its arc than Anthe's and Methone's masses are.

Following the discovery that asteroid (3753) Cruithne was a coorbital companion of the Earth, a new theory of coorbital motion has been developed whereby planets or satellites can maintain companion objects in the same orbit as themselves. This has led to the prediction of hitherto unknown types of stable motion, all of which are seen in the evolution of specific near-Earth asteroids. The slow diffusion of such objects through the Earth's coorbital region is shown to lead to temporary capture, suggesting the existence of undiscovered retrograde moons of the Earth.