arXiv:1804.08351v2 [astro-ph.SR] 1 Jul 2018 ∗ cl eatutr fti yehsbe siae tas at estimated as been little has type this of megastructure scale 2016 status II type thereby star, towards host scale their ( Kardashev of the power radiation climbing the frac- significant of a into tion tap advanced principle in extremely could civilizations planets, dismantled from material ohsasadselrrmat aeas enconsid- been around also built have (e.g. are remnants ered stellar shades and large stars single both or radiation. sparse rings stellar where the swarms, designs, of astro-engineering fraction ambitious small Less satel- a each spherical with absorbing dense, satellites, lite a absorbing as of rather shroud but or swarm shell, solid a as visioned -al [email protected] E-mail: adse 1964 Kardashev rpittpstuigL using 2018 typeset 3, Preprint July version Draft ycntutn yo pee( sphere Dyson a constructing By EIWT AA H BEVTOA INTRSO ERYCOMP NEARLY OF SIGNATURES OBSERVATIONAL THE GAIA: WITH SETI ; rgte l 2016 al. et Wright sta h srmti ouinhsbe opoie ya unseen an by compromised ( been dwarf has white massive solution rather astrometric discrepan the the for para that explanation Gaia param is plausible its stellar a RAVE, the from by confirms derived inferred largely that distance analysis our twice TYC While about dwarf excess. is p F-type infrared we object late population, this the outlier stars, of this these distance find in of We objects one 5. of for spectroscopy properties Release f the expected Data Dyson-sphe on type RAVE light the for and of shed search 1 discrepancies distance limited Release display a G Data indeed stars combining out Gaia by to carry common time, and stars current technique this the He of at surveys. ings possible ground-based also from f is distances allow spectrophotometric type years signat this coming waste-heat in of on will based search searches mission to Gaia complementary are the may that criterion, fraction s this covering estimate Using distance high spectrophotometric with distance. a sphere with Dyson object a subluminous in enshrouded star A Keywords: releases. data ∼ uen1977 Suffern 10 2 r( yr .TeDsnshr stpclynten- not typically is sphere Dyson The ). 1. eateto hsc n srnm,UpaaUiest,B University, Uppsala Astronomy, and Physics of Department rsrn adeg2013 Sandberg & Armstrong rsoy–ds,etnto tr:idvda TC6111-1162-1 (TYC individual stars: – extinction dust, – troscopy xrtretilitliec tr:dsacs–mtos data methods: – distances stars: – intelligence Extraterrestrial INTRODUCTION A T E ; ,adtetm obidasmall- a build to time the and ), tl ATX .1.0 v. AASTeX6 style X rgte l 2014a al. et Wright rkZackrisson Erik yo 1960 Dyson ≈ 1 M ⊙ ; nvriyo edleg Germany Heidelberg, of University 1 .Ti cnrocnb ute etdtruhucmn Gaia upcoming through tested further be can scenario This ). Osmanov ∗ nra .Korn J. Andreas , u of out ) .Pro- ). oansReiter Johannes ABSTRACT ue nteMlyWy( Way Milky with- struc- the such out, for in individual carried evidence tures compelling for been any searches already unraveling out have of spheres number Dyson a (SETI), ligence asadln-ag omncto (e.g. communication long-range supercom- habi- and artificial include of tats powering energy the collected propulsion, vessel the puting, for uses posed edn nteszsaddsrbto fteobscuring the of distribution and sizes the on any astro-engineering pending 2015 up of al. et turn type Zackrisson to this failed of ( Dyson also signatures by have tell-tale powered technology colonization sphere intergalactic even galaxy-spanning or in engaging civilizations for searches lar 2000 2013 1999 Sandberg ni 1999 Annis ntefil fSacigfrEtaersra Intel- Extraterrestrial for Searching of field the In 1 ; ). 2 ngrWehrhahn Ansgar , uau&Nsiua2004 Nishimura & Jugaku e edsustemrt n shortcom- and merits the discuss we re, rsa nrrdwvlnts limited A wavelengths. infrared at ures x55 E712 psl,Sweden Uppsala, 20 SE-751 515, ox 1116-.Tespectrophotometric The 6111-1162-1. rnal opeeDsnshrs To spheres. Dyson complete nearly or ; gicnl necs fisparallax its of excess in ignificantly rgte l 2014b al. et Wright itneetmtso hsobject this of estimates distance t tr n h spectrophotometric the and eters ; eetflo-phigh-resolution follow-up resent lx n hr sn detectable no is there and llax, ecniae ntesml of sample the in candidates re rod2005 Arnold iaycmain osbya possibly companion, binary aietisl sa optically an as itself manifest i aalxdsacswith distances parallax aia nlss–mtos spec- methods: – analysis ; rDsnshr searches Dyson-sphere or ) art 2015 Garrett htasalfato of fraction small a that EEDSNSPHERES DYSON LETE 1 ls 1985 Slysh ; rsrn Sandberg & Armstrong ; , argn2009 Carrigan a ; rffihe l 2015 al. et Griffith ; lo 2017 Olson ; ioeve al. et Timofeev uen1977 Suffern .Simi- ). .De- ). ; ; 2 Zackrisson et al. satellites, these structures may give rise to conspicu- /far-infrared excess through a comparison of predicted ous temporal variations in the apparent brightness of and observed fluxes at mid/far-IR wavelengths from in- the affected star (e.g. Arnold 2005; Wright et al. 2016; frared telescopes like WISE, IRAS, Spitzer, Akari and Osmanov 2016). However, most searches so far have fo- Herschel. cused on the overall, long-term changes in the spectral In section 2, we describe the method and its limita- energy distribution of the star that the satellites may tions. In section 3 we apply this method to the ≈ 2×105 induce. stars in common between Gaia DR1 and RAVE DR5 Following the AGENT formalism of Wright et al. in search of objects with signatures that match those (2014a), the energy budget of a Dyson sphere can be expected for partial Dyson spheres with high covering described as: fractions. As an example of the type of object that is α + ǫ = γ + ν, (1) identified by this method, we discuss the properties of TYC 6111-1162-1 in section 4, and present follow-up ob- where α represents the radiation energy collected by servations to confirm the stellar parameters derived by the Dyson sphere, ǫ the energy produced by other RAVE. Section 5 discusses the prospects of applying the means, γ the thermal waste heat and ν the energy lost method to much larger data sets in the future. Section 6 through non-photonic means (e.g. neutrinos, gravita- summarizes our findings. tional waves). Under the assumption that the collected starlight overwhelms all other forms of energy generated 2. DISCREPANCIES BETWEEN and that non-thermal losses are negligible, equation (1) SPECTROPHOTOMETRIC AND PARALLAX simplifies to α ≈ γ. In observational terms, even a minor DISTANCES FOR PARTIAL DYSON SPHERES fraction of collected starlight can then generate a very Through the use of stellar-atmosphere modeling and pronounced infrared excess in the spectral energy distri- theoretical stellar isochrones in conjunction with ob- bution of the star, as long as the Dyson sphere operates served photometric fluxes and spectroscopic data at at a temperature in the 50–1000K range. Consequently, optical/near-infrared wavelength, it is possible to ac- most Dyson-sphere searches have focused on detecting curately determine stellar parameters, dust-attenuation this infrared signature. However, in case where losses factors and spectrophotometric distances Dspec, or con- due to non-thermal photons are significant (i.e. where ν versely spectrophotometric parallaxes πspec, for various is non-negligible), where the energy is somehow stored types of stars (e.g Binney et al. 2014). This procedure rather than radiated away or where the waste heat radi- basically assumes that once the observed optical flux ation is anisotropic (and not emitted in our direction), fopt has been corrected for dust, the relation between the infrared signature may be weaker than usually as- the intrinsic optical luminosity Lopt and its spectropho- sumed. Here, we outline a method to search for partial tometric distance is: Dyson spheres with weak or even absent waste heat sig- Lopt natures. Fopt = , (2) 4πD2 Since a star surrounded by a partial Dyson sphere with spec large (but non-unity) covering fraction is expected to ap- which implies pear unusually faint for its type at optical/near-infrared 1/2 Lopt wavelengths, its spectrophotometric distance is expected Dspec = (3) to be overestimated. A trigonometric parallax measure- 4πFopt  ment, as provided by Hipparcos or Gaia would, on the However, consider a star enshrouded by a Dyson sphere other hand, still be able to provide an accurate dis- which covers a fraction fcov of the stellar surface tance estimate. Hence, it should be possible single out area. If this obscuring structure blocks the emerging tentative Dyson-sphere candidates by comparing paral- optical/near-IR light like a grey absorber (i.e. all wave- lax distances to their spectrophotometric counterparts. lengths throughout this interval are equally affected), With future Gaia data releases, both distance estimates the only impact on the observed flux would be a general can be obtained from Gaia data, but even with Gaia dimming, whereas the detailed shape of the spectrum at Data Release (DR) 1 (Gaia Collaboration et al. 2016) it these wavelengths would be unaffected. The true dis- is possible to carry out a search of this kind by com- tance would then be given by D′: bining the parallax distances provided by the Tycho- L (1 − f ) 1/2 Gaia Astrometric Solution (TGAS) with spectrophoto- D′ = opt cov . (4) metric distances from ground-based surveys like the Ra-  4πFopt  dial Velocity Experiment (RAVE; Kunder et al. 2017) Combining eqs. 3 and 4, we get: or GALAH (Martell et al. 2017). Candidates for Dyson ′ 1/2 spheres identified this way can also be tested for mid- D = Dspec(1 − fcov) (5) SETI with Gaia 3

