We report on our observing campaign of the compact binary merger GW190814, detected by the Advanced LIGO and Advanced Virgo detectors on August 14th, 2019. This signal has the best localisation of any observed gravitational wave (GW) source, with a 90% probability area of 18.5 deg$^2$, and an estimated distance of ~ 240 Mpc. We obtained wide-field observations with the Deca-Degree Optical Transient Imager (DDOTI) covering 88% of the probability area down to a limiting magnitude of $w$ = 19.9 AB. Nearby galaxies within the high probability region were targeted with the Lowell Discovery Telescope (LDT), whereas promising candidate counterparts were characterized through multi-colour photometry with the Reionization and Transients InfraRed (RATIR) and spectroscopy with the Gran Telescopio de Canarias (GTC). We use our optical and near-infrared limits in conjunction with the upper limits obtained by the community to constrain the possible electromagnetic counterparts associated with the merger. A gamma-ray burst seen along its jet's axis is disfavoured by the multi-wavelength dataset, whereas the presence of a burst seen at larger viewing angles is not well constrained. Although our observations are not sensitive to a kilonova similar to AT2017gfo, we can rule out high-mass (> 0.1 M$_{\odot}$) fast-moving (mean velocity >= 0.3c) wind ejecta for a possible kilonova associated with this merger.

We derive limits on any electromagnetic counterpart to the compact binary merger S190814bv, whose parameters are consistent with the merger of a black hole and a neutron star. We present observations with the new wide-field optical imager DDOTI and also consider Swift/BAT observations reported by Palmer et al. (2019). We show that Swift/BAT would have detected a counterpart with similar properties to a typical on-axis short GRB at the 98 per cent confidence level, whereas our DDOTI observations only rule out such a counterpart at the 27 per cent confidence level. Neither have sufficient sensitivity to rule out an off-axis counterpart like GW 170817. We compare the efficiency of Swift/BAT and DDOTI for future observations, and show that DDOTI is likely to be about twice as efficient as Swift/BAT for off-axis events up to about 100 Mpc.

We present the results from four stellar occultations by (486958) Arrokoth, the flyby target of the New Horizons extended mission. Three of the four efforts led to positive detections of the body, and all constrained the presence of rings and other debris, finding none. Twenty-five mobile stations were deployed for 2017 June 3 and augmented by fixed telescopes. There were no positive detections from this effort. The event on 2017 July 10 was observed by SOFIA with one very short chord. Twenty-four deployed stations on 2017 July 17 resulted in five chords that clearly showed a complicated shape consistent with a contact binary with rough dimensions of 20 by 30 km for the overall outline. A visible albedo of 10% was derived from these data. Twenty-two systems were deployed for the fourth event on 2018 Aug 4 and resulted in two chords. The combination of the occultation data and the flyby results provides a significant refinement of the rotation period, now estimated to be 15.9380 $\pm$ 0.0005 hours. The occultation data also provided high-precision astrometric constraints on the position of the object that were crucial for supporting the navigation for the New Horizons flyby. This work demonstrates an effective method for obtaining detailed size and shape information and probing for rings and dust on distant Kuiper Belt objects as well as being an important source of positional data that can aid in spacecraft navigation that is particularly useful for small and distant bodies.

Inverse problems play a key role in modern image/signal processing methods. However, since they are generally ill-conditioned or ill-posed due to lack of observations, their solutions may have significant intrinsic uncertainty. Analysing and quantifying this uncertainty is very challenging, particularly in high-dimensional problems and problems with non-smooth objective functionals (e.g. sparsity-promoting priors). In this article, a series of strategies to visualise this uncertainty are presented, e.g. highest posterior density credible regions, and local credible intervals (cf. error bars) for individual pixels and superpixels. Our methods support non-smooth priors for inverse problems and can be scaled to high-dimensional settings. Moreover, we present strategies to automatically set regularisation parameters so that the proposed uncertainty quantification (UQ) strategies become much easier to use. Also, different kinds of dictionaries (complete and over-complete) are used to represent the image/signal and their performance in the proposed UQ methodology is investigated.

Until recently mass-mapping techniques for weak gravitational lensing convergence reconstruction have lacked a principled statistical framework upon which to quantify reconstruction uncertainties, without making strong assumptions of Gaussianity. In previous work we presented a sparse hierarchical Bayesian formalism for convergence reconstruction that addresses this shortcoming. Here, we draw on the concept of local credible intervals (cf. Bayesian error bars) as an extension of the uncertainty quantification techniques previously detailed. These uncertainty quantification techniques are benchmarked against those recovered via Px-MALA - a state of the art proximal Markov Chain Monte Carlo (MCMC) algorithm. We find that typically our recovered uncertainties are everywhere conservative, of similar magnitude and highly correlated (Pearson correlation coefficient $\geq 0.85$) with those recovered via Px-MALA. Moreover, we demonstrate an increase in computational efficiency of $\mathcal{O}(10^6)$ when using our sparse Bayesian approach over MCMC techniques. This computational saving is critical for the application of Bayesian uncertainty quantification to large-scale stage IV surveys such as LSST and Euclid.

