The Sun-like star Kepler-160 (KOI-456) has been known to host two transiting planets, Kepler-160 b and c, of which planet c shows substantial transit-timing variations (TTVs). We used the archival Kepler photometry of Kepler-160 to search for additional transiting planets using a combination of our Wotan detrending algorithm and our transit least-squares (TLS) detection algorithm. We also used the Mercury N-body gravity code to study the orbital dynamics of the system. First, we recovered the known transit series of planets Kepler-160 b and c. Then we found a new transiting candidate with a radius of 1.91 (+0.17, -0.14) Earth radii (R_ear), an orbital period of 378.417 (+0.028, -0.025) d, and Earth-like insolation. The vespa software predicts that this signal has an astrophysical false-positive probability of FPP_3 = 1.8e-3 when the multiplicity of the system is taken into account. Kepler vetting diagnostics yield a multiple event statistic of MES = 10.7, which corresponds to an ~85 % reliability against false alarms due to instrumental artifacts such as rolling bands. We are also able to explain the observed TTVs of planet c with the presence of a previously unknown planet. The period and mass of this new planet, however, do not match the period and mass of the new transit candidate. Our Markov chain Monte Carlo simulations of the TTVs of Kepler-160 c can be conclusively explained by a new nontransiting planet with a mass between about 1 and 100 Earth masses and an orbital period between about 7 and 50 d. We conclude that Kepler-160 has at least three planets, one of which is the nontransiting planet Kepler-160 d. The expected stellar radial velocity amplitude caused by this new planet ranges between about 1 and 20 m/s. We also find the super-Earth-sized transiting planet candidate KOI-456.04 in the habitable zone of this system, which could be the fourth planet.

It has recently been suggested in this journal by Benford (2019) that "Lurkers" in the form of interstellar exploration probes could be present in the solar system. Similarly, extraterrestrial intelligence could send long-lived probes to many other stellar systems, to report back science and surveillance. If probes and planets with technological species exist in more than a handful of systems in our galaxy, it is beneficial to use a coordinated communication scheme. Due to the inverse square law, data rates decrease strongly for direct connections over long distances. The network bandwidth could be increased by orders of magnitude if repeater stations (nodes) are used in an optimized fashion. This introduction to a series of papers makes the assumptions of the communication scheme explicit. Subsequent papers will discuss technical aspects such as transmitters, repeaters, wavelengths, and power levels. The overall purpose is to gain insight into the physical characteristics of an interstellar communication network, allowing us to describe the most likely sizes and locations of nodes and probes.

The detection of transiting exoplanets in time-series photometry requires the removal or modeling of instrumental and stellar noise. While instrumental systematics can be reduced using methods such as pixel level decorrelation, removing stellar trends while preserving transit signals proves challenging. Due to vast archives of light curves from recent transit surveys, there is a strong need for accurate automatic detrending, without human intervention. A large variety of detrending algorithms are in active use, but their comparative performance for transit discovery is unexplored. We benchmark all commonly used detrending methods against hundreds of Kepler, K2, and TESS planets, selected to represent the most difficult cases for systems with small planet-to-star radius ratios. The full parameter range is explored for each method to determine the best choices for planet discovery. We conclude that the ideal method is a time-windowed slider with an iterative robust location estimator based on Tukey's biweight. This method recovers 99% and 94% of the shallowest Kepler and K2 planets, respectively. We include an additional analysis for young stars with extreme variability and conclude they are best treated using a spline-based method with a robust Huber estimator. All stellar detrending methods explored are available for public use in wotan, an open-source Python package on GitHub (see https://github.com/hippke/wotan).

The extended Kepler mission (K2) has revealed more than 500 transiting planets in roughly 500,000 stellar light curves. All of these were found either with the box least-squares algorithm or by visual inspection. Here we use our new transit least-squares (TLS) algorithm to search for additional planets around all K2 stars that are currently known to host at least one planet. We discover and statistically validate 17 new planets with radii ranging from about 0.7 Earth radii to roughly 2.2 Earth radii and a median radius of 1.18 Earth radii. EPIC201497682.03, with a radius of 0.692 (-0.048, +0.059) Earth radii, is the second smallest planet ever discovered with K2. The transit signatures of these 17 planets are typically 200 ppm deep (ranging from 100 ppm to 2000 ppm), and their orbital periods extend from about 0.7 d to 34 d with a median value of about 4 d. Fourteen of these 17 systems only had one known planet before, and they now join the growing number of multi-planet systems. Most stars in our sample have subsolar masses and radii. The small planetary radii in our sample are a direct result of the higher signal detection efficiency that TLS has compared to box-fitting algorithms in the shallow-transit regime. Our findings help in populating the period-radius diagram with small planets. Our discovery rate of about 3.7 % within the group of previously known K2 systems suggests that TLS can find over 100 additional Earth-sized planets in the data of the Kepler primary mission.

