Wind is the process that connects Mars' climate system. Measurements of Mars atmospheric winds from orbit would dramatically advance our understanding of Mars and help prepare for human exploration of the Red Planet. Multiple instrument candidates are in development and will be ready for flight in the next decade. We urge the Decadal Survey to make these measurements a priority for 2023-2032.

Efficient operation of a submillimeter interferometer requires remote (preferably automated) control of mechanically tuned local oscillators, phase-lock loops, mixers, optics, calibration vanes and cryostats. The present control system for these aspects of the Submillimeter Array (SMA) will be described. Distributed processing forms the underlying architecture and the software is split between hardware platforms in a leader/follower arrangement. In each antenna cabin, a serial network of up to ten independent 80C196 microcontroller boards attaches to the real-time PowerPC computer (running LynxOS). A multi-threaded, gcc-compiled leader program on the PowerPC accepts top-level requests via remote procedure calls (RPC), subsequently dispatches tuning commands to the relevant follower microcontrollers, and regularly reports the system status to optical-fiber-based reflective memory for common access by the telescope monitor and error reporting system. All serial communication occurs asynchronously via encoded, variable-length packets. The microcontrollers respond to the requested commands and queries by accessing non-volatile, rewriteable lookup-tables (when appropriate) and executing embedded software that operates additional electronic devices (DACs, ADCs, etc.). Since various receiver hardware components require linear or rotary motion, each microcontroller also implements a position servo via a one-millisecond interrupt service routine which drives a DC-motor/encoder combination that remains standard across each subsystem.

This paper documents the sixteenth data release (DR16) from the Sloan Digital Sky Surveys; the fourth and penultimate from the fourth phase (SDSS-IV). This is the first release of data from the southern hemisphere survey of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2); new data from APOGEE-2 North are also included. DR16 is also notable as the final data release for the main cosmological program of the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), and all raw and reduced spectra from that project are released here. DR16 also includes all the data from the Time Domain Spectroscopic Survey (TDSS) and new data from the SPectroscopic IDentification of ERosita Survey (SPIDERS) programs, both of which were co-observed on eBOSS plates. DR16 has no new data from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey (or the MaNGA Stellar Library "MaStar"). We also preview future SDSS-V operations (due to start in 2020), and summarize plans for the final SDSS-IV data release (DR17).

The performance of tomographic adaptive optics systems is intrinsically linked to the vertical profile of optical turbulence. Firstly, a sufficient number of discrete turbulent layers must be reconstructed to model the true continuous turbulence profile. Secondly over the course of an observation, the profile as seen by the telescope changes and the tomographic reconstructor must be updated. These changes can be due to the unpredictable evolution of turbulent layers on meteorological timescales as short as minutes. Here we investigate the effect of changing atmospheric conditions on the quality of tomographic reconstruction by coupling fast analytical adaptive optics simulation to a large database of 10 691 high resolution turbulence profiles measured over two years by the Stereo-SCIDAR instrument at ESO Paranal, Chile. This work represents the first investigation of these effects with a large, statistically significant sample of turbulence profiles. The statistical nature of the study allows us to assess not only the degradation and variability in tomographic error with a set of system parameters (e.g. number of layers, temporal update period) but also the required parameters to meet some error threshold. In the most challenging conditions where the profile is rapidly changing, these parameters must be far more tightly constrained in order to meet this threshold. By providing estimates of these constraints for a wide range of system geometries as well as the impact of different temporal optimisation strategies we may assist the designers of tomographic AO for the ELT to dimension their systems.

We describe the implementation of a robotic SLODAR instrument at the Cerro Paranal observatory. The instrument measures the vertical profile of the optical atmospheric turbulence strength, in 8 resolution elements, to a maximum altitude ranging between 100 m and 500 m. We present statistical results of measurements of the turbulence profile on a total of 875 nights between 2014 and 2018. The vertical profile of the ground layer of turbulence is very varied, but in the median case most of the turbulence strength in the ground layer is concentrated within the first 50 m altitude, with relatively weak turbulence at higher altitudes up to 500 m. We find good agreement between measurements of the seeing angle from the SLODAR and from the Paranal DIMM seeing monitor, and also for seeing values extracted from the Shack-Hartmann active optics sensor of VLT UT1, adjusting for the height of each instrument above ground level. The SLODAR data suggest that a median improvement in the seeing angle from 0.689 arcsec to 0.481 arcsec at wavelength 500 nm would be obtained by fully correcting the ground-layer turbulence between the height of the UTs (taken as 10 m) and altitude 500 m.

Gaia16aye was a binary microlensing event discovered in the direction towards the northern Galactic disc and was one of the first microlensing events detected and alerted to by the Gaia space mission. Its light curve exhibited five distinct brightening episodes, reaching up to I=12 mag, and it was covered in great detail with almost 25,000 data points gathered by a network of telescopes. We present the photometric and spectroscopic follow-up covering 500 days of the event evolution. We employed a full Keplerian binary orbit microlensing model combined with the motion of Earth and Gaia around the Sun to reproduce the complex light curve. The photometric data allowed us to solve the microlensing event entirely and to derive the complete and unique set of orbital parameters of the binary lensing system. We also report on the detection of the first-ever microlensing space-parallax between the Earth and Gaia located at L2. The properties of the binary system were derived from microlensing parameters, and we found that the system is composed of two main-sequence stars with masses 0.57$\pm$0.05 $M_\odot$ and 0.36$\pm$0.03 $M_\odot$ at 780 pc, with an orbital period of 2.88 years and an eccentricity of 0.30. We also predict the astrometric microlensing signal for this binary lens as it will be seen by Gaia as well as the radial velocity curve for the binary system. Events such as Gaia16aye indicate the potential for the microlensing method of probing the mass function of dark objects, including black holes, in directions other than that of the Galactic bulge. This case also emphasises the importance of long-term time-domain coordinated observations that can be made with a network of heterogeneous telescopes.

Atmospheric scintillation caused by optical turbulence in the Earth's atmosphere can be the dominant source of noise in ground-based photometric observations of bright targets, which is a particular concern for ground-based exoplanet transit photometry. We demonstrate the implications of atmospheric scintillation for exoplanet transit photometry through contemporaneous turbulence profiling and transit observations. We find a strong correlation between measured intensity variations and scintillation determined through optical turbulence profiling. This correlation indicates that turbulence profiling can be used to accurately model the amount of scintillation noise present in photometric observations on another telescope at the same site. We examine the conditions under which scintillation correction would be beneficial for transit photometry through turbulence profiling, and find that for the atmosphere of La Palma, scintillation dominates for bright targets of magnitude above $V\sim10.1$ mag for a 0.5~m telescope, and at $V\sim11.7$ mag for a 4.2 m telescope under median atmospheric conditions. Through Markov-chain Monte Carlo methods we examine the effect of scintillation noise on the uncertainty of the measured exoplanet parameters, and determine the regimes where scintillation correction is especially beneficial. The ability to model the amount of noise in observations due to scintillation, given an understanding of the atmosphere, is a crucial test for our understanding of scintillation and the overall noise budget of our observations.

Global dust storms are the most thermodynamically significant dust events on Mars. They are produced from the combination of multiple local and regional lifting events and maintained by positive radiative-dynamic feedbacks. The most recent of these events, which began in June 2018, was monitored by several spacecraft in orbit and on the surface, but many questions remain regarding its onset, expansion and decay. We model the 2018 global dust storm with the NASA Ames Mars Global Climate Model to better understand the evolution of the storm and how the general circulation and finite surface dust reservoirs impact it. The global dust storm is characterized by the rapid eastward transport of dust in the equatorial regions and subsequent lifting. We highlight the rapid transfer of dust between western and eastern hemispheres reservoirs, which may play an important role in the storm development through the replenishment of surface dust. Both the Hadley cell circulation and the diurnal cycle of atmospheric heating increase in intensity with increasing dustiness. Large dust plumes are predicted during the mature stage of the storm, injecting dust up to 80 km. The water ice cloud condensation level migrates to higher altitudes, leading to the enrichment of water vapor in the upper atmosphere. In our simulations, the intensity of the Hadley cell is significantly stronger than that of non-dusty conditions. This feedback is strongly sensitive to the radiative properties of dust, which depends on the effective size of the lifted dust distribution.

For extremely large telescopes, adaptive optics will be required to correct the Earth's turbulent atmosphere. The performance of tomographic adaptive optics is strongly dependent on the vertical distribution (profile) of this turbulence. An important way in which this manifests is the tomographic error, arising from imperfect measurement and reconstruction of the turbulent phase at altitude. Conventionally, a small number of reference profiles are used to obtain this error in simulation however these profiles are not constructed to be representative in terms of tomographic error. It is therefore unknown whether these simulations are providing realistic performance estimates. Here, we employ analytical adaptive optics simulation that drastically reduces computation times to compute tomographic error for 10 691 measurements of the turbulence profile gathered by the Stereo-SCIDAR instrument at ESO Paranal. We assess for the first time the impact of the profile on tomographic error in a statistical manner. We find, in agreement with previous work, that the tomographic error is most directly linked with the distribution of turbulence into discrete, stratified layers. Reference profiles are found to provide mostly higher tomographic error than expected, which we attribute to the fact that these profiles are primarily composed of averages of many measurements resulting in unrealistic, continuous distributions of turbulence. We propose that a representative profile should be defined with respect to a particular system, and that as such simulations with a large statistical sample of profiles must be an important step in the design process.