Figure 1. Distribution of fcov for a) the whole Gaia DR1-RAVE DR5 sample; b) the subset of the Gaia DR1-RAVE DR5 sample with the smallest parallax errors (trigonometric parallax error ≤ 10%, spectroscopic parallax error ≤ 20%).

Hence, for a Dyson sphere which obscures a fraction spectrophotometric fit also allows for a prediction of the 1 > fcov > 0 of the intrinsic optical radiation, the intrinsic stellar flux at infrared wavelengths that can true distance D′ would be smaller than that inferred by be directly compared to observed fluxes with mid/far- standard spectrophotometric distance-estimation tech- infrared survey data in search of any potential infrared niques. Even in this case, a measurement of the trigono- excess. metric parallax distance Dtrig should still be able to re- It should be noted that the covering fraction fcov used trieve the true distance, since this technique is not di- in the above formalism does not need to represent the rectly affected by the brightness of the star. If we adopt global covering fraction of the star, but simply the ob- ′ Dtrig ≈ D and use the relation between trigonometric scured fraction of the face of the star facing the observer. parallax πtrig and distance D (π = 1/D, as long as the Hence, the principle not only applies to isotropic Dyson fractional parallax error is ≤ 0.2, Bailer-Jones 2015), spheres, but also different types of stellar engines (e.g. then we can express fcov as: Badescu & Cathcart 2000) and large transiting megas- 2 2 tructures (Wright et al. 2016) that provide anisotropic Dtrig πspec fcov =1 − =1 − (6) coverage of the star. However, stellar engines that re- Dspec   πtrig  flect some fraction of the light back at the host star – However, the square-root dependence in eq. (5) makes like the Shkadov thruster (Shkadov 1987; Forgan 2013)– the expected difference between the two distances very would likely affect its thermal equilibrium and may pre- small unless the star is close to being completely en- vent stellar parameters from being accurately retrieved shrouded, i.e. fcov is close to unity. For example, using standard spectral-analysis techniques. fcov ≈ 0.1 would correspond to a distance discrepancy of just ≈ 5%, which is undetectably small. In fact, it 3. A SEARCH FOR PARTIAL DYSON-SPHERE takes fcov > 0.75 to produce a factor of > 2 discrepancy, CANDIDATES IN GAIA DATA RELEASE 1 which may be more in the range of what one can hope For this pilot project, we have combined Gaia DR1 to achieve given current errorbars. with the RAVE DR5 dataset. Gaia DR1 contains In practice, to single out fcov > 0.75 candidates from 5-parameter astrometric data for 2 million objects, a large sample of objects requires that distances or par- whereas RAVE DR51 has spectroscopic data for over allaxes can be measured with relative errors of no more 4.8×105 stars. Of these, about 2.3×105 objects feature than ≈ 10–20%. This level of parallax precision is rou- in both catalogs. tinely achieved with space-based missions like Hipparcos If one attempts to infer fcov using eq. (6) for the full or Gaia, and the Gaia mission is eventually expected to data set, one arrives at the fcov distribution shown in 8 provide this level of parallax precision for ∼ 10 main the left panel of Figure 1, with extremely large num- sequence stars (Luri et al. 2014). The precision of spec- bers of stars at both very high and very low covering trophotometric parallaxes is typically somewhat lower, fractions. The reason for this is simply that many have with relative errors at the ≈ 20–30% level (Binney et al. distance errors much too large to allow meaningful con- 2014) in the case of the RAVE survey, although the sub- set of stars with the very best estimates are still likely to be useful. 1 Stellar ages and masses are not provided in RAVE DR5, there- In cases where a high-fcov candidate is identified, the fore these quantities are extracted from RAVE DR4 4 Zackrisson et al.

