Abstract
CaII K-line (3933 Å) absorption profiles observed towards 15 A-type and two B-type stars with known (or suspected) debris disks, in addition to archival spectral data for three A-type stars, are presented. Inspection of the extracted CaII K-line absorption lines has determined that one late B-type (HD 58647) and four new A-type (HD 56537, HD 64145, HD 108767, and HD 109573) stellar systems exhibit short-term (night-to-night) absorption variation within these profiles. This variability is due to the liberation of gas from falling evaporating bodies (the so-called FEBs activity) on the grazing approach of exocomets towards their parent star. Also we present archival spectra of the FeI λ3860 Å and CaII K-lines towards β Pictoris (HD 39060) that, for the first time, reveal the absorption signatures of FEB activity at similar velocities in both line profiles. FEB-hosting stars seem to be of an earlier type (on average spectral type Al.8 as compared with A3.1) than those where no FEB activity has been observed. This could be due to stellar activity levels (such as chromospheric activity or nonradial pulsations) being of a higher level within the atmospheres of these hotter A-type stars.
1. Introduction
The existence of reservoirs of small rocky bodies (i.e., asteroids/planetesimals) in orbits around many stellar systems is now well established. Their presence is required by current planetary formation theories [1, 2] and also through direct and indirect observations of stars with associated dusty and gaseous circumstellar disks [3–5]. Their appearance in large numbers close to the host star is likely episodic and triggered by the formation of large planets [6]. Since dust and gas-rich systems will transition to optically thin and gas-free debris disks, observation of asteroids/planetesimals is best performed in young or debris-rich systems, which may have associated planetary bodies in the early stages of formation. In such systems, the gravitation fields can potentially enable the disruption of large numbers of these kilometer-sized icy bodies into trajectories directed towards the young central star. Indeed the detection of evaporating gas liberated from the comae of swarms of these infalling comet-like bodies (or exocomets) can be interpreted as being a signature of the early stages of the formation of a planetary system. To date, the search for exocomets has focused on young (<50 Myr) A-type stars, since roughly one-third of normal A-type main sequence stars exhibit infrared (IR) excesses due to the presence of dusty debris disks [7, 8]. Since the radiation pressure from A-type stars removes gas and dust particles from the disk on timescales ≪1 Myr [9], then any dust and/or gas found close to the star must have arisen recently from the collision and evaporation events amongst planetesimals.
The evaporation of gaseous products from infalling exocomets occurring close to a star can potentially cause small disruptions in the ambient circumstellar disk plasma. For circumstellar disks that are viewed “edge-on” this evaporating material may be directly observed through transient gas absorption features seen at rapidly changing redshifts. This short-term (night-to-night) variability in absorption was first detected towards the A6V star β Pictoris, through repeated observations of the CaII K-line at 3933 Å [10, 11]. The relatively warm ( K–2500 K) gas liberated during these evaporative events generally lies within 1 AU of the star, whereas the cold and optically thick circumstellar gas disk (as traced by CO emission) is located at further distances >20 AU from the central star [12].
Over the last few years renewed interest in the detection of exocomet activity through observation of short-term (night-to-night) absorption variability of the circumstellar CaII 3933 Å K-line towards A-type stars has resulted in the addition of over a dozen new stellar systems thought to be somewhat similar to that of β Pictoris [3, 14, 15]. The release of evaporating gas from these kilometer-sized “falling evaporating bodies” (FEBs or exocomets) generally requires their orbits to be perturbed by an outer planet onto star-grazing orbits [16]. For the case of β Pictoris this prompted a twenty-year search for an associated exoplanet, which culminated in the detection of the giant planet, β Pictoris b, by Lagrange et al. [17]. The widely used techniques of exoplanet detection towards (older) stars of late spectral types G to M, via radial velocity variation or the observation of transits, are generally ill-suited to young A-type stars. Hence, the detection of short-term exocomet activity in such young (<50 Myr old) systems may currently be the best indication of the presence of associated exoplanets, whose existence can then be confirmed using very high contrast direct imaging techniques. This is of particular importance since recent high contrast imaging observations suggest that fewer than 10% of B- and A-type stars may possess giant exoplanets (of mass > 10 ) at distances between 40 and 650 AU [18].
An analysis of ~30 rapidly rotating A-type stars that have been repeatedly observed in absorption at the CaII K-line suggests that those (ten) stars that exhibit FEB activity are younger (<100 Myr) and are more likely to exhibit chemical peculiarities than those A-type stars in which FEB activity has not been detected [3]. Interestingly, FEB activity does not appear to be directly associated with a strong mid-IR excess. This is most probably linked to the viewing angle of the circumstellar gas disk. Even if a dust disk or infalling comets are present, they will be very difficult to detect if the viewing angle is not within a few degrees of edge-on [16]. We note that β Pictoris, the system in which FEB activity is most frequently observed, has a disk with an inclination angle of 88.5° [29]. In addition, the orbital inclination of the inflating comets is also a prime factor in determining the observed velocity of an FEB event, and for the case of β Pictoris the angle between the axis of the cometary orbit and the line of sight is tightly constrained to be within 10° of l50° [16]. In the absence of direct imaging, the only way we have of selecting near edge-on disk systems is by observing stars that rapidly rotate ( km s−1). Given that , an imperfect proxy for viewing angle, it is perhaps not surprising that the (10) stars in which FEB activity has thus far been detected are not characterized by significantly higher rotation velocities when compared with those (21) stars in which FEB activity was not detected by Welsh and Montgomery [3].
β Pictoris, an A6V star, has subsequently become the de facto canonical young (~12 Myr) A-type star with an IR excess and a debris disk containing planetesimal belts containing myriads of small rocky bodies [6, 30]. In addition to showing absorption variability at the circumstellar CaII K-line on a night-to-night (and sometimes hourly) basis, β Pictoris also exhibits absorption from the metastable 3d levels of the CaII IR triplet at 8542 Å [31, 32]. Such absorption, which is due to radiative pumping, can only occur within a circumstellar environment and cannotbe of interstellar origin. These observations indicate the location of the warm circumstellar gas around the central star to be <10 AU. Such infalling gas would appear to possess a high electron density ( cm−3) and also a high electron temperature ( K). β Pictoris also exhibits two other types of circumstellar absorption behavior at visible wavelengths. Firstly, the main circumstellar CaII K-line component arises close to (or at) the stellar radial velocity (i.e., km s−1) which is located at a slight red shift of km s−1 from that of the far weaker circumstellar NaI D-line absorption [33]. Fernández et al. [34] have explained this velocity difference by a gas braking mechanism in which neutral circumstellar species are rapidly ionized by the strong stellar radiation force and thus only reach limited (and slightly different) velocities.
Typical column density ratios of N(CaII)/N(NaI) > 30 : 1 have been determined for the gas at the central velocity of the CaII K-line. These values are at least an order of magnitude higher than those found for the general interstellar medium. A second phenomenon that seems peculiar only to the β Pictoris system is the unusually high level of circumstellar absorption associated with the rarely observed FeI 3860 Å line. This spectral feature is generally found to be very weak in the interstellar medium, being found mainly in cold and dense regions [35]. However, for the case of β Pictoris, it has an equivalent width of ~6 mÅ with the circumstellar absorption being centered near the radial velocity of the star (see Section 3.1.8 of this paper).
Now since FEB activity has been detected towards several other nearly A-type stars, it is instructive to pose the question as to whether β Pictoris really is a typical example of an A-type stellar system with a debris disk that exhibits FEB activity due to the evaporation of exocomets. We note that, of 16 debris discs around A-type stars that have been resolved by imaging in the IR, eleven show evidence of multiple dust components [36]. These dust disks exhibit a wide diversity in physical appearance, with radii from a few AU to hundreds of AU. Some of these disks are confined dust rings with clear inner and outer edges as found around HR 4796 by Schneider et al. [37], whereas some have very extended disks like that of β Pictoris [38]. Additionally, several authors have posited that the origin of inner dust (and gas) disks around A-type stars might be associated with a stellar phenomenon, even though their mass-loss rates are generally very low and should not give rise to significant near-infrared emission. However, we note that recent work has shown a possible correlation between the K-band IR excess and the rotational velocity of such stars, which would hint at a possible stellar mass-loss origin for the circumstellar dust (and gas) [39].
Other spectral types of young A-type stars have also been found to exhibit temporal absorption variability of the CaII K-line that can be associated with a surrounding circumstellar disk. These are the young (but metal poor) λ Bootis-type stars that are surrounded by a thick cloud of dusty gas [40, 41]. In addition, the premain sequence Herbig Ae/Be stars typically also have large IR excesses with associated circumstellar gas/dust disks [42–44]. Since it is widely accepted that circumstellar disks are a common feature of both star formation and planetary formation processes, then comparing the variability behavior of circumstellar disk gas surrounding all types of young A-type stars of differing ages (1–100 Myr) may provide insights into the early stages of the evolution of planetary bodies.
In this paper we present nightly spectral observations of 15 A-type stars (plus two late B-type stars) with known (or suspected) debris disks, in addition to presenting similar spectral data for three A-type stars retrieved from the European Southern Observatory data archive. In Section 2 we describe how we have derived the residual intensity absorption profiles of the circumstellar CaII K-line from each observation to measure the corresponding line equivalent width value. These measurements are discussed in Section 3 for each target star, and we have determined that one late B-type (HD 58647) and four new A-type (HD 56537, HD 64145, HD 108767, and HD 109573) stars exhibit short-term (night-to-night) variation in the equivalent width values of their circumstellar CaII K-line absorption line profiles. Based on the similar behavior of CaII absorption exhibited by the β Pictoris stellar system, we associate this absorption variability with the evaporation of gas from comet-like bodies on their grazing approach to a parent star and note that this brings the total number of systems thought to harbor exocomets to 16 (see Table 4), 13 of which have been discovered by the present authors.
2. Observations and Data Extraction
High spectral resolution observations of a total of 17 stars listed in Table 1 are presented in this paper. With the exception of the radial velocities [25], the astrophysical data for each star shown in Table 1 have been drawn from the SIMBAD data retrieval system. Target selection was based on the physical properties shared with other circumstellar gas systems already known to exhibit variable absorption at the CaII wavelengths. Thus, most of the program targets are similarly young, nearby ( < 100 pc) A-type main sequence stars that rotate rapidly and exhibit a mid-infrared excess. Such A-type stars often possess a known (or suspected) associated debris disk [26, 45], which may have the potential to harbor exocomet activity.
Spectral observations of the four stars listed in the uppermost section of Table 1 were performed in the service mode by staff on the 2.1 m Bernard Lyot telescope at the Pic du Midi Observatory in France. Spectra for three stars were recorded on four nights over the time period 10-28- to 11-1-2014, and the star HD 184006 was similarly observed on three occasions over the time period 09-21- to 10-7-2013. The intent of these (and all of our) observations was to observe each of the targets (weather permitting) on a nightly basis to search for absorption variability of the CaII K-line at 3933 Å. Each observation consisted typically of an integration time of ~60 minutes, such that the spectra were well exposed with typical ratios of ~200 : 1 in the stellar continuum around the CaII K-line. Medium spectral resolution observations of (4.6 km s−1) were obtained for the 3 stars over the entire 3705–9485 Å range using the EsPaDOns echelle spectrograph used in the “star only” mode. Spectra were recorded with a EEV CCD of 2000 × 4500 pixels and the data were processed at the Observatory using the Libre ESpRIT data reduction software tool, which uses many of the same data extraction techniques outlined in the last paragraph of this section.
In the second section of Table 1 we present information regarding repeated spectroscopic observations of 11 A- or B-type stars using the Sandiford Echelle Spectrograph mounted at the Cassegrain focus of the 2.1 m telescope at the McDonald Observatory, Texas. Medium spectral resolution observations of 60,000 (5 km s−1) were obtained over the 3800–5000 Å region and the data was recorded with a Reticon CCD detector of 1200 × 400 pixels. The first nine of these stars were observed over the nights of 03-1- to 03-7-2013. The tenth star, HD 158352, was observed in an identical fashion on four occasions over the time period 08-18- to 08-23-2010. We also present CaII K-line absorption observations of the star HD 9672 (49 Ceti) recorded on the nights of 08-11-12-2013 using the same instrumental set-up. This star was previously known to exhibit FEB absorption activity similar to that of β Pictoris [15] and these observations were chosen to simultaneously coincide with ultraviolet absorption observations recorded with the STIS instrument aboard the Hubble Space Telescope. Preliminary results from observations of the neutral ultraviolet absorption species detected towards 49 Ceti have been presented by Roberge et al. [46].
Finally spectral observations of the two stars listed in the lowest section of Table 1 were performed during the nights of 05-8- to 05-13-2014 using the GIRAFFE fiber-fed echelle spectrograph mounted on the 1.9 m Radcliffe telescope of the South African Astronomical Observatory. The spectral data which encompassed a wavelength range of 3850 Å–5500 Å were recorded with a Tektronix 1024 × 1024 CCD detector at a resolving power of ~40,000 (7.5 km s−1).
All of the spectral data recorded at the various observatories were extracted and reduced in a similar manner. The recorded images were bias-subtracted and flat-fielded using standard IRAF echelle data reduction routines. Cosmic rays and bad pixels were removed, and the spectral orders were then extracted using the standard optimal spectrum extraction algorithm of Horne [47]. The wavelength dispersion calibration of the stellar spectra was obtained by cross-reference to Th-Ar emission spectra recorded at the beginning and end of each night. This resulted in a wavelength accuracy of ~0.015 Å (1 km s−1) for all of the stellar spectra. These wavelengths were then transformed into the heliocentric frame of reference for all future discussion in this paper.
2.1. Archival Data
We have also made a cursory inspection of some of the numerous observations of A-type stars recorded with the FEROS and/or HARPS high-resolution spectrographs that are deposited within the European Southern Observatory (ESO) data archive (http://archive.eso.org/cms.html). In Table 2 we list three A-type stars (together with their astrophysical information) for which we have downloaded spectral data that have been processed using the standard ESO data pipeline. We list the dates of each of these observations together with the spectrograph that recorded these spectra. We note that the HARPS spectral data are recorded with a resolving power of , whereas the FEROS spectral data have a lower resolving power of . In addition to the CaII and NaI absorption lines, the FeI 3860 Å circumstellar absorption line was also extracted from the β Pictoris spectra. None of the other targets listed in Tables 1 and 2 had any significant absorption feature in their spectra near this wavelength to a 3σ upper limit of EW < 1.5 mÅ.
2.2. Absorption Line Analysis
All of the observed spectra contained a rotationally broadened stellar (CaII-K) absorption line centered at 3933 Å, most of which also had a very narrow circumstellar absorption line superposed at the bottom of the stellar absorption trough. A typical example of this type of absorption profile can be seen in Figure 1 of Montgomery and Welsh [15]. This sort of spectral signature is typical of A-type stars that possess circumstellar gaseous debris disks corotating near to their stellar radial velocity [3, 15, 27, 28].
A stellar continuum level for each target was determined by fitting the broad CaII K-line photospheric absorption trough with a 6th (or higher) order polynomial over a velocity range of ±150 km s−1 from the central absorption. This resulted in a residual intensity absorption profile for the far weaker circumstellar CaII K-line. The absorption equivalent width of this circumstellar line was then determined from a straightforward integration of the residual intensity spectrum summed over the line profile. The continuum placement software assigns an rms error to each of the spectral data points, which is subsequently adopted as the error for these points as defined in Vallerga et al. [48].
All of the spectra of our target stars may have some (small) contribution due to interstellar absorption along their sightline path lengths. Hence, in Table 3 we also list the projected velocity vector of the local interstellar cloud (LIC) for each sightline using the respective cloud properties listed in Table of Redfield and Linsky [13]. In many cases these velocities can be used to discriminate between circumstellar absorption features (which should form at velocities close to the radial velocity of the star) and any local interstellar features.
3. Discussion of Individual Target Spectra
Previous observations of several local A-type stars with debris disks have shown two main types of circumstellar absorption variability. First, the absorption strength or shape of the absorption profile of the main circumstellar CaII K-line (usually centered at the radial velocity of the star) may change over time periods on the order of days (or in some cases on an hourly basis). Secondly, transient weak absorption features (FEBs) may occasionally appear at velocities that are red or blue-shifted by tens of km s−1 away from the main circumstellar K-line, which itself may often vary in absorption strength. An assessment of a significant change in the main circumstellar line can be obtained by comparison of its night-to-night equivalent width (EW) values, as integrated over the same velocity interval.
In the following sections we discuss the types of CaII absorption variability observed towards the 20 stars listed in Tables 1 and 2.
3.1. Stars Exhibiting CaII Circumstellar Variability or the Presence of FEBs
3.1.1. HD 64145
HD 64145 (λ Gem) is an A3V type star (of, as yet, undetermined age) located at a distance of 68 pc with a warm ( K) debris disk located close to the central star that exhibits a small IR excess at 18 μm [49]. In Figure 1 we see that each of the 4 nights’ observations of this star has revealed an absorption component at km s−1, which is close to that of the stellar radial velocity. We note that the projected velocity of the local interstellar cloud according to Redfield and Linsky [13] in this galactic direction is +18.5 km s−1, and thus we identify the nightly absorption at km s−1 as being of circumstellar origin. The measured EW of this line was ~3.3 mÅ on three of the four observations, whereas the observation of 10-29-2014 shows the absorption strength to increase to EW = 4.6 mÅ which (on inspection of the associated measurement errors) is change in value. On the night of 10-31-2014 an extra absorption feature appeared at km s−1 with an EW = mÅ and on night 10-30-2014 a weaker absorption feature at km s−1 is present with an EW = mÅ. We identify the extra absorption on 10-31-2014 as being due to an FEB event due to the evaporation of exocomet gas, and although the extra absorption seen on 10-30-2014 is only just above the continuum noise level (and is of lower statistical significance), this may also be of an FEB origin.
For comparison purposes, the main circumstellar CaII K-line detected at km s−1 towards β Pictoris typically has an EW > 100 mÅ. In addition, its FEB absorption components appear on a regular (< nightly) basis at velocities up to +150 km s−1 that have an EW > 20 mÅ [50]. See Section 3.1.8 for a more detailed discussion of exocomet activity in the canonical β Pictoris system.
3.1.2. HD 56537
This A3V star is located at a distance of 31 pc and has a suspected debris disk based on its IR excess flux at 24 μm [20, 51]. It has one of the largest 24 μm IR excesses amongst older A-type stars, which has been associated with emission from an extended circumstellar (disk) structure [20, 39]. It has an estimated age of ~550 Myr and is the oldest disk studied in our present sample of stars. Thus, we might expect that the rate of occurrence of any detectable exocomet activity would be very low. However, in Figure 2 we show that absorption variability was (unexpectedly) detected over the velocity range −20 to +15 km s−1 on two of the six nights of observation. This velocity range encompasses that of the radial velocity of the star ( km s−1), thus suggesting that at least some of this absorption profile is of a circumstellar origin. Figure 2 shows four of the five nights’ observations of the CaII K-line in which the nights of 03-01-2013, 03-04-2013 (not shown), and 03-06-2013 show no measurable absorption features. However, on the nights of 03-20-2013 and 03-07-2013 an extended absorption feature centered close to the radial velocity of the star ( = −7.4 km s−1) was present. We note that during one night’s observation of this star by Holweger et al. [52] they failed to detect any circumstellar CaII K-line absorption.
3.1.3. HD 58647
HD 58647 is a very young Herbig Ae/Be star of age < 1 Myr that has an associated dusty debris disk identified by a significant mid-IR excess [7]. However, we also note that Manoj et al. [53] have argued that it is a classical Be star. Changes in the broad photospheric NaI and OI absorption line profiles have been attributed to a time-variable circumstellar contribution [54], and Grady et al. [55] have observed similar circumstellar absorption variability (due to gas accretion) in the UV lines of MgII. Spitzer Space Telescope IR spectra recorded towards HD 58647 have revealed molecular emission from polycyclic aromatic hydrocarbons (PAHs) that are known to often populate the planet-forming disks of gas and dust around young stars [56], and CO () emission has been detected by Brittain et al. [57] at a velocity of +16 km s−1.
Our CaII K-line observations shown in Figure 3 reveal a strong absorption line centered at km s−1. At a distance of 320 pc to HD 58647 this absorption component must certainly be of interstellar origin, since it has a similar absorption strength and velocity to that measured towards the nearby stars of HD 58343 and HD 57139 which are both at a comparable distance based on the NaI observations of Welsh et al. [58]. However, on the night of 03-01-2013 an absorption event of EW = 3.5 mÅ at km s−1 was observed which was not seen on subsequent nights. We attribute this temporal absorption event to that of an FEB origin due to the evaporation of exocomet gas.
3.1.4. HD 108767
This nearby AOIV star located at a distance of 27 pc has only a marginal mid-IR excess [20, 22] and has an estimated age of ~260 Myr. Photometric H- and K-band imaging observations have revealed an extended emission profile consistent with the presence of a disk of exozodiacal dust around the star [59]. Previous (single night) CaII K-line observations of this star have revealed a weak absorption centered at km s−1 with an EW = 3.1 mA [60], whose origin has been associated with the local Gemini interstellar cloud by Redfield and Linsky [13]. Our CaII K-line observations also reveal this interstellar absorption with a very similar EW centered at km s−1 (see Figure 4, night 03-02-2013 data). However, all four night’s observations also show broad absorption events in the velocity range = +70 to +130 km s−1 which we associated with FEB events. Over the course of our observations the velocity of some of this circumstellar gas appears to move towards km s−1 and reveals itself as a time-variable absorption that affects the red wing of the interstellar absorption profile (whose total EW increases to 7.4 mÅ) by the night of 03-07-2013. Such variable absorption behavior is very similar to that observed towards stars like β Pictoris and 49 Ceti, both of which have been shown to harbor exocomets in their debris disks [15].
3.1.5. HD 9672
The circumstellar CaII K-line absorption activity associated with HD 9672 (49 Ceti) has been well described by Montgomery and Welsh [15], in which this nearby 40 Myr old AlV star has exhibited significant nightly FEB activity presumably due to the evaporation of gas from infalling exocomets. It has a well-observed dust disk [61] and is one of the few debris disks that have associated sub-mm CO emission [62]. Due to its many physical similarities to the canonical β Pictoris debris disk system, HD 9672 was chosen for subsequent ultraviolet absorption observations using the HST-STIS instrument over the period 08-11–16-2013. Preliminary results from the ultraviolet observations of the neutral absorption lines of CI and OI have been presented by Roberge et al. [46], in which (based on the derived C/Fe ratio) it appears that 49 Ceti has a volatile-rich gas disk similar in some respects to that of β Pictoris.
In Figure 5 we show our observations of the circumstellar CaII K-line recorded three times during the night of 8-11-2013 and one observation on the following night of 8-12-2013. During the 2.5 hours of observations on the first night we detected an FEB event at km s−1, which disappeared by the following night only to be replaced by a new FEB absorption event at km s−1. The EW of the main CaII K-line at km s−1 remained unchanged (within the measurement errors) throughout all of these observations. We note that the projected velocity of the local interstellar cloud of = 11.0 km s−1 derived by Redfield and Linsky [13] is very similar to that of the radial velocity of the star, and thus the main absorption profile consists of both circumstellar and interstellar contributions. The evolution of the FEB event over the 8-11-2013 to 8-12-2013 time period in relation to the circumstellar UV absorption lines with a similar ionization potential to that of CaII will be discussed in a forthcoming paper by Roberge et al. [63].
3.1.6. HD 80007
This A1 III star, located at a distance of 35 pc, is often referred to as being “Vega-like,” but Redfield et al. [64] have failed to detect a significant mid-IR excess that is normally associated with emission from a surrounding dusty debris disk. Previous observations of this star by Hempel and Schmitt [65] have shown a variation in both the EW of the CaII K-line and its absorption velocity recorded over a six-year time period. In addition, very high spectral resolution observations of the CaII K-line by Redfield et al. [64] revealed EW values varying from 0.5 to 3.1 mÅ recorded over a two-month period. Both of these sets of data are in agreement with our present observations of HD 80007 shown in Figure 6. Spectra recorded on the first two nights show a broad absorption profile with a quasi two-component structure, with one component centered near the radial velocity of the star at km s−1. The latter two night’s observations show only a one-component absorption with overall reduced EW, centered close to the stellar radial velocity.
This type of variable absorption behavior is unlike that of the “classic” FEB variable absorption phenomenon, such as that routinely observed towards β Pictoris and 49 Ceti. This may be related to the advanced age of 260 Myr for this system. We note that the Balmer line observed towards HD 80007 has been shown to be variable [66], which might suggest a stellar origin for the absorption variability we have observed. One possible scenario could be associated with mass loss accompanied by diffusion/accretion of matter with a deficit of metals from the surrounding gas/dust disk. Since the convective zone of an A-type star is extremely thin, this may lead to contamination by circumstellar and/or interstellar matter into the stellar photospheric layers. Clearly further observations of this bright and nearby star are required in order to establish whether the observed absorption variability is due to a stellar or exocomet evaporation phenomenon.
3.1.7. HD 109573
In Figure 7 we show archival spectra of the AO-type binary star HD 109573 (HR 4796A) recorded with the ESO FEROS spectrograph on 3 nights within a two-month period in 2007, together with one additional observation in 2010. This star has long been known to possess a significant IR excess [67] with the emission originating in a ring-like dusty circumstellar debris disk thought to be a protocometary cloud [68]. The star has an estimated age of only ~8 Myr [69] and Koerner et al. [70] have suggested that the observed morphology of the dust emission is consistent with its disk being in a transitional planet-forming stage, between that of massive gaseous protostellar disks and more tenuous debris disks such as the one detected around Vega. More recent high spatial resolution imaging of this system by Schneider et al. [71] has found that the spatial distribution of the disk grains is consistent with the possibility of one or more (as yet) unseen exoplanets, and the colors of the dust grains within the disk are consistent with a collisional evolved population of debris.
Far UV absorption observations of the gas associated with the HD 109573 disk have failed to detect H2, and this nondetection implies that the debris disk is almost depleted of molecular gas [72]. Observations of the (weak) CaII K-line over a period of 24 hours by Welsh et al. [27] in 1997 showed no change in its absorption characteristics, which is consistent with higher resolution observations that revealed the central CaII absorption to consist of two absorption components at = −4 and −14 km s−1, both well displaced from the = +7.1 km s−1 radial velocity of the star [65]. The projected velocity of the local G-cloud derived by Redfield and Linsky [13] is = −14.9 km s−1 and thus the observed CaII K-line (unchanging) absorption is not likely to be of circumstellar origin but likely to be of interstellar origin. However, the two observations recorded on the night of 2007-05-04 clearly show a broad and shallow absorption feature at km s−1, which was not detected on the other two nights’ observations. This is most probably an FEB absorption event, since HD 109573 is young, has an associated debris disk, and is thought to harbor exoplanet-like bodies.
3.1.8. Pictoris
HD 39060 (β Pictoris) is probably the most well-observed A-type star with an associated debris disk, exocomet activity, and at least one exoplanet. It is both young (20 Myr) and nearby ( = 19 pc) and its debris disk is viewed nearly as edge-on. A recent summary of some of the properties of the gas disk in this system can be found in present Introduction and in Kiefer et al. [73]. The vast majority of observations of the circumstellar gas in the system have focused on temporal variations of the CaII K-line [50, 74] in order to probe exocomet activity associated with planetesimal material in the debris disk. However, we note that studies of the disk gas surrounding β Pictoris have also focused on the circumstellar absorption associated with NaI [33], CO [75], and CI [76].
Observations of the spatial structure of the optical emission from the debris disk gas by Brandeker et al. [77] and Nillson et al. [78] have revealed over 80 extended emission lines, with the brightest being FeI, NaI, Ti, and CaII. These data have been used in braking gas models by Fernández et al. [34] in order to explain the continued presence of metallic gas in the disk at fixed circumstellar velocities. Such models work best if the gas disk has an enhanced carbon FeI λ3859.91 Å emission from the β Pictoris disk that prompted us to search the ESO archive for any observations of this spectral feature in absorption. This spectral line, which arises in cold ( < 500 K) and neutral (ionization potential < 7.9 eV) gas, is very rarely seen along interstellar sightlines [79, 80] and thus its presence in the β Pictoris gas disk, as first mentioned by Lagrange et al. [81], is potentially intriguing. We note that this absorption line, when present in the spectra of other A-type stars with debris disks that have been thus far observed by us, is much weaker than that found towards β Pictoris.
In Figure 8 we show CaII K-line and FeI λ3860 Å absorption spectra of β Pictoris simultaneously recorded with the HARPS spectrograph at ESO on three nights over the 6-day period 12-23-2003 to 12-29-2003. All 3 spectra show the strong circumstellar CaII K-line at km s−1, each accompanied by an FEB event of varying absorption equivalent widths with velocities in the = +40 to +60 km s−1 range. We also detect FeI absorption, centered at the circumstellar velocity of km s−1, with an (average) absorption equivalent width of ~6.3 mÅ, corresponding to a FeI column density of N(FeI) = 2.5 × 1011 cm−2. Since β Pictoris lies well within the highly ionized region of the Local Cavity as mapped by Welsh et al. [58], this gas column density must be of circumstellar rather than of interstellar origin. We also note the appearance of what appears to be an FEB event centered at km s−1 in the FeI spectrum recorded on 12-28-2003. This is mirrored by an FEB event at a similar velocity simultaneously recorded in the CaII K-line spectrum of that night. To our knowledge, this is the first detection of an FEB absorption event that accompanies the FeI circumstellar line. FEB events were also recorded in the CaII K-line spectra on each of the other two nights, but these were not accompanied by similar absorption events in the corresponding FeI spectra. The FEB event in the CaII K-line spectrum of 12-28-2003 was of a far larger EW value than the other 2 nights, which may suggest that there is a CaII absorption strength threshold above which FEB events must have in order to be detected in the FeI line.
A recent analysis of over 1000 archival spectra of the circumstellar CaII H and K-lines detected towards β Pictoris by Kiefer et al. [73] has revealed the existence of two families of exocomets with different physical properties. The CaII K-line FEB absorption event recorded on 12-28-2003 and shown in Figure 8 is of type “D” (i.e., deep absorption with a narrow FWHM) as opposed to the type “S” FEB events (i.e., shallow absorption with a broad FWHM). The FeI FEB absorption event of 12-28-2003 covers a broad range of velocities (30 to 55 km s−1) but is of a shallow absorption depth with an EW = 4.0 mÅ (i.e., N(FeI) = 1.3 × 1011 cm−2). Its broad profile suggests that it is composed of several absorption components closely spaced in velocity.
This type of absorption profile could be interpreted as originating in the evaporation of gas from several swarms of exocomets, each with similar velocities. Clearly more observations are required in order to investigate any physical relationship between FeI, CaII, and NaI circumstellar absorption. This will be accomplished in a future publication by Montgomery et al. [82] that will utilize all of the 1000+ spectra of β Pictoris that currently reside in the ESO data archive.
3.2. Stars Showing No Variation in CaII-K Absorption
In Figure 9 we show CaII K-line absorption spectra of six stars (HD 88195, HD 106591, HD 135742, HD 158352, HD 161868, and HD 184006) that showed no significant change in the absorption strength between multiple observations typically spanning several nights. On each plot we show the measured EW value of the CaII K-line and in order to help to determine whether the observed line is of interstellar or circumstellar origin we indicate the position of the projected velocity of the most probable cloud component of local interstellar cloud complex as derived by Redfield and Linsky [13]. These data are summarized in Table 3, the last column of which indicates whether the measured absorption is most probably of interstellar (ISM) or circumstellar (C/S) origin.
For the case of HD 88195, with a distance of 183 pc, both of the components in the observed CaII K-line profile are probably dominated by absorption arising along the line-of-sight interstellar medium. However, we note that the (interstellar) absorption component at = +21.8 km s−1 has a very similar velocity to that of the radial velocity of the star and hence a small part of this absorption could be of circumstellar origin. The stars HD 106591, HD 135742, and HD 161868 all have distances < 60 pc such that their sightline absorption may well be dominated by the local interstellar cloud complex. A comparison of the stellar radial velocities of these three stars, as listed in Table 3, indicates that the measured absorption can indeed be associated with that of the local interstellar clouds. Our observations of HD 106591 are consistent with those of Frisch et al. [83] who detected an (interstellar) absorption component of similar strength at km s−1. The present CaII K-line profile of HD 135742 is centered at a very similar velocity to the observations of Redfield and Linsky [13] who detected UV lines (of interstellar origin) at km s−1. We note that our observations of HD 161868 are consistent with previous (far higher resolution) observations by Crawford et al. [84] in which the CaII K-line profile was composed of two components (at = −30 and −33 km s−1) with a total EW = mÅ.
The CaII K-line profile of HD 158352 shown in Figure 9 shows a two-component absorption structure with the main component centered at km s−1, which is identical to the projected velocity of the local interstellar cloud complex listed in Table 3. However, the blue wing of the absorption profile has a component at km s−1 which is very close to the stellar radial velocity value of −36.1 km s−1. Therefore, we identify this component as being of circumstellar origin. However, we observed no change in the strength of this circumstellar component over six nights of observation, and also no FEB absorption activity was detected.
The star HD 184006 ( Cyg) was observed four times over a two-week period in 2013, with all four spectra showing a two-component absorption structure (see Figure 9). The stronger of the two components has a velocity of km s−1, which is close to the projected velocity of the local interstellar cloud complex as listed in Table 3. However, the second (and weaker) absorption component has a velocity of km s−1, which is close to the radial velocity, = −19.5 km s−1, of the star. Thus, we deem it highly likely that this gas component is circumstellar in nature. Neither of the two components changed their respective EW values over the present period of observations. However previous observations of this star with a similar spectral resolving power by us in 2010 [15] suggested that the = −17.4 km s−1 circumstellar component possessed a higher EW value than that recorded in 2013. In fact, the 2010 data showed that the circumstellar component was stronger than the interstellar component over the period 8-19-2010 to 8-23-2010. Those observations now fully support our contention that HD 184006 does indeed possess a circumstellar gas disk that varies in density, albeit slowly, over time.
Finally we note that although we have found no temporal variation in the absorption towards these stars (which would mostly appear to be of interstellar origin) because of the limited number and time coverage of these data we cannot totally rule out the presence of any circumstellar gas and associated FEB activity towards each of these stars.
3.3. Stars with No Measurable Circumstellar CaII K-Line Absorption
Six of the targets listed in Tables 1 and 2 (HD 31295, HD 38678, HD 71155, HD 95418, HD 112429, and HD 125161) had no measurable circumstellar CaII K-line absorption, all typically to a level of < 1.2 mA. All six stars have distances < 38 pc, thus placing them all within the rarefied Local Cavity, which is known to be deficient in neutral gas [58]. We assume that these nondetections are due to either a low level of circumstellar gas surrounding these stars or the fact that the line-of-sight inclination angle of the debris disk is large.
4. Discussion
Our present research has revealed one new late B-type (HD 58647) and four new A-type stellar systems (HD 64145, HD 56537, HD 108767, and HD 109573) that exhibit short-term (night to-night) variation in their circumstellar CaII K-line absorption profiles (i.e., FEB activity). We associate such variability with the evaporation of gas from comet-like bodies on their grazing approach to a parent star. This evaporation of gas, at least for the case of β Pictoris, takes place within of the star. Although the CaII K-line profiles for HD 80007 presented in Figure 6 show circumstellar absorption variability over a six-night time frame, we are unable at present to confidently associate this absorption change with FEB activity. Additionally we have detected circumstellar CaII K-line absorption towards five other stars (HD 88195, HD 106591, HD 135742, HD 158352, and HD 161868) that showed no significant change in absorption strength between several observations. As mentioned previously, although the present (short-term) observations of HD 184006 showed no circumstellar absorption variability, a comparison with previous observations made in 2010 does indicate variability over the long term. All of our stellar targets were selected for observation on the basis of their relative youth, rapid stellar rotation, and mid-IR excess such that they possessed a known (or suspected) associated debris disk. With such similar physical properties it is therefore instructive to determine why some A-type systems with dust disks exhibit variable circumstellar gas absorption and others do not. Our previous work on this subject indicated that FEB activity tends to occur in significantly younger stellar systems that also exhibit chemical peculiarities [3]. We now update these findings using the additional data in the present paper.
In the upper half of Table 4 we list 13 A-type stars (with relevant references) that have observed circumstellar CaII K-line absorption but do not exhibit any significant short-term absorption variability when measured over several nights’ observations. Since circumstellar CaII K-line absorption can often be confused with (local) interstellar absorption with a similar velocity, we deem an absorption component to be of circumstellar origin if it lies within ±6 km s−1 of the radial velocity of the central star (i.e., ~one spectral resolution element) and its velocity is inconsistent with that of the projected velocity of the local interstellar cloud complex as derived by Redfield and Linsky [13]. In the lower part of Table 4 we list (16) stars that do possess short-term (night-to-night) circumstellar absorption variability of the CaII K-line profile that has been associated with FEB activity. Note that our list of stars does not include the A0-type star HD 32297 which has displayed variable circumstellar NaI absorption but no CaII K-line variability [85] nor does it include the AOV-type star HD 158643 (51 Oph) in which Hempel and Schmitt [65] suspect circumstellar CaII K-line variability but has yet to show FEB activity.
In addition to the spectral type, age, and mid-IR excess data in Table 4, column 6 lists whether the star was identified as chemically peculiar (CP) in the work of Zorec and Royer [86]. In column 7 we identify stars labeled as λ Bootis [87, 88], and in column 8 we list the equivalent width value of the λ4481 Å MgII absorption line as measured by Abt and Morrell [89]. These stars are characterized by an underabundance of elements such as Mg, Ca, and Fe, such that a weakness of the MgII line is generally adopted as a primary criterion for the classification of λ Bootis stars.
Firstly all of the 29 targets listed in Table 4 possess a circumstellar gas disk, as indicated by the presence of the CaII K-line near the radial velocity of the central star. We note that 12 of the 29 stars in Table 4 have been classified as shell-type [90], which by definition implies the presence of circumstellar gas surrounding the central star. The majority of those stars with known ages of <100 Myr also possess debris (or transitional) disks, as determined either by their observed mid-IR excess or by direct high-resolution imaging studies in the IR. Hence, just based on these two factors alone we might expect that such systems would in addition to a gas disk that is presumably located close to the central star also possess primordial objects such as planetesimals, comets, and asteroids that are associated with the more distant debris disk.
A statistical analysis of the (unfortunately incomplete) data listed in Table 4 reveals the following results:(1)Although the FEB-hosting stars are on average ~70 Myr younger than the non-FEB-hosting stars, there are only 78% chances that this difference is actually real. Also deriving accurate ages for A-type stars is problematic, with several of the ages listed in Table 4 being obtained from cluster membership rather than by individual stellar observations.(2)There is no statistical difference between the two sets of stars based on either their mid-IR excess or their rotational velocities.(3)As opposed to our previous finding there is no highly statistical difference between the two sets of stars based on either their metallicity or their chemical peculiarity (even if we include λ Bootis stars).(4)One interesting difference between the two sets of stars seems to be that FEB-hosting stars are of an earlier type (on average A1.8 as compared with A3.1). A -test reveals that there is an 88% chance that the FEB-hosting stars are of an earlier type. This could be due to stellar activity levels (such as chromospheric activity or nonradial pulsations) being higher within the atmospheres of hotter A-type stars.
Finally we note that the detection of FEB events is critically dependent on two observational factors that are independent of any of the stellar parameters listed in Table 4. Firstly the viewing angle of the gas disk must be close to being edge-on for an observer to detect a strong absorption signal from the circumstellar gas and also any changes in the CaII K-line absorption profile due to FEB activity. Secondly, our observations generally span timeframes of only ~ one week, and thus A-type debris disk systems that have lower levels of FEB activity would be incorrectly categorized in our Table 4.
5. Conclusion
We have presented medium resolution () spectral observations of the CaII K-line (3933 Å) recorded in absorption towards 15 A-type stars (plus two late B-type stars) with known (or suspected) surrounding debris disk. In addition, we have presented similar spectral absorption data for three A-type stars that were retrieved from the European Southern Observatory on-line data archive. Nightly changes in the CaII K-line absorption profile were observed for five stars (HD56537, HD 58647, HD 64145, HD 108767, and HD 109573). This type of absorption variability that is observed in A-type stars with associated debris disks is thought to originate in a circumstellar gas disk surrounding such young stellar systems. The weak absorption features that sporadically appear within ±100 km s−1 of the central circumstellar CaII K-line absorption profile are most probably associated with the presence of “falling evaporating bodies” (FEBs or exocomets) that liberate evaporative gas on their grazing trajectory towards and around the central star. It is widely believed that planets form through a long period of impacts, mergers, and the accretion of kilometer-sized planetesimals, which themselves have been formed from the accretion of dust grains formed from heavy elements condensing out gaseous circumstellar disks [9, 91]. Our present detection of gas liberated from such bodies that are thought to be the building blocks of planetary systems fully supports current theoretical models of (exo)planetary formation. We note that the detection of such short-term exocomet activity around these young (<100 MYr) stars may currently be the best indication of the presence of associated (but as yet unobserved) exoplanets orbiting these young stellar systems. Hence, we urge further observations of the 16 stars exhibiting FEB activity listed in Table using high precision imaging techniques such as that used by Lagrange et al. [17] in the successful detection of the β Pictoris b exoplanet.
We have also presented archival spectra of the circumstellar FeI λ3860 Å and CaII K-lines observed towards β Pictoris (HD 39060) that, for the first time, reveal the absorption signatures of FEB activity at similar velocities in both of the line profiles. Although we have only made a cursory inspection of the 1000+ spectra residing in the European Southern Observatory (ESO) data archive, our preliminary finding is that an FEB absorption event is only detected in the FeI line whenever a similar FEB absorption event is observed with an equivalent width value of >50 mÅ in the CaII K-line profile. Recent findings, based on observations of the circumstellar CaII K-line, suggest that there are at least two families of exocomets with different physical properties in orbit around β Pictoris [73]. Thus, it would be interesting to investigate any physical relationship between FeI, CaII, and NaI circumstellar absorption in this young disk system. This will be accomplished in a future publication by Montgomery et al. [82] that will use all of the 1000+ spectra of β Pictoris that currently reside in the ESO data archive.
Finally we attempted to determine which, if any, of the stellar physical parameters of young debris disk systems determine whether FEB activity (as defined by the sporadic presence of weak absorption features in the circumstellar CaII K-line profile) is present. None of the parameters, such as rotational velocity, mid-IR excess, or chemical peculiarity, seem to differentiate between stars in which FEB activity is present and those in which it is not observed. Our data suggests that FEB-hosting stars are of an earlier type (on average spectral type Al.8 as compared with A3.l) than those in which no FEB activity has been observed. This could be due to stellar activity levels (such as chromospheric activity or nonradial pulsations) being of a higher level within the atmospheres of these hotter A-type stars. Finally it would appear that stellar age does not seem to be a significant factor in determining the presence of exocomet activity around A-type stars. This is in conflict with most predictions of planetary system formation and may well be due to the inherent uncertainty in the derivation of the ages we have listed for the stars in Table 4.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The authors particularly acknowledge the assistance of the dedicated staff at the McDonald Observatory, the Pie du Midi Observatory, and the South African Astronomical Observatory. Both Ben Bukoski and Sarah Strausbaugh helped with the data extraction and the authors appreciate the work carried out by Rosine Lallement in accessing the NARVAL data for them. Finally they acknowledge the unerring support of the amazing C2.