Search SciRate

We have used solar oscillation frequencies and frequency splittings obtained over solar cycles 23, 24 and the rising phase of solar cycle 25 to investigate whether the tachocline properties (jump i.e., the change in the rotation rate across the tachocline, width and position) show any time variation. We confirm that the change in rotation rate across the tachocline changes substantially, however, the change does not show a simple correlation with solar cycle unlike, for instance, changes in mode frequencies. The change during the ascending phase of solar cycle 25 is almost a mirror image of the change during the descending part of solar cycle 24, tempting us to speculate that the tachocline has a much longer period than either the sunspot or the magnetic cycle. We also find that the position of the tachocline, defined as the mid-point of the change in rotation rate, showed significant changes during solar cycle 24. The width of the tachocline, on the other hand, has showed significant changes during solar cycle 23, but not later. The change in the tachocline becomes more visible if we look at the upper and lower extents of the tachocline, defined as (position +/- width). We find that for epochs around solar maxima and minima, the extent decreases before increasing again - a few more years of data should clarify this trend. Our results reinforce the need to continue helioseismic monitoring of the Sun to understand solar activity and its evolution.

We present a new iterative rotation inversion technique based on the Simultaneous Algebraic Reconstruction Technique developed for image reconstruction. We describe in detail our algorithmic implementation and compare it to the classical inversion techniques like the Regularized Least Squares (RLS) and the Optimally Localized Averages (OLA) methods. In our implementation, we are able to estimate the formal uncertainty on the inferred solution using standard error propagation, and derive the averaging kernels without recourse to any Monte-Carlo simulation. We present the potential of this new technique using simulated rotational frequency splittings. We use noiseless sets that cover the range of observed modes and associate to these artificial splittings observational uncertainties. We also add random noise to present the noise magnification immunity of the method. Since the technique is iterative we also show its potential when using an apriori solution. With the right regularization this new method can outperform our RLS implementation in precision, scope and resolution. Since it results in very different averaging kernels where the solution is poorly constrained, this technique infers different values. Adding such a technique to our compendium of inversion methods will allow us to improve the robustness of our inferences when inverting real observations and better understand where they might be biased and/or unreliable, as we push our techniques to maximize the diagnostic potential of our observations.

We have used solar oscillation frequencies and frequency splittings obtained over solar cycles 23 and 24 to investigate whether the base of the solar convection zone shows any departure from spherical symmetry. We used the even-order splitting coefficients, $a_2$-$a_8$, and estimated the contributions from each one separately. The average asphericity over the two solar cycles was determined using frequencies and splittings obtained with a 9216-day time-series. We find that evidence of asphericity is, \em at best, marginal: the $a_2$ component is consistent with no asphericity, the $a_4$ and $a_6$ components yield results at a level a little greater than $1\,\sigma$, while the $a_8$ component shows a signature below $1\,\sigma$. The combined results indicate that the time average of the departure from the spherically symmetric position of the base of the convection zone is $\lesssim 0.0001R_\odot$. We have also used helioseismic data obtained from time-series of lengths 360 days, 576 days, 1152 days, and 2304 days in order to examine the consistency of the results and evaluate whether there is any time variation. We find that the evidence for time variation is statistically marginal in all cases, except for the $a_6$ component, for which tests consistently yield $p$ values of less than $0.05$.

The HMI project recently started processing the continuum intensity images following global helioseismology procedures similar to those used to process the velocity images. The spatial decomposition of these images has produced time series of spherical harmonic coefficients for degrees up to $\ell=300$, using a different apodization than the one used for velocity observations. The first 360 days of observations were processed and made available. In this paper I present initial results from fitting these time series using my state of the art fitting methodology and compare the derived mode characteristics to those estimated using co-eval velocity observations.