We apply, for the first time, the Transit Least Squares (TLS) algorithm to search for new transiting exoplanets. TLS is a successor to the Box Least Squares (BLS) algorithm, which has served as a standard tool for the detection of periodic transits. In this proof-of-concept paper, we demonstrate how TLS finds small planets that have previously been missed. We showcase TLS' capabilities using the K2 EVEREST-detrended light curve of the star K2-32 (EPIC205071984) that was known to have three transiting planets. TLS detects these known Neptune-sized planets K2-32b, d, and c in an iterative search and finds an additional transit signal with a high signal detection efficiency (SDE_TLS) of 26.1 at a period of 4.34882 (-0.00075, +0.00069) d. We show that this signal remains detectable (SDE_TLS = 13.2) with TLS in the K2SFF light curve of K2-32, which includes a less optimal detrending of the systematic trends. The signal is below common detection thresholds, however, if searched with BLS in the K2SFF light curve (SDE_BLS = 8.9) as in previous searches. Markov Chain Monte Carlo sampling shows that the radius of this candidate is 1.01 (-0.09, +0.10) Earth radii. We analyze its phase-folded transit light curve using the vespa software and calculate a false positive probability FPP = 3.1e-3, formally validating K2-32e as a planet. Taking into account the multiplicity boost of the system, FPP < 3.1e-4. K2-32 now hosts at least four planets that are very close to a 1:2:5:7 mean motion resonance chain. The offset of the orbital periods of K2-32e and b from a 1:2 mean motion resonance is in very good agreement with the sample of transiting multi-planet systems from Kepler, lending further credence to the planetary nature of K2-32e. We expect that TLS can find many more transits of Earth-sized and smaller planets in the Kepler data that have hitherto remained undetected with BLS and similar algorithms.

We present a new method to detect planetary transits from time-series photometry, the Transit Least Squares (TLS) algorithm. TLS searches for transit-like features while taking the stellar limb darkening and planetary ingress and egress into account. We have optimized TLS for both signal detection efficiency (SDE) of small planets and computational speed. TLS analyses the entire, unbinned phase-folded light curve. We compensate for the higher computational load by (i.) using algorithms like "Mergesort" (for the trial orbital phases) and by (ii.) restricting the trial transit durations to a smaller range that encompasses all known planets, and using stellar density priors where available. A typical K2 light curve, including 80d of observations at a cadence of 30min, can be searched with TLS in ~10s real time on a standard laptop computer, as fast as the widely used Box Least Squares (BLS) algorithm. We perform a transit injection-retrieval experiment of Earth-sized planets around sun-like stars using synthetic light curves with 110ppm white noise per 30min cadence, corresponding to a photometrically quiet KP=12 star observed with Kepler. We determine the SDE thresholds for both BLS and TLS to reach a false positive rate of 1% to be SDE~7 in both cases. The resulting true positive (or recovery) rates are ~93% for TLS and ~76% for BLS, implying more reliable detections with TLS. We also test TLS with the K2 light curve of the TRAPPIST-1 system and find six of seven Earth-sized planets using an iterative search for increasingly lower signal detection efficiency, the phase-folded transit of the seventh planet being affected by a stellar flare. TLS is more reliable than BLS in finding any kind of transiting planet but it is particularly suited for the detection of small planets in long time series from Kepler, TESS, and PLATO. We make our Python implementation of TLS publicly available.

PDS 110 is a young disk-hosting star in the Orion OB1A association. Two dimming events of similar depth and duration were seen in 2008 (WASP) and 2011 (KELT), consistent with an object in a closed periodic orbit. In this paper we present data from a ground-based observing campaign designed to measure the star both photometrically and spectroscopically during the time of predicted eclipse in September 2017. Despite high-quality photometry, the predicted eclipse did not occur, although coherent structure is present suggesting variable amounts of stellar flux or dust obscuration. We also searched for RV oscillations caused by any hypothetical companion and can rule out close binaries to 0.1 $M_\odot$. A search of Sonneberg plate archive data also enabled us to extend the photometric baseline of this star back more than 50 years, and similarly does not re-detect any deep eclipses. Taken together, they suggest that the eclipses seen in WASP and KELT photometry were due to aperiodic events. It would seem that PDS 110 undergoes stochastic dimmings that are shallower and shorter-duration than those of UX Ori variables, but may have a similar mechanism.

In recent publications, Borra (2013); Borra & Trottier (2016); Borra (2017) claimed the discovery of ultra-short ($10^{-12}\,$s) optical pulses originating from stars and galaxies, asserted to be sent by extraterrestrial intelligence. I show that these signals are not astrophysical or instrumental in nature, but originate from the non-random spacings of spectral absorption lines. They can be shown to arise in their clearest form in synthetic solar spectra, as these do not suffer from noise.

