Search

We report the discovery by the HATSouth network of HATS-18 b: a 1.980 +/- 0.077 Mj, 1.337 +0.102 -0.049 Rj planet in a 0.8378 day orbit, around a solar analog star (mass 1.037 +/- 0.047 Msun, and radius 1.020 +0.057 -0.031 Rsun) with V=14.067 +/- 0.040 mag. The high planet mass, combined with its short orbital period, implies strong tidal coupling between the planetary orbit and the star. In fact, given its inferred age, HATS-18 shows evidence of significant tidal spin up, which together with WASP-19 (a very similar system) allows us to constrain the tidal quality factor for Sun-like stars to be in the range 6.5 <= lg(Q*/k_2) <= 7 even after allowing for extremely pessimistic model uncertainties. In addition, the HATS-18 system is among the best systems (and often the best system) for testing a multitude of star--planet interactions, be they gravitational, magnetic or radiative, as well as planet formation and migration theories., Comment: Submitted. 12 pages, 9 figures, 5 tables

We report the discovery and characterization of four transiting exoplanets by the HATNet survey. The planet HAT-P-50b has a mass of 1.35 M_J and a radius of 1.29 R_J, and orbits a bright (V = 11.8 mag) M = 1.27 M_sun, R = 1.70 R_sun star every P = 3.1220 days. The planet HAT-P-51b has a mass of 0.31 M_J and a radius of 1.29 R_J, and orbits a V = 13.4 mag, M = 0.98 M_sun, R = 1.04 R_sun star with a period of P = 4.2180 days. The planet HAT-P-52b has a mass of 0.82 M_J and a radius of 1.01 R_J, and orbits a V = 14.1 mag, M = 0.89 M_sun, R = 0.89 R_sun star with a period of P = 2.7536 days. The planet HAT-P-53b has a mass of 1.48 M_J and a radius of 1.32 R_J, and orbits a V = 13.7 mag, M = 1.09 M_sun, R = 1.21 R_sun star with a period of P = 1.9616 days. All four planets are consistent with having circular orbits and have masses and radii measured to better than 10% precision. The low stellar jitter and favorable R_P/R_star ratio for HAT-P-51 make it a promising target for measuring the Rossiter-McLaughlin effect for a Saturn-mass planet., Comment: Submitted to AJ. 20 pages, 9 figures, 5 tables. Data available at http://hatnet.org/

A number of planet-host stars have been observed to rotate with a period equal to an integer multiple of the orbital period of their close planet. We expand this list by analyzing Kepler data of HAT-P-11 and finding a period ratio of 6:1. In particular, we present evidence for a long-lived spot on the stellar surface that is eclipsed by the planet in the same position four times, every sixth transit. We also identify minima in the out-of-transit lightcurve and confirm that their phase with respect to the stellar rotation is mostly stationary for the 48-month timeframe of the observations, confirming the proposed rotation period. For comparison, we apply our methods to Kepler-17 and confirm the findings of Bonomo & Lanza (2012) that the period ratio is not exactly 8:1 in that system. Finally, we provide a hypothesis on how interactions between a star and its planet could possibly result in an observed commensurability for systems where the stellar differential rotation profile happens to include a period at some latitude which is commensurable to the planetary orbit., Comment: Accepted for publication in ApJ. 10 pages, 7 figures

We present HST STIS observations of two occultations of the transiting exoplanet HAT-P-1b. By measuring the planet to star flux ratio near opposition, we constrain the geometric albedo of the planet, which is strongly linked to its atmospheric temperature gradient. An advantage of HAT-P-1 as a target is its binary companion ADS 16402 A, which provides an excellent photometric reference, simplifying the usual steps in removing instrumental artifacts from HST time-series photometry. We find that without this reference star, we would need to detrend the lightcurve with the time of the exposures as well as the first three powers of HST orbital phase, and this would introduce a strong bias in the results for the albedo. However, with this reference star, we only need to detrend the data with the time of the exposures to achieve the same per-point scatter, therefore we can avoid most of the bias associated with detrending. Our final result is a 2 sigma upper limit of 0.64 for the geometric albedo of HAT-P-1b between 577 and 947 nm., Comment: 8 pages, 2 figures, 3 tables

