Stellar Photosphere Research Articles

The convective photosphere of the red supergiant CE Tauri

Context. Red supergiant stars are one of the latest stages in the evolution of massive stars. Their photospheric convection may play an important role in the launching mechanism of their mass loss; however, its characteristics and dynamics are still poorly constrained. Aims. By observing red supergiant stars with near infrared interferometry at different epochs, we expect to reveal the evolution of bright convective features on their stellar surface. Methods. We observed the M2Iab-Ib red supergiant star CE Tau with the VLTI/PIONIER instrument in the H band at two different epochs separated by one month. Results. We derive the angular diameter of the star and basic stellar parameters, and reconstruct two reliable images of its H-band photosphere. The contrast of the convective pattern of the reconstructed images is 5 ± 1% and 6 ± 1% for our two epochs of observation. Conclusions. The stellar photosphere shows few changes between the two epochs. The contrast of the convective pattern is below the average contrast variations obtained on 30 randomly chosen snapshots of the best matching 3D radiative hydrodynamics simulation: 23 ± 1% for the original simulation images and 16 ± 1% for the maps degraded to the reconstruction resolution. We offer two hypotheses to explain this observation. CE Tau may be experiencing a quiet convective activity episode or it could be a consequence of its warmer effective temperature (hence its smaller radius) compared to the simulation.

New Circumstellar Sources of PO and PN: The Increasing Role of Phosphorus Chemistry in Oxygen-rich Stars

Abstract PO and PN have been newly identified in several oxygen-rich circumstellar envelopes, using the Submillimeter Telescope of the Arizona Radio Observatory. The J = 5 → 4 and J = 6 → 5 transitions of PN near 235 and 282 GHz, and the lambda doublets originating in the J = 5.5 → 4.5 and J = 6.5 → 5.5 lines of PO at 240 and 284 GHz, have been detected toward the shells of asymptotic giant branch (AGB) stars TX Cam and R Cas. A similar set of lines has been observed toward the supergiant NML Cyg, and new transitions of these two molecules were also measured toward the AGB star IK Tau. Along with the previous data from VY Canis Majoris (VY CMa), these spectral lines were analyzed using the non-local thermodynamic equilibrium (non-LTE) circumstellar modeling code, ESCAPADE. For the AGB stars, peak abundances found for PN and PO were f ∼ (1–2) × 10−8 and (0.5–1) × 10−7, respectively, while those for the supergiants were f(PN) ∼ (0.3–0.7) × 10−8 and f(PO) ∼ (5–7) × 10−8. PN was well modeled with a spherical radial distribution, suggesting formation near the stellar photosphere, perhaps enhanced by shocks. PO was best reproduced by a shell model, indicating a photochemical origin, except for VY CMa. Overall, the abundance of PO is a factor of 5–20 greater than that of PN. This study suggests that phosphorus-bearing molecules are common in O-rich envelopes, and that a significant amount of phosphorus (>20%) remains in the gas phase.

A physical model of mass ejection in failed supernovae

During the core collapse of massive stars, the formation of the protoneutron star is accompanied by the emission of a significant amount of mass-energy ($\sim 0.3 \, M_{\odot}$) in the form of neutrinos. This mass-energy loss generates an outward-propagating pressure wave that steepens into a shock near the stellar surface, potentially powering a weak transient associated with an otherwise-failed supernova. We analytically investigate this mass-loss-induced wave generation and propagation. Heuristic arguments provide an accurate estimate of the amount of energy contained in the outgoing sound pulse. We then develop a general formalism for analyzing the response of the star to centrally concentrated mass loss in linear perturbation theory. To build intuition, we apply this formalism to polytropic stellar models, finding qualitative and quantitative agreement with simulations and heuristic arguments. We also apply our results to realistic pre-collapse massive star progenitors (both giants and compact stars). Our analytic results for the sound pulse energy, excitation radius, and steepening in the stellar envelope are in good agreement with full time-dependent hydrodynamic simulations. We show that {prior} to the sound pulses arrival at the stellar photosphere, the photosphere has already reached velocities $\sim 20-100 \%$ of the local sound speed, thus likely modestly decreasing the stellar effective temperature prior to the star disappearing. Our results provide important constraints on the physical properties and observational appearance of failed supernovae.

