Parallax Distance Research Articles

A Bayesian approach to calibrating period-luminosity relations of RR Lyrae stars in the mid-infrared

A Bayesian approach to calibrating period-luminosity (PL) relations has substantial benefits over generic least-squares fits. In particular, the Bayesian approach takes into account the full prior distribution of the model parameters, such as the a priori distances, and refits these parameters as part of the process of settling on the most highly-constrained final fit. Additionally, the Bayesian approach can naturally ingest data from multiple wavebands and simultaneously fit the parameters of PL relations for each waveband in a procedure that constrains the parameter posterior distributions so as to minimize the scatter of the final fits appropriately in all wavebands. Here we describe the generalized approach to Bayesian model fitting and then specialize to a detailed description of applying Bayesian linear model fitting to the mid-infrared PL relations of RR Lyrae variable stars. For this example application we quantify the improvement afforded by using a Bayesian model fit. We also compare distances previously predicted in our example application to recently published parallax distances measured with the Hubble Space Telescope and find their agreement to be a vindication of our methodology. Our intent with this article is to spread awareness of the benefits and applicability of this Bayesian approach and encourage future PL relation investigations to consider employing this powerful analysis method.

VLBA DETERMINATION OF THE DISTANCE TO NEARBY STAR-FORMING REGIONS. V. DYNAMICAL MASS, DISTANCE, AND RADIO STRUCTURE OF V773 Tau A

(ABRIDGED) We present multi-epoch Very Long Baseline Array (VLBA) observations of V773 Tau A, the 51-day binary subsystem in the multiple young stellar system V773 Tau. Combined with previous interferometric and radial velocity measurements, these new data enable us to improve the characterization of the physical orbit of the A subsystem. In particular, we infer updated dynamical masses for the primary and the secondary components of 1.55 pm 0.11 Msun, and 1.293 pm 0.068 Msun, respectively, and an updated orbital parallax distance to the system of 135.7 pm 3.2 pc, all consistent with previous estimates. Using the improved orbit, we can calculate the absolute coordinates of the barycenter of V773 Tau A at each epoch of our VLBA observations, and fit for its trigonometric parallax and proper motion. This provides a direct measurement of the distance to the system almost entirely independent of the orbit modeling. The best fit yields a distance of 129.9 pm 3.2 pc, in good agreement (i.e. within 1 sigma) with the distance estimate based on the orbital fit. Taking the mean value of the orbital and trigonometric parallaxes, we conclude that V773 Tau is located at d =132.8 pm 2.3 pc. The accuracy of this determination is nearly one order of magnitude better than that of previous estimates. In projection, V773 Tau and two other young stars (Hubble 4 and HDE 283572) recently observed with the VLBA are located toward the dark cloud Lynds 1495, in the central region of Taurus. These three stars appear to have similar trigonometric parallaxes, radial velocities, and proper motions, and we argue that the weighted mean and dispersion of their distances (d = 131.4 pc and sigma_d = 2.4 pc) provide a good estimate of the distance to and depth of Lynds 1495 and its associated stellar population.

XMM-NEWTONOBSERVATIONS OF THE DWARF NOVA RU Peg IN QUIESCENCE: PROBE OF THE BOUNDARY LAYER

We present an analysis of X-ray and UV data obtained with the XMM-Newton Observatory of the long period dwarf nova RU Peg. RU Peg contains a massive white dwarf, possibly the hottest white dwarf in a dwarf nova, it has a low inclination, thus optimally exposing its X-ray emitting boundary layer, and has an excellent trigonometric parallax distance. We modeled the X-ray data using XSPEC assuming a multi-temperature plasma emission model built from the MEKAL code. We obtained a maximum temperature of 31.7 keV, based on the EPIC MOS1, 2 and pn data, indicating that RU Peg has an X-ray spectrum harder than most dwarf novae, except U Gem. This result is consistent with and indirectly confirms the large mass of the white dwarf in RU Peg. The X-ray luminosity we computed corresponds to a boundary layer luminosity for a mass accretion rate of 2.E-11 Msun/yr (assuming Mwd=1.3Msun), in agreement with an expected quiescent accretion rate. The modeling of the O VIII emission line at 19A as observed by the RGS implies a projected stellar rotational velocity of 695 km/s, i.e. the line is emitted from material rotating at about 936-1245 km/s (for i about 34-48deg) or about 1/6 of the Keplerian speed; this velocity is much larger than the rotation speed of the white dwarf inferred from the FUSE spectrum. Cross-correlation analysis yielded an undelayed component and a delayed component of 116 +/- 17 sec where the X-ray variations/fluctuations lagged the UV variations. This indicates that the UV fluctuations in the inner disk are propagated into the X-ray emitting region in about 116 sec. The undelayed component may be related to irradiation effects.