We present the first accurate characterization of high-degree modes, derived using the best MDI full-disk full-resolution data set available. A ninety day long time series of full-disk two arc-second per pixel resolution dopplergrams was acquired in 2001. These dopplergrams were spatially decomposed using our best estimate of the image scale and the known components of MDI's image distortion. A multi-taper power spectrum estimator was used to generate power spectra up to l = 1000, with a large number of tapers to reduce the realization noise, while the blending at high degrees negates the need for high spectral resolution. These power spectra were fitted for all degrees and all azimuthal orders, between l = 100 and l = 1000. This fitting generated in excess of 6x10^6 individual estimates of ridge frequencies, line-widths, amplitudes and asymmetries (singlets: l,n,m), corresponding to some 6,000 multiplets (l,n). Fitting at high degrees generates ridge characteristics, characteristics that do not correspond to the underlying mode characteristics. We used a sophisticated forward modeling to recover the best possible estimate of the underlying mode characteristics (mode frequencies, as well as line-widths, amplitudes and asymmetries). We describe in detail this modeling and its validation. The modeling has been extensively reviewed and refined, by including an iterative process to improve its input parameters to better match the observations. Also, the contribution of the leakage matrix on the accuracy of the procedure has been carefully assessed. We present the derived set of corrected mode characteristics, discuss their uncertainties and the precision of the ridge to mode correction schemes. In our conclusions, we address how to further improve these estimates, and the implications for other data sets, like GONG+ and HMI.

The dynamics of the solar radiative interior are still poorly constrained by comparison to the convective zone. This disparity is even more marked when we attempt to derive meaningful temporal variations. Many data sets contain a small number of modes that are sensitive to the inner layers of the Sun, but we found that the estimates of their uncertainties are often inaccurate. As a result, these data sets allow us to obtain, at best, a low resolution estimate of the solar core rotation rate down to approximately 0.2R. We present inferences based on mode determination resulting from an alternate peak-fitting methodology aimed at increasing the amount of observed modes that are sensitive to the radiative zone, while special care was taken in the determination of their uncertainties. This methodology has been applied to MDI and GONG data, for the whole Solar Cycle 23, and to the newly available HMI data. The numerical inversions of all these data sets result in the best inferences to date of the rotation in the radiative region. These results and the method used to obtain them are discussed. The resulting profiles are shown and analyzed, and the significance of the detected changes discussed.

Improved wavelength calibrators for high-resolution astrophysical spectrographs will be essential for precision radial velocity (RV) detection of Earth-like exoplanets and direct observation of cosmological deceleration. The astro-comb is a combination of an octave-spanning femtosecond laser frequency comb and a Fabry-Pérot cavity used to achieve calibrator line spacings that can be resolved by an astrophysical spectrograph. Systematic spectral shifts associated with the cavity can be 0.1-1 MHz, corresponding to RV errors of 10-100 cm/s, due to the dispersive properties of the cavity mirrors over broad spectral widths. Although these systematic shifts are very stable, their correction is crucial to high accuracy astrophysical spectroscopy. Here, we demonstrate an \emphin-situ technique to determine the systematic shifts of astro-comb lines due to finite Fabry-Pérot cavity dispersion. The technique is practical for implementation at a telescope-based spectrograph to enable wavelength calibration accuracy better than 10 cm/s.

We present preliminary asteroseismic results from Kepler on three G-type stars. The observations, made at one-minute cadence during the first 33.5d of science operations, reveal high signal-to-noise solar-like oscillation spectra in all three stars: About 20 modes of oscillation may be clearly distinguished in each star. We discuss the appearance of the oscillation spectra, use the frequencies and frequency separations to provide first results on the radii, masses and ages of the stars, and comment in the light of these results on prospects for inference on other solar-type stars that Kepler will observe.

