PHOEBE 2 is a Python package for modeling the observables of eclipsing star systems, but until now has focused entirely on the forward-model -- that is, generating a synthetic model given fixed values of a large number of parameters describing the system and the observations. The inverse problem, obtaining orbital and stellar parameters given observational data, is more complicated and computationally expensive as it requires generating a large set of forward-models to determine which set of parameters and uncertainties best represent the available observational data. The process of determining the best solution and also of obtaining reliable and robust uncertainties on those parameters often requires the use of multiple algorithms, including both optimizers and samplers. Furthermore, the forward-model of PHOEBE has been designed to be as physically robust as possible, but is computationally expensive compared to other codes. It is useful, therefore, to use whichever code is most efficient given the reasonable assumptions for a specific system, but learning the intricacies of multiple codes presents a barrier to doing this in practice. Here we present the 2.3 release of PHOEBE (publicly available from http://phoebe-project.org) which introduces a general framework for defining and handling distributions on parameters, and utilizing multiple different estimation, optimization, and sampling algorithms. The presented framework supports multiple forward-models, including the robust model built into PHOEBE itself.

Traditionally, the effects of interstellar extinction on binary star light curves have been treated as a uniform reduction in the observed brightness of the system that is independent of orbital phase. However, unless the orbital plane of the system coincides with the plane of the sky, or if the two stars are completely identical and present with minimal mutual irradiation and tidal/rotational distortions, then this is unlikely to be an accurate representation of the effect of interstellar extinction. Here, we present an updated treatment of interstellar extinction as incorporated in the PHOEBE 2.2 release (publicly available from http://phoebe-project.org) and assess the importance of using such an approach in the modeling of different types of binary systems. We also present the incorporation of PHOENIX model atmospheres into the PHOEBE 2.2 release, providing increased fidelity on computed observables down to lower temperatures than previously available. The importance of these new code developments is then highlighted via an extincted toy model of the eclipsing white-dwarf-subdwarf binary SDSS~J235524.29+044855.7 -- demonstrating that, in the age of LSST as well as complementary space-based photometric missions, a proper accounting for extinction and as well as the use of realistic model atmospheres will be essential in deriving accurate binary parameters.

Binary systems where the axis of rotation (spin) of one or both components is tilted w.r.t. the axis of revolution are called misaligned binary systems. The angle of misalignment, obliquity, has been measured for a handful of stars and extrasolar planets to date. Here we present a mathematical framework for a complete and rigorous treatment of misalignment and introduce an extension to the public PHOEBE code that implements this framework. We discuss misalignment for the Roche geometry and introduce methods for computing stellar shapes, equilibrium (generalized Lagrange) points of the potential and minimal requirements for lobe existence. Efficient parametrization of misalignment is proposed in the plane-of-sky coordinates and implementation details in PHOEBE are given alongside the proof-of-concept toy model, comparison with a known misaligned binary DI Her, and comparison with a misaligned planetary system Kepler-13. We provide important mathematical details of the model in the Appendix. This paper accompanies the release of PHOEBE 2.1, which will be available soon from its website http://phoebe-project.org.

The precision of photometric and spectroscopic observations has been systematically improved in the last decade, mostly thanks to space-borne photometric missions and ground-based spectrographs dedicated to finding exoplanets. The field of eclipsing binary stars strongly benefited from this development. Eclipsing binaries serve as critical tools for determining fundamental stellar properties (masses, radii, temperatures and luminosities), yet the models are not capable of reproducing observed data well either because of the missing physics or because of insufficient precision. This led to a predicament where radiative and dynamical effects, insofar buried in noise, started showing up routinely in the data, but were not accounted for in the models. PHOEBE (PHysics Of Eclipsing BinariEs; http://phoebe-project.org) is an open source modeling code for computing theoretical light and radial velocity curves that addresses both problems by incorporating missing physics and by increasing the computational fidelity. In particular, we discuss triangulation as a superior surface discretization algorithm, meshing of rotating single stars, light time travel effect, advanced phase computation, volume conservation in eccentric orbits, and improved computation of local intensity across the stellar surfaces that includes photon-weighted mode, enhanced limb darkening treatment, better reflection treatment and Doppler boosting. Here we present the concepts on which PHOEBE is built on and proofs of concept that demonstrate the increased model fidelity.

We have analyzed the first 3.75 years of data from TAOS, the Taiwanese American Occultation Survey. TAOS monitors bright stars to search for occultations by Kuiper Belt Objects (KBOs). This dataset comprises 5e5 star-hours of multi-telescope photometric data taken at 4 or 5 Hz. No events consistent with KBO occultations were found in this dataset. We compute the number of events expected for the Kuiper Belt formation and evolution models of Pan & Sari (2005), Kenyon & Bromley (2004), Benavidez & Campo Bagatin (2009), and Fraser (2009). A comparison with the upper limits we derive from our data constrains the parameter space of these models. This is the first detailed comparison of models of the KBO size distribution with data from an occultation survey. Our results suggest that the KBO population is comprised of objects with low internal strength and that planetary migration played a role in the shaping of the size distribution.

The Taiwanese-American Occultation Survey (TAOS) operates four fully automatic telescopes to search for occultations of stars by Kuiper Belt Objects. It is a versatile facility that is also useful for the study of initial optical GRB afterglows. This paper provides a detailed description of the TAOS multi-telescope system, control software, and high-speed imaging.

Achieving maximum scientific results from the overwhelming volume of astronomical data to be acquired over the next few decades will demand novel, fully automatic methods of data analysis. Artificial intelligence approaches hold great promise in contributing to this goal. Here we apply neural network learning technology to the specific domain of eclipsing binary (EB) stars, of which only some hundreds have been rigorously analyzed, but whose numbers will reach millions in a decade. Well-analyzed EBs are a prime source of astrophysical information whose growth rate is at present limited by the need for human interaction with each EB data-set, principally in determining a starting solution for subsequent rigorous analysis. We describe the artificial neural network (ANN) approach which is able to surmount this human bottleneck and permit EB-based astrophysical information to keep pace with future data rates. The ANN, following training on a sample of 33,235 model light curves, outputs a set of approximate model parameters (T2/T1, (R1+R2)/a, e sin(omega), e cos(omega), and sin i) for each input light curve data-set. The whole sample is processed in just a few seconds on a single 2GHz CPU. The obtained parameters can then be readily passed to sophisticated modeling engines. We also describe a novel method polyfit for pre-processing observational light curves before inputting their data to the ANN and present the results and analysis of testing the approach on synthetic data and on real data including fifty binaries from the Catalog and Atlas of Eclipsing Binaries (CALEB) database and 2580 light curves from OGLE survey data. [abridged]

The serendipitous detection of stellar occultations by Outer Solar System objects is a powerful method for ascertaining the small end ($r \lesssim 15$ km) of the size distribution of Kuiper Belt Objects and may potentially allow the exploration of objects as far out as the Oort Cloud. The design and implementation of an occultation survey is aided by a detailed understanding of how diffraction and observational parameters affect the detection of occultation events. In this study, stellar occultations are simulated, accounting for diffraction effects, finite source sizes, finite bandwidths, stellar spectra, sampling, and signal-to-noise. Finally, the possibility of detecting small Outer Solar System objects from the Kuiper Belt all the way out to the Oort Cloud is explored for three photometric systems: a proposed space telescope, Whipple (Kaplan et al. 2003), the Taiwanese-American Occultation Survey (Lehner et al. 2006), and the Multi Mirror Telescope (Bianco 2007).

The Taiwan-America Occultation Survey (TAOS) aims to determine the number of small icy bodies in the outer reach of the Solar System by means of stellar occultation. An array of 4 robotic small (D=0.5 m), wide-field (f/1.9) telescopes have been installed at Lulin Observatory in Taiwan to simultaneously monitor some thousand of stars for such rare occultation events. Because a typical occultation event by a TNO a few km across will last for only a fraction of a second, fast photometry is necessary. A special CCD readout scheme has been devised to allow for stellar photometry taken a few times per second. Effective analysis pipelines have been developed to process stellar light curves and to correlate any possible flux changes among all telescopes. A few billion photometric measurements have been collected since the routine survey began in early 2005. Our preliminary result of a very low detection rate suggests a deficit of small TNOs down to a few km size, consistent with the extrapolation of some recent studies of larger (30--100 km) TNOs.

We present a methodology to discover outliers in catalogs of periodic light-curves. We use cross-correlation as measure of ``similarity'' between two individual light-curves and then classify light-curves with lowest average ``similarity'' as outliers. We performed the analysis on catalogs of variable stars of known type from the MACHO and OGLE projects and established that our method correctly identifies light-curves that do not belong to those catalogs as outliers. We show how our method can scale to large datasets that will be available in the near future such as those anticipated from Pan-STARRS and LSST.