We report the discovery of HAT-P-31b, a transiting exoplanet orbiting the V=11.660 dwarf star GSC 2099-00908. HAT-P-31b is the first HAT planet discovered without any follow-up photometry, demonstrating the feasibility of a new mode of operation for the HATNet project. The 2.17 Mj, 1.1Rj planet has a period P = 5.0054 days and maintains an unusually high eccentricity of e = 0.2450+/-0.0045, determined through Keck, FIES and Subaru high precision radial velocities. Detailed modeling of the radial velocities indicates an additional quadratic residual trend in the data detected to very high confidence. We interpret this trend as a long-period outer companion, HAT-P-31c, of minimum mass 3.4Mj and period >2.8 years. Since current RVs span less than half an orbital period, we are unable to determine the properties of HAT-P-31c to high confidence. However, dynamical simulations of two possible configurations show that orbital stability is to be expected. Further, if HAT-P-31c has non-zero eccentricity, our simulations show that the eccentricity of HAT-P-31b is actively driven by the presence of c, making HAT-P-31 a potentially intriguing dynamical laboratory., Comment: 12 pages, 7 figures, 7 tables. Accepted in AJ. Minor changes over previous version

We report the discovery of HAT-P-27b, an exoplanet transiting the moderately bright G8 dwarf star GSC 0333-00351 (V=12.214). The orbital period is 3.039586 +/- 0.000012 d, the reference epoch of transit is 2455186.01879 +/- 0.00054 (BJD), and the transit duration is 0.0705 +/- 0.0019 d. The host star with its effective temperature 5300 +/- 90 K is somewhat cooler than the Sun, and is more metal-rich with a metallicity of +0.29 +/- 0.10. Its mass is 0.94 +/- 0.04 Msun and radius is 0.90 +/- 0.04 Rsun. For the planetary companion we determine a mass of 0.660 +/- 0.033 MJ and radius of 1.038 +0.077 -0.058 RJ. For the 30 known transiting exoplanets between 0.3 MJ and 0.8 MJ, a negative correlation between host star metallicity and planetary radius, and an additional dependence of planetary radius on equilibrium temperature are confirmed at a high level of statistical significance., Comment: Submitted to ApJ on 2011-01-18. 12 pages, 7 figures, 7 tables

Using data from the HATNet survey for transiting exoplanets we measure photometric rotation periods for 368 Pleiades stars with 0.4 Msun < M < 1.3 Msun. We detect periodic variability for 74% of the cluster members in this mass range that are within our field-of-view, and 93% of the members with 0.7 Msun < M < 1.0 Msun. This increases, by a factor of five, the number of Pleiades members with measured periods. Included in our sample are 14 newly identified probable cluster members which have proper motions, photometry, and rotation periods consistent with membership. We compare this data to the rich sample of spectroscopically determined projected equatorial rotation velocities (vsini) available in the literature for this cluster. For stars with M > 0.85 Msun the rotation periods, vsini and radius estimates are consistent with the stars having an isotropic distribution of rotation axes, if a moderate differential rotation law is assumed. For stars with M < 0.85 Msun the inferred sini values are systematically larger than 1.0. These observations imply that the combination of measured parameters P(vsini)/R is too large by \sim 24% for low-mass stars in this cluster. By comparing our new mass-period relation for the Pleiades to the slightly older cluster M35, we confirm previous indications that the spin-down stalls at \sim 100 Myr for the slowest rotating stars with 0.7 Msun < M < 1.1 Msun a fact which may indicate that the internal transport of angular momentum is inefficient in slowly rotating solar mass stars., Comment: 17 pages, 15 figures, 4 tables. Accepted for publication in MNRAS. Data available at http://www.cfa.harvard.edu/~jhartman/Pleiades/catalog.tgz

