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.

The regularity of pulsar emissions becomes apparent once we reference the pulses' times of arrivals to the inertial rest frame of the solar system. It follows that errors in the determination of Earth's position with respect to the solar-system barycenter can appear as a time-correlated bias in pulsar-timing residual time series, affecting the searches for low-frequency gravitational waves performed with pulsar timing arrays. Indeed, recent array datasets yield different gravitational-wave background upper limits and detection statistics when analyzed with different solar-system ephemerides. Crucially, the ephemerides do not generally provide usable error representations. In this article we describe the motivation, construction, and application of a physical model of solar-system ephemeris uncertainties, which focuses on the degrees of freedom (Jupiter's orbital elements) most relevant to gravitational-wave searches with pulsar timing arrays. This model, BayesEphem, was used to derive ephemeris-robust results in NANOGrav's 11-yr stochastic-background search, and it provides a foundation for future searches by NANOGrav and other consortia. The analysis and simulations reported here suggest that ephemeris modeling reduces the gravitational-wave sensitivity of the 11-yr dataset; and that this degeneracy will vanish with improved ephemerides and with the longer pulsar timing datasets that will become available in the near future.

We report a search for astronomical neutrinos in the energy region from several GeV to TeV in the direction of the blazar TXS0506+056 using the Super-Kamiokande detector following the detection of a 100 TeV neutrino from the same location by the IceCube collaboration. Using Super-Kamiokande neutrino data across several data samples observed from April 1996 to February 2018 we have searched for both a total excess above known backgrounds across the entire period as well as localized excesses on smaller time scales in that interval. No significant excess nor significant variation in the observed event rate are found in the blazar direction. Upper limits are placed on the electron and muon neutrino fluxes at 90\% confidence level as $6.03 \times 10^{-7}$ and $4.52 \times 10^{-7}$ to $9.26 \times 10^{-10}$ [${\rm erg}/{\rm cm}^2/{\rm s}$], respectively.

Supernova detection is a major objective of the Super-Kamiokande (SK) experiment. In the next stage of SK (SK-Gd), gadolinium (Gd) sulfate will be added to the detector, which will improve the ability of the detector to identify neutrons. A core-collapse supernova will be preceded by an increasing flux of neutrinos and anti-neutrinos, from thermal and weak nuclear processes in the star, over a timescale of hours; some of which may be detected at SK-Gd. This could provide an early warning of an imminent core-collapse supernova, hours earlier than the detection of the neutrinos from core collapse. Electron anti-neutrino detection will rely on inverse beta decay events below the usual analysis energy threshold of SK, so Gd loading is vital to reduce backgrounds while maximising detection efficiency. Assuming normal neutrino mass ordering, more than 200 events could be detected in the final 12 hours before core collapse for a 15-25 solar mass star at around 200 pc, which is representative of the nearest red supergiant to Earth, $\mathrm{\alpha}$Ori (Betelgeuse). At a statistical false alarm rate of 1 per century, detection could be up to 10 hours before core collapse, and a pre-supernova star could be detected by SK-Gd up to 600 pc away. A pre-supernova alert could be provided to the astrophysics community following gadolinium loading.

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.

We search for an isotropic stochastic gravitational-wave background (GWB) in the newly released $11$-year dataset from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). While we find no significant evidence for a GWB, we place constraints on a GWB from a population of supermassive black-hole binaries, cosmic strings, and a primordial GWB. For the first time, we find that the GWB upper limits and detection statistics are sensitive to the Solar System ephemeris (SSE) model used, and that SSE errors can mimic a GWB signal. We developed an approach that bridges systematic SSE differences, producing the first PTA constraints that are robust against SSE uncertainties. We thus place a $95\%$ upper limit on the GW strain amplitude of $A_\mathrm{GWB}<1.45\times 10^{-15}$ at a frequency of $f=1$ yr$^{-1}$ for a fiducial $f^{-2/3}$ power-law spectrum, and with inter-pulsar correlations modeled. This is a factor of $\sim 2$ improvement over the NANOGrav $9$-year limit, calculated using the same procedure. Previous PTA upper limits on the GWB will need revision in light of SSE systematic uncertainties. We use our constraints to characterize the combined influence on the GWB of the stellar mass-density in galactic cores, the eccentricity of SMBH binaries, and SMBH--galactic-bulge scaling relationships. We constrain cosmic-string tension using recent simulations, yielding an SSE-marginalized $95\%$ upper limit on the cosmic string tension of $G\mu < 5.3\times 10^{-11}$---a factor of $\sim 2$ better than the published NANOGrav $9$-year constraints. Our SSE-marginalized $95\%$ upper limit on the energy density of a primordial GWB (for a radiation-dominated post-inflation Universe) is $\Omega_\mathrm{GWB}(f)h^2<3.4\times10^{-10}$.