The Magellan Evolution of Galaxies Spectroscopic and Ultraviolet Reference Atlas (MegaSaura). II. Stacked Spectra

Abstract We stack the rest-frame ultraviolet spectra of N = 14 highly magnified gravitationally lensed galaxies at redshifts . The resulting new composite spans Å, with a peak signal-to-noise ratio (S/N) of 103 per spectral resolution element (∼100 km s−1). It is the highest S/N, highest spectral resolution composite spectrum of z ∼ 2–3 galaxies yet published. The composite reveals numerous weak nebular emission lines and stellar photospheric absorption lines that can serve as new physical diagnostics, particularly at high redshift with the James Webb Space Telescope (JWST). We report equivalent widths to aid in proposing for and interpreting JWST spectra. We examine the velocity profiles of strong absorption features in the composite, and in a matched composite of COS/HST galaxy spectra. We find remarkable similarity in the velocity profiles at and , suggesting that similar physical processes control the outflows across cosmic time. While the maximum outflow velocity depends strongly on ionization potential, the absorption-weighted mean velocity does not. As such, the bulk of the high-ionization absorption traces the low-ionization gas, with an additional blueshifted absorption tail extending to at least −2000 km s−1. We interpret this tail as arising from the stellar wind and photospheres of massive stars. Starburst99 models are able to replicate this high-velocity absorption tail. However, these theoretical models poorly reproduce several of the photospheric absorption features, indicating that improvements are needed to match observational constraints on the massive stellar content of star-forming galaxies at . We publicly release our composite spectra.

The First Post-Kepler Brightness Dips of KIC 8462852

Abstract We present a photometric detection of the first brightness dips of the unique variable star KIC 8462852 since the end of the Kepler space mission in 2013 May. Our regular photometric surveillance started in 2015 October, and a sequence of dipping began in 2017 May continuing on through the end of 2017, when the star was no longer visible from Earth. We distinguish four main 1%–2.5% dips, named “Elsie,” “Celeste,” “Skara Brae,” and “Angkor,” which persist on timescales from several days to weeks. Our main results so far are as follows: (i) there are no apparent changes of the stellar spectrum or polarization during the dips and (ii) the multiband photometry of the dips shows differential reddening favoring non-gray extinction. Therefore, our data are inconsistent with dip models that invoke optically thick material, but rather they are in-line with predictions for an occulter consisting primarily of ordinary dust, where much of the material must be optically thin with a size scale ≪1 μm, and may also be consistent with models invoking variations intrinsic to the stellar photosphere. Notably, our data do not place constraints on the color of the longer-term “secular” dimming, which may be caused by independent processes, or probe different regimes of a single process.

Transition Probabilities of Co ii Weak Lines to the Ground and Low Metastable Levels.

New branching fraction (BF) measurements based primarily on data from a cross-dispersed echelle spectrometer are reported for 84 lines of Co ii. The BFs for 82 lines are converted to absolute atomic transition probabilities using radiative lifetimes from laser-induced fluorescence (LIF) measurements on 19 upper levels of the lines. A lifetime of 3.3(2) ns for the z5D0 level is used based on LIF measurements for lifetimes of the four other levels in the z5D term. Twelve of the eighty-four lines are weak transitions connecting to the ground and low metastable levels of Co+. Another 46 lines are strong transitions connecting to the ground and low metastable levels of Co+. For these lines, log(gf) values were measured in earlier studies and, with a few exceptions, are confirmed in this study. Such lines, if unblended in stellar spectra, have the potential to yield Co abundance values unaffected by any breakdown of the local thermodynamic equilibrium approximation in stellar photospheres because the ground and low metastable levels of Co+ are the primary population reservoirs of Co in the photospheres of interest. Weak lines, if unblended, are useful in photospheres with high Co abundance, and strong lines are useful in metal-poor photospheres. New hyperfine structure A constants for 28 levels of ionized Co from least-squares fits to Fourier transform spectra line profiles are reported. These laboratory data are applied to re-determine the Co abundance in the metal-poor halo star HD 84937. BFs and transition probabilities for 19 lines are reported for the first time.

