We report velocity-delay maps for prominent broad emission lines, Ly_alpha, CIV, HeII and H_beta, in the spectrum of NGC5548. The emission-line responses inhabit the interior of a virial envelope. The velocity-delay maps reveal stratified ionization structure. The HeII response inside 5-10 light-days has a broad single-peaked velocity profile. The Ly_alpha, CIV, and H_beta responses peak inside 10 light-days, extend outside 20 light-days, and exhibit a velocity profile with two peaks separated by 5000 km/s in the 10 to 20 light-day delay range. The velocity-delay maps show that the M-shaped lag vs velocity structure found in previous cross-correlation analysis is the signature of a Keplerian disk with a well-defined outer edge at R=20 light-days. The outer wings of the M arise from the virial envelope, and the U-shaped interior of the M is the lower half of an ellipse in the velocity-delay plane. The far-side response is weaker than that from the near side, so that we see clearly the lower half, but only faintly the upper half, of the velocity--delay ellipse. The delay tau=(R/c)(1-sin(i))=5 light-days at line center is from the near edge of the inclined ring, giving the inclination i=45 deg. A black hole mass of M=7x10^7 Msun is consistent with the velocity-delay structure. A barber-pole pattern with stripes moving from red to blue across the CIV and possibly Ly_alpha line profiles suggests the presence of azimuthal structure rotating around the far side of the broad-line region and may be the signature of precession or orbital motion of structures in the inner disk. Further HST observations of NGC 5548 over a multi-year timespan but with a cadence of perhaps 10 days rather than 1 day could help to clarify the nature of this new AGN phenomenon.

This is the final report submitted to NASA for a Probe-class concept study of the "Farside Array for Radio Science Investigations of the Dark ages and Exoplanets" (FARSIDE), a low radio frequency interferometric array on the farside of the Moon. The design study focused on the instrument, a deployment rover, the lander and base station, and delivered an architecture broadly consistent with the requirements for a Probe mission. This notional architecture consists of 128 dipole antennas deployed across a 10 km area by a rover, and tethered to a base station for central processing, power and data transmission to the Lunar Gateway, or an alternative relay satellite. FARSIDE would provide the capability to image the entire sky each minute in 1400 channels spanning frequencies from 150 kHz to 40 MHz, extending down two orders of magnitude below bands accessible to ground-based radio astronomy. The lunar farside can simultaneously provide isolation from terrestrial radio frequency interference, auroral kilometric radiation, and plasma noise from the solar wind. This would enable near-continuous monitoring of the nearest stellar systems in the search for the radio signatures of coronal mass ejections and energetic particle events, and would also detect the magnetospheres for the nearest candidate habitable exoplanets. Simultaneously, FARSIDE would be used to characterize similar activity in our own solar system, from the Sun to the outer planets. Through precision calibration via an orbiting beacon, and exquisite foreground characterization, FARSIDE would also measure the Dark Ages global 21-cm signal at redshifts from 50-100. It will also be a pathfinder for a larger 21-cm power spectrum instrument by carefully measuring the foreground with high dynamic range.

We present the results of the first transient survey from the Owens Valley Radio Observatory Long Wavelength Array (OVRO-LWA) using 31 hr of data, in which we place the most constraining limits on the instantaneous transient surface density at timescales of 13 s to a few minutes and at frequencies below 100 MHz. The OVRO-LWA is a dipole array that images the entire viewable hemisphere with 58 MHz of bandwidth from 27 to 84 MHz at 13 s cadence. No transients are detected above a 6.5$\sigma$ flux density limit of 10.5 Jy, implying an upper limit to the transient surface density of $2.5\times10^{-8}$ deg$^{-2}$ at the shortest timescales probed, which is orders of magnitude deeper than has been achieved at sub-100 MHz frequencies and comparable flux densities to date. The nondetection of transients in the OVRO-LWA survey, particularly at minutes-long timescales, allows us to place further constraints on the rate of the potential population of transients uncovered by Stewart et al. (2016). From their transient rate, we expect a detection of $8.4^{+31.8}_{-8.0}$ events, and the probability of our null detection is $1.9^{+644}_{-1.9}\times10^{-3}$, ruling out a transient rate $>1.4\times10^{-4} \text{days}^{-1} \text{deg}^{-2}$ with 95% confidence at a flux density limit of 18.1 Jy, under the assumption of a flat spectrum and wide bandwidth. We discuss the implications of our nondetection for this population and further constraints that can be made on the source spectral index, intrinsic emission bandwidth, and resulting luminosity distribution.

We present the discovery of a nuclear transient with the Caltech-NRAO Stripe 82 Survey (CNSS), a dedicated radio transient survey carried out with the Karl G. Jansky Very Large Array (VLA). This transient, CNSS J001947.3+003527, exhibited a turn-on over a timescale of $\lesssim$1 yr, increasing in flux density at 3 GHz from $<0.14$ mJy in February 2014 to $4.4\pm0.1$ mJy in March 2015, reaching a peak luminosity of $5\times10^{28} \text{erg s}^{-1} \text{Hz}^{-1}$ around October 2015. The association of CNSS J0019+00 with the nucleus (Gaia and our VLBI positions are consistent to within 1 pc) of a nearby S0 Seyfert galaxy at 77 Mpc, together with the radio spectral evolution, implies that this transient is most likely a tidal disruption event (TDE). Our equipartition analysis indicates the presence of a $\sim$15,000 km s$^{-1}$ outflow, having energy $\sim$10$^{49}$ erg. We derive the radial density profile for the circum-nuclear material in the host galaxy to be proportional to R$^{-2.5}$. All of these properties suggest resemblance with radio-detected thermal TDEs like ASASSN-14li and XMMSL1 J0740-85. No significant X-ray or optical emission is detected from CNSS J0019+00, although this may simply be due to the thermal emission being weak during our late-time follow up observations. From the CNSS survey we have obtained the first unbiased measurement of the rate of radio TDEs, $R(>500 \mu{\rm Jy})$ of about $2\times10^{-3}$ deg$^{-2}$, or equivalently a volumetric rate of about 10 Gpc$^{-3}$ yr$^{-1}$. This rate implies that all-sky radio surveys such as the VLA Sky Survey and those planned with ASKAP, will find many tens of radio TDEs over the next few years.

The Spectral Imaging of the Coronal Environment (SPICE) instrument is a high-resolution imaging spectrometer operating at extreme ultraviolet (EUV) wavelengths. In this paper, we present the concept, design, and pre-launch performance of this facility instrument on the ESA/NASA Solar Orbiter mission. The goal of this paper is to give prospective users a better understanding of the possible types of observations, the data acquisition, and the sources that contribute to the instrument's signal. The paper discusses the science objectives, with a focus on the SPICE-specific aspects, before presenting the instrument's design, including optical, mechanical, thermal, and electronics aspects. This is followed by a characterisation and calibration of the instrument's performance. The paper concludes with descriptions of the operations concept and data processing. The performance measurements of the various instrument parameters meet the requirements derived from the mission's science objectives. The SPICE instrument is ready to perform measurements that will provide vital contributions to the scientific success of the Solar Orbiter mission.

A successful ground array Radio Frequency (RF)-only self-trigger on 10 high-energy cosmic ray events is demonstrated with 256 dual-polarization antennas of the Owens Valley Radio Observatory Long Wavelength Array (OVRO-LWA). This RF-only capability is predicated on novel techniques for Radio Frequency Interference (RFI) identification and mitigation with an analysis efficiency of 45\% for shower-driven events with a Signal-to-noise ratio $\gtrsim$ 5 against the galactic background noise power of individual antennas. This technique enables more efficient detection of cosmic rays over a wider range of zenith angles than possible via triggers from in-situ particle detectors and can be easily adapted by neutrino experiments relying on RF-only detection. This paper discusses the system design, RFI characterization and mitigation techniques, and initial results from 10 cosmic ray events identified within a 40-hour observing window. A design for a future optimized commensal cosmic-ray detector for the OVRO-LWA is presented, as well as recommendations for developing a similar capability for other experiments -- these designs either reduce data-rate or increase sensitivity by an order of magnitude for many configurations of radio instruments.

We model the ultraviolet spectra of the Seyfert 1 galaxy NGC~5548 obtained with the Hubble Space Telescope during the 6-month reverberation-mapping campaign in 2014. Our model of the emission from NGC 5548 corrects for overlying absorption and deblends the individual emission lines. Using the modeled spectra, we measure the response to continuum variations for the deblended and absorption-corrected individual broad emission lines, the velocity-dependent profiles of Ly$\alpha$ and C IV, and the narrow and broad intrinsic absorption features. We find that the time lags for the corrected emission lines are comparable to those for the original data. The velocity-binned lag profiles of Ly$\alpha$ and C IV have a double-peaked structure indicative of a truncated Keplerian disk. The narrow absorption lines show delayed response to continuum variations corresponding to recombination in gas with a density of $\sim 10^5~\rm cm^{-3}$. The high-ionization narrow absorption lines decorrelate from continuum variations during the same period as the broad emission lines. Analyzing the response of these absorption lines during this period shows that the ionizing flux is diminished in strength relative to the far-ultraviolet continuum. The broad absorption lines associated with the X-ray obscurer decrease in strength during this same time interval. The appearance of X-ray obscuration in $\sim\,2012$ corresponds with an increase in the luminosity of NGC 5548 following an extended low state. We suggest that the obscurer is a disk wind triggered by the brightening of NGC 5548 following the decrease in size of the broad-line region during the preceding low-luminosity state.