We present refined parameters for the extrasolar planetary system HAT-P-2 (also known as HD 147506), based on new radial velocity and photometric data. HAT-P-2b is a transiting extrasolar planet that exhibits an eccentric orbit. We present a detailed analysis of the planetary and stellar parameters, yielding consistent results for the mass and radius of the star, better constraints on the orbital eccentricity, and refined planetary parameters. The improved parameters for the host star are M_star = 1.36 +/- 0.04 M_sun and R_star = 1.64 +/- 0.08 R_sun, while the planet has a mass of M_p = 9.09 +/- 0.24 M_Jup and radius of R_p = 1.16 +/- 0.08 R_Jup. The refined transit epoch and period for the planet are E = 2,454,387.49375 +/- 0.00074 (BJD) and P = 5.6334729 +/- 0.0000061 (days), and the orbital eccentricity and argument of periastron are e = 0.5171 +/- 0.0033 and omega = 185.22 +/- 0.95 degrees. These orbital elements allow us to predict the timings of secondary eclipses with a reasonable accuracy of ~15 minutes. We also discuss the effects of this significant eccentricity including the characterization of the asymmetry in the transit light curve. Simple formulae are presented for the above, and these, in turn, can be used to constrain the orbital eccentricity using purely photometric data. These will be particularly useful for very high precision, space-borne observations of transiting planets., Comment: Revised version, accepted for publication in MNRAS, 11 pages, 6 figures

We report the identification of 32 transiting-planet candidates in HATNet field G205. We describe the procedures that we have used to follow up these candidates with spectroscopic and photometric observations, and we present a status report on our interpretation of the 28 candidates for which we have follow-up observations. Eight are eclipsing binaries with orbital solutions whose periods are consistent with their photometric ephemerides; two of these spectroscopic orbits are singled-lined and six are double-lined. For one of the candidates, a nearby but fainter eclipsing binary proved to be the source for the HATNet light curve, due to blending in the HATNet images. Four of the candidates were found to be rotating more rapidly than vsini = 50 km/s and were not pursued further. Thirteen of the candidates showed no significant velocity variation at the level of 0.5 to 1.0 km/s . Seven of these were eventually withdrawn as photometric false alarms based on an independent reanalysis using more sophisticated tools. Of the remaining six, one was put aside because a close visual companion proved to be a spectroscopic binary, and two were not followed up because the host stars were judged to be too large. Two of the remaining candidates are members of a visual binary, one of which was previously confirmed as the first HATNet transiting planet, HAT-P-1b. In this paper we confirm that the last of this set of candidates is also a a transiting planet, which we designate HAT-P-8b, with mass Mp = 1.52 +/- 0.18/0.16 Mjup, radius Rp = 1.50 +/- 0.08/0.06 Rjup, and photometric period P = 3.076320 +/- 0.000004 days. HAT-P-8b has an inflated radius for its mass, and a large mass for its period. The host star is a solar-metallicity F dwarf, with mass M* = 1.28 +/- 0.04 Msun and Rp = 1.58 +/- 0.08/0.06 Rsun., Comment: 16 pages, 6 figures, 13 tables

The recently discovered transiting very hot Jupiter, HAT-P-7b, a planet detected by the telescopes of HATNet, turned out to be among the ones subjected to the highest irradiation from the parent star. As known, the combination of photometric and spectroscopic data for such an object yields the stellar, orbital and planetary parameters. In order to best characterize this particular planet, we carried out a complex analysis based on a complete and simultaneous Monte-Carlo solution using all available data. We included the discovery light curves, partial follow-up light curves, the radial velocity data, and we used the stellar evolution models to infer the stellar properties. This self-consistent way of modeling provides the most precise estimate of the a posteriori distributions of all of the system parameters of interest, and avoids making assumptions on the values and uncertainties of any of the internally derived variables describing the system. This analysis demonstrates that even partial light curve information can be valuable. This may become very important for future discoveries of planets with longer periods -- and therefore longer transit durations -- where the chance of observing a full event is small., Comment: 4 pages, 2 figures. To appear in the Proceedings of IAU Symposium 253, "Transiting Planets", May 2008, Cambridge, MA, USA