Observing the linked depletion of dust and CO gas at 0.1–10 au in disks of intermediate-mass stars

We report on the discovery of correlations between dust and CO gas tracers of the 0.1–10 au region in planet-forming disks around young intermediate-mass stars. The abundance of refractory elements on stellar photospheres decreases as the location of hot CO gas emission recedes to larger disk radii, and as the near-infrared excess emission from hot dust in the inner disk decreases. The linked behavior between these observables demonstrates that the recession of infrared CO emission to larger disk radii traces an inner disk region where dust is being depleted. We also find that Herbig disk cavities have either low (~5–10%) or high (~20–35%) near-infrared excess, a dichotomy that has not been captured by the classic definition of “pre-transitional” disks.

Orbital misalignment of the Neptune-mass exoplanet GJ 436b with the spin of its cool star.

The angle between the spin of a star and the orbital planes of its planets traces the history of the planetary system. Exoplanets orbiting close to cool stars are expected to be on circular, aligned orbits because of strong tidal interactions with the stellar convective envelope. Spin-orbit alignment can be measured when the planet transits its star, but such ground-based spectroscopic measurements are challenging for cool, slowly rotating stars. Here we report the three-dimensional characterization of the trajectory of an exoplanet around an M dwarf star, derived by mapping the spectrum of the stellar photosphere along the chord transited by the planet. We find that the eccentric orbit of the Neptune-mass exoplanet GJ 436b is nearly perpendicular to the stellar equator. Both eccentricity and misalignment, surprising around a cool star, can result from dynamical interactions (via Kozai migration) with a yet-undetected outer companion. This inward migration of GJ 436b could have triggered the atmospheric escape that now sustains its giant exosphere.

Upper stellar mass limit by radiative feedback at low-metallicities: metallicity and accretion rate dependence

We investigate the upper stellar mass limit set by radiative feedback by the forming star with various accretion rates and metallicities. To this end, we numerically solve the structures of both a protostar and its surrounding accretion envelope assuming a spherical symmetric and steady flow. The optical depth of the dust cocoon, a dusty part of the accretion envelope, differs among the direct light from the stellar photosphere and the diffuse light re-emitted as dust thermal emission. As a result, varying the metallicity qualitatively changes the way that the radiative feedback suppresses the accretion flow. With a fixed accretion rate of $10^{-3}M_{\odot} {\rm yr^{-1}}$, the both direct and diffuse lights jointly operate to prevent the mass accretion at $Z \gtrsim 10^{-1}Z_{\odot}$. At $Z \lesssim 10^{-1}Z_{\odot}$, the diffuse light is no longer effective, and the direct light solely limits the mass accretion. At $Z \lesssim 10^{-3}Z_{\odot}$, the HII region formation plays an important role in terminating the accretion. The resultant upper mass limit increases with decreasing metallicity, from a few $\times~10~M_\odot$ to $\sim 10^3~M_\odot$ over $Z = 1Z_{\odot} - 10^{-4}~Z_\odot$. We also illustrate how the radiation spectrum of massive star-forming cores changes with decreasing metallicity. First, the peak wavelength of the spectrum, which is located around $30 \mu {\rm m}$ at $1Z_{\odot}$, shifts to $< 3 \mu {\rm m}$ at $Z \lesssim 0.1 Z_\odot$. Second, a characteristic feature at $10 \mu \rm{m}$ due to the amorphous silicate band appears as a dip at $1Z_{\odot}$, but changes to a bump at $Z \lesssim 0.1 Z_{\odot}$. Using these spectral signatures, we can search massive accreting protostars in nearby low-metallicity environments with up-coming observations.