We present extensive optical photometric and spectroscopic observations, from 4 to 482 days after explosion, of the Type II-plateau (II-P) supernova (SN) 2017eaw in NGC 6946. SN 2017eaw is a normal SN II-P intermediate in properties between, for example, SN 1999em and SN 2012aw and the more luminous SN 2004et, also in NGC 6946. We have determined that the extinction to SN 2017eaw is primarily due to the Galactic foreground and that the SN site metallicity is likely subsolar. We have also independently confirmed a tip-of-the-red-giant-branch (TRGB) distance to NGC 6946 of 7.73+/-0.78 Mpc. The distances to the SN that we have also estimated via both the standardized candle method and expanding photosphere method corroborate the TRGB distance. We confirm the SN progenitor identity in pre-explosion archival Hubble Space Telescope (HST) and Spitzer Space Telescope images, via imaging of the SN through our HST Target of Opportunity program. Detailed modeling of the progenitor's spectral energy distribution indicates that the star was a dusty, luminous red supergiant consistent with an initial mass of ~15 Msuns.

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.

Advanced adaptive optics (AO) instruments on ground-based telescopes require accurate knowledge of the atmospheric turbulence strength as a function of altitude. This information assists point spread function reconstruction, AO temporal control techniques and is required by wide-field AO systems to optimize the reconstruction of an observed wavefront. The variability of the atmosphere makes it important to have a measure of the optical turbulence profile in real time. This measurement can be performed by fitting an analytically generated covariance matrix to the cross-covariance of Shack-Hartmann wavefront sensor (SHWFS) centroids. In this study we explore the benefits of reducing cross-covariance data points to a covariance map region of interest (ROI). A technique for using the covariance map ROI to measure and compensate for SHWFS misalignments is also introduced. We compare the accuracy of covariance matrix and map ROI optical turbulence profiling using both simulated and on-sky data from CANARY, an AO demonstrator on the 4.2 m William Herschel telescope, La Palma. On-sky CANARY results are compared to contemporaneous profiles from Stereo-SCIDAR - a dedicated high-resolution optical turbulence profiler. It is shown that the covariance map ROI optimizes the accuracy of AO telemetry optical turbulence profiling. In addition, we show that the covariance map ROI reduces the fitting time for an extremely large telescope-scale system by a factor of 72. The software package we developed to collect all of the presented results is now open source.

An HST COS Far UV spectrum (1170 A to 1800 A) was obtained for the short orbital period recurrent novae (T Pyxidis subclass), CI Aquilae. CI Aql is the only classical CV known to have two eclipses of sensible depth per orbit cycle and also have pre- and post-outburst light curves that are steady enough to allow estimates of mass and orbital period changes. Our FUV spectral analysis with model accretion disks and NLTE high gravity photospheres, together with the Gaia parallax, reveal CI Aql's FUV light is dominated by an optically thick accretion disk with an accretion rate of the order of $4\times 10^{-8}$ $M_{\odot}/yr$. Its database of light curves, radial velocity curves, and eclipse timings is among the best for any CV. Its orbit period ($P$), $dP/dt$, and reference time are re-derived via simultaneous analysis of the three data types, giving a dimensionless post-outburst $dP/dt$ of $-2.49\pm 0.95\times 10^{-10}$. Lack of information on loss of orbital to rotational angular momentum leads to some uncertainty in the translation of $dP/dt$ to white dwarf mass change rate, $dM_1/dt$, but within the modest range of $+4.8\times 10^{-8}$ to $+7.8\times 10^{-8}$ $M_{\odot} /yr$. The estimated white dwarf mass change through outburst for CI Aql, based on simple differencing of its pre- and post outburst orbit period, is unchanged from the previously published $+5.3 \times 10^{-6} M_{\odot}$. At the WD's estimated mass increase rate, it will terminate as a Type Ia supernova within 10 million years.

Knowledge of the optical turbulence profile is important in adaptive optics (AO) systems, particularly tomographic AO systems such as those to be employed by the next generation of 40 m class extremely large telescopes (ELTs). Site characterisation and monitoring campaigns have produced large quantities of turbulence profiling data for sites around the world. However AO system design and performance characterisation is dependent on Monte-Carlo simulations that cannot make use of these large datasets due to long computation times. Here we address the question of how to reduce these large datasets into small sets of profiles that can feasibly be used in such Monte-Carlo simulations, whilst minimising the loss of information inherent in this effective compression of the data. We propose hierarchical clustering to partition the dataset according to the structure of the turbulence profiles and extract a single profile from each cluster. This method is applied to the Stereo-SCIDAR dataset from ESO Paranal containing over 10000 measurements of the turbulence profile from 83 nights. We present two methods of extracting turbulence profiles from the clusters, resulting in two sets of 18 profiles providing subtly different descriptions of the variability across the entire dataset. For generality we choose integrated parameters of the turbulence to measure the representativeness of our profiles and compare to others. Using these criterion we also show that such variability is difficult to capture with small sets of profiles associated with integrated turbulence parameters such as seeing.

HiPERCAM is a quintuple-beam imager that saw first light on the 4.2m William Herschel Telescope (WHT) in October 2017 and on the 10.4m Gran Telescopio Canarias (GTC) in February 2018. The instrument uses re-imaging optics and 4 dichroic beamsplitters to record ugriz (300-1000nm) images simultaneously on its five CCD cameras. The detectors in HiPERCAM are frame-transfer devices cooled thermo-electrically to -90degC, thereby allowing both long-exposure, deep imaging of faint targets, as well as high-speed (over 1000 windowed frames per second) imaging of rapidly varying targets. In this paper, we report on the as-built design of HiPERCAM, its first-light performance on the GTC, and some of the planned future enhancements.

Using spectroscopic radial velocities with the APOGEE instrument and Gaia distance estimates, we demonstrate that Kepler-503b, currently considered a validated Kepler planet, is in fact a brown-dwarf/low-mass star in a nearly circular 7.2-day orbit around a subgiant star. Using a mass estimate for the primary star derived from stellar models, we derive a companion mass and radius of $0.075\pm0.003 \ M_{\odot}$ ($78.6\pm3.1 \ M_{Jup}$) and $0.099^{+0.006}_{-0.004}\ R_{\odot}$ ($0.96^{+0.06}_{-0.04}\ R_{Jup}$), respectively. Assuming the system is coeval, the evolutionary state of the primary indicates the age is $\sim6.7$ Gyr. Kepler-503b sits right at the hydrogen burning mass limit, straddling the boundary between brown dwarfs and very low-mass stars. More precise radial velocities and secondary eclipse spectroscopy with James Webb Space Telescope will provide improved measurements of the physical parameters and age of this important system to better constrain and understand the physics of these objects and their spectra. This system emphasizes the value of radial velocity observations to distinguish a genuine planet from astrophysical false positives, and is the first result from the SDSS-IV monitoring of Kepler planet candidates with the multi-object APOGEE instrument.

Knowledge of the Earth's atmospheric optical turbulence is critical for astronomical instrumentation. Not only does it enable performance verification and optimisation of existing systems but it is required for the design of future instruments. As a minimum this includes integrated astro-atmospheric parameters such as seeing, coherence time and isoplanatic angle, but for more sophisticated systems such as wide field adaptive optics enabled instrumentation the vertical structure of the turbulence is also required. Stereo-SCIDAR is a technique specifically designed to characterise the Earth's atmospheric turbulence with high altitude resolution and high sensitivity. Together with ESO, Durham University has commissioned a Stereo-SCIDAR instrument at Cerro Paranal, Chile, the site of the Very Large Telescope (VLT), and only 20~km from the site of the future Extremely Large Telescope (ELT). Here we provide results from the first 18 months of operation at ESO Paranal including instrument comparisons and atmospheric statistics. Based on a sample of 83 nights spread over 22 months covering all seasons, we find the median seeing to be 0.64" with 50% of the turbulence confined to an altitude below 2 km and 40% below 600 m. The median coherence time and isoplanatic angle are found as 4.18 ms and 1.75" respectively. A substantial campaign of inter-instrument comparison was also undertaken to assure the validity of the data. The Stereo-SCIDAR profiles (optical turbulence strength and velocity as a function of altitude) have been compared with the Surface-Layer SLODAR, MASS-DIMM and the ECMWF weather forecast model. The correlation coefficients are between 0.61 (isoplanatic angle) and 0.84 (seeing).

Future remote sensing of exoplanets will be enhanced by a thorough investigation of our solar system Ice Giants (Neptune-size planets). What can the configuration of the magnetic field tell us (remotely) about the interior, and what implications does that field have for the structure of the magnetosphere; energy input into the atmosphere, and surface geophysics (for example surface weathering of satellites that might harbour sub-surface oceans). How can monitoring of auroral emission help inform future remote observations of emission from exoplanets? Our Solar System provides the only laboratory in which we can perform in-situ experiments to understand exoplanet formation, dynamos, systems and magnetospheres.

The Apache Point Observatory Galactic Evolution Experiment (APOGEE) has observed $\sim$600 transiting exoplanets and exoplanet candidates from \textit{Kepler} (Kepler Objects of Interest, KOIs), most with $\geq$18 epochs. The combined multi-epoch spectra are of high signal-to-noise (typically $\geq$100) and yield precise stellar parameters and chemical abundances. We first confirm the ability of the APOGEE abundance pipeline, ASPCAP, to derive reliable [Fe/H] and effective temperatures for FGK dwarf stars -- the primary \textit{Kepler} host stellar type -- by comparing the ASPCAP-derived stellar parameters to those from independent high-resolution spectroscopic characterizations for 221 dwarf stars in the literature. With a sample of 282 close-in ($P<100$ days) KOIs observed in the APOGEE KOI goal program, we find a correlation between orbital period and host star [Fe/H] characterized by a critical period, $P_\mathrm{crit}$= $8.3^{+0.1}_{-4.1}$ days, below which small exoplanets orbit statistically more metal-enriched host stars. This effect may trace a metallicity dependence of the protoplanetary disk inner-radius at the time of planet formation or may be a result of rocky planet ingestion driven by inward planetary migration. We also consider that this may trace a metallicity dependence of the dust sublimation radius, but find no statistically significant correlation with host $T_\mathrm{eff}$ and orbital period to support such a claim.

We report the discovery of four Fast Radio Bursts (FRBs) in the ongoing SUrvey for Pulsars and Extragalactic Radio Bursts (SUPERB) at the Parkes Radio Telescope: FRBs 150610, 151206, 151230 and 160102. Our real-time discoveries have enabled us to conduct extensive, rapid multi-messenger follow-up at 12 major facilities sensitive to radio, optical, X-ray, gamma-ray photons and neutrinos on time scales ranging from an hour to a few months post-burst. No counterparts to the FRBs were found and we provide upper limits on afterglow luminosities. None of the FRBs were seen to repeat. Formal fits to all FRBs show hints of scattering while their intrinsic widths are unresolved in time. FRB 151206 is at low Galactic latitude, FRB 151230 shows a sharp spectral cutoff, and FRB 160102 has the highest dispersion measure (DM = $2596.1\pm0.3$ pc cm$^{-3}$) detected to date. Three of the FRBs have high dispersion measures (DM >$1500$ pc cm$^{-3}$), favouring a scenario where the DM is dominated by contributions from the Intergalactic Medium. The slope of the Parkes FRB source counts distribution with fluences $>2$ Jyms is $\alpha=-2.2^{+0.6}_{-1.2}$ and still consistent with a Euclidean distribution ($\alpha=-3/2$). We also find that the all-sky rate is $1.7^{+1.5}_{-0.9}\times10^3$FRBs/($4\pi$ sr)/day above $\sim2$ Jyms and there is currently no strong evidence for a latitude-dependent FRB sky-rate.

As telescopes become larger, into the era of ~40 m Extremely Large Telescopes, the high- resolution vertical profile of the optical turbulence strength is critical for the validation, optimization and operation of optical systems. The velocity of atmospheric optical turbulence is an important parameter for several applications including astronomical adaptive optics systems. Here, we compare the vertical profile of the velocity of the atmospheric wind above La Palma by means of a comparison of Stereo-SCIntillation Detection And Ranging (Stereo- SCIDAR) with the Global Forecast System models and nearby balloon-borne radiosondes. We use these data to validate the automated optical turbulence velocity identification from the Stereo-SCIDAR instrument mounted on the 2.5 m Isaac Newton Telescope, La Palma. By comparing these data we infer that the turbulence velocity and the wind velocity are consistent and that the automated turbulence velocity identification of the Stereo-SCIDAR is precise. The turbulence velocities can be used to increase the sensitivity of the turbulence strength profiles, as weaker turbulence that may be misinterpreted as noise can be detected with a velocity vector. The turbulence velocities can also be used to increase the altitude resolution of a detected layer, as the altitude of the velocity vectors can be identified to a greater precision than the native resolution of the system. We also show examples of complex velocity structure within a turbulent layer caused by wind shear at the interface of atmospheric zones.

We present new time-resolved photometry of 74 cataclysmic variables (CVs), 47 of which are eclipsing. 13 of these eclipsing systems are newly discovered. For all 47 eclipsing systems we show high cadence (1-20 seconds) light curves obtained with the high-speed cameras ultracam and ultraspec. We provide new or refined ephemerides, and supply mid-eclipse times for all observed eclipses. We assess the potential for light curve modelling of all 47 eclipsing systems to determine their system parameters, finding 20 systems which appear to be suitable for future study.

A 32 GHz bandwidth VLBI capable correlator and phased array has been designed and deployed at the Smithsonian Astrophysical Observatory's Submillimeter Array (SMA). The SMA Wideband Astronomical ROACH2 Machine (SWARM) integrates two instruments: a correlator with 140 kHz spectral resolution across its full 32 GHz band, used for connected interferometric observations, and a phased array summer used when the SMA participates as a station in the Event Horizon Telescope (EHT) Very Long Baseline Interferometry (VLBI) array. For each SWARM quadrant, Reconfigurable Open Architecture Computing Hardware (ROACH2) units shared under open source from the Collaboration for Astronomy Signal Processing and Electronics Research (CASPER) are equipped with a pair of ultra-fast Analog-to- Digital Converters (ADCs), a Field Programmable Gate Array (FPGA) processor, and eight 10 Gigabit Ethernet ports. A VLBI data recorder interface designated the SWARM Digital Back End, or SDBE, is implemented with a ninth ROACH2 per quadrant, feeding four Mark6 VLBI recorders with an aggregate recording rate of 64 Gbps. This paper describes the design and implementation of SWARM, as well as its deployment at SMA with reference to verification and science data.

After more than 26 years in quiescence, the black hole transient V404 Cyg went into a luminous outburst in June 2015, and additional activity was detected in late December of the same year. Here, we present an optical spectroscopic follow-up of the December mini-outburst, together with X-ray, optical and radio monitoring that spanned more than a month. Strong flares with gradually increasing intensity are detected in the three spectral ranges during the ~10 days following the Swift trigger. Our optical spectra reveal the presence of a fast outflowing wind, as implied by the detection of a P-Cyg profile (He I - 5876 A) with a terminal velocity of ~2500 km/s. Nebular-like spectra -- with an H_alpha equivalent width of ~500 A -- are also observed. All these features are similar to those seen during the main June 2015 outburst. Thus, the fast optical wind simultaneous with the radio jet is most likely present in every V404 Cyg outburst. Finally, we report on the detection of a strong radio flare in late January 2016, when X-ray and optical monitoring had stopped due to Sun constraints.

We present long-term photometric observations of the young open cluster IC 348 with a baseline time-scale of 2.4 yr. Our study was conducted with several telescopes from the Young Exoplanet Transit Initiative (YETI) network in the Bessel $R$ band to find periodic variability of young stars. We identified 87 stars in IC 348 to be periodically variable; 33 of them were unreported before. Additionally, we detected 61 periodic non-members of which 41 are new discoveries. Our wide field of view was the key to those numerous newly found variable stars. The distribution of rotation periods in IC 348 has always been of special interest. We investigate it further with our newly detected periods but we cannot find a statistically significant bimodality. We also report the detection of a close eclipsing binary in IC 348 composed of a low-mass stellar component ($M \gtrsim 0.09\,\mathrm{M}_{\odot}$) and a K0 pre-main sequence star ($M \approx 2.7\,\mathrm{M}_{\odot}$). Furthermore, we discovered three detached binaries among the background stars in our field of view and confirmed the period of a fourth one.

HiPERCAM is a high-speed camera for the study of rapid variability in the Universe. The project is funded by a 3.5MEuro European Research Council Advanced Grant. HiPERCAM builds on the success of our previous instrument, ULTRACAM, with very significant improvements in performance thanks to the use of the latest technologies. HiPERCAM will use 4 dichroic beamsplitters to image simultaneously in 5 optical channels covering the u'g'r'i'z' bands. Frame rates of over 1000 per second will be achievable using an ESO CCD controller (NGC), with every frame GPS timestamped. The detectors are custom-made, frame-transfer CCDs from e2v, with 4 low-noise (2.5e-) outputs, mounted in small thermoelectrically-cooled heads operated at 180 K, resulting in virtually no dark current. The two reddest CCDs will be deep-depletion devices with anti-etaloning, providing high quantum efficiencies across the red part of the spectrum with no fringing. The instrument will also incorporate scintillation noise correction via the conjugate-plane photometry technique. The opto-mechanical chassis will make use of additive manufacturing techniques in metal to make a light-weight, rigid and temperature-invariant structure. First light is expected on the 4.2m William Herschel Telescope on La Palma in 2017 (on which the field of view will be 10' with a 0.3"/pixel scale), with subsequent use planned on the 10.4m Gran Telescopio Canarias on La Palma (on which the field of view will be 4' with a 0.11"/pixel scale) and the 3.5m New Technology Telescope in Chile.

We present a long-term programme for timing the eclipses of white dwarfs in close binaries to measure apparent and/or real variations in their orbital periods. Our programme includes 67 close binaries, both detached and semi-detached and with M-dwarfs, K-dwarfs, brown dwarfs or white dwarfs secondaries. In total, we have observed more than 650 white dwarf eclipses. We use this sample to search for orbital period variations and aim to identify the underlying cause of these variations. We find that the probability of observing orbital period variations increases significantly with the observational baseline. In particular, all binaries with baselines exceeding 10 yrs, with secondaries of spectral type K2 -- M5.5, show variations in the eclipse arrival times that in most cases amount to several minutes. In addition, among those with baselines shorter than 10 yrs, binaries with late spectral type (>M6), brown dwarf or white dwarf secondaries appear to show no orbital period variations. This is in agreement with the so-called Applegate mechanism, which proposes that magnetic cycles in the secondary stars can drive variability in the binary orbits. We also present new eclipse times of NN Ser, which are still compatible with the previously published circumbinary planetary system model, although only with the addition of a quadratic term to the ephemeris. Finally, we conclude that we are limited by the relatively short observational baseline for many of the binaries in the eclipse timing programme, and therefore cannot yet draw robust conclusions about the cause of orbital period variations in evolved, white dwarf binaries.

We report results of the first five observing campaigns for the open stellar cluster NGC 7243 in the frame of project Young Exoplanet Transit Initiative (YETI). The project focuses on the monitoring of young and nearby stellar clusters, with the aim to detect young transiting exoplanets, and to study other variability phenomena on time-scales from minutes to years. After five observing campaigns and additional observations during 2013 and 2014, a clear and repeating transit-like signal was detected in the light curve of J221550.6+495611. Furthermore, we detected and analysed 37 new eclipsing binary stars in the studied region. The best fit parameters and light curves of all systems are given. Finally, we detected and analysed 26 new, presumably pulsating variable stars in the studied region. The follow-up investigation of these objects, including spectroscopic measurements of the exoplanet candidate, is currently planned.

pt5m is a 0.5m robotic telescope located on the roof of the 4.2m William Herschel Telescope (WHT) building, at the Roque de los Muchachos Observatory, La Palma. Using a 5-position filter wheel and CCD detector, and bespoke control software, pt5m provides a high quality robotic observing facility. The telescope first began robotic observing in 2012, and is now contributing to transient follow-up and time-resolved astronomical studies. In this paper we present the scientific motivation behind pt5m, as well as the specifications and unique features of the facility. We also present an example of the science we have performed with pt5m, where we measure the radius of the transiting exoplanet WASP-33b. We find a planetary radius of 1.603 +/- 0.014 R(J).

We report the discovery and characterisation of a deeply eclipsing AM CVn-system, Gaia14aae (= ASSASN-14cn). Gaia14aae was identified independently by the All-Sky Automated Survey for Supernovae (ASAS-SN; Shappee et al. 2014) and by the Gaia Science Alerts project, during two separate outbursts. A third outburst is seen in archival Pan-STARRS-1 (PS1; Schlafly et al. 2012; Tonry et al. 2012; Magnier et al. 2013) and ASAS-SN data. Spectroscopy reveals a hot, hydrogen-deficient spectrum with clear double-peaked emission lines, consistent with an accreting double degenerate classification. We use follow-up photometry to constrain the orbital parameters of the system. We find an orbital period of 49.71 min, which places Gaia14aae at the long period extremum of the outbursting AM CVn period distribution. Gaia14aae is dominated by the light from its accreting white dwarf. Assuming an orbital inclination of 90 degrees for the binary system, the contact phases of the white dwarf lead to lower limits of 0.78 M solar and 0.015 M solar on the masses of the accretor and donor respectively and a lower limit on the mass ratio of 0.019. Gaia14aae is only the third eclipsing AM CVn star known, and the first in which the WD is totally eclipsed. Using a helium WD model, we estimate the accretor's effective temperature to be 12900+-200 K. The three out-burst events occurred within 4 months of each other, while no other outburst activity is seen in the previous 8 years of Catalina Real-time Transient Survey (CRTS; Drake et al. 2009), Pan-STARRS-1 and ASAS-SN data. This suggests that these events might be rebrightenings of the first outburst rather than individual events.

Scintillation noise due to the Earth's turbulent atmosphere can be a dominant noise source in high-precision astronomical photometry when observing bright targets from the ground. Here we describe the phenomenon of scintillation from its physical origins to its effect on photometry. We show that Young's (1967) scintillation-noise approximation used by many astronomers tends to underestimate the median scintillation noise at several major observatories around the world. We show that using median atmospheric optical turbulence profiles, which are now available for most sites, provides a better estimate of the expected scintillation noise and that real-time turbulence profiles can be used to precisely characterise the scintillation noise component of contemporaneous photometric measurements. This will enable a better understanding and calibration of photometric noise sources and the effectiveness of scintillation correction techniques. We also provide new equations for calculating scintillation noise, including for extremely large telescopes where the scintillation noise will actually be lower than previously thought. These equations highlight the fact that scintillation noise and shot noise have the same dependence on exposure time and so if an observation is scintillation limited, it will be scintillation limited for all exposure times. The ratio of scintillation noise to shot noise is also only weakly dependent on telescope diameter and so a bigger telescope may not yield a reduction in fractional scintillation noise.

V471 Tauri, a white dwarf--red dwarf eclipsing binary in the Hyades, is well known for stimulating development of common envelope theory, whereby novae and other cataclysmic variables form from much wider binaries by catastrophic orbit shrinkage. Our evaluation of a recent imaging search that reported negative results for a much postulated third body shows that the object could have escaped detection or may have actually been seen. The balance of evidence continues to favor a brown dwarf companion about 12 AU from the eclipsing binary. A recently developed algorithm finds unified solutions from three datatypes. New radial velocities (RVs) of the red dwarf and BV RCIC light curves are solved simultaneously along with white dwarf and red dwarf RVs from the literature, uvby data, the MOST mission light curve, and 40 years of eclipse timings. Precision-based weighting is the key to proper information balance among the various datasets. Timewise variation of modeled starspots allows unified solution of multiple data eras. Light curve amplitudes strongly suggest decreasing spottedness from 1976 to about 1980, followed by approximately constant spot coverage from 1981 to 2005. An explanation is proposed for lack of noticeable variation in 1981 light curves, in terms of competition between spot and tidal variations. Photometric spectroscopic distance is estimated. The red dwarf mass comes out larger than normal for a K2V star, and even larger than adopted in several structure and evolution papers. An identified cause for this result is that much improved red dwarf RVs curves now exist.

A very large lobe overflow event is suggested to explain the $0.^m4$ brightening observed in K band at pericenter passage of the star known as S2 that orbits the Galaxy's supermassive black hole (SMBH). Known observed properties of S2 that contribute to lobe filling are 1) the enormous mass ratio, $M_{SMBH}/M_{S2}$, 2) S2's fast rotation, and 3) S2's large orbital eccentricity. Published estimates have given limiting lobe sizes of order 100 to 300 $R_\odot$ but, with S2's fast rotation taken into account, the computed lobe size is much smaller, being compatible with either a main sequence OB star or a stripped evolved star. An important evolutionary consideration that predicts very large pericenter overflows is envelope expansion following mass loss that is characteristic of highly evolved stars. Material removed by lobe overflow at pericenter is replenished by envelope expansion as an evolved star awaits its next pericenter passage. An observational signature of lobe overflow for upcoming pericenter passages would be appearance of emission lines as the ejected gas expands and becomes optically thin.

We have produced a multiannual climatology of airborne dust from Martian year 24 to 31 using multiple datasets of retrieved or estimated column optical depths. The datasets are based on observations of the Martian atmosphere from April 1999 to July 2013 made by different orbiting instruments: the Thermal Emission Spectrometer (TES) aboard Mars Global Surveyor, the Thermal Emission Imaging System (THEMIS) aboard Mars Odyssey, and the Mars Climate Sounder (MCS) aboard Mars Reconnaissance Orbiter (MRO). The procedure we have adopted consists of gridding the available retrievals of column dust optical depth (CDOD) from TES and THEMIS nadir observations, as well as the estimates of this quantity from MCS limb observations. Our gridding method calculates averages and uncertainties on a regularly spaced, but possibly incomplete, spatio-temporal grid, using an iterative procedure weighted in space, time, and retrieval uncertainty. In order to evaluate strengths and weaknesses of the resulting gridded maps, we validate them with independent observations of CDOD. We have statistically analyzed the irregularly gridded maps to provide an overview of the dust climatology on Mars over eight years, specifically in relation to its interseasonal and interannual variability. Finally, we have produced multiannual, regular daily maps of CDOD by spatially interpolating the irregularly gridded maps using a kriging method. These synoptic maps are used as dust scenarios in the Mars Climate Database version 5, and are useful in many modelling applications in addition to forming a basis for instrument intercomparisons. The derived dust maps for the eight available Martian years are publicly available and distributed with open access.

The next generation of adaptive optics (AO) systems will require tomographic reconstruction techniques to map the optical refractive index fluctuations, generated by the atmospheric turbulence, along the line of sight to the astronomical target. These systems can be enhanced with data from an external atmospheric profiler. This is important for Extremely Large Telescope scale tomography. Here we propose a new instrument which utilises the generalised SCIntillation Detection And Ranging (SCIDAR) technique to allow high sensitivity vertical profiles of the atmospheric optical turbulence and wind velocity profile above astronomical observatories. The new approach, which we refer to as 'Stereo-SCIDAR', uses a stereoscopic system with the scintillation pattern from each star of a double-star target incident on a separate detector. Separating the pupil images for each star has several advantages including: increased magnitude difference tolerance for the target stars; negating the need for re-calibration due to the normalisation errors usually associated with SCIDAR; an increase of at least a factor of two in the signal-to-noise ratio of the cross-covariance function and hence the profile for equal magnitude target stars and up to a factor of 16 improvement for targets of 3 magnitudes difference; and easier real-time reconstruction of the wind-velocity profile. Theoretical response functions are calculated for the instrument, and the performance is investigated using a Monte-Carlo simulation. The technique is demonstrated using data recorded at the 2.5 m Nordic Optical Telescope and the 1.0 m Jacobus Kapteyn Telescope, both on La Palma.

We present results from the first year of a campaign to characterize and monitor the optical turbulence profile at the SAAO Sutherland observing station in South Africa. A MASS-DIMM (MultiAperture Scintillation Sensor - Differential Image Motion Monitor) was commissioned in March 2010 to provide continuous monitoring of the seeing conditions. Over the first month of the campaign, a SLODAR (SLOpe Detection And Ranging) from Durham University was also installed allowing an independent verification of the performance of the MASS-DIMM device. After the first year of data collection, the overall median seeing value is found to be 1.32" as measured at ground level. The ground layer which includes all layers below 1 km accounts for 84% of the turbulence, while the free atmosphere above 1 km accounts for 16% with a median value of 0.41". The median isoplanatic angle value is 1.92", which is similar to other major astronomical sites. The median coherence time, calculated from corrected MASS measurements, is 2.85 ms. The seeing conditions at the site do show a strong correlation with wind direction with bad seeing conditions being associated with winds from the South-East.

We started a new project which aims to find compact hot subdwarf binaries at low Galactic latitudes. Targets are selected from several photometric surveys and a spectroscopic follow-up campaign to find radial velocity variations on timescales as short as tens of minutes has been started. Once radial variations are detected phase-resolved spectroscopy is obtained to measure the radial velocity curve and the mass function of the system. The observing strategy is described and the discovery of two short period hot subdwarf binaries is presented. UVEXJ032855.25+503529.8 contains a hot subdwarf B star (sdB) orbited by a cool M-dwarf in a P=0.11017 days orbit. The lightcurve shows a strong reflection effect but no eclipses are visible. HS 1741+2133 is a short period (P=0.20 days) sdB most likely with a white dwarf (WD) companion.

We present z'-band secondary eclipse photometry of the highly irradiated hot Jupiter WASP-12b using ULTRACAM on the 4.2m William Herschel Telescope. We measure a decrease in flux of \delta = 0.130+/-0.013% during the passage of the planet behind the star, which is significantly deeper than the previous measurement at this wavelength (0.082+/-0.015%, L\'opez-Morales et al. 2010). Our secondary eclipse is best fit with a mid-eclipse phase, \phi, that is compatible with a circular orbit \phi = 0.501+/-0.002, in agreement with previous results (Croll et al. 2011). In combination with existing data, our eclipse depth measurement allows us to constrain the characteristics of the planet's atmosphere, which is consistent with a carbon-rich model, with no evidence for a strong thermal inversion. If the difference in eclipse depth reported here compared to that of L\'opez-Morales et al. (2010) is of physical origin, as opposed to due to systematics, it may be caused by temporal variability in the flux, due to atmospheric dynamics.

We describe the design and performance of the digital servo and motion control system for the 6-meter diameter parabolic antennas of the Submillimeter Array (SMA) on Mauna Kea, Hawaii. The system is divided into three nested layers operating at a different, appropriate bandwidth. (1) A rack-mounted, real-time Unix system runs the position loop which reads the high resolution azimuth and elevation encoders and sends velocity and acceleration commands at 100 Hz to a custom-designed servo control board (SCB). (2) The microcontroller-based SCB reads the motor axis tachometers and implements the velocity loop by sending torque commands to the motor amplifiers at 558 Hz. (3) The motor amplifiers implement the torque loop by monitoring and sending current to the three-phase brushless drive motors at 20 kHz. The velocity loop uses a traditional proportional-integral-derivative (PID) control algorithm, while the position loop uses only a proportional term and implements a command shaper based on the Gauss error function. Calibration factors and software filters are applied to the tachometer feedback prior to the application of the servo gains in the torque computations. All of these parameters are remotely adjustable in software. The three layers of the control system monitor each other and are capable of shutting down the system safely if a failure or anomaly occurs. The Unix system continuously relays antenna status to the central observatory computer via reflective memory. In each antenna, a Palm Vx hand controller displays system status and allows full local control of the drives in an intuitive touchscreen user interface. It can also be connected outside the cabin for convenience during antenna reconfigurations. Excellent tracking performance (0.3 arcsec rms) is achieved with this system. It has been in reliable operation on 8 antennas for over 10 years and has required minimal maintenance.

We present time-resolved J-band spectroscopy of the short period cataclysmic variable SDSS J143317.78+101123.3. We detect absorption lines from the sub-stellar donor star in this system, which contributes 38 +/- 5% to the J-band light. From the relative strengths of the absorption lines in the J-band, we estimate the spectral type of the donor star to be L2 +/- 1. These data are the first spectroscopic detection of a confirmed sub-stellar donor in a cataclysmic variable, and the spectral type is consistent with that expected from semi-empirical evolutionary models. Using skew mapping, we have been able to derive an estimate for the radial velocity of the donor of Kd = 520 +/- 60 km/s. This value is consistent with, though much less precise than, predictions from mass determinations found via photometric fitting of the eclipse light curves.

We have constructed a five station 12 GHz atmospheric phase interferometer (API) for the Submillimeter Array (SMA) located near the summit of Mauna Kea, Hawaii. Operating at the base of unoccupied SMA antenna pads, each station employs a commercial low noise mixing block coupled to a 0.7 m off-axis satellite dish which receives a broadband, white noise-like signal from a geostationary satellite. The signals are processed by an analog correlator to produce the phase delays between all pairs of stations with projected baselines ranging from 33 to 261 m. Each baseline's amplitude and phase is measured continuously at a rate of 8 kHz, processed, averaged and output at 10 Hz. Further signal processing and data reduction is accomplished with a Linux computer, including the removal of the diurnal motion of the target satellite. The placement of the stations below ground level with an environmental shield combined with the use of low temperature coefficient, buried fiber optic cables provides excellent system stability. The sensitivity in terms of rms path length is 1.3 microns which corresponds to phase deviations of about 1 degree of phase at the highest operating frequency of the SMA. The two primary data products are: (1) standard deviations of observed phase over various time scales, and (2) phase structure functions. These real-time statistical data measured by the API in the direction of the satellite provide an estimate of the phase front distortion experienced by the concurrent SMA astronomical observations. The API data also play an important role, along with the local opacity measurements and weather predictions, in helping to plan the scheduling of science observations on the telescope.

The digital phase lock loop described in this paper has been in use on the Submillimeter Array (SMA) front-end receivers for over a decade and has been a key element in obtaining all of the high quality images that have been published from this telescope over the years. The technical achievements enabled by these devices include the first phase closure observations in the 690 GHz band, the first attempts at band-to-band phase transfer at submillimeter wavelengths, and the first successful demonstration of interferometry using a fully photonic millimeter-wave local oscillator. Copies of these devices are also in place at the Caltech Submillimeter Observatory and the James Clerk Maxwell Telescope in support of the eSMA project and submillimeter VLBI experiments. Additional units of this design were used by the Princeton Millimeter Interferometer and the Microwave Anisotropy Telescope. In total, over three dozen units have been constructed and used in astronomical studies. In this paper, we briefly describe the background theory, design, performance, and calibration steps, and provide useful testing and repair information.

We report a spectral-line survey of the extreme carbon star IRC+10216 carried out between 293.9 and 354.8 GHz with the Submillimeter Array. A total of 442 lines were detected, more than 200 for the first time; 149 are unassigned. Maps at an angular resolution of ~3" were obtained for each line. A substantial new population of narrow lines with an expansion velocity of ~4 km/s (i.e. ~ 30% of the terminal velocity) was detected. Most of these are attributed to rotational transitions within vibrationally excited states, emitted from energy levels above the v=0, J=0 ground state with excitation energy of 1000-3000 K. Emission from these lines appears to be centered on the star with an angular extent of <1". We use multiple transitions detected in several molecules to derive physical conditions in this inner envelope of IRC+10216.

High precision fast photometry from ground-based observatories is a challenge due to intensity fluctuations (scintillation) produced by the Earth's atmosphere. Here we describe a method to reduce the effects of scintillation by a combination of pupil reconjugation and calibration using a comparison star. Because scintillation is produced by high altitude turbulence, the range of angles over which the scintillation is correlated is small and therefore simple correction by a comparison star is normally impossible. We propose reconjugating the telescope pupil to a high dominant layer of turbulence, then apodizing it before calibration with a comparison star. We find by simulation that given a simple atmosphere with a single high altitude turbulent layer and a strong surface layer a reduction in the intensity variance by a factor of ~30 is possible. Given a more realistic atmosphere as measured by SCIDAR at San Pedro M\'artir we find that on a night with a strong high altitude layer we can expect the median variance to be reduced by a factor of 11. By reducing the scintillation noise we will be able to detect much smaller changes in brightness. If we assume a 2 m telescope and an exposure time of 30 seconds a reduction in the scintillation noise from 0.78 mmag to 0.21 mmag is possible, which will enable the routine detection of, for example, the secondary transits of extrasolar planets from the ground.

Using a model of force balance in Saturn's disk-like magnetosphere, we show that variations in hot plasma pressure can change the magnetic field configuration. This effect changes (i) the location of the magnetopause, even at fixed solar wind dynamic pressure, and (ii) the magnetic mapping between ionosphere and disk. The model uses equatorial observations as a boundary condition-we test its predictions over a wide latitude range by comparison with a Cassini high-inclination orbit of magnetic field and hot plasma pressure data. We find reasonable agreement over time scales larger than the period of Saturn kilometric radiation (also known as the camshaft period).

A distance is found for the W UMa type binary U Pegasi, with a newly modified version of the Wilson-Devinney program (W-D) that makes use of the direct distance estimation (DDE) algorithm. The reported distance of d = 123.6pc is an average based on solutions for B and V data and a primary star temperature of 5800K. Standardized light curves (not differential), radial velocities, and a spectroscopic primary star temperature are input to the pro- gram. Differential corrections were performed for each light curve band along with the velocities for two primary temperatures that span 100K. Log10d is a model parameter like many others that are adjustable in W-D. The eclipsing binary distance agrees with the Hipparcos parallax distance and is more precise.

The Gamma-Ray Burst Monitor (GBM) will significantly augment the science return from the Fermi Observatory in the study of Gamma-Ray Bursts (GRBs). The primary objective of GBM is to extend the energy range over which bursts are observed downward from the energy range of the Large Area Telescope (LAT) on Fermi into the hard X-ray range where extensive previous data exist. A secondary objective is to compute burst locations on-board to allow re-orientiong the spacecraft so that the LAT can observe delayed emission from bright bursts. GBM uses an array of twelve sodium iodide scintillators and two bismuth germanate scintillators to detect gamma rays from ~8 keV to ~40 MeV over the full unocculted sky. The on-board trigger threshold is ~0.7 photons/cm2/s (50-300 keV, 1 s peak). GBM generates on-board triggers for ~250 GRBs per year.

As shown by Mennickent, et al(2003), a subset of the blue variable stars in the Large Magellanic Cloud exhibit brightness variability of small amplitude in the period range 2.4 to 16 days as well as larger amplitude variability with periods of 140 to 600 days, with a remarkably tight relation between the long and the short periods. Our re-examination of these objects has led to the discovery of additional variability. The Fourier spectra of 11 of their 30 objects have 3 or 4 peaks above the noise level and a linear relation of the form f_a = 2(f_b - f_L) among three of the frequencies. An explanation of this relation requires an interplay between the binary motion and that of a third object. The two frequency relations together with the Fourier amplitude ratios pose a challenging modeling problem.

Observational astronomy has changed drastically in the last decade: manually driven target-by-target instruments have been replaced by fully automated robotic telescopes. Data acquisition methods have advanced to the point that terabytes of data are flowing in and being stored on a daily basis. At the same time, the vast majority of analysis tools in stellar astrophysics still rely on manual expert interaction. To bridge this gap, we foresee that the next decade will witness a fundamental shift in the approaches to data analysis: case-by-case methods will be replaced by fully automated pipelines that will process the data from their reduction stage, through analysis, to storage. While major effort has been invested in data reduction automation, automated data analysis has mostly been neglected despite the urgent need. Scientific data mining will face serious challenges to identify, understand and eliminate the sources of systematic errors that will arise from this automation. As a special case, we present an artificial intelligence (AI) driven pipeline that is prototyped in the domain of stellar astrophysics (eclipsing binaries in particular), current results and the challenges still ahead.

We present molecular line mapping of the Giant Molecular Cloud G1.6-0.025, which is located at the high longitude end of the Central Molecular Zone of our Galaxy. We assess the degree of star formation activity in that region using several tracers and find very little. We made a large scale, medium (2') resolution map in the J = 2-1 transition of SiO for which we find clumpy emission over a ~0.8 x 0.3 degree-sized region stretching along the Galactic plane. Toward selected positions we also took spectra in the easy to excite J_k=2_k-1_k quartet of CH3OH and the CS 2-1 line. Throughout the cloud these \meth lines are, remarkably, several times stronger than, both, the CS and the SiO lines. The large widths of all the observed lines, similar to values generally found in the Galactic center, indicate a high degree of turbulence. Several high LSR velocity clumps that have 0-80 km/s higher velocities than the bulk of the molecular cloud appear at the same projected position as "normal" velocity material; this may indicate cloud-cloud collisions. Statistical equilibrium modeling of the CH3OH lines observed by us and others yield relatively high densities and moderate temperatures for a representative dual velocity position. We find 8 10^4 cm-3/30 K for material in the G1.6-0.025 cloud and a higher temperature (190 K), but a 50% lower density in a high velocity clump projected on the same location. Several scenarios are discussed in which shock chemistry might enhance the CH3OH and SiO abundances in G1.6-0.025 and elsewhere in the Central Molecular Zone.

Using the Submillimeter Array we have detected the J=3-2 and 2-1 rotational transitions from within the first vibrationally excited state of CO toward the extreme carbon star IRC+10216 (CW Leo). The emission remains spatially unresolved with an angular resolution of ~2" and, given that the lines originate from energy levels that are ~3100 K above the ground state, almost certainly originates from a much smaller (~10^{14} cm) sized region close to the stellar photosphere. Thermal excitation of the lines requires a gas density of ~10^{9} cm^{-3}, about an order of magnitude higher than the expected gas density based previous infrared observations and models of the inner dust shell of IRC+10216.

A spectral-line survey of IRC+10216 in the 345 GHz band has been undertaken with the Submillimeter Array. Although not yet completed, it has already yielded a fairly large sample of narrow molecular emission lines with line-widths indicating expansion velocities of ~4 km/s, less than 3 times the well-known value of the terminal expansion velocity (14.5 km/s) of the outer envelope. Five of these narrow lines have now been identified as rotational transitions in vibrationally excited states of previously detected molecules: the v=1, J=17--16 and J=19--18 lines of Si34S and 29SiS and the v=2, J=7--6 line of CS. Maps of these lines show that the emission is confined to a region within ~60 AU of the star, indicating that the narrow-line emission is probing the region of dust-formation where the stellar wind is still being accelerated.

Observation of the CMB is central to observational cosmology, and the Antarctic Plateau is an exceptionally good site for this work. The first attempt at CMB observations from the Plateau was an expedition to the South Pole in December 1986 by the Radio Physics Research group at Bell Laboratories. Sky noise and opacity were measured. The results were sufficiently encouraging that in the Austral summer of 1988-1989, three CMB groups participated in the "Cucumber" campaign, where a temporary site dedicated to CMB anisotropy measurements was set up 2 km from South Pole Station. Winter-time observations became possible with the establishment in 1990 of the Center for Astrophysical Research in Antarctica (CARA), a National Science Foundation Science and Technology Center. CARA developed year-round observing facilities in the "Dark Sector", a section of Amundsen-Scott South Pole Station dedicated to astronomical observations. CARA scientists fielded several astronomical instruments: AST/RO, SPIREX, White Dish, Python, Viper, ACBAR, and DASI. By 2001, data from CARA, together with BOOMERANG, a CMB experiment on a long-duration balloon launched from McMurdo Station on the coast of Antarctica, showed clear evidence that the overall geometry of the Universe is flat, as opposed to being positively or negatively curved. In 2002, the DASI group reported the detection of polarization in the CMB. These observations strongly support the concordance model of cosmology, where the dynamics of a flat Universe are dominated by forces exerted by the Dark Energy and Dark Matter. The South Pole Telescope (SPT) is a newly-operational 10 m diameter offset telescope designed to rapidly measure anisotropies on scales much smaller than 1 degree.

We report the development and first results of an instrument called Low Layer Scidar (LOLAS) which is aimed at the measurement of optical-turbulence profiles in the atmospheric boundary layer with high altitude-resolution. The method is based on the Generalized Scidar (GS) concept, but unlike the GS instruments which need a 1- m or larger telescope, LOLAS is implemented on a dedicated 40-cm telescope, making it an independent instrument. The system is designed for widely separated double-star targets, which enables the high altitude-resolution. Using a 20000-separation double- star, we have obtained turbulence profiles with unprecedented 12-m resolution. The system incorporates necessary novel algorithms for autoguiding, autofocus and image stabilisation. The results presented here were obtained at Mauna Kea Observatory. They show LOLAS capabilities but cannot be considered as representative of the site. A forthcoming paper will be devoted to the site characterisation. The instrument was built as part of the Ground Layer Turbulence Monitoring Campaign on Mauna Kea for Gemini Observatory.

Temporal fluctuations of the atmospheric piston are critical for interferometers as they determine their sensitivity. We characterize an instrumental set-up, termed the piston scope, that aims at measuring the atmospheric time constant, tau0, through the image motion in the focal plane of a Fizeau interferometer. High-resolution piston scope measurements have been obtained at the observatory of Paranal, Chile, in April 2006. The derived atmospheric parameters are shown to be consistent with data from the astronomical site monitor, provided that the atmospheric turbulence is displaced along a single direction. Piston scope measurements, of lower temporal and spatial resolution, were for the first time recorded in February 2005 at the Antarctic site of DomeC. Their re-analysis in terms of the new data calibration sharpens the conclusions of a first qualitative examination.

Adaptive optics systems are essential on all large telescopes where image quality is important. These are complex systems with many design parameters requiring optimisation before good performance can be achieved. The simulation of adaptive optics systems is therefore necessary to categorise the expected performance. This paper describes an adaptive optics simulation platform, developed at Durham University, which can be used to simulate adaptive optics systems on the largest proposed future extremely large telescopes (ELTs) as well as current systems. This platform is modular, object oriented and has the benefit of hardware application acceleration which can be used to improve the simulation performance, essential for ensuring that the run time of a given simulation is acceptable. The simulation platform described here can be highly parallelised using parallelisation techniques suited for adaptive optics simulation, whilst still offering the user complete control while the simulation is running. Results from the simulation of a ground layer adaptive optics system are provided as an example to demonstrate the flexibility of this simulation platform.

Using five independent analytic and Monte Carlo simulation codes, we have studied the performance of wide field ground layer adaptive optics (GLAO), which can use a single, relatively low order deformable mirror to correct the wavefront errors from the lowest altitude turbulence. GLAO concentrates more light from a point source in a smaller area on the science detector, but unlike traditional adaptive optics, images do not become diffraction-limited. Rather the GLAO point spread function (PSF) has the same functional form as a seeing-limited PSF, and can be characterized by familiar performance metrics such as Full-Width Half-Max (FWHM). The FWHM of a GLAO PSF is reduced by 0.1" or more for optical and near-infrared wavelengths over different atmospheric conditions. For the Cerro Pachon atmospheric model this correction is even greater when the image quality is worst, which effectively eliminates "bad-seeing" nights; the best seeing-limited image quality, available only 20% of the time, can be achieved 60 to 80% of the time with GLAO. This concentration of energy in the PSF will reduce required exposure times and improve the efficiency of an observatory up to 30 to 40%. These performance gains are relatively insensitive to a number of trades including the exact field of view of a wide field GLAO system, the conjugate altitude and actuator density of the deformable mirror, and the number and configuration of the guide stars.

We report the measurement results and compensation of the antenna elevation angle dependences of the Sub-millimeter Array (SMA) antenna characteristics. Without optimizing the subreflector (focus) positions as a function of the antenna elevation angle, antenna beam patterns show lopsided sidelobes, and antenna efficiencies show degradations. The sidelobe level increases and the antenna efficiencies decrease about 1% and a few %, respectively, for every 10 degrees change in the elevation angle at the measured frequency of 237 GHz. We therefore obtained the optimized subreflector positions for X (azimuth), Y (elevation), and Z (radio optics) focus axes at various elevation angles for all the eight SMA antennas. The X axis position does not depend on the elevation angle. The Y and Z axes positions depend on the elevation angles, and are well fitted with a simple function for each axis with including a gravity term (cosine and sine of elevation, respectively). In the optimized subreflector positions, the antenna beam patterns show low level symmetric sidelobe of at most a few %, and the antenna efficiencies stay constant at any antenna elevation angles. Using one set of fitted functions for all antennas, the SMA is now operating with real-time focusing, and showing constant antenna characteristics at any given elevation angle.

A method for producing a laser guide star wavefront sensor for adaptive optics with reduced focal anisoplanatism is presented. A theoretical analysis and numerical simulations have been carried out and the results are presented. The technique, named SPLASH (Sky-Projected Laser Array Shack-Hartmann), is shown to suffer considerably less from focal anisoplanatism than a conventional laser guide star system. The method is potentially suitable for large telescope apertures (~8m), and possibly for extremely large telescopes.

Numerical Simulation is an essential part of the design and optimisation of astronomical adaptive optics systems. Simulations of adaptive optics are computationally expensive and the problem scales rapidly with telescope aperture size, as the required spatial order of the correcting system increases. Practical realistic simulations of AO systems for extremely large telescopes are beyond the capabilities of all but the largest of modern parallel supercomputers. Here we describe a more cost effective approach through the use of hardware acceleration using field programmable gate arrays. By transferring key parts of the simulation into programmable logic, large increases in computational bandwidth can be expected. We show that the calculation of wavefront sensor image centroids can be accelerated by a factor of four by transferring the algorithm into hardware. Implementing more demanding parts of the adaptive optics simulation in hardware will lead to much greater performance improvements, of up to 1000 times.

We report here on serendipitous observations of the intense gamma-ray flare from SGR 1806-20 that occured on 27 December 2004. Unique data from the Cluster and Double Star-2 satellites, designed to study the Earth's magnetosphere, provide the first observational evidence of three separate timescales within the early (first 100ms) phases of this class of events. These observations reveal that, in addition to the initial very steep (<0.25ms) X-ray onset, there is firstly a 4.9ms exponential rise timescale followed by a continued exponential rise in intensity on a timescale of 70ms. These three timescales are a prominent feature of current theoretical models including the timescale (several ms) for fracture propagation in the crust of the neutron star.

The next large NASA mission in the field of gamma-ray astronomy, GLAST, is scheduled for launch in 2007. Aside from the main instrument LAT (Large-Area Telescope), a gamma-ray telescope for the energy range between 20 MeV and > 100 GeV, a secondary instrument, the GLAST burst monitor (GBM), is foreseen. With this monitor one of the key scientific objectives of the mission, the determination of the high-energy behaviour of gamma-ray bursts and transients can be ensured. Its task is to increase the detection rate of gamma-ray bursts for the LAT and to extend the energy range to lower energies (from ~10 keV to \~30 MeV). It will provide real-time burst locations over a wide FoV with sufficient accuracy to allow repointing the GLAST spacecraft. Time-resolved spectra of many bursts recorded with LAT and the burst monitor will allow the investigation of the relation between the keV and the MeV-GeV emission from GRBs over unprecedented seven decades of energy. This will help to advance our understanding of the mechanisms by which gamma-rays are generated in gamma-ray bursts.

A rationale is presented for the use of a relatively low-altitude Rayleigh Laser Guide Star to provide partial adaptive optics correction across a large fraction of the sky on the 4.2m William Herschel Telescope. The scientific motivation is highlighted and supported by model calculations. An overview the technical implementation of the system is presented.

Published measurements of distances to detached eclipsing binaries in the Large Magellanic Cloud combine stellar surface areas (in absolute units) determined from photometric light and radial velocity curves with surface brightnesses from model atmospheres and observed spectra. The method does not require the stars to be normal or undistorted, and is not limited in its applicability to the well detached systems that have traditionally been considered. With this in mind we carry out semi-detached light curve solutions for SMC binaries discovered by the OGLE collaboration, identify candidates for distance estimation, and tabulate OGLE light curve solutions. We point out important advantages for semi-detached binaries as standard candles; (1) light curve solutions can be strengthened by exploiting lobe-filling configurations, (2) only single-lined spectra may be needed because the mass ratio can often be determined from photometry, and (3) nearly all semi-detached binaries have sensibly circular orbits. We carry out simulations with synthetic data to show that semi-detached binaries can often be reliably identified and to quantify the accuracy of solutions. The simulations demonstrate two additional advantages for semi-detached distance determination candidates; (4) the well-known difficulty in distinguishing solutions with interchanged radii (aliasing) is much less severe for semi-detached than for detached binaries, and (5) the desirable condition of complete eclipse is identified with improved reliability. In many cases we find that detached and semi-detached systems can be distinguished, and we select 36 candidate semi-detached systems from the OGLE SMC eclipsing binary catalog. We also note the near absence of Algol-like light curves in the OGLE sample and discuss possible explanations.

Detached eclipsing binary stars provide a robust one-step distance determination to nearby galaxies. As a by-product of Galactic microlensing searches, catalogs of thousands of variable stars including eclipsing binaries have been produced by the OGLE, MACHO and EROS collaborations. We present photometric solutions for detached eclipsing binaries in the Small Magellanic Cloud (SMC) discovered by the OGLE collaboration. The solutions were obtained with an automated version of the Wilson-Devinney program. By fitting mock catalogs of eclipsing binaries we find that the normalized stellar radii (particularly their sum) and the surface brightness ratio are accurately described by the fitted parameters and estimated standard errors, despite various systematic uncertainties. In many cases these parameters are well constrained. In addition we find that systems exhibiting complete eclipses can be reliably identified where the fractional standard errors in the radii are small. We present two quantitatively selected sub-samples of eclipsing binaries that will be excellent distance indicators. These can be used both for computation of the distance to the SMC and to probe its structure. One particularly interesting binary has a very well determined solution, exhibits complete eclipses, and is comprised of well detached G-type, class $II$ giants.

A quantitative assessment is presented of diffraction-limited stellar images with Strehl ratios of 0.25-0.30 obtained by selection of short-exposure CCD images of stars brighter than +6m at 810nm with the Nordic Optical Telescope.

AST/RO, a 1.7 m diameter telescope for astronomy and aeronomy studies at wavelengths between 200 and 2000 microns, was installed at the South Pole during the 1994-1995 Austral summer. The telescope operates continuously through the Austral winter, and is being used primarily for spectroscopic studies of neutral atomic carbon and carbon monoxide in the interstellar medium of the Milky Way and the Magellanic Clouds. The South Pole environment is unique among observatory sites for unusually low wind speeds, low absolute humidity, and the consistent clarity of the submillimeter sky. Four heterodyne receivers, an array receiver, three acousto-optical spectrometers, and an array spectrometer are installed. A Fabry-Perot spectrometer using a bolometric array and a Terahertz receiver are in development. Telescope pointing, focus, and calibration methods as well as the unique working environment and logistical requirements of the South Pole are described.

Ginga and Rossi X-ray Timing Explorer (RXTE) observations have allowed an unprecedented view of the recurrent systematic pulse shape changes associated with the 35-day cycle of Hercules X-1, a phenomenon currently unique among the known accretion-powered pulsars. We present observations of the pulse shape evolution. An explanation for the pulse evolution in terms of a freely precessing neutron star is reviewed and shown to have several major difficulties in explaining the observed pulse evolution pattern. Instead, we propose a phenomenological model for the pulse evolution based upon an occultation of the pulse emitting region by the tilted, inner edge of a precessing accretion disk. The systematic and repeating pulse shape changes require a resolved occultation of the pulse emission region. The observed pulse profile motivates the need for a pulsar beam consisting of a composite coaxial pencil and fan beam but the observed evolution pattern requires the fan beam to be focused around the neutron star and beamed in the antipodal direction. The spectral hardness of the pencil beam component suggests an origin at the magnetic polar cap, with the relatively softer fan beam emission produced by backscattering from within the accretion column, qualitatively consistent with several theoretical models for X-ray emission from the accretion column of an accreting neutron star.

We present observations of a series of 10 outbursts of pulsed hard X-ray flux from the transient 10.6 mHz accreting pulsar GS 1843-02, using the Burst and Transient Source Experiment on the Compton Gamma Ray Observatory. These outbursts occurred regularly every 242 days, coincident with the ephemeris of the periodic transient GRO J1849-03 (Zhang et al. 1996), which has recently been identified with the SAS 3 source 2S 1845-024 (Soffitta et al. 1998). Our pulsed detection provides the first clear identification of GS 1843-02 with 2S 1845-024. We present a pulse timing analysis which shows that the 2S 1845-024 outbursts occur near the periastron passage of the neutron star's highly eccentric (e = 0.88+-0.01) 242.18+-0.01 day period binary orbit about a high mass (M > 7 solar masses) companion. The orbit and transient outburst pattern strongly suggest the pulsar is in a binary system with a Be star. Our observations show a long-term spin-up trend, with most of the spin-up occurring during the outbursts. From the measured spin-up rates and inferred luminosities we conclude that an accretion disk is present during the outbursts.

Dramatic torque reversals between spin up and spin down have been observed in half of the persistent X-ray pulsars monitored by the BATSE all-sky monitor on CGRO. Theoretical models developed to explain early pulsar timing data can explain spin down torques via a disk-magnetosphere interaction if the star nearly corotates with the inner accretion disk. To produce the observed BATSE torque reversals, however, these equilibrium models require the disk to alternate between two mass accretion rates, with $\dot M_{\pm}$ producing accretion torques of similar magnitude, but always of opposite sign. Moreover, in at least one pulsar (GX 1+4) undergoing secular spin down the neutron star spins down faster during brief ($\sim 20$ day) hard X-ray flares -- this is opposite the correlation expected from standard theory, assuming BATSE pulsed flux increases with mass accretion rate. The $10$ day to 10 yr intervals between torque reversals in these systems are much longer than any characteristic magnetic or viscous time scale near the inner disk boundary and are more suggestive of a global disk phenomenon. We discuss possible explanations of the observed torque behavior. Despite the preferred sense of rotation defined by the binary orbit, the BATSE observations are surprisingly consistent with an earlier suggestion by Makishima \etal (1988) for GX~1+4: the disks in these systems somehow alternate between episodes of prograde and retrograde rotation. We are unaware of any mechanism that could produce a stable retrograde disk in a binary undergoing Roche-lobe overflow, but such flip-flop behavior does occur in numerical simulations of wind-fed systems. One possibility is that the disks in some of these binaries are fed by an X-ray excited wind.

We summarize five years of continuous monitoring of accretion-powered pulsars with the Burst and Transient Source Experiment (BATSE) on the Compton Gamma Ray Observatory. Our 20-70 keV observations have determined or refined the orbital parameters of 13 binaries, discovered 5 new transient accreting pulsars, measured the pulsed flux history during outbursts of 12 transients (GRO J1744-28, 4U 0115+634, GRO J1750-27, GS 0834-430, 2S 1417-624, GRO J1948+32, EXO 2030+375, GRO J1008-57, A 0535+26, GRO J2058+42, 4U 1145-619 and A 1118-616), and also measured the accretion torque history of during outbursts of 6 of those transients whose orbital parameters were also known. We have also continuously measured the pulsed flux and spin frequency for eight persistently accreting pulsars (Her X-1, Cen X-3, Vela X-1, OAO 1657-415, GX 301-2, 4U 1626-67, 4U 1538-52, and GX 1+4). Because of their continuity and uniformity over a long baseline, BATSE observations have provided new insights into the long-term behavior of accreting magnetic stars. We have found that all accreting pulsars show stochastic variations in their spin frequencies and luminosities, including those displaying secular spin-up or spin-down on long time scales, blurring the conventional distinction between disk-fed and wind-fed binaries. Pulsed flux and accretion torque are strongly correlated in outbursts of transient accreting pulsars, but uncorrelated, or even anticorrelated, in persistent sources.

Over five years of daily hard X-ray (>20 keV) monitoring of the 2-min accretion-powered pulsar GX 1+4 with the Compton Gamma Ray Observatory/BATSE large-area detectors has found nearly continuous rapid spin-down, interrupted by a bright 200-d spin-up episode. During spin-down, the torque becomes more negative as the luminosity increases (assuming that the 20-60 keV pulsed flux traces bolometric luminosity), the opposite of what is predicted by standard accretion torque theory. No changes in the shape of the 20-100 keV pulsed energy spectrum were detected, so that a very drastic change in the spectrum below 20 keV or the pulsed fraction would be required to make the 20-60 keV pulsed flux a poor luminosity tracer. These are the first observations which flatly contradict standard magnetic disk accretion theory, and they may have important implications for understanding the spin evolution of X-ray binaries, cataclysmic variables, and protostars. We briefly discuss the possibility that GX 1+4 may be accreting from a retrograde disk during spin-down, as previously suggested.

We report the results of a CGRO 3-week observation of the binary system containing the 47 ms pulsar PSR B1259-63 orbiting around a Be star companion in a very eccentric orbit. The PSR B1259-63 binary is a unique system for the study of the interaction of a rapidly rotating pulsar with time-variable nebular surroundings. CGRO observed the PSR B1259-63 system in coincidence with its most recent periastron passage (January 3-23, 1994). Unpulsed and non-thermal hard X-ray emission was detected up to 200 keV, with a photon index $1.8 \pm 0.2$ and a flux of ~4 mCrab, corresponding to a luminosity of a few 10^{34} erg/s at the distance of 2 kpc. The hard X-ray flux and spectrum detected by CGRO agrees with the X-ray emission measured by simultaneous ASCA observations. EGRET upper limits are significant, and exclude strong inverse Compton cooling in the PSR B1259-63 system. We interpret the observed non-thermal emission as synchrotron radiation of shocked electron/positron pairs of the relativistic pulsar wind interacting with the mass outflow from the Be star. Our results clearly indicate, for the first time in a binary pulsar, that high energy emission can be shock-powered rather than caused by accretion. The lack of X-ray/gamma-ray pulsations constrains models of high-energy emission from rapidly rotating pulsars.

The cyclotron line in the spectrum of the accretion-powered pulsar Her X-1 offers an opportunity to assess the ability of the BATSE Spectroscopy Detectors (SDs) to detect lines like those seen in some GRBs. Preliminary analysis of an initial SD pulsar mode observation of Her X-1 indicated a cyclotron line at an energy of approximately 44 keV, rather than at the expected energy of approximately 36 keV. Our analysis of four SD pulsar mode observations of Her X-1 made during high-states of its 35 day cycle confirms this result. We consider a number of phenomenological models for the continuum spectrum and the cyclotron line. This ensures that we use the simplest models that adequately describe the data, and that our results are robust. We find modest evidence (significance Q ~ 10^-4-10^-2) for a line at approximately 44 keV in the data of the first observation. Joint fits to the four observations provide stronger evidence (Q ~ 10^-7-10^-4) for the line. Such a shift in the cyclotron line energy of an accretion-powered pulsar is unprecedented.

We investigate spectral evolution in 37 bright, long gamma-ray bursts observed with the BATSE Spectroscopy Detectors. High resolution spectra are characterized by the energy of the peak of \nfn~and the evolution of this quantity is examined relative to the emission intensity. In most cases it is found that this peak energy either rises with or slightly precedes major intensity increases and softens for the remainder of the pulse. Inter-pulse emission is generally harder early in the burst. For bursts with multiple intensity pulses, later spikes tend to be softer than earlier ones indicating that the energy of the peak of \nfn~is bounded by an envelope which decays with time. Evidence is found that bursts in which the bulk of the flux comes well after the event which triggers the instrument tend to show less peak energy variability and are not as hard as several bursts in which the emission occurs promptly after the trigger. Several recently proposed burst models are examined in light of these results and no qualitative conflicts with the observations presented here are found.