Interferometric observations of rapidly rotating stars

Optical interferometry provides us with a unique opportunity to improve our understanding of stellar structure and evolution. Through direct observation of rotationally distorted photospheres at sub-milliarcsecond scales, we are now able to characterize latitude dependencies of stellar radius, temperature structure, and even energy transport. These detailed new views of stars are leading to revised thinking in a broad array of associated topics, such as spectroscopy, stellar evolution, and exoplanet detection. As newly advanced techniques and instrumentation mature, this topic in astronomy is poised to greatly expand in depth and influence.


Introduction
One of the most fundamental stellar characteristics that is most frequently taken for granted is the shape of a star. This is perhaps somewhat unsurprising-in our direct experience of seeing stellar disks, there is only the sun, which is very nearly a perfect sphere: on an average radius of 959 .28 ± 0 .15 (Kuhn et al. 2004) there is only a variation of only 9.0 ± 1.8 mas (Rozelot et al. 2003) from equator to pole, indicating an oblateness (b/a − 1) of less than 10 −5 . However, for a surprisingly non-trivial number of stars, this degree of oblateness is in excess of 20% and, in certain cases, even 30%. Hints of this have been seen for decades in stellar spectroscopy, but this phenomenon has eluded direct observation until only recently.
The basis for departures from true sphericity lie in the rapid rotation of a star. Our sun turns about its axis at a stately rate of slightly more than once per month (Snodgrass and Ulrich 1990), which for a gaseous body of its size and mass leads to the ∼10 −5 oblateness value seen above. However, as we shall see below, certain objects with greater mass are spinning at rates that are 10-60× more rapid, leading to distortions of their shapes.
Direct observation of these oblate spheroids has been a tantalizingly close possibility for decades, but it has only been within the last 10 years that there has been a sufficient convergence of theory, capability and technique for imaging at submilliarcsecond levels of resolution to be possible. The early spectroscopic results a century ago guided theoretical expectations, which then unfortunately had to lay dormant while the technology caught up. The advent of modern long-baseline interferometry in the visible and near-infrared has opened the door to directly characterize and even image these rapidly rotating stars-and in particular, the most recent developments in this field of ultra-high-resolution astronomy have enabled rapid strides to be made in probing the surprising stellar structure of these objects.
In Sect. 2 we will explore the history of these developments, starting with the spectroscopic background (Sect. 2.1) which led to the early expectations and tests of interferometry (Sect. 2.2). The basic physics will be described in Sect. 3, including the simple Roche shape of a rigidly rotating gaseous spheroid (Sect. 3.1) and the resulting latitude dependence of flux known as the von Zeipel effect (Sect. 3.2). The specifics of interferometric observations of such targets is addressed in Sect. 4, including a general discussion of observational quantities (Sect. 4.1), and discussion of extensions of the technique to improved image reconstruction (Sect. 4.2) and greater spectral resolution (Sect. 4.3). A review of observational results to date is given in Sect. 5, followed by a discussion of the broad impact of these results to date Sect. 6; a summary of 191 future targets are suggested in Sect. 7.