Calibrating the Cepheid period-luminosity relation from the infrared surface brightness technique

We determine Period-Luminosity relations for Milky Way Cepheids in the optical and near-IR bands. These relations can be used directly as reference for extra-galactic distance determination to Cepheid populations with solar metallicity, and they form the basis for a direct comparison with relations obtained in exactly the same manner for stars in the Magellanic Clouds, presented in an accompanying paper. In that paper we show that the metallicity effect is very small and consistent with a null effect, particularly in the near-IR bands, and we combine here all 111 Cepheids from the Milky Way, the LMC and SMC to form a best relation. We employ the near-IR surface brightness (IRSB) method to determine direct distances to the individual Cepheids after we have recalibrated the projection factor using the recent parallax measurements to ten Galactic Cepheids and the constraint that Cepheid distances to the LMC should be independent of pulsation period. We confirm our earlier finding that the projection factor for converting radial velocity to pulsational velocity depends quite steeply on pulsation period, p=1.550-0.186*log(P) in disagrement with recent theoretical predictions. We delineate the Cepheid PL relation using 111 Cepheids with direct distances from the IRSB analysis. The relations are by construction in agreement with the recent HST parallax distances to Cepheids and slopes are in excellent agreement with the slopes of apparent magnitudes versus period observed in the LMC.

Trigonometric parallax measurements of post-AGB stars with VERA

AbstractA trigonometric parallax distance is measurable for a planetary nebula if it harbors H2O (or OH) maser spots. We have demonstrated it for H2O maser sources, K3-35, IRAS 19312+1950 and IRAS 18286-0959 using the VLBI Exploration of Radio Astrometry (VERA). They are post-AGB stars and exhibit peculiar morpho-kinematical structures like a bipolar planetary nebula or a bipolar molecular jet (water fountain). The luminosities and secular motions of these sources derived from the trigonometric distances and proper motions of the H2O masers suggest that their parental stars may be intermediate-mass (2M⊙ < M* < 8 M⊙) evolved stars.

THE EVOLUTION OF CATACLYSMIC VARIABLES AS REVEALED BY THEIR DONOR STARS

ABSTRACT We present an attempt to reconstruct the complete evolutionary path followed by cataclysmic variables (CVs), based on the observed mass–radius relationship of their donor stars. Along the way, we update the semi-empirical CV donor sequence presented previously by one of us, present a comprehensive review of the connection between CV evolution and the secondary stars in these systems, and reexamine most of the commonly used magnetic braking (MB) recipes, finding that even conceptually similar ones can differ greatly in both magnitude and functional form. The great advantage of using donor radii to infer mass-transfer and angular-momentum-loss (AML) rates is that they sample the longest accessible timescales and are most likely to represent the true secular (evolutionary average) rates. We show explicitly that if CVs exhibit long-term mass-transfer-rate fluctuations, as is often assumed, the expected variability timescales are so long that other tracers of the mass-transfer rate—including white dwarf (WD) temperatures—become unreliable. We carefully explore how much of the radius difference between CV donors and models of isolated main-sequence stars may be due to mechanisms other than mass loss. The tidal and rotational deformation of Roche-lobe-filling stars produces ≃ 4.5% radius inflation below the period gap and ≃ 7.9% above. A comparison of stellar models to mass–radius data for non-interacting stars suggests a real offset of ≃ 1.5% for fully convective stars (i.e., donors below the gap) and ≃ 4.9% for partially radiative ones (donors above the gap). We also show that donor bloating due to irradiation is probably smaller than, and at most comparable to, these effects. After calibrating our models to account for these issues, we fit self-consistent evolution sequences to our compilation of donor masses and radii. In the standard model of CV evolution, AMLs below the period gap are assumed to be driven solely by gravitational radiation (GR), while AMLs above the gap are usually described by an MB law first suggested by Rappaport et al. We adopt simple scaled versions of these AML recipes and find that these are able to match the data quite well. The optimal scaling factors turn out to be f GR = 2.47 ± 0.22 below the gap and f MB = 0.66 ± 0.05 above (the errors here are purely statistical, and the standard model corresponds to f GR = f MB = 1). This revised model describes the mass–radius data significantly better than the standard model. Some of the most important implications and applications of our results are as follows. (1) The revised evolution sequence yields correct locations for the minimum period and the upper edge of the period gap; the standard sequence does not. (2) The observed spectral types of CV donors are compatible with both standard and revised models. (3) A direct comparison of predicted and observed WD temperatures suggests an even higher value for f GR, but this comparison is sensitive to the assumed mean WD mass and the possible existence of mass-transfer-rate fluctuations. (4) The predicted absolute magnitudes of donor stars in the near-infrared form a lower envelope around the observed absolute magnitudes for systems with parallax distances. This is true for all of our sequences, so any of them can be used to set firm lower limits on (or obtain rough estimates of) the distances toward CVs based only on P orb and single epoch near-IR measurements. (5) Both standard and revised sequences predict that short-period CVs should be susceptible to dwarf nova (DN) eruptions, consistent with observations. However, both sequences also predict that the fraction of DNe among long-period CVs should decline with P orb above the period gap. Observations suggest the opposite behavior, and we discuss the possible explanations for this discrepancy. (6) Approximate orbital period distributions constructed from our evolution sequences suggest that the ratio of long-period CVs to short-period, pre-bounce CVs is about 3 × higher for the revised sequence than the standard one. This may resolve a long-standing problem in CV evolution. Tables describing our donor and evolution sequences are provided in electronically readable form.

