Realistic synthetic observations of theoretical source models are essential for our understanding of real observational data. In using synthetic data, one can verify the extent to which source parameters can be recovered and evaluate how various data corruption effects can be calibrated. These studies are important when proposing observations of new sources, in the characterization of the capabilities of new or upgraded instruments, and when verifying model-based theoretical predictions in a comparison with observational data. We present the SYnthetic Measurement creator for long Baseline Arrays (SYMBA), a novel synthetic data generation pipeline for Very Long Baseline Interferometry (VLBI) observations. SYMBA takes into account several realistic atmospheric, instrumental, and calibration effects. We used SYMBA to create synthetic observations for the Event Horizon Telescope (EHT), a mm VLBI array, which has recently captured the first image of a black hole shadow. After testing SYMBA with simple source and corruption models, we study the importance of including all corruption and calibration effects. Based on two example general relativistic magnetohydrodynamics (GRMHD) model images of M87, we performed case studies to assess the attainable image quality with the current and future EHT array for different weather conditions. The results show that the effects of atmospheric and instrumental corruptions on the measured visibilities are significant. Despite these effects, we demonstrate how the overall structure of the input models can be recovered robustly after performing calibration steps. With the planned addition of new stations to the EHT array, images could be reconstructed with higher angular resolution and dynamic range. In our case study, these improvements allowed for a distinction between a thermal and a non-thermal GRMHD model based on salient features in reconstructed images.

The determination of the neutrino mass is one of the major challenges in astroparticle physics today. Direct neutrino mass experiments, based solely on the kinematics of beta-decay, provide a largely model-independent probe to the neutrino mass scale. The Karlsruhe Tritium Neutrino (KATRIN) experiment is designed to directly measure the effective electron antineutrino mass with a sensitivity of 0.2 eV 90% CL. In this work we report on the first operation of KATRIN with tritium which took place in 2018. During this commissioning phase of the tritium circulation system, excellent agreement of the theoretical prediction with the recorded spectra was found and stable conditions over a time period of 13 days could be established. These results are an essential prerequisite for the subsequent neutrino mass measurements with KATRIN in 2019.

We report on the neutrino mass measurement result from the first four-week science run of the Karlsruhe Tritium Neutrino experiment KATRIN in spring 2019. Beta-decay electrons from a high-purity gaseous molecular tritium source are energy analyzed by a high-resolution MAC-E filter. A fit of the integrated electron spectrum over a narrow interval around the kinematic endpoint at 18.57 keV gives an effective neutrino mass square value of $(-1.0^{+0.9}_{-1.1})$ eV$^2$. From this we derive an upper limit of 1.1 eV (90$\%$ confidence level) on the absolute mass scale of neutrinos. This value coincides with the KATRIN sensitivity. It improves upon previous mass limits from kinematic measurements by almost a factor of two and provides model-independent input to cosmological studies of structure formation.

Gravitational-wave observations of double compact object (DCO) mergers are providing new insights into the physics of massive stars and the evolution of binary systems. Making the most of expected near-future observations for understanding stellar physics will rely on comparisons with binary population synthesis models. However, the vast majority of simulated binaries never produce DCOs, which makes calculating such populations computationally inefficient. We present an importance sampling algorithm, STROOPWAFEL, that improves the computational efficiency of population studies of rare events, by focusing the simulation around regions of the initial parameter space found to produce outputs of interest. We implement the algorithm in the binary population synthesis code COMPAS, and compare the efficiency of our implementation to the standard method of Monte Carlo sampling from the birth probability distributions. STROOPWAFEL finds $\sim$25-200 times more DCO mergers than the standard sampling method with the same simulation size, and so speeds up simulations by up to two orders of magnitude. Finding more DCO mergers automatically maps the parameter space with far higher resolution than when using the traditional sampling. This increase in efficiency also leads to a decrease of a factor $\sim$3-10 in statistical sampling uncertainty for the predictions from the simulations. This is particularly notable for the distribution functions of observable quantities such as the black hole and neutron star chirp mass distribution, including in the tails of the distribution functions where predictions using standard sampling can be dominated by sampling noise.