We present an analysis of three years of precision radial velocity measurements of 160 metal-poor stars observed with HIRES on the Keck 1 telescope. We report on variability and long-term velocity trends for each star in our sample. We identify several long-term, low-amplitude radial-velocity variables worthy of follow-up with direct imaging techniques. We place lower limits on the detectable companion mass as a function of orbital period. Our survey would have detected, with a 99.5% confidence level, over 95% of all companions on low-eccentricity orbits with velocity semi-amplitude K > 100 m/s, or M_p*sin(i) > 3.0 M_JUP*(P/yr)^(1/3), for orbital periods P< 3 yr. None of the stars in our sample exhibits radial-velocity variations compatible with the presence of Jovian planets with periods shorter than the survey duration. The resulting average frequency of gas giants orbiting metal-poor dwarfs with -2.0 < [Fe/H] < -0.6 is f_p<0.67% (at the 1-sigma confidence level). We examine the implications of this null result in the context of the observed correlation between the rate of occurrence of giant planets and the metallicity of their main-sequence solar-type stellar hosts. By combining our dataset with the Fischer & Valenti (2005) uniform sample, we confirm that the likelihood of a star to harbor a planet more massive than Jupiter within 2 AU is a steeply rising function of the host's metallicity. However, the data for stars with -1.0 < [Fe/H] < 0.0 are compatible, in a statistical sense, with a constant occurrence rate f_p~1%. Our results can usefully inform theoretical studies of the process of giant planet formation across two orders of magnitude in metallicity.

We used the m-averaged spectrum technique ("collapsogram") to extract the low-frequency solar p-mode parameters of low- and intermediate-angular degrees (l ≤35) in long time series of GONG and MDI observations. Rotational splittings and central frequencies have been measured down to ~850\muHz, including predicted modes which have not been measured previously. Both GONG and MDI frequency splitting data sets were numerically inverted to extract the internal solar rotation rate. The impact of the very low-frequency observables and the differences between GONG and MDI data sets on the inversion results are also analyzed.

The detection of the signature of dipole gravity modes has opened the path to study the solar inner radiative zone. Indeed, g modes should be the best probes to infer the properties of the solar nuclear core that represents more than half of the total mass of the Sun. Concerning the dynamics of the solar core, we can study how future observations of individual g modes could enhance our knowledge of the rotation profile of the deep radiative zone. Applying inversions on a set of real p-mode splittings coupled with either one or several g modes, we have checked the improvement of the inferred rotation profile when different error bars are considered for the g modes. Moreover, using a new methodology based on the analysis of the almost constant separation of the dipole gravity modes, we can introduce new constraints on solar models. For that purpose, we can compare g-mode predictions computed from several models including different physical inputs with the g-mode asymptotic signature detected in GOLF data and calculate the correlation. This work shows the great consistency between the signature of dipole gravity modes and our knowledge of p-modes: incompatibility of data with a present standard model including the Asplund composition

We report new spectroscopic and photometric observations of the parent stars of the recently discovered transiting planets TrES-3 and TrES-4. A detailed abundance analysis based on high-resolution spectra yields [Fe/H] $= -0.19\pm 0.08$, $T_\mathrm{eff} = 5650\pm 75$ K, and $\log g = 4.4\pm 0.1$ for TrES-3, and [Fe/H] $= +0.14\pm 0.09$, $T_\mathrm{eff} = 6200\pm 75$ K, and $\log g = 4.0\pm0.1$ for TrES-4. The accuracy of the effective temperatures is supported by a number of independent consistency checks. The spectroscopic orbital solution for TrES-3 is improved with our new radial-velocity measurements of that system, as are the light-curve parameters for both systems based on newly acquired photometry for TrES-3 and a reanalysis of existing photometry for TrES-4. We have redetermined the stellar parameters taking advantage of the strong constraint provided by the light curves in the form of the normalized separation $a/R_\star$ (related to the stellar density) in conjunction with our new temperatures and metallicities. The masses and radii we derive are $M_\star=0.928_{-0.048}^{+0.028} M_{\sun}$,$R_\star = 0.829_{-0.022}^{+0.015} R_{\sun}$, and $M_\star = 1.404_{-0.134}^{+0.066} M_{\sun}$, $R_\star=1.846_{-0.087}^{+0.096} R_{\sun}$ for TrES-3 and TrES-4, respectively. With these revised stellar parameters we obtain improved values for the planetary masses and radii. We find $M_p = 1.910_{-0.080}^{+0.075} M_\mathrm{Jup}$, $R_p=1.336_{-0.036}^{+0.031} R_\mathrm{Jup}$ for TrES-3, and $M_p=0.925 \pm 0.082 M_\mathrm{Jup}$, $R_p=1.783_{-0.086}^{+0.093} R_\mathrm{Jup}$ for TrES-4. We confirm TrES-4 as the planet with the largest radius among the currently known transiting hot Jupiters.