We report the discovery of a planet transiting a moderately faint (V=12.3 mag) late F star, with an orbital period of 3.92289 +/- 0.00004 days. From the transit light curve and radial velocity measurements we determine that the radius of the planet is R_p = 1.40 +/- 0.06 R_Jup and that the mass is M_p = 0.78 +/- 0.09 M_Jup. The density of the new planet, rho = 0.35 +/- 0.06 g cm^{-3}, fits to the low-density tail of the currently known transiting planets. We find that the center of transit is at T_c = 2454417.9077 +/- 0.0003 (HJD), and the total transit duration is 0.143 +/- 0.004 days. The host star has M_s = 1.28 +/- 0.13 M_Sun and R_s = 1.32 +/- 0.07 R_Sun., Comment: Submitted to ApJ; V2: Replaced with accepted version

We use multi-band photometry to refine estimates for the planetary radius and orbital inclination of the transiting planet system HD 209458. We gathered 1066 spectra over four distinct transits with the STIS spectrometer on the Hubble Space Telescope using two gratings with a resolution R=1500 and a combined wavelength range of 290-1030 nm. We divide the spectra into ten spectrophotometric bandpasses, five for each grating, of equal wavelength span within each grating, and fit a transit curve over all bandpasses simultaneously. In our fit we use theoretical values for the stellar limb-darkening to further constrain the planetary radius. We find that the radius of HD 209458b is 1.320 +/- 0.025 R_Jup, which is a factor of two more precise than current estimates. We also obtain improved estimates for the orbital period P and time of center of transit T_C. Although in principle the photon-limited precision of the STIS data should allow us to measure the timing of individual transits to a precision of 2-7 s, we find that systematic instrumental offsets in the measured flux from one orbit of the spacecraft to the next degrade these measurements to a typical precision of +/- 14 s. Within this level of error, we find no significant variations in the timing of the eight events examined in this work., Comment: final accepted version of paper

We show that the very close-by (19 pc) K0 star HD 189733, already found to be orbited by a transiting giant planet, is the primary of a double-star system, with the secondary being a mid-M dwarf with projected separation of about 216 AU from the primary. This conclusion is based on astrometry, proper motion and radial velocity measurements, spectral type determination and photometry. We also detect differential proper motion of the secondary. The data appear consistent with the secondary orbiting the primary in a clockwise orbit, lying nearly in the plane of the sky (that is, nearly perpendicular to the orbital plane of the transiting planet), and with period about 3200 years., Comment: 4 pages, submitted to APJ Letters

We report g, V, and r photometric time series of HD 149026 spanning predicted times of transit of the Saturn-mass planetary companion, which was recently discovered by Sato and collaborators. We present a joint analysis of our observations and the previously reported photometry and radial velocities of the central star. We refine the estimate of the transit ephemeris to Tc [HJD] = 2453527.87455^{+0.00085}_{-0.00091} + N * 2.87598^{+0.00012}_{-0.00017}. Assuming that the star has a radius of 1.45 +/- 0.10 R_Sun and a mass of 1.30 +/- 0.10 M_Sun, we estimate the planet radius to be 0.726 +/- 0.064 R_Jup, which implies a mean density of 1.07^{+0.42}_{-0.30} g/cm^3. This density is significantly greater than that predicted for models which include the effects of stellar insolation and for which the planet has only a small core of solid material. Thus we confirm that this planet likely contains a large core, and that the ratio of core mass to total planet mass is more akin to that of Uranus and Neptune than that of either Jupiter or Saturn., Comment: 18 pages, 4 figures, submitted to ApJ

We determine the stellar, planetary, and orbital properties of the transiting planetary system HD 209458, through a joint analysis of high-precision radial velocities, photometry, and timing of the secondary eclipse. Of primary interest is the strong detection of the Rossiter-McLaughlin effect, the alteration of photospheric line profiles that occurs because the planet occults part of the rotating surface of the star. We develop a new technique for modeling this effect, and use it to determine the inclination of the planetary orbit relative to the apparent stellar equator [lambda = (-4.4 +/- 1.4) degrees], and the line-of-sight rotation speed of the star [v*sin(I) = (4.70 +/- 0.16) km/s]. The uncertainty in these quantities has been reduced by an order of magnitude relative to the pioneering measurements by Queloz and collaborators. The small but nonzero misalignment is probably a relic of the planet formation epoch, because the expected timescale for tidal coplanarization is larger than the age of the star. Our determination of v*sin(I) is a rare case in which rotational line broadening has been isolated from other broadening mechanisms., Comment: 21 pages, 6 figures. To match the version accepted by ApJ. Errors have been corrected in Table 1. (In previous versions, some confidence limits were mistakenly 80% rather than 90% limits.)

