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.

High-resolution spectroscopy and speckle interferometry reveal the young star HD 86588 as a quadruple system with a 3-tier hierarchy. The 0.3" resolved binary A,B with an estimated period around 300 years contains the 8-year pair Aa,Abc (also potentially resolvable), where Ab,Ac is a double-lined binary with equal components, for which we compute the spectroscopic orbit. Despite the short period of 2.4058 day, the orbit of Ab,Ac is eccentric (e=0.086+-0.003). It has a large inclination, but there are no eclipses; only a 4.4 mmag light modulation apparently caused by star spots on the components of this binary is detected with Evryscope. Assuming a moderate extinction of A_V = 0.5 mag and a parallax of 5.2 mas, we find that the stars are on or close to the main sequence (age >10 Myr) and their masses are from 1 to 1.3 solar. We measure the strength of the Lithium line in the visual secondary B which, together with rotation, suggests that the system is younger than 150 Myr. This object is located behind the extension of the Chamaeleon I dark cloud (which explains extinction and interstellar Sodium absorption), but apparently does not belong to it. We propose a scenario where the inner orbit has recently acquired its high eccentricity through dynamical interaction with the outer two components; it is now undergoing rapid tidal circularization on a time scale of ~1 Myr. Alternatively, the eccentricity could be excited quasi-stationary by the outer component Aa.

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.

The gamma-ray binary LS I +61 303 is composed of a Be star and a compact companion orbiting in an eccentric orbit. Variable flux modulated with the orbital period of ~26.5 d has been detected from radio to very high-energy gamma rays. In addition, the system presents a superorbital variability of the phase and amplitude of the radio outburst with a period of ~4.6 yr. We present optical photometric observations of LS I +61 303 spanning ~1.5 yr and contemporaneous Halpha equivalent width (EW Halpha) data. The optical photometry shows, for the first time, that the known orbital modulation suffers a positive orbital phase shift and an increase in flux for data obtained 1-yr apart. This behavior is similar to that already known at radio wavelengths, indicating that the optical flux follows the superorbital variability as well. The orbital modulation of the EW Halpha presents the already known superorbital flux variability but shows, also for the first time, a positive orbital phase shift. In addition, the optical photometry exhibits a lag of ~0.1-0.2 in orbital phase with respect to the EW Halpha measurements at similar superorbital phases, and presents a lag of ~0.1 and ~0.3 orbital phases with respect noncontemperaneous radio and X-ray outbursts, respectively. The phase shifts detected in the orbital modulation of thermal indicators, such as the optical flux and the EW Halpha, are in line with the observed behavior for nonthermal indicators, such as X-ray or radio emission. This shows that there is a strong coupling between the thermal and nonthermal emission processes in the gamma-ray binary LS I +61 303. The orbital phase lag between the optical flux and the EW Halpha is naturally explained considering different emitting regions in the circumstellar disk, whereas the secular evolution might be caused by the presence of a moving one-armed spiral density wave in the disk.

Low-cost mass-produced sensors and optics have recently made it feasible to build telescope arrays which observe the entire accessible sky simultaneously. In this article we discuss the scientific motivation for these telescopes, including exoplanets, stellar variability and extragalactic transients. To provide a concrete example we detail the goals and expectations for the Evryscope, an under-construction 780 MPix telescope which covers 8,660 square degrees in each two-minute exposure; each night, 18,400 square degrees will be continuously observed for an average of approximately 6 hours. Despite its small 61mm aperture, the system's large field of view provides an etendue which is ~10% of LSST. The Evryscope, which places 27 separate individual telescopes into a common mount which tracks the entire accessible sky with only one moving part, will return 1%-precision, many-year-length, high-cadence light curves for every accessible star brighter than mV=16.5, with brighter stars having few-millimagnitude photometric precision in long-term light curves. It will be capable of searching for transiting giant planets around the brightest and most nearby stars, where the planets are much easier to characterize; it will also search for small planets nearby M-dwarfs, for planetary occultations of white dwarfs, and will perform comprehensive nearby microlensing and eclipse-timing searches for exoplanets inaccessible to other planet-finding methods. The Evryscope will also monitor outbursting young stars, white dwarf activity, and stellar activity of all types, along with finding a large sample of very-long-period M-dwarf eclipsing binaries. When relatively rare transients events occur, such as gamma-ray bursts (GRBs), nearby supernovae, or even gravitational wave detections, the array will return minute-by-minute light curves without needing pointing towards the event as it occurs. (abridged)

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).