10 1% of stellar luminosity reprocessed into 300 K black body Dyson sphere 12 AB mag 14 Observation Model 16 5 10 15 Wavelength [ m] Figure 2. Density of fcov as a function of Teff for all stars Figure 3. Comparison of observed and predicted broadband after the cut. The position of TYC 6111-1162-1 is indicated optical-to-IR fluxes for TYC-6111-1162-1, based on the stel- by a red cross. The numbers show the included fraction of lar parameters and spectrophotometric parallax presented in stars for each contour. RAVE DR5. Blue circles indicate the observed optical (BV), near-IR (JHKs) and mid-IR (WISE W1, W2, W3) fluxes (ex- pressed in AB magnitudes, with error bars smaller than the straints on fcov. By pruning the sample, and only keep- symbols used) and red diamonds the corresponding fluxes ing objects with parallax errors at ≤ 10% for Gaia and predicted by the PARSEC 1.1 isochrones given the best- ≤ 20% for RAVE (a good compromise between error fitting model (stellar parameters, interstellar extinction and spectrophotometric parallax) presented for this star in RAVE and sample size, given the current data), one arrives at DR5. The black solid curve is a black body fit based on the the much narrower distribution (featuring 8441 stars) bolometric luminosity and effective temperature of this stel- shown in the right panel of Figure 1. Once this cut lar model. As seen, there is reasonable agreement between ∼ observed and modelled fluxes at all wavelengths plotted, and has been applied, a small group of 75 stars (i.e. 1% no indication of any significant IR excess. The dashed curve of the full sample) at fcov > 0.9 clearly stands out in illustrates the strong IR excess expected even in a case where the distribution. Under the assumption that parallax just 1% of the bolometric stellar luminosity would be repro- errors are Gaussian, one would naively expect only a cessed into black body radiation by a Dyson sphere operating at 300 K. fraction . 0.04% of the sample (i.e. . 4 objects) to spu- riously appear at fcov > 0.9. However, the probability distribution for the spectrophotometric distance mod- 2012) isochrones on which the RAVE fits were based. ulus is in reality far more complicated and frequently Occasionally, these stellar parameters yield nonsensical features multiple peaks (e.g. Binney et al. 2014). All results when compared to the isochrones (e.g. a mass types of stars are moreover not equally well-described too high for the isochrone), which suggests some sort of by current stellar-evolution tracks or stellar-atmosphere glitch in the fitting. Once these are also removed, we models, which makes it extremely difficult to reliably es- are left with 6 outliers at fcov > 0.7 out of 8365 stars. timate the expected number of spurious outliers in the To investigate these stars further, using additional ob- fcov distribution. In the following, we therefore take a servations from APASS (optical), 2MASS (near-IR) and closer look at the stars that end up in this high-fcov tail WISE (mid-IR), we compare the observed Spectral En- in an attempt to further constrain their nature. ergy Distribution (SED), to the theoretical predictions Since the most interesting outliers would be those that from the PARSEC isochrones given the best-fitting stel- belong to a spectral type for which current models are lar parameter. Four stars show conspicuous inconsis- generally expected to produce highly reliable results, tencies between the two, either indicating an incorrect we are mainly interested in main-sequence FGK dwarfs. RAVE distance modulus or a problem with the data However, the vast majority of objects in the conspicu- in the observed catalogs. The only two stars remain- ous fcov > 0.9 peak in the right panel of Figure 1 have ing are TYC 7169-1532-1 and TYC 6111-1162-1, nei- inferred values for the surface gravity log g < 2 (g in ther of which show any obvious mid-IR excess above cgs units) and are consequently classified as low-gravity the predictions from the stellar modeling. As discussed giants. Since it has been shown that the RAVE spec- by McMillan et al. (2017), discrepancies between spec- trophotometric distances systematically tend to then trophotometric and astrometric parallaxes exhibit a sig- overestimate the distances to such stars (Binney et al. nificant effective-temperature dependence, likely reflect- 2014), albeit not to the level seen here, we will not con- ing systematic problems and degeneracies in the spec- sider these further. As an additional consistency check trophotometric fitting procedure. In Figure 2 we plot of the remaining outliers, we check the best-fitting mass, the inferred fcov for our sample as a function of Teff . age and metallicity against the PARSEC (Bressan et al. With its Teff of ≈6200 K, TYC-6111-1162-1 lies in the SETI with Gaia 5 effective-temperature regime for which distance agree- The model atmospheres used and the line-synthesis code ments tend to be very good, whereas TYC 7169-1532-1 assume LTE (local thermodynamic equilibrium). Only with Teff of 7500 K is found in a region where distance very small departure from LTE (at most a few 0.01dex overestimates seem to be far more common. We there- in the log of the abundance) are to be expected in dif- fore take TYC 6111-1162-1 to be the most interesting ferential studies of stars this close to Solar parameters. individual candidate for further investigation in the cur- This technique has proven very successful for local F- rent sample. and G-type stars, with distances compatible with HIP- PARCOS parallaxes to within 5% (rms). Significantly 4. CASE STUDY: TYC 6111-1162-1 larger offsets are a clear sign of stellar anomalies. In Figure 3, we compare the observed apparent flux of Starting from a longer list of iron lines that was short- TYC 6111-1162-1 at optical (BV ), near-IR (JHKs) and ened by an iterative 2σ-clipping procedure, 26 lines mid-IR (WISE W 1, W 2, W 3) wavelengths to the corre- with line strengths (equivalent widths) between 28 and sponding fluxes provided by the PARSEC 1.1 isochrones 103mA˚ were used in the derivation of the final param- given the best-fitting model (stellar parameters, inter- eters. The resulting stellar parameters are: Teff = 6150 stellar extinction and spectrophotometric parallax) for ± 80 K, log g = 4.24 ± 0.1, [Fe/H] = +0.32 ± 0.03, ξ this star in RAVE DR5. As seen, this model provides = 1.2 ± 0.1 km/s (all uncertainties given as 1σ). As il- an adequate fit to the data throughout the whole plot- lustrated by Figure 4, these stellar parameters describe ted wavelength range. The observed flux shows no de- the line spectra of selected diagnostic regions very well. tectable mid-IR flux excess compared to this model – Furthermore, the final log g (of importance for the spec- in stark contrast to the expectation for waste heat ra- trophotometric distance, see below) is fully supported diation from a Dyson sphere (Dyson 1960), even in the by the pressure-sensitive wings of the Mg Ib lines, with case where just a tiny fraction of the stellar luminosity is the magnesium abundance determined from the weaker captured/reprocessed into waste heat. We exemplify the Mg I 4730 line. In summary, TYC 6111-1162-1 seems to latter argument by showing the mid-IR excess (dashed be a perfectly normal, metal-rich main-sequence star of line in Figure 3) in a case where 1% of the bolometric late F-type. luminosity of the star is re-radiated as black body ra- We can now compare our analysis with the results diation at 300 K. This scenario would raise the flux at published by the RAVE consortium. Kunder et al. 12µm (the WISE W 3 band) by more than two magni- (2017) list the DR5 results based on a SNR≈45 tudes, in direct conflict with the data, and can therefore spectrum as 6169±54 K/4.54±0.11/+0.06±0.1 be ruled out. Hence, the RAVE DR5 stellar model and (Teff /log g/[m/H]), while the data-driven RAVE- spectrophotometric parallax provides a decent fit to the on analysis returns 6039±91K/4.29±0.13/+0.06±0.07. observed SED of TYC 6111-1162-1, without any con- The agreement with both data set is good for Teff spicuous anomalies. and log g, less so for the metallicity. Other metallicity To verify or refute the identified distance discrep- indicators provided in RAVE DR5 (“Fe abund” and ancy between RAVE and Gaia, we obtained a higher- “calib met”) agree better at +0.18 and +0.16, respec- quality spectrum of this star. We took a 1h exposure tively, and it seems justified to compare to one of these of TYC 6111-1162-1 on February 9th, 2017 with the instead, as we are actually measuring iron. However, high-resolution spectrograph FIES (Telting et al. 2014) we find [Mg/Fe] = [Ca/Fe] = 0.00 ± 0.05 (based on on the Nordic Optical Telescope (NOT). FIES provides Mg I 4730 and Ca I 6572) and thus see no reason why continuous wavelength coverage from 370-830 nm (to the Ca II triplet region used by RAVE should lead to a 910nm with some gaps) at a resolving power R = 67,000. systematically underestimated metallicity. Our analysis Data reduction was performed on-site using FIEStools based on better data may indicate that the parameter (Stempels & Telting 2017). In the wavelength region space of metal-rich stars is not (yet) well-calibrated in used for the spectroscopic analysis (mostly λ ≥ 5000 A),˚ RAVE. the signal-to-noise ratio (SNR) is approximately 70 per Having determined reliable spectroscopic stellar pa- pixel. rameters, we can now use evolutionary tracks to deter- We perform a line-by-line differential analysis relative mine the stellar mass. Interpolating among MIST stellar to the Sun using tools very similar to those employed by tracks (Choi et al. 2016), the star’s mass is constrained Fuhrmann (2004). In short, Hα is used to constrain the to M = 1.34 ± 0.05 M⊙ (RAVE DR4 lists 1.17 ± 0.1 M⊙ effective temperature Teff and the ionization equilibrium as the mass estimate for essentially the same set of stel- of neutral and ionized iron (Fe I/II) to constrain the sur- lar parameters as in DR5). Using AV =0.10 following face gravity log g and the iron abundance [Fe/H]. The RAVE and interpolating in the tables of Alonso et al. microturbulence ξ is determined by requiring weak and (1995) to estimate an appropriate bolometric correction strong iron lines to return the same abundance value. (BCV =−0.07), we arrive at a spectrophotometric dis- 6 Zackrisson et al. TYC6111 Halpha region TYC6111 MgIb region

1.0 1.0

0.8 0.8

0.6 0.6

normalized flux 0.4 normalized flux 0.4

0.2 0.2

6540 6550 6560 6570 6580 5150 5160 5170 5180 5190 5200 λ [Å] λ [Å] Figure 4. Spectral regions of diagnostic power for the spectroscopic analysis of TYC 6111-1162-1: Hα, a sensitive indicator of the effective temperature (left) and the Mg Ib lines (right) that react sensitively to the surface gravity. Bullets mark the observations, the full-drawn line is the best synthetic fit to the data. The overall fit quality is very satisfactory (with the exception of the Hα line core which is affected by the chromosphere and departures from LTE) indicative of a normal, metal-rich main-sequence star.

+48 tance estimate of Dspec = 238 −34 pc (πspec = 4.2 ± culations to check whether the binary nature of this sys- 0.7 mas). In RAVE DR5, the spectrophotometric par- tem by itself can explain the discrepant distance infor- allax is listed as πRAVE = 4.51 ± 0.8 mas which cor- mation. +50 responds to a distance of 222 −35 pc, while Gaia DR1 To approximately double the parallax from 4.2 mas tabulates πGaia = 8.72 ± 0.53 mas, i.e. implying a dis- (our best estimate) to 8.7 mas (Gaia’s best DR1 esti- tance only half as large as the one derived by us. The mate) requires an additional ≈4 mas wobble that is not distance discrepancy is thus confirmed. If interpreted caused by the distance-related trigonometric parallax. as caused by a partial Dyson sphere, this would convert This assumes that the spectroscopic analysis of the SB1 into fcov ≈ 0.77 and a dimming of the star by ≈ 1.6 mag spectrum is not affected by the binary companion. At at optical/near-IR wavelength. 250pc, 4mas corresponds to 1AU, that is, the secondary What can cause this ≈100% distance offset between has to move the primary by this distance as the two Gaia DR1 and the spectroscopic analyses? In practi- stars circle their common center of mass. Assuming cal terms, the spectrophotometric parallax (log πspec = circular orbits, this can be achieved by a ≈1 M⊙ star 0.5([g] − [M]) − 2[Teff] − 0.2(V − AV +BCV +0.25) with at a distance of ≈1 AU. This star cannot be a main- [X]:=log(X/X⊙)) reacts most sensitively to changes in sequence star, as it would contribute ≈25% of the pho- log g, as this parameter has the largest relative uncer- tons at optical wavelengths which would surely be visi- tainty in the logarithm. However, an unrealistically ble in the spectrum. Instead, it has to be a white dwarf large log g correction of ≈0.6 dex (a factor of 4 in g) at the high-mass end of the white-dwarf mass distri- would have to be applied to align the spectrophotomet- bution. This scenario thus explains a) the measured ric parallax with the astrometric one. By tweaking the radial-velocity amplitude (FIES vs MIKE) and b) the spectroscopic analysis the discrepancy can thus be di- difference between spectrophotometric and trigonomet- minished by a few tens of percent only. ric distance (FIES vs Gaia). It is more or less compat- Additional clues come from radial-velocity follow-up ible with the short-term radial-velocity variation which we have conducted since February 2017. The FIES spec- with an orbital velocity of the primary of 20 km/s (lower trum points to a heliocentric radial velocity of −4.66 limit) would be 0.33 km/s (0.49 km/s as measured by ± 0.5 km/s, fully compatible with the value published HERMES). All of these are modulo orbital ellipticities by RAVE (−4.728 ± 1.097km/s). However, a Magel- and projection effects due to the unknown inclination of lan/MIKE spectrum taken on 2017-08-13 gives −25.3 the orbital plane of the system with respect to the line ± 1.0 km/s instead. More recently, two measurements of sight. with Mercator/HERMES taken six days apart show a re- While this scenario is astrophysically plausible and at- turn to values closer to zero, but also indicate short-term tractive for explaining an extreme outlier as an extreme variations: −5.40 ± 0.05km/s (2018-03-17) and −5.89 binary system, it does not explain why and how the ± 0.05km/s (2018-03-23). Clearly, TYC 6111-1162-1 repeat observations compiled in Gaia DR1 (14 months is a single-lined spectroscopic binary (SB1). We don’t of data with a few dozen individual measurements) re- have nearly enough data to constrain the binary-system sulted in a consistent answer πGaia = 8.72 ± 0.53 mas parameters, and thus only make order-of-magnitude cal- in the presence of a time-variable error source of order SETI with Gaia 7