We propose a Bayesian uncertainty quantification method for large-scale imaging inverse problems. Our method applies to all Bayesian models that are log-concave, where maximum-a-posteriori (MAP) estimation is a convex optimization problem. The method is a framework to analyse the confidence in specific structures observed in MAP estimates (e.g., lesions in medical imaging, celestial sources in astronomical imaging), to enable using them as evidence to inform decisions and conclusions. Precisely, following Bayesian decision theory, we seek to assert the structures under scrutiny by performing a Bayesian hypothesis test that proceeds as follows: firstly, it postulates that the structures are not present in the true image, and then seeks to use the data and prior knowledge to reject this null hypothesis with high probability. Computing such tests for imaging problems is generally very difficult because of the high dimensionality involved. A main feature of this work is to leverage probability concentration phenomena and the underlying convex geometry to formulate the Bayesian hypothesis test as a convex problem, that we then efficiently solve by using scalable optimization algorithms. This allows scaling to high-resolution and high-sensitivity imaging problems that are computationally unaffordable for other Bayesian computation approaches. We illustrate our methodology, dubbed BUQO (Bayesian Uncertainty Quantification by Optimization), on a range of challenging Fourier imaging problems arising in astronomy and medicine.

Uncertainty quantification is a critical missing component in radio interferometric imaging that will only become increasingly important as the big-data era of radio interferometry emerges. Statistical sampling approaches to perform Bayesian inference, like Markov Chain Monte Carlo (MCMC) sampling, can in principle recover the full posterior distribution of the image, from which uncertainties can then be quantified. However, for massive data sizes, like those anticipated from the Square Kilometre Array (SKA), it will be difficult if not impossible to apply any MCMC technique due to its inherent computational cost. We formulate Bayesian inference problems with sparsity-promoting priors (motivated by compressive sensing), for which we recover maximum a posteriori (MAP) point estimators of radio interferometric images by convex optimisation. Exploiting recent developments in the theory of probability concentration, we quantify uncertainties by post-processing the recovered MAP estimate. Three strategies to quantify uncertainties are developed: (i) highest posterior density credible regions; (ii) local credible intervals (cf. error bars) for individual pixels and superpixels; and (iii) hypothesis testing of image structure. These forms of uncertainty quantification provide rich information for analysing radio interferometric observations in a statistically robust manner. Our MAP-based methods are approximately $10^5$ times faster computationally than state-of-the-art MCMC methods and, in addition, support highly distributed and parallelised algorithmic structures. For the first time, our MAP-based techniques provide a means of quantifying uncertainties for radio interferometric imaging for realistic data volumes and practical use, and scale to the emerging big-data era of radio astronomy.

Uncertainty quantification is a critical missing component in radio interferometric imaging that will only become increasingly important as the big-data era of radio interferometry emerges. Since radio interferometric imaging requires solving a high-dimensional, ill-posed inverse problem, uncertainty quantification is difficult but also critical to the accurate scientific interpretation of radio observations. Statistical sampling approaches to perform Bayesian inference, like Markov Chain Monte Carlo (MCMC) sampling, can in principle recover the full posterior distribution of the image, from which uncertainties can then be quantified. However, traditional high-dimensional sampling methods are generally limited to smooth (e.g. Gaussian) priors and cannot be used with sparsity-promoting priors. Sparse priors, motivated by the theory of compressive sensing, have been shown to be highly effective for radio interferometric imaging. In this article proximal MCMC methods are developed for radio interferometric imaging, leveraging proximal calculus to support non-differential priors, such as sparse priors, in a Bayesian framework. Furthermore, three strategies to quantify uncertainties using the recovered posterior distribution are developed: (i) local (pixel-wise) credible intervals to provide error bars for each individual pixel; (ii) highest posterior density credible regions; and (iii) hypothesis testing of image structure. These forms of uncertainty quantification provide rich information for analysing radio interferometric observations in a statistically robust manner.

We performed an analysis of all RXTE observations of the Low Mass X-ray Binary and Black Hole Candidate IGR J17091-3624 during the 2011-2013 outburst of the source. By creating lightcurves, hardness-intensity diagrams and power density spectra of each observation, we have created a set of 9 variability `classes' that phenomenologically describe the range of types of variability seen in this object. We compare our set of variability classes to those established by Belloni et al. (2000) to describe the similar behaviour of the LMXB GRS 1915+105, finding that some types of variability seen in IGR J17091-3624 are not represented in data of GRS 1915+105. We also use all available X-ray data of the 2011-2013 outburst of IGR J17091-3624 to analyse its long-term evolution, presenting the first detection of IGR J17091-3624 above 150 keV as well as noting the presence of `re-flares' during the latter stages of the outburst. Using our results we place new constraints on the mass and distance of the object, and find that it accretes at <33% of its Eddington limit. As such, we conclude that Eddington-limited accretion can no longer be considered a sufficient or necessary criterion for GRS 1915+105-like variability to occur in Low Mass X-Ray Binaries.