Additional pulsation modes have been discovered in many Cepheids, RR Lyrae, and other variable stars. Fourier transforms are used to find, fit and subtract the main pulsation period and its harmonics to reveal additional modes. Commonly, for every star, the strongest of these modes is presented in a "Petersen diagram", where the shorter-to-longer period ratio is plotted against the longer period. This diagram discards the information about temporal variations, multiple pulsation modes, and signals which are below some chosen signal to noise threshold. I here present a new tool to map pulsation modes in variable stars, dubbed the "spectral Petersen diagram". It shows all signals, irrespective of their multiplicity or significance. Many (thousands) of light curves can be stacked to improve the signal to noise ratio and reveal unprecedented detail about the pulsation mode space. This tool is useful to constraint the parameters of variable star models.

Transit photometry of the exoplanet candidate Kepler-1625b has recently been interpreted to show hints of a moon. We aim to clarify whether the exomoon-like signal is really caused by a large object in orbit around Kepler-1625b. We explore several detrending procedures, i.e. polynomials and the Cosine Filtering with Autocorrelation Minimization (CoFiAM). We then supply a light curve simulator with the co-planar orbital dynamics of the system and fit the resulting planet-moon transit light curves to the Kepler data. We employ the Bayesian Information Criterion (BIC) to assess whether a single planet or a planet-moon system is a more likely interpretation of the light curve variations. We carry out a blind hare-and-hounds exercise using many noise realizations by injecting simulated transits into different out-of-transit parts of the original Kepler-1625 data: 100 sequences with 3 synthetic transits of a Kepler-1625b-like planet and 100 sequences with 3 synthetic transits of this planet with a Neptune-sized moon. The statistical significance and characteristics of the exomoon-like signal strongly depend on the detrending method, and the data chosen for detrending, and on the treatment of gaps in the light curve. Our injection-retrieval experiment shows evidence for moons in about 10% of those light curves that do not contain an injected moon. Strikingly, many of these false-positive moons resemble the exomoon candidate. We recover up to about half of the injected moons, depending on the detrending method, with radii and orbital distances broadly corresponding to the injected values. A $\Delta$BIC of -4.9 for the CoFiAM-based detrending indicates an exomoon around Kepler-1625b. This solution, however, is only one out of many and we find very different solutions depending on the details of the detrending method. It is worrying that the detrending is key to the interpretation of the data.

We have explored the optimal frequency of interstellar photon communications and benchmarked other particles as information carriers in previous papers of this series. We now compare the latency and bandwidth of sending probes with inscribed matter. Durability requirements such as shields against dust and radiation, as well as data duplication, add negligible weight overhead at velocities <0.2c. Probes may arrive in full, while most of a photon beam is lost to diffraction. Probes can be more energy efficient per bit, and can have higher bandwidth, compared to classical communication, unless a photon receiver is placed in a stellar gravitational lens. The probe's advantage dominates by order of magnitude for long distances (kpc) and low velocities (<0.1c) at the cost of higher latency.

Many rocky exoplanets are heavier and larger than the Earth, and have higher surface gravity. This makes space-flight on these worlds very challenging, because the required fuel mass for a given payload is an exponential function of planetary surface gravity. We find that chemical rockets still allow for escape velocities on Super-Earths up to 10x Earth mass. More massive rocky worlds, if they exist, would require other means to leave the planet, such as nuclear propulsion.

The search for extraterrestrial communication has mainly focused on microwave photons since the 1950s. We compare other high speed information carriers to photons, such as electrons, protons, and neutrinos, gravitational waves, inscribed matter, and artificial megastructures such as occulters. The performance card includes the speed of exchange, information per energy and machine sizes, lensing performance, cost, and complexity. In fast point-to-point communications, photons are superior to other carriers by orders of magnitude. Sending probes with inscribed matter requires less energy, but has higher latency. For isotropic beacons with low data rates, our current technological level is insufficient to determine the best choice. We discuss cases where our initial assumptions do not apply, and describe the required properties of hypothetical particles to win over photons.

Many rocky exoplanets are heavier and larger than Earth, so-called "Super-Earths". Some of these may be habitable, and a few may be inhabited by Super-Earthlings. Due to the higher surface gravity on these worlds, space-flight is much more challenging. We find that chemical rockets still allow for escape velocities on Super-Earths up to 10x Earth mass. Much heavier rocky worlds, if they exist, will require using up most of the planet as chemical fuel for the (one) launch, a rather risky undertaking. We also briefly discuss launching rockets from water worlds, which requires Alien megastructures.