Redshifted X-rays from the material accreting onto TW Hydrae: Evidence of a low-latitude accretion spot

High resolution spectroscopy, providing constraints on plasma motions and temperatures, is a powerful means to investigate the structure of accretion streams in CTTS. In particular, the accretion shock region, where the accreting material is heated to temperatures of a few MK as it continues its inward bulk motion, can be probed by X-ray spectroscopy. To attempt to detect for the first time the motion of this X-ray-emitting post-shock material, we searched for a Doppler shift in the deep Chandra/HETGS observation of the CTTS TW Hya. This test should unveil the nature of this X-ray emitting plasma component in CTTS, and constrain the accretion stream geometry. We searched for a Doppler shift in the X-ray emission from TW Hya with two different methods, by measuring the position of a selected sample of emission lines, and by fitting the whole TW Hya X-ray spectrum, allowing the line-of-sight velocity to vary. We found that the plasma at T~2-4 MK has a line-of-sight velocity of 38.3+/-5.1 km/s with respect to the stellar photosphere. This result definitively confirms that this X-ray-emitting material originates in the post-shock region, at the base of the accretion stream, and not in coronal structures. The comparison of the observed velocity along the line of sight, 38.3+/-5.1 km/s, with the inferred intrinsic velocity of the post shock of TW Hya, v_post~110-120 km/s, indicates that the footpoints of the accretion streams on TW Hya are located at low latitudes on the stellar surface. Our results indicate that complex magnetic field geometries, such as that of TW Hya, permit low-latitude accretion spots. Moreover, since on TW Hya the redshift of the soft X-ray emission is very similar to that of the narrow component of the CIV resonance doublet at 1550 Ang, as found by Ardila et al. (2013), then the plasma at 2-4 MK and that at 0.1 MK likely originate in the same post-shock regions.

A differential least-squares deconvolution method for high precision spectroscopy of stars and exoplanets – I. Application to obliquity measurements of HARPS observations of HD189733b

High precision measurements of stellar spectroscopic line profiles and their changes over time contain very valuable information about the physics of the stellar photosphere (stellar activity) and can be used to characterize extrasolar planets via the Rossiter–McLaughlin effect or from reflected light from the planet. In this paper, we present a new method for measuring small changes in the mean line profile of a spectrum by performing what we call differential least-squares deconvolution (dLSD). The method consists in finding the convolution function (or kernel) required to transform a high-signal-to-noise ratio template of the star into each observed spectrum. Compared to similar techniques, the method presented here does not require any assumptions on the template spectrum (e.g. no line-list or cross-correlation mask required). We show that our implementation of dLSD is able to perform at least as good as other techniques by applying it to star–planet obliquity measurements of exoplanet HD183799 during its transit.1 Among other things, the method should enable model-independent detection of light reflected by an exoplanet.

Forecasting the Impact of Stellar Activity on Transiting Exoplanet Spectra

Abstract Exoplanet host star activity, in the form of unocculted starspots or faculae, alters the observed transmission and emission spectra of the exoplanet. This effect can be exacerbated when combining data from different epochs if the stellar photosphere varies between observations due to activity. Here, we present a method to characterize and correct for relative changes due to stellar activity by exploiting multi-epoch ( visits/transits) observations to place them in a consistent reference frame. Using measurements from portions of the planet’s orbit where negligible planet transmission or emission can be assumed, we determine changes to the stellar spectral amplitude. With the analytical methods described here, we predict the impact of stellar variability on transit observations. Supplementing these forecasts with Kepler-measured stellar variabilities for F-, G-, K-, and M-dwarfs, and predicted transit precisions by the James Webb Space Telescope’s (JWST) NIRISS, NIRCam, and MIRI, we conclude that stellar activity does not impact infrared transiting exoplanet observations of most presently known or predicted TESS targets by current or near-future platforms, such as JWST, as activity-induced spectral changes are below the measurement precision.