We present the design and implementation of an automated data calibration and reduction pipeline for very-long-baseline interferometric (VLBI) observations taken at millimeter wavelengths. These short radio-wavelengths provide the best imaging resolution available from ground-based VLBI networks such as the Event Horizon Telescope (EHT) and the Global Millimeter VLBI Array (GMVA), but require specialized processing due to the strong effects from atmospheric opacity and turbulence as well as the heterogeneous nature of existing global arrays. The pipeline builds upon a calibration suite (HOPS) originally designed for precision geodetic VLBI. To support the reduction of data for astronomical observations, we have developed an additional framework for global phase and amplitude calibration which provides output in a standard data format for astronomical imaging and analysis. We test the pipeline on observations taken at 3.5 mm (86 GHz) by the GMVA joined by the phased Atacama Large Millimeter/submillimeter Array in April 2017, and demonstrate the benefits from the specialized processing of high frequency VLBI data with respect to classical analysis techniques.

We investigate the impact of uncertainty in the metallicity-specific star formation rate over cosmic time on predictions of the rates and masses of double compact object mergers observable through gravitational waves. We find that this uncertainty can change the predicted detectable merger rate by more than an order of magnitude, comparable to contributions from uncertain physical assumptions regarding binary evolution, such as mass transfer efficiency or supernova kicks. We statistically compare the results produced by the COMPAS population synthesis suite against a catalog of gravitational-wave detections from the first two Advanced LIGO and Virgo observing runs. We find that the rate and chirp mass of observed binary black hole mergers can be well matched under our default evolutionary model with a star formation metallicity spread of $0.39$ dex around a mean metallicity $\left<Z\right>$ that scales with redshift $z$ as $\left<Z\right>=0.035 \times 10^{-0.23 z}$, assuming a star formation rate of $0.01 \times (1+z)^{2.77} / (1+((1+z)/2.9)^{4.7}) \, \rm{M}_\odot$ Mpc$^{-3}$ yr$^{-1}$. Intriguingly, this default model predicts that 80\% of the approximately one binary black hole merger per day that will be detectable at design sensitivity will have formed through isolated binary evolution with only dynamically stable mass transfer, i.e., without experiencing a common-envelope event.

We report a framework for spectroscopic follow-up design for optimizing supernova photometric classification. The strategy accounts for the unavoidable mismatch between spectroscopic and photometric samples, and can be used even in the beginning of a new survey -- without any initial training set. The framework falls under the umbrella of active learning (AL), a class of algorithms that aims to minimize labelling costs by identifying a few, carefully chosen, objects which have high potential in improving the classifier predictions. As a proof of concept, we use the simulated data released after the Supernova Photometric Classification Challenge (SNPCC) and a random forest classifier. Our results show that, using only 12\% the number of training objects in the SNPCC spectroscopic sample, this approach is able to double purity results. Moreover, in order to take into account multiple spectroscopic observations in the same night, we propose a semi-supervised batch-mode AL algorithm which selects a set of $N=5$ most informative objects at each night. In comparison with the initial state using the traditional approach, our method achieves 2.3 times higher purity and comparable figure of merit results after only 180 days of observation, or 800 queries (73% of the SNPCC spectroscopic sample size). Such results were obtained using the same amount of spectroscopic time necessary to observe the original SNPCC spectroscopic sample, showing that this type of strategy is feasible with current available spectroscopic resources. The code used in this work is available in the COINtoolbox: https://github.com/COINtoolbox/ActSNClass .

Double neutron stars (DNSs) have been observed as Galactic radio pulsars, and the recent discovery of gravitational waves from the DNS merger GW170817 adds to the known DNS population. We perform rapid population synthesis of massive binary stars and discuss model predictions, including formation rates, mass distributions, and delay time distributions. We vary assumptions and parameters of physical processes such as mass transfer stability criteria, supernova kick distributions, remnant mass distributions and common-envelope energetics. We compute the likelihood of observing the orbital period-eccentricity distribution of the Galactic DNS population under each of our population synthesis models, allowing us to quantitatively compare the models. We find that mass transfer from a stripped post-helium-burning secondary (case BB) onto a neutron star is most likely dynamically stable. We also find that a natal kick distribution composed of both low (Maxwellian $\sigma=30\rm~km~s^{-1}$) and high ($\sigma=265\rm~km~s^{-1}$) components is preferred over a single high-kick component. We find that the observed DNS mass distribution can place strong constraints on model assumptions.