Previous and ongoing searches for extraterrestrial optical and infrared nanosecond laser pulses and narrow line-width continuous emissions have so far returned null results. At the commonly used observation cadence of $\sim 10^{-9}\,$s, sky-integrated starlight is a relevant noise source for large field-of-view surveys. This can be reduced with narrow bandwidth filters, multipixel detectors, or a shorter observation cadence. We examine the limits of short pulses set by the uncertainty principle, interstellar scattering, atmospheric scintillation, refraction, dispersion and receiver technology. We find that optimal laser pulses are time-bandwidth limited Gaussians with a duration of $\Delta t \approx\,10^{-12}\,$s at a wavelength $\lambda_{0}\approx1\,\mu$m, and a spectral width of $\Delta \lambda \approx 1.5\,$nm. Shorter pulses are too strongly affected through Earth's atmosphere. Given certain technological advances, survey speed can be increased by three orders of magnitude when moving from ns to ps pulses. Faster (and/or parallel) signal processing would allow for an all-sky-at-once survey of lasers targeted at Earth.

It has recently been argued from a laser engineering point of view that there are only a few magic colors for optical SETI. These are primarily the Nd:YAG line at 1064 nm and its second harmonic 532.1 nm. Next best choices would be the sum frequency and/or second harmonic generation of Nd:YAG and Nd:YLF laser lines, 393.8 nm (near Fraunhofer CaK), 656.5 nm (H$\alpha$) and 589.1 nm (NaD2). In this paper, we examine the interstellar extinction, atmospheric transparency and scintillation, as well as noise conditions for these laser lines. For strong signals, we find that optical wavelengths are optimal for distances $d\lesssim\,$kpc. Nd:YAG at $\lambda=1{,}064\,$nm is a similarly good choice, within a factor of two, under most conditions and out to $d\lesssim3\,$kpc. For weaker transmitters, where the signal-to-noise ratio with respect to the blended host star is relevant, the optimal wavelength depends on the background source, such as the stellar type. Fraunhofer spectral lines, while providing lower stellar background noise, are irrelevant in most use cases, as they are overpowered by other factors. Laser-pushed spaceflight concepts, such as "Breakthrough Starshot", would produce brighter and tighter beams than ever assumed for OSETI. Such beamers would appear as naked eye stars out to kpc distances. If laser physics has already matured and converged on the most efficient technology, the laser line of choice for a given scenario (e.g., Nd:YAG for strong signals) can be observed with a narrow filter to dramatically reduce background noise, allowing for large field-of-view observations in fast surveys.

A complex message from space may require the use of computers to display, analyze and understand. Such a message cannot be decontaminated with certainty, and technical risks remain which can pose an existential threat. Complex messages would need to be destroyed in the risk averse case.

We combine and calibrate publicly available data for Boyajian's star including photometry from ASAS (SN, V, I), Kepler, Gaia, SuperWASP, and citizen scientist observations (AAVSO, HAO and Burke-Gaffney). Precise (mmag) photometry covers the years 2006-2017. We show that the year-long flux variations with an amplitude of ~4% can not be explained with cyclical symmetric or asymmetric models with periods shorter than ten years. If the dips are transits, their period must exceed ten years, or their structure must evolve significantly during each 4-year long cycle.

A photon can encode several bits of information based on an alphabet of its time of arrival, energy, and polarization. Heisenberg's uncertainty principle places a limit on measuring pairs of physical properties of a particle, limiting the maximal information efficiency to <59 bits per photon in practice, and <171 bits per photon at Planck energy, at a data rate of one photon per second.

We have previously argued that targeted interstellar communication has a physical optimum at narrowband X-ray wavelengths $\lambda\approx1\,$nm, limited by the surface roughness of focusing devices at the atomic level (arXiv:1711.05761). We search 24,247 archival X-ray spectra (of 6,454 unique objects) for such features and present 19 sources with monochromatic signals. Close examination reveals that these are most likely of natural origin. The ratio of artificial to natural sources must be <0.01%. This first limit can be improved in future X-ray surveys.

How many bits of information can a single photon carry? Intuition says "one", but this is incorrect. With an alphabet based on the photon's time of arrival, energy, and polarization, several bits can be encoded. In this introduction to photon information efficiency, we explain how to calculate the maximum number of bits per photon depending on the number of encoding modes, noise, and losses.

Recent technological advances could make interstellar travel possible, using ultra-lightweight sails pushed by lasers or solar photon pressure, at speeds of a few percent the speed of light. Obtaining remote observational data from such probes is not trivial because of their minimal instrumentation (gram scale) and large distances (pc). We derive the optimal communication scheme to maximize the data rate between a remote probe and home-base. he framework requires coronagraphic suppression of the stellar background at the level of $10^{-9}$ within a few tenths of an arcsecond of the bright star. Our work includes models for the loss of photons from diffraction, technological limitations, interstellar extinction, and atmospheric transmission. Major noise sources are atmospheric, zodiacal, stellar and instrumental. We examine the maximum capacity using the "Holevo bound" which gives an upper limit to the amount of information (bits) that can be encoded through a quantum state (photons), which is a few bits per photon for optimistic signal and noise levels. This allows for data rates of order bits per second per Watt from a transmitter of size 1 m at a distance of $\alpha\,$Centauri (1.3 pc) to an earth-based large receiving telescope (E-ELT, 39 m). The optimal wavelength for this distance is 300 nm (space-based receiver) to 400 nm (earth-based) and increases with distance, due to extinction, to a maximum of $\approx3\,\mu$m to the center of the galaxy at 8 kpc.