Annual Parallax Distance and Kinematical Property of H2O Masers in IRAS 19312+1950

Abstract We report on results of astrometric observations of H$_{2}$O masers in the pre-planetary nebula candidate source IRAS 19312$+$1950 with the VLBI Exploration of Radio Astrometry (VERA), yielding an annual parallax, $\pi$$=$ 263$\ \pm\ $47$\ \mu$as, and a secular motion, ($\mu_{\alpha}$, $\mu_{\delta}$) $=$ ($-$2.61$\ \pm\ $0.47, $-$6.73$\ \pm\ $0.14) [mas yr$^{-1}$], for IRAS 19312$+$1950. We then derived a heliocentric distance of $D$$=$ 3.80$^{+0.83}_{\-0.58}\ $kpc, a location in the Galaxy ($R$, $z$) $=$ (7.07$\ \pm\ $0.12 kpc, 28$\ \pm\ $3 pc), and a three-dimensional velocity vector ($V_{R}$, $V_{\theta}$, $V_{z}$) $=$ (33$\ \pm\ $28, 214$\ \pm\ $4, $-$14$\ \pm\ $8) [km s$^{-1}$] in galactic cylindrical coordinates. Based on these parameters, we discuss the intrinsic properties of IRAS 19312$+$1950. We estimated the mass of the central star to be $\sim\ $8$\ M_{\odot}$ and the deviation of the source motion from the circular galactic rotation to be $\sim\ $40 km s$^{-1}$. These are consistent with this source being an intermediate-mass late AGB or early post-AGB star.

ROTATIONALLY MODULATEDg-MODES IN THE RAPIDLY ROTATING δ SCUTI STAR RASALHAGUE (α OPHIUCHI)

Despite a century of remarkable progress in understanding stellar interiors, we know surprisingly little about the inner workings of stars spinning near their critical limit. New interferometric imaging of these so-called ``rapid rotators'' combined with breakthroughs in asteroseismology promise to lift this veil and probe the strongly latitude-dependent photospheric characteristics and even reveal the internal angular momentum distribution of these luminous objects. Here, we report the first high precision photometry on the low-amplitude delta cuti variable star Rasalhague (alpha Oph, A5IV, 2.18 Msun, omega/omega_c~0.88) based on 30 continuous days of monitoring using the MOST satellite. We have identified 57+/-1 distinct pulsation modes above a stochastic granulation spectrum with a cutoff of ~26 cycles per day. Remarkably, we have also discovered that the fast rotation period of 14.5~hours modulates low-frequency modes (1-10 day periods) that we identify as a rich family of g-modes (|m| up to 7). The spacing of the g-modes is surprisingly linear considering Coriolis forces are expected to strongly distort the mode spectrum, suggesting we are seeing prograde ``equatorial Kelvin'' waves (modes l=m). We emphasize the unique aspects of Rasalhague motivating future detailed asteroseismic modeling -- a source with a precisely measured parallax distance, photospheric oblateness, latitude temperature structure, and whose low-mass companion provides an astrometric orbit for precise mass determinations.