The properties of the population of merging binary black holes encode some of the uncertain physics of the evolution of massive stars in binaries. The binary black hole merger rate and chirp mass distribution are being measured by ground-based gravitational-wave detectors. We consider isolated binary evolution and explore how accurately the physical model can be constrained with such observations by applying the Fisher information matrix to the merging black hole population simulated with the rapid binary population synthesis code COMPAS. We investigate variations in four COMPAS parameters: common envelope efficiency, kick velocity dispersion, and mass loss rates during the luminous blue variable and Wolf--Rayet stellar evolutionary phases. We find that 1000 observations would constrain these model parameters to a fractional accuracy of a few percent. Given the empirically determined binary black hole merger rate, we can expect gravitational-wave observations alone to place strong constraints on the physics of stellar and binary evolution within a few years.

The nearby red giant Aldebaran is known to host a gas giant planetary companion from decades of ground-based spectroscopic radial velocity measurements. Using Gaussian Process-based Continuous Auto-Regressive Moving Average (CARMA) models, we show that these historic data also contain evidence of acoustic oscillations in the star itself, and verify this result with further dedicated ground-based spectroscopy and space-based photometry with the Kepler Space Telescope. From the frequency of these oscillations we determine the mass of Aldebaran to be $1.16 \pm 0.07 \, M_\odot$, and note that this implies its planet will have been subject to insolation comparable to the Earth for some of the star's main sequence lifetime. Our approach to sparse, irregularly sampled time series astronomical observations has the potential to unlock asteroseismic measurements for thousands of stars in archival data, and push to lower-mass planets around red giant stars.

As we enter the era of gravitational wave astronomy, we are beginning to collect observations which will enable us to explore aspects of astrophysics of massive stellar binaries which were previously beyond reach. In this paper we describe COMPAS (Compact Object Mergers: Population Astrophysics and Statistics), a new platform to allow us to deepen our understanding of isolated binary evolution and the formation of gravitational-wave sources. We describe the computational challenges associated with their exploration, and present preliminary results on overcoming them using Gaussian process regression as a simulation emulation technique.

During its first 4 months of taking data, Advanced LIGO has detected gravitational waves from two binary black hole mergers, GW150914 and GW151226, along with the statistically less significant binary black hole merger candidate LVT151012. We use our rapid binary population synthesis code COMPAS to show that all three events can be explained by a single evolutionary channel -- classical isolated binary evolution via mass transfer including a common envelope phase. We show all three events could have formed in low-metallicity environments (Z = 0.001) from progenitor binaries with typical total masses $\gtrsim 160 M_\odot$, $\gtrsim 60 M_\odot$ and $\gtrsim 90 M_\odot$, for GW150914, GW151226, and LVT151012, respectively.

Molecules with hyperfine splitting of their rotational line spectra are useful probes of optical depth, via the relative line strengths of their hyperfine components.The hyperfine splitting is particularly advantageous in interpreting the physical conditions of the emitting gas because with a second rotational transition, both gas density and temperature can be derived. For HCN however, the relative strengths of the hyperfine lines are anomalous. They appear in ratios which can vary significantly from source to source, and are inconsistent with local thermodynamic equilibrium. This is the HCN hyperfine anomaly, and it prevents the use of simple LTE models of HCN emission to derive reliable optical depths. In this paper we demonstrate how to model HCN hyperfine line emission, and derive accurate line ratios, spectral line shapes and optical depths. We show that by carrying out radiative transfer calculations over each hyperfine level individually, as opposed to summing them over each rotational level, the anomalous hyperfine emission emerges naturally. To do this requires not only accurate radiative rates between hyperfine states, but also accurate collisional rates. We investigate the effects of different sets of hyperfine collisional rates, derived via the 'proportional method' and through direct recoupling calculations. Through an extensive parameter sweep over typical low mass star forming conditions, we show the HCN line ratios to be highly variable to optical depth. We also reproduce an observed effect whereby the red-blue asymmetry of the hyperfine lines (an infall signature) switches sense within a single rotational transition.

