Closure statistics in interferometric data

Interferometric visibilities, reflecting the complex correlations between signals recorded at antennas in an interferometric array, carry information about the angular structure of a distant source. While unknown antenna gains in both amplitude and phase can prevent direct interpretation of these measurements, certain combinations of visibilities called closure phases and closure amplitudes are independent of antenna gains and provide a convenient set of robust observables. However, these closure quantities have subtle noise properties and are generally both linearly and statistically dependent. These complications have obstructed the proper use of closure quantities in interferometric analysis, and they have obscured the relationship between analysis with closure quantities and other analysis techniques such as self calibration. We review the statistics of closure quantities, noting common pitfalls that arise when approaching low signal-to-noise due to the nonlinear propagation of statistical errors. We then develop a strategy for isolating and fitting to the independent degrees of freedom captured by the closure quantities through explicit construction of linearly independent sets of quantities along with their noise covariance in the Gaussian limit, valid for moderate signal-to-noise, and we demonstrate that model fits have biased posteriors when this covariance is ignored. Finally, we introduce a unified procedure for fitting to both closure information and partially calibrated visibilities, and we demonstrate both analytically and numerically the direct equivalence of inference based on closure quantities to that based on self calibration of complex visibilities with unconstrained antenna gains.Comment: 31 pages, 17 figure

The Megamaser Cosmology Project. X. High Resolution Maps and Mass Constraint for SMBHs

We present high resolution (sub-mas) VLBI maps of nuclear H2O megamasers for seven galaxies. In UGC6093, the well-aligned systemic masers and high-velocity masers originate in an edge-on, flat disk and we determine the mass of the central SMBH to be M_SMBH = 2.58*10^7Msun(+-7%). For J1346+5228, the distribution of masers is consistent with a disk, but the faint high-velocity masers are only marginally detected, and we constrain the mass of the SMBH to be in the range 1.5-2.0*10^7Msun. The origin of the masers in Mrk1210 is less clear, as the systemic and high-velocity masers are misaligned and show a disorganized velocity structure. We present one possible model in which the masers originate in a tilted, warped disk, but we do not rule out the possibility of other explanations including outflow masers. In NGC6926, we detect a set of redshifted masers, clustered within a pc of each other, and a single blueshifted maser about 4.4pc away, an offset that would be unusually large for a maser disk system. Nevertheless, if it is a disk system, we estimate the enclosed mass to be M_SMBH<4.8*10^7 Msun . For NGC5793, we detect redshifted masers spaced about 1.4pc from a clustered set of blueshifted features. The orientation of the structure supports a disk scenario as suggested by Hagiwara et al.(2001). We estimate the enclosed mass to be M SMBH<1.3*10^7 Msun. For NGC2824 and J0350-0127, the masers may be associated with pc or sub-pc scale jets or outflows.Comment: Accepted by Ap

Measuring the Polarization of a Rapidly Precessing Deuteron Beam

This paper describes a time-marking system that enables a measurement of the in-plane (horizontal) polarization of a 0.97-GeV/c deuteron beam circulating in the Cooler Synchrotron (COSY) at the Forschungszentrum J\"ulich. The clock time of each polarimeter event is used to unfold the 120-kHz spin precession and assign events to bins according to the direction of the horizontal polarization. After accumulation for one or more seconds, the down-up scattering asymmetry can be calculated for each direction and matched to a sinusoidal function whose magnitude is proportional to the horizontal polarization. This requires prior knowledge of the spin tune or polarization precession rate. An initial estimate is refined by re-sorting the events as the spin tune is adjusted across a narrow range and searching for the maximum polarization magnitude. The result is biased toward polarization values that are too large, in part because of statistical fluctuations but also because sinusoidal fits to even random data will produce sizeable magnitudes when the phase is left free to vary. An analysis procedure is described that matches the time dependence of the horizontal polarization to templates based on emittance-driven polarization loss while correcting for the positive bias. This information will be used to study ways to extend the horizontal polarization lifetime by correcting spin tune spread using ring sextupole fields and thereby to support the feasibility of searching for an intrinsic electric dipole moment using polarized beams in a storage ring. This paper is a combined effort of the Storage Ring EDM Collaboration and the JEDI Collaboration.Comment: 28 pages, 15 figures, prepared for Physical Review ST - Accelerators and Beam

The HST Survey of BL Lacertae Objects. II. Host Galaxies

We have used the HST WFPC2 camera to survey 132 BL Lac objects comprising seven complete radio-, X-ray-, and optically-selected samples. We obtained useful images for 110 targets spanning the redshift range 0 < z < 1.3. In two thirds of the BL Lac images, host galaxies are detected, including nearly all for z < 0.5 (58 of 63). The highest redshift host galaxy detected is in a BL Lac object at z=0.664. In 58 of the 72 resolved host galaxies, a de Vaucouleurs profile is significantly preferred, at >99% confidence, over a pure exponential disk; the two fits are comparable in the remaining 14 cases. These results limit the number of disk systems to at most 8% of BL Lacs (at 99% confidence), and are consistent with all BL~Lac host galaxies being ellipticals. The detected host galaxies are luminous ellipticals with a median absolute K-corrected magnitude of M_R= -23.7 +- 0.6 mag, at least one magnitude brighter than M* and comparable to brightest cluster galaxies. The galaxy morphologies are generally smooth and undisturbed, with small or negligible ellipticities (<0.2). There is no correlation between host galaxy and observed nuclear magnitude or estimated jet power corrected for beaming. If black hole mass is correlated linearly with bulge mass in general, this implies a large range in Eddington ratio. Present data strongly support the unification picture with FR I galaxies constituting the bulk of the parent population of BL Lac objects.Comment: Accepted for publication on ApJ. 38 pages, 8 figure

Observing Ultra High Energy Cosmic Particles from Space: SEUSO, the Super Extreme Universe Space Observatory Mission

The experimental search for ultra high energy cosmic messengers, from $E\sim 10^{19}$ eV to beyond $E\sim 10^{20}$ eV, at the very end of the known energy spectrum, constitutes an extraordinary opportunity to explore a largely unknown aspect of our universe. Key scientific goals are the identification of the sources of ultra high energy particles, the measurement of their spectra and the study of galactic and local intergalactic magnetic fields. Ultra high energy particles might, also, carry evidence of unknown physics or of exotic particles relics of the early universe. To meet this challenge a significant increase in the integrated exposure is required. This implies a new class of experiments with larger acceptances and good understanding of the systematic uncertainties. Space based observatories can reach the instantaneous aperture and the integrated exposure necessary to systematically explore the ultra high energy universe. In this paper, after briefly summarising the science case of the mission, we describe the scientific goals and requirements of the SEUSO concept. We then introduce the SEUSO observational approach and describe the main instrument and mission features. We conclude discussing the expected performance of the mission

Experimental Implementation of the Deutsch-Jozsa Algorithm for Three-Qubit Functions using Pure Coherent Molecular Superpositions

The Deutsch-Jozsa algorithm is experimentally demonstrated for three-qubit functions using pure coherent superpositions of Li$_{2}$ rovibrational eigenstates. The function's character, either constant or balanced, is evaluated by first imprinting the function, using a phase-shaped femtosecond pulse, on a coherent superposition of the molecular states, and then projecting the superposition onto an ionic final state, using a second femtosecond pulse at a specific time delay