Spatial distribution of interstellar dust in the Sun's vicinity

3D tomography of the interstellar dust and gas may be useful in many respects, from the physical and chemical evolution of the ISM itself to foreground decontamination of the CMB, or various studies of the environments of specific objects. Our goal here is to bring more precise information on the distance to nearby interstellar dust and gas clouds within 250 pc. We apply the best available calibration methods to a carefully screened set of stellar Stromgren photometry data for targets possessing a Hipparcos parallax and spectral type classification. We combine the derived interstellar extinctions and the parallax distances for about 6,000 stars to build a 3D tomography of the local dust. We use an inversion method based on a regularized Bayesian approach and a least squares criterion. We obtain 3D maps of the opacity and the distance to the main dust-bearing clouds with 250 pc. We calculate the integrated extinction between the Sun and the cube boundary and compare with the total galactic extinction derived from infrared 2D maps. We also compare the opacity distribution with the 3D distribution of interstellar neutral sodium resulting from the inversion of sodium columns towards about 1,700 target stars. Our maps show a larger high latitude dust opacity in the North compared to the South, while gas maps do not show the same asymmetry, suggesting a polar asymmetry of the dust to gas ratio at small distances.

TAMING THE INVISIBLE MONSTER: SYSTEM PARAMETER CONSTRAINTS FOR ϵ AURIGAE FROM THE FAR-ULTRAVIOLET TO THE MID-INFRARED

We have assembled new Spitzer Space Telescope Infrared Array Camera observations of the mysterious binary star Epsilon Aurigae, along with archival far-ultraviolet to mid-infrared data, to form an unprecedented spectral energy distribution spanning three orders of magnitude in wavelength from 0.1 microns to 100 microns. The observed spectral energy distribution can be reproduced using a three component model consisting of a 2.2+0.9/-0.8 Msun F type post-asymptotic giant branch star, and a 5.9+/-0.8 Msun B5+/-1 type main sequence star that is surrounded by a geometrically thick, but partially transparent, disk of gas and dust. At the nominal HIPPARCOS parallax distance of 625 pc, the model normalization yields a radius of 135+/-5 Rsun for the F star, consistent with published interferometric observations. The dusty disk is constrained to be viewed at an inclination of i > 87 deg, and has effective temperature of 550+/-50 K with an outer radius of 3.8 AU and a thickness of 0.95 AU. The dust content of the disk must be largely confined to grains larger than ~10 microns in order to produce the observed gray optical-infrared eclipses and the lack of broad dust emission features in the archival Spitzer mid-infrared spectra. The total mass of the disk, even considering a potential gaseous contribution in addition to the dust that produces the observed infrared excess, is << 1 Msun. We discuss evolutionary scenarios for this system that could lead to the current status of the stellar components and suggests possibilities for its future evolution, as well as potential observational tests of our model.

Far Ultraviolet Spectroscopy of Seven Nova-Like Variables

ABSTRACTWe present the results of a multicomponent synthetic spectral analysis of the archival far ultraviolet spectra of several key nova-like variables including members of the SW Sex, RW Tri, UX UMa, and VY Scl subclasses: KR Aur, RW Tri, V825 Her, V795 Her, BP Lyn, V425 Cas, and HL Aqr. Accretion rates as well as the possible flux contribution of the accreting white dwarf are included in our analysis. Except for RW Tri, which has a reliable trigonometric parallax, we computed the distances to the nova-like systems using the method of Knigge. Our analysis of seven archival IUE spectra of RW Tri at its parallax distance of 341 pc consistently indicates a low mass (∼0.4 M⊙) white dwarf and an average accretion rate, . For KR Aur, we estimate that the white dwarf has Teff = 29,000 ± 2000 K, log g = 8.4, and contributes 18% of the far-UV flux, while an accretion disk with accretion rate at an inclination of 41° contributes the remainder. We find that an accretion disk dominates the far-UV spectrum of V425 Cas but a white dwarf contributes nonnegligibly with approximately 18% of the far-UV flux. For the two high state nova-likes, HL Aqr and V825 Her, their accretion disks totally dominate with and 3 × 10-9 M⊙ yr-1, respectively. For BP Lyn we find while for V795 Her, we find an accretion rate of . We discuss the implications of our results for the evolutionary status of nova-like variables.