The optimal frequency for interstellar communication, using "Earth 2017" technology, was derived in papers I and II of this series (arXiv:1706.03795, arXiv:1706.05570). The framework included models for the loss of photons from diffraction (free space), interstellar extinction, and atmospheric transmission. A major limit of current technology is the focusing of wavelengths $\lambda<300\,$nm (UV). When this technological constraint is dropped, a physical bound is found at $\lambda\approx1\,$nm ($E\approx\,$keV) for distances out to kpc. While shorter wavelengths may produce tighter beams and thus higher data rates, the physical limit comes from surface roughness of focusing devices at the atomic level. This limit can be surpassed by beam-forming with electromagnetic fields, e.g. using a free electron laser, but such methods are not energetically competitive. Current lasers are not yet cost efficient at nm wavelengths, with a gap of two orders of magnitude, but future technological progress may converge on the physical optimum. We recommend expanding SETI efforts towards targeted (at us) monochromatic (or narrow band) X-ray emission at 0.5-2 keV energies.

Upon entering a star system, light sails are subject to both gravitational forces and radiation pressure, and can use both in concert to modify their trajectory. Moreover, stars possess significant magnetic fields, and if the sail is in any way charged, it will feel the Lorentz force also. We investigate the dynamics of so-called "photogravimagnetic assists" of sailcraft around $\alpha$ Centauri A, a potential first destination en route to Proxima Centauri (the goal of the Breakthrough Starshot program). We find that a 10m$^2$ sail with a charge-to-mass-ratio of around 10 $\mu$C/g or higher will need to take account of magnetic field effects during orbital maneouvres. The magnetic field can provide an extra source of deceleration and deflection, and allow capture onto closer orbits around a target star. However, flipping the sign of the sailcraft's charge can radically change resulting trajectories, resulting in complex loop-de-loops around magnetic field lines and essentially random ejection from the star system. Even on well-behaved trajectories, the field can generate off-axis deflections at $\alpha$ Centauri that, while minor, can result in very poor targetting of the final destination (Proxima) post-assist. Fortunately for Breakthrough Starshot, nanosails are less prone to charging en route than their heavier counterparts, but can still accrue relatively high charge at both the origin and destination, when travelling at low speeds. Photogravimagnetic assists are highly non-trivial, and require careful course correction to mitigate against unwanted changes in trajectory.

We have shown in paper I of this series (arXiv:1706.03795) that interstellar communication to nearby (pc) stars is possible at data rates of bits per second per Watt between a 1 m sized probe and a large receiving telescope (E-ELT, 39 m), when optimizing all parameters such as frequency at 300-400 nm. We now apply our framework of interstellar extinction and quantum state calculations for photon encoding to the solar gravitational lens (SGL), which enlarges the aperture (and thus the photon flux) of the receiving telescope by a factor of $>10^9$. For the first time, we show that the use of the SGL for communication purposes is possible. This was previously unclear because the Einstein ring is placed inside the solar coronal noise, and contributing factors are difficult to determine. We calculate point-spread functions, aperture sizes, heliocentric distance, and optimum communication frequency. The best wavelength for nearby (<100 pc) interstellar communication is limited by current technology to the UV and optical band. To suppress coronal noise, an advanced coronograph is required, alternatively an occulter could be used which would require a second spacecraft in formation flight 78 km from the receiver, and ~10 m in size. Data rates scale approximately linear with the SGL telescope size and with heliocentric distance. Achievable (receiving) data rates from Alpha Cen are 1-10 Mbits per second per Watt for a pair of meter-sized telescopes, an improvement of $10^6$ compared to using the same receiving telescope without the SGL. A 1 m telescope in the SGL can receive data at rates comparable to a km-class "normal" telescope.

New means of interstellar travel are now being considered by various research teams, assuming lightweight spaceships to be accelerated via either laser or solar radiation to a significant fraction of the speed of light (c). We recently showed that gravitational assists can be combined with the stellar photon pressure to decelerate an incoming lightsail from Earth and fling it around a star or bring it to rest. Here, we demonstrate that photogravitational assists are more effective when the star is used as a bumper (i.e. the sail passes "in front of" the star) rather than as a catapult (i.e. the sail passes "behind" or "around" the star). This increases the maximum deceleration at $\alpha$ Cen A and B and reduces the travel time of a nominal graphene-class sail (mass-to-surface ratio 8.6e-4 gram m$^{-2}$) from 95 to 75 yr. The maximum possible velocity reduction upon arrival depends on the required deflection angle from $\alpha$ Cen A to B and therefore on the binary's orbital phase. Here, we calculate the variation of the minimum travel times from Earth into a bound orbit around Proxima for the next 300 yr and then extend our calculations to roughly 22,000 stars within about 300 ly. Although $\alpha$ Cen is the most nearby star system, we find that Sirius A offers the shortest possible travel times into a bound orbit: 69 yr assuming 12.5% c can be obtained at departure from the solar system. Sirius A thus offers the opportunity of flyby exploration plus deceleration into a bound orbit of the companion white dwarf after relatively short times of interstellar travel.