Here we present a detailed analysis of solar acoustic mode frequencies and their rotational splittings for modes with degree up to 900. They were obtained by applying spherical harmonic decomposition to full-disk solar images observed by the Michelson Doppler Imager onboard the Solar and Heliospheric Observatory spacecraft. Global helioseismology analysis of high-degree modes is complicated by the fact that the individual modes cannot be isolated, which has limited so far the use of high-degree data for structure inversion of the near-surface layers (r > 0.97 R). In this work, we took great care to recover the actual mode characteristics using a physically motivated model which included a complete leakage matrix. We included in our analysis the following instrumental characteristics: the correct instantaneous image scale, the radial and non-radial image distortions, the effective position angle of the solar rotation axis and a correction to the Carrington elements. We also present variations of the mode frequencies caused by the solar activity cycle. We have analyzed seven observational periods from 1999 to 2005 and correlated their frequency shift with four different solar indices. The frequency shift scaled by the relative mode inertia is a function of frequency alone and follows a simple power law, where the exponent obtained for the p modes is twice the value obtained for the f modes. The different solar indices present the same result.

Accurate determination of the rotation rate in the radiative zone of the sun from helioseismic observations requires rotational frequency splittings of exceptional quality as well as reliable inversion techniques. We present here inferences based on mode parameters calculated from 2088-days long MDI, GONG and GOLF time series that were fitted to estimate very low frequency rotational splittings (nu < 1.7 mHz). These low frequency modes provide data of exceptional quality, since the width of the mode peaks is much smaller than the rotational splitting and hence it is much easier to separate the rotational splittings from the effects caused by the finite lifetime and the stochastic excitation of the modes. We also have implemented a new inversion methodology that allows us to infer the rotation rate of the radiative interior from mode sets that span l=1 to 25. Our results are compatible with the sun rotating like a rigid solid in most of the radiative zone and slowing down in the core (R_sun < 0.2). A resolution analysis of the inversion was carried out for the solar rotation inverse problem. This analysis effectively establishes a direct relationship between the mode set included in the inversion and the sensitivity and information content of the resulting inferences. We show that such an approach allows us to determine the effect of adding low frequency and low degree p-modes, high frequency and low degree p-modes, as well as some g-modes on the derived rotation rate in the solar radiative zone, and in particular the solar core. We conclude that the level of uncertainties that is needed to infer the dynamical conditions in the core when only p-modes are included is unlikely to be reached in the near future, and hence sustained efforts are needed towards the detection and characterization of g-modes.

The solar rotation profile is well constrained down to about 0.25 R thanks to the study of acoustic modes. Since the radius of the inner turning point of a resonant acoustic mode is inversely proportional to the ratio of its frequency to its degree, only the low-degree p modes reach the core. The higher the order of these modes, the deeper they penetrate into the Sun and thus they carry more diagnostic information on the inner regions. Unfortunately, the estimates of frequency splittings at high frequency from Sun-as-a-star measurements have higher observational errors due to mode blending, resulting in weaker constraints on the rotation profile in the inner core. Therefore inversions for the solar internal rotation use only modes below 2.4 mHz for l < 4. In the work presented here, we used an 11.5 year-long time series to compute the rotational frequency splittings for modes l < 4 using velocities measured with the GOLF instrument. We carried out a theoretical study of the influence of the low-degree modes in the region 2 to 3.5 mHz on the inferred rotation profile as a function of their error bars.