THE EXPANDING BIPOLAR SHELL OF THE HELIUM NOVA V445 PUPPIS

From multi-epoch adaptive optics imaging and integral field unit spectroscopy we report the discovery of an expanding and narrowly confined bipolar shell surrounding the helium nova V445 Puppis (Nova Puppis 2000). An equatorial dust disc obscures the nova remnant, and the outflow is characterised by a large polar outflow velocity of 6720 +/- 650 km/s and knots moving at even larger velocities of 8450 +/- 570 km/s. We derive an expansion parallax distance of 8.2 +/- 0.5 kpc and deduce a pre-outburst luminosity of the underlying binary of log L/L_Sun = 4.34 +/- 0.36. The derived luminosity suggests that V445 Puppis probably contains a massive white dwarf accreting at high rate from a helium star companion making it part of a population of binary stars that potentially lead to supernova Ia explosions due to accumulation of helium-rich material on the surface of a massive white dwarf.

ELLIPSOIDAL VARIABLE V1197 ORIONIS: ABSOLUTE LIGHT-VELOCITY ANALYSIS FOR KNOWN DISTANCE

V1197 Orionis light curves from a long-term observing program for red giant binaries show ellipsoidal variation of small amplitude in the V and RC bands, although not clearly in U and B. Eclipses are not detected. All four bands show large irregular intrinsic variations, including fleeting quasi-periodicities identified by power spectra, that degrade analysis and may be caused by dynamical tides generated by orbital eccentricity. To deal with the absence of eclipses and consequent lack of astrophysical and geometrical information, direct use is made of the Hipparcos parallax distance while the V and RC light curves and (older) radial velocity curves are analyzed simultaneously in terms of absolute flux. The red giant's temperature is estimated from new spectra. This type of analysis, called Inverse Distance Estimation for brevity, is new and can also be applied to other ellipsoidal variables. Advantages gained by utilization of definite distance and temperature are discussed in regard to how radius, fractional lobe filling, and mass ratio information are expressed in the observations. The advantages were tested in solutions of noisy synthetic data. Also discussed and tested by simulations are ideas on the optimal number of light curves to be solved simultaneously under various conditions. The dim companion has not been observed or discussed in the literature but most solutions find its mass to be well below that of the red giant. Solutions show red giant masses that are too low for evolution to the red giant stage within the age of the Galaxy, although that result is probably an artifact of the intrinsic brightness fluctuations.

PRECISION SOUTHERN HEMISPHERE VLBI PULSAR ASTROMETRY. II. MEASUREMENT OF SEVEN PARALLAXES

Accurate measurement of pulsar distances via astrometry using very long baseline interferometry enables the improvement of Galactic electron density distribution models, improving distance estimates for the vast majority of pulsars for which parallax measurements are unavailable. However, pulsars at southern declinations have been under-represented in previous interferometric astrometry campaigns. In order to redress this imbalance, we have conducted a two-year astrometric campaign targeting eight southern pulsars with the Australian Long Baseline Array. The program summarized in this paper has resulted in the measurement of seven new pulsar parallaxes, with success on objects down to a mean flux density of 0.8 mJy at 1600 MHz. Our results highlight the substantial uncertainties that remain when utilizing free electron density models for individual pulsar distances. Until this study, PSR J0630-2834 was believed to convert 16% of its spin-down energy into x-rays, but our measured parallax distance of 332 (+52 -40) pc has revised this value to <1%. In contrast, PSR J0108-1431 was found to be almost a factor of two more distant than previously thought, making its conversion of spin-down energy to x-rays the most efficient known (>1%). The 8.5 second radio pulsar J2144-3933 was found to be closer than previously predicted, making its apparent 1400 MHz radio luminosity the lowest of any known pulsar (20 microJy kpc^2). We have examined the growing population of neutron stars with accurate parallaxes to determine the effect of distance errors on the underlying neutron star velocity distribution, and find that typical distance errors may be biasing the estimated mean pulsar velocity upwards by 5%, and are likely to exaggerate the distribution's high-velocity tail.