4 mas. Future Gaia data releases (DR2 still treats stars One important assumption in the proposed method is as single objects, while DR3 lifts this assumption) will that a Dyson sphere would block a fraction fcov of the likely shed light on this question2. stellar light at optical/near-IR wavelengths, but other- wise not distort the shape of the spectrum in this wave- 5. DISCUSSION length range. A violation of this assumption (for in- While the number of stars for which errors on spec- stance through light reflected off Dyson-sphere elements trophotometric and parallax distances are sufficiently into our line of sight, or some non-thermal emission pro- small to allow Dyson-sphere candidates to be reliably cess) would affect the model fit in unpredictable ways, identified remains relatively small at the current time and could, for instance, make bona-fide Dyson spheres (∼ 104 objects), the sample is expected to increase sig- in a sample of stars end up in the part of the distribution nificantly in the coming years. Gaia DR2 will provide with very poor spectroscopic fits. Because such objects parallaxes for 1.3 billion objects (compared to 2 mil- would supposedly be assigned very large distance errors, lion for DR1) which will immediately increase the over- our current approach would bias against them, since we lap with ground-based surveys like RAVE and GALAH have argued that high-fcov candidates can only be sin- and thereby boost the number of objects for which the gled out (using the technique outlined in this paper) technique can be applied, albeit likely only by factors from the subset of objects for which distance errors are of a few. Gaia DR2 will also feature luminosity and small. ≈ × 7 extinction estimates for 7.7 10 objects. However, 5.2. Grey dust since the stellar parameters will be based on photometry and not on spectroscopy, strong degeneracies between One effect that potentially could give rise to spurious effective temperature and dust reddening are expected, high-fcov Dyson-sphere candidates is grey dust, i.e. line- and it is not clear that this will increase the number of of-sight material that obscures starlight without causing useful objects for our purposes. By Gaia DR33, these any significant reddening of the spectrum. Any dust degeneracies will be significantly reduced through the component with large, micrometer-sized grains would use of data from the Gaia Radial Velocity Spectrometer effectively be grey at optical/near-IR wavelengths, and (RVS) data, and based on the projected performance could – at sufficiently high optical depth – reduce the (Recio-Blanco et al. 2016; Dafonte et al. 2016), we esti- apparent brightness of a star in a way that would be mate that ∼ 106 stars will have sufficiently small spec- very similar to the expected signatures of a Dyson trophotometric distance errors to be included in searches sphere at these wavelengths. Such large dust grains like ours. If we assume the same errors as in the cuts have been seen forming in the circumstellar mate- used in the current paper, a sample of ∼ 106 objects rial around supernovae (Gall et al. 2014; Nielsen et al. 3 2018) and have also been reported in the interstellar could give ∼ 10 spurious detections at fcov > 0.7. How- ever, multi-epoch photometry (Gaia and Pan-starrs) medium (e.g. Wang et al. 2015). High levels of grey combined with archival multi-wavelength data can likely extinction have also been claimed along certain sight prune the sample further, so that only a minor fraction lines (Gorbikov & Brosch 2010), occasionally amount- of these need to be targeted by spectroscopic observ- ing to several magnitudes of optical extinction (e.g. ing campaigns. Additionally, Gaia DR3 will identify Sk´orzy´nski et al. 2003; Krelowski et al. 2016). However, non-single stars, thereby making it easier to weed out it has been argued that grey extinction correlates with binaries like TYC 6111-1162-1. the strengths of certain diffuse interstellar bands in the We can, however, foresee a couple of complications for optical region (e.g. Sk´orzy´nski et al. 2003). High signal- our Dyson-sphere search that are likely to linger despite to-noise spectroscopy of the spectral regions featuring upcoming improvements in Gaia data, which we discuss these bands may therefore be a suitable way to further in the following. scrutinize potential Dyson-sphere candidates. 6. SUMMARY 5.1. The grey-absorber assumption We have presented a method that, through a com- parison of trigonometric and spectrophotometric paral- laxes/distances, makes it possible to single out candi- 2 Note added in proof: Gaia DR2, which was released while this paper was in review, confirms that the Gaia DR1 parallax is flawed dates for nearly complete Dyson spheres, even in cases and claims a significantly smaller parallax of πGaia = 5.75 ± 0.17. where the Dyson-sphere IR-excess signature is much However, this DR2 parallax is still based on the assumption of a single-star solution. We therefore expect the agreement with smaller than usually assumed. The Gaia mission will our spectrophotometric estimate to become even better with Gaia eventually provide both of the required distances (with DR3. highly variable accuracy) and we estimate that Gaia 3 scheduled for late 2020 DR 3 (currently scheduled for late 2020) should allow 8 Zackrisson et al. the technique to be applied to samples of ∼106 stars. A.J.K. acknowledges long-term support by the Swedish A small-scale version of this search can also be imple- National Space Board for work as part of the Gaia Data mented at the current time by combining Gaia paral- Processing and Analysis Consortium (DPAC). This laxes with spectrophotometric distances from a ground- work is based on observations made with the Nordic based survey. By combining parallax data from Gaia Optical Telescope, operated by the Nordic Optical DR1 with spectrophotometric distances from the RAVE Telescope Scientific Association at the Observatorio survey, we find that ≈ 8000 stars (≈ 4% of the ob- del Roque de los Muchachos, La Palma, Spain, of the jects in common between both surveys) have sufficient Instituto de Astrofisica de Canarias. data quality to allow potential Dyson-sphere candidates We would like to thank Alex Ji, Anna Frebel, Peter de to be efficiently singled out. One of the candidates Cat, Sophie van Eck and Alain Jorissen for providing identified this way is TYC-6111-1162-1, an F-type star additional radial-velocity measurements. with an apparent factor-of-2 distance discrepancy (as This work has made use of data from would be expected for a Dyson sphere with covering the European Space Agency (ESA) mission fraction fcov ≈ 0.75), but no detectable mid-IR ex- Gaia (https://www.cosmos.esa.int/gaia), cess. Through new observations of TYC-6111-1162-1 processed by the Gaia DPAC with the NOT/FIES high-resolution spectrograph, we (https://www.cosmos.esa.int/web/gaia/dpac/consortium). largely confirm the stellar parameters found by RAVE Funding for the DPAC has been provided by national (with the exception of the metallicity), i.e. the star re- institutions, in particular the institutions participating mains a marked outlier. Additional observations with in the Gaia Multilateral Agreement. the Magellan/MIKE high-resolution spectrograph reveal Funding for RAVE has been provided by: the Aus- significant radial-velocity differences (≈ 20km/s) with tralian Astronomical Observatory; the Leibniz-Institut respect to the RAVE and NOT/FIES observations. This fuer Astrophysik Potsdam (AIP); the Australian Na- indicates that TYC-6111-1162-1 is part of a binary sys- tional University; the Australian Research Council; tem with an unseen companion. We find it plausible that the French National Research Agency; the German the astrometric wobble resulting from this fact has inter- Research Foundation (SPP 1177 and SFB 881); the fered with the Gaia DR1 astrometric solution and that, European Research Council (ERC-StG 240271 Galac- for this system, the astrometric parallax has been over- tica); the Istituto Nazionale di Astrofisica at Padova; estimated. Taken at face value (given the poor sampling The Johns Hopkins University; the National Science of the radial-velocity curve, assuming circular orbits and Foundation of the USA (AST-0908326); the W. M. Keck ignoring unknown projection effects), a white dwarf of foundation; the Macquarie University; the Netherlands ≈1 M⊙ is needed to double the parallax. Later Gaia data Research School for Astronomy; the Natural Sciences releases may shed further clues on this scenario. Overall, and Engineering Research Council of Canada; the weeding out this type of interloper should become easier Slovenian Research Agency; the Swiss National Science with Gaia DR3 which will feature astrometric solutions Foundation; the Science & Technology Facilities Council for binary stars. of the UK; Opticon; Strasbourg Observatory; and the Universities of Groningen, Heidelberg and Sydney. The ACKNOWLEDGMENTS RAVE web site is at https://www.rave-survey.org. E.Z. acknowledges funding from the Magnus Bergvall foundation and Nordenskj¨oldska Swendeborgsfonden.

REFERENCES

Alonso, A., Arribas, S., & Martinez-Roger, C. 1995, A&A, 297, Choi, J., Dotter, A., Conroy, C., et al. 2016, ApJ, 823, 102 197 Dafonte, C., Fustes, D., Manteiga, M., et al. 2016, A&A, 594, Annis, J. 1999, Journal of the British Interplanetary Society, 52, A68 33 Dyson, F. J. 1960, Science, 131, 1667 Armstrong, S., & Sandberg, A. 2013, Acta Astronautica, 89, 1 Forgan, D. H. 2013, Journal of the British Interplanetary Arnold, L. F. A. 2005, ApJ, 627, 534 Society, 66, 144 Badescu, V., & Cathcart, R. B. 2000, Journal of the British Fuhrmann, K. 2004, Astronomische Nachrichten, 325, 3 Interplanetary Society, 53, 297 Gaia Collaboration, Prusti, T., de Bruijne, J. H. J., et al. 2016, Bailer-Jones, C. A. L. 2015, PASP, 127, 994 A&A, 595, A1 Binney, J., Burnett, B., Kordopatis, G., et al. 2014, MNRAS, Gall, C., Hjorth, J., Watson, D., et al. 2014, Nature, 511, 326 437, 351 Garrett, M. A. 2015, A&A, 581, L5 Bressan, A., Marigo, P., Girardi, L., et al. 2012, MNRAS, 427, Gorbikov, E., & Brosch, N. 2010, MNRAS, 401, 231 127 Griffith, R. L., Wright, J. T., Maldonado, J., et al. 2015, ApJS, Carrigan, Jr., R. A. 2009, ApJ, 698, 2075 217, 25 SETI with Gaia 9

Jugaku, J., & Nishimura, S. 2004, in IAU Symposium, Vol. 213, Slysh, V. I. 1985, in IAU Symposium, Vol. 112, The Search for Bioastronomy 2002: Life Among the Stars, ed. R. Norris & Extraterrestrial Life: Recent Developments, ed. M. D. F. Stootman, 437 Papagiannis, 315–319 Kardashev, N. S. 1964, Soviet Ast., 8, 217 Stempels, E., & Telting, J. 2017, FIEStool: Automated data Kre lowski, J., Galazutdinov, G. A., Strobel, A., & Mulas, G. reduction for FIber-fed Echelle Spectrograph (FIES), 2016, AcA, 66, 469 Astrophysics Source Code Library, , , ascl:1708.009 Kunder, A., Kordopatis, G., Steinmetz, M., et al. 2017, AJ, 153, Suffern, K. G. 1977, Proceedings of the Astronomical Society of 75 Australia, 3 Luri, X., Palmer, M., Arenou, F., et al. 2014, A&A, 566, A119 Telting, J. H., Avila, G., Buchhave, L., et al. 2014, Martell, S. L., Sharma, S., Buder, S., et al. 2017, MNRAS, 465, Astronomische Nachrichten, 335, 41 3203 Timofeev, M. Y., Kardashev, N. S., & Promyslov, V. G. 2000, McMillan, P. J., Kordopatis, G., Kunder, A., et al. 2017, ArXiv Acta Astronautica, 46, 655 e-prints, arXiv:1707.04554 Wang, S., Li, A., & Jiang, B. W. 2015, ApJ, 811, 38 Nielsen, A.-S. B., Hjorth, J., & Gall, C. 2018, A&A, 611, A67 Wright, J. T., Cartier, K. M. S., Zhao, M., Jontof-Hutter, D., & Olson, S. J. 2017, International Journal of Astrobiology, 16, 176 Ford, E. B. 2016, ApJ, 816, 17 Osmanov, Z. 2016, International Journal of Astrobiology, 15, 127 Wright, J. T., Griffith, R. L., Sigurdsson, S., Povich, M. S., & Recio-Blanco, A., de Laverny, P., Allende Prieto, C., et al. 2016, Mullan, B. 2014a, ApJ, 792, 27 A&A, 585, A93 Wright, J. T., Mullan, B., Sigurdsson, S., & Povich, M. S. 2014b, Sandberg, A. 1999, Journal of Evolution and Technology, 5, 1 ApJ, 792, 26 Shkadov, L. M., ed. 1987, Possibility of controlling solar system Zackrisson, E., Calissendorff, P., Asadi, S., & Nyholm, A. 2015, motion in the Galaxy ApJ, 810, 23 Sk´orzy´nski, W., Strobel, A., & Galazutdinov, G. A. 2003, A&A, 408, 297