Using full-disk observations obtained with the Michelson Doppler Imager (MDI) on board the Solar and Heliospheric Observatory (SOHO) spacecraft, we present variations of the solar acoustic mode frequencies caused by the solar activity cycle. High-degree (100 < l < 900) solar acoustic modes were analyzed using global helioseismology analysis techniques over most of solar cycle 23. We followed the methodology described in details in Korzennik, Rabello-Soares and Schou (2004) to infer unbiased estimates of high-degree mode parameters (see also Rabello-Soares, Korzennik and Schou, 2006). We have removed most of the known instrumental and observational effects that affect specifically high-degree modes. We show that the high-degree changes are in good agreement with the medium-degree results, except for years when the instrument was highly defocused. We analyzed and discuss the effect of defocusing on high degree estimation. Our results for high-degree modes confirm that the frequency shift scaled by the relative mode inertia is a function of frequency and it is independent of degree.

We describe a high-precision Doppler search for giant planets orbiting a well-defined sample of metal-poor dwarfs in the field. This experiment constitutes a fundamental test of theoretical predictions which will help discriminate between proposed giant planet formation and migration models. We present here details on the survey as well as an overall assessment of the quality of our measurements, making use of the results for the stars that show no significant velocity variation.

I describe and present the results of a newly developed fitting methodology optimized for very long time series. The development of this new methodology was motivated by the fact that we now have more than half a decade of nearly uninterrupted observations by GONG and MDI, with fill factors as high as 89.8% and 82.2% respectively. It was recently prompted by the availability of a 2088-day-long time series of spherical harmonic coefficients produced by the MDI team. The fitting procedure uses an optimal sine-multi-taper spectral estimator -- whith the number of tapers based on the mode linewidth, the complete leakage matrix (i.e., horizontal as well as vertical components), and an asymmetric mode profile to fit simultaneously all the azimuthal orders with individually parametrized profiles. This method was applied to 2088-day-long time series of MDI and GONG observations, as well as 728-day-long subsets, and for spherical harmonic degrees between 1 and 25. The values resulting from these fits are inter-compared (MDI versus GONG) and compared to equivalent estimates from the MDI team and the GONG project. I also compare the results from fitting the 728-day-long subsets to the result of the 2088-day-long time series. This comparison shows the well known change of frequencies with solar activity -- and how it scales with a nearly constant pattern in frequency and m/l. This comparison also shows some changes in the mode linewidth and the constancy of the mode asymmetry.

As a part of an on-going program to explore the signature of p-modes in solar-like stars by means of high-resolution absorption lines pectroscopy, we have studied four stars (alfaCMi, etaCas A, zetaHer A and betaVir). We present here new results from two-site observations of Procyon A acquired over twelve nights in 1999. Oscillation frequencies for l=1 and l=0 (or 2) p-modes are detected in the power spectra of these Doppler shift measurements. A frequency analysis points out the dificulties of the classical asymptotic theory in representing the p-mode spectrum of Procyon A.

We present radial-velocity measurements obtained with the ELODIE and AFOE spectrographs for GJ 777 A (HD 190360), a metal-rich ([Fe/H]=0.25) nearby (d=15.9 pc) star in a stellar binary system. A long-period low radial-velocity amplitude variation is detected revealing the presence of a Jovian planetary companion. Some of the orbital elements remain weakly constrained because of the smallness of the signal compared to our instrumental precision. The detailed orbital shape is therefore not well established. We present our best fitted orbital solution: an eccentric (e=0.48) 10.7--year orbit. The minimum mass of the companion is 1.33 M_Jup.

ESPI has been proposed for direct imaging and spectral analysis of giant planets orbiting solar-type stars. ESPI extends the concept suggested by Nisenson and Papaliolios (2001) for a square aperture apodized telescope that has sufficient dynamic range to directly detect exo-planets. With a 1.5 M square mirror, ESPI can deliver high dynamic range imagery as close as 0.3 arcseconds to bright sources, permitting a sensitive search for exoplanets around nearby stars and a study of their characteristics in reflected light.