PULSED GAMMA RAYS FROM THE MILLISECOND PULSAR J0030+0451 WITH THEFERMILARGE AREA TELESCOPE

We report the discovery of gamma-ray pulsations from the nearby isolated millisecond pulsar PSR J0030+0451 with the Large Area Telescope (LAT) on the \emph{Fermi} Gamma-ray Space Telescope (formerly GLAST). This discovery makes PSR J0030+0451 the second millisecond pulsar to be detected in gamma-rays after PSR J0218+4232, observed by the EGRET instrument on the Compton Gamma Ray Observatory. The spin-down power $\dot E = $ 3.5 $\times$ 10$^{33}$ ergs s$^{-1}$ is an order of magnitude lower than the empirical lower bound of previously known gamma-ray pulsars. The emission profile is characterized by two narrow peaks, respectively 0.07 $\pm$ 0.01 and 0.08 $\pm$ 0.02 wide, separated by 0.44 $\pm$ 0.02 in phase. The first gamma-ray peak falls 0.15 $\pm$ 0.01 after the main radio peak. The pulse shape is similar to that of the normal gamma-ray pulsars. An exponentially cut-off power-law fit of the emission spectrum leads to an integral photon flux above 100 MeV of (6.76 $\pm$ 1.05 $\pm$ 1.35) $\times 10^{-8}$ cm$^{-2}$ s$^{-1}$ with cut-off energy (1.7 $\pm$ 0.4 $\pm$ 0.5) GeV. Based on its parallax distance of $(300 \pm 90)$ pc, we obtain a gamma-ray efficiency $L_\gamma / \dot{E} \simeq 15%$ for the conversion of spin-down energy rate into gamma-ray radiation, assuming isotropic emission.

DISTANCES TO FOUR SOLAR NEIGHBORHOOD ECLIPSING BINARIES FROM ABSOLUTE FLUXES

Eclipsing binary (EB)-based distances are estimated for four solar neighborhood EBs by means of the Direct Distance Estimation (DDE) algorithm. Results are part of a project to map the solar neighborhood EBs in three dimensions, independently of parallaxes, and provide statistical comparisons between EB and parallax distances. Apart from judgments on adopted temperature and interstellar extinction, DDE's simultaneous light-velocity solutions are essentially objective and work as well for semidetached (SD) and overcontact binaries as for detached systems. Here, we analyze two detached and two SD binaries, all double lined. RS Chamaeleontis is a pre-main-sequence (MS), detached EB with weak ? Scuti variations. WW Aurigae is detached and uncomplicated, except for having high metallicity. RZ Cassiopeiae is SD and has very clear ? Scuti variations and several peculiarities. R Canis Majoris (R CMa) is an apparently simple but historically problematic SD system, also with weak ? Scuti variations. Discussions include solution rules and strategies, weighting, convergence, and third light problems. So far there is no indication of systematic band dependence among the derived distances, so the adopted band-calibration ratios seem consistent. Agreement of EB-based and parallax distances is typically within the overlapped uncertainties, with minor exceptions. We also suggest an explanation for the long-standing undermassiveness problem of R CMa's hotter component, in terms of a fortuitous combination of low metallicity and evolution slightly beyond the MS.

PRECISION ASTROMETRY WITH THE VERY LONG BASELINE ARRAY: PARALLAXES AND PROPER MOTIONS FOR 14 PULSARS

Astrometry can bring powerful constraints to bear on a variety of scientific questions about neutron stars, including their origins, astrophysics, evolution, and environments. Using phase-referenced observations at the VLBA, in conjunction with pulsar gating and in-beam calibration, we have measured the parallaxes and proper motions for 14 pulsars. The smallest measured parallax in our sample is 0.13+-0.02 mas for PSR B1541+09, which has a most probable distance of 7.2+1.3-1.1 kpc. We detail our methods, including initial VLA surveys to select candidates and find in-beam calibrators, VLBA phase-referencing, pulsar gating, calibration, and data reduction. The use of the bootstrap method to estimate astrometric uncertainties in the presence of unmodeled systematic errors is also described. Based on our new model-independent estimates for distance and transverse velocity, we investigate the kinematics and birth sites of the pulsars and revisit models of the Galactic electron density distribution. We find that young pulsars are moving away from the Galactic plane, as expected, and that age estimates from kinematics and pulsar spindown are generally in agreement, with certain notable exceptions. Given its present trajectory, the pulsar B2045-16 was plausibly born in the open cluster NGC 6604. For several high-latitude pulsars, the NE2001 electron density model underestimates the parallax distances by a factor of two, while in others the estimates agree with or are larger than the parallax distances, suggesting that the interstellar medium is irregular on relevant length scales. The VLBA astrometric results for the recycled pulsar J1713+0747 are consistent with two independent estimates from pulse timing, enabling a consistency check between the different reference frames.

