The nearby ultracool dwarf TRAPPIST-1 possesses several Earth-sized terrestrial planets, three of which have equilibrium temperatures that may support liquid surface water, making it a compelling target for exoplanet characterization. TRAPPIST-1 is an active star with frequent flaring, with implications for the habitability of its planets. Superflares (stellar flares whose energy exceeds 10^33 erg) can completely destroy the atmospheres of a cool star's planets, allowing ultraviolet radiation and high-energy particles to bombard their surfaces. However, ultracool dwarfs emit little ultraviolet flux when quiescent, raising the possibility of frequent flares being necessary for prebiotic chemistry that requires ultraviolet light. We combine Evryscope and Kepler observations to characterize the high-energy flare rate of TRAPPIST-1. The Evryscope is an array of 22 small telescopes imaging the entire Southern sky in g' every two minutes. Evryscope observations, spanning 170 nights over 2 years, complement the 80-day continuous short-cadence K2 observations by sampling TRAPPIST-1's long-term flare activity. We update TRAPPIST-1's superflare rate, finding a cumulative rate of 4.2 (+1.9 -0.2) superflares per year. We calculate the flare rate necessary to deplete ozone in the habitable-zone planets' atmospheres, and find that TRAPPIST-1's flare rate is insufficient to deplete ozone if present on its planets. In addition, we calculate the flare rate needed to provide enough ultraviolet flux to power prebiotic chemistry. We find TRAPPIST-1's flare rate is likely insufficient to catalyze some of the Earthlike chemical pathways thought to lead to RNA synthesis, and flux due to flares in the biologically relevant UV-B band is orders of magnitude less for any TRAPPIST-1 planet than has been experienced by Earth at any time in its history.

We measure rotation periods and sinusoidal amplitudes in Evryscope light curves for 122 two-minute K5-M4 TESS targets selected for strong flaring. The Evryscope array of telescopes has observed all bright nearby stars in the South, producing two-minute cadence light curves since 2016. Long-term, high-cadence observations of rotating flare stars probe the complex relationship between stellar rotation, starspots, and superflares. We detect periods from 0.3487 to 104 d, and observe amplitudes from 0.008 to 0.216 g' mag. We find the Evryscope amplitudes are larger than those in TESS with the effect correlated to stellar mass (p-value=0.01). We compute the Rossby number (Ro), and find our sample selected for flaring has twice as many intermediate rotators (0.04<Ro<0.4) as fast (Ro<0.04) or slow (Ro>0.44) rotators; this may be astrophysical or a result of period-detection sensitivity. We discover 30 fast, 59 intermediate, and 33 slow rotators. We measure a median starspot coverage of 13% of the stellar hemisphere and constrain the minimum magnetic field strength consistent with our flare energies and spot coverage to be 500 G, with later-type stars exhibiting lower values than earlier-types. We observe a possible change in superflare rates at intermediate periods. However, we do not conclusively confirm the increased activity of intermediate rotators seen in previous studies. We split all rotators at Ro~0.2 into Prot<10 d and Prot>10 d bins to confirm short-period rotators exhibit higher superflare rates, larger flare energies, and higher starspot coverage than do long-period rotators, at p-values of 3.2 X 10^-5, 1.0 X 10^-5, and 0.01, respectively.

Camera lenses are increasingly used in wide-field astronomical surveys due to their high performance, wide field-of-view (FOV) unreachable from traditional telescope optics, and modest cost. The machining and assembly tolerances for commercially available optical systems cause a slight misalignment (tilt) between the lens and CCD, resulting in PSF degradation. We have built an automated alignment system (Robotilters) to solve this challenge, optimizing 4 degrees of freedom - 2 tilt axes, a separation axis (the distance between the CCD and lens), and the lens focus (the built-in focus of the lens by turning the lens focusing ring which moves the optical elements relative to one another) in a compact and low-cost package. The Robotilters remove tilt and optimize focus at the sub 10 micron level, are completely automated, take 2 hours to run, and remain stable for multiple years once aligned. The Robotilters were built for the Evryscope telescope (a 780 MPix 22-camera array with an 8150 sq.deg. field of view and continuous 2-minute cadence) designed to detect short timescale events across extremely large sky areas simultaneously. Variance in quality across the image field, especially the corners and edges compared to the center, is a significant challenge in wide-field astronomical surveys like the Evryscope. The individual star PSFs (which typically extend only a few pixels) are highly susceptible to slight increases in optical aberrations in this situation. The Robotilter solution resulted in a limiting magnitude improvement of .5 mag in the center of the image and 1.0 mag in the corners for typical Evryscope cameras, with less distorted and smaller PSFs (half the extent in the corners and edges in many cases). In this paper we describe the Robotilter mechanical and software design, camera alignment results, long term stability, and image improvement.