Can we really use chemical properties of red-giant stars as age indicators?

AbstractThe cornerstone mission of the European Space Agency, Gaia, together with forthcoming complementary surveys (CoRoT, Kepler, K2, APOGEE and Gaia-ESO), will revolutionize our understanding of the formation and history of our Galaxy, providing accurate stellar masses, radii, ages, distances, as well as chemical properties for a very large sample of stars across different Galactic stellar populations. Using improved population synthesis approach and new stellar evolution models we attempt to evaluate the possibility to derive ages of clump stars from their chemical properties. A new version of the Besançon Galaxy models (BGM) is used in which new stellar evolutionary tracks are computed from the stellar evolution code STAREVOL. The effects of mixing on chemical composition of the stellar photosphere has a significant impact on the determined stellar age from the observed [C/N] ratio. We clearly show that transport processes occurring in red-giant stars should be taken into account in the determination of ages for future Galactic archaeology studies.

The Effects of Magnetic Fields on the Dynamics of Radiation Pressure–dominated Massive Star Envelopes

Abstract We use three-dimensional radiation magnetohydrodynamic simulations to study the effects of magnetic fields on the energy transport and structure of radiation pressure–dominated main sequence massive star envelopes at the region of the iron opacity peak. We focus on the regime where the local thermal timescale is shorter than the dynamical timescale, corresponding to inefficient convective energy transport. We begin with initially weak magnetic fields relative to the thermal pressure, from 100 to 1000 G in differing geometries. The unstable density inversion amplifies the magnetic field, increasing the magnetic energy density to values close to equipartition with the turbulent kinetic energy density. By providing pressure support, the magnetic field’s presence significantly increases the density fluctuations in the turbulent envelope, thereby enhancing the radiative energy transport by allowing photons to diffuse out through low-density regions. Magnetic buoyancy brings small-scale magnetic fields to the photosphere and increases the vertical energy transport, with the energy advection velocity proportional to the Alfvén velocity, although in all cases we study, photon diffusion still dominates the energy transport. The increased radiative and advective energy transport causes the stellar envelope to shrink by several scale heights. We also find larger turbulent velocity fluctuations compared with the purely hydrodynamic case, reaching at the stellar photosphere. The photosphere also shows vertical oscillations with similar averaged velocities and periods of a few hours. The increased turbulent velocity and oscillations will have strong impacts on the line broadening and periodic signals in massive stars.

SONS: The JCMT legacy survey of debris discs in the submillimetre

Debris discs are evidence of the ongoing destructive collisions between planetesimals, and their presence around stars also suggests that planets exist in these systems. In this paper, we present submillimetre images of the thermal emission from debris discs that formed the SCUBA-2 Observations of Nearby Stars (SONS) survey, one of seven legacy surveys undertaken on the James Clerk Maxwell telescope between 2012 and 2015. The overall results of the survey are presented in the form of 850 microns (and 450 microns, where possible) images and fluxes for the observed fields. Excess thermal emission, over that expected from the stellar photosphere, is detected around 49 stars out of the 100 observed fields. The discs are characterised in terms of their flux density, size (radial distribution of the dust) and derived dust properties from their spectral energy distributions. The results show discs over a range of sizes, typically 1-10 times the diameter of the Edgeworth-Kuiper Belt in our Solar System. The mass of a disc, for particles up to a few millimetres in size, is uniquely obtainable with submillimetre observations and this quantity is presented as a function of the host stars' age, showing a tentative decline in mass with age. Having doubled the number of imaged discs at submillimetre wavelengths from ground-based, single dish telescope observations, one of the key legacy products from the SONS survey is to provide a comprehensive target list to observe at high angular resolution using submillimetre/millimetre interferometers (e.g., ALMA, SMA).

VLTI/PIONIER images the Achernar disk swell