Spectroscopic underpinnings
The earliest investigations into stellar rotation have their roots in Galileo's observations of sunspots (Drake 1957). 1 The technique of spectroscopy was developed in the late 19th century, and was the necessary tool to quantitatively observe rotational effects on stars, rather than just our sun. However, prior to that development, there is clearly evidence for consideration of the effect of rotation in the intervening years upon other astronomical observables, such as photometry (cf. the work of Bouillaud in 1667 on Mira, and later by Cassini, Fontenelle, and Miraldi as presented in Brunet 1931). Thus, roughly a dozen generations after Galileo, Abney (1877) was the first to suggest axial rotation of stars could be observed from spectral line broadening. This suggestion was swiftly rebuked by Vogel (1877), who pointed out that broadened hydrogen lines suggested by Abney are frequently accompanied by other lines which themselves appear narrow-it was only later established that there are many stars in which all lines are broadened. An actual measurement of the rotation effect on spectral lines was carried out first by Schlesinger (1909Schlesinger ( , 1911 on the eclipsing binaries λ Tauri and δ Librae. In Schlesinger's observations, the less luminous companion occulted varying parts of the rapidly rotating primary, allowing measurement of variations in apparent radial velocity. This 'Rossiter-McLaughlin' effect (Rossiter 1924;McLaughlin 1924) is now a phenomenon commonly observed with transiting extrasolar planets, and is a useful tool for probing the alignment of the orbital plane relative to the stellar rotation axis (see, e.g., Winn et al. 2005).
The specific line shapes expected from rotational Doppler broadening were first predicted by Shajn and Struve (1929). Although their work emphasized once again the observation of this effect in binaries, they included relevant predictions for single stars. These line contours were used by Elvey (1930) in publishing the first list of rotational velocities (although there had been a few efforts on specific individual stars a few years before), and the rest of the decade saw rapid observational progress. Struve and Elvey (1931) linked rotation rate to spectral type, finding that A-type stars were the most likely fast rotators; Westgate (1933aWestgate ( , 1933bWestgate ( , 1934 published extensive observational catalogs for hundreds of stars between spectral types O and F. These efforts solidified this new cornerstone of observational astronomy. Slettebak, in a dozen papers between 1949 and1956, discovered that the most rapid rotators were to be found among Be stars, and established a relationship between rotation and mass. Contemporaneously, von Zeipel (1924a, 1924b demonstrated that (under the assumption of rigid body rotation) the local surface brightness at any point on a star is proportional to the local effective gravity, and as such, the temperature at the poles would be greater than at the equator for a rotating star. Slettebak (1949) went on to take the implications of this 'von Zeipel effect' and computed the first modifications for expected spectral line shapes of rapidly rotating, bright stars. These implications were developed in detail in Collins (1963Collins ( , 1965 for continuum emission, and Collins and Harrington (1966), who incorporated shape distortion, aspect effects, gravity and limb darkening, and latitude variation in calculating H-β profiles. Harrington and Collins (1968) went on to characterize the intrinsic polarization expected for rapidly rotating early-type stars, and Collins and Sonneborn (1977) demonstrated that the spread of the main sequence (e.g. M V versus b − y coordinates) due to rotation alone was 2 to 3 times larger than previously expected. This latter result had significant implications for inference of both ages and distance moduli for clusters of young stars.
Overall, the apparent position of stars on the HR diagram is significantly affected (e.g. 2-3 subtypes) in considering rapidly rotating stars to their nonrotating counterparts (see Slettebak et al. 1980;Collins and Smith 1985, and references therein). For the interested reader, the detailed history of line profile analysis is given a substantively more thorough treatment in Slettebak (1985). On the broad subject of rotation the multiple editions of Tassoul (1978Tassoul ( , 2000 and proceedings of the two dedicated IAU meetings (Colloquium 4 in 1969 organized by Arne Slettebak, Symposium 215 in 2002 organized by Andre Maeder & Philippe Eenens, both titled simply "Stellar Rotation") are indispensable resources.

Interferometric observations
It is unclear as to when the possibility of directly measuring stellar rotational distortion with optical interferometry 2 was first considered as a plausible exercise. The pioneering angular diameter measurements of Michelson and Pease (1921) with the 20-foot beam interferometer on the Hooker 100 was the technology gate that opened up the possibility of such measurements. However, it is clear that the idea of potential observations did not fully develop until the guiding spectroscopic rotational velocities were themselves surveyed between 1930 and 1960, and the implications of the extremes of those velocities then evaluated. Furthermore, Pease's observational experience of operating the 50-foot beam interferometer, the successor to the 20-foot, was of extreme difficultly and highly limited success. Single measurements of individual stars were few in number, none were ever published in refereed journals (see the commentary in Townes 1999), and the collection of the sufficient data density necessary to establish a detection of rapid rotation on a given object was not forthcoming.
The entire endeavor of interferometric stellar angular diameter measurements lapsed into a state of dormancy for more than three decades, until the innovative proposal by Hanbury Brown and Twiss (1956) to pursue this task with intensity interferometry began to produce results. These efforts led to the construction and operation of the Narrabri Intensity Interferometer (NI 2 ), which produced the seminal results on single-star diameters (Hanbury Brown et al. 1974) and binary star orbits (Herbison-Evans et al. 1971). The success of NI 2 in this regard, combined with the maturity of the underlying theory of rapid rotators that predicted observable effects, led to considerations of the possibility of using NI 2 to observe the oblateness of Altair, as described in the PhD dissertations of Jordahl (1972) and Lake (1975). Unfortunately, the northern hemisphere location (δ = +8 • 52 ) of Altair is at odds with the southern location of Narrabri (latitude = −30 • 19 ), which limited NI 2 's ability to collect a sufficient range of baseline projections upon the object. Simple diameter measurements from NI 2 were published, but no detection of oblateness was made. (No mention of attempts to observe the more favorable declination, bright rapid rotator Acherner with NI 2 is obvious in the literature.) Another contributing factor in this missed opportunity may have been the novel audacity of the intensity interferometer-and the degree to which its underlying physical principles were (and still are) poorly understood by the majority of other astronomers in the field. Gaining general acceptance of the simple single-star diameter measurements was difficult as it was, and a more esoteric result like stellar oblateness was perhaps viewed as too ambitious. The completion of NI 2 operations brought this second era of optical interferometry to a close with no detection of this phenomenon.
The third, and current, generation of optical interferometers has met with resounding success in probing more than just 1-dimensional parameterizations of stellar sizes. Many of the facilities in use today are characterized by multiple apertures, some even being relocatable, delivering dense data sets that permit sophisticated image reconstructions. An early attempt was made by Robert R. Thompson and myself using the IOTA interferometer in 1998 to detect the rotationally induced oblateness of S Cep, spending two weeks of observing using multiple baselines to characterize the shape of the object. Unfortunately, data on our check star SS Cep indicated spurious results, and the data had to be discarded.
Stellar oblateness was first detected when a team I led observed Altair with the Palomar Testbed Interferometer (PTI) on two separate baselines , leading to disparate uniform-disk size measurements, particularly in comparison to the two contemporaneous size measurements on the check star Vega, which did agree. Extraction of a multi-parameter solution though use of Monte Carlo techniques, patterned after the similar approach applied to Keplerian orbit solving for interferometric data , indicated a v sin i = 210 ± 13 km s −1 , in agreement with spectroscopic values. The data were insufficient to further constrain the additional parameters of inclination or gravity darkening, but a new sub-field in optical interferometry was opened up by that study.
Subsequent observations of rapidly rotating stars have been carried out with the more capable facilities VLTI, NOI, and (especially) the CHARA Array. These ensuing data sets have been sufficiently rich to allow for detailed parameterizations of the observed objects, including constraints on inclination and gravity darkening. A full discussion on the seven objects studied in-depth to date will be reviewed in Sect. 5, and is summarized in Table 1.

Description of the geometry and basic physics
3.1 Roche model for stellar shape For a non-rotating star of mass M in radiative equilibrium, its uniformly spherical shape of radius R is trivially defined in terms of an equipotential surface (Eddington 1926): However, once rotation is imparted onto the star, a term must be added to account for the rotational potential:  'Oblateness' is defined as R eq /R pole − 1. In certain cases, equation (10) has been invoked to convert between ω/ω crit and v/v crit Fig. 1 Illustration of relevant rapid rotator geometry: observer is in XZ plane, with stellar polar axis of radius R p at inclination i to observer. Stellar angular rotation rate ω (and stellar mass M) sets equatorial radius R e . Stellar surface point at colatitude θ , longitude φ has radius R and local effective gravity vector g, with μ being the cosine of the angle between g and the observer. Adopted to be consistent with Fig. 2 in Collins (1965) where θ is the colatitude, R pole is the polar radius, Ω is the angular velocity (since our basic model here assumes uniform rotation, Ω = Ω(θ)), and R(θ) is the radius at a given colatitude; this 'Roche model' (Roche 1837) geometry is illustrated in Fig. 1. The last equality in (2) comes from maintaining that equality across the stellar surface at the non-rotating pole. Another important caveat of this simple treatment is that polar radius will be treated as a constant, regardless of rotation speed. A more accurate examination of this overall phenomenon accounts for a decrease in polar radius as the rotating star approaches maximum speedthis effect is slight for most stars ( 2%), although for 1M objects this effect is thought to increase rapidly with decreasing mass to a maximum of ∼8% (Sackmann 1970; Deupree 2011). A important implication of advancing from (1) to (2) is that stellar models need to graduate from one-dimensional treatments to twodimensional. The maximum "critical" angular velocity, Ω crit , where the outward centrifugal acceleration of rotation equals inward gravitational acceleration, can then be derived as At this rotational speed the oblateness of the object is at its greatest, 3 with R eq,critical = 3/2R pole . The fractional angular velocity ω is then defined as and solving for the cubic equipotential Φ(R pole ) = Φ(R, θ) with trigonometric methods (see Collins 1963;Collins and Harrington 1966;Jordahl 1972), 4 we arrive at an expression for the stellar radius at any colatitude θ for a given fractional rotational velocity ω: Owocki et al. (1994) identify the convenient identities where w 0 ≡ v 2 eq R pole /2GM and v eq is the equatorial rotational velocity, presumably derived from spectroscopy; these identities allow use of the solution for R(ω, θ) in (5). Figure 2 illustrates the expected distortion imparted upon a (Rocheapproximated) stellar surface as it progresses from a non-rotating to a rapidly rotating situation. For this object, an increase of rotational speed from ω = 0 to ω = 0.92 induces a ∼15% increase in the equatorial radius. Owocki et al. (1994) point out that using the point-mass form for the gravitational potential found in (2) ignores high order multipole components that might arise from the rotationally distorted stellar mass distribution. However, Orlov (1961) shows that polytropic stellar structure models using the correct potential yield a change of less than 1% from the oblate surface radii predicted by (5).
To compute the surface gravity at a given colatitude, we compute the negative gradient of the effective potential in (2). The two components of the local effective gravity g at a given stellar colatitude θ in spherical polar coordinates are Thus g, directed inwards along the local surface normal, is equal in magnitude to g 2 r + g 2 θ (Cranmer and Owocki 1995; Aufdenberg et al. 2006). The further details of the mathematics describing rotating Roche equipotential surfaces and it application to the stars can be found in Collins (1963), Collins and  Harrington (1966), Jordahl (1972), Tassoul (1978), Cranmer andOwocki (1995), van Belle et al. (2001), and Aufdenberg et al. (2006). Interesting related discussions regarding the oblateness of a planetary objects can also be found in the appropriate literature (e.g. Baron et al. 1989 wherein a discussion of Uranus' oblateness is considered).
For the benefit of the reader who often encounters rotational velocities from spectroscopy expressed in terms of v sin i (rather than angular velocity expressions), it is useful to note here a few relationships. Starting with (3), the critical linear velocity at the equator can be derived as Also, there is a subtle difference between fractional angular velocity and fractional linear velocity, since the linear radius that relates angular velocities to linear velocities is different between actual speed and the breakup speed: 3.2 The von Zeipel effect von Zeipel (1924a) and later Chandrasekhar (1933) showed that a rotationally oblate star has a local surface flux that is proportional to the local gravity, F ∝ g. Thus, the Stefan-Boltzmann law, F = σ SB T 4 eff implies that T eff ∝ g 0.25 . This particular proportionality is applied in the case of radiative equilibrium; Lucy (1967) pointed out that for stars with convective envelopes, a weaker dependency is expected, with T eff ∝ g 0.08 given as 'representative' (although the value of 0.08 is weakly dependent upon mass, radius, luminosity, and other parameters). A closer examination of the Lucy (1967) discussion points out that progressing from slow to rapid stellar rotation should lead to pressure gradients that need to be balanced by the Coriolis force, with the consequence that the value of 0.08 quickly becoming less valid. A general parametrization expresses the relationship as T eff ∝ g β , with β either being a model predicted value or a parameter to fit as part of a larger overall solution for some set of observed astronomical data.
Thus, when considering the total flux from a rapidly rotating star, it is important to consider integrating over the surface accounting for the flux as it is tied to the local effective temperature. This colatitude-dependent temperature can be written as allowing for a modification of the Stefan-Boltzmann law for characterization of the local stellar radiative flux, The von Zeipel parameter β has been tabulated for a range of stellar models varying in mass and age by Claret (1998Claret ( , 2000. 5 It is worth noting that the validity of (11) and (12) has been questioned by Connon Smith and Worley (1974) when the effects of differential rotation are considered. For the case where the concerns regarding differential rotation may be neglected, Aufdenberg et al. (2006) use the preceding formalism to derive the temperature difference between the equator and pole: It should be noted that use of a single 'β' term in describing a star perhaps oversimplifies the true nature of the stellar surface. As noted above, surfaces that are fully convective or fully radiative are expected to have values that are 0.08 and 0.25, respectively. However, as the temperature ranges from hot to cool between a stellar pole and equator, the dominant mode of heat transportation may shift from radiative to convective, depending upon the value of those temperatures. One may easily conceive of an object whose temperature profile ranges pole to equator, straddling a 'sweet spot' value where such a shift occurs-the low latitudes would be convective, the high latitudes radiative. As such, a colatitude-dependent value of β(θ) may be appropriate to consider in constructing the most valid stellar model; the surface appearance of such an object will be most intriguing to image.
In light of the recent observational results that will be reviewed in Sect. 5, Espinosa Lara and Rieutord (2011) re-examine the underlying assumptions of the von Zeipel law that lead to the T EFF ∝ g 1/4 EFF formulation. In particular, the 'strong hypothesis' of barotropicity (pressure only depends on the density) is already poor for slowly rotating stars (first noticed by Eddington 1925, see review in Rieutord 2006, and breaks down in noticeable ways for rapid rotators. These authors propose a new gravity darkening model, based on energy flux being a divergence-free vector antiparallel to the effective gravity, which appears to have good agreement with the interferometric results on Altair and Regulus.

Observational quantities, predictions
Optical interferometers are, by nature, telescopes that can rightly be considered rather 'non-traditional'. Elements of an optical interferometer would be eminently familiar to the non-specialist (e.g. the individual telescopes, for example), whereas there are also elements that would be rather alien (e.g. the beam-combining back end). This is also true of the elements of the data products-familiar concepts, such as photometry, quickly give way to unfamiliar ones, such as fringe visibility. For the purposes of this review, we will briefly introduce the relevant observational quantities in a concise manner, directing the reader to more thorough discussions should she or he be interested.
To begin on familiar territory, consider a 'normal' filled-aperture telescope. If it were to observe an unresolved point source, one would see on the telescope's image plane an Airy pattern with the characteristic width of 1.22λ/D radians, where λ is the wavelength of observation, and D is the telescope aperture size. For a more complex object with on-sky brightness distribution O(α, δ) being observed by a telescope with point-spread function P(α, δ), the mathematical description of this process is For the simple case of a point source, O(α , δ ) reduces to a Dirac delta function, P(α − α , δ − δ ) is the point source function associated with a filled-aperture telescope, and I(α, δ) returns the Airy pattern. The Fourier transform of (15) gives which takes us directly to the points of interest for this review. Specifically, as a curiosity of their construction and measurement processes, optical interferometers tend to provide measurements that sample the Fourier information associated with an object, I (u, v), rather than image plane information (I(α, δ)). Secondly, the transfer function associated with an optical interferometer, T (u, v), can be well-characterized. Thus, since the information associated with I (u, v) and T (u, v) can be collected, the nature of O (u, v) can be computed-and by extension, the original object brightness distribution O(α, δ). For further expansive details on the underlying theory of optical interferometers, reference the van Cittert-Zernike theorem in a general optics textbook such as Born and Wolf (1980), or books specific to this topic (e.g. Goodman 2005). The proceedings of recent summer schools on optical interferometry are also very instructive in considering these instruments further (in particular, Lawson 2000;Haniff 2007aHaniff , 2007b, as is Bracewell (2000)'s book on the Fourier transform itself. In practice, an optical interferometer takes light from two or more telescopes and recombines in such a way that the light from the apertures combines coherently, or interferes with itself (hence the name). Typically the interference criteria include maintaining wavefront quality, polarization, and path length from the point of light collection at the individual telescopes, through the system, to the point of beam recombination. In the simplest case of two telescopes, one may consider the interference of light as it is recombined and falls on a photodetector. Moving a mirror in the beam train that brought the light to the recombination point, the path length through the system may be swept through the point of equal path between the two arms of the interferometer, producing constructive and destructive interference. This characteristic pattern is referred to as a 'fringe', with its contrast (the intensity as measured by the photodetector between maximum constructive and maximum destructive interference) being the measured visibility: As related to our discussion of (16), the visibility V is the amplitude of a single point of complex image information I (u, v). Fringe visibility is the basic observable of an optical interferometer and is directly related to the angular size of objects being observed by the instrument. A second observable, closure phase, is also frequently measured by optical interferometers. Transmission of starlight through the turbulent atmosphere corrupts the phase information at optical wavelengths, so an absolute measurement of the phase information associated with I (u, v) is not possible. However, with interferometers consisting of three or more telescopes, these corrupted phases can be summed about triangles to recover partial, uncorrupted phase information. If we consider a triangle of points from three telescopes {l, m, n} that form baselines sampling the complex image I (u, v), the measure phased for pair (m, n) is ψ mn , which consists of the source phase φ mn and phase errors (ξ m − ξ n ). Thus, summing around a triangle of {u, v} points, we see that the measured sum CP is directly related to the original uncorrupted φ values, since the ξ errors cancel. Closure phase is an observable that is sensitive to the degree of asymmetry in an image I (u, v). A good example of a striking closure phase signature is in the signal found in resolved stellar disk observations during planet transits, which can be considered 'perfect' star spots (van Belle 2008). Additionally, higher order such constructions, such as the closure amplitude (CV lmno = V lm V no /V lo V mn ) are also possible but will not be considered herein.
To take this discussion of interferometric observables and place it in a context of actual measurements, we shall consider the representative, pioneering case of observing the rapid rotator Altair, which we shall see in Sect. 5 is a perennial favorite for demonstrating the prowess of an interferometer or instrument. In Fig. 3 we can see the basic appearance of Altair in the {u, v} plane; the left panel is a representation of the visibility amplitude, and the right panel illustrates the visibility phase. These points were generated by construction of a Roche model based on the rotation parameters reported in Monnier et al. (2007) (R pole = 1.634R and π = 194.45 mas, ω = 0.923, i = 57.2 • , α = −61.8 • , β = 0.19, T pole = 8450 K), and using the prescription found in Sect. 4 of Aufdenberg et al. (2006) to generate colatitude dependent values for radius R(θ), local surface gravity g(θ), and local effective temperature T EFF (θ ). (Limb darkening was ignored in this example, which has a small ∼1% effect on visibility amplitude measurements.) T EFF was then converted to 1.6 μm flux under the simplifying assumption of blackbody radiation (BBR). For a qualitative comparison of the effect rapid rotation has on these observables, a second nearly identical set of {u, v} points was generated, with the exception that the rotation rate was artificially reduced from ω = 0.923 to ω = 0.01, in Fig. 4. The general deviation of these two figures are most clearly seen in the right-hand phase plots: the uniform source of Fig. 4 has a clear plane associated with a φ = 0 value, while the asymmetric appearance of the gravity-darkened rapid rotator in Fig. 3 shows a continuum of interesting phase values.
These differences, and the associated visibility amplitude excursions associated with the stellar shape produced by rapid rotation, are sharply highlighted in Fig. 5, which plots the difference of Figs. 3 and 4. Not only are the obvious phase excursions seen on the right, but previously hidden ±10% visibility amplitude differences are made obvious. The deviations in visibility amplitude on the left a due primarily to the increased angular size on the sky in the stellar equatorial region; the closure phase  . We see that the observed visibility deviates for a rapid rotator from the values expected for a non-rotating object at the ±10% level; the phase values vary by many tens, if not hundreds, of degrees-and as such, the resultant closure phase values would deviate on similar scales. Instruments such as CHARA Array MIRC and VLTI AMBER are able to measure visibility amplitude and closure phase; for the former, visibility amplitude precision is at the ∼σ V = 5% level, and closure phase precision is as good as ∼σ CP = 1 • for single-measurement accuracy with substantial √ N improvements possible (Zhao et al. 2011). Figures 7, 8, 9, 10, and 11 illustrate differences in visibil-   ity amplitude and phase (e.g. Fig. 3) as the solution parameters {ω, R pole , i, α, β} found in Monnier et al. (2007) are given 1-σ excursions from their best-fit values.  The practicalities of χ 2 minimization of interferometric data are highlighted in Sect. 6 of van Belle et al. (2006) (henceforth referred to as VB06), which is illustrated pictorially in Fig. 12. VB06 used a Monte Carlo approach which began by constructing models of a rapidly rotating star based upon rotation ω and polar radius R pole . Model surfaces were constructed for the star at some sufficiently dense interval (for VB06, the spacing was 0.8 • ) in both colatitude and longitude across the whole volume. Flux for a given surface area can then computed using the appropriate influence of gravity darkening as described in Sect. 3.2, with the β parameter quantifying the effect, anchored to pole effective temperature T pole . In addition to the gravity darkening, limb darkening appropriate for these model stars as indicated by the quadratic laws can be generated by Claret and Hauschildt (2003). These models are then mapped onto the sky, through the use of two additional free parameters describing inclination i and on-sky rotational orientation α.
Thus, for a given set of six randomized free parameters {ω, R pole , i, α, β, T pole }, a ∼100,000 element volume surface is generated, projected upon the sky, rotated and the resultant image Fourier transformed for comparison to each of the observed visi-  (2006) bility and closure phase data points, from which a χ 2 /DOF is calculated. As cleverly noted in Aufdenberg et al. (2006), the full Fourier transform of the model star's image does not need to be performed-only visibilities and closure phases associated with a small, specific set of {u, v} points corresponding to the interferometer sampling need to be generated. The resulting computational load is substantially lighter and allows for faster and more complete Monte Carlo searches. A detailed review of optical interferometry in astronomy can be found in Monnier (2003). A detailed discussion of phases in particular is continued in Monnier (2007), including expansive discussions on closure phase, and more esoteric observables possible with optical interferometers, such as differential closure phase, closure differential phase, and closure amplitudes.

Image reconstruction
Relative to the VB06 study, there are two specific extensions that can be made to this technique. First, inclusion of closure phase data can significantly expand the capability of the Monte Carlo routine to quickly recover the six characterizing parameters of the original star. The sensitivity of closure phase to asymmetries in brightness distribution speeds and improves the fidelity of the parametric fitting approach. The code written for VB06 in fact already includes generation of the closure phase data, but no comparison was made in that investigation since the CHARA Array's Classic pairwise beam combiner (from which the VB06 data were obtained) produces no measures of closure phase. However, both the CHARA Array's MIRC combiner and the VLTI AMBER instruments provide measurements of closure phase, and data from these instruments could be utilized for the parameter fitting; an example of this can be found in Sect. 4 of Zhao et al. (2009).
Second, and most instructively, optical interferometry instrumentation and techniques is advancing beyond simple pairwise telescope recombination to now include 3-, 4-, and even 6-way telescope recombination, which allows model-independent maps to be recovered from the data. These 'imaging' techniques are advancing rapidly (Cotton et al. 2008;Baron et al. 2010), and their great strength lies in the fact that they use little or no a priori assumptions, making them powerful, independent checks on the model-dependent approach to interpretation of interferometric data.
An excellent example of image reconstruction can be found in Sect. 3 of Zhao et al. (2009), who use the Maximum Entropy Method of radio synthesis imaging (Narayan and Nityananda 1986) as developed by Ireland et al. (2006) into the "Markov-Chain Imager for Optical Interferometry" (MACIM) code. Blind tests of the results from MACIM (Lawson et al. 2006) have shown it to be a robust tool for image reconstruction with optical interferometric data. For the application to rapid rotators, Zhao et al. (2009) incorporated an ellipse prior, justified by the reasonable expectation that the stellar photosphere has a sharp cutoff at the edge. This prior was generated through MACIM imaging of a grid of uniform surface brightness ellipses. Additionally, a Gaussian beam deweighting (as is common in radio synthesis imaging) was also applied, for image smoothness. Following these considerations, the final image was the one with the global maximum entropy. As closure phase precision progresses to levels much better than a fraction of a degree , the fidelity of these image reconstructions will progress as well.

Spectroscopy and spectro-interferometry
Another ongoing, extremely powerful extension of these interferometric techniques is the combination of them with unresolved spectroscopy, and extension of the interferometric observables through spectral dispersion as well. For example, the rapidly rotating model generated to fit the interferometrically observed stellar surface brightness distribution of Vega ) was upgraded to generate high dispersion spectra (Aufdenberg et al. 2007) and found that the high rotational velocitynecessary to reproduce the surface brightness distribution-was incompatible with the lower rotational velocity preferred by the line fits. These observations were taken to suggest that a simple von Zeipel gravity darkening law and uniform rotation are insufficient to describe the stellar photosphere. Yoon et al. (2010) demonstrate a further extension of the approach that attempts to unify high-resolution spectroscopy with interferometry, deducing that Vega's metallicity and mass are both lower than previously expected.
Additionally, the interferometric beam combiner may itself spectrally disperse the recombined light, permitting a wavelength-dependent analysis of the stellar surface brightness distribution. Sufficiently high-resolution dispersion, such as the R ∼ 12,000 available with the VLTI AMBER instrument (Petrov et al. 2007), make it possible to examine spectral line profiles in a spatially resolved or even imaged sense. AMBER has already been used for spectro-interferometric observations of the Herbig Be star MWC 297 as this high spectral resolution (Weigelt et al. 2011).

Stars studied to date
An increasing number of facilities have been successfully employed to make a measurement of rotationally induced oblateness, starting in 2001 with the Palomar Testbed Interferometer (PTI, Colavita et al. 1999), the Very Large Telescope Interferometer (VLTI, Glindemann et al. 2003), the Navy Prototype Optical Interferometer (NOI, Armstrong et al. 1998), and by 2005 the Center for High Angular Resolution Astronomy (CHARA) Array (ten ) was demonstrating the technique as well. These interferometric observations were interpreted with parametric fits similar to the example described in Sect. 4.1; recently, interferometric imaging (Sect. 4.2) has been demonstrated as a powerful tool in investigating these objects in a significantly more model-independent approach. Below we will review the studies to date; a timeline of these studies and their achievements is presented in Table 2.

Altair
Of the rapidly rotating stars, Altair (α Aql, HR7557, HD187642; α = 19 h 50 m 47.00 s , δ = +08 • 52 m 05.96 s ) is one of the nearest, brightest prototypes, and located close enough to the celestial equator to be observed by both northern-and southernhemisphere facilities (d = 5.143 ± 0.024 pc, m V = 0.77, m K = 0.10, α = 19 h 50 , δ = +08 • 52 ; Perryman et al. 1997;Cutri et al. 2003). Altair is an A7IV-V star (Johnson and Morgan 1953), the 12th brightest of all stars, and has been known to be a rapid rotator for over half a century (240 km s −1 from Slettebak 1955). These characteristics made it an attractive target for the first interferometric observations of rapid rotation, and continue to make it appealing for instrument commissioning.
The interferometric observational implications were first considered theoretically in the PhD theses of Jordahl (1972) and Lake (1975), who investigated the effects that rapid rotation would have upon the shape, brightness distribution, and spectral features of that particular star. The observations of the object's size by the Narrabri Intensity Interferometer in 1965 were published (Hanbury Brown et al. 1974), but with only a uniform disk fit to the data of 2.78 ± 0.13 mas.
More than 25 years would pass between the publication of Hanbury Brown et al. (1974) and the first direct observation of rotational oblateness by van , with Altair as the target. During that time, the possibility of making such a measurement was consciously known by investigators in the field, and attempts were even made at such a measurement on various targets, without success. However, using the Palomar Testbed Interferometer (PTI), we were able to take advantage of a number of aspects of the instrument to establish departures from spherical symmetry. First, PTI provided at the time of the observations two baselines of sufficient length (>80 m) to resolve the ∼3.2 mas disk of Altair in the K-band. The baselines also provided sufficient {u, v} coverage on the sky, such that the angle between those baselines when projected onto the sky was on order ∼50 • . Second, the visibility data quality was such that the night-to-night repeatability was (under appropriate weather conditions) on order σ V 2 0.018, which was necessary to compare the differing baseline's data with sufficient precision to indicate departures from spherical symmetry. PTI's highly automated operations also allowed for rapid collection of a large body of data points with each baseline as well. Finally, the nature of the experimental design showed the apparent change in interferometric visibility with sky projection angle was not seen with data collected contemporaneously on the check star, Vega, as expected from that star's small v sin i value (21.9 ± 0.1 km s −1 ; Hill et al. 2004).
Collection of the observational data was only the first necessary step in achieving the confirmed detection of rotational oblateness. The proper direction for interpretation of the data was, at first, unclear, given the novel nature of the investigation at the time. Many of the publications of the time that discussed stellar rotation concentrated upon spectral line profile shapes but not the apparent shapes of the stellar disks. Indeed, some of the best early guidance was found in solar system literature discussing the shapes of the rapidly rotating gas giants (e.g., Baron et al. 1989). However, as recommended to us by George Collins (private communication), we found guidance in using the dissertation of Jordahl (1972) set against a computational methodology similar to that employed for reducing interferometric binary star data (e.g., Boden et al. 1999). This framework made it possible to extend the result from a mere toy model for a star with an elliptical on-sky appearance, to a physically meaningful family of dynamic Roche models of various inclinations with a single characteristic value of v sin i = 210 ± 13 km s −1 , which was then compared, successfully, against the spectroscopic values.
This 2001 investigation contains one notable error, 6 which fortunately does not alter the physics or the projected rotational velocity result: the {u, v} coordinates were inadvertently swapped, resulting in an incorrect on-sky orientation angle (for example, Fig. 6 in that manuscript should be rotated about the α = δ diagonal). The reported value in the manuscript of −21.6 ± 6.2 deg corrects to −68.4 ± 6.2 deg when this is taken into account.
This result on Altair quickly set the stage for further investigations of the phenomenon, with a variety of instruments and technique coming to bear on this particular object, along with other additional rapid rotators, beginning with Achernar the following year (as will be discussed in the next section, Sect. 5.2). This result was noted in the general astrophysical review literature as being a significant development (Trimble and Aschwanden 2003). With Altair, Ohishi et al. (2004) followed the initial result with observations of the star with the Navy Optical Interferometer (NOI). These investigators not only confirmed the oblate nature of the stellar photosphere, but reported a detection of a non-zero closure phase, consistent with a bright polar region produced by gravity darkening. This result was expanded upon with the addition of observations from the Very Large Telescope Interferometer (VLTI) by Domiciano de Souza et al. (2005), who synthesized together the PTI, NOI, and VLTI data sets into a unified solution. The following NOI observations and analysis by Peterson et al. (2006a) was able to unambiguously establish a stellar inclination for Altair.
Ultimately, these studies collectively set the stage for the ground-breaking study by Monnier et al. (2007), which was the first to directly image a main sequence star. The data from the MIRC instrument on the CHARA Array directly confirmed the presence of gravity darkening using an image reconstruction independent of the previous parametric modeling approaches, while at the same time providing evidence that departures from the standard formulation of that effect were insufficient to explain the observations. Using a value of β = 0.25, the investigators were able to fit their data in a reasonably satisfactory way, but their goodness-of-fit criteria significantly improved when this was allowed to be a free parameter, which converged on a value of β = 0.190 ± 0.012. Differential rotation, opacity, and convection-related phenomena are all cited as possible reasons for this value of β; detailed spectrointerferometric line profile analysis (as noted in Sect. 4.3) is suggested as the preferred tool to disentangle which of these effects is the dominant one.
Given its size and brightness, one obvious further application of interferometric observations of Altair would be an attempt to directly observe stellar rotation with time-resolved imaging. The ∼10 hour rotational period of star would have to be accounted for in such observations-no one exposure could be overly long, lest it risk smearing the image; 'snapshot' imaging would be a must. However, this could in fact be an advantage: during a single evening of observing, the full face of the star would be seen over time, allowing surface morphology maps to be built on a night-to-night basis. Such imaging would require a non-homogeneous surface, but that may be possible (even likely) given the known magnetic activity of Altair (Robrade and Schmitt 2009).

Achernar (α Eri)
Following on the heels of the Altair observations, Domiciano de Souza et al. (2002) presented theoretical details on modeling rapid rotators which were soon applied to VLTI observations of the rapidly rotating Be star, Achernar (α Eri, HR472, HD10144; α = 01 h 37 m 42.85 s , δ = −57 • 14 m 12.33 s ; Domiciano de Souza et al. 2003). This object is the southern-hemisphere complement to the northern-hemisphere object Altair: it is bright and presents an remarkably large variation of angular size as a function of sky orientation. In fact, in evaluating their VLTI data as a solitary, rapidly rotating photosphere, the investigators found an extreme oblateness ratio of a/b = 1.56 ± 0.05-a degree of photospheric distortion that could not be explained by Roche approximation. The authors claimed that the lack of Hγ line emission during their observations was indicative of no disk contamination of the interferometric visibility data, and that "the classical assumption of the Roche approximation becomes questionable". Investigations of more exotic photospheric models were spurred onwards by this claim (e.g., Jackson et al. 2004).
However, ensuing investigations of this object (Kervella and Domiciano de Souza 2006;Carciofi et al. 2008) branched out to add contributions from a circumstellar environment (CSE) to model, which lessen the derived oblateness ratios; Kervella and Domiciano de Souza (2006) thereby derive the photospheric oblateness ratio dropped to a/b = 1.41 ± 0.053 with a photosphere-to-CSE relative flux ratio of 4.7 ± 0.3% Additionally, there has been detection of a close-in companion by Kervella et al. (2008); however, inclusion of the companion and a rotationally distorted star alone in the visibility fits cannot explain the interferometry data, and a CSE contribution is still required (Kanaan et al. 2008). Companion thermal-IR interferometric observations with VLTI-MIDI has provided data consistent with extended emission associated with a fast wind ejected along the hot polar caps of star (Kervella et al. 2009).
Clearly this particular object has found a vigorous existence in the realm of interferometric study. However, in contrast to Altair, the complications of disks and winds have presented additional challenges (and rewards) for investigators attempting to fully fit models for Achernar. As the brightest Be star prototype, however, this object is uniquely situated to allow in-depth interferometric investigations to connect the effects of stellar rotation to the overall Be star phenomenon.

Regulus (α Leo)
As with Altair, Regulus (α Leo, HR 3982, HD 87901, ADS 7654A; α = 10 h 08 m 22.31 s , δ = +11 • 54 m 01.95 s ) was observed by the Narrabri Intensity Interferometer (Hanbury Brown et al. 1974), with similar result: insufficient {u, v} coverage prevented an in-depth characterization of this rapidly rotating star, although a gross size characterization was reported. Regulus had been recognized by Slettebak (1954)   . In addition to demonstrating the capabilities of the new array, this investigation was notable in the specific context of rapid rotator study for the following reasons: First, it marked the first time a combined interferometric and spectroscopic approach was applied to the data reduction methodology. Second, it was the first such investigation that resulted in a reasonably specific constraint on the inclination of the rotation axis for such an object. As of the publication of the Regulus results, previous investigations had painted broad suggestions on what various objects' true inclinations would be, but no specific values had be offered. The Regulus error was still quite generous (i = 90 +0 −15 deg) and yet specific, representing one of the ultimate benefits of this manner of investigation.
The resulting solution set for Regulus is a rather rich characterization of the object, including values for inclination (i), on-sky orientation (α), gravity darkening (β), rotational velocity (v), fractional rotational velocity (v e /v c ), polar and equatorial radii (R p , R e ), polar and equatorial temperatures (T p , T e ), mass (M), luminosity (L), and interstellar extinction A V . Such a rich data set allows for in-depth evaluation of the object on a variety of fronts. For example, comparison of the on-sky orientation with its known proper motion shows that Regulus's space motion is roughly along the axis of its rotational spin. This observation is thought-provoking from the context of considering the formation history of the star and the dynamics of its protostellar cloud.
More specifically, providing a the true inclination solution had far-reaching consequences. From the inclination, the true rotation velocity and fractional rotational velocity could be established. The degree to which such objects are rotating close to their breakup velocity has significant implications regarding the nature of mass loss in Be stars (e.g., Porter and Rivinius 2003). The range in temperature from the pole to the equator has cast doubt on the 'young' age estimate of 150 Myr (Gerbaldi et al. 2001) and a significantly older estimate of 1 Gyr has been suggested (Rappaport et al. 2009). The recent discovery of close companion by Gies et al. (2008) to Regulus 7 has led to consideration of the object's 'huge' quadrupole mass moment Q (Iorio 2008); estimates of Q indicate that detection of the correction to the Keplerian period due to Q would be possible with period measurements precise to ∼20 seconds; however, the current state of the art for the period measurement is about two orders of magnitude beyond that at P = 40.11 ± 0.02 days (Gies et al. 2008).
Also as with Altair, Regulus has now been imaged by CHARA-MIRC ( (Table 1), but further rich detail has been added from the imaging. In particular, a much tighter constraint on inclination is now available (i = 7 In addition to its close companion, Regulus has a wide companion α Leo B at a separation of ∼175 , also a binary. However, at such a distance, it will not have ever interacted with Regulus (McAlister et al. 2005). 86.3 +1.0 −1.6 ), along with tighter constraints on pole and equator temperature, and further strong evidence of departures from the 'standard' von Zeipel gravity darkening, with β = 0.188 +0.012 −0.029 .

Alderamin (α Cep)
A second CHARA Array commissioning target, Alderamin (α Cep, HR 8162, HD 203280; α = 21 h 18 m 34.77 s , δ = +62 • 35 m 08.06 s ), was also observed with the array, also benefiting from the large amount of {u, v} coverage afforded by the facility's multiple baselines (van Belle et al. 2006). Alderamin is an A7IV star (Johnson and Morgan 1953), making it an attractive observing candidate given its similarities to Altair. Alderamin is 3× further than Altair but, strangely, at roughly the same X-ray luminosity (Robrade and Schmitt 2009); this suggests that coronal X-ray emission sets in at about A7 for main sequence stars (at least at the sensitivity levels of current instrumentation). Alderamin has been known to be a rapid rotator for many decades (Slettebak 1955), with spectroscopic estimates rotational velocity that range from v sin i = 180-200 km s −1 (Gray 1980;Abt and Morrell 1995) to 245-265 km s −1 (Abt and Moyd 1973;Bernacca and Perinotto 1970). This object was recently reobserved (Zhao et al. 2009) with the array using the newer multi-way MIRC combiner, which not only provided visibility amplitude data on baseline pairs but closure phase data from baseline triplets. Both the van Belle et al. (2006) and Zhao et al. (2009) investigations attempted to fully characterize the fundamental parameters associated with the rapidly rotating stellar photosphere, including inclination i, sky position angle α, polar and equatorial radii and temperatures, and most significantly, gravity darkening β. As noted in Sect. 3.2, a single value of β is perhaps an incomplete characterization of the true physical nature of the star's surface, particularly given the range of temperatures, from T pole = 8600 ± 300 K to T equator = 6600 ± 200 K; following the canonical temperatures of static models found in Cox (2000), we see this ranges from spectral types ∼A3V to ∼F5V, respectively. This range of temperatures brackets the radiativeto-convective transition temperature of ∼8300 K (Christensen-Dalsgaard 2000; Neff and Simon 2008), allowing for equatorial convective activity but no similar polar activity (Rachford and Foight 2009).
Contrasting the results from the earlier van Belle et al. (2006) study with that of Zhao et al. (2009) is instructive in illustrating the power of various interferometric data sets. The former, having only visibility amplitude data available (only the 2-way CHARA Classic combiner was online at the time of the study), came up with an inclination much larger than the latter (88.2 +1.8 −13.3 versus 55.70 ± 6.23 degrees). Simply put, the former was rather insensitive to the asymmetric stellar surface brightness distribution, which manifests itself in the poor value for i and (probably unreasonably) low value for β; the latter, using CHARA-MIRC data which produce data sets not just with size-sensitive visibility amplitude data, but also with asymmetry-sensitive closure phase data, was able to produce a significantly more precise-and presumably accurate-values for i and β. This in turn affected the entire ensemble of the values in the solution set; it is interesting that the best value for gravity darkening still deviates from the 'von Zeipel ideal' with β = 0.216 ± 0.021.
It is also noteworthy that the similar visibility-only solution for Regulus in McAlister et al. (2005) also found a lower value for β than the spectroscopy/visibility solution found in the same paper (β = 0.13 ± 0.05 versus 0.25 ± 0.11). Investigators may wish to consider if such data sets are insufficiently rich to providing constraints upon the stellar surface morphology to characterize gravity darkening: an excellent example of this can be seen in Fig. 12 of Zhao et al. (2009), where for certain onsky orientations, degeneracies may form between various parameters (in this case, β and i, when i is large). In particular, for those investigators wishing to establish interferometry-only data sets to determine absolute rotation velocities for comparison with spectroscopic v sin i values, closure phase characterizations of asymmetric surface bright distributions will be highly sensitive to rotation axis inclination values.
Extending the implications of the observations, in van Belle et al. (2006), a simple analysis of the rotation history of the star via the simple application of conservation of angular momentum was considered, and implied that its present fractional breakup velocity (v e /v c = 0.83 ± 0.05-similar to the Zhao et al. (2009) value, so this line of reasoning is unaffected by the poor values of i and β in that study) would have been even higher in its zero-age main sequence history, with v e /v c 0.92-0.98. Evidence of such an large earlier rotation rate has significant implications for the star's formation history.
In their study of the onset of convective zones in main sequence stars, both Altair and Alderamin are cited by Neff and Simon (2008) as cases that demonstrate the lack of fidelity in categorizing stars solely by spectral type or mass. These same investigators found no FUSE UV evidence for deepening of convective zones due to the rapid rotation of these objects, as suggested not only by the interferometry data but also by the models of MacGregor et al. (2007).

Vega (α Lyr)
Vega (α Lyr, HR7001, HD172167; α = 18 h 36 m 56.34 s , δ = +38 • 47 m 01.29 s ) is a particularly interesting special case of an observed rapid rotator in that it is viewed nearly pole-on relative to the other examples cited thus far. The long-standing use of Vega as a spectrophotometric standard (e.g., Hayes 1967;Straizys et al. 1976) makes understanding the surface brightness distribution of this object at the highest levels of detail more than an abstract scientific exercise, and one of fundamental utility. The discovery of Vega's debris disk (Aumann et al. 1984) due to its mid-infrared excess flux called into question its validity as a fundamental standard. This situation has been further complicated by the fact that the star is a rapid rotator, as first suggested by Gray (1988). This result, based upon detailed spectroscopic line modeling, has been confirmed by interferometry, as we shall see below.
As with Altair and Regulus, Hanbury Brown et al. (1974) observed Vega and obtained a limb-darkened angular size of 3.24 ± 0.07 mas (and a uniform-disk size of 3.08 ± 0.07 mas) through a 443 ± 5 nm filter. As with all of the objects in this review, the intervening generation of optical interferometer prototypes after NII-GI2T, IRMA, IOTA (to name a few)-did not observe Vega as a scientific target, due to limited spatial resolution. Further scientifically significant data from interferometer on this object came with the longer baselines of PTI when Ciardi et al. (2001) observed it and found puzzling signs of residuals in the K-band fits, consistent with a debris disk signal contaminating the stellar photospheric signal (as discussed in more detail in Barnes 2009;Lawler et al. 2009;Akeson et al. 2009). This finding was consistent with the further investigations with the CHARA Array ) and IOTA (Defrère et al. 2011). These PTI and CHARA Array interferometric studies set the stage for two further, more detailed, studies of the star's photosphere itself.
Both NOI (Peterson et al. 2006b) and the CHARA Array  published studies of Vega early in 2006, in the visible (∼500-800 nm) and K -band, respectively. Both studies found strong evidence for rapid rotation in a nearly pole-on star, arriving at remarkable agreement in inclination (4.54 ± 0.33 • and 4.7 ± 0.3), equatorial rotational velocity (274 ± 14 and 270 ± 15), and other parameters (see Table 3) despite completely independent development of their methodologies and differing observational wavelengths. Agreement between these two disparate approaches is reassuring; however, further investigations by the NOI group (Yoon et al. 2008(Yoon et al. , 2010 indicate a lower rotation speed (v eq = 175 ± 33) coupled with lower than previously expected mass and sub-solar metallicity (2.14 ± 0.08M , Z = 0.0080 ± 0.0033).
Interferometric results indicating a 'fast' rotation speed of ≈270 km s −1 are disputed in the spectroscopic analysis by Takeda et al. (2008), who favor a more 'moderate' speed of ≈175 km s −1 . However, using the parameters from Peterson et al. (2006b), Yoon et al. (2008) evaluate ELODIE archival spectra of Vega and invoke turbulence on large scales ('cyclones', even, to use their language) in their analysis to achieve a match between the data and their calculated spectral lines. The implications of such a result are intriguing, to say the least: these results are themselves rather dramatic, but the underlying approach of interferometric observations to guide spectroscopic analysis could be a far-reaching in application to more than just studies of rapidly rotating stars.
It is interesting to note that there has been no detection of X-rays by CHANDRA even after 29 ks of observing (Pease et al. 2006), even though (as noted above) they have been detected for Altair and Alderamin (Robrade and Schmitt 2009). The overall set of 'problems' with Vega are reviewed in the article by Gray (2007), of which the interferometrically detected phenomena of IR excess and rapid rotation play a large role. Ultimately, all of these various pathologies may lead astronomers to use stars other than Vega for photometric standards, as already proposed by Engelke et al. (2010).

Rasalhague (α Oph)
Rasalhague (α Oph A, HR6556, HD 159561; α = 17 h 34 m 56.07 s , δ = +12 • 33 m 36.13 s ), due to a number of unique circumstances, is a particularly interesting rapid rotator and is quickly becoming a fundamental laboratory for exploring physics of rapid rotation. As with many of these objects, it was originally observed interferometrically by Hanbury Brown et al. (1974).
An A5IV (Gray et al. 2001) δ Scuti variable star, α Oph A has been observed to be a nearly edge-on rapid rotator with ω/ω c ∼ 0.88 with the imaging study of  Zhao et al. (2009). As noted above in the discussion for Alderamin (Sect. 5.4), the edge-on orientation makes solving for β problematic, so Zhao et al. (2009) proceeded with a fixed value of β = 0.25. Monnier et al. (2010b) observed α Oph A with the MOST astroseismology satellite, noting that this inclination is particularly favorable (odd-parity l − |m| modes are suppressed, simplifying mode identifications) and detected rotationally modulated g-modes in the object. The (unexpected) linear relationships between g-mode spacings were provisionally explained as dispersion-free Kelvin waves. Hinkley et al. (2011) present an improved astrometric orbit of α Oph A and B from eight years of AO imaging data, and point out that a significantly improved determination of the 7.9-year astrometric orbit will be possible during the ∼50 mas April 2012 periastron, using optical interferometry. This in turn will permit mass determination of α Oph A at the few percent level. Collectively, these studies allow Rasalhague to be a prototype for probing the physics of stellar interiors, and understanding the effects of rapid rotation in that regime.

β Cas (Caph)
The most recent addition to the stable of rapid rotators observed with optical interferometry is β Cas . Its spectral type of F2III-IV (Rhee et al. 2007) indicates it is a 'retired' A-type star, evolving off the main sequence, and makes this observation a noteworthy one, in that it is the first object that is not (currently) 'A' or 'B' spectral type. As a cooler object, it represents an expansion of the tests derived from this technique into new areas of discovery space. The rotational velocity of ∼70 km s −1 (Rachford and Foight 2009) did not, on its face, give an a priori suggestion that it was an rapid rotator, although the small value of P / sin i reported therein (2.48 d ) certainly hinted at the possibility. Observationally, β Cas presents itself nearly pole-on, with i = 19.9 ± 1.9 • , indicating a significant rotation rate (v = 206 km s −1 , P = 1.12 +0.03 −0.04 days). As with Regulus, Che et al. (2011) also fit for β, finding a value of 0.146 +0.013 −0.007 . This lower value is consistent with trends toward a lower expected value for a star with lower surface temperature (see Fig. 9 of Che et al. 2011, andClaret 2000), although it is not nearly as small as would be expected for a fully radiative photosphere (β = 0.08).

Broad impact in astronomy
Although there is only a small number of results thus far, there are already broad and significant impacts throughout astronomy as the implications of such detailed new views of stars diffuse through the field. For example, the particulars of the analysis of the spectral energy distribution of Vega already noted above ) significantly revises our view of this fundamental standard (Rieke et al. 2008).
The interaction between spectroscopy and rapid rotator analysis has a long history, notably from the work by Collins (1963Collins ( , 1965, Collins and Harrington (1966). The interferometric observations have begun to guide spectroscopy-for example, discrepancies in T EFF indicated in Hα and Hβ Balmer line profiles in Altair found in Smalley et al. (2002) and Smalley and Kupka (2003) are resolved when the lower effective temperatures indicated by interferometry are used. Similarly, as predicted in Smalley et al. (2002), subsequent observations of Rasalhague also indicate a lower average T EFF , particularly given the edge-on orientation of this object, which resolve similar temperature anomalies. Such agreements are contrasted with the example of v sin i estimates for Vega, which even when account for pole-on rapid rotation (Hill et al. 2004) fall short of the true magnitude of the star's extreme rotational speed (Aufdenberg et al. 2007). Similarly, a line profile analysis of Altair incorporating not just the Balmer lines but a larger 650-line list (Reiners and Royer 2004) does not quite fall into agreement is the most recent values for inclination and rotation speed ). Analyses are now being forced to cast a wider net and consider such things as revised metallicities (Yoon et al. 2010) in order to reconcile spectroscopy and interferometry, a move which perhaps mirrors other upheavals in astronomy related to metallicity (Asplund et al. 2009).
Folding details of rapid rotation from these observations into spectroscopic modeling has not been prevented by these challenges, however. The line profiles of Huang and Gies (2006) are built using hemispheric averages for {T EFF , log g, v sin i} in their analysis using line profiles to investigate evolution of stellar rotation in young clusters. Neff and Simon (2008) examined O VI emission lines in detail in the context of the rapid rotation reported for Altair and Alderamin, to explore the onset of convection zones as made possible by rotational equatorial cooling-their detection of such emission in a variety of stars indicates that magnetic activity without substantial convective zones needs to be possible in current models. Khalack (2005) also point out that magnetic dipole strength can be overestimated if stellar oblateness is not taken into account. Combining interferometric and spectroscopic observations is suggested as a way to probe the properties of the rotation law of the external layers of rapidly rotating stars (Zorec et al. 2011).
Underlying spectroscopy, observational characterization of rapid rotation has also started to make significant inroads with stellar modeling (Domiciano de Souza et al. 2002). Comparison of the results of observed rapid rotators to stellar models is often difficult, in that typical models are constructed for non-rotating stars. Che et al. (2011) take the approach of applying corrections from predictions such as Sackmann (1970) to facilitate comparison to Y 2 models (Yi et al. 2001), with some success, although they note further work is warranted.
From an examination of collective results to date, Che et al. (2011) suggests that for general use in modeling rapidly rotating stars, a gravity darkening value of β = 0.19 is more appropriate. The interferometric results seem to validate the suggestion in Tassoul (2000) that solid-body rotation is impossible for a rapid rotator: temperature and pressure consistency across the stellar surface is disrupted by such rotation. This disruption leads to temperature and pressure gradients between the poles and equator, which in turn induces meridional circulation and breaks down strict radiative equilibrium. The differential angular momentum of the meridional circulation's matter flow-the higher latitude material carries less angular momentum-will then cause differential rotation.
Stellar rotation affects abundance, mass loss rates, and overall evolution (for a discussion of the implications for the most massive stars, see the 13-part series that be-gins with Meynet and Maeder 1997). New grids are in the process of being produced which incorporate the observational interferometric results, including smoother transition values for gravity darkening (Claret 2005). Alternatively, the basic parameters being discovered through interferometry are being used as anchor points in stellar evolution models to 'run the clock backwards', and explore the history of parameters such as angular momentum (van Belle et al. 2006).
Beyond of spectroscopy and modeling, there are a broad number of areas where these observational results are being considered as noteworthy. The naturally complementary techniques of asteroseismology and interferometry are seen as deeply connected on the point of rotation by Cunha et al. (2007) and Aerts et al. (2009); observation extension of these two techniques in a conjoined fashion has already been demonstrated in Monnier et al. (2010b) In considering transiting extrasolar planets, Seager and Hui (2002) noted the stellar oblateness measured by van Belle et al. (2001), and uses that as a starting point to discuss constraining the rotation rate of those planets through measurements of the planetary oblateness. Such measurements are expected to be possible with the high-precision data being produced by the Kepler mission. Barnes (2009) notes that the Rossiter-McLaughlin measurements that have been employed to date to probe the spin-orbit alignments of transiting systems are difficult about rapid rotators, but suggests that the unusual and distinctive transiting lightcurves associated with such systems could instead be utilized to explore those alignments.
Evaluation of rapid rotation is seen as a necessary part of understanding microlensing signals (Han and Chang 2006), with shapes of caustics being affected by source oblateness; Rattenbury et al. (2005) attempt to derive from microlensing event photometry the shape of the source star. Space climate is seen as directly connected to the irradiance properties of our sun, which in turn are connected to its shape (Lefebvre et al. 2007;Lefebvre et al. 2005). On a fundamental level, measurement of stellar oblateness is now acknowledged as one of the basic techniques for measurement of stellar rotation (Royer 2005).

Future target list
The seminal work of Hanbury Brown et al. (1974) providing angular diameter sizes on 32 bright, hot, nearby stars may be reviewed in considering the importance of expanding our knowledge base of rapidly rotating stars. The error estimates on the limb-darkened angular diameters presented in this paper ranged from 2%-14%, with a median value of 6.3%. Cross-referencing these results against the rotational velocities readily found in the catalog of Glebocki and Gnacinski (2005), we find 13 objects in common. We can estimate average linear radii R and mass M from the appropriate sections in Cox (2000), and thereby estimate oblateness using the prescription found in the appendix of van Belle et al. (2006).
What is instructive is to not consider the overall degree of oblateness predicted for these 13 objects, but to compare that degree of non-sphericity with the error estimates of spherical size given in Hanbury Brown et al. (1974). In two cases (α Aql, α Eri), the predicted oblateness is significantly in excess (∼2-3× greater) than the spherical size quoted error; for three objects (β Car, ζ Pup, Car) oblateness is comparable (∼0.9-1.0×) to the quoted error, with a fourth also being reasonably large in comparison (α Vir, at 0.6×). Better than 45% of the objects available for such a line of inquiry show strong motivation for developing a deeper understanding of their true size, particularly since the Hanbury Brown et al. (1974) values have propagated widely in the literature.
To date, only 7 targets have been investigated with optical/near-infrared interferometry, which may leave the impression that a substantially larger sample is beyond the capabilities of current facilities. This is simply not true. To probe the possibilities of this manner of study, we created a list of targets that are appealing to presentday facilities. Our assumptions were modest: angular sizes in excess of 0.50 mas (roughly the current limiting angular size of modern optical interferometers), and size ratios in excess of R b /R a − 1 > 0.05 (as dictated by the visibility precisions achievable with those facilities.) Starting with the rotational velocity catalog of Glebocki and Gnacinski (2005), we cross-referenced it against Hipparcos and 2MASS (Perryman et al. 1997;Cutri et al. 2003) to establish V and K magnitudes. From these values, using the rough angular size predictor found in van Belle (1999), we estimated angular sizes and cut those stars with <0.50 mas. Using the quoted spectral types of Hipparcos to estimate average linear radius R and mass M from the nominal values found in Cox (2000), we then were able to estimate oblateness, cutting for The resultant catalog of 354 objects is illustrated in Fig. 13, of which almost a third are 'Altair-class' rotators with R b /R a − 1 > 0.15; the complete catalog is provided in an Appendix. This catalog is, of course, incomplete, due to the imperfection of the assumptions involved, and it further reflects the incompleteness of Glebocki and Gnacinski (2005) (although the rather uniform sky distribution suggests this degree of incompleteness is low-e.g. no bias of northern hemisphere targets over those in the south). This approach also completely misses low-inclination rapid rotators with low values of v sin i, as illustrated in the case of Caph (Sect. 5.7). The limitations on use of v sin i as a predictive parameter for oblateness is clearly seen in Fig. 13: each of the targets that has been observed to date has an actual value in excess (in some cases, well in excess) of the prediction from v sin i. This of course illustrates the impact of the sin i term-and portends many fruitful future observations of oblateness with interferometry!

Conclusion
Direct observations of rapid rotators by long-baseline interferometers have matured rapidly over the past 10 years. Simple parametric models have given way to detailed models that probe the underlying stellar structure, and stunning images of photospheres that confirm the input physics in a dramatic, model-independent fashion (Fig. 14). The lessons learned and confidence collectively gained by the specialists carrying out these investigations is now being applied to a greater variety of objects. As the observed sample is expanded to objects of lower temperatures, some surprises are expected and have already been found . Matching this expansion of target selection has been a growth of instrumental capabilities. The 4-telescope CHARA-MIRC instrumentation that provided the first images ) is being expanded to full 6-way combination with sensitivity improvements through external fringe tracking (Monnier et al. 2010a); the first results from 4-telescope beam combination at VLTI with the PIONIER visitor instrument are very impressive (Berger et al. 2010); a new NOI combiner ('VI-SION') will add true 6-way combination to this facility as well and is slated to be on-sky in mid-2012. All of these developments bode well for further investigation into rapid rotators, possibly with complementary data sets spanning from the 500 nm through 2.4 µm, along with to increased 'snapshot' capability from the richer data sets.
What has been especially gratifying to see has been the rapid ingestion of these results into a wide variety of other endeavors in the field (Sect. 6). This has taken this activity from being just a 'cottage industry' to being one of deep implications for many of the foundational ideas in astronomy.