FARSIDE (Farside Array for Radio Science Investigations of the Dark ages and Exoplanets) is a Probe-class concept to place a low radio frequency interferometric array on the farside of the Moon. A NASA-funded design study, focused on the instrument, a deployment rover, the lander and base station, delivered an architecture broadly consistent with the requirements for a Probe mission. This notional architecture consists of 128 dual polarization antennas deployed across a 10 km area by a rover, and tethered to a base station for central processing, power and data transmission to the Lunar Gateway. FARSIDE would provide the capability to image the entire sky each minute in 1400 channels spanning frequencies from 100 kHz to 40 MHz, extending down two orders of magnitude below bands accessible to ground-based radio astronomy. The lunar farside can simultaneously provide isolation from terrestrial radio frequency interference, auroral kilometric radiation, and plasma noise from the solar wind. This would enable near-continuous monitoring of the nearest stellar systems in the search for the radio signatures of coronal mass ejections and energetic particle events, and would also detect the magnetospheres for the nearest candidate habitable exoplanets. Simultaneously, FARSIDE would be used to characterize similar activity in our own solar system, from the Sun to the outer planets, including the hypothetical Planet Nine. Through precision calibration via an orbiting beacon, and exquisite foreground characterization, FARSIDE would also measure the Dark Ages global 21-cm signal at redshifts z=50-100. The unique observational window offered by FARSIDE would enable an abundance of additional science ranging from sounding of the lunar subsurface to characterization of the interstellar medium in the solar system neighborhood.