Context.The mechanism of disk formation around fast-rotating Be stars is not well understood. In particular, it is not clear which mechanisms operate, in addition to fast rotation, to produce the observed variable ejection of matter. The star Achernar is a privileged laboratory to probe these additional mechanisms because it is close, presents B ⇌ Be phase variations on timescales ranging from ~6 yr to ~15 yr, a companion star was discovered around it, and probably presents a polar wind or jet.Aims.Despite all these previous studies, the disk around Achernar was never directly imaged. Therefore we seek to produce an image of the photosphere and close environment of the star.Methods.We used infrared long-baseline interferometry with the PIONIER instrument at the Very Large Telescope Interferometer (VLTI) to produce reconstructed images of the photosphere and close environment of the star over four years of observations. To study the disk formation, we compared the observations and reconstructed images to previously computed models of both the stellar photosphere alone (normal B phase) and the star presenting a circumstellar disk (Be phase).Results.The observations taken in 2011 and 2012, during the quiescent phase of Achernar, do not exhibit a disk at the detection limit of the instrument. In 2014, on the other hand, a disk was already formed and our reconstructed image reveals an extendedH-band continuum excess flux. Our results from interferometric imaging are also supported by several Hαline profiles showing that Achernar started an emission-line phase sometime in the beginning of 2013. The analysis of our reconstructed images shows that the 2014 near-IR flux extends to ~1.7–2.3 equatorial radii. Our model-independent size estimation of theH-band continuum contribution is compatible with the presence of a circumstellar disk, which is in good agreement with predictions from Be-disk models.

K2 Ultracool Dwarfs Survey. I. Photometry of an L Dwarf Superflare

Abstract We report on K2 Campaign 8 measurements of a huge white light flare on the L1 dwarf SDSSp J005406.55-003101.8 (EPIC 220186653). The source is a typical L1 dwarf at a distance of ∼50 pc, probably an old hydrogen-burning star rather than a young brown dwarf. In the long (30-minute) cadence photometry, the flare peak is 21 times the flux of the stellar photosphere in the broad optical Kepler filter, which we estimate corresponds to ΔV ≈ −7.1. The total equivalent duration of the flare is 15.4 hr. We estimate that the total bolometric energy of the flare was 4 × 1033 erg, more powerful than the previously reported Kepler white light flares for the L1 dwarf WISEP J190648.47+401106.8, but weaker than the ΔV = −11 L0 dwarf superflare ASASSN-16ae. The initial (impulsive) cooling phase is too rapid to resolve with our 30-minute cadence data, but after 1 hour the gradual cooling phase has an exponential time constant of 1.8 hr. We use template fitting to estimate that the full time-width-at-half-amplitude of the light curve is &lt;10 minutes and that the true flare maximum reached ∼70 times the stellar photosphere, or ΔV ≈ −8. This flare is comparable to the most powerful Kepler flares observed on the active M4 dwarf GJ 1243.

An observational study of dust nucleation in Mira (o Ceti)

The formation of silicate dust in oxygen-rich envelopes of evolved stars is thought to be initiated by formation of seed particles that can withstand the high temperatures close to the stellar photosphere and act as condensation cores farther away from the star. Candidate species considered as first condensates are TiO and TiO$_2$. We aim to identify and characterize the circumstellar gas-phase chemistry of titanium that leads to the formation of solid titanium compounds in the envelope of $o$ Cet, the prototypical Mira, and seek an observational verification of whether titanium oxides play a major role in the onset of dust formation in M-type AGB stars. We present high angular-resolution ALMA observations at submillimeter (submm) wavelengths supplemented by APEX and Herschel spectra of the rotational features of TiO and TiO$_2$. In addition, circumstellar features of TiO and TiI are identified in optical spectra which cover multiple pulsation cycles of $o$ Cet. The submm ALMA data reveal TiO and TiO$_2$ bearing gas within the extended atmosphere of Mira. While TiO is traceable up to a radius (FWHM/2) of 4.0 R$_{\star}$, TiO$_2$ extends as far as 5.5 R$_{\star}$ and unlike TiO appears to be anisotropically distributed. Optical spectra display variable emission of TiI and TiO from inner parts of the extended atmosphere (<3 R$_{\star}$). Chemical models which include shocks are in general agreement with the observations of gas-phase titanium-bearing molecules. It is unlikely that substantial amounts of titanium is locked up in solids because the abundance of the gaseous titanium species is very high. In particular, formation of hot titanium-rich condensates is very improbable because we find no traces of their hot precursor species in the gas phase. It therefore appears unlikely that the formation of dust in Mira, and possibly other M-type AGB stars, is initiated by titanium oxides.