We report high precision spectrophotometric observations of four planetary transits of HD 209458, in the region of the sodium resonance doublet at 589.3 nm. We find that the photometric dimming during transit in a bandpass centered on the sodium feature is deeper by (2.32 +/- 0.57) x 10^{-4} relative to simultaneous observations of the transit in adjacent bands. We interpret this additional dimming as absorption from sodium in the planetary atmosphere, as recently predicted from several theoretical modeling efforts. Our model for a cloudless planetary atmosphere with a solar abundance of sodium in atomic form predicts more sodium absorption than we observe. There are several possibilities that may account for this reduced amplitude, including reaction of atomic sodium into molecular gases and/or condensates, photoionization of sodium by the stellar flux, a low primordial abundance of sodium, or the presence of clouds high in the atmosphere., Comment: 26 pages, 8 figures, accepted by ApJ 2001 November 19

We have observed 4 transits of the planet of HD 209458 using the STIS spectrograph on HST. Summing the recorded counts over wavelength between 582 nm and 638 nm yields a photometric time series with 80 s time sampling and relative precision of about 1.1E-4 per sample. The folded light curve can be fit within observational errors using a model consisting of an opaque circular planet transiting a limb-darkened stellar disk. In this way we estimate the planetary radius R_p = 1.347 +/- 0.060 R_Jup, the orbital inclination i = 86.68 +/- 0.14 degrees, the stellar radius R_* = 1.146 +/- 0.050 R_solar, and one parameter describing the stellar limb darkening. Our estimated radius is smaller than those from earlier studies, but is consistent within measurement errors, and is also consistent with theoretical estimates of the radii of irradiated Jupiter-like planets. Satellites or rings orbiting the planet would, if large enough, be apparent from distortions of the light curve or from irregularities in the transit timings. We find no evidence for either satellites or rings, with upper limits on satellite radius and mass of 1.2 Earth radii and 3 Earth masses, respectively. Opaque rings, if present, must be smaller than 1.8 planetary radii in radial extent. The high level of photometric precision attained in this experiment confirms the feasibility of photometric detection of Earth-sized planets circling Sun-like stars., Comment: 28 pages, 9 figures. To be published in the Astrophysical Journal

Close-in extrasolar giant planets may be directly detectable by their reflected light, due to the proximity of the planet to the illuminating star. The spectrum of the system will contain a reflected light component that varies in amplitude and Doppler shift as the planet orbits the star. Intensive searches for this effect have been carried out for only one extrasolar planet system, tau Boo. There exist several other attractive targets, including the transiting planet system HD 209458., Comment: 6 pages, 2 figures, to appear in "Disks, Planetesimals, and Planets", ASP Conference Series, F. Garzon and T. J. Mahoney, eds

The planet orbiting tau Boo at a separation of 0.046 AU could produce a reflected light flux as bright as 1e-4 relative to that of the star. A spectrum of the system will contain a reflected light component which varies in amplitude and Doppler-shift as the planet orbits the star. Assuming the secondary spectrum is primarily the reflected stellar spectrum, we can limit the relative reflected light flux to be less than 5e-5. This implies an upper limit of 0.3 for the planetary geometric albedo near 480 nm, assuming a planetary radius of 1.2 R_Jup. This albedo is significantly less than that of any of the giant planets of the solar system, and is not consistent with certain published theoretical predictions., Comment: 5 pages, 1 figure, accepted by ApJ Letters

We discuss the interpretation of distortions to stellar spectral lines with particular attention to line bisectors in the presence of an orbiting planetary companion. We present a simple model whereby light reflected by the companion can cause temporal variations to the observed line profiles. These distortions have a characteristic signature depending on the inclination angle of the system. For the known close-in extrasolar giant planets, the expected amplitude of the effect might not be far from current detection capabilities. This method could be used to directly detect the presence of the companion, yielding the orbital inclination and hence the planetary mass. Futhermore, a detection would measure a combination of the planetary radius and albedo, from which a minimum radius may be deduced., Comment: 8 pages, 4 figures, uses emulateapj.sty, accepted by ApJ Letters