We analyze the internal kinematics of 26 Planetary Nebulae (PNe) with low metallicity that appear to derive from progenitor stars of the lowest masses, including the halo PN population. Based upon spatially-resolved, long-slit, echelle spectroscopy drawn from the San Pedro M\'artir Kinematic Catalogue of PNe (L\'opez et al. 2012), we characterize the kinematics of these PNe measuring their global expansion velocities based upon the largest sample used to date for this purpose. We find kinematics that follow the trends observed and predicted in other studies, but also find that most of the PNe studied here tend to have expansion velocities less than 20 km/s in all of the emission lines considered. The low expansion velocities that we observe in this sample of low metallicity planetary nebulae with low mass progenitors are most likely a consequence of a weak central star wind driving the kinematics of the nebular shell. This study complements previous results (Pereyra et al. 2013, and references therein) that link the expansion velocities of the PN shells with the characteristics of the central star.

We have selected a group of 100 evolved planetary nebulae (PNe) and study their kinematics based upon spatially-resolved, long-slit, echelle spectroscopy. The data have been drawn from the San Pedro M\'artir Kinematic Catalogue of PNe (L\'opez et al. 2012). The aim is to characterize in detail the global kinematics of PNe at advanced stages of evolution with the largest sample of homogenous data used to date for this purpose. The results reveal two groups that share kinematics, morphology, and photo-ionization characteristics of the nebular shell and central star luminosities at the different late stages under study.The typical flow velocities we measure are usually larger than seen in earlier evolutionary stages, with the largest velocities occurring in objects with very weak or absent [N II] \lambda 6584 line emission, by all indications the least evolved objects in our sample. The most evolved objects expand more slowly. This apparent deceleration during the final stage of PNe evolution is predicted by hydrodynamical models, but other explanations are also possible. These results provide a template for comparison with the predictions of theoretical models.

We study the line widths in the [\ion{O}{3}]$\lambda$5007 and H$\alpha$ lines for two groups of planetary nebulae in the Milky Way bulge based upon spectroscopy obtained at the Observatorio Astron\'omico Nacional in the Sierra San Pedro M\'artir (OAN-SPM) using the Manchester Echelle Spectrograph. The first sample includes objects early in their evolution, having high H$\beta$ luminosities, but [\ion{O}{3}]$\lambda 5007/\mathrm H\beta < 3$. The second sample comprises objects late in their evolution, with \ion{He}{2} $\lambda 4686/\mathrm H\beta > 0.5$. These planetary nebulae represent evolutionary phases preceeding and following those of the objects studied by Richer et al. (2008). Our sample of planetary nebulae with weak [\ion{O}{3}]$\lambda$5007 has a line width distribution similar to that of the expansion velocities of the envelopes of AGB stars, and shifted to systematically lower values as compared to the less evolved objects studied by Richer et al. (2008). The sample with strong \ion{He}{2} $\lambda 4686$ has a line width distribution indistinguishable from that of the more evolved objects from Richer et al. (2008), but a distribution in angular size that is systematically larger and so they are clearly more evolved. These data and those of Richer et al. (2008) form a homogeneous sample from a single Galactic population of planetary nebulae, from the earliest evolutionary stages until the cessation of nuclear burning in the central star. They confirm the long-standing predictions of hydrodynamical models of planetary nebulae, where the kinematics of the nebular shell are driven by the evolution of the central star.

We present a systematic study of line widths in the [\ion{O}{3}]$\lambda$5007 and H$\alpha$ lines for a sample of 86 planetary nebulae in the Milky Way bulge based upon spectroscopy obtained at the \facility{Observatorio Astron\'omico Nacional in the Sierra San Pedro M\'artir (OAN-SPM)} using the Manchester Echelle Spectrograph. The planetary nebulae were selected with the intention of simulating samples of bright extragalactic planetary nebulae. We separate the planetary nebulae into two samples containing cooler and hotter central stars, defined by the absence or presence, respectively, of the \ion{He}{2} $\lambda$6560 line in the H$\alpha$ spectra. This division separates samples of younger and more evolved planetary nebulae. The sample of planetary nebulae with hotter central stars has systematically larger line widths, larger radii, lower electron densities, and lower H$\beta$ luminosities. The distributions of these parameters in the two samples all differ at significance levels exceeding 99%. These differences are all in agreement with the expectations from hydrodynamical models, but for the first time confirmed for a homogeneous and statistically significant sample of galactic planetary nebulae. We interpret these differences as evidence for the acceleration of the nebular shells during the early evolution of these intrinsically bright planetary nebulae. As is the case for planetary nebulae in the Magellanic Clouds, the acceleration of the nebular shells appears to be the direct result of the evolution of the central stars.