We report the results of a neutrino search in Super-Kamiokande for coincident signals with the first detected gravitational wave produced by a binary neutron star merger, GW170817, which was followed by a short gamma-ray burst, GRB170817A, and a kilonova/macronova. We searched for coincident neutrino events in the range from 3.5 MeV to $\sim$100 PeV, in a time window $\pm$500 seconds around the gravitational wave detection time, as well as during a 14-day period after the detection. No significant neutrino signal was observed for either time window. We calculated 90% confidence level upper limits on the neutrino fluence for GW170817. From the upward-going-muon events in the energy region above 1.6 GeV, the neutrino fluence limit is $16.0^{+0.7}_{-0.6}$ ($21.3^{+1.1}_{-0.8}$) cm$^{-2}$ for muon neutrinos (muon antineutrinos), with an error range of $\pm5^{\circ}$ around the zenith angle of NGC4993, and the energy spectrum is under the assumption of an index of $-2$. The fluence limit for neutrino energies less than 100 MeV, for which the emission mechanism would be different than for higher-energy neutrinos, is also calculated. It is $6.6 \times 10^7$ cm$^{-2}$ for anti-electron neutrinos under the assumption of a Fermi-Dirac spectrum with average energy of 20 MeV.

This paper reports the first OH 18-cm line observation of the first detected interstellar object 1I/2017 U1 (`Oumuamua) using the Green Bank Telescope. We have observed the OH lines at 1665.402 MHz, 1667.359, and 1720.53 MHz frequencies with a spectral resolution of 357 Hz (approximately 0.06 km-s^{-1}). At the time of the observation, `Oumuamua was at topocentric distance and velocity of 1.07 au and 63.4 km-s^{-1}, respectively, or at heliocentric distance and velocity of 1.8 au and 39 km-s^{-1}, respectively. Based on a detailed data reduction and an analogy-based inversion, our final results confirm the asteroidal origin of `Oumuamua (as discussed in Meech et al., 2017) with an upper bound of OH production of Q[OH] < 0.17 x 10^{28} s^{-1}.

We have conducted a search for radio emission consistent with an artificial source targeting 1I/'Oumuamua with the Robert C. Byrd Green Bank Telescope (GBT) between 1.1 and 11.6 GHz. We searched the data for narrowband signals and found none. Given the close proximity to this interstellar object, we can place limits to putative transmitters with extremely low power (0.08 W).

In-spiraling supermassive black holes should emit gravitational waves, which would produce characteristic distortions in the time of arrival residuals from millisecond pulsars. Multiple national and regional consortia have constructed pulsar timing arrays by precise timing of different sets of millisecond pulsars. An essential aspect of precision timing is the transfer of the times of arrival to a (quasi-)inertial frame, conventionally the solar system barycenter. The barycenter is determined from the knowledge of the planetary masses and orbits, which has been refined over the past 50 years by multiple spacecraft. Within the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), uncertainties on the solar system barycenter are emerging as an important element of the NANOGrav noise budget. We describe what is known about the solar system barycenter, touch upon how uncertainties in it affect gravitational wave studies with pulsar timing arrays, and consider future trends in spacecraft navigation.

We present the results of a search in the Super-Kamiokande (SK) detector for excesses of neutrinos with energies above a few GeV that are in the direction of the track events reported in IceCube. Data from all SK phases (SK-I through SK-IV) were used, spanning a period from April 1996 to April 2016 and corresponding to an exposure of 225 kilotonne-years . We considered the 14 IceCube track events from a data set with 1347 livetime days taken from 2010 to 2014. We use Poisson counting to determine if there is an excess of neutrinos detected in SK in a 10 degree search cone (5 degrees for the highest energy data set) around the reconstructed direction of the IceCube event. No significant excess was found in any of the search directions we examined. We also looked for coincidences with a recently reported IceCube multiplet event. No events were detected within a $\pm$ 500 s time window around the first detected event, and no significant excess was seen from that direction over the lifetime of SK.

A comprehensive study on the atmospheric neutrino flux in the energy region from sub-GeV up to several TeV using the Super-Kamiokande water Cherenkov detector is presented in this paper. The energy and azimuthal spectra of the atmospheric ${\nu}_e+{\bar{\nu}}_e$ and ${\nu}_{\mu}+{\bar{\nu}}_{\mu}$ fluxes are measured. The energy spectra are obtained using an iterative unfolding method by combining various event topologies with differing energy responses. The azimuthal spectra depending on energy and zenith angle, and their modulation by geomagnetic effects, are also studied. A predicted east-west asymmetry is observed in both the ${\nu}_e$ and ${\nu}_{\mu}$ samples at 8.0 {\sigma} and 6.0 {\sigma} significance, respectively, and an indication that the asymmetry dipole angle changes depending on the zenith angle was seen at the 2.2 {\sigma} level. The measured energy and azimuthal spectra are consistent with the current flux models within the estimated systematic uncertainties. A study of the long-term correlation between the atmospheric neutrino flux and the solar magnetic activity cycle is also performed, and a weak indication of a correlation was seen at the 1.1 {\sigma} level, using SK I-IV data spanning a 20 year period. For particularly strong solar activity periods known as Forbush decreases, no theoretical prediction is available, but a deviation below the typical neutrino event rate is seen at the 2.4 {\sigma} level.

We report the results from a search in Super-Kamiokande for neutrino signals coincident with the first detected gravitational wave events, GW150914 and GW151226, using a neutrino energy range from 3.5 MeV to 100 PeV. We searched for coincident neutrino events within a time window of $\pm$500 seconds around the gravitational wave detection time. Four neutrino candidates are found for GW150914 and no candidates are found for GW151226. The remaining neutrino candidates are consistent with the expected background events. We calculated the 90\% confidence level upper limits on the combined neutrino fluence for both gravitational wave events, which depends on event energy and topologies. Considering the upward going muon data set (1.6 GeV - 100 PeV) the neutrino fluence limit for each gravitational wave event is 14 - 37 (19 - 50) cm$^{-2}$ for muon neutrinos (muon antineutrinos), depending on the zenith angle of the event. In the other data sets, the combined fluence limits for both gravitational wave events range from 2.4$\times 10^{4}$ to 7.0$\times 10^{9}$ cm$^{-2}$.

Previous observations suggested that Ceres has active but possibly sporadic water outgassing, and possibly varying spectral characteristics in a time scale of months. We used all available data of Ceres collected in the past three decades from the ground and the Hubble Space Telescope, and the newly acquired images by Dawn Framing Camera to search for spectral and albedo variability on Ceres, in both a global scale and local regions, particularly the bright spots inside Occator crater, over time scales of a few months to decades. Our analysis has placed an upper limit on the possible temporal albedo variation on Ceres. Sporadic water vapor venting, or any possibly ongoing activity on Ceres, is not significant enough to change the albedo or the area of the bright features in Occator crater by >15%, or the global albedo by >3% over various time scales that we searched. Recently reported spectral slope variations can be explained by changing Sun-Ceres-Earth geometry. The active area on Ceres is less than 1 km$^2$, too small to cause global albedo and spectral variations detectable in our data. Impact ejecta due to impacting projectiles of tens of meters in size like those known to cause observable changes to the surface albedo on Asteroid Scheila cannot cause detectable albedo change on Ceres due to its relatively large size and strong gravity. The water vapor activity on Ceres is independent of Ceres' heliocentric distance, rulling out the possibility of comet-like sublimation process as a possible mechanism driving the activity.

The MOdified Newtonian Dynamics (MOND) is an attempt to modify the gravitation theory to solve the Dark Matter problem. This phenomenology is very successful at the galactic level. The main effect produced by MOND in the Solar System is called the External Field Effect parametrized by the parameter $Q_2$. We have used 9 years of Cassini range and Doppler measurements to constrain $Q_2$. Our estimate of this parameter based on Cassini data is given by $Q_2=(3 \pm 3)\times 10^{-27} \ \rm{s^{-2}}$ which shows no deviation from General Relativity and excludes a large part of the relativistic MOND theories. This limit can also be interpreted as a limit on a external tidal potential acting on the Solar System coming from the internal mass of our galaxy (including Dark Matter) or from a new hypothetical body.

In this communication, the current tests of gravitation available at Solar System scales are recalled. These tests rely mainly on two frameworks: the PPN framework and the search for a fifth force. Some motivations are given to look for deviations from General Relativity in other frameworks than the two extensively considered. A recent analysis of Cassini data in a MOND framework is presented. Furthermore, possibilities to constrain Standard Model Extension parameters using Solar System data are developed.

Observations of asteroid 4 Vesta by NASA's Dawn spacecraft are interesting because its surface has the largest range of albedo, color and composition of any other asteroid visited by spacecraft to date. These hemispherical and rotational variations in surface brightness and composition have been attributed to impact processes since Vesta's formation. Prior to Dawn's arrival at Vesta, its surface properties were the focus of intense telescopic investigations for nearly a hundred years. Ground-based photometric and spectroscopic observations first revealed these variations followed later by those using Hubble Space Telescope. Here we compare interpretations of Vesta's rotation period, pole, albedo, topographic, color, and compositional properties from ground-based telescopes and HST with those from Dawn. Rotational spectral variations observed from ground-based studies are also consistent with those observed by Dawn. While the interpretation of some of these features was tenuous from past data, the interpretations were reasonable given the limitations set by spatial resolution and our knowledge of Vesta and HED meteorites at that time. Our analysis shows that ground-based and HST observations are critical for our understanding of small bodies and provide valuable support for ongoing and future spacecraft missions.

We report laboratory results of a coronagraphic test bench to assess the intensity reduction differences between a "Gaussian" tapered focal plane coronagraphic mask and a classical hard-edged "Top Hat" function mask at Extreme Adaptive Optics (ExAO) Strehl ratios of ~94%. However, unlike a traditional coronagraph design, we insert a reflective focal plane mask at 45 degree to the optical axis. We also used an intermediate secondary mask ("Mask_2") before a final image in order to block additional mask-edge diffracted light. The test bench simulates the 8.1m Gemini North telescope. It includes one spider vane, different mask radii (r= 1.9, 3.7, 7.4 lambda/D) and two types of reflective focal plane masks (hard-edged "Top Hat" and "Gaussian" tapered profiles). In order to investigate the relative performance of these competing coronagraphic designs with regard to extra-solar planet detection sensitivity, we utilize the simulation of realistic extra-solar planet populations (Nielson et al. 2006). With an appropriate translation of our laboratory results to expected telescope performance, a "Gaussian" tapered mask radius of 3.7 lambda/D with an additional mask ("Mask_2") performs best (highest planet detection sensitivity). For a full survey with this optimal design, the simulation predicts ~30% more planets detected compared to a similar sized "Top Hat" function mask with "Mask_2." Using the best design, the point contrast ratio between the stellar PSF peak and the coronagraphic PSF at 10 lambda/D (0.4" in H band if D = 8.1m) is ~10 times higher than a classical Lyot "Top Hat" coronagraph. Hence, we find a Gaussian apodized mask with an additional blocking mask is a superior (~10x higher contrast) than use of a classical Lyot coronagraph for ExAO-like Strehls.