Current time-domain wide-field sky surveys generally operate with few-degree-sized fields and take many individual images to cover large sky areas each night. We present the design and project status of the Evryscope ("wide-seer"), which takes a different approach: using an array of 7cm telescopes to form a single wide-field-of-view pointed at every part of the accessible sky simultaneously and continuously. The Evryscope is a gigapixel-scale imager with a 9060 sq. deg. field of view and has an etendue three times larger than the Pan-STARRS sky survey. The system will search for transiting exoplanets around bright stars, M-dwarfs and white dwarfs, as well as detecting microlensing events, nearby supernovae, and gamma-ray burst afterglows. We present the current project status, including an update on the Evryscope prototype telescopes we have been operating for the last three years in the Canadian High Arctic.

The Qatar Exoplanet Survey (QES) is discovering hot Jupiters and aims to discover hot Saturns and hot Neptunes that transit in front of relatively bright host stars. QES currently operates a robotic wide-angle camera system to identify promising transiting exoplanet candidates among which are the confirmed exoplanets Qatar 1b and 2b. This paper describes the first generation QES instrument, observing strategy, data reduction techniques, and follow-up procedures. The QES cameras in New Mexico complement the SuperWASP cameras in the Canary Islands and South Africa, and we have developed tools to enable the QES images and light curves to be archived and analysed using the same methods developed for the SuperWASP datasets. With its larger aperture, finer pixel scale, and comparable field of view, and with plans to deploy similar systems at two further sites, the QES, in collaboration with SuperWASP, should help to speed the discovery of smaller radius planets transiting bright stars in northern skies.

We report on 13 sub-arcsecond binaries, detected by means of lunar occultations in the near-infrared at the ESO Very Large Telescope. They are all first-time detections, except for the visual binary HD~158122 which we resolved for the first time in the near infrared. The primaries have magnitudes in the range $K$=4.5 to 10.0, and the companions in the range $K$=6.8 to 11.1. The magnitude differences have a median value of 2.4, with the largest being 4.6. The projected separations are in the range 4 to 168 milliarcseconds with a median of 13 milliarcseconds. We discuss and compare our results with the available literature. With this paper, we conclude the mining for binary star detections in the volume of 1226 occultations recorded at the VLT with the ISAAC instrument. We expect that the majority of these binaries may be unresolvable by adaptive optics on current telescopes, and they might be challenging for long-baseline interferometry. However they constitute an interesting sample for future larger telescopes and for astrometric missions such as GAIA.

We report on 25 sub-arcsecond binaries, detected for the first time by means of lunar occultations in the near-infrared as part of a long-term program using the ISAAC instrument at the ESO Very Large Telescope. The primaries have magnitudes in the range K=3.8 to 10.4, and the companions in the range K=6.4 to 12.1. The magnitude differences have a median value of 2.8, with the largest being 5.4. The projected separations are in the range 6 to 748 milliarcseconds and with a median of 18 milliarcseconds, or about 3 times less than the diffraction limit of the telescope. Among our binary detections are a pre-main sequence star and an enigmatic Mira-like variable previously suspected to have a companion. Additionally, we quote an accurate first-time near-IR detection of a previously known wider binary. We discuss our findings on an individual basis as far as made possible by the available literature, and we examine them from a statistical point of view. We derive a typical frequency of binarity among field stars of ~10%, in the resolution and sensitivity range afforded by the technique (~0.003" to ~0.5", and K~12mag, respectively). This is in line with previous results by the same technique but we point out interesting differences that we can trace up to sensitivity, time sampling, and average distance of the targets. Finally, we discuss the prospects for further follow-up studies.