Because of our relatively low spectral resolution, we compare our observations with Gray's line bisector data by fitting observed line profiles to an expansion in terms of orthogonal (Hermite) functions. To obtain an accurate comparison, we model the emergent line profiles from rotating and pulsating stars, taking the instrumental point spread function into account. We describe this modeling process in detail. We find no evidence for line profile or strength variations at the radial velocity period in either 51 Peg or in Tau Boo. For 51 Peg, our upper limit for line shape variations with 4.23-day periodicity is small enough to exclude with 10 sigma confidence the bisector curvature signal reported by Gray & Hatzes; the bisector span and relative line depth signals reported by Gray (1997) are also not seen, but in this case with marginal (2 sigma) confidence. We cannot, however, exclude pulsations as the source of 51 Peg's radial velocity variation, because our models imply that line shape variations associated with pulsations should be much smaller than those computed by Gray & Hatzes; these smaller signals are below the detection limits both for Gray & Hatzes' data and for our own. Tau Boo's large radial velocity amplitude and v*sin(i) make it easier to test for pulsations in this star. Again we find no evidence for periodic line-shape changes, at a level that rules out pulsations as the source of the radial velocity variability. We conclude that the planet hypothesis remains the most likely explanation for the existing data., Comment: 44 pages, 19 figures, plain TeX, accepted to ApJS (companion to letter astro-ph/9712279)

We develop a formulation of the Baade-Wesselink method which uses the Fourier coefficients of the observables. We derive an explicit, analytic expression to determine the mean radius from each Fourier order. The simplicity of this method allows us to derive the uncertainty in the mean radius due to measurement errors. Using simulations and a recent dataset we demonstrate that the precision of the radius measurement with optical magnitudes is in most cases limited by the accuracy of the measurement of the phase difference between the light and the color index curve. In this case it is advantageous to determine the inverse radius, because it has normal errors., Comment: 18 pages, postscript, accepted for publication in ApJ

This is the final report for this NASA grant. Work under the first three years of the Grant (from May 1 2001 through April 30 2004) has been described in previous annual reports. Here we briefly summarize that work, and then focus on activities between May 1 2004 (the start of a 1-year no-cost extension period) and April 30 2005, the end of the Grant period.

The upgrade to the AFOE spectrograph at the 1.5-m Tillinghast telescope at FLWO Observatory, described in detail in the previous progress report, was completed. Work included installing a new CCD detector (using SA0 internal funds) to replace the previous detector, which had failed after many thermal cyclings during the previous nearly ten years of operation. The new detector (a back-illuminated Marconi 2K x 2K detector) shows excellent quantum efficiency, and employs dual readout channels to keep the readout time as short as 20 seconds. Tests of the attainable radial velocity precision of the upgraded spectrograph with the new detector show an improvement in radial velocity detection efficiency by a factor of about 10 over the previous AFOE configuration, allowing us to push the photon-noise-limited doppler precision of 7 m/s to correspondingly fainter stars. This work was carried out by Co-investigator Sylvain Korzennik, along with the late Peter Nisenson.

We initiated a significant hardware upgrade to the AFOE, to increase its efficiency for precise radial velocity studies to the level where we can continue to contribute usefully to extrasolar planet research on relatively bright stars. The AFOE, at a 1.5-m telescope, will of course not have the sensitivity of radial velocity instruments at larger telescopes, such as the HIRES on Keck or the Hectochelle on the MMT telescope (about to come on line). However, it has been possible to increase its efficiency for precise radial velocity studies by a factor of 4 to 5, which-combined with the large amount of telescope time available at the 1.5-m telescope-will permit us to do intensive follow-up observations of stars brighter than about 8 magnitude. The AFOE was originally designed primarily for asteroseismology using a ThAr reference. This provided useful wavelength stability over tens of minutes as required for asteroseismology, but we were unable to get a long-term (month-to-month) velocity precision better than about 15 m/s with that setup. Hence, we implemented an iodine cell as a wavelength reference for extrasolar planet studies. However, the optical design of the original AFOE did not completely span the wavelength range covered by the iodine absorption spectrum, and furthermore the optics suffered significant light loss through optical obscuration in the camera secondary. To remedy this, we replaced the AFOE grating with a new one that covered the entire iodine spectral range at somewhat lower spectral resolution, and replaced the camera with a transmitting lens. (The use of a lens was made possible by restricting the spectral range covered by the upgraded AFOE to only the iodine region.) These upgrades were successfully completed, and the instrument was tested for three nights in fall of 2002. The expected improvement in sensitivity by a factor of 4 to 5 was observed: that is, the same velocity precision as previously attained (of order 5 to 7 m/s) was now obtained on a stars that are about 4 to 5 times fainter. Unfortunately, just after installing the new AFOE optics, there was a catastrophic failure of the old Tektronix CCD. A careful investigation was made, and it was determined that the failure was probably due to too many thermal cyclings over the previous 8 years. The only solution was to purchase a new CCD. Internal funds of $16,000 were found for this purpose, and the new CCD (a thinned, back-illuminated Marconi 2K x 2K device with 13.5 micron pixels, and employing dual readout channels to keep the readout time as short as 20 seconds) was ordered, has now arrived, and is currently under test before 1 being installed in the AFOE liquid nitrogen dewar. At the same time, we expect to implement new CCD control electronics. This system should be better in many ways than the previous detector (including eliminating the need for UV flashing to maximize the quantum efficiency). We plan to install the detector on the telescope and begin operations again in April of this year.

This is the final report of NASA Grant NAG5-7505, for 'Study of Extra-solar Planets with the Advanced Fiber Optic Echelle'. This program was funded in response to our proposal submitted under NASA NRA 97-OSS-06, with a total period of performance from June 1, 1998 through Feb 28 2002. Principal Investigator is Robert W. Noyes; co-Investigators are Sylvain G. Korzennik (SAO), Peter Niserison (SAO), and Timothy M. Brown (High Altitude Observatory). Since the start of this program we have carried out more than 30 observing runs, typically of 5 to 7 days duration. We obtained a total of around 2000 usable observations of about 150 stars, where a typical observation consists of 3 exposures of 10 minutes each. Using this data base we detected thc two additional planetary companions to the star Upsilon Andromedae. This detection was made independently of, and essentially simultaneously with, a similar detection by the Berkeley group (Marcy et al): the fact that two data sets were completely independent and gave essentially the same orbital parameters for this three-planet system gave a strong confirmation of this important result. We also extended our previous detection of the planet orbiting Rho Coronae Borealis to get a better determination of its orbital eccentricity: e=0.13 +/- 0.05. We detected a new planet in orbit around the star HD 89744, with orbital period 256 days, semi-major axis 0.88 AU, eccentricity 0.70, and minimum mass m sini = 7.2 m(sub Jup). This discovery is significant because of the very high orbital eccentricity, arid also because HD 89744 has both high metallicity [Fe/H] and at the same time a low [C/Fe] abundance ratio.

This grant report highlights activity in the following areas: (1) Improvement in Precise Radial Velocity (PRV) analysis code; (2) Reanalysis of previous data; (3) Improvements to the AFOE (Advanced Fiber Optic Echelle) spectrograph; (4) Development of PRV capabilities for the Hectochelle; (5) Extra-solar planet studies; (6) Longer-term plans for the AFOE; (7) Completion and publication of the analysis of the transiting gas-giant planet HD 209458b.

This grant funded an observational and theoretical program to study the structure and dynamics of the solar photosphere and low chromosphere, and the spectral signatures that result. The overall goal is to learn about mechanisms that cause heating of the overlying atmosphere, and produce variability of solar emission in spectral regions important for astrophysics and space physics. The program exploited two new ground-based observational capabilities: one using the Swedish Solar Telescope on La Palma for very high angular resolution observations of the photospheric intensity field (granulation) and proxies of the magnetic field (G-band images); and the other using the Near Infrared Magnetograph at the McMath-Pierce Solar Facility to map the spatial variation and dynamic behavior of the solar temperature minimum region using infrared CO lines. We have interpreted these data using a variety of theoretical and modelling approaches, some developed especially for this project. Previous annual reports cover the work done up to 31 May 1997. This final report summarizes our work for the entire period, including the period of no-cost extension from 1 June 1997 through September 30 1997. In Section 2 we discuss observations and modelling of the photospheric flowfields and their consequences for heating of the overlying atmosphere, and in Section 3 we discuss imaging spectroscopy of the CO lines at 4.67 mu.

We continued our program of CO observations with the McMath-Pierce facility at Kitt Peak National Solar Observatory. Uitenbroek has developed a two-and three dimensional radiative transfer code that now includes chemical equilibrium calculations. This code allows us to compute a CO spectrum from for instance a snapshot of a solar granulation simulation (e.g. Stein & Nordlund 1989, Apj 342, L95) and compare these theoretical spectra with our spatially resolved CO spectroscopy. Van Ballegooijen and Uitenbroek have started calculations of two-dimensional fluxtube models that account consistently for hydrogen ionization in the calculation of the electron density. To this end we solve radiative transfer for hydrogen (bound-bound and bound-free transitions) in the two-dimensional models, including the effect of partial frequency redistribution (PRD) in the Lyman (alpha) and (beta) lines. From our internally consistent models we will calculate emergent spectra and the way these vary with location of some well-known spectral diagnostics and compare our results with observed line profiles. We can readily compare theoretical CO profiles from out models with our spatially resolved CO observations. Also we can compare with spatially resolved Ca II (ground based) and Mg II (we have observations done with the UVSP/SMM instrument), and Lyman (alpha) observations that should be available from SUMER/SOHO.

This progress report covers the first year of NASA Grant NAGw-2545, a study of magnetic structure in the solar photosphere and chromosphere. We have made significant progress in three areas: (1) analysis of vorticity in photospheric convection, which probably affects solar atmospheric heating through the stresses it imposes on photospheric magnetic fields; (2) modelling of the horizontal motions of magnetic footpoints in the solar photosphere using an assumed relation between brightness and vertical motion as well as continuity of flow; and (3) observations and analysis of infrared CO lines formed near the solar temperature minimum, whose structure and dynamics also yield important clues to the nature of heating of the upper atmosphere. Each of these areas are summarized in this report, with copies of those papers prepared or published this year included.

We have studied an outstanding sequence of continuum images of the solar granulation from Pic du Midi Observatory. We have calculated the horizontal vector flow field using a correlation tracking algorithm, and from this determined three scalar field: the vertical component of the curl; the horizontal divergence; and the horizontal flow speed. The divergence field has substantially longer coherence time and more power than does the curl field. Statistically, curl is better correlated with regions of negative divergence - that is, the vertical vorticity is higher in downflow regions, suggesting excess vorticity in intergranular lanes. The average value of the divergence is largest (i.e., outflow is largest) where the horizontal speed is large; we associate these regions with exploding granules. A numerical simulation of general convection also shows similar statistical differences between curl and divergence. Some individual small bright points in the granulation pattern show large local vorticities.

We discuss recent observations of the fundamental vibration-rotation transitions of carbon monoxide (CO) in the solar infrared spectrum. Employing a new array detector at the McMath-Pierce facility on Kitt Peak we find that the CO lines sketch a rich picture of the dynamics of the solar temperature minimum region, the lower boundary of the chromosphere. In a spectra-spectroheliogram and a time-sequence of the slit-spectra obtained during exceptional seeing conditions we observe small-scale bright, ring shaped, blueshifted features. We speculate that they are the signature of granular overshoot into the convectively stable temperature minimum. The centers of the rings are among the coolest elements seen in strong CO-line heliograms on the disk, and may be instrumental to the low temperature observed in CO close to the solar limb.

This is the final report on a project to develop a theoretical basis for interpreting solar oscillation data in terms of the interior dynamics and structure of the Sun. The topics covered include the following: (1) studies of the helioseismic signatures of differential rotation and convection in the solar interior; (2) wave generation by turbulent convection; and (3) the study of antipodal sunspot imaging of an active region tomography.