The 21\,cm transition of neutral hydrogen is opening an observational window into the cosmic dawn of the universe---the epoch of first star formation. We use 28\,hr of data from the Owens Valley Radio Observatory Long Wavelength Array (OVRO-LWA) to place upper limits on the spatial power spectrum of 21\,cm emission at $z \approx 18.4$ ($\Delta_{21} \lesssim 10^4\,\text{mK}$), and within the absorption feature reported by the EDGES experiment (Bowman et al. 2018). In the process we demonstrate the first application of the double Karhunen-Lo\`{e}ve transform for foreground filtering, and diagnose the systematic errors that are currently limiting the measurement. We also provide an updated model for the angular power spectrum of low-frequency foreground emission measured from the northern hemisphere, which can be used to refine sensitivity forecasts for next-generation experiments.

New CCD photometry has been combined with published and unpublished earlier observations to study the three Cepheid variables in M13: V1, V2 and V6. The light curve characteristics in $B$, $V$ and $I_{\rm C}$ have been determined and the periods updated. A period change analysis shows all three stars have increasing periods but for V1 and V2 the rate of period increase does not appear to be constant over the 118 years of observation. The observed rates of period increase are in good agreement with the predictions of the Pisa theoretical models with helium abundance $Y = 0.25$. Theory suggests V1 and V6 have masses of $\sim0.57\,M_\odot$ and are in the redward-evolving final stage of the ``blue loop'' evolutionary phase that is produced when helium-shell ignition occurs. The larger period and period change rate for V2 indicate it has a mass of $\sim0.52\,M_\odot$. A study of eighteen metal-poor BL Her stars shows the observed period changes for such objects in general can be reasonably well explained using the predictions from horizontal branch evolutionary tracks. BL Her stars with periods less than $\sim$3 d and relatively large secular period change rates ($dP/dt\approx5-15$ d/Myr) are in the evolutionary stage before He-shell ignition; the remaining cases are stars that have already experienced He-shell ignition. Moreover, an analysis of crossing time through the instability strip indicates that it is likely that few, if any, BL Her stars have a He abundance as large as $Y = 0.33$.

Multimessenger observations of the binary neutron star merger GW170817 have enabled the discovery of a diverse array of electromagnetic counterparts to compact binary mergers, including an unambiguous kilonova, a short gamma-ray burst, and a late-time radio jet. Beyond these counterparts, compact binary mergers are additionally predicted to be accompanied by prompt low-frequency radio emission. The successful observation of a prompt radio counterpart would be immensely valuable, but is made difficult by the short delay between the gravitational-wave and prompt electromagnetic signals as well as the poor localization of gravitational-wave sources. Here, we present the first search for prompt radio emission accompanying a gravitational-wave event, targeting the binary black hole merger GW170104 detected by the Advanced LIGO and Virgo gravitational-wave observatories during their second (O2) observing run. Using the Owens Valley Radio Observatory Long Wavelength Array (OVRO-LWA), we search a $\sim900\,\mathrm{deg}^2$ region for transient radio emission within approximately one hour of GW170104, obtaining an upper limit of $2.5\times10^{41}\,\mathrm{erg}\,\mathrm{s}^{-1}$ on its equivalent isotropic luminosity between 27-84 MHz. We additionally discuss plans to target binary neutron star mergers in Advanced LIGO and Virgo's upcoming O3 observing run.

The Space Technology 7 Disturbance Reduction System (ST7-DRS) is a NASA technology demonstration payload that operated from January 2016 through July of 2017 on the European Space Agency's LISA Pathfinder spacecraft. The joint goal of the NASA and ESA missions was to validate key technologies for a future space-based gravitational wave observatory targeting the source-rich milliHertz band. The two primary components of ST7-DRS are a micropropulsion system based on colloidal micro-Newton thrusters (CMNTs) and a control system that simultaneously controls the attitude and position of the spacecraft and the two free-flying test masses (TMs). This paper presents our main experimental results and summarizes the overall the performance of the CMNTs and control laws. We find that the CMNT performance to be consistent with pre-flight predictions, with a measured system thrust noise on the order of $100\,\textrm{nN}/\sqrt{\textrm{Hz}}$ in the $1\,\textrm{mHz}\leq f \leq 30\,\textrm{mHz}$ band. The control system maintained the TM-spacecraft separation with an RMS error of less than 2$\,$nm and a noise spectral density of less than $3\,\textrm{nm}/\sqrt{\textrm{Hz}}$ in the same band. Thruster calibration measurements yield thrust values consistent with the performance model and ground-based thrust-stand measurements, to within a few percent. We also report a differential acceleration noise between the two test masses with a spectral density of roughly $3\,\textrm{fm}/\textrm{s}^2/\sqrt{\textrm{Hz}}$ in the $1\,\textrm{mHz}\leq f \leq 30\,\textrm{mHz}$ band, slightly less than twice as large as the best performance reported with the baseline LISA Pathfinder configuration and below the current requirements for the Laser Interferometer Space Antenna (LISA) mission.

We present new HST-COS G130M spectroscopy for a sightline toward a filament projected 1.9 kpc from the nucleus of M87. The combination of the sensitivity of COS and the proximity of M87 allows us to study the structure of this filament in unparalleled detail. We propose that the filament is composed of many cold clumps, each surrounded by an FUV-emitting boundary layer, with the filament having a radius $r_c \sim 10$ pc and the clumps filling the cylinder with a low volume filling factor. The observed velocity dispersion in emission lines from the filament results from the random motions of these clumps within the filament. We measure fluxes and kinematics for emission lines of Ly$\alpha$, C II $\lambda$1335, and N V $\lambda1238$. We associate these three lines, as well as archival measurements of H$\alpha$, C IV $\lambda$1549, and He II $\lambda$1640, with a multitemperature boundary layer around clumps which are moving with supersonic random motions in the filament. This boundary layer is a significant coolant of the hot gas. We show that the [C II] $\lambda$158$\mu$m flux observed by Herschel-PACS from this region implies the existence of a massive cold ($T \sim 10^3$ K) component in the filament which contains significantly more mass than the FUV-emitting boundary layer. We also detect [Fe XXI] $\lambda$1354 in emission at $4-5\sigma$. This line is emitted from 1 keV ($T \approx 10^7$ K) plasma, and we use it to measure the bulk radial velocity and velocity dispersion of the plasma at this temperature. In contrast to the intermediate-temperature FUV lines, [\ion{Fe}{xxi}] is blueshifted relative to M87 and matches the bulk velocity of a nearby filament to the south. We hypothesize that this line arises from the approaching face of the radio bubble expanding through this sightline, while the filament lies on the receding side of the bubble. (abstract is abridged)

The Large-Aperture Experiment to Detect the Dark Age (LEDA) was designed to detect the predicted O(100)mK sky-averaged absorption of the Cosmic Microwave Background by Hydrogen in the neutral pre- and intergalactic medium just after the cosmological Dark Age. The spectral signature would be associated with emergence of a diffuse Ly$\alpha$ background from starlight during 'Cosmic Dawn'. Recently, Bowman et al. (2018) have reported detection of this predicted absorption feature, with an unexpectedly large amplitude of 530 mK, centered at 78 MHz. Verification of this result by an independent experiment, such as LEDA, is pressing. In this paper, we detail design and characterization of the LEDA radiometer systems, and a first-generation pipeline that instantiates a signal path model. Sited at the Owens Valley Radio Observatory Long Wavelength Array, LEDA systems include the station correlator, five well-separated redundant dual polarization radiometers and backend electronics. The radiometers deliver a 30-85MHz band (16<z<34) and operate as part of the larger interferometric array, for purposes ultimately of in situ calibration. Here, we report on the LEDA system design, calibration approach, and progress in characterization as of January 2016. The LEDA systems are currently being modified to improve performance near 78 MHz in order to verify the purported absorption feature.

We summarize and reanalyze observations bearing upon missing galactic baryons, where we propose a consistent picture for halo gas in L >~ L* galaxies. The hot X-ray emitting halos are detected to 50-70 kpc, where typically, M_hot(<50 kpc) ~ 5E9 Msun, and with density n \propto r^-3/2. When extrapolated to R200, the gas mass is comparable to the stellar mass, but about half of the baryons are still missing from the hot phase. If extrapolated to 1.9-3 R200, the baryon to dark matter ratio approaches the cosmic value. Significantly flatter density profiles are unlikely for R < 50 kpc and they are disfavored but not ruled out for R > 50 kpc. For the Milky Way, the hot halo metallicity lies in the range 0.3-1 solar for R < 50 kpc. Planck measurements of the thermal Sunyaev-Zeldovich effect toward stacked luminous galaxies (primarily early-type) indicate that most of their baryons are hot, near the virial temperature, and extend beyond R200. This stacked SZ signal is nearly an order of magnitude larger than that inferred from the X-ray observations of individual (mostly spiral) galaxies with M_* > 10^11.3 Msun. This difference suggests that the hot halo properties are distinct for early and late type galaxies, possibly due to different evolutionary histories. For the cooler gas detected in UV absorption line studies, we argue that there are two absorption populations: extended halos; and disks extending to ~50 kpc, containing most of this gas, and with masses a few times lower than the stellar masses. Such extended disks are also seen in 21 cm HI observations and in simulations.

The baryon content around local galaxies is observed to be much less than is needed in Big Bang nucleosynthesis. Simulations indicate that a significant fraction of these "missing baryons" may be stored in a hot tenuous circum-galactic medium (CGM) around massive galaxies extending to or even beyond the virial radius of their dark matter halos. Previous observations in X-ray and Sunyaev-Zel'dovich (SZ) signal claimed that $\sim(1-50)\%$ of the expected baryons are stored in a hot CGM within the virial radius. The large scatter is mainly caused by the very uncertain extrapolation of the hot gas density profile based on the detection in a small radial range (typically within 10\%-20\% of the virial radius). Here we report stacking X-ray observations of six local isolated massive spiral galaxies from the CGM-MASS sample. We find that the mean density profile can be characterized by a single power law out to a galactocentric radius of $\approx 200\rm~kpc$ (or $\approx130\rm~kpc$ above the 1~$\sigma$ background uncertainty), about half the virial radius of the dark matter halo. We can now estimate that the hot CGM within the virial radius accounts for $(8\pm4)\%$ of the baryonic mass expected for the halos. Including the stars, the baryon fraction is $(27\pm16)\%$, or $(39\pm20)\%$ by assuming a flattened density profile at $r\gtrsim130\rm~kpc$. We conclude that the hot baryons within the virial radius of massive galaxy halos are insufficient to explain the "missing baryons".

One of the essential prerequisites for detection of Earth-like extra-solar planets or direct measurements of the cosmological expansion is the accurate and precise wavelength calibration of astronomical spectrometers. It has already been realized that the large number of exactly known optical frequencies provided by laser frequency combs ('astrocombs') can significantly surpass conventionally used hollow-cathode lamps as calibration light sources. A remaining challenge, however, is generation of frequency combs with lines resolvable by astronomical spectrometers. Here we demonstrate an astrocomb generated via soliton formation in an on-chip microphotonic resonator ('microresonator') with a resolvable line spacing of 23.7 GHz. This comb is providing wavelength calibration on the 10 cm/s radial velocity level on the GIANO-B high-resolution near-infrared spectrometer. As such, microresonator frequency combs have the potential of providing broadband wavelength calibration for the next-generation of astronomical instruments in planet-hunting and cosmological research.

Double-lobed radio galaxies (DLRGs) often have radio lobes which subtend an angle of less than 180 degrees, and these bent DLRGs have been shown to associate preferentially with galaxy clusters and groups. In this study, we utilize a catalog of DLRGs in SDSS quasars with radio lobes visible in VLA FIRST 20 cm radio data. We cross-match this catalog against three catalogs of galaxies over the redshift range $0 < z < 0.70$, obtaining 81 tentative matches. We visually examine each match and apply a number of selection criteria, eventually obtaining a sample of 44 securely detected DLRGs which are paired to a nearby massive galaxy, galaxy group, or galaxy cluster. Most of the DLRGs identified in this manner are not central galaxies in the systems to which they are matched. Using this sample, we quantify the projected density of these matches as a function of projected separation from the central galaxy, finding a very steep decrease in matches as the impact parameter increases (for $\Sigma \propto b^{-m}$ we find $m = 2.5^{+0.4}_{-0.3}$) out to $b \sim 2$ Mpc. In addition, we show that the fraction of DLRGs with bent lobes also decreases with radius, so that if we exclude DLRGs associated with the central galaxy in the system the bent fraction is 78\% within 1 Mpc and 56\% within 2 Mpc, compared to just 29\% in the field; these differences are significant at $3.6\sigma$ and $2.8\sigma$ respectively. This behavior is consistent with ram pressure being the mechanism that causes the lobes to bend.

We have conducted the most sensitive low frequency (below 100 MHz) search to date for prompt, low-frequency radio emission associated with short-duration gamma-ray bursts (GRBs), using the Owens Valley Radio Observatory Long Wavelength Array (OVRO-LWA). The OVRO-LWA's nearly full-hemisphere field-of-view ($\sim20$,$000$ square degrees) allows us to search for low-frequency (sub-$100$ MHz) counterparts for a large sample of the subset of GRB events for which prompt radio emission has been predicted. Following the detection of short GRB 170112A by Swift, we used all-sky OVRO-LWA images spanning one hour prior to and two hours following the GRB event to search for a transient source coincident with the position of GRB 170112A. We detect no transient source, with our most constraining $1\sigma$ flux density limit of $650~\text{mJy}$ for frequencies spanning $27~\text{MHz}-84~\text{MHz}$. We place constraints on a number of models predicting prompt, low-frequency radio emission accompanying short GRBs and their potential binary neutron star merger progenitors, and place an upper limit of $L_\text{radio}/L_\gamma \lesssim 7\times10^{-16}$ on the fraction of energy released in the prompt radio emission. These observations serve as a pilot effort for a program targeting a wider sample of both short and long GRBs with the OVRO-LWA, including bursts with confirmed redshift measurements which are critical to placing the most constraining limits on prompt radio emission models, as well as a program for the follow-up of gravitational wave compact binary coalescence events detected by advanced LIGO and Virgo.

A host of new low-frequency radio telescopes seek to measure the 21-cm transition of neutral hydrogen from the early universe. These telescopes have the potential to directly probe star and galaxy formation at redshifts $20 \gtrsim z \gtrsim 7$, but are limited by the dynamic range they can achieve against foreground sources of low-frequency radio emission. Consequently, there is a growing demand for modern, high-fidelity maps of the sky at frequencies below 200 MHz for use in foreground modeling and removal. We describe a new widefield imaging technique for drift-scanning interferometers, Tikhonov-regularized $m$-mode analysis imaging. This technique constructs images of the entire sky in a single synthesis imaging step with exact treatment of widefield effects. We describe how the CLEAN algorithm can be adapted to deconvolve maps generated by $m$-mode analysis imaging. We demonstrate Tikhonov-regularized $m$-mode analysis imaging using the Owens Valley Long Wavelength Array (OVRO-LWA) by generating 8 new maps of the sky north of $\delta=-30^\circ$ with 15 arcmin angular resolution, at frequencies evenly spaced between 36.528 MHz and 73.152 MHz, and $\sim$800 mJy/beam thermal noise. These maps are a 10-fold improvement in angular resolution over existing full-sky maps at comparable frequencies, which have angular resolutions $\ge 2^\circ$. Each map is constructed exclusively from interferometric observations and does not represent the globally averaged sky brightness. Future improvements will incorporate total power radiometry, improved thermal noise, and improved angular resolution -- due to the planned expansion of the OVRO-LWA to 2.6 km baselines. These maps serve as a first step on the path to the use of more sophisticated foreground filters in 21-cm cosmology incorporating the measured angular and frequency structure of all foreground contaminants.

We present the analysis of the XMM-Newton data of the Circum-Galactic Medium of MASsive Spirals (CGM-MASS) sample of six extremely massive spiral galaxies in the local Universe. All the CGM-MASS galaxies have diffuse X-ray emission from hot gas detected above the background extending $\sim(30-100)\rm~kpc$ from the galactic center. This doubles the existing detection of such extended hot CGM around massive spiral galaxies. The radial soft X-ray intensity profile of hot gas can be fitted with a $\beta$-function with the slope typically in the range of $\beta=0.35-0.55$. This range, as well as those $\beta$ values measured for other massive spiral galaxies, including the Milky Way (MW), are in general consistent with X-ray luminous elliptical galaxies of similar hot gas luminosity and temperature, and with those predicted from a hydrostatic isothermal gaseous halo. Hot gas in such massive spiral galaxy tends to have temperature comparable to its virial value, indicating the importance of gravitational heating. This is in contrast to lower mass galaxies where hot gas temperature tends to be systematically higher than the virial one. The ratio of the radiative cooling to free fall timescales of hot gas is much larger than the critical value of $\sim10$ throughout the entire halos of all the CGM-MASS galaxies, indicating the inefficiency of gas cooling and precipitation in the CGM. The hot CGM in these massive spiral galaxies is thus most likely in a hydrostatic state, with the feedback material mixed with the CGM, instead of escaping out of the halo or falling back to the disk. We also homogenize and compare the halo X-ray luminosity measured for the CGM-MASS galaxies and other galaxy samples and discuss the "missing" galactic feedback detected in these massive spiral galaxies.

During the Space Telescope and Optical Reverberation Mapping Project (STORM) observations of NGC 5548, the continuum and emission-line variability became de-correlated during the second half of the 6-month long observing campaign. Here we present Swift and Chandra X-ray spectra of NGC 5548 obtained as a part of the campaign. The Swift spectra show that excess flux (relative to a power-law continuum) in the soft X-ray band appears before the start of the anomalous emission-line behavior, peaks during the period of the anomaly, and then declines. This is a model-independent result suggesting that the soft excess is related to the anomaly. We divide the Swift data into on- and off-anomaly spectra to characterize the soft excess via spectral fitting. The cause of the spectral differences is likely due to a change in the intrinsic spectrum rather than being due to variable obscuration or partial covering. The Chandra spectra have lower signal-to-noise ratios, but are consistent with Swift data. Our preferred model of the soft excess is emission from an optically thick, warm Comptonizing corona, the effective optical depth of which increases during the anomaly. This model simultaneously explains all the three observations: the UV emission line flux decrease, the soft-excess increase, and the emission line anomaly.

Supermassive binary black holes (BBHs) on sub-parsec scales are prime targets for gravitational wave experiments. They also provide insights on close binary evolution and hierarchical structure formation. Sub-parsec BBHs cannot be spatially resolved but indirect methods can identify candidates. In 2015 Liu et al. reported an optical-continuum periodicity in the quasar PSO J334.2028+01.4075, with the estimated mass and rest-frame period suggesting an orbital separation of about 0.006 pc (0.7 micro-arcsec). The persistence of the quasar's optical periodicity has recently been disfavored over an extended baseline. However, if a radio jet is launched from a sub-parsec BBH, the binary's properties can influence the radio structure on larger scales. Here, we use the Very Long Baseline Array (VLBA) and Karl G. Jansky Very Large Array (VLA) to study the parsec- and kiloparsec-scale emission energized by the quasar's putative BBH. We find two VLBA components separated by 3.6 mas (30 pc), tentatively identifying one as the VLBA "core" from which the other was ejected. The VLBA components contribute to a point-like, time-variable VLA source that is straddled by lobes spanning 8 arcsec (66 kpc). We classify PSO J334.2028+01.4075 as a lobe-dominated quasar, albeit with an atypically large twist of 39 deg between its elongation position angles on parsec- and kiloparsec-scales. By analogy with 3C 207, a well-studied lobe-dominated quasar with a similarly-rare twist, we speculate that PSO J334.2028+01.4075 could be ejecting jet components over an inner cone that traces a precessing jet in a BBH system.

We present the results of an optical spectroscopic monitoring program targeting NGC 5548 as part of a larger multi-wavelength reverberation mapping campaign. The campaign spanned six months and achieved an almost daily cadence with observations from five ground-based telescopes. The H$\beta$ and He II $\lambda$4686 broad emission-line light curves lag that of the 5100 $\AA$ optical continuum by $4.17^{+0.36}_{-0.36}$ days and $0.79^{+0.35}_{-0.34}$ days, respectively. The H$\beta$ lag relative to the 1158 $\AA$ ultraviolet continuum light curve measured by the Hubble Space Telescope is roughly $\sim$50% longer than that measured against the optical continuum, and the lag difference is consistent with the observed lag between the optical and ultraviolet continua. This suggests that the characteristic radius of the broad-line region is $\sim$50% larger than the value inferred from optical data alone. We also measured velocity-resolved emission-line lags for H$\beta$ and found a complex velocity-lag structure with shorter lags in the line wings, indicative of a broad-line region dominated by Keplerian motion. The responses of both the H$\beta$ and He II $\lambda$4686 emission lines to the driving continuum changed significantly halfway through the campaign, a phenomenon also observed for C IV, Ly $\alpha$, He II(+O III]), and Si IV(+O IV]) during the same monitoring period. Finally, given the optical luminosity of NGC 5548 during our campaign, the measured H$\beta$ lag is a factor of five shorter than the expected value implied by the $R_\mathrm{BLR} - L_\mathrm{AGN}$ relation based on the past behavior of NGC 5548.

We conduct a multiwavelength continuum variability study of the Seyfert 1 galaxy NGC 5548 to investigate the temperature structure of its accretion disk. The 19 overlapping continuum light curves (1158 to 9157 angstroms) combine simultaneous HST , Swift , and ground-based observations over a 180 day period from 2014 January to July. Light-curve variability is interpreted as the reverberation response of the accretion disk to irradiation by a central time-varying point source. Our model yields the disk inclination, i, temperature T1 at 1 light day from the black hole, and a temperature-radius slope, alpha. We also infer the driving light curve and find that it correlates poorly with both the hard and soft X-ray light curves, suggesting that the X-rays alone may not drive the ultraviolet and optical variability over the observing period. We also decompose the light curves into bright, faint, and mean accretion-disk spectra. These spectra lie below that expected for a standard blackbody accretion disk accreting at L/LEdd = 0.1

We present the results of a recent reverberation mapping campaign for UGC 06728, a nearby low-luminosity Seyfert 1 in a late-type galaxy. Nightly monitoring in the spring of 2015 allowed us to determine an H$\beta$ time delay of $\tau = 1.4 \pm 0.8$ days. Combined with the width of the variable H$\beta$ line profile, we determine a black hole mass of $M_{\rm BH} = (7.1 \pm 4.0) \times 10^5$ M$_{\odot}$. We also constrain the bulge stellar velocity dispersion from higher-resolution long slit spectroscopy along the galaxy minor axis and find $\sigma_{\star} = 51.6 \pm 4.9$ km s$^{-1}$. The measurements presented here are in good agreement with both the $R_{\rm BLR} - L$ relationship and the $M_{\rm BH}-\sigma_{\star}$ relationship for AGNs. Combined with a previously published spin measurement, our mass determination for UGC 06728 makes it the lowest-mass black hole that has been fully characterized, and thus an important object to help anchor the low-mass end of black hole evolutionary models.

The Circum-Galactic Medium of MASsive Spirals (CGM-MASS) is a project studying the overall content, physical and chemical properties, and spatial distributions of the multi-phase circum-galactic medium (CGM) around a small sample of the most massive isolated spiral galaxies in the local Universe. We introduce the sample and present a detailed case study of the XMM-Newton observation of the hot gas halo of NGC5908. After data calibration, point source removal, and background analysis, we find that the diffuse soft X-ray emission of NGC5908 is significantly more extended than the stellar light in the vertical direction. The 0.5-1.25keV radial intensity profile tracing hot gas emission can be detected above the background out to about $2^\prime$, or $30\rm~kpc$ from the nucleus. The unresolved soft X-ray emission can be characterized with a $\beta$-model with a slope of $\beta\approx0.68$. The unresolved 0.5-2keV luminosity within $1^\prime$ is $6.8\times10^{39}\rm~ergs~s^{-1}$, but young stellar sources may contribute significantly to it. Assuming a metallicity of 0.2solar, an upper limit (without subtracting the very uncertain young stellar contribution) to the mass of hot gas within this radius is $2.3\times10^9\rm~M_\odot$. The cooling radius is $r_{\rm cool}\approx25\rm~kpc$ or $\approx0.06r_{\rm 200}$, within which the hot gas could cool radiatively in less than 10Gyr, and the cooling of hot gas could significantly contribute in replenishing the gas consumed in star formation. The hot gas accounts for $\approx1.9\%$ of the baryon detected within the cooling radius. By comparing NGC5908 to other galaxies, we find that its X-ray luminosity per stellar mass is consistent with lower-mass non-starburst field spiral galaxies. However, a large scatter in hot gas soft X-ray emissivity is indicated for spiral galaxies with $M_*\gtrsim2\times10^{11}\rm~M_\odot$.

Non-thermal pressure from turbulence and bulk flows is a fundamental ingredient in hot gaseous halos, and in the intracluster medium it will be measured through emission line kinematics with calorimeters on future X-ray spacecraft. In this paper we present a complementary method for measuring these effects, using forbidden FUV emission lines of highly ionized Iron which trace $10^7$ K gas. The brightest of these is [Fe XXI] $\lambda$1354.1. We search for these lines in archival HST-COS spectra from the well-known elliptical galaxies M87 and NGC4696, which harbor large reservoirs of $10^7$ K gas. We report a 2.2$\sigma$ feature which we attribute to [Fe XXI] from a filament in M87, and positive residuals in the nuclei of M87 and NGC4696, for which the 90\% upper limits on the line flux are close to the predicted fluxes based on X-ray observations. In a newer reduction of the data from the Hubble Spectroscopic Legacy Archive, these limits become tighter and the [Fe XXI] feature reaches a formal significance of 5.3$\sigma$, neglecting uncertainty in fitting the continuum. Using our constraints, we perform emission measure analysis, constraining the characteristic path length and column density of the $\sim10^7$ K gas. We also examine several sightlines towards filaments or cooling flows in other galaxy clusters, for which the fraction of gas at $10^7$ K is unknown, and place upper limits on its emission measure in each case. A medium-resolution HST-COS observation of the M87 filament for $\sim$10 orbits would confirm our detection of [Fe XXI] and measure its width.

During an intensive Hubble Space Telescope (HST) Cosmic Origins Spectrograph (COS) UV monitoring campaign of the Seyfert~1 galaxy NGC 5548 performed from 2014 February to July, the normally highly correlated far-UV continuum and broad emission-line variations decorrelated for ~60 to 70 days, starting ~75 days after the first HST/COS observation. Following this anomalous state, the flux and variability of the broad emission lines returned to a more normal state. This transient behavior, characterised by significant deficits in flux and equivalent width of the strong broad UV emission lines, is the first of its kind to be unambiguously identified in an active galactic nucleus reverberation mapping campaign. The largest corresponding emission-line flux deficits occurred for the high-ionization collisionally excited lines, C IV and Si IV(+O IV]), and also He II(+O III]), while the anomaly in Ly-alpha was substantially smaller. This pattern of behavior indicates a depletion in the flux of photons with E_{\rm ph} > 54 eV, relative to those near 13.6 eV. We suggest two plausible mechanisms for the observed behavior: (i) temporary obscuration of the ionizing continuum incident upon BLR clouds by a moving veil of material lying between the inner accretion disk and inner BLR, perhaps resulting from an episodic ejection of material from the disk, or (ii) a temporary change in the intrinsic ionizing continuum spectral energy distribution resulting in a deficit of ionizing photons with energies > 54 eV, possibly due to a transient restructuring of the Comptonizing atmosphere above the disk. Current evidence appears to favor the latter explanation.

We study the merger of binary neutron stars with different mass ratios adopting three different realistic, microphysical nuclear equations of state, as well as incorporating neutrino cooling effects. In particular, we concentrate on the influence of the equation of state on the gravitational wave signature and also on its role, in combination with neutrino cooling, in determining the properties of the resulting hypermassive neutron star, of the neutrinos produced, and of the ejected material. The ejecta we find are consistent with other recent studies that find that small mass ratios produce more ejecta than equal mass cases (up to some limit) and this ejecta is more neutron rich. This trend indicates the importance with future kilonovae observations of measuring the individual masses of an associated binary neutron star system, presumably from concurrent gravitational wave observations, in order to be able to extract information about the nuclear equation of state

We present ground-based optical photometric monitoring data for NGC 5548, part of an extended multi-wavelength reverberation mapping campaign. The light curves have nearly daily cadence from 2014 January to July in nine filters (\emph{BVRI} and \emph{ugriz}). Combined with ultraviolet data from the \emph{Hubble Space Telescope} and \emph{Swift}, we confirm significant time delays between the continuum bands as a function of wavelength, extending the wavelength coverage from 1158\,\AA\ to the $z$ band ($\sim\!9160$\,\AA). We find that the lags at wavelengths longer than the {\it V} band are equal to or greater than the lags of high-ionization-state emission lines (such as He\,{\sc ii}\,$\lambda 1640$ and $\lambda 4686$), suggesting that the continuum-emitting source is of a physical size comparable to the inner broad-line region (BLR). The trend of lag with wavelength is broadly consistent with the prediction for continuum reprocessing by an accretion disk with $\tau \propto \lambda^{4/3}$. However, the lags also imply a disk radius that is 3 times larger than the prediction from standard thin-disk theory, assuming that the bolometric luminosity is 10\% of the Eddington luminosity ($L = 0.1L_{\rm Edd}$). Using optical spectra from the Large Binocular Telescope, we estimate the bias of the interband continuum lags due to BLR emission observed in the filters. We find that the bias for filters with high levels of BLR contamination ($\sim\! 20\%$) can be important for the shortest continuum lags, and likely has a significant impact on the {\it u} and {\it U} bands owing to Balmer continuum emission.

Methods to solve the relativistic hydrodynamic equations are a key computational kernel in a large number of astrophysics simulations and are crucial to understanding the electromagnetic signals that originate from the merger of astrophysical compact objects. Because of the many physical length scales present when simulating such mergers, these methods must be highly adaptive and capable of automatically resolving numerous localized features and instabilities that emerge throughout the computational domain across many temporal scales. While this has been historically accomplished with adaptive mesh refinement (AMR) based methods, alternatives based on wavelet bases and the wavelet transformation have recently achieved significant success in adaptive representation for advanced engineering applications. This work presents a new method for the integration of the relativistic hydrodynamic equations using iterated interpolating wavelets and introduces a highly adaptive implementation for multidimensional simulation. The wavelet coefficients provide a direct measure of the local approximation error for the solution and place collocation points that naturally adapt to the fluid flow while providing good conservation of fluid quantities. The resulting implementation, OAHU, is applied to a series of demanding one- and two-dimensional problems which explore high Lorentz factor outflows and the formation of several instabilities, including the Kelvin-Helmholtz instability and the Rayleigh-Taylor instability.

We use weak gravitational lensing to measure mean mass profiles around Locally Brightest Galaxies (LBGs). These are selected from the SDSS/DR7 spectroscopic and photometric catalogues to be brighter than any neighbour projected within 1.0 Mpc and differing in redshift by $<1000$ km/s. Most ($> 83\%$) are expected to be the central galaxies of their dark matter halos. Previous stacking analyses have used this LBG sample to measure mean Sunyaev-Zeldovich flux and mean X-ray luminosity as a function of LBG stellar mass. In both cases, a simulation of the formation of the galaxy population was used to estimate effective halo mass for LBGs of given stellar mass, allowing the derivation of scaling relations between the gas properties of halos and their mass. By comparing results from a variety of simulations to our lensing data, we show that this procedure has significant model dependence reflecting: (i) the failure of any given simulation to reproduce observed galaxy abundances exactly; (ii) a dependence on the cosmology underlying the simulation; and (iii) a dependence on the details of how galaxies populate halos. We use our lensing results to recalibrate the scaling relations, eliminating most of this model dependence and explicitly accounting both for residual modelling uncertainties and for observational uncertainties in the lensing results. The resulting scaling relations link the mean gas properties of dark halos to their mass over an unprecedentedly wide range, $10^{12.5}<M_{500}/ \mathrm{M_\odot}<10^{14.5}$, and should fairly and robustly represent the full halo population.

We examine 11 XMM-Newton observations of the giant spiral galaxy NGC 1961, allowing us to study the hot gaseous halo of a spiral galaxy in unprecedented detail. We perform a spatial and a spectral analysis; with the former, the hot halo is detected to at least 80 kpc and with the latter its properties can be measured in detail up to 42 kpc. We find evidence for a negative gradient in the temperature profile as is common for elliptical galaxies. We measure a rough metallicity profile, which is consistent with being flat at $Z \sim 0.2 Z_{\odot}$. Converting to this metallicity, the deprojected density profile is consistent with previous parametric fits, with no evidence for a break within 42 kpc ($\sim$0.1R$_{\text{vir}}$). Extrapolating to the virial radius, we infer a hot halo mass comparable to the stellar mass of the galaxy, and a baryon fraction from the stars and hot gas of around 30%. The cooling time of the hot gas is orders of magnitude longer than the dynamical time, making the hot halo stable against cooling instabilities, and we argue that an extended stream of neutral Hydrogen seen to the NW of this galaxy is instead likely due to accretion from the intergalactic medium. The low metallicity of the hot halo suggests it too was likely accreted. We compare the hot halo of NGC 1961 to hot halos around isolated elliptical galaxies, and show that the total mass determines the hot halo properties better than the stellar mass.

We conduct a comprehensive search for X-ray emission lines from sterile neutrino dark matter, motivated by recent claims of unidentified emission lines in the stacked X-ray spectra of galaxy clusters and the centers of the Milky Way and M31. Since the claimed emission lines lie around 3.5 keV, we focus on galaxies and galaxy groups (masking the central regions), since these objects emit very little radiation above $\sim 2$ keV and offer a clean background against which to detect emission lines. We develop a formalism for maximizing the signal-to-noise of decaying dark matter emission lines by weighing each X-ray event according to the expected dark matter profile. In total, we examine 81 and 89 galaxies with Chandra and XMM-Newton respectively, totaling 15.0 and 14.6 Ms of integration time. We find no significant evidence of any emission lines, placing strong constraints on the mixing angle of sterile neutrinos with masses between 4.8-12.4 keV. In particular, if the 3.57 keV feature from Bulbul et al. (2014) were due to 7.1 keV sterile neutrino emission, we would have detected it at $4.4\sigma$ and $11.8\sigma$ in our two samples. The most conservative estimates of the systematic uncertainties reduce these constraints to $4.4\sigma$ and 7.8$\sigma$, or letting the line energy vary between 3.50 and 3.60 keV reduces these constraints to $2.7\sigma$ and $11.0\sigma$ respectively. Unlike previous constraints, our measurements do not depend on the model of the X-ray background or on the assumed logarithmic slope of the center of the dark matter profile.

We study the merger of binary neutron stars using different realistic, microphysical nuclear equations of state, as well as incorporating magnetic field and neutrino cooling effects. In particular, we concentrate on the influence of the equation of state on the gravitational wave signature and also on its role, in combination with cooling and electromagnetic effects, in determining the properties of the hypermassive neutron star resulting from the merger, the production of neutrinos, and the characteristics of ejecta from the system. The ejecta we find are consistent with other recent studies that find soft equations of state produce more ejecta than stiffer equations of state. Moreover, the degree of neutron richness increases for softer equations of state. In light of reported kilonova observations (associated to GRB~130603B and GRB~060614) and the discovery of relatively low abundances of heavy, radioactive elements in deep sea deposits (with respect to possible production via supernovae), we speculate that a soft EoS might be preferred---because of its significant production of sufficiently neutron rich ejecta---if such events are driven by binary neutron star mergers. We also find that realistic magnetic field strengths, obtained with a sub-grid model tuned to capture magnetic amplification via the Kelvin-Helmholtz instability at merger, are generally too weak to affect the gravitational wave signature post-merger within a time scale of $\approx 10$~ms but can have subtle effects on the post-merger dynamics.

We examine a sample of $\sim 250 000$ "locally brightest galaxies" selected from the Sloan Digital Sky Survey to be central galaxies within their dark matter halos. These were previously stacked by the Planck Collaboration to measure the Sunyaev-Zel'dovich signal as a function of central galaxy stellar mass. Here, we stack the X-ray emission from these halos using data from the ROSAT All-Sky Survey. We detect emission across almost our entire sample, including emission which we attribute to hot gas around galaxies spanning a range of 1.2 dex in stellar mass (corresponding to nearly two orders of magnitude in halo mass) down to $M* = 10^{10.8} M_{\odot}$ ($M_{500} \approx 10^{12.6} M_{\odot}$). Over this range, the X-ray luminosity can be fit by a power-law, either of stellar mass or of halo mass. A single unified scaling relation between mass and $L_X$ applies for galaxies, groups, and clusters. This relation has a steeper slope than expected for self-similarity, in contrast to the $Y_{SZ}$-$M_{500}$ relation, showing the importance of non-gravitational heating. If this heating is predominantly due to AGN feedback, the lack of a break in our relation suggests that AGN feedback is tightly self-regulated and fairly gentle, in agreement with recent high-resolution simulations. Our results are consistent with earlier measurements of the $L_X$-$L_K$ relation for elliptical galaxies and of the $L_X$-$M_{500}$ relation for optically-selected galaxy clusters. However, our $L_X$-$M_{500}$ relation lies more than a factor of two below most previous relations based on X-ray-selected cluster samples. We argue that optical selection offers a less biased view of the $L_X$-$M_{500}$ relation.

We present the results of a reverberation-mapping program targeting NGC5273, a nearby early-type galaxy with a broad-lined active galactic nucleus. Over the course of the monitoring program, NGC5273 showed strong variability that allowed us to measure time delays in the responses of the broad optical recombination lines to changes in the continuum flux. A weighted average of these measurements results in a black hole mass determination of $M_{\rm BH} = (4.7 \pm 1.6) \times 10^6$ M$_{\odot}$. An estimate of the size of the black hole sphere of influence in NGC5273 puts it just at the limit of the resolution achievable with current ground-based large aperture telescopes. NGC5273 is therefore an important future target for a black hole mass determination from stellar dynamical modeling, especially because it is the only nearby early-type galaxy hosting an AGN with a reverberation-based mass, allowing the best comparison for the masses determined from these two techniques.

Two red supergiants of the Per OB1 association, RS Per and T Per, have been observed in H band using the MIRC instrument at the CHARA array. The data show clear evidence of departure from circular symmetry. We present here new techniques specially developed to analyze such cases, based on state-of-the-art statistical frameworks. The stellar surfaces are first modeled as limb-darkened discs based on SATLAS models that fit both MIRC interferometric data and publicly available spectrophotometric data. Bayesian model selection is then used to determine the most probable number of spots. The effective surface temperatures are also determined and give further support to the recently derived hotter temperature scales of red su- pergiants. The stellar surfaces are reconstructed by our model-independent imaging code SQUEEZE, making use of its novel regularizer based on Compressed Sensing theory. We find excellent agreement between the model-selection results and the reconstructions. Our results provide evidence for the presence of near-infrared spots representing about 3-5% of the stellar flux.

We incorporate realistic, tabulated equations of state into fully relativistic simulations of magnetized neutron stars along with a neutrino leakage scheme which accounts for cooling via neutrino emission. Both these improvements utilize open-source code (GR1D) and tables from http://stellarcollapse.org. Our implementation makes use of a novel method for the calculation of the optical depth which simplifies its use with distributed adaptive mesh refinement. We present various tests with and without magnetization and preliminary results both from single stars and from the merger of a binary system.

Extended X-ray emission can be studied either spatially (through its surface brightness profile) or spectrally (by analyzing the spectrum at various locations in the field). Both techniques have advantages and disadvantages, and when the emission becomes particularly faint and/or extended, the two methods can disagree. We argue that an ideal approach would be to model the events file directly, and therefore to use both the spectral and spatial information which are simultaneously available for each event. In this work we propose a first step in this direction, introducing a method for spatial analysis which can be extended to leverage spectral information simultaneously. We construct a model for the entire X-ray image in a given energy band, and generate a likelihood function to compare the model to the data. A critical goal of this modeling is disentangling vignetted and unvignetted backgrounds through their different spatial distributions. Employing either maximum likelihood or Markov Chain Monte Carlo, we can derive probability distribution functions for the source and background parameters together, or we can fit and subtract the background, leaving the description of the source non-parametric. We calibrate and demonstrate this method against a variety of simulated images, and then apply it to Chandra observations of the hot gaseous halo around the elliptical galaxy NGC 720. We are able to follow the X-ray emission below a tenth of the background, and to infer a hot gas mass within 35 kpc of 4-5x10^9 Msun, with some indication that the profile continues to at least 50 kpc and that it steepens as the radius increases. We derive much stronger constraints on the surface brightness profile than previous studies, which employed the spectral method, and we show that the density profiles inferred from these studies are in conflict with the observed surface brightness profile. (abridged)

We expand on our study of the gravitational and electromagnetic emissions from the late stage of an inspiraling neutron star binary as presented in Ref. \cite{Palenzuela:2013hu}. Interactions between the stellar magnetospheres, driven by the extreme dynamics of the merger, can yield considerable outflows. We study the gravitational and electromagnetic waves produced during the inspiral and merger of a binary neutron star system using a full relativistic, resistive MHD evolution code. We show that the interaction between the stellar magnetospheres extracts kinetic energy from the system and powers radiative Poynting flux and heat dissipation. These features depend strongly on the configuration of the initial stellar magnetic moments. Our results indicate that this power can strongly outshine pulsars in binaries and have a distinctive angular and time-dependent pattern. Our discussion provides more detail than Ref. \cite{Palenzuela:2013hu}, showing clear evidence of the different effects taking place during the inspiral. Our simulations include a few milliseconds after the actual merger and study the dynamics of the magnetic fields during the formation of the hypermassive neutron star. We also briefly discuss the possibility of observing such emissions.

The late stage of an inspiraling neutron star binary gives rise to strong gravitational wave emission due to its highly dynamic, strong gravity. Moreover, interactions between the stellar magnetospheres can produce considerable electromagnetic radiation. We study this scenario using fully general relativistic, resistive magneto-hydrodynamics simulations. We show that these interactions extract kinetic energy from the system, dissipate heat, and power radiative Poynting flux, as well as develop current sheets. Our results indicate that this power can: (i) outshine pulsars in binaries, (ii) display a distinctive angular- and time-dependent pattern, and (iii) radiate within large opening angles. These properties suggest that some binary neutron star mergers are ideal candidates for multimessenger astronomy.

We place general constraints on the luminosity and mass of hot X-ray emitting gas residing in extended "hot halos" around nearby massive galaxies. We examine stacked images of 2165 galaxies from the 2MASS Very Isolated Galaxy Catalog (2MVIG), as well as subsets of this sample based on galaxy morphology and K-band luminosity. We detect X-ray emission at high confidence (ranging up to nearly 10\sigma) for each subsample of galaxies. The average L_X within 50 kpc is 1.0\pm0.1 (statistical) \pm0.2 (systematic) x10^40 erg/s, although the early-type galaxies are more than twice as luminous as the late-type galaxies. Using a spatial analysis, we also find evidence for extended emission around five out of seven subsamples (the full sample, the luminous galaxies, early-type galaxies, luminous late-type galaxies, and luminous early-type galaxies) at 92.7%, 99.3%, 89.3%, 98.7%, and 92.1% confidence, respectively. Several additional lines of evidence also support this conclusion and suggest that about 1/2 of the total emission is extended, and about 1/3 of the extended emission comes from hot gas. For the sample of luminous galaxies, which has the strongest evidence for extended emission, the average hot gas mass is 4x10^9 Msun within 50 kpc and the implied accretion rate is 0.4 Msun/yr.

While the importance of dusty asymptotic giant branch (AGB) stars to galactic chemical enrichment is widely recognised, a sophisticated understanding of the dust formation and wind-driving mechanisms has proven elusive due in part to the difficulty in spatially-resolving the dust formation regions themselves. We have observed twenty dust-enshrouded AGB stars as part of the Keck Aperture Masking Experiment, resolving all of them in multiple near-infrared bands between 1.5 microns and 3.1 microns. We find 45% of the targets to show measurable elongations that, when correcting for the greater distances of the targets, would correspond to significantly asymmetric dust shells on par with the well-known cases of IRC+10216 or CIT6. Using radiative transfer models, we find the sublimation temperature of 1130 +- 90 K and 1170 +- 60 K for silicates and amorphous carbon respectively, both somewhat lower than expected from laboratory measurements and vastly below temperatures inferred from the inner edge of YSO disks. The fact that O-rich and C-rich dust types showed the same sublimation temperature was surprising as well. For the most optically-thick shells (tau > 2 at 2.2 microns), the temperature profile of the inner dust shell is observed to change substantially, an effect we suggest could arise when individual dust clumps become optically-thick at the highest mass-loss rates.

We present our XMM-Newton observation of the fastest rotating spiral galaxy UGC 12591. We detect hot gas halo emission out to 110 kpc from the galaxy center, and constrain the halo gas mass to be smaller than 3.5e11 solar masses. We also measure the temperature of the hot gas as T=0.64\pm0.03 keV. Combining our X-ray constraints and the near-infrared and radio measurements in the literature, we find a baryon mass fraction of 0.03-0.04 in UGC 12591, suggesting a missing baryon mass of 75% compared with the cosmological mean value. Combined with another recent measurement in NGC 1961, the result strongly argues that the majority of missing baryons in spiral galaxies does not reside in their hot halos. We also find that UGC 12591 lies significantly below the baryonic Tully-Fisher relationship. Finally, we find that the baryon fractions of massive spiral galaxies are similar to those of galaxy groups with similar masses, indicating that the baryon loss is ultimately controlled by the gravitational potential well. The cooling radius of this gas halo is small, similar to NGC 1961, which argues that the majority of stellar mass of this galaxy is not assembled as a result of cooling of this gas halo.

The Allen Telescope Array was used to monitor Mars between 9 March and 2 June 2010, over a total of approximately 30 hours, for radio emission indicative of electrostatic discharge. The search was motivated by the report from Ruf et al. (2009) of the detection of non-thermal microwave radiation from Mars characterized by peaks in the power spectrum of the kurtosis, or kurtstrum, at 10 Hz, coinciding with a large dust storm event on 8 June 2006. For these observations, we developed a wideband signal processor at the Center for Astronomy Signal Processing and Electronics Research (CASPER). This 1024-channel spectrometer calculates the accumulated power and power-squared, from which the spectral kurtosis is calculated post-observation. Variations in the kurtosis are indicative of non-Gaussianity in the signal, which can be used to detect variable cosmic signals as well as radio frequency interference (RFI). During the three month period of observations, dust activity occurred on Mars in the form of small-scale dust storms, however no signals indicating lightning discharge were detected. Frequent signals in the kurtstrum that contain spectral peaks with an approximate 10 Hz fundamental were seen at both 3.2 and 8.0 GHz, but were the result of narrowband RFI with harmonics spread over a broad frequency range.

Hot gaseous halos are predicted around all large galaxies and are critically important for our understanding of galaxy formation, but they have never been detected at distances beyond a few kpc around a spiral galaxy. We used the Chandra ACIS-I instrument to search for diffuse X-ray emission around an ideal candidate galaxy: the isolated giant spiral NGC 1961. We observed four quadrants around the galaxy for 30 ks each, carefully subtracting background and point source emission, and found diffuse emission that appears to extend to 40-50 kpc. We fit $\beta$-models to the emission, and estimate a hot halo mass within 50 kpc of $5\times10^9 M_{\odot}$. When this profile is extrapolated to 500 kpc (the approximate virial radius), the implied hot halo mass is $1-3\times10^{11} M_{\odot}$. These mass estimates assume a gas metallicity of $Z = 0.5 Z_{\odot}$. This galaxy's hot halo is a large reservoir of gas, but falls significantly below observational upper limits set by pervious searches, and suggests that NGC 1961 is missing 75% of its baryons relative to the cosmic mean, which would tentatively place it below an extrapolation of the baryon Tully-Fisher relationship of less massive galaxies. The cooling rate of the gas is no more than 0.4 $M_{\odot}$/year, more than an order of magnitude below the gas consumption rate through star formation. We discuss the implications of this halo for galaxy formation models.

We report an experimental observation of an instability in gas of constant density (air) with an initial non-uniform seeding of small droplets that develops as a planar shock wave passes through the gas-droplet mix. The seeding non-uniformity is produced by vertical injection of a slow-moving jet of air pre-mixed with glycol droplets into the test section of a shock tube, with the plane of the shock parallel to the axis of the jet. After the shock passage, we observe development of two counter-rotating vortices in the plane normal to that axis. The physical mechanism of the instability we observe is peculiar to multiphase flow, where the shock acceleration causes the second (embedded) phase to move with respect to the embedding medium. With sufficient seeding concentration, this leads to entrainment of the embedding phase that acquires a relative velocity dependent on the initial seeding, resulting in vortex formation in the flow.

We investigate the merger of a neutron star (of compaction ratio $0.1$) in orbit about a spinning black hole in full general relativity with a mass ratio of $5:1$, allowing for the star to have an initial magnetization of $10^{12} {\rm Gauss}$. We present the resulting gravitational waveform and analyze the fallback accretion as the star is disrupted. The evolutions suggest no significant effects from the initial magnetization. We find that only a negligible amount of matter becomes unbound; $99%$ of the neutron star material has a fallback time of 10 seconds or shorter to reach the region of the central engine and that $99.99%$ of the star will interact with the central disk and black hole within 3 hours.

The metallicity in galaxy clusters is expected to originate from the stars in galaxies, with a population dominated by high mass stars likely being the most important stellar component, especially in rich clusters. We examine the relationship between the metallicity and the prominence of galaxies as measured by the star to baryon ratio, M$_*$/M$_{bary}$. Counter to expectations, we rule out a metallicity that is proportional to M$_*$/M$_{bary}$, where the best fit has the gas phase metallicity decreasing with M$_*$/M$_{bary}$, or the metallicity of the gas plus the stars being independent of M$_*$/M$_{bary}$. This implies that the population of stars responsible for the metals is largely proportional to the total baryonic mass of the cluster, not to the galaxy mass within the cluster. If generally applicable, most of the heavy elements in the universe were not produced within galaxies.

Light curves in the B, V, and I_c passbands have been obtained for the type II Cepheids V154 in M3 and V42 and V84 in M5. Alternating cycle behavior, similar to that seen among RV Tauri variables, is confirmed for V84. Old and new observations, spanning more than a century, show that V154 has increased in period while V42 has decreased in period. V84, on the other hand, has shown large, erratic changes in period that do not appear to reflect the long term evolution of V84 through the HR diagram.

Galaxies are missing most of their baryons, and many models predict these baryons lie in a hot halo around galaxies. We establish observationally motivated constraints on the mass and radii of these haloes using a variety of independent arguments. First, the observed dispersion measure of pulsars in the Large Magellanic Cloud allows us to constrain the hot halo around the Milky Way: if it obeys the standard NFW profile, it must contain less than 4-5% of the missing baryons from the Galaxy. This is similar to other upper limits on the Galactic hot halo, such as the soft X-ray background and the pressure around high velocity clouds. Second, we note that the X-ray surface brightness of hot haloes with NFW profiles around large isolated galaxies is high enough that such emission should be observed, unless their haloes contain less than 10-25% of their missing baryons. Third, we place constraints on the column density of hot haloes using nondetections of OVII absorption along AGN sightlines: in general they must contain less than 70% of the missing baryons or extend to no more than 40 kpc. Flattening the density profile of galactic hot haloes weakens the surface brightness constraint so that a typical L$_*$ galaxy may hold half its missing baryons in its halo, but the OVII constraint remains unchanged, and around the Milky Way a flattened profile may only hold $6-13%$ of the missing baryons from the Galaxy ($2-4 \times 10^{10} M_{\odot}$). We also show that AGN and supernovae at low to moderate redshift - the theoretical sources of winds responsible for driving out the missing baryons - do not produce the expected correlations with the baryonic Tully-Fisher relationship and so are insufficient to explain the missing baryons from galaxies. We conclude that most of missing baryons from galaxies do not lie in hot haloes around the galaxies, and that the missing baryons never fell into the potential wells of protogalaxies in the first place. They may have been expelled from the galaxies as part of the process of galaxy formation.

We simulate the possible emission from a disk perturbed by a recoiling super-massive black hole. To this end, we study radiation transfer from the system incorporating bremsstrahlung emission from a Maxwellian plasma and absorption given by Kramer's opacity law modified to incorporate blackbody effects. We employ this model in the radiation transfer integration to compute the luminosity at several frequencies, and compare with previous bremsstrahlung luminosity estimations from a transparent limit (in which the emissivity is integrated over the computational domain and over all frequencies) and with a simple thermal emission model. We find close agreement between the radiation transfer results and the estimated bremsstrahlung luminosity from previous work for electromagnetic signals above $10^{14}$ Hz. For lower frequencies, we find a self-eclipsing behavior in the disk, resulting in a strong intensity variability connected to the orbital period of the disk.

In addition to producing gravitational waves (GW), the dynamics of a binary black hole system could induce emission of electromagnetic (EM) radiation by affecting the behavior of plasmas and electromagnetic fields in their vicinity. We here study how the electromagnetic fields are affected by a pair of orbiting black holes through the merger. In particular, we show how the binary's dynamics induce a variability in possible electromagnetically induced emissions as well as a possible enhancement of electromagnetic fields during the late-merge and merger epochs. These time dependent features will likely leave their imprint in processes generating detectable emissions and can be exploited in the detection of electromagnetic counterparts of gravitational waves.

We present the radio detection of 18 BL Lac objects from our survey of over 575 square degrees of sky. These 18 objects are located within 20 arcsec of the X-ray position, of which 11 have a measured red-shift. All candidates are radio emitters above ~1 mJy and fall within the range of existing samples on the two colour, alpha_ro vs alpha_ox, diagram with a transitional population of three (3) evident. Two unusual sources have been identified, a candidate radio quiet BL Lac, RX J0140.9-4130, and an extreme HBL, RX J0109.9-4020, with Log(nu_peak)~19.2. The BL Lac Log(N)-Log(S) relation is consistent with other samples and indicates the ROSAT All Sky Survey (RASS) could contain (2000+-400) BL Lac objects.

The merger process of a binary black hole system can have a strong impact on a circumbinary disk. In the present work we study the effect of both central mass reduction (due to the energy loss through gravitational waves) and a possible black hole recoil (due to asymmetric emission of gravitational radiation). For the mass reduction case and recoil directed along the disk's angular momentum, oscillations are induced in the disk which then modulate the internal energy and bremsstrahlung luminosities. On the other hand, when the recoil direction has a component orthogonal to the disk's angular momentum, the disk's dynamics are strongly impacted, giving rise to relativistic shocks. The shock heating leaves its signature in our proxies for radiation, the total internal energy and bremsstrahlung luminosity. Interestingly, for cases where the kick velocity is below the smallest orbital velocity in the disk (a likely scenario in real AGN), we observe a common, characteristic pattern in the internal energy of the disk. Variations in kick velocity simply provide a phase offset in the characteristic pattern implying that observations of such a signature could yield a measure of the kick velocity through electromagnetic signals alone.

We explore the connection between different classes of active galactic nuclei (AGNs) and the evolution of their host galaxies, by deriving host galaxy properties, clustering, and Eddington ratios of AGNs selected in the radio, X-ray, and infrared. We study a sample of 585 AGNs at 0.25 < z < 0.8 using redshifts from the AGN and Galaxy Evolution Survey (AGES) and data in the radio (WSRT 1.4 GHz), X-rays (Chandra XBootes), and mid-IR (IRAC Shallow Survey). The radio, X-ray, and IR AGN samples show modest overlap, indicating that to the flux limits of the survey, they represent largely distinct classes of AGNs. We derive host galaxy colors and luminosities, as well as Eddington ratios (lambda), for obscured or optically faint AGNs. We also measure the two-point cross-correlation between AGNs and galaxies on scales of 0.3-10 h^-1 Mpc, and derive typical dark matter halo masses. We find that: (1) radio AGNs are mainly found in luminous red galaxies, are strongly clustered (with M_halo ~ 3x10^13 h^-1 M_sun), and have very low lambda <~ 10^-3; (2) X-ray-selected AGNs are preferentially found in galaxies in the "green valley" of color-magnitude space and are clustered similarly to typical AGES galaxies (M_halo ~ 10^13 h^-1 M_sun), with 10^-3 <~ lambda <~ 1; (3) IR AGNs reside in slightly bluer, less luminous galaxies than X-ray AGNs, are weakly clustered (M_halo <~ 10^12 h^-1 M_sun), and have lambda > 10^-2. We interpret these results in terms of a simple model of AGN and galaxy evolution, whereby a "quasar" phase and the growth of the stellar bulge occurs when a galaxy's dark matter halo reaches a critical mass between ~10^12 and 10^13 M_sun. Subsequently, star formation ceases and AGN accretion shifts from radiatively efficient (optical- and IR- bright) to radiatively inefficient (optically-faint, radio-bright) modes.

Clusters of galaxies provide a closed box within which one can determine the chemical evolution of the gaseous baryons with cosmic time. We studied this metallicity evolution in the hot X-ray emitting baryons through an analysis of XMM-Newton observations of 29 galaxy clusters in the redshift range 0.3 < z < 1.3. Taken alone, this data set does not show evidence for significant evolution. However, when we also include a comparable sample of 115 clusters observed with Chandra (Maughan et al. 2008) and a lower redshift sample of 70 clusters observed with XMM at z < 0.3 (Snowden et al. 2008), there is definitive evidence for a decrease in the metallicity. This decrease is approximately a factor of two from z = 0 to z \approx 1, over which we find a least-squares best-fit line Z(z) / Z_{\odot} = (0.46 \pm 0.05) - (0.38 \pm 0.03)z. The greatest uncertainty in the evolution comes from poorly constrained metallicities in the highest redshift bin.

We investigate the influence of magnetic fields upon the dynamics of and resulting gravitational waves from a binary neutron star merger in full general relativity coupled to ideal magnetohydrodynamics (MHD). We consider two merger scenarios, one where the stars begin with initially aligned poloidal magnetic fields and one with no magnetic field. Both mergers result in a strongly differentially rotating object. In comparison to the non-magnetized scenario, the aligned magnetic fields delay the final merger of the two stars. During and after merger we observe phenomena driven by the magnetic field, including Kelvin-Helmholtz instabilities in shear layers, winding of the field lines, and transition from poloidal to toroidal fields. These effects not only produce electromagnetic radiation, but also can have a strong influence on the gravitational waves. Thus, there are promising prospects for studying such systems with both types of waves.

We model two mergers of orbiting binary neutron stars, the first forming a black hole and the second a differentially rotating neutron star. We extract gravitational waveforms in the wave zone. Comparisons to a post-Newtonian analysis allow us to compute the orbital kinematics, including trajectories and orbital eccentricities. We verify our code by evolving single stars and extracting radial perturbative modes, which compare very well to results from perturbation theory. The Einstein equations are solved in a first order reduction of the generalized harmonic formulation, and the fluid equations are solved using a modified convex essentially non-oscillatory method. All calculations are done in three spatial dimensions without symmetry assumptions. We use the \had computational infrastructure for distributed adaptive mesh refinement.

Recent developments in ``Einstein Dehn filling'' allow the construction of infinitely many Einstein manifolds that have different topologies but are geometrically close to each other. Using these results, we show that for many spatial topologies, the Hartle-Hawking wave function for a spacetime with a negative cosmological constant develops sharp peaks at certain calculable geometries. The peaks we find are all centered on spatial metrics of constant negative curvature, suggesting a new mechanism for obtaining local homogeneity in quantum cosmology.

We examine the leading order corrections to the Nambu effective action for the motion of a cosmic string, which appear at fourth order in the ratio of the width to radius of curvature of the string. We determine the numerical coefficients of these extrinsic curvature corrections, and derive the equations of motion of the worldsheet. Using these equations, we calculate the corrections to the motion of a collapsing loop, a travelling wave, and a helical breather. From the numerical coefficients we have calculated, we discuss whether the string motion can be labelled as `rigid' or `antirigid,' and hence whether cusp or kink formation might be suppressed or enhanced.

We compare the radio and soft X-ray brightness as a function of position within the young supernova remnant Cassiopeia A. A moderately strong correlation (r = 0.7) was found between the X-ray emission (corrected for interstellar absorption) and radio emission, showing that the thermal and relativistic plasmas occupy the same volumes and are regulated by common underlying parameters. The logarithmic slope of the relationship, ln(Sx-ray) = 1.2 x Sradio + ln(k) implies that the variations in brightness are primarily due to path length differences. The X-ray and radio emissivities are both high in the same general locations, but their more detailed relationship is poorly constrained and probably shows significant scatter. The strongest radio and X-ray absorption is found at the western boundary of Cas A. Based on the properties of Cas A and the absorbing molecular cloud, we argue that they are physically interacting. We also compare ASCA derived column densities with 21 cm H I and 18 cm OH optical depths in the direction of Cas A, in order to provide an independent estimate of ISM properties. We derive an average value for the H I spin temperature of about 40 K and measure the ratio of OH to molecular hydrogen to be nominally larger than previous estimates. Keywords: Cas A, Cassiopeia A, interstellar medium, molecular clouds, radio astronomy, supernova remnants, X-ray astronomy