A Baker-Nunn Camera (BNC), originally installed at the Real Instituto y Observatorio de la Armada (ROA) in 1958, was refurbished and robotized. The new facility, called Telescope Fabra ROA Montsec (TFRM), was installed at the Observatori Astron\`omic del Montsec (OAdM). The process of refurbishment is described in detail. Most of the steps of the refurbishment project were accomplished by purchasing commercial components, which involve little posterior engineering assembling work. The TFRM is a 0.5m aperture f/0.96 optically modified BNC, which offers a unique combination of instrumental specifications: fully robotic and remote operation, wide-field of view (4.4 deg x 4.4 deg), moderate limiting magnitude (V~19.5mag), ability of tracking at arbitrary right ascension and declination rates, as well as opening and closing CCD shutter at will during an exposure. Nearly all kind of image survey programs can benefit from those specifications. Apart from other less time consuming programs, since the beginning of science TFRM operations we have been conducting two specific and distinct surveys: super-Earths transiting around M-type dwarfs stars, and geostationary debris in the context of Space Situational Awareness / Space Surveillance and Tracking (SSA/SST) programs. Preliminary results for both cases will be shown.

Four gamma-ray binaries, namely PSR B1259-63, HESS J0632+057, HD 215227 and LS I +61 303, contain compact objects orbiting around massive Be stars. The nature of the compact object is only known in the case of PSR B1259-63, but the other systems could also contain young non-accreting pulsars with relativistic winds. Around periastron passage the compact objects should produce significant changes in the structure of the Be discs due to gravitational forces and eventually by ram pressure from the putative pulsar wind. Indeed, variability in the Halpha emission line has been detected in all these systems, and periodic variability in the optical photometry has been detected in two of them. However, there is lack of a systematic monitoring with accurate photometry, which could be used to constrain the shape of the disc during the periastron passage. This information is important to build accurate physical models to explain the broadband spectral energy distribution of these sources. Here we present an ongoing program to monitor the optical photometry of gamma-ray binaries and we show preliminary results for the case of HD 215227.

Calibration is one of the long-standing problems in optical interferometric measurements, particularly with long baselines which demand stars with angular sizes on the milliarcsecond scale and no detectable companions. While systems of calibrators have been generally established for the near-infrared in the bright source regime (K$\la 3$\,mag), modern large interferometers are sensitive to significantly fainter magnitudes. We aim at providing a list of sources found unresolved from direct observations with high angular resolution and dynamic range, which can be used to choose interferometric calibrators. To this purpose, we have used a large number of lunar occultations recorded with the ISAAC instrument at the VLT to select sources found to be unresolved and without close companions. An algorithm has been used to determine the limiting angular resolution achieved for each source, taking into account a noise model built from occulted and unocculted portions of the light curves. We have obtained upper limits on the angular sizes of 556 sources, with magnitudes ranging from K$_{\rm s} \approx$4 to 10, with a median of 7.2\,mag. The upper limits on possible undetected companions (within $\approx 0\farcs5$) range from K$_{\rm s} \approx$8 to 13, with a median of 11.5\,mag. One-third of the sources have angular sizes $\le 1$, and two-thirds $\le 2$ milliarcseconds. This list of unresolved sources matches well the capabilities of current large interferometric facilities. We also provide available cross-identifications, magnitudes, spectral types, and other auxiliary information. A fraction of the sources are found to be potentially variable. The list covers parts of the Galactic Bulge and in particular the vicinity of the Galactic Center, where extinction is very significant and traditional lists of calibrators are often insufficient.

The telescope Fabra ROA at Montsec (TFRM) is a 0.5m f/1 refurbished Baker-Nunn Camera (BNC) operated by a collaboration between the Fabra Observatory - Royal Academy of Arts and Sciences of Barcelona and the Spanish Navy Observatory (ROA), and installed at Montsec Astronomical Observatory (Spain). Among other capabilities, its CCD FoV (4.4{\deg}x4.4{\deg}), the telescope tracking at arbitrary RA and DEC rates, and the CCD shutter commanding at will during the exposure are specially remarkable for Space Surveillance and Tracking (SST) observational programs. On Feb 2011, the TFRM participated, in the CO-VI third run satellite tracking campaign of the ESA SST/Space Surveillance Awareness Preparatory Program (SST/SSA-PP). During this multi-asset 7-day campaign the TFRM conducted systematic observations of artificial satellites which yielded to the determination of 1137 accurate position measurements. Since Feb 2011, the TFRM is observing in remote and fully unattended robotic modes under commissioning status. A summary of the results of the ESA CO-VI SST optical observational campaign and insights of other SST-like observations in process will be presented.

Since the beginning of the Space Age optical sensors have been one of the main instruments for positioning and tracking known space objects. Nowadays, the unrelenting growth of man-made objects together with the overcrowding of the useful satellite orbits, and the real space debris and NEO hazards, has made necessary to carry out surveys of the space looking for uncatalogued objects. Optical telescopes play a key role in the Space Surveillance and Tracking (SST) as a primary Space Situational Awareness element and, it is known, that the best instrument for this task is a fully robotic wide-field telescope with a minimum aperture of 40cm. The Baker-Nunn Cameras (BNCs) were produced by the Smithsonian Institution during the late 50s as an optical tracking system for artificial satellites. These wide-field telescopes of 50cm of aperture were manufactured by Perkin- Elmer (optics) and Boller & Chivens (mechanics) with the highest quality specifications. The TFRM is a fully robotic refurbished BNC that exploits the excellent mechanical and optical original design to obtain an equatorial telescope with a useful 4.4{\deg}x4.4{\deg} CCD field of view. TFRM is therefore a European asset very well suited for SST (satellites and space debris) and NEO observations. Moreover, its control system allows tracking of all kinds of orbits, including LEOs.

Lunar occultations (LO) provide a unique combination of high angular resolution and sensitivity at near-infrared wavelenghts. At the ESO Very Large Telescope, it is possible to achieve about 1 milliarcsecond (mas) resolution and detect sources as faint as K$\approx$12\,mag. We have taken advantage of a passage of the Moon over two crowded and reddened regions in the direction of the inner part of the galactic bulge, in order to obtain a high number of occultation light curves over two half nights. Our goal was to detect and characterize new binary systems, and to investigate highly extincted and relatively unknown infrared sources in search of circumstellar shells and similar peculiarities. Our target list included a significant number of very late-type stars, but in fact the majority of the sources was without spectral classification. A LO event requires the sampling of the light curve at millisecond rates in order to permit a detailed study of the diffraction fringes. For this, we used the so-called burst mode of the ISAAC instrument at the Melipal telescope. Our observing efficiency was ultimately limited by overheads for telescope pointing and data storage, to about one event every three minutes. We could record useful light curves for 184 sources. Of these, 24 were found to be binaries or multiples, all previously unknown. The projected separations are as small as 7.5\,mas, and the magnitude differences as high as $\Delta$K=6.5\,mag. Additionally we could establish, also for the first time, the resolved nature of at least two more stars, with indication of circumstellar emission. We could also put upper limits on the angular size of about 165 unresolved stars, an information that combined with previous and future observations will be very helpful in establishing a list of reliable calibrators for long baseline interferometers.

(Abridged) Lunar occultations (LO) are a very efficient and powerful technique, that achieves the best combination of high angular resolution and sensitivity possible today at near-infrared wavelengths. Given that the events are fixed in time, that the sources are occulted randomly, and that the telescope use is minimal, the technique is very well suited for service mode observations. We have established a program of routine LO observations at the VLT observatory, especially designed to take advantage of short breaks available in-between other programs. We have used the ISAAC instrument in burst mode, capable of producing continuous read-outs at millisecond rates on a suitable subwindow. Given the random nature of the source selection, our aim has been primarily the investigation of a large number of stellar sources at the highest angular resolution in order to detect new binaries. Serendipitous results such as resolved sources and detection of circumstellar components were also anticipated. We have recorded the signal from background stars for a few seconds, around the predicted time of occultation by the Moon's dark limb. At millisecond time resolution, a characteristic diffraction pattern can be observed. Patterns for two or more sources superimpose linearly, and this property is used for the detection of binary stars. The detailed analysis of the diffraction fringes can be used to measure specific properties such as the stellar angular size and the presence of extended light sources such as a circumstellar shell. We present a list of 191 stars for which LO data could be recorded and analyzed. Results include the detection of 16 binary and 2 triple stars, all but one of which were previously unknown. The projected angular separations are as small as 4 milliarcseconds and magnitude differences as high as ?K=5.8 mag...

A number of Baker-Nunn Camera (BNC) were manufactured by Smithsonian Institution during the 60s as optical tracking systems for artificial satellites with optimal optical and mechanical specifications. One of them was installed at the Real Instituto y Observatorio de la Armada (ROA). We have conducted a profound refurbishment project of the telescope to be installed at Observatori Astron\`omic del Montsec (OAdM). As a result, the BNC offers the largest combination of a huge FOV (4.4$\deg$x4.4$\deg$) and aperture (leading to a limiting magnitude of V$\sim$20). These specifications, together with their remote and robotic natures, allows this instrument to face an observational program of exoplanets detection by means of transit technique with high signal-to-noise ratio in the appropiate magnitude range.

For the first time, the lunar occultation technique has been employed on a very large telescope in the near-IR with the aim of achieving systematically milliarcsecond resolution on stellar sources. We have demonstrated the burst mode of the ISAAC instrument, using a fast read-out on a small area of the detector to record many tens of seconds of data at a time on fields of few squared arcsec. We have used the opportunity to record a large number of LO events during a passage of the Moon close to the Galactic Center in March 2006. We have developed a data pipeline for the treatment of LO data, including the automated estimation of the main data analysis parameters using a wavelet-based method, and the preliminary fitting and plotting of all light curves. We recorded 51 LO events over about four hours. Of these, 30 resulted of sufficient quality to enable a detailed fitting. We detected two binaries with subarcsec projected separation and three stars with a marginally resolved angular diameter of about 2 mas. Two more SiO masers, were found to be resolved and in one case we could recover the brightness profile of the extended emission, which is well consistent with an optically thin shell. The remaining unresolved stars were used to characterize the performance of the method. The LO technique at a very large telescope is a powerful and efficient method to achieve angular resolution, sensitivity, and dynamic range that are among the best possible today with any technique. The selection of targets is naturally limited and LOs are fixed-time events, however each observation requires only a few minutes including overheads. As such, LOs are ideally suited to fill small gaps of idle time between standard observations.

We present lunar occultation observations obtained in August 2006 with the recently demonstrated burst mode of the ISAAC instrument at the ESO VLT. The results presented here follow the previously reported observations carried out in March 2006 on a similar but unrelated set of sources. Interstellar extinction in the inner regions of the galactic bulge amounts to tens of magnitudes at visual wavelengths. As a consequence, the majority of sources in that area are poorly studied and large numbers of potentially interesting sources such as late-type giants with circumstellar shells, stellar masers, infrared stars, remain excluded from the typical investigations which are carried out in less problematic regions. Also undetected are a large numbers of binaries. By observing LO events in this region, we gain the means to investigate at least a selected number of sources with an unprecedented combination of sensitivity and angular resolution. The LO technique permits to achieve mas resolution with a sensitivity of K~12mag at a very large telescope. We have used the opportunity of a favorable passage of the Moon over a crowded region in the general direction of the Galactic Center to observe 78 LO events of heavily reddened stellar sources. We have detected six new binary and one triple star, with typical projected separation of ~0.01". We have also detected the compact circumstellar emission around one maser and one central star of a planetary nebula. Additionally we have measured the diameter and/or circumstellar shell of two carbon stars and other IR sources. In agreement with our previous result, we conclude that lunar occultations in fast read-out mode on a detector subwindow at an 8m-class telescope can achieve an angular resolution close to 0.001" with a sensitivity K~12mag.

The introduction of infrared arrays for lunar occultations (LO) work and the improvement of predictions based on new deep IR catalogues have resulted in a large increase in the number of observable occultations. We provide the means for an automated reduction of large sets of LO data. This frees the user from the tedious task of estimating first-guess parameters for the fit of each LO lightcurve. At the end of the process, ready-made plots and statistics enable the user to identify sources which appear to be resolved or binary and to initiate their detailed interactive analysis. The pipeline is tailored to array data, including the extraction of the lightcurves from FITS cubes. Because of its robustness and efficiency, the wavelet transform has been chosen to compute the initial guess of the parameters of the lightcurve fit. We illustrate and discuss our automatic reduction pipeline by analyzing a large volume of novel occultation data recorded at Calar Alto Observatory. The automated pipeline package is available from the authors.

The recent implementation of the high-speed burst mode at the ISAAC instrument on UT1, and its propagation to other ESO instruments, has opened the door to observational capabilities which hold the potential for a wealth of novel results. In the ELT era, when the accent will likely be on lengthy programs aimed at the best sensitivity and angular resolution, the VLT telescopes could continue to play a significant and largely unique role by performing routinely observations of transient events at high temporal resolution. In our contribution, we provide details on two such kinds of observations, namely lunar occultations of stars and of asteroids. For the first ones, we report on two passages of the Moon in regions with high stellar density as the Galactic Center. The VLT-UT1 telescope was used for the first time to record successfully 53 and 71 occultations on March 22 and August 6, 2006, with an angular resolution of 0.5-1 milliarcsecond and $K\sim12.5$ limiting magnitude. We note that the angular resolution is superior to that achieved at present by Adaptive Optics on any telescope, and also superior to that foreseen for the ELT at the same wavelength. LO are also very efficient in terms of telescope time. We present some of the results, including the discovery of close binaries, and the detection and study of compact circumstellar components of cool giants, AGB stars and embedded IR sources. Rest of the abstract follows at the paper

(Abridged) In the first part of this thesis, a general methodology for applying image deconvolution to wide-field CCD imagery. Results show that wavelet-based deconvolution can increase limiting magnitude up to 0.6 mag and improve limiting resolution 1 pixel with respect to original image with no astrometric accuracy degradation. In the second part, a new observational technique based on CCD fast drift scanning has been proposed for lunar occultations (LO) and speckle interferometry. This enables all kind of professional and high-end amateur observatories to perform such kind of observations. For LO, 16 new binaries up to 2mas of projected separation were detected and stellar diameters measurements in the 7 mas regime were obtained with that CCD and IR subarray based techniques. A new wavelet-based LO reduction pipeline was implemented. For speckle, CCD fast drift scanning technique was validated with the observation of four binary systems with well determined orbits. The results of separation, position angle and magnitude difference are in accordance with published measurements. A new approach for calibrating speckle transfer function from the binary power spectrum itself has been introduced. It does not require point source observations, which gives a more effective use of observation time.

We present an update of the lunar occultation program which is routinely carried out in the near-IR at the Calar Alto Observatory. A total of 350 events were recorded since our last report (Fors et al. 2004). In the course of eight runs we have observed, among others, late-type giants, T-Tauri stars, and infrared sources. Noteworthy was a passage of the Moon close to the galactic center, which produced a large number of events during just a few hours in July 2004. Results include the determinations of the angular diameter of RZ Ari, and the projected separations and brightness ratios for one triple and 13 binary stars, almost all of which representing first time detections. Projected separations range from 0.09arcsec to 0.007arcsec. We provide a quantitative analysis of the performance achieved in our observations in terms of angular resolution and sensitivity, which reach about 0.003arcsec and K~8.5mag, respectively. We also present a statistical discussion of our sample, and in particular of the frequency of detection of binaries among field stars.

A new application of a fast CCD drift scanning technique that allows us to perform speckle imaging of binary stars is presented. For each observation, an arbitrary number of speckle frames is periodically stored on a computer disk, each with an appropriate exposure time given both atmospheric and instrumental considerations. The CCD charge is shifted towards the serial register and read out sufficiently rapidly to avoid an excessive amount of interframe dead time. Four well-known binary systems (ADS 755, ADS 2616, ADS 3711 and ADS 16836) are observed in to show the feasibility of the proposed technique. Bispectral data analysis and power spectrum fitting is carried out for each observation, yielding relative astrometry and photometry. A new approach for self-calibrating this analysis is also presented and validated. The proposed scheme does not require any additional electronic or optical hardware, so it should allow most small professional observatories and advanced amateurs to enjoy the benefits of diffraction-limited imaging.

We present a program of routine lunar occultations, at optical and near-IR wavelengths, recently started at the 1.5m Spanish telescope at the Calar Alto Observatory. Both a CCD and an infrared array detector are used. The program is aimed mainly at the detection and investigation of binary systems, although results in other areas of stellar research are also anticipated. Occultations are reported for a total of 40 stars. Among these, SAO 164567, SAO 78258 and AG+24 788 have been discovered to be binaries, with projected separations as small as 0.006". Furthermore, binarity is suspected in the case of SAO 78119 and SAO 79251. Additionally, the angular diameter of the late-type giant 30 Psc and of the infrared star V349 Gem have been accurately measured, this latter for the first time. We finally evaluate the instrumentation performance in terms of limiting magnitude and angular resolution, and discuss applications to larger telescopes.

We have compiled optical and radio astrometric data of the microquasar LS 5039 and derived its proper motion. This, together with the distance and radial velocity of the system, allows us to state that this source is escaping from its own regional standard of rest, with a total systemic velocity of about 150 km/s and a component perpendicular to the galactic plane larger than 100 km/s. This is probably the result of an acceleration obtained during the supernova event that created the compact object in this binary system. We have computed the trajectory of LS 5039 in the past, and searched for OB associations and supernova remnants in its path. In particular, we have studied the possible association between LS 5039 and the supernova remnant G016.8-01.1, which, despite our efforts, remains dubious. We have also discovered and studied an HI cavity in the ISM, which could have been created by the stellar wind of LS 5039 or by the progenitor of the compact object in the system. Finally, in the symmetric supernova explosion scenario, we estimate that at least 17 solar masses were lost in order to produce the high eccentricity observed. Such a mass loss could also explain the observed runaway velocity of the microquasar.

A CCD drift-scanning technique for lunar occultation (LO) observations is presented. While this approach has been addressed before by Sturmann (1994) for the case of large telescopes, the technical validity has never been discussed for sub-meter class telescopes. In contrast to Sturmann's scheme, the proposed technique places the CCD in the image plane of the telescope. This does not represent a problem in the case of small telescopes, where the practical angular resolution attainable by LO is not limited by aperture smoothing. Photon-generated charge is read out at millisecond rate on a column-per-column basis, as the diffraction pattern of the occulted star is being tracked. Two LO events (SAO 79031 and SAO 77911) were observed to demonstrate the feasibility of the method. Data analysis was carried out and no angular diameter the observed stars were resolved. We show, however, that the technique could be useful for close binary detections with small telescopes. A discussion of the limiting resolution and magnitude imposed by our instrumentation is carried out, showing that drift-scanning technique could be extended to 1-2m telescopes for stellar diameter determination purposes. Finally, we point out that the technical demands required by this technique can be easily met by most small professional observatories and advanced amateurs.