We revisit the nature of large dips in flux from extinction by dusty circumstellar material that is observed by Kepler for many young stars in the Upper Sco and $\rho$ Oph star formation regions. These young, low-mass "dipper" stars are known to have low accretion rates and primarily host moderately evolved dusty circumstellar disks. Young low mass stars often exhibit rotating star spots that cause quasi-periodic photometric variations. We found no evidence for periods associated with the dips that are different from the star spot rotation period in spectrograms constructed from the light curves. The material causing the dips in most of these light curves must be approximately corotating with the star. We find that disk temperatures computed at the disk corotation radius are cool enough that dust should not sublime. Crude estimates for stellar magnetic field strengths and accretion rates are consistent with magnetospheric truncation near the corotation radius. Magnetospheric truncation models can explain why the dips are associated with material near corotation and how dusty material is lifted out of the midplane to obscure the star which would account for the large fraction of young low mass stars that are dippers. We propose that variations in disk orientation angle, stellar magnetic field dipole tilt axis, and disk accretion rate are underlying parameters accounting for differences in the dipper light curves.

The F3 main sequence star KIC 8462852 (Boyajian's Star) showed deep (up to 20%) day-long brightness dips of unknown cause during the 4 years of the Kepler mission. A 0.164 mag (16%) dimming between 1890 and 1990 was claimed, based on the analysis of photographic plates from the Harvard Observatory. We have gathered an independent set of historic plates from Sonneberg Observatory, Germany, covering the years 1934-1995. With 861 magnitudes in B, and 397 magnitudes in V, we find the star to be of constant brightness within 0.03 mag per century (3%). Consistent outcomes are found using by-eye estimates of the best 119 plates. Results are supported by data from Sternberg Observatory, Moscow, which show the star as constant between 1895 and 1995. The previously claimed century-long dimming is inconsistent with our results at the 5-sigma level, however the recently reported modest dimming of 3% in the Kepler data is not inconsistent with our data. We find no periodicities or shorter trends within our limits of 5% per 5-year bin, but note a possible dimming event on 24 Oct 1978.

At a distance of about 4.22 lightyears, it would take about 100,000 years for humans to visit our closest stellar neighbor Proxima Centauri using modern chemical thrusters. New technologies are now being developed that involve high-power lasers firing at 1 gram solar sails in near-Earth orbits, accelerating them to 20% the speed of light (c) within minutes. Although such an interstellar probe could reach Proxima 20 years after launch, without propellant to slow it down it would traverse the system within hours. Here we demonstrate how the stellar photon pressures of the stellar triple $\alpha$ Cen A, B, and C (Proxima) can be used together with gravity assists to decelerate incoming solar sails from Earth. The maximum injection speed at $\alpha$ Cen A to park a sail with a mass-to-surface ratio ($\sigma$) similar to graphene (7.6e-4 gram/m$^{-2}$) in orbit around Proxima is about 13,800 km/s (4.6% c), implying travel times from Earth to $\alpha$ Cen A and B of about 95 years and another 46 years (with a residual velocity of 1280 km/s) to Proxima. The size of such a low-$\sigma$ sail required to carry a payload of 10 grams is about 10$^5$ m$^2$ = (316 m)$^2$. Such a sail could use solar photons instead of an expensive laser system to gain interstellar velocities at departure. Photogravitational assists allow visits of three stellar systems and an Earth-sized potentially habitable planet in one shot, promising extremely high scientific yields.

Distance estimates from Gaia parallax and expected luminosities are compared for KIC 8462852. Gaia DR1 yields a parallax of $2.55\pm0.31$mas, that is a distance of $391.4\substack{+53.6 \\ -42.0}$pc, or $391.4\substack{+122.1 \\ -75.2}$pc including systematic uncertainty. The distance estimate based on the absolute magnitude of an F3V star and measured reddening is $\sim454\pm35$pc. Both estimates agree within $<1\sigma$, which only excludes some of the most extreme theorized scenarios for KIC 8462852. Future Gaia data releases will determine the distance to within 1% and thus allow for the determination a precise absolute luminosity.

We present new ways to identify single and multiple moons around extrasolar planets using planetary transit timing variations (TTVs) and transit duration variations (TDVs). For planets with one moon, measurements from successive transits exhibit a hitherto undescribed pattern in the TTV-TDV diagram, originating from the stroboscopic sampling of the planet's orbit around the planet-moon barycenter. This pattern is fully determined and analytically predictable after three consecutive transits. The more measurements become available, the more the TTV-TDV diagram approaches an ellipse. For planets with multiple moons in orbital mean motion resonance (MMR), like the Galilean moons, the pattern is much more complex and addressed numerically in this report. Exomoons in MMR can also form closed, predictable TTV-TDV figures if the drift of the moons' pericenters is sufficiently slow. We find that MMR exomoons produce loops in the TTV-TDV diagram and that the number of these loops is equal to the order of the MMR, or the largest integer in the MMR ratio. We use a Bayesian model and Monte Carlo simulations to test the discoverability of exomoons using TTV-TDV diagrams with current and near-future technology. In a blind test, two of us (BP, DA) successfully retrieved a large moon from simulated TTV-TDV by co-authors MH and RH, which resembled data from a known Kepler planet candidate. Single exomoons with a 10% moon-to-planet mass ratio, like to Pluto-Charon binary, can be detectable in the archival data of the Kepler primary mission. Multi-exomoon systems, however, require either larger telescopes or brighter target stars. Complementary detection methods invoking a moon's own photometric transit or its orbital sampling effect can be used for validation or falsification. A combination of TESS, CHEOPS, and PLATO data would offer a compelling opportunity for exomoon discoveries around bright stars.

The century-long photometric record of the DASCH project provides a unique window into the variability of stars normally considered to be photometrically inactive. In this paper, we look for long-term trends in the brightness of F stars, with particular attention to KIC 8462852,an F3 main sequence star that has been identified as significant short-term variability according to Kepler observations. Although a simple search for variability suggests long-term dimming of a number of F stars, we find that such trends are artifacts of the 'Menzel Gap' in the DASCH data. That includes the behavior of KIC 8462852, which we believe is consistent with constant flux over the full duration of observations. We do, however, present a selection of F stars that do have significant photometric trends, even after systematics are taken into account.

We present a statistical analysis of the accuracy of the digitized magnitudes of photometric plates on the time scale of decades. In our examination of archival Johnson B photometry from the Harvard DASCH archive, we find a median RMS scatter of lightcurves of order 0.15mag over the range B~9-17 for all calibrations. Slight underlying systematics (trends or flux discontinuities) are on a level of ~<=0.2mag per century (1889-1990) for the majority of constant stars. These historic data can be unambiguously used for processes that happen on scales of magnitudes, and need to be carefully examined in cases approaching the noise floor. The characterization of these limits in photometric stability may guide future studies in their use of plate archives. We explain these limitations for the example case of KIC8462852, which has been claimed to dim by 0.16mag per century, and show that this trend cannot be considered as significant.

The orbital sampling effect (OSE) appears in phase-folded transit light curves of extrasolar planets with moons. Analytical OSE models have hitherto neglected stellar limb darkening and non-zero transit impact parameters and assumed that the moon is on a circular, co-planar orbit around the planet. Here, we present an analytical OSE model for eccentric moon orbits, which we implement in a numerical simulator with stellar limb darkening that allows for arbitrary transit impact parameters. We also describe and publicly release a fully numerical OSE simulator (PyOSE) that can model arbitrary inclinations of the transiting moon orbit. Both our analytical solution for the OSE and PyOSE can be used to search for exomoons in long-term stellar light curves such as those by Kepler and the upcoming PLATO mission. Our updated OSE model offers an independent method for the verification of possible future exomoon claims via transit timing variations and transit duration variations. Photometrically quiet K and M dwarf stars are particularly promising targets for an exomoon discovery using the OSE.

Dramatically improved data from observatories like the CoRoT and Kepler spacecraft have recently facilitated nonlinear time series analysis and phenomenological modeling of variable stars, including the search for strange (aka fractal) or chaotic dynamics. We recently argued [Lindner et al., Phys. Rev. Lett. 114 (2015) 054101] that the Kepler data includes "golden" stars, whose luminosities vary quasiperiodically with two frequencies nearly in the golden ratio, and whose secondary frequencies exhibit power-law scaling with exponent near -1.5, suggesting strange nonchaotic dynamics and singular spectra. Here we use a series of phenomenological models to make plausible the connection between golden stars and fractal spectra. We thereby suggest that at least some features of variable star dynamics reflect universal nonlinear phenomena common to even simple systems.

Trojans are small bodies in planetary Lagrangian points. In our solar system, Jupiter has the largest number of such companions. Their existence is assumed for exoplanetary systems as well, but none has been found so far. We present an analysis by super-stacking $\sim4\times10^4$ Kepler planets with a total of $\sim9\times10^5$ transits, searching for an average trojan transit dip. Our result gives an upper limit to the average Trojan transiting area (per planet) corresponding to one body of radius $<460$km at $2\sigma$ confidence. We find a significant Trojan-like signal in a sub-sample for planets with more (or larger) Trojans for periods $>$60 days. Our tentative results can and should be checked with improved data from future missions like PLATO2.0, and can guide planetary formation theories.

Time-series transit photometry from the Kepler space telescope has allowed for the discovery of thousands of exoplanets. We explore the potential of yet improved future missions such as PLATO 2.0 in detecting solar system analogues. We use real-world solar data and end-to-end simulations to explore the stellar and instrumental noise properties. By injecting and retrieving planets, rings and moons of our own solar system, we show that the discovery of Venus- and Earth-analogues transiting G-dwarfs like our Sun is feasible at high S/N after collecting 6yrs of data, but Mars and Mercury will be difficult to detect due to stellar noise. In the best cases, Saturn's rings and Jupiter's moons will be detectable even in single transit observations. Through the high number (>1 bn) of observed stars by PLATO 2.0, it will become possible to detect thousands of single-transit events by cold gas giants, analogue to our Jupiter, Saturn, Uranus and Neptune. Our own solar system aside, we also show, through signal injection and retrieval, that PLATO 2.0-class photometry will allow for the secure detection of exomoons transiting quiet M-dwarfs. This is the first study analyzing in-depth the potential of future missions, and the ultimate limits of photometry, using realistic case examples.

Despite the discovery of thousands of exoplanets, no exomoons have been detected so far. We test a recently developed method for exomoon search, the orbital sampling effect (OSE), using the full exoplanet photometry from the Kepler Space Telescope. The OSE is applied to phase-folded transits, for which we present a framework to detect false positives, and discuss four candidates which pass several of our tests. Using numerical simulations, we inject exomoon signals into real Kepler data and retrieve them, showing that under favorable conditions, exomoons can be found with Kepler and the OSE method. In addition, we super-stack a large sample of Kepler planets to search for the average exomoon OSE and the accompanying increase in noise, the scatter peak. We find a significant OSE-like signal, which might indicate the presence of moons, for planets with 35d<P<80d, having an average dip per planet of 6+-2ppm, corresponding to a moon radius of 2120 +330 -370km for the average star radius of 1.24R_Sun in this sample.

Fast Radios Bursts (FRBs) show large dispersion measures (DMs), suggesting an extragalactic location. We analyze the DMs of the 11 known FRBs in detail and identify steps as integer multiples of half the lowest DM found, 187.5cm$^{-3}$ pc, so that DMs occur in groups centered at 375, 562, 750, 937, 1125cm$^{-3}$ pc, with errors observed <5%. We estimate the likelhood of a coincidence as 5:10,000. We speculate that this could originate from a Galaxy population of FRBs, with Milky Way DM contribution as model deviations, and an underlying generator process that produces FRBs with DMs in discrete steps. However, we find that FRBs tend to arrive at close to the full integer second, like man-made perytons. If this holds, FRBs would also be man-made. This can be verified, or refuted, with new FRBs to be detected.

The unprecedented light curves of the Kepler space telescope document how the brightness of some stars pulsates at primary and secondary frequencies whose ratios are near the golden mean, the most irrational number. A nonlinear dynamical system driven by an irrational ratio of frequencies generically exhibits a strange but nonchaotic attractor. For Kepler's "golden" stars, we present evidence of the first observation of strange nonchaotic dynamics in nature outside the laboratory. This discovery could aid the classification and detailed modeling of variable stars.

Learned et. al. proposed that a sufficiently advanced extra-terrestrial civilization may tickle Cepheid and RR Lyrae variable stars with a neutrino beam at the right time, thus causing them to trigger early and jogging the otherwise very regular phase of their expansion and contraction. This would turn these stars into beacons to transmit information throughout the galaxy and beyond. The idea is to search for signs of phase modulation (in the regime of short pulse duration) and patterns, which could be indicative of intentional, omnidirectional signaling. We have performed such a search among variable stars using photometric data from the Kepler space telescope. In the RRc Lyrae star KIC 5520878, we have found two such regimes of long and short pulse durations. The sequence of period lengths, expressed as time series data, is strongly auto correlated, with correlation coefficients of prime numbers being significantly higher ($p=99.8$\%). Our analysis of this candidate star shows that the prime number oddity originates from two simultaneous pulsation periods and is likely of natural origin. Simple physical models elucidate the frequency content and asymmetries of the KIC 5520878 light curve. Despite this SETI null result, we encourage testing other archival and future time-series photometry for signs of modulated stars. This can be done as a by-product to the standard analysis, and even partly automated.

We document the presence of a few Cepheid and RR Lyrae variable stars with previously unrecognized characteristics. These stars exhibit the property of a period ratio of main pulsation divided by secondary pulsation P1/P2 very close to sqrt(2). Other stars of these types have period ratios which do not show clustering with a close association and a single remarkable non-harmonic number. Close examination reveals a deviation of multiples of a few times ~0.06% for these stars. This deviation seems to be present in discrete steps on the order of ~0.000390(4), indicating the possible presence of a sort of fine structure in this oscillation.