We report development of a simple and affordable radio interferometer suitable as an educational laboratory experiment. With the increasing importance of interferometry in astronomy, the lack of educational interferometers is an obstacle to training the future generation of astronomers. This interferometer provides the hands-on experience needed to fully understand the basic concepts of interferometry. The design of this interferometer is based on the Michelson & Pease stellar optical interferometer, but operates at a radio wavelength (~11 GHz; ~2.7cm); thus the requirement for optical accuracy is much less stringent. We utilize a commercial broadcast satellite dish and feedhorn. Two flat side mirrors slide on a ladder, providing baseline coverage. This interferometer resolves and measures the diameter of the Sun, a nice daytime experiment which can be carried out even in marginal weather (i.e., partial cloud cover). Commercial broadcast satellites provide convenient point sources for comparison to the Sun's extended disk. We describe the mathematical background of the adding interferometer, the design and development of the telescope and receiver system, and measurements of the Sun. We present results from a students' laboratory report.

The Submillimeter Array (SMA), a collaboration between the Smithsonian Astrophysical Observatory and the Academica Sinica Institute for Astronomy and Astrophysics of Taiwan, is an eight-element radio-interferometer designed to operate throughout the major atmospheric windows from about 180 to 900 GHz. In an effort to mitigate the effects of atmospheric instabilities which limit the phase coherence of the array especially in the higher frequency bands, the array was designed to allow simultaneous operation of a low frequency receiver (<350 GHz) with a high frequency receiver (>330 GHz). The overlap region of 330-350 GHz was included to facilitate dual polarization measurements in the frequency range considered to offer the highest sensitivity for continuum observations with the array. So far, the array is equipped with working SIS receivers covering the frequency ranges 176-256 GHz, 260-350 GHz, and 600-700 GHz, and single frequency operation has been routine in the lower two frequency bands for the past year. More recently, with the completion of IF hardware required to make full use of the SMA cross-correlator, dual receiver operation became possible. We have since made a number of Galactic and extra-galactic astronomical observations in dual-band mode with the hopes of using the 230 GHz receiver as a phase reference to enable improved interferometry in the 650 GHz band. We will present the current antenna and receiver performance, some of the first interferometric images in the 650 GHz receiver band, and our initial attempts at phase transfer.

We present CO observations of a large sample of ultraluminous IR galaxies out to z = 0.3. Most of the galaxies are interacting, but not completed mergers. All but one have high CO(1-0) luminosities, log(Lco [K-km/s-pc^2]) = 9.92 +/- 0.12. The dispersion in Lco is only 30%, less than that in the FIR luminosity. The integrated CO intensity correlates Strongly with the 100 micron flux density, as expected for a black body model in which the mid and far IR radiation are optically thick. We use this model to derive sizes of the FIR and CO emitting regions and the enclosed dynamical masses. Both the IR and CO emission originate in regions a few hundred parsecs in radius. The median value of Lfir/Lco = 160 Lsun/(K-km/s-pc^2), within a factor of two of the black body limit for the observed FIR temperatures. The entire ISM is a scaled up version of a normal galactic disk with densities a factor of 100 higher, making even the intercloud medium a molecular region. Using three different techniques of H2 mass estimation, we conclude that the ratio of gas mass to Lco is about a factor of four lower than for Galactic molecular clouds, but that the gas mass is a large fraction of the dynamical mass. Our analysis of CO emission reduces the H2 mass from previous estimates of 2-5e10 Msun to 0.4-1.5e10 Msun, which is in the range found for molecular gas rich spiral galaxies. A collision involving a molecular gas rich spiral could lead to an ultraluminous galaxy powered by central starbursts triggered by the compression of infalling preexisting GMC's.