TRIGONOMETRIC PARALLAXES OF MASSIVE STAR-FORMING REGIONS. I. S 252 &amp; G232.6+1.0

We are conducting a large program with the NRAO Very Long Baseline Array (VLBA) to measure trigonometric parallaxes of massive star-forming regions across the Milky Way. Here we report measurement of the parallax and proper motion of methanol masers in S 252 and G232.6+1.0. The parallax of S 252 is 0.476 +/- 0.006 mas (2.10 [+0.027/-0.026] kpc), placing it in the Perseus spiral arm. The parallax of G232.6+1.0 is 0.596 +/- 0.035 mas (1.68 [+0.11/-0.09] kpc), placing it between the Carina-Sagittarius and Perseus arms, possibly in a Local (Orion) spur of the Carina-Sagittarius arm. For both sources, kinematic distances are significantly greater than their parallax distances. Our parallaxes and proper motions yield full space motions accurate to about 1 km/s. Both sources orbit the Galaxy about 13 km/s slower than circular rotation.

Extremely High Precision VLBI Astrometry of PSR J0437-4715 and Implications for Theories of Gravity

Using the recently upgraded Long Baseline Array, we have measured the trigonometric parallax of PSR J0437-4715 to better than 1% precision, the most precise pulsar distance determination made to date. Comparing this VLBI distance measurement to the kinematic distance obtained from pulsar timing, which is calculated from the pulsar's proper motion and apparent rate of change of orbital period, gives a precise limit on the unmodeled relative acceleration between the Solar System and PSR J0437-4715, which can be used in a variety of applications. Firstly, it shows that Newton's gravitational constant G is stable with time (\dot{G}/G = (-5 +- 26) x 10^{-13} yr^{-1}, 95% confidence). Secondly, if a stochastic gravitational wave background existed at the currently quoted limit, this null result would fail ~50% of the time. Thirdly, it excludes Jupiter-mass planets within 226 AU of the Sun in 50% of the sky (95% confidence). Finally, the ~1% agreement of the parallax and orbital period derivative distances provides a fundamental confirmation of the parallax distance method upon which all astronomical distances are based.

Suzaku Detection of an Intense X-Ray Flare from an A-Type Star, HD161084

Abstract We report on a serendipitous detection of an intense X-ray flare from the Tycho reference source on HD161084 during a Suzaku observation of the galactic center region for $\sim$20ks. The X-ray Imaging Spectrometer recorded a flare from this A1-type dwarf or subgiant star with a flux of $\sim$1.4$\times$10$^{-12}$erg s$^{-1}$ cm$^{-2}$ (0.5-10keV) and a decay time scale of $\sim$0.5hr. The spectrum is hard with a prominent FeXXV K$\alpha$ emission line at 6.7keV, which is explained by a $\sim$5keV thin-thermal plasma model attenuated by a $\sim$1.4$\times$10$^{21}$cm$^{-2}$ extinction. The low extinction, which is consistent with the optical reddening, indicates that the source is a foreground star toward the galactic center region. Based on a spectroscopic parallax distance of $\sim$530pc, the peak X-ray luminosity amounts to $\sim$1$\times$10$^{32}$erg s$^{-1}$ (0.5-10keV). This is much larger than the X-ray luminosity of ordinary late-type main-sequence stars, and the X-ray emission is unattributable to a hidden late-type companion that comprises a wide binary system with the A star. We discuss possible nature of HD161084, and suggest that it is most likely an interacting binary with elevated magnetic activity in the companion, such as the Algol-type system. The flux detected by Suzaku during the burst is $\sim$100-times larger than the quiescent level measured using the archived XMM-Newton and Chandra data. The large flux amplification makes this star a unique example among sources of this class.