Search

We describe for the first time in the scientific literature the Planetary Ephemeris Program (PEP), an open-source general-purpose astrometric data analysis program. We discuss, in particular, the implementation of pulsar timing analysis, which was recently upgraded in PEP to handle more options. This implementation was done independently of other pulsar programs, with minor exceptions that we discuss. We illustrate the implementation of this capability by comparing the post-fit residuals from the analyses of time-of-arrival observations by both PEP and Tempo2. The comparison shows substantial agreement: 22 ns rms differences for 1,065 pulse time-of-arrival measurements for the millisecond pulsar in a binary system, PSR J1909-3744 (pulse period 2.947108 ms; full-width half-maximum of pulse 43 $\mu$s) for epochs in the interval from December 2002 to February 2011., Comment: 9 pages, 1 figure, submitted to Astronomical Journal (AJ)

A long standing goal in astrophysics is to directly observe the immediate environment of a black hole with angular resolution comparable to the event horizon. Realizing this goal would open a new window on the study of General Relativity in the strong field regime, accretion and outflow processes at the edge of a black hole, the existence of an event horizon, and fundamental black hole physics (e.g., spin). Steady long-term progress on improving the capability of Very Long Baseline Interferometry (VLBI) at short wavelengths has now made it extremely likely that this goal will be achieved within the next decade. The most compelling evidence for this is the recent observation by 1.3mm VLBI of Schwarzschild radius scale structure in SgrA*, the compact source of radio, submm, NIR and xrays at the center of the Milky Way. SgrA* is thought to mark the position of a ~4 million solar mass black hole, and because of its proximity and estimated mass presents the largest apparent event horizon size of any black hole candidate in the Universe. Over the next decade, existing and planned mm/submm facilities will be combined into a high sensitivity, high angular resolution "Event Horizon Telescope" that will bring us as close to the edge of black hole as we will come for decades. This white paper describes the science case for mm/submm VLBI observations of both SgrA* and M87 (a radio loud AGN of a much more luminous class that SgrA*). We emphasize that while there is development and procurement involved, the technical path forward is clear, and the recent successful observations have removed much of the risk that would normally be associated with such an ambitious project., Comment: 8 pages, Science White Paper submitted to the ASTRO2010 Decadal Review Panels

The recent discovery of "dark energy" has challenged Einstein's general theory of relativity as a complete model for our macroscopic universe. From a theoretical view, the challenge is even stronger: general relativity clearly does not extend to the very small, where quantum mechanics holds sway. Fundamental physics models thus require some major revisions. We must explore deeper to both constrain and inspire this needed new physics. In the realm of the solar-system, we can effectively probe for small deviations from the predictions of general relativity: Technology now offers a wide range of opportunities to pursue experiments with accuracies orders of magnitude better than yet achieved. We describe both the relevant theoretical backgrounds and the opportunities for far more accurate solar system experiments., Comment: 8 pages, LaTeX, Astro2010 Science White Paper

In 1981, a faint radio source (G') was detected near the center of the lensing galaxy of the famous "twin quasar" Q0957+561. It is still unknown whether this central radio source is a third quasar image or an active nucleus of the lensing galaxy, or a combination of both. In an attempt to resolve this ambiguity, we observed Q0957+561 at radio wavelengths of 13cm, 18cm, and 21cm, using the Very Long Baseline Array in combination with the phased Very Large Array and the Green Bank Telescope. We measured the spectrum of G' for the first time and found it to be significantly different from the spectra of the two bright quasar images. This finding suggests that the central component is primarily or entirely emission from the foreground lens galaxy, but the spectrum is also consistent with the hypothesis of a central quasar image suffering free-free absorption. In addition, we confirm the previously-reported VLBI position of G' just north of the optical center of the lens galaxy. The position slightly favors the hypothesis that G' originates in the lens, but is not conclusive. We discuss the prospects for further clarification of this issue., Comment: 18 pages, accepted for publication in AJ

We examine models of the mass distribution for the first known case of gravitational lensing. Several new sets of constraints are used, based on recent observations. We remodel the VLBI observations of the radio jets in the two images of Q0957+561, showing that the previously derived positions and uncertainties were incorrect. We use as additional constraints the candidate lensed pairs of galaxies identified recently with the Hubble Space Telescope. We explore a wider range of lens models than before, and find that the Hubble constant is not tightly constrained once elliptical lens models are considered. We also discuss the systematic uncertainties caused by the cluster containing the lens galaxy. We conclude that additional observations of the candidate lensed galaxies as well as direct measurements of the cluster mass profile are needed before an accurate value of the Hubble constant can be derived from this lens., Comment: 28 pages, Latex + 4 figures, accepted for ApJ. Final version, with minor changes. Model parameters in Table 4 corrected

We present long-slit LRIS/Keck spectroscopic observations of the gravitational lens system 0957+561. Averaged over all of our data, the rest-frame velocity dispersion sigma_v of the central lens galaxy G1 is sigma_v = 279 +- 12 km/s. However, there appears to be a significant decrease in sigma_v as a function of distance from the center of G1 that is not typical of brightest cluster galaxies. Within 0.2" of the center of G1, we find the average sigma_v = 316+- 14 km/s, whereas for positions > 0.2" from the center of G1, we find the average sigma_v= 266+-12 km/s. A plausible explanation is that G1 contains a central massive dark object of mass M_MDO ~ 4x10^9 h_100^{-1} M_sun (h_100 = H_0/100km/s/Mpc), which contributes to the central velocity dispersion, and that the outer value of sigma_v is the appropriate measure of the depth of the potential well of G1. The determination of a luminosity-weighted estimate of sigma_v is essential for a determination of H_0 from 0957+561; our accurate measurements remove one of the chief uncertainties in the assumed form of the mass distribution of the lens. Thus, with the recent apparent reduction in the uncertainty in the measurement of the time delay for the images A and B of 0957+561, Delta tau_BA = 417+-3 d (Kundic et al. 1996), we obtain an estimate for the Hubble constant: H_0 = 62+-7 km/s/Mpc. If for some reason the trend of sigma_v with slit position is spurious and we should use the dispersion averaged along the slit, then the estimate of H_0 increases to 67+-8 km/s/Mpc. These standard errors do not, however, include any contribution from any errors in the assumed form of the mass distribution of the lens. In particular, we used the mass model described by Falco, Gorenstein & Shapiro (1991), as updated by Grogin & Narayan (1996a, b). The reduced chi^2, Comment: 29 pages, 12 .eps figs, AAS Latex, ApJ, in press

Quadruple gravitational lens systems offer the possibility of measuring time delays for image pairs, microlensing effects, and extinction in distant galaxies. Observations of these systems may be used to obtain estimates of H_o and to study the various mass components of lens galaxies at high redshifts. With HST, we have observed the reddest known gravitational lens system, MG 0414+0534. We used WFPC2/PC1 to obtain deep, high-resolution images with two filters, F675W and F814W. We present a detailed analysis of all of the components, as well as macrolens models. Our main results are: (1) confirmation that MG 0414+0534 is inescapably a gravitational lens system; (2) discovery of a blue arc connecting the 3 brightest images of the QSO central core; (3) accurate positions and apparent brightnesses for all 4 known images of the QSO central core and for the lens galaxy G; (4) a good representation of the brightness distribution of G by elliptical isophotes with a De Vaucouleurs profile, characteristic of an elliptical galaxy; (5) models that consist of simple elliptical potentials and account qualitatively, not quantitatively, for the HST image positions, arc morphology and radio flux ratios for the images of the QSO central core; (6) a possible new test to distinguish between reddening in the host galaxy of the QSO and in the lens galaxy, based on future accurate measurements of spatial variations in the color of the arc; and (7) the suggestion that microlensing is a plausible cause for the differences between the radio and optical flux ratios for the brightest images, A1 and A2. Further observations and measurements such as of the redshift of the lens galaxy, can be used fruitfully to study microlensing for this system., Comment: 27 pages, 8 .ps figs, AAS Latex, AJ, in press, Feb 1997

This research is aimed at testing gravitational theory, primarily on an interplanetary scale and using mainly observations of objects in the solar system. Our goal is either to detect departures from the standard model (general relativity) - if any exist within the level of sensitivity of our data - or to support this model by placing tighter bounds on any departure from it. For this project, we have analyzed a combination of observational data with our model of the solar system, including planetary radar ranging, lunar laser ranging, and spacecraft tracking, as well as pulsar timing and pulsar VLBI measurements.

This Annual Report illustrates the work carried out during the last grant-year activity on the Test of the Equivalence Principle in an Einstein Elevator. The activity focused on the following main topics: (1) analysis and conceptual design of a detector configuration suitable for the flight tests; (2) development of techniques for extracting a small signal from data strings with colored and white noise; (3) design of the mechanism that spins and releases the instrument package inside the cryostat; and (4) experimental activity carried out by our non-US partners (a summary is shown in this report). The analysis and conceptual design of the flight-detector (point 1) was focused on studying the response of the differential accelerometer during free fall, in the presence of errors and precession dynamics, for various detector's configurations. The goal was to devise a detector configuration in which an Equivalence Principle violation (EPV) signal at the sensitivity threshold level can be successfully measured and resolved out of a much stronger dynamics-related noise and gravity gradient. A detailed analysis and comprehensive simulation effort led us to a detector's design that can accomplish that goal successfully.

The scientific goal of the experiment is to test the equality of gravitational and inertial mass (i.e., to test the Principle of Equivalence) by measuring the independence of the rate of fall of bodies from the composition of the falling body. The measurement is accomplished by measuring the relative displacement (or equivalently acceleration) of two falling bodies of different materials which are the proof masses of a differential accelerometer. The goal of the experiment is to measure the Eoetvoes ratio sigma g/g (differential acceleration/common acceleration) with an accuracy goal of several parts in 10(exp 15). The estimated accuracy is about two orders of magnitude better than the present state of the art. The main goal of the study to be carried out under this grant is part of the flight definition of the experiment and laboratory testing of key components. The project involves an international cooperation in which the responsibility of the US side is the flight definition of the experimental facility while the responsibility of the non-US partners is the flight definition and laboratory prototyping of the differential acceleration detector.In summary, the experiment to be designed is for taking differential acceleration measurements with a high-sensitivity detector (the sensor) during free fall conditions lasting up to 30 s in a disturbance-free acceleration environment. The experiment strategy consists in letting the sensor free fall inside a few meters long (in the vertical direction) evacuated capsule that is falling simultaneously in the rarefied atmosphere after release from a helium balloon flying at a stratospheric altitude.

The present report covers the first year of a grant which represents a direct continuation of NASA NAG5-4050, with the same title as before. It is dedicated as before to the discovery and characterization of new astrophysical molecules. This year, has been extremely productive, yielding many new discoveries of astronomical interest at both radio and optical wavelengths, and the publication or submission of the 15 papers listed below. Nearly all of these articles have or will soon appear in the leading refereed journals of astrophysics, chemical physics, physics, or molecular spectroscopy. One is a major invited review for Molecular Physics. One of our other invited reviews published in Spectrochimica Acta in 2001 was recently awarded the Sir Harold Thompson Memorial Award, annually given to the best paper in that journal. During the past year significant advances have been made by our group in the laboratory study of exotic silicon and carbon molecules of astronomical interest. The most exciting discoveries include the pure silicon cluster Si3, several novel silicon hydrides, and the detection of phenyl radical, C6H5, a fundamental reactive organic ring. In addition, the rotational spectra of many carbon chains terminated with Si, N, O, and other heteroatoms have also been detected for the first time. The laboratory astrophysics of the whole set is complete in the sense that the entire radio spectrum of each species has now been measured or can be calculated to very high accuracy. Nearly all of these newly found molecules are plausible candidates for the detection by radio astronomers in the interstellar gas or in circumstellar sources because they are similar in structure and composition to known astronomical species, and because most are calculated to possess large permanent dipole moments.

We are engaged in testing gravitational theory, mainly using observations of objects in the solar system and mainly on the interplanetary scale. Our goal is either to detect departures from the standard model (general relativity) - if any exist within the level of sensitivity of our data - or to place tighter bounds on such departures. For this project, we have analyzed a combination of observational data with our model of the solar system, including primarily planetary radar ranging, lunar laser ranging, and spacecraft tracking, but also including both pulsar timing and pulsar VLBI measurements. In the past year, we have included new data in the analysis, primarily tracking data from the Mars Pathfinder mission. Although these data are relatively few in number, they extend the time span of high-precision tracking on the surface of Mars from six years to over 20. As a result, the statistical standard deviation of our estimate of Mars precession rate has nearly halved, and the rest of the parameters in our solar-system model have experienced a corresponding, albeit smaller, improvement (about 20% for t,he relevant asteroid masses, 10% for the semimajor axis of Mars orbit, and smaller amounts for most other parameters). In the coming year, we plan to continue adding data to our set, as available. Ne 2 expect to use these data and improved models to obtain estimates of the gravitational- theory parameters and to publish these results.