The First 40 Million Years of Circumstellar Disk Evolution: The Signature of Terrestrial Planet Formation

Abstract We characterize the first 40 Myr of evolution of circumstellar disks through a unified study of the infrared properties of members of young clusters and associations with ages from 2 Myr up to ∼40 Myr: NGC 1333, NGC 1960, NGC 2232, NGC 2244, NGC 2362, NGC 2547, IC 348, IC 2395, IC 4665, Chamaeleon I, Orion OB1a and OB1b, Taurus, the β Pictoris Moving Group, ρ Ophiuchi, and the associations of Argus, Carina, Columba, Scorpius–Centaurus, and Tucana–Horologium. Our work features: (1) a filtering technique to flag noisy backgrounds; (2) a method based on the probability distribution of deflections, P(D), to obtain statistically valid photometry for faint sources; and (3) use of the evolutionary trend of transitional disks to constrain the overall behavior of bright disks. We find that the fraction of disks three or more times brighter than the stellar photospheres at 24 μm decays relatively slowly initially and then much more rapidly by ∼10 Myr. However, there is a continuing component until ∼35 Myr, probably due primarily to massive clouds of debris generated in giant impacts during the oligarchic/chaotic growth phases of terrestrial planets. If the contribution from primordial disks is excluded, the evolution of the incidence of these oligarchic/chaotic debris disks can be described empirically by a log-normal function with the peak at 12–20 Myr, including ∼13% of the original population, and with a post-peak mean duration of 10–20 Myr.

ACCESS I. AN OPTICAL TRANSMISSION SPECTRUM OF GJ 1214b REVEALS A HETEROGENEOUS STELLAR PHOTOSPHERE

ABSTRACT GJ 1214b is the most studied sub-Neptune exoplanet to date. Recent measurements have shown its near-infrared transmission spectrum to be flat, pointing to a high-altitude opacity source in the exoplanet's atmosphere, either equilibrium condensate clouds or photochemical hazes. Many photometric observations have been reported in the optical by different groups, though simultaneous measurements spanning the entire optical regime are lacking. We present an optical transmission spectrum (4500–9260 Å) of GJ 1214b in 14 bins, measured with Magellan/IMACS repeatedly over three transits. We measure a mean planet-to-star radius ratio of R p / R s = 0.1146 ± 2 × 10 − 4 and mean uncertainty of σ ( R p / R s ) = 8.7 × 10 − 4 in the spectral bins. The optical transit depths are shallower on average than observed in the near-infrared. We present a model for jointly incorporating the effects of a composite photosphere and atmospheric transmission through the exoplanet's limb (the CPAT model), and use it to examine the cases of absorber and temperature heterogeneities in the stellar photosphere. We find the optical and near-infrared measurements are best explained by the combination of (1) photochemical haze in the exoplanetary atmosphere with a mode particle size r = 0.1 μm and haze-forming efficiency f haze = 10 % and (2) faculae in the unocculted stellar disk with a temperature contrast Δ T = 354 − 46 + 46 K, assuming 3.2% surface coverage. The CPAT model can be used to assess potential contributions of heterogeneous stellar photospheres to observations of exoplanet transmission spectra, which will be important for searches for spectral features in the optical.