Over the duration of the Kepler mission KIC 8462852 was observed to undergo irregularly shaped aperiodic dips in flux of up to 20 per cent. The dipping activity can last for between 5 and 80 d. We characterize the object with highresolution spectroscopy spectral energy distribution fitting radial velocity measurements highresolution imaging and Fourier analyses of the Kepler light curve. We determine that KIC 8462852 is a typical mainsequence F3 V star that exhibits no significant IR excess and has no very close interacting companions. In this paper we describe various scenarios to explain the dipping events observed in the Kepler light curve. We confirm that the dipping signals in the data are not caused by any instrumental or data processing artefact and thus are astrophysical in origin. We construct scenarioindependent constraints on the size and location of a body in the system that are needed to reproduce the observations. We deliberate over several assorted stellar and circumstellar astrophysical scenarios most of which have problems explaining the data in hand. By considering the observational constraints on dust clumps in orbit around a normal mainsequence star we conclude that the scenario most consistent with the data in hand is the passage of a family of exocomet or planetesimal fragments all of which are associated with a single previous breakup event possibly caused by tidal disruption or thermal processing. The minimum total mass associated with these fragments likely exceeds 10SUP6SUP MSUBSUB corresponding to an original rocky body of gt100 km in diameter. We discuss the necessity of future observations to help interpret the system.
['comets: general', 'planets and satellites: dynamical evolution and stability', 'stars: activity', 'stars: individual: KIC 8462852', 'stars: peculiar', 'Astrophysics - Solar and Stellar Astrophysics', 'Astrophysics - Earth and Planetary Astrophysics']
Planet Hunters IX. KIC 8462852  wheres the flux
{"KIC8462852 - Where's the flux? †": "- T. S. Boyajian 1 , D. M. LaCourse 2 , S. A. Rappaport 3 ,  \nD. Fabrycky 4 , D. A. Fischer 1 , D. Gandolfi 5 , 6 , G. M. Kennedy 7 , H. Korhonen 8 , 9 , M. C. Liu 10 , A. Moor 11 , K. Olah 11 , K. Vida 11 , M. C. Wyatt 7 , W. M. J. Best 10 , J. Brewer 1 , F. Ciesla 12 , B. Cs'ak 13 , H. J. Deeg 14 , 15 , T. J. Dupuy 16 , G. Handler 17 , K. Heng 18 , S. B. Howell 19 , S. T. Ishikawa 20 , J. Kov'acs 13 , T. Kozakis 21 , L. Kriskovics 11 , J. Lehtinen 22 , C. Lintott 23 , S. Lynn 24 , D. Nespral 14 , 15 , S. Nikbakhsh 22 , 25 , K. Schawinski 26 , J. R. Schmitt 1 , A. M. Smith 27 , Gy. Szabo 11 , 13 , 28 , R. Szabo 11 , J. Viuho 22 , J. Wang 1 , 29 , A. Weiksnar 20 , M. Bosch 2 , J. L. Connors 2 , S. Goodman 2 , G. Green 2 , A. J. Hoekstra 2 , T. Jebson 2 , K. J. Jek 2 , 2 2 2", 'M. R. Omohundro , H. M. Schwengeler , A. Szewczyk': "1 Department of Astronomy, Yale University, New Haven, CT 06511, USA  \n- 2 Amateur Astronomer\n- 3 Department of Physics, and Kavli Institute for Astrophysics and Space Research, Massachusetts Institute of Technology, Cambridge, MA 02139, USA\n- 4 Department of Astronomy and Astrophysics, University of Chicago, 5640 South Ellis Avenue, Chicago, IL 60637, USA\n- 5 Dipartimento di Fisica, Universit'a di Torino, via P. Giuria 1, I-10125, Torino, Italy\n- 6 Landessternwarte Konigstuhl, Zentrum fur Astronomie der Universitat Heidelberg, Konigstuhl 12, D-69117 Heidelberg, Germany\n- 7 Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge CB3 0HA, UK\n- 8 Finnish Centre for Astronomy with ESO (FINCA), University of Turku, Vaisalantie 20, FI-21500 Piikkio, Finland\n- 9 Centre for Star and Planet Formation, Niels Bohr Institute, University of Copenhagen, Øster Voldgade 5-7, DK-1350, København K, Denmark\n- 10 Institute for Astronomy, University of Hawaii, 2680 Woodlawn Drive, Honolulu HI 96822, USA\n- 11 Konkoly Observatory, Research Centre of Astronomy and Earth Sciences, Hungarian Academy of Sciences, H-1121 Budapest, Konkoly Th. M. 'ut 15 -17, Hungary\n- 12\n- Department of the Geophysical Sciences, The University of Chicago, 5734 South Ellis Avenue, Chicago, IL 60637\n- 13 ELTE Gothard Astrophysical Observatory, H-9704 Szombathely, Szent Imre herceg ut 112, Hungary\n- 14 Instituto de Astrof'ısica de Canarias, C. V'ıa L'actea S/N, E-38205 La Laguna, Tenerife, Spain\n- 15 Departamento de Astrof'ısica, Universidad de La Laguna, E-38200 La Laguna, Tenerife, Spain\n- 16 The University of Texas at Austin, Department of Astronomy, 2515 Speedway C1400, Austin, TX 78712, USA\n- 17 Copernicus Astronomical Center, Bartycka 18, 00-716 Warsaw, Poland\n- 18 University of Bern, Center for Space and Habitability, Sidlerstrasse 5, CH-3012, Bern, Switzerland\n- 19 NASA Ames Research Center, Moffett Field, CA 94035, USA\n- 20 Adler Planetarium, Department of Citizen Science, 1300 S Lake Shore Dr, Chicago, IL 60605\n- 21 Carl Sagan Institute, Cornell University, Ithaca, NY 14853, USA\n- 22 Department of Physics, PO Box 64, 00014 University of Helsinki, Finland\n- 23 Department of Physics, University of Oxford, Denys Wilkinson Building, Keble Road, Oxford, OX1 3RH, UK\n- 24 CartoDB, 247 Centre Street, New York, NY 10013, USA\n- 25 Finnish Meteorological Institute, Post Office Box 503, FI-00101 Helsinki, Finland\n- 26 Institute for Astronomy, Department of Physics, ETH Zurich, Wolfgang-Pauli-Strasse 27, CH-8093 Zurich, Switzerland\n- 27 GitHub Inc, 88 Colin P Kelly Jr St, San Francisco, CA 94107 USA\n- 28 Gothard-Lendulet Research Team, H-9704 Szombathely, Szent Imre herceg 'ut 112, Hungary\n- 29 California Institute of Technology, Pasadena, CA 91109, USA  \n26 January 2016", 'ABSTRACT': 'Over the duration of the Kepler mission, KIC 8462852 was observed to undergo irregularly shaped, aperiodic dips in flux of up to ∼ 20 %. The dipping activity can last for between 5 and 80 days. We characterize the object with high-resolution spectroscopy, spectral energy distribution fitting, radial velocity measurements, high-resolution imaging, and Fourier analyses of the Kepler light curve. We determine that KIC 8462852 is a typical main-sequence F3 V star that exhibits no significant IR excess, and has no very close interacting companions. In this paper, we describe various scenarios to explain the dipping events observed in the Kepler light curve. We confirm that the dipping signals in the data are not caused by any instrumental or data processing artifact, and thus are astrophysical in origin. We construct scenario-independent constraints on the size and location of a body in the system that is needed to reproduce the observations. We deliberate over several assorted stellar and circumstellar astrophysical scenarios, most of which have problems explaining the data in hand. By considering the observational constraints on dust clumps in orbit around a normal mainsequence star, we conclude that the scenario most consistent with the data in hand is the passage of a family of exocomet or planetesimal fragments, all of which are associated with a single previous break-up event, possibly caused by tidal disruption or thermal processing. The minimum total mass associated with these fragments likely exceeds 10 -6 M ⊕ , corresponding to an original rocky body of > 100 km in diameter. We discuss the necessity of future observations to help interpret the system.  \nKey words: stars: individual (KIC 8462852), stars: peculiar, stars: activity, comets: general, planets and satellites: dynamical evolution and stability', '1 INTRODUCTION': "For over four years, NASA's Kepler mission measured the brightness of objects within a ∼ 100 square-degree patch of sky in the direction of the constellations Cygnus and Lyrae. The program's targets were primarily selected to address the Kepler mission goals of discovering Earth-like planets orbiting other stars. Kepler targeted over > 150 , 000 stars, primarily with a 30-minute observing cadence, leading to over 2.5-billion data points per year ( > 10 billion data points over the nominal mission lifetime).  \nThe Kepler mission's data processing and identification of transiting planet candidates was done in an automated manner through sophisticated computer algorithms (e.g., Jenkins et al. 2010). Complementary to this analysis, the Zooniverse citizen science network provided the means to crowd source the review of light curves with the Planet Hunters project 1 (e.g., Fischer et al. 2012). In this framework, Planet Hunter volunteers view 30 day segments of light curves in the 'Classify' web interface. A volunteer's main task is to identify signals of transiting planets by harnessing the human eye's unique ability for pattern recognition. This process has shown to have a detection efficiency to identify planetary transits > 85 % using the first Quarter of Kepler data (Schwamb et al. 2012). The Planet Hunters project has now discovered almost a hundred exoplanet candidates, including several con-  \nfirmed systems (Fischer et al. 2012; Lintott et al. 2013; Schwamb et al. 2013; Wang et al. 2013; Schmitt et al. 2014).  \nBecause Planet Hunter volunteers look at every light curve by eye, serendipitous discoveries are inevitable, especially in rich data sets such as that which Kepler has provided. As such, a key aspect of the Planet Hunters project is the 'Talk' interface. 'Talk' is a community discussion board/site where volunteers can discuss light curves and present further analysis on objects viewed in the main 'Classify' interface. In a handful of cases, such as the discovery of the unusual cataclysmic variable, KIC 9406652 (Gies et al. 2013), the default aperture mask used to generate the Kepler light curve was not perfectly centered on the object of interest. Because of this, interesting events in the Kepler light curve would appear to come and go as a result of the shifting orientation of the aperture mask when the spacecraft underwent a quarterly rotation. Events such as these are tagged and discussed on 'Talk', making it possible to return to the raw data target pixel files (TPF) to extract improved light curves with modified aperture masks, for example.  \nThis paper presents the discovery of a highly unusual dipping source, KIC 8462852, from the Planet Hunters project. In just the first quarter of Kepler data, Planet Hunter volunteers identified KIC 8462852's light curve as a 'bizarre', 'interesting', 'giant transit' (Q1 event depth was 0.5% with a duration of 4 days). As new Kepler data were released in subsequent quarters, discussions continued on 'Talk' about KIC 8462852's light curve peculiarities, particularly ramping up pace in the final observations of the Kepler mission.  \nIn this work we examine the full 4 years of Kepler observations of KIC 8462852 as well as supplemental data provided by additional ground- and space-based observations. In Section 2, we characterize KIC 8462852 using Kepler photometry, spectroscopic analysis, AO imaging, and spectral energy distribution analysis. We discover a wide M-dwarf companion to the system and argue that with the data sets we have in-hand, we can exclude the presence of an additional massive gravitationally bound companion nearby.  \nIn Section 3, we visit possible explanations for the peculiar observations of KIC 8462852, including instrumental artifacts, intrinsic/extrinsic variability, and a variety of scenarios invoking lightblocking events. We formulate a variety of scenario independent constraints in Section 4, and elaborate on specific occultation scenarios in Section 5. In Section 6, we conclude by discussing future observations needed to constrain the nature of the object.", '2 DATA': 'KIC8462852, also known as TYC 3162-665-1 and 2MASS J20061546+4427248, is a V ∼ 12 mag star in the Kepler field of view. Its light curve was identified serendipitously by the Planet Hunters project, and was deemed an interesting object that was worthy of further investigation. In the following sections, we characterize the system with data from Kepler as well as additional data from various targeted and archived programs.', '2.1 Kepler photometry': "The Kepler mission was launched on 2009 March 7, and it started science observations on 2009 May 13. The nominal mission was finished almost 4 years later, on 2013 May 12, after the failure of the second reaction wheel. KIC 8462852 was observed throughout the main Kepler mission (divided into Quarters 0 - 17) under longcadence (30-minute) observations yielding an ultra-precise, nearly uninterrupted, light curve during this time. Kepler data files provide both the 'uncorrected' Simple Aperture Photometry (SAP) and the 'corrected' Pre-search Data Conditioning (PDCSAP) fluxes (for details, see Christiansen et al. 2012). In this work, our analysis uses the normalized, PDCSAP data. Note that we have thoroughly validated the data to ensure that any flux variations represent physical events in or near the star (and they do); these processes are described in detail within Section 3.1, and we do not repeat them here.  \nIn Figure 1, we present a montage of plots capturing much of the interesting flux variations observed in the Kepler timeseries data. The top two panels, '(a)' and '(b)', show the flux time series for the entire Kepler mission, but with different vertical flux scales. These show that the flux is relatively constant for most of that time, but is punctuated by a number of substantial dips in flux, including a 15% drop near day 800, and a whole sequence of dips (with one reaching a depth of 22%) after day 1500. Panel '(b)' marks the occurrence of 10 discrete dips (see Table 1). For convenience, we hereafter refer to the two main dip structures between day 788 and 795 and between day 1510 and 1570, as events 'D800' and 'D1500', respectively. Panel '(c)' is a zoom in on the dip D800. The remaining three panels are progressively zoomed in around the exotic complex of dips at D1500.  \nThe D800 dip feature is clean, sharp, and asymmetric in shape. It possesses a gradual dimming lasting almost a week, and transitions back to its nominal brightness in just a couple of days. The D1500 complex consists of many dips, with variable shape and duration, often occurring concurrently as if several independent occultation events were superimposed upon each other. The D1500 dips persist for ∼ 100 days until the Kepler mission's end, and only for a small part of this time does it appear 'quiescent'. There are also other smaller ∼ 0 . 5 % dips, including three earlier in the mission around day 140, day 260, and day 359, and another after the D800 event, around day 1205 (dips #1, 2, 3 and 6, respectively; Figure 1 '(b)', Table 1). Several more 0 . 5 -1 % dips appear in  \nFigure 2. Fourier transform for KIC 8462852. The peaks are labeled with the harmonic numbers starting with 1 for the base frequency. Refer to Section 2.1 for details.  \n<!-- image -->  \nand around the two deep D1500 features, including a ∼ 3 % dip around day 1540. Two small dips occurring at day 1205 and day 1540 have shapes with a similar distinctive, 'triple-dip', symmetric profile, however, they differ in duration by a factor of 3 and in degree of dimming by a factor of 5. All of the fluctuations in intensity visible on these plots are real, i.e., not due to statistical or instrumental variations (Section 3.1).  \nThere are also modulations in the raw flux data at the ∼ 500 ppm level which are visible by eye. To further explore whether any of these modulations are periodic, or have a periodic component, we generated a Fourier transform (FT) of the data with the dips excised from the data train. Figure 2 shows the FT of the Kepler photometry and one can see a clear periodicity of 0.88 day (1.14 cycles/day) and its next two higher harmonics.  \nThis 0.88-day signal is a broad feature that resembles typical FTs of Kepler targets for early type stars (Balona 2013, see their figure 6). If this is a rotation period, then the projected rotational velocity (from Section 2.2) of 84 ± 4 km s -1 represents a minimum stellar radius of ∼ 1 . 46 R glyph[circledot] , consistent with the radius of an F-type star (also see Section 2.2). Also seen in Figure 2 just to the left of the base frequency is a broad collection of smaller peaks. This suggests that something more complicated than a single rotating surface inhomogeneity is producing the observed signal.  \nWe investigate the stability of the frequencies observed in the FT by performing a Short-Term Fourier Transform (STFT), again clipping the data in the dipping regions. In the STFT method, the data are broken up into 'short' segments of 43 d. This segment duration has been selected to optimize both time and frequency resolution. The FT is computed and displayed vertically on the plot, and this is repeated as a function of time, with overlap in time segments to gain back some temporal resolution.  \nThe STFT is presented in Figure 3. This shows that the 0.88 day signal is present in most of the Kepler time series, with the strongest presence occurring around day 1200. Interestingly however, around day 400 and day 1400, we see major contributions at different frequencies, corresponding to ∼ 0 . 96 days and ∼ 0 . 90 days, respectively. We conclude that these are the source of the broad collection of peaks to the left of the base frequency noted above. These low-frequency side-bands could possibly be due to  \nFigure 1. Montage of flux time series for KIC 8462852 showing different portions of the 4-year Kepler observations with different vertical scalings. The top two panels show the entire Kepler observation time interval. The starting time of each Kepler quarter Q is marked and labeled in red in the top panel '(a)'. Dip numbers corresponding to the 10 discrete dips listed in Table 1 are marked and labeled in blue in panel '(b)'. Panel '(c)' is a blowup of the dip # 5 near day 793 (D800). The remaining three panels, '(d)', '(e)', and '(f)', explore the dips (labeled in blue) which occur during the 90-day interval from day 1490 to day 1580 (D1500). Refer to Section 2.1 for details.  \n<!-- image -->  \n1.01  \nDip #  \n1  \n2  \n3  \n4  \n5  \n6  \n7 8 9 10  \nb)  \nNormalized flux  \nNormalized flux  \nNormalized flux  \n1.00  \n0.99  \n0.98  \n0.97  \n0.96  \n0  \n500  \n1000  \n1500  \nc)  \n786  \n788  \n790  \n792  \n794  \n796  \n798  \n800  \n7  \n8  \n9  \n10  \ne)  \n1500  \n1520  \n1540  \n1560  \n1580  \nTime (BJD-2454833)  \nd)  \n1500  \n1520  \n1540  \n1560  \n1580  \n7  \n8  \n9  \n10  \nf)  \n1500  \n1520  \n1540  \n1560  \n1580  \nTime (BJD-2454833)  \nTime (BJD-2454833)  \nTime (BJD-2454833)  \nDip 5 (D800)  \nTime (BJD-2454833)  \n1.00  \nNormalized flux  \nNormalized flux  \n0.95  \n0.90  \n0.85  \n0.80  \n1.010  \n1.005  \n1.000  \n0.995  \n0.990  \n0.985  \n0.980  \n1.00  \n0.95  \n0.90  \n0.85  \n0.80  \n1.01  \n1.00  \n0.99  \n0.98  \n0.97  \n0.96  \n0.95  \nDips 7-10 (D1500)  \nFigure 3. The STFT for the Kepler flux time series. The main base period of ∼ 0 . 88 days is present throughout the span of observations. We identify (at least) two additional frequencies appearing around day 400 and 1400, corresponding to periods of 0.96 to 0.90 days, which we attribute to differential rotation. Refer to Section 2.1 for details.  \n<!-- image -->  \nregions contrasted in flux (e.g., starspots, chemically peculiar regions) appearing at higher latitudes coupled with differential rotation. This is consistent with the differential rotation (or inferred fractional frequency difference of ∼ 10 %) for F-type stars (Reinhold et al. 2013). We would like to note however, that we cannot completely discount the possibility that these periods are due to pulsations. The position of KIC 8462852 is within the Gamma Doradus ( γ Dor) region of the instability strip, where pulsations are observed at < 5 cycles d -1 (e.g., Uytterhoeven et al. 2011). To investigate this, we then compared the STFTs of known γ Dor pulsators to the STFT of KIC 8462852. We found that the dominant frequencies in STFTs for known γ Dor stars do not evolve with time, contrary to the STFT for KIC 8462852. This supports the interpretation that the ∼ 0 . 88 d signal is due to the star's rotational period.  \nWe also report on the presence of variability on the timescale of 10 - 20 days (Figure 2), which, when present, is visible by eye in the light curve 2 . We illustrate this in Figure 4, showing zoomed in regions of the Kepler light curve. The star's 0.88 d period is also evident in each panel as the higher-frequency flux variations. The panel second from the bottom '(c)' shows no low-frequency (10 - 20 day) variations, but the rest do. While the largest of the dipping structures within the D1500 events could also be described as having a periodic structure close to 20 days, the magnitude of the variability and the temporal behavior are much different than these low-amplitude variations described here. Thus, we cannot suggest any connection between the D1500 features and the 10-20 day variability. Finally, we note that the 10 - 20 day variability may actually arise on a faint neighboring star (see Sect. 2.3).  \nThere is another possible periodicity that is worth discussing briefly. In Table 1, we summarize the times and depths of 10 discrete dips present in the Kepler light curve, also labeled in panel '(b)' - '(e)' of Figure 1. If we examine the two most prominent dips (D1568 and D1520; also see panel '(d)' in Figure 1), we see that they have a separation of ∼ 48 . 8 days. We can also see that the D800 dip (dip #5 in Table 1) is separated from the D1520 dip by 15 of these intervals, if the interval is more precisely defined to be 48.4 days. Furthermore, the very shallow dips early in the Kepler time series at D260 and D360 are very close to 26 and 24 of these 48.4-day cycles from the D1520 dip. The other five identified discrete dips (four of which are very shallow), also listed in Table 1, are about a half cycle out of phase with this period to within  \nFigure 4. Stacked plots showing a zoomed-in portion of the Kepler light curve. The star's rotation period of 0.88 d is seen in each panel as the highfrequency modulation in flux. With the exception of panel 'c)', a longer term (10 -20 day) brightness variation is observed, also present in the FT shown in Figure 2. Refer to Section 2.1 for details.  \n<!-- image -->  \nTable 1. Principal Dip Times of KIC 8462852 vs. 48.4-day Period  \n∼ ± 5 %of a cycle. In this exercise, we have neglected the fact that the three most prominent dips in the D1500 region are quite highly structured, and they also have additional minima whose times could have been tagged and included in the analysis. At this time we do not ascribe any particular significance to this period, but it is something to bear in mind as more data on this object become available.", '2.2 Spectroscopy': "We obtained four high resolution ( R = 47000 ) spectra of KIC8462852 with the FIES spectrograph (Frandsen & Lindberg 1999; Telting et al. 2014) mounted at the 2.56-m Nordic Optical Telescope (NOT) of Roque de los Muchachos Observatory in La Palma, Spain. The observations were performed on 11 August 2014, 5 November 2014, 20 November 2015, and 26 November 2015. The data were reduced using standard procedures, which include bias subtraction, flat fielding, order tracing and extraction, and wavelength calibration. The extracted spectra have a S/N ratio of 45-55 per pixel at 5500 ˚ A.  \nFollowing the same spectral analysis procedure described in Rappaport et al. (2015), we use the SPECTRUM code to calculate a grid of synthetic spectra using ATLAS9 models. We then use the co-added FIES spectrum to determine the stellar effective temperature T eff , surface gravity log g , projected rotational velocity v sin i , metal abundance [M/H], and spectral type of KIC 8462852 (Table 3). The plots in Figure 5 show selected regions of the observed spectrum (black) along with the best fit model (red). The temperature we derive ( T eff = 6750 ± 140 K) is consistent with the photometric estimate of T eff = 6584 +178 -279 K from the revised Kepler Input Catalog properties (Huber et al. 2014), as well as with T eff = 6780 K derived from the empirical ( V -K ) colortemperature relation from Boyajian et al. (2013). The projected rotational velocity we measure v sin i = 84 ± 4 km s -1 is also well in line with the one predicted from rotation in Section 2.1, if the 0.88 d signal is in fact the rotation period. Overall, the star's spectrum is unremarkable, as it looks like an ordinary early F-star with no signs of any emission lines or P-Cygni profiles. Finally, we use the stellar properties derived from our spectroscopic analysis to estimate a stellar mass M = 1 . 43 M glyph[circledot] , luminosity L = 4 . 68 L glyph[circledot] , and radius R = 1 . 58 R glyph[circledot] , corresponding to a main-sequence F3 V star based on the empirical calibrations from Pecaut & Mamajek (2013) 3 . Combining the radius (assuming a conservative value of 20% for the radius error), projected rotational velocity, and rotation period (Section 2.1), we determine a stellar rotation axis inclination of 68 ± 29 degrees.  \nWhile interstellar medium features are not typically related to indicators of astrophysically interesting happenings in stars, we note the presence of stellar and interstellar Na D lines in our spectra. In the bottom panel of Figure 5, we show a close up of the region containing the Na D lines ( λλ 5890 , 5896 ˚ A). Within the two broad stellar features, there are two very deep and narrow Na D lines with split line profiles, indicating the presence of two discrete ISM clouds with different velocities of ∼ 20 km s -1 .", '2.3 Imaging': "Figure 6 shows the UKIRT image of KIC 8462852 as well as a similarly bright source ∼ 40 '' away. The PSF of KIC 8462852 is asym-  \nFigure 5. NOT spectrum closeups for KIC 8462852, the best fit stellar model shown in red. Panels show region near H α , H β , Mg, and Na D (top to bottom). The bottom panel shows both the stellar (broad) and interstellar (narrow) counterparts of the Na D lines. Refer to Section 2.2 for details.  \n<!-- image -->  \nby comparison, leading us to speculate that KIC 8462852 has a faint companion star about 1 . 5 -2 '' away.  \nWe observed KIC 8462852 on UT 16 Oct 2014 using the natural guide star adaptive optics (AO) system (Wizinowich et al. 2000) of the 10-meter Keck II Telescope on Mauna Kea, Hawaii. We used the facility IR camera NIRC2 and the J (1.25 µ m), H (1.64 µ m), and K (2.20 µ m) filters from the Mauna Kea Observatories (MKO) filter consortium (Simons & Tokunaga 2002; Tokunaga et al. 2002). We used NIRC2's narrow camera, which produces a  \n0.00994 '' pixel -1 scale and a 10.2 '' field of view. Conditions were cloudy with variable seeing, around 1 '' FWHM. KIC8462852 was observed over an airmass range of 1.26-1.28.  \nThe AO-corrected images have full widths at half maxima (FWHMs) of 39 mas, 43 mas, and 51 mas at JHK , respectively, with RMS variations of about 1-3%. We obtained a series of nine images in each filter. The total on-source integration time was 65 seconds per filter. The images were reduced in a standard fashion using custom scripts written in the Interactive Data Language (IDL). We constructed flat fields from the differences of images of the telescope dome interior with and without lamp illumination. We subtracted an average bias from the images and divided by the flat-field. Then we created a master sky frame from the median average of the bias-subtracted, flat-fielded images and subtracted it from the individual reduced images. The individual reduced images were registered and stacked to form a final mosaic (Figure 7).  \nAs suspected from the asymmetric UKIRT image, the Keck AO images reveal an obvious faint companion at a separation of 1.95 '' and position angle of 96.6 · . To measure the flux ratios and relative positions of the two components, we used an analytic model of the point spread function (PSF) as the sum of two elliptical Gaussian components, a narrow component for the PSF core and a broad component for the PSF halo, as we have done for other binaries (Liu et al. 2008). For the individual images obtained with each filter, we fitted for the flux ratio, separation, and position angle of the binary. To correct for optical distortions in NIRC2, we used the calibration of Yelda et al. (2010). The averages of the results were adopted as the final measurements and the standard deviations as the errors (Table 3).  \nIt is unclear whether this is a physical or visual binary, though given the delta magnitude and separation, the chance alignment of the companion being a background or foreground object is only ∼ 1 %(Rappaport et al. 2014). At ∼ 2 %of the flux of the brighter star, this would be a ∼ 0 . 4 M glyph[circledot] M2V star, if it is indeed at the same distance as our target F star (Kraus & Hillenbrand 2007). The JHK colors are also consistent with the companion being a dwarf, not a giant (Bessell & Brett 1988). If we take the magnitude of KIC8462852 as V = 11 . 705 , and the absolute visual magnitude of an F3V star to be V = 3 . 08 (Pecaut & Mamajek 2013), then the (reddened) distance modulus is 8.625. We derive a de-reddened distance of ∼ 454 pc using E ( B -V ) = 0 . 11 (Section 2.4; corresponding to a V -band extinction of A V = 0 . 341 ). Assuming the fainter star is associated with the main F-star target, the angular separation of ∼ 1 . 95 '' translates to a distance of ∼ 885 AU. At this separation, the second star cannot currently be physically affecting the behavior of the Kepler target star, though could be affecting bodies in orbit around it via long term perturbations (see Kaib et al. 2013). If such a star is unbound from KIC 8462852, but traveling through the system perpendicular to our line of sight, it would take only 400 years to double its separation if traveling at 10 km sec -1 . So, the passage would be relatively short-lived in astronomical terms.  \nWealso obtained Speckle observations of KIC 8462852 on the night of UT 22 Oct 2015 using the DSSI instrument on the WIYN telescope located on Kitt Peak (Howell et al. 2011). Observations were made simultaneously in two filters with central wavelengths at 692 and 880 nm. Both filters show the source to be single, with no visible companion observed to within 0.08 arcsec and brighter than a delta magnitude of 3.8 and 4.2 magnitudes (for the 692 and 880 nm filters, respectively). The companion star seen in the Keck NIRC2 image was not detected, favoring the conclusion that it is an M-dwarf, which would be too faint to be detected in the reddest  \nFigure 6. UKIRT image for KIC 8462852 and another bright star for comparison, showing that it has a distinct protrusion to the left (east). For reference, the grid lines in the image are 10 '' × 10 '' . The color coding is logarithmically scaled. Refer to Section 2.3 for details.  \n<!-- image -->  \nFigure 7. Keck AO H -band image for KIC 8462852 showing the companion was detected with a 2 '' separation and a magnitude difference ∆ H = 3 . 8 . The color coding is logarithmically scaled. Refer to Section 2.3 for details.  \n<!-- image -->  \nDSSI filter (880 nm). However, it is important to note that these speckle results provide an independent confirmation of the results from Keck AO: KIC 8462852 has no additional companion down to a separation of ∼ 20 AU detectable within the relative brightness limits with each instrument.  \nFinally, we speculate that the 10-20 day periodicity discussed in Sect. 2.1 might actually arise on the neighboring faint M star. The amplitude of those modulations are ∼ 500 ppm of the total target flux. If they arise on the M-star, then their fractional modulation of that star would be as high as 3%, which would not be unusual for an M star.", '2.4 Spectral energy distribution': 'The spectral energy distribution (SED) of KIC 8462852 including optical, 2MASS (Skrutskie et al. 2006), (ALL)WISE (Wright et al. 2010), and Galex NUV (Morrissey et al. 2007) flux densities is shown in Figure 8. Optical photometry in BV ( RI ) C filters was obtained by the 90 cm Schmidt telescope of the Konkoly Observatory at Piszk\'estet"o Mountain Station. For standard magnitudes  \nFigure 8. SED for KIC 8462852. The Black solid line is a model for a star with T eff = 6750 K and E ( B -V ) = 0 . 11 . Flux calibrated photometry are plotted in red, where the extent of the \'error-bars\' in the X-direction indicate the wavelength range of each bandpass and the Y-direction shows the error of the flux measurement. Refer to Section 2.4 for details.  \n<!-- image -->  \nGD391ABCE photometric standard stars were used as comparison (Landolt 2013). Photometric magnitudes are listed in Table 3.  \nIn order to study whether the system exhibits excess at midinfrared wavelengths, we first fitted an ATLAS9 atmosphere model (Castelli & Kurucz 2004) to the photometric data points between 0.15 and 3.6 µ m. From the grid of model atmospheres we selected the one that has the closest metallicity, surface gravity, and effective temperature to those derived from our spectroscopic study. Thus we fixed T eff , log g , and [Fe/H] parameters to 6750 K, 4.0, and 0.0, respectively, and only the amplitude of the model and the reddening were fitted. The best fitted photospheric model is displayed in Figure 8. We derive a reddening of 0 . 11 ± 0 . 03 mags. By comparing the measured W 2 and W 3 WISE flux densities at 4.6 and 11.6 µ m(at 22 µ mwehave only an upper limit) with the predicted fluxes derived from the photosphere model we found them to be consistent, i.e., no excess emission can be detected at mid-infrared wavelengths. This lack of significant IR excess is independently confirmed using warm Spitzer/IRAC data by Marengo et al. (2015).  \nHowever, this does not exclude the existence of a colder debris disk or a warmer, but relatively tenuous disk. Assuming that the emitting grains act like a blackbody, we can derive their characteristic temperature at a specific stellar-centric distance. Using this approach, we compute the SED of a narrow dust belt located at a distance of 1, 2, 3, 5, and 10 AU from a star with a luminosity of 4.7 L glyph[circledot] , corresponding to the main-sequence stage (Pecaut & Mamajek 2013). The W 3 and W 4 band photometry were then used as upper limits to set the amplitude of the excess. Figure 8 shows the result of these computations and summarizes the fundamental disk properties (dust temperature, upper limits for fractional luminosity) of the dust belts at different radii. It is worth noting that this very simple model accounts only for large blackbody grains, smaller ( µ m-sized) grains are ineffective emitters and may be heated to higher temperatures compared to larger grains at the same location. We revisit this analysis in more detail later in Section 4 (also see Figure 11).', '2.5 Ground-based photometric surveys': 'We reviewed the ∼ 700 photometric intensities from the years 1900 - 2000 from the Digital Access to a Sky Century Harvard  \n(DASCH) project 4 (Grindlay et al. 2012). The error bars on the photometry are about ∼ 10 %. At this level, we found the star did not do anything spectacular over the past 100 years. However, if it underwent several ∼ 20 %dips in flux lasting for several days each during that period, the chances are high that there were no plates exposed at those times.  \nSuperWASP data (Butters et al. 2010) are unremarkable for KIC8462852. We note that there is a 0.2 magnitude offset between the available SuperWASP data sets. However, we see the same offset when comparing its photometry with a similarly bright source nearby KIC 8462852. Thus, we reject this being real (e.g., due to a flaring event, etc.).  \nUnfortunately, KIC 8462852 falls outside the area covered by the KELT network (T. Beatty, private communication).', '2.6 Limits on a close companion': "We use the four FIES spectra (Section 2.2) to measure the presence of any Doppler shifts induced by a companion. We traced the radial velocity (RV) drift of the instrument by taking longexposed ThAr spectra in a bracketed sequence, i.e., right before and after each target observation. RV measurements were derived by cross-correlating the target spectra with the rotationally broadened best fitting Kurucz model. The RV measurements are listed in Table 2 along with the error bars and the barycentric Julian dates in barycentric dynamical time. To within the ∼ 400 m s -1 uncertainties in the RV measurements, the four values we measure are quite consistent with no change at all over the 470 day observation interval.  \nIn order to quantify what limits we can set on the mass of an hypothetical close companion star, we carried out the following analysis. We assumed a circular orbit because there are insufficient data points to fit for the parameters in an eccentric orbit. Then, for each in a sequence of 4 × 10 5 trial orbital periods, P , in the range of 0.5 to 3000 days, we fit the four RV points with a sine and cosine term to represent the orbit and a systemic γ velocity. From this fit we computed the velocity semi-amplitude K and added its 2σ uncertainty to establish a conservative upper limit to K. We then used the upper limit on K to compute the corresponding upper limit on the mass function. Finally, we solved for the upper limit on the mass of the hypothetical close companion by taking the mass of the F star to be 1.4 M glyph[circledot] , and assuming three different orbital inclination angles ( 30 · , 60 · , and 90 · ). The results are shown in Figure 9. The spikes are at values of P where the epochs of the four RV measurements are commensurate with being at orbital phase 0, and the mass constraints are weaker at these periods. For longer periods, the density of these spikes diminishes greatly and the lower locus of points can be taken as a likely upper limit on the mass of any companion. Therefore, we conclude that for periods between ∼ 30 and 300 days, the mass of any companion is very unlikely to exceed that of a brown dwarf.  \nAnother diagnostic to constrain the nature of the companion uses the FT in Figure 2, which shows no sharp, narrow peaks without harmonics (Section 2.1). With this information, a very basic limit can be set on a companion from the lack of observed ellipsoidal light variations (ELVs). The ELV amplitude A ELV is expressed as:  \nA ELV glyph[similarequal] 1 . 5( M c /M ∗ )( R ∗ /a ) 3 sin 2 i (1)  \nTable 2. FIES RVs of KIC 8462852  \nSee Section 2.6 for details.  \nFigure 9. Upper limits (2σ confidence) to a hypothetical companion mass from the lack of ELVs (red curves) and lack of RV variations on four occasions (blue curves). Each type of constraint is shown for three different assumed orbital inclination angles ( 30 · , 60 · , 90 · ); these are marked directly on the ELV constraint curves, and can be inferred from the dashed, dotted, and solid curves, respectively, for the RV constraints. These results indicate that there are no objects heavier than a super-Jupiter in close orbits with P orb glyph[lessorsimilar] 2 days, and likely no heavier in mass than a brown dwarf for P orb glyph[lessorsimilar] 300 days. Refer to Section 2.6 for details.  \n<!-- image -->  \n(e.g., Kopal 1959; Carter et al. 2011) where M ∗ and R ∗ are the mass and radius of the primary, a and i are the semimajor axis and orbital inclination, and M c is the mass of a putative companion. Rearranging to express a as the orbital period P using Kepler's third law, this equation simplifies to:  \nA ELV glyph[similarequal] 3 . 3 × 10 -5 ( M c / M J )(1d /P ) 2 sin 2 i (2)  \nwhere now the companion mass M c is expressed in Jupiter masses M J and P is in days. If ELVs were present, we would have seen a peak glyph[greaterorsimilar] 50 ppm for all periods shorter than 4 days ( glyph[greaterorsimilar] 0 . 25 cycles day -1 ) in the FT (Figure 2).  \nThe limits on the companion mass that we can set from the lack of ELVs, as a function of orbital period, are illustrated in Figure 9. They are plotted as red lines for three different assumed inclination angles. Note that an angle of 90 · is not allowed or we would have seen (regular) transits; it is shown in this figure for instructive purposes only. These ELV mass constraints are superposed on those discussed above based on the lack of RV differences among our four measurements (blue curves). Taken together, these results indicate that there is not likely to be a close companion to the F star more massive than a super-Jupiter with P orb glyph[lessorsimilar] 2 days, nor more massive than a brown dwarf for P orb glyph[lessorsimilar] 300 days.  \nTable 3. Properties of KIC 8462852  \na Upper limit.", '2.7 Space motion and age': 'Using our distance estimate of 454 pc (Section 2.3), the radial velocity obtained from the FIES spectrum (Section 2.6), and proper motions and positions from the UCAC4 catalogue we computed the Galactic space motion of the target, yielding +31 . 5 , -2 . 5 , and +10 . 2 km s -1 for the U (defined as positive toward the Galactic center), V, and W velocity components, respectively. Young disk population stars have low velocity dispersion and they occupy a special region within the velocity space. Based on the studies of Eggen (1989), Leggett (1992) defined a box by -50 < U < +20 km s -1 , -30 < V < 0 km s -1 , and -25 < W < 10 km s -1 , which includes most of the young disk stars in our neighborhood. Our target lies outside of this box. In fact, its galactic space motion - especially the U component - deviates significantly from the characteristic space motion of any nearby young ( < 100 Myr) kinematic groups, open clusters, and star forming regions (Makarov 2007; Mamajek 2015). Altogether, it implies that KIC 8462852 likely does not belong to the youngest stellar population.  \nIn making this distance estimate, we assumed that KIC8462852 is a main-sequence star (Section 2.3). We note that assuming a pre-main or post-main sequence phase does not change our previous conclusion. These evolutionary stages would be accompanied by larger luminosities and thereby larger distances. This would result in a galactic space motion that deviates even more significantly from that of typical young disk stars. Unfortunately, our star falls outside the region where empirically calibrated age diagnostics such as chromospheric activity or stellar rotation period can be used (e.g., Mamajek & Hillenbrand 2008).', '2.8 Similar dippers in the Kepler field?': "The anomalous dips in KIC 8462852 were serendipitously found by the Planet Hunter citizen science group. Due to its aperiodic nature, it likely never would have been flagged/recovered by most searches for transits, eclipsing binaries, or asteroseismologically interesting stars. However, knowing the existence of KIC 8462852's light curve, we naturally wondered if there are, in fact, numerous other such objects in the main-field Kepler data base. We therefore applied a simple algorithm to search the data base for other systems similar to KIC 8462852. The algorithm consisted of searching for dips with depths of greater than 10% (i.e., normalized fluxes of < 0 . 9 ) that consist of 5 or more consecutive Kepler long-cadence samples (i.e, lasting more than ∼ 2 . 5 hours). In all, this search turned up more than a thousand targets with this signature. The vast majority of them, however, were due to (1) eclipsing binaries, (2) the rotation signature of large amplitude starspots, and (3) some obvious Kepler data artifacts. We carefully examined the remaining small number of systems by eye, but could identify none that was reminiscent of KIC 8462852. We also lowered the threshold for dips to 5%, but the search likewise turned up no candidates that one would believe closely resemble KIC 8462852. Of course, some of the visual comparison work is necessarily qualitative, but we were satisfied that there are at most a few similar systems to be found in the main Kepler field.", '3 POSSIBLE EXPLANATIONS OF THE OBSERVED DIPPING EVENTS OBSERVED IN KIC8462852': 'The main issue in explaining the peculiar light curve for KIC8462852 is related to the presence of multiple dimming events, that are not periodic and of which the D800 single event has a smooth, yet highly asymmetric, profile, and the D1500 events are the deepest and most complex. Here, we introduce several scenarios to explain KIC 8462852 and discuss how the observational data do and do not support each theory.', '3.1 Instrumental effects or data reduction artifacts?': "The Kepler light curve for KIC 8462852 is unique, and we have thoroughly explored the raw data for defects/instrumental effects, which could cause the observed variations in KIC 8462852's flux. We use the PYKE software tools for Kepler data analysis to check the data for instrumental effects. We check the following possibilities:  \n- · We checked that the same flux variations, i.e., the 'dips', are present in the SAP FLUX data set.\n- · We verified that data gaps and cosmic rays events 5 do not coincide with the dipping events, as they are prone to produce glitches in the corrected fluxes.\n- · Weverified at the pixel-level that there are no signs of peculiar photometric masks used in making the light curves.\n- · We verified at the pixel level that the image light centroid does not shift during the 'dipping' events\n- · We inspected the Kepler light curves of neighboring sources and find that they do not show similar variability patterns in their light curves.  \n- · We determined that CCD cross talk, reflection, and column anomaly cannot be the cause (Coughlin et al. 2014).\n- · We verified with the Kepler team mission scientists that the data were of good quality.  \nThis analysis concludes that instrumental effects or artifacts in the data reduction are not the cause of the observed dipping events, and thus the nature of KIC 8462852's light curve is astrophysical in origin.", '3.2 Intrinsic variability?': "An example of a class of stars which display intrinsic variability are the R Coronae Borealis (RCB) type variables. These are highly evolved F-G supergiants (e.g., Clayton 1996) that have light curves which show pulsations (on the order of months) and irregular deep dips (lasting weeks to months). Their 'dipping' variability is associated with formation of clouds that obscure the photosphere, and is often observed as a sharp decrease in flux followed by a more gradual, and sometimes staggered, recovery. In the case of KIC 8462852 the time scales of the dips are different than those of a RCB variable. Likewise, the ingress at D800 has a gradual decrease in flux, which is inverse to what is expected in a RCB, and the dip shapes at D1500 are also non-characteristic of a RCB. Lastly, the spectroscopic log g and v sin i are far from those of a supergiant. These items together strongly rule out the possibility of KIC 8462852 being a RCB variable.  \nAnother possibility is the self-emission of disk material from the star itself, as in the case of Be-stars. Be stars are rapidly rotating (almost near breakup) stars that are usually of spectral class O and B, but sometimes A, and exhibit irregular episodic outbursts. Usually these outbursts are in emission, but in some cases it can also result in dimming (see Hubert & Floquet 1998). Be stars also often exhibit quasi-periodic oscillations in the range of ∼ 0 . 5 -1 . 5 days. This also fits the bill for what we see in the FT of KIC 8462852 ( § 2.1). It has been hypothesized (e.g., Rappaport & van den Heuvel 1982) that most, if not all, Be stars have a binary companion which originally transferred mass to the current Be star to spin it up to near breakup (the remnant of that star is sometimes found to be a neutron star). The periods of these binaries range from a couple of weeks to thousands of days (perhaps longer). If KIC 8462852 is a Be star, we would get an unprecedented look into the inner disk behavior. In such as case, the broad peak in the FT at frequencies just below the 0.88 d periodicity could be explained by ejected material in a so-called 'excretion disk' that is moving outward but with roughly Keplerian velocity.  \nThe lack of observed IR excess does not support the existence of an excretion disk. There is also an absence of H α emission in the star's spectrum (although, as noted above, Be star H α emission is known to be variable and turn off and on with timescales from days to years). Furthermore, the temperature of KIC 8462852, T eff = 6750 K, is too cool to be a Be star. It is also unlikely to have been spun-up by a donor star whose remnant is still orbiting the F star because of the constraints set by the four RV measurements and the limits on any ELVs (see Section 2.6). Though, we cannot rule out remnants orbiting with P glyph[greaterorsimilar] a few years.", '3.3 Extrinsic variability?': "3.3.1 Related to the secondary star  \nWe first consider whether KIC 8462852's flux is contaminated by the nearby M-dwarf detected with high-resolution images ( § 2.3).  \nWhether or not the system is bound, the faint companion contributes light in the Kepler photometric aperture, which in turn affects the observed signal in the light curve. Our observations show that the flux ratio in the infrared is ∼ 30 , which translates to a factor of several hundred in the Kepler bandpass. Thus, the maximum imprint that the M-dwarf has on the light curve variability is ∼ 30 mmags; this is insufficient to make an impression on KIC8462852's light curve at anything greater than ∼ 3%, and, in particular, it couldn't possibly explain any of the large dips.", '3.3.2 Occultation by circumstellar dust clumps': "The dips could be readily explained in terms of occultation by an inhomogeneous circumstellar dust distribution. However, this does not mean that the dust distribution that would be required to explain the observations is physically plausible or would necessarily apply to KIC 8462852.  \nInhomogeneous dust distributions have been invoked to explain dips seen towards some young stars such as UX Orionis or AA Tau-like 'dipper' systems (Herbst et al. 1994; Herbst & Shevchenko 1999; Morales et al. 2009; Cody et al. 2014; Ansdell et al. 2015). At an age of only a few tens of Myrs, these dipper stars have V -band light curves characterized by sporadic photometric minima with amplitudes of 2 - 3 magnitudes and with durations of days to many weeks. These objects also generally exhibit strong infrared excess, starting at ∼ 2 -5 µ mand show signs of accretion (emission) in their spectra. However, in contrast to such systems, KIC8462852 has no detectable IR excess or accretion signature to suggest that it is a young T Tauri star (Sections 2.2, 2.4). Thus a scenario in which material in a gas-dominated protoplanetary disk occults the star due either to accretion columns or non-axisymmetric azimuthal or vertical structure in the inner disk (e.g. Herbst et al. 1994; Herbst & Shevchenko 1999; Bouvier et al. 1999; McGinnis et al. 2015) is strongly disfavoured.  \nWe therefore are left to consider scenarios that could arise around a main-sequence or weak-line T Tauri star that has dispersed its protoplanetary disk, but still hosts a gas-poor planetary system that may include planets, asteroids, and comets. The 'clumps' of dust passing in front of the star could perhaps lie within an optically thin asteroid belt analogue that is otherwise undetected, or be more isolated objects such as remnants of a broken up comet. As in the above scenarios, the typical minimum sizes of the dust grains are ∼ µ m(e.g. Backman & Paresce 1993), which are able to cause stellar variation by absorbing and scattering starlight at optical wavelengths. Before considering such scenarios in more detail, we start with some scenario-independent constraints that can be gleaned from the observations.", '4 SCENARIO-INDEPENDENT CONSTRAINTS': "To understand what could be the origin of the clumps it would help to know where they are located in the system, how big they are, and how long they last. To aid with this discussion, Figure 10 shows some scenario-independent constraints on the size and orbital distance of the clumps that are discussed further below. The only assumption for now is that the clumps are on circular orbits, but this assumption is relaxed later in Section 5.4. Some of the constraints also assume the clumps to be opaque, but again this assumption is relaxed later.  \nDip duration: The timescale t dip for the transit of a clump of radius s with transverse velocity v t across the equator of a star with  \nFigure 10. Size vs. semi-major axis parameter space for optically thick, spherical dust clumps on circular orbits around a star of M ∗ = 1 . 43 M glyph[circledot] and R ∗ = 1 . 58 R glyph[circledot] . The solid lines represent dips of equal duration (as labelled). Dotted lines show minimum clump sizes for dips of different depths. Vertical dashed lines show where the orbital period is 1500 days, and where the light curve gradient for an optically thick 'knife edge' could be as high as 0.5 d -1 . Diagonal dashed lines show Hill radii of planetesimals of different sizes, assuming a density of 3 g cm -3 . Combined, the period, gradient, and duration constraints in the circular orbit scenario suggests the clumps lie between ∼ 3 -10 AU, and have sizes similar to the star.  \n<!-- image -->  \nradius R ∗ is t dip = 2( s + R ∗ ) /v t . If the clump is on a circular orbit around a star of mass M ∗ with semi-major axis a , and is much less massive than the star, then  \ns ≈ 1 . 85 t dip ( M ∗ a ) 1 / 2 -R ∗ , (3)  \nfor a is in units of AU, M ∗ in M glyph[circledot] , s and R ∗ in R glyph[circledot] , t dip in days. Thus, the several-day duration of the events for KIC 8462852 suggests that the clumps are either close-in and large compared to the star, or far-away from the star and small. However, clumps that are too distant move too slowly across the stellar disk to explain the observed duration regardless of their size; e.g., a 3-day duration dip cannot arise from a clump beyond ∼ 15 AU.  \nDip depth: A minimum clump size is set by the depth of the dimming events, which we characterise as 1 minus the normalised flux, which we call τ . For example, even if the clump is completely opaque, the maximum dip depth is max( τ ) = ( s/R ∗ ) 2 . The deepest τ = 20 % dimming event at D1500 thus implies that at least some clumps are a sizeable fraction of the stellar size. A dip caused by a fully optically thick symmetrical clump would also have a characteristic symmetrical shape which does not resemble those observed (e.g., panel 'c' in Figure 1), so this can be regarded as a strong lower limit. While there appear to be a range of event durations, the duration of the deepest events is at most about 3 days. The middle solid line in Figure 10 (for t dip = 3 d and a depth of τ = 20 %) therefore decreases the outer limit on the clump locations mentioned above to closer to 8 AU.  \nLight-curve gradient: A similar, but independently derived, outer constraint on the clump location can be set by examining the gradients in the light-curve, which are at most half of the total stellar flux per day (i.e. 0.5 d -1 when the light curve is normalised to 1). Orbiting material can change the light-curve most rapidly when it is optically thick and passing the stellar equator (i.e., the  \nFigure 11. Fractional luminosity limits (blue lines) and an estimate of the system dust content from the light curve (green line). The dust level is constrained to lie below the blue line by the WISE photometry (4.6 µ m, 12 µ m, and 22 µ m). The green line integrates the optical depth in the light curve assuming that clumps are similar in size to the star and on circular orbits. If the clumps lie beyond about 0.2 AU the IR non-detection of the dust is unsurprising, although many scenarios require more emission than that from dust seen to pass along our line-of-sight to the star. Refer to Section 4 for details.  \n<!-- image -->  \n'knife edge' model of van Werkhoven et al. 2014). The high rate of change in the KIC 8462852 light curve translates to a lower limit on the transverse velocity of the orbiting material of about 9 km s -1 , which corresponds to an upper limit of 13 AU for material on circular orbits, although as discussed in Section 5.4, this upper limit is closer to the star if the clump is optically thin.  \nNon-periodicity: The lack of evidence for a clear periodicity in the dips in the observed light-curve excludes orbital periods shorter than ∼ 1500 days, which thus constrains the location to lie beyond about 3 AU. This constraint could be broken if the clumps disperse within a single orbit. Likewise, if the two deep dipping events at D800 and D1500 are from the same orbiting body (or bodies), a period of 700 - 800 days remains a possibility.  \nGravitational binding: To address the survival of the clumps, we note that in any scenario where the clumps are not selfgravitating, they cannot be long-lived in the face of orbital shear (e.g. Kenyon & Bromley 2005) and their internal velocity dispersion (e.g. Jackson & Wyatt 2012). Figure 10 therefore shows planetesimal sizes required to retain dust clouds within their Hill sphere, R Hill = a ( M pl / [3 M ∗ ]) 1 / 3 , as one way of ensuring long-lived clumps.  \nThus, under the assumption of circular orbits, the depth, duration and lack of periodicity of the dimming events constrains their location to a region roughly corresponding to that occupied by the giant planets in the Solar System (i.e., between the green dashed lines). Clump sizes would thus be comparable to, but larger than, the star (i.e., above the uppermost horizontal dotted purple line), and they would have to have high, but not necessarily unity optical depth. It might be possible to explain the clumps as dust bound to planetesimals larger than around 1000 km, which means such planetesimals are not necessarily large enough for direct transit detection (the lack of which could provide another constraint).  \nInfrared excess: Another constraint on the origin of the clumps comes from the lack of infrared emission (Section 2.4).  \nFigure 12. Inverted light curve for KIC 8462852 portraying the blocking factors needed to reproduce the light curve as a function of time. Refer to Section 4 for details.  \n<!-- image -->  \nAssuming the clumps are larger than the star, the Kepler light curve provides blocking factors needed as a function of time, ln(normalized flux) , where ln(normalised flux) ≈ τ for small τ , as shown in Figure 12. This optical depth and the assumption that the clump crosses the star at its orbital velocity allows conversion to optical depth as a function of distance along the clump. The dimming events therefore allow an estimate of the minimum possible cross-sectional area σ tot of dust in orbit around the star. That is,  \nσ tot = v t h ∫ τ ( t ) dt, (4)  \nwhere the light-curve yields ∫ τ ( t ) dt ≈ 0 . 86 days, v t is the velocity of the clumps (assumed to be uniform at circular velocity for a given semi-major axis), and h the 'height' of the clumps (i.e. their size along the dimension perpendicular to their velocity). The height of the clumps is assumed to be 2 R ∗ , though it could be higher if not all of the clump crosses the stellar disk (e.g., this could be assumed to be πs/ 2 for large spherical clumps passing directly across the star). This calculation gives the minimum possible crosssectional area as  \nσ tot = 2 . 6 × 10 -4 a -1 / 2 AU 2 (5)  \nwhere a is in AU, the dependence on which arises from the velocity at which the clump crosses the star.  \nThis cross-sectional area can then be converted to fractional luminosity at a given distance from the star using f = σ tot / (4 πa 2 ) . The blue lines in Figure 11 show the limits on the dust fractional luminosity f = L dust /L ∗ derived from the SED (Section 2.4). These can be thought of as the maximum luminosity of blackbodies at a range of dust temperatures (or stellocentric radii) that fit under the WISE photometry. The dust estimate from Equation 4 is shown as a green line, and the fact that it lies below the blue line at all radii beyond 0.2 AU indicates that it is perhaps not particularly surprising that no mid-IR excess was seen.  \nHowever, this dust area estimate is only a lower limit since it only includes the dust which passed in front of the star during the lifetime of the Kepler mission. The true area would be larger if there are more clumps further along the orbit which have yet to pass in front of the star, and could also be larger if the dips do not capture all of the cross-sectional area in their clumps. Furthermore, for some specific scenarios discussed in the following sections, the  \npresence of clumps that pass in front of the star requires the existence of other clumps that do not pass along our line-of-sight. The lack of infrared emission thus places constraints on how many such clumps there are in the system. For example, Figure 11 shows that for clumps at a few AU the cross-sectional area can only be increased by 3 orders of magnitude before it is detectable by WISE. The calculation is further complicated should the clumps be considered to be short-lived, or on non-circular orbits.  \nMass estimates: The minimum possible cross-sectional area required to cause the observed dips, σ tot (Equation 5), can also be used to determine a minimum possible dust mass, m tot . If the dust all has the same diameter D and density ρ then m tot /σ tot = 2 ρD/ 3 , resulting in a total mass of 6 . 7 × 10 18 g for 1 µ mdiameter dust of density 3 g cm -3 orbiting at 3 AU (and scaling as a -1 / 2 ρD for different assumptions). If all of this mass were put in a single body of the same density ρ , this would have a diameter of 16 km. This illustrates that the minimum mass of the parent body required to cause this phenomenon is approaching the mass of comet HaleBopp. However, this calculation has a few caveats. For one, the value derived for σ tot is an absolute minimum given that it only accounts for the material which passed in front of the star during the observations. It also does not account for the possibility that the dust in the clump has a range of sizes. For example, for dust with a power law size distribution with index of 3.5 (Dohnanyi 1969) extending from D min to D max , the ratio of mass to cross-sectional area scales ∝ ρ √ D min D max . Thus, this estimate would be 100 times larger than that derived above if the size distribution extended from 1 µ mup to 1 cm.  \nGiven these basic constraints we now consider several scenarios that may explain the observations. The first two are related to collisions within an asteroid belt (Section 5.1) or unstable planetary system (Section 5.2), the third considers dust that orbits within the Hill spheres of large planetesimals which may reside in an asteroid belt but are not required to collide (Section 5.3), and the fourth is that the dips are the passage of a series of fragments from a brokenup comet or asteroid on a highly elliptical orbit (Section 5.4).", '5.1 Aftermath of catastrophic collisions in asteroid belt': "One possibility is that the dimming events are caused by dust thrown off in collisions between planetesimals in an otherwise unseen asteroid belt analogue (e.g., Wyatt & Dent 2002; Zeegers et al. 2014). The dust clouds created in these destructive collisions expand at roughly the planetesimals' escape velocity from the colliding bodies, eventually spreading and shearing out to form a smooth dust component in which the clumps reside. Such a scenario is a promising explanation for the star RZ Psc (de Wit et al. 2013), though in that case evidence that the underlying asteroid belt exists is given by a strong IR excess.  \nThere are several problems with this scenario as applied to KIC8462852 however. Probably the most fundamental of these is the absence of an IR excess from the smooth component. This is because for every clump we see, remembering that these were inferred to be slightly larger than the star, there should be many more that have spread out. The infrared emission from the dispersed clumps would likely sum up to a detectable level, even before counting dust produced in non-dip forming events. Moreover we should see dips from the clumps in the middle of being dispersed (i.e., dips with longer duration albeit lower optical depth), as well as dips  \nwith a continuum of depths and durations from the many different scales of planetesimal impacts that would occur. The clustering of dips at D1500 also points to these events being correlated which is hard to reconcile with this scenario, though the planetesimals in the belts could be shepherded by planets into confined azimuthal regions (e.g., Wyatt 2003; Nesvorn'y et al. 2013).", '5.2 Aftermath of giant impact in planetary system': 'A possible way around the issues in Section 5.1 is to invoke dust thrown off in a single collision, perhaps analogous to the EarthMoon system forming event (Jackson & Wyatt 2012). In this case there need not be an underlying asteroid belt, as the collision could be between planets whose orbits recently became unstable, or between growing planetary embryos. Such events are expected to result in strong IR excesses (e.g. Jackson & Wyatt 2012; Genda et al. 2015), and are indeed seen in systems such as HD 172555 where giant impacts are the favored explanation (Lisse et al. 2009). In this scenario, the putative collision would need to have occurred between the WISE observation taken in Kepler Q5 and the first large dip at D800. The dip at D1500 is then interpreted as the same material seen one orbit later, with the ∼ 750 day period implying an orbit at ∼ 1 . 6 AU. The difference in the dip structure from D800 to D1500 could arise because the clump(s) created in the original impact are expanding and shearing out. This scenario therefore predicts that KIC 8462852 may now have a large mid-IR excess, but the most recent IR observations taken in 2015 January with Spitzer IRAC show no significant excess for KIC 8462852 (Marengo et al. 2015). However, non-detection of an excess would not necessarily rule this scenario out, as the dust levels derived in Section 4 (which account for the dust seen passing in front of the star) were shown to be consistent with a non-detection. A more robust prediction is that future dimming events should occur roughly every 750 days, with one in 2015 April and another in 2017 May.  \nTwo new issues arise with this scenario however. Firstly, if the period of the orbiting material is a few years, what is the origin of the two small 0 . 5 % dips seen in the first few hundred days (D140 and D260; Table 1), and why did they not repeat 750 days later? It is a concern that these could require the existence of an outer planetesimal belt, which may contradict the lack of infrared emission to this star. Perhaps more problematic is the probability that this star (of unknown age) should suffer such an event that occurs within a few-year window between the WISE observation and the end of the prime Kepler mission, and that the geometry of the system is such that material orbiting at ∼ 1.6 AU lies almost exactly between us and the star. Taking this few year window, the main sequence lifetime, and an optimistic estimate for the scale height of giant impact debris, and the number of Kepler stars observed, this suggests that every star would have to undergo 10 4 such impacts throughout its lifetime for us to be likely to witness one in the Kepler field. Thus, while this scenario is attractive because it is predictive, the periodicity argument may be inconsistent, and the probability of witnessing such an event may be very low (though of course difficult to estimate).', '5.3 Dust-enshrouded planetesimals': "Scenarios in which the clumps can be long-lived are attractive because they suffer less from being improbable. Thus, one possibility is that the clumps are held together because they are in fact themselves orbiting within the Hill sphere of large planetesimals.  \nThey can therefore be thought of as planetesimals enshrouded by near-spherical swarms of irregular satellites, which are themselves colliding to produce the observed dust. This scenario is therefore analogous to that suggested for the enigmatic exoplanet Fomalhaut b (Kalas et al. 2008; Kennedy & Wyatt 2011), which borrows from the irregular satellites seen in the Solar System (e.g. Jewitt & Haghighipour 2007; Bottke et al. 2010). This scenario suffers from several problems. First, the observed dips already require multiple large planetesimals. Unless these all orbit within the same plane to a high degree (i.e., to within a few stellar radii), there must be many more large planetesimals which never (or have yet to) pass in front of the star. Debris disks with low levels of stirring are theoretically possible (Heng & Tremaine 2010; Krivov et al. 2013). However, these low stirring levels require the absence of large planetesimals which through mutual interactions would stir the relative velocities to their escape speeds. This is in addition to the problem of filling the Hill sphere of such planetesimals almost completely with dust. This may be reasonable if the planetesimals are embedded in a belt of debris. However, that would incur the problem of the lack of infrared excess. The question also remains why the D1500 events are so clustered, and why there are several deep dimming events and no intermediate ones. A population of planetesimals should have a variety of inclinations with respect to our line of sight, so they should pass in front of the star at a range of impact parameters and cause a range of dip depths.  \nA related scenario is that the planetesimals are surrounded by large ring systems, similar to that invoked to explain the ∼ 50 day dimming event seen for 1SWASP J140747.93-394542.6 (normally called 'J1407', Mamajek et al. 2012; van Werkhoven et al. 2014; Kenworthy & Mamajek 2015). In that case however, a single relatively time-symmetric dimming event was seen, whereas KIC8462852 has multiple asymmetric events. Thus, a single ringed planet(esimal) would not reproduce the observed light-curve, and a scenario with multiple ringed-planetesimals would be essentially the same as the irregular satellite scenario above.", '5.4 A family of objects on a comet-like trajectory': "One of the scenario independent constraints considered in Section 4 was the presence of light-curve gradients as large as 0.5 d -1 , which results in an upper limit of 13 AU for the clumps' semi-major axis assuming optically thick clumps (Figure 10). However, the star is never completely occulted, so this estimate should be corrected for the optical depth of the clump τ . That is, the steepness of the gradient is diluted either by flux transmitted through a large optically thin clump (or by unocculted parts of the star for an optically thick small clump). Assuming τ = 0 . 2 the velocity estimate given by the gradients is then 5 times higher than assumed in Section 4; this would predict a more realistic minimum transverse velocity of ∼ 50 km s -1 to cause the observed gradient, which for a circular orbit yields a maximum semimajor axis of a = 0 . 5 AU. While this estimate is uncertain, for example because of the unknown optical depth structure of the different clumps, this highlights the possibility that the material may be moving so fast that the velocity for a circular orbit is inconsistent with the non-repetition of the events.  \nOne solution to this problem is that the orbits need not be circular. That is, we could be seeing material close to the pericenter of a highly eccentric orbit, reminiscent of comets seen in the inner Solar System at pericenter (Marsden 1967; Sekanina 1984). Comets around other stars have also been detected, the first of these being found around Beta Pictoris (Smith & Terrile 1984; Lagrange-Henri et al. 1989; Beust et al. 1990). We therefore envision a scenario in  \nFigure 13. Size vs. pericenter parameter space for high eccentricity cometlike orbits. Dotted lines show lower limits on the clump sizes from the dip depths. The dashed line is the outer limit set by the light curve gradient, noting that this limit decreases with decreasing optical depth, e.g., the limit would be at a pericentre that is 25 times smaller than that plotted if the clumps have optical depth of 0.2 (line not shown in figure). The dot-dashed line is where the clump radius equals the pericenter distance, though the clumps could exist above here if they are elongated along the orbital direction. The solid lines are of constant dip duration.  \n<!-- image -->  \nwhich the dimming events are caused by the passage of a series of chunks of a broken-up planetesimal on a comet-like orbit. That planetesimal may have been analogous to what we refer to in the Solar System as a 'comet', in which case it could be volatile-rich and may have broken up as a result of thermal processes. However, it may alternatively have closer analogy with Solar System asteroids in having a more refractory composition, which might require non-thermal processes such as tidal disruption to break it up. The disruption mechanism and composition of the planetesimal are not defined for this scenario, just its orbit which is comet-like, and so we refer to it here-on as a 'comet-like' without bias to their origin or physical make-up. Regardless of its disruption process, the resulting chunks would have to have since spread around the orbit, and may be continuing to fragment to cause the erratic nature of the observed dips.  \nTo assess this scenario, Figure 13 revisits the clump - orbit parameter space of Figure 10 (discussed in Section 4), but now uses the pericenter of the clump's orbit instead of its semi-major axis. The orbits are assumed to be highly eccentric ( e ≈ 1 ), with the dips arising from material close to pericenter, so that their orbital velocity is roughly √ 2 times the circular Keplerian velocity at that distance. The limits from the dip depths and light-curve gradient are again shown, as are lines of constant dip duration. The planetesimal Hill radius lines are not shown, because they are not applicable to the cometary scenario considered here, though these would be slightly modified versions of those in Figure 10 (see eq. B5 of Pearce & Wyatt 2014). In general, the main change compared with Figure 10 is that the higher orbital velocity relaxes the constraints on how far out the clumps can be orbiting. However, as mentioned above, if the clumps are optically thin (as opposed to optically thick as assumed in Figure 13) the constraint from the light curve gradient may be more stringent. For example, decreasing the optical depth to 0.2 would result in a transverse velocity of 50 km s -1 (see above), thereby moving the light curve gradient constraint on the upper limit from 26 AU closer to 1 AU.  \nThe proximity of the comet-like clump to the star when causing the dip does not present a problem for this scenario, as it did when the clump was on a circular orbit. This is because the pericenter distance does not necessarily bear any relation to the period with which the comet-like fragments return to pass in front of the star. That period is set by the semimajor axis which has the same constraint as shown on Figure 10, and there is no such constraint on the pericenter in Figure 13. Thus the point of note from Figure 13 is that the pericenter could be significantly within 1 AU. Closer pericenters are favored both because this geometry results in a higher probability of the clumps occulting the star along our line-of-sight, and because of the greater opportunities for fragmentation of the bodies. The temperatures of comets (i.e., with volatiles) at such close proximity to the star ( > 410 K) would render them susceptible to thermal stresses. The existence of multiple super-Earth planets orbiting < 1 AU from many main sequence stars also points to the possibility that the body could have been tidally disrupted in a close encounter with one such planet. It is even possible that the body came close enough to the star for tidal disruption in the absence of other considerations; e.g., a comet similar to Halley's comet would fall apart by tidal forces on approach to within 3-7 stellar radii (0.02 - 0.05 AU). By contrast, a rocky body would require a closer encounter to tidally disrupt.  \nFor close pericenters it is important to point out that while the constraint is discussed in terms of the clump's radius, the clump can not in fact be spherical at that size. Figure 13 shows a blue dotdashed line where the 'clump radius' is the same as the pericenter distance. At such proximity, the clump could not be elongated in the radial direction, but could only be elongated azimuthally along the orbit. In fact, this mostly linear clump structure is the correct way to visualize debris from the breakup of a comet or planetesimal. The small velocity kicks (from fragmentation or tidal disruption) would cause a small dispersion in semimajor axis for material in the clump, and the resulting differential orbital motion causes the material to spread around the orbit. These small kicks do not significantly change the periastron distance or the orbital inclination angle.  \nThis scenario is attractive, because comets are known in the Solar System to have highly eccentric orbits and disrupt for various reasons near pericenter, and infalling comets are the most robust explanation for the falling evaporating body (FEB) phenomenon seen around many nearby A-type stars (e.g. Kondo & Bruhweiler 1985; Beust et al. 1990; Welsh & Montgomery 2013; Kiefer et al. 2014). Also, since fragments of the comet family would all have very similar orbits, this mitigates the problem noted in Section 5.1 that the detection of multiple transits may require orders of magnitude more clumps to be present in the system. Instead, the observed clumps may be essentially in a single orbit which is that of the progenitor, and that orbit happens to be preferentially aligned for its transit detection. That is, it is not excluded that we have observed all the clumps present in the system. While a quick look at Figure 11 suggests that the lack of infrared excess might still be problematic for the closest pericenters (noting that σ tot also needs to be increased by √ 2 due to the higher transverse velocity at pericenter in Equation 5), in fact that is not necessarily the case. Rather, in that figure we assumed that the clumps were present at the given distance at all times, whereas the clumps in the comet-like group scenario were at much larger separation from the star at the time of the WISE observations. The total mass of the fragmented body was considered in Section 4, but since the clumps can be closer to the star in this scenario, and are moving faster than circular Keplerian velocity, a better minimum mass estimate for clumps seen at a  \npericenter of 0.1 AU is ∼ 3 × 10 19 g. Again, the size distribution and any material not contributing to the observed dips will increase this minimum mass, perhaps by a factor of 100, leading to a more realistic parent body mass of 3 × 10 21 g , consistent with a rocky body ∼ 100 km in diameter.  \nIt remains to be shown that this model can explain the more detailed structure of the light-curves. Some potential positives are that the clustered nature of the dips could be explained by subsequent fragmentation of a large fragment from an earlier breakup. The smaller dips could also potentially be explained by smaller fragments which may also be expected to receive larger kicks during fragmentation. However, the structure of individual clumps may be problematic. For example, a fairly generic prediction of transits of comet-like bodies may be that their light-curves show signs of their tails. The light-curve expected for a typical event then has a relatively fast ingress as the head of the comet passes in front of the star, but a slower egress as the tail passes (e.g. Lecavelier Des Etangs et al. 1999; Rappaport et al. 2012). However, the D800 event shows the opposite (see panel 'c' in Figure 1). Possible resolutions of this issue are that the D800 comet fragment received a large kick with an orientation that sheared it out in such a way to form a 'forward tail'. Such forward comet tails produced by the fragments being kicked toward the star have been studied in the literature, but require the grains in the tail to be large enough to overcome the effects of radiation pressure (Sanchis-Ojeda et al. 2015). Alternatively, this event could be comprised of two dips superimposed to have the appearance of a forward tail. While several issues remain to be explored, of the scenarios considered we conclude that a cometary-like group of bodies seems most consistent with the data at hand.", '6 SUMMARY AND CONCLUSIONS': "In this paper, we have shown that KIC 8462852 is a unique source in the Kepler field. This otherwise seemingly normal F star undergoes erratic and completely unpredictable dips in flux ranging from glyph[lessorsimilar] 1% to more than 20%. Most of the approximately 7 dips observed before D1500 have fairly smooth, but unexplained, dip profiles that are each several days long. The D1500 sequence lasts continuously for at least 80 days, but the majority of that time is spent with the flux depressed by less than ∼ 2%.  \nWe have conducted numerous follow-up investigations of the star and its environment, including spectroscopy, adaptive optics imaging, construction of a spectral energy distribution, generation of a Fourier transform and a sonogram using the Kepler time series, and examination of ground-based photometry. Our analysis characterizes the object as both remarkable (e.g., the 'dipping' events in the Kepler light curve) and unremarkable (ground-based data reveal no deviation from a normal F-type star) at the same time.  \nAn extensive set of scenarios has been presented to explain the occurrence of the dips, most of which are unsuccessful in explaining the observations in their entirety. Of these, the scenarios invoking intrinsic variability, such as the Be star framework, were deemed unlikely, but they are not entirely ruled out as a plausible option to explain the dips. However, we pointed out that the relatively low T eff and lack of H α emission and IR excess in KIC 8462852 are not suggestive of Be-star activity.  \nA broad range of scenarios for the dipping behavior that involve occultation by circumstellar dust clumps was considered. Among these, we find that the break-up of one or more massive exocomets (or planetesimals on comet-like orbits) provides the  \nmost compelling explanation consistent with the data in hand. The required mass of the original body may have been in excess of 3 × 10 21 grams (only ∼ 0.3% the mass of Ceres, and perhaps ∼ 100 km in diameter).  \nWe can envision a scenario in which a barrage of bodies, such as described above, could be triggered by the passage of a field star through the system. And, in fact, as discussed in Sect. 2.3, there is a small star nearby ( ∼ 1000 AU; Section 2.3) which, if moving near to KIC 8462852, but not bound to it, could trigger such a barrage into the vicinity of the host star. On the other hand, if the companion star is bound, it could be pumping up comet eccentricities through the Kozai mechanism. Measuring the motion/orbit of the companion star with respect to KIC 8462852 would be telling as to whether or not they are physically associated, and we could then be better able to make assessments about the timescale and repeatability of comet showers based on bound or unbound star-comet perturbing models.  \nContinuing observations of KIC 8462852 should aid in unraveling the peculiar dips in its light curve. First and foremost, longterm photometric monitoring is imperative in order to catch future dipping events. It will be helpful to know whether such observations reveal continued, possibly periodic dips, or no further dips. If the dips continue, it will be important to search for a clear periodicity, and to look for changes in depth or shape. To completely solidify the hypothesis that the dips are due to dust, observations should study the wavelength dependence of the obscuration soon after a new dip is discovered. In the case of a family of giant comet-like bodies there presumably should be at least a few events similar to those seen with Kepler over the next decade. However, if the cometlike objects actually populate a very long eccentric orbit (i.e., that of the original planetesimal), the material may be spread out around that orbit, and future dippings events could continue to appear over hundreds of years.  \nSeveral of the proposed scenarios are ruled out by the lack of observed IR excess (Section 2.4), but the comet/planetesimal fragments scenario has the least stringent IR constraints. In the comet scenario, the level of emission could vary quite rapidly in the nearIR as clumps pass through pericenter (close to the time they are transiting) and are shedding new material. If the system is currently in the aftermath of a giant impact, there could be a semi-steady increase in IR flux over years/decades. The WISE observations were made in Q5, and assuming that an impact occurred in Q8 (D800, Section 5.2), detecting the IR emission from such an impact is still a possibility in the future. The only Spitzer IRAC observation of KIC8462852, taken in January 2015, showed a marginal, but below 3σ , excess at 4 . 5 µ m, disfavoring the impact scenario (Marengo et al. 2015). Continued monitoring in the IR will allow us to firmly distinguish between the giant-impact and cometary-group scenarios.  \nIn summary, it will require some observational skill and patience to find the next dipping event from this object using groundbased observations. As we pointed out, the source spent a rather small fraction of its time during the 4-year Kepler mission with dips of greater than 2%. Nonetheless, the key to unraveling the mysterious dips will require such observations.", 'ACKNOWLEDGMENTS': "We thank Jason Wright and Jason Curtis for fruitful discussions on the object. We further acknowledge Planet Hunter user 'Exoplanet1' for their contributions to the discussion of this object. We  \nare grateful to Sherry Guo and Bhaskar Balaji for running an automated search through the Kepler set to find other similar dippers. We acknowledge Mike Jura for very insightful comments about the required mass of the body that is the origin of the obscuring material. We thank Josh Carter for pointing out the possible 48-day periodicity. We appreciate the efforts of Jeff Coughlin, Jon Jenkins, and Jeffrey Smith for taking a careful look at the raw Kepler photometry to decided if it was all good, i.e., not artifacts. We thank Mark Everett and Lea Hirsch for making the DSSI observations. We thank Huan Meng, Massimo Marengo, and Casey Lisse, for insightful comments related to the IR excess. We are grateful for thoughtful discussions with members of the Kepler Eclipsing Binary Working Group, and attendees of the K2 SciCon 2015. Last but not least, we are grateful for the anonymous referee's comments to help improve the paper.  \nTSB acknowledges support provided through NASA grant ADAP12-0172 and ADAP14-0245. MCW and GMK acknowledge the support of the European Union through ERC grant number 279973. The authors acknowledge support from the Hungarian Research Grants OTKA K-109276, OTKA K-113117, the Lendulet2009 and Lendulet-2012 Program (LP2012-31) of the Hungarian Academy of Sciences, the Hungarian National Research, Development and Innovation Office - NKFIH K-115709, and the ESA PECS Contract No. 4000110889/14/NL/NDe. This work was supported by the Momentum grant of the MTA CSFK Lendulet Disk Research Group. GH acknowledges support by the Polish NCN grant 2011/01/B/ST9/05448. Based on observations made with the Nordic Optical Telescope, operated by the Nordic Optical Telescope Scientific Association at the Observatorio del Roque de los Muchachos, La Palma, Spain, of the Instituto de Astrofisica de Canarias. This research made use of The Digital Access to a Sky Century at Harvard (DASCH) project, which is grateful for partial support from NSF grants AST-0407380, AST-0909073, and AST-1313370. The research leading to these results has received funding from the European Community's Seventh Framework Programme (FP7/2007-2013) under grant agreements no. 269194 (IRSES/ASK) and no. 312844 (SPACEINN). We thank Scott Dahm, Julie Rivera, and the Keck Observatory staff for their assistance with these observations. This research was supported in part by NSF grant AST-0909222 awarded to M. Liu. The authors wish to recognize and acknowledge the very significant cultural role and reverence that the summit of Mauna Kea has always had within the indigenous Hawaiian community. We are most fortunate to have the opportunity to conduct observations from this mountain. KS gratefully acknowledges support from Swiss National Science Foundation Grant PP00P2 138979/1. HJD and DN acknowledge support by grant AYA2012-39346-C02-02 of the Spanish Secretary of State for R&D&i (MINECO). This paper makes use of data from the first public release of the WASP data (Butters et al. 2010) as provided by the WASP consortium and services at the NASA Exoplanet Archive, which is operated by the California Institute of Technology, under contract with the National Aeronautics and Space Administration under the Exoplanet Exploration Program. This publication makes use of data products from the Widefield Infrared Survey Explorer, which is a joint project of the University of California, Los Angeles, and the Jet Propulsion Laboratory/California Institute of Technology, and NEOWISE, which is a project of the Jet Propulsion Laboratory/California Institute of Technology. WISE and NEOWISE are funded by the National Aeronautics and Space Administration. This research made use of the SIMBAD and VIZIER Astronomical Databases, operated  \nat CDS, Strasbourg, France (http://cdsweb.u-strasbg.fr/), and of NASA's Astrophysics Data System."}