The characteristic of the solar acoustic spectrum is such that mode lifetimes get shorter and spatial leaks get closer in frequency as the degree of a mode increases for a given order. A direct consequence of this property is that individual p-modes are only resolved at low and intermediate degrees, and that at high degrees, individual modes blend into ridges. Once modes have blended into ridges, the power distribution of the ridge defines the ridge central frequency and it will mask the true underlying mode frequency. An accurate model of the amplitude of the peaks that contribute to the ridge power distribution is needed to recover the underlying mode frequency from fitting the ridge. We present the results of fitting high degree power ridges (up to l = 900) computed from several two to three-month-long time-series of full-disk observations taken with the Michelson Doppler Imager (MDI) on-board the Solar and Heliospheric Observatory between 1996 and 1999. We also present a detailed discussion of the modeling of the ridge power distribution, and the contribution of the various observational and instrumental effects on the spatial leakage, in the context of the MDI instrument. We have constructed a physically motivated model (rather than some ad hoc correction scheme) resulting in a methodology that can produce an unbiased determination of high-degree modes, once the instrumental characteristics are well understood. Finally, we present changes in high degree mode parameters with epoch and thus solar activity level and discuss their significance.

We have analyzed changes in the acoustic oscillation eigenfrequencies measured over the past 7 years by the GONG, MDI and LOWL instruments. The observations span the period from 1994 to 2001 that corresponds to half a solar cycle, from minimum to maximum solar activity. These data were inverted to look for a signature of the activity cycle on the solar stratification. A one-dimensional structure inversion was carried out to map the temporal variation of the radial distribution of the sound speed at the boundary between the radiative and convective zones. Such variation could indicate the presence of a toroidal magnetic field anchored in this region. We found no systematic variation with time of the stratification at the base of the convection zone. However we can set an upper limit to any fractional change of the sound speed at the level of $3 \times 10^{-5}$.

HD 89744 is an F7 V star with mass 1.4 M, effective temperature 6166 K, age 2.0 Gy and metallicity [Fe/H]= 0.18. The radial velocity of the star has been monitored with the AFOE spectrograph at the Whipple Observatory since 1996, and evidence has been found for a low mass companion. The data were complemented by additional data from the Hamilton spectrograph at Lick Observatory during the companion's periastron passage in fall 1999. As a result, we have determined the star's orbital wobble to have period P = 256 d, orbital amplitude K = 257 m/s, and eccentricity e = 0.7. From the stellar mass we infer that the companion has minimum mass m2 sin i = 7.2 MJup in an orbit with semi-major axis a2 = 0.88 AU. The eccentricity of the orbit, among the highest known for extra-solar planets, continues the trend that extra-solar planets with semi-major axes greater than about 0.15 AU tend to have much higher eccentricities than are found in our solar system. The high metallicity of the parent star reinforces the trend that parent stars of extra-solar planets tend to have high metallicity

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.

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.

Spectroscopic observations of 51 Pegasi and tau Bootis show no periodic changes in the shapes of their line profiles; these results for 51 Peg are in significant conflict with those reported by Gray and Hatzes (1997). Our detection limits are small enough to rule out nonradial pulsations as the cause of the variability in tau Boo, but not in 51 Peg. The absence of line shape changes is consistent with these stars' radial velocity variability arising from planetary-mass companions.

We report the discovery of near-sinusoidal radial velocity variations of the G0V star rhoCrB, with period 39.6 days and amplitude 67 m/s. These variations are consistent with the existence of an orbital companion in a circular orbit. Adopting a mass of 1.0 M(Sun) for the primary, the companion has minimum mass about 1.1 Jupiter masses, and orbital radius about 0.23 AU. Such an orbital radius is too large for tidal circularization of an initially eccentric orbit during the lifetime of the star, and hence we suggest that the low eccentricity is primordial, as would be expected for a planet formed in a dissipative circumstellar disk.