We search for superflares from 4,068 cool stars in 2+ years of Evryscope photometry, focusing on those with high-cadence data from both Evryscope and TESS. The Evryscope array of small telescopes observed 575 flares from 284 stars, with a median energy of 10^34.0 erg. Since 2016, Evryscope has enabled the detection of rare events from all stars observed by TESS through multi-year, high-cadence continuous observing. We report ~2X the previous largest number of 10^34 erg high-cadence flares from nearby cool stars. We find 8 flares with amplitudes of 3+ g' magnitudes, with the largest reaching 5.6 magnitudes and releasing 10^36.2 erg. We observe a 10^34 erg superflare from TOI-455 (LTT 1445), a mid-M with a rocky planet candidate. We measure the superflare rate per flare-star and quantify the average flaring of active stars as a function of spectral type, including superflare rates, FFDs, and typical flare amplitudes in g'. We confirm superflare morphology is broadly consistent with magnetic re-connection. We estimate starspot coverage necessary to produce superflares, and hypothesize maximum-allowed superflare energies and waiting-times between flares corresponding to 100% coverage of the stellar hemisphere. We observe decreased flaring at high galactic latitudes. We explore the effects of superflares on ozone loss to planetary atmospheres: we observe 1 superflare with sufficient energy to photo-dissociate all ozone in an Earth-like atmosphere in one event. We find 17 stars that may deplete an Earth-like atmosphere via repeated flaring. Of the 1822 stars around which TESS may discover temperate rocky planets, we observe 14.6% +/- 2% emit large flares.

The regions around the celestial poles offer the ability to find and characterize long-term variables from ground-based observatories. We used multi-year Evryscope data to search for high-amplitude (~5% or greater) variable objects among 160,000 bright stars (Mv < 14.5) near the South Celestial Pole. We developed a machine learning based spectral classifier to identify eclipse and transit candidates with M-dwarf or K-dwarf host stars - and potential low-mass secondary stars or gas giant planets. The large amplitude transit signals from low-mass companions of smaller dwarf host stars lessens the photometric precision and systematics removal requirements necessary for detection, and increases the discoveries from long-term observations with modest light curve precision. The Evryscope is a robotic telescope array that observes the Southern sky continuously at 2-minute cadence, searching for stellar variability, transients, transits around exotic stars and other observationally challenging astrophysical variables. In this study, covering all stars 9 < Mv < 14.5, in declinations -75 to -90 deg, we recover 346 known variables and discover 303 new variables, including 168 eclipsing binaries. We characterize the discoveries and provide the amplitudes, periods, and variability type. A 1.7 Jupiter radius planet candidate with a late K-dwarf primary was found and the transit signal was verified with the PROMPT telescope network. Further followup revealed this object to be a likely grazing eclipsing binary system with nearly identical primary and secondary K5 stars. Radial velocity measurements from the Goodman Spectrograph on the 4.1 meter SOAR telescope of the likely-lowest-mass targets reveal that six of the eclipsing binary discoveries are low-mass (.06 - .37 solar mass) secondaries with K-dwarf primaries, strong candidates for precision mass-radius measurements.

The Evryscope is a telescope array designed to open a new parameter space in optical astronomy, detecting short timescale events across extremely large sky areas simultaneously. The system consists of a 780 MPix 22-camera array with an 8150 sq. deg. field of view, 13" per pixel sampling, and the ability to detect objects down to Mg=16 in each 2 minute dark-sky exposure. The Evryscope, covering 18,400 sq.deg. with hours of high-cadence exposure time each night, is designed to find the rare events that require all-sky monitoring, including transiting exoplanets around exotic stars like white dwarfs and hot subdwarfs, stellar activity of all types within our galaxy, nearby supernovae, and other transient events such as gamma ray bursts and gravitational-wave electromagnetic counterparts. The system averages 5000 images per night with ~300,000 sources per image, and to date has taken over 3.0M images, totaling 250TB of raw data. The resulting light curve database has light curves for 9.3M targets, averaging 32,600 epochs per target through 2018. This paper summarizes the hardware and performance of the Evryscope, including the lessons learned during telescope design, electronics design, a procedure for the precision polar alignment of mounts for Evryscope-like systems, robotic control and operations, and safety and performance-optimization systems. We measure the on-sky performance of the Evryscope, discuss its data-analysis pipelines, and present some example variable star and eclipsing binary discoveries from the telescope. We also discuss new discoveries of very rare objects including 2 hot subdwarf eclipsing binaries with late M-dwarf secondaries (HW Vir systems), 2 white dwarf / hot subdwarf short-period binaries, and 4 hot subdwarf reflection binaries. We conclude with the status of our transit surveys, M-dwarf flare survey, and transient detection.

We report improved masses, radii, and densities for four planets in two bright M-dwarf systems, K2-3 and GJ3470, derived from a combination of new radial velocity and transit observations. Supplementing K2 photometry with follow-up Spitzer transit observations refined the transit ephemerides of K2-3 b, c, and d by over a factor of 10. We analyze ground-based photometry from the Evryscope and Fairborn Observatory to determine the characteristic stellar activity timescales for our Gaussian Process fit, including the stellar rotation period and activity region decay timescale. The stellar rotation signals for both stars are evident in the radial velocity data and are included in our fit using a Gaussian process trained on the photometry. We find the masses of K2-3 b, K2-3 c and GJ3470 b to be 6.48$^{+0.99}_{-0.93}$, 2.14$^{+1.08}_{-1.04}$, and 12.58$^{+1.31}_{-1.28}$ M$_\oplus$ respectively. K2-3 d was not significantly detected and has a 3-$\sigma$ upper limit of 2.80 M$_\oplus$. These two systems are training cases for future TESS systems; due to the low planet densities ($\rho$ $<$ 3.7 g cm$^{-3}$) and bright host stars (K $<$ 9 mag), they are among the best candidates for transmission spectroscopy in order to characterize the atmospheric compositions of small planets.

There have been many efforts to correct systematic effects in astronomical light curves to improve the detection and characterization of planetary transits and astrophysical variability. Algorithms like the Trend Filtering Algorithm (TFA) use simultaneously-observed stars to remove systematic effects, and binning is used to reduce high-frequency random noise. We present TFAW, a wavelet-based modified version of TFA. TFAW aims to increase the periodic signal detection and to return a detrended and denoised signal without modifying its intrinsic characteristics. We modify TFA's frequency analysis step adding a Stationary Wavelet Transform filter to perform an initial noise and outlier removal and increase the detection of variable signals. A wavelet filter is added to TFA's signal reconstruction to perform an adaptive characterization of the noise- and trend-free signal and the noise contribution at each iteration while preserving astrophysical signals. We carried out tests over simulated sinusoidal and transit-like signals to assess the effectiveness of the method and applied TFAW to real light curves from TFRM. We also studied TFAW's application to simulated multiperiodic signals, improving their characterization. TFAW improves the signal detection rate by increasing the signal detection efficiency (SDE) up to a factor ~2.5x for low SNR light curves. For simulated transits, the transit detection rate improves by a factor ~2-5x in the low-SNR regime compared to TFA. TFAW signal approximation performs up to a factor ~2x better than bin averaging for planetary transits. The standard deviations of simulated and real TFAW light curves are ~40x better than TFA. TFAW yields better MCMC posterior distributions and returns lower uncertainties, less biased transit parameters and narrower (~10x) credibility intervals for simulated transits. We present a newly-discovered variable star from TFRM.

Proxima b is a terrestrial-mass planet in the habitable-zone of Proxima Centauri. Proxima Centauri's high stellar activity however casts doubt on the habitability of Proxima b: sufficiently bright and frequent flares and any associated proton events may destroy the planet's ozone layer, allowing lethal levels of UV flux to reach its surface. In March 2016, the Evryscope observed the first naked-eye-brightness superflare detected from Proxima Centauri. Proxima increased in optical flux by a factor of ~68 during the superflare and released a bolometric energy of 10^33.5 erg, ~10X larger than any previously-detected flare from Proxima. Over the last two years the Evryscope has recorded 23 other large Proxima flares ranging in bolometric energy from 10^30.6 erg to 10^32.4 erg; coupling those rates with the single superflare detection, we predict at least five superflares occur each year. Simultaneous high-resolution HARPS spectroscopy during the Evryscope superflare constrains the superflare's UV spectrum and any associated coronal mass ejections. We use these results and the Evryscope flare rates to model the photochemical effects of NOx atmospheric species generated by particle events from this extreme stellar activity, and show that the repeated flaring may be sufficient to reduce the ozone of an Earth-like atmosphere by 90% within five years; complete depletion may occur within several hundred kyr. The UV light produced by the Evryscope superflare would therefore have reached the surface with ~100X the intensity required to kill simple UV-hardy microorganisms, suggesting that life would have to undergo extreme adaptations to survive in the surface areas of Proxima b exposed to these flares.

Certain instrumental effects and data reduction anomalies introduce systematic errors in photometric time-series. Detrending algorithms such as the Trend Filtering Algorithm (TFA) (Kov\'{a}cs et al. 2004) have played a key role in minimizing the effects caused by these systematics. Here we present the results obtained after applying the TFA, Savitszky-Golay (Savitzky & Golay 1964) detrending algorithms and the Box Least Square phase folding algorithm (Kov\'{a}cs et al. 2002) to the TFRM-PSES data (Fors et al. 2013). Tests performed on this data show that by applying these two filtering methods together, the photometric RMS is on average improved by a factor of 3-4, with better efficiency towards brighter magnitudes, while applying TFA alone yields an improvement of a factor 1-2. As a result of this improvement, we are able to detect and analyze a large number of stars per TFRM-PSES field which present some kind of variability. Also, after porting these algorithms to Python and parallelizing them, we have improved, even for large data samples, the computing performance of the overall detrending+BLS algorithm by a factor of $\sim$10 with respect to Kov\'{a}cs et al. (2004).