We implement the adaptive step size scheme from the optimization methods AdaGrad and Adam in a novel variant of the Proximal Gradient Method (PGM). Our algorithm, dubbed AdaProx, avoids the need for explicit computation of the Lipschitz constants or additional line searches and thus reduces per-iteration cost. In test cases for Constrained Matrix Factorization we demonstrate the advantages of AdaProx in fidelity and performance over PGM, while still allowing for arbitrary penalty functions. The python implementation of the algorithm presented here is available as an open-source package at https://github.com/pmelchior/proxmin.

We present a Bayesian machine learning architecture that combines a physically motivated parametrization and an analytic error model for the likelihood with a deep generative model providing a powerful data-driven prior for complex signals. This combination yields an interpretable and differentiable generative model, allows the incorporation of prior knowledge, and can be utilized for observations with different data quality without having to retrain the deep network. We demonstrate our approach with an example of astronomical source separation in current imaging data, yielding a physical and interpretable model of astronomical scenes.

The Large Synoptic Survey Telescope (LSST) is an ambitious astronomical survey with a similarly ambitious Data Management component. Data Management for LSST includes processing on both nightly and yearly cadences to generate transient alerts, deep catalogs of the static sky, and forced photometry light-curves for billions of objects at hundreds of epochs, spanning at least a decade. The algorithms running in these pipelines are individually sophisticated and interact in subtle ways. This paper provides an overview of those pipelines, focusing more on those interactions than the details of any individual algorithm.

We present the source separation framework SCARLET for multi-band images, which is based on a generalization of the Non-negative Matrix Factorization to alternative and several simultaneous constraints. Our approach describes the observed scene as a mixture of components with compact spatial support and uniform spectra over their support. We present the algorithm to perform the matrix factorization and introduce constraints that are useful for optical images of stars and distinct stellar populations in galaxies, in particular symmetry and monotonicity with respect to the source peak position. We also derive the treatment of correlated noise and convolutions with band-dependent point spread functions, rendering our approach applicable to coadded images observed under variable seeing conditions. SCARLET thus yields a PSF-matched photometry measurement with an optimally chosen weight function given by the mean morphology in all available bands. We demonstrate the performance of SCARLET for deblending crowded extragalactic scenes and on an AGN jet -- host galaxy separation problem in deep 5-band imaging from the Hyper Suprime-Cam Stategic Survey Program. Using simulations with prominent crowding we show that SCARLET yields superior results to the HSC-SDSS deblender for the recovery of total fluxes, colors, and morphologies. Due to its non-parametric nature, a conceptual limitation of SCARLET is its sensitivity to undetected sources or multiple stellar population within detected sources, but an iterative strategy that adds components at the location of significant residuals appears promising. The code is implemented in Python with C++ extensions and is available at https://github.com/fred3m/scarlet

Using images taken with the Dark Energy Camera (DECam), the first extensive survey of low mass and substellar objects is made in the 15-20 Myr Upper Centaurus Lupus (UCL) and Lower Centaurus Crux (LCC) subgroups of the Scorpius Centaurus OB Association (Sco-Cen). Due to the size of our dataset (>2Tb) we developed an extensive open source set of python libraries to reduce our images, including astrometry, coaddition, and PSF photometry. Our survey consists of 29$\times$3 deg$^2$ fields in the UCL and LCC subgroups of Sco-Cen and the creation of a catalog with over 11 million point sources. We create a prioritized list of candidate for members in UCL and LCC, with 118 \emph{best} and another 348 \emph{good} candidates. We show that the luminosity and mass functions of our low mass and substellar candidates are consistent with measurements for the younger Upper Scorpius subgroup and estimates of a universal IMF, with spectral types ranging from M1 down to L1.

As the size of images and data products derived from astronomical data continues to increase, new tools are needed to visualize and interact with that data in a meaningful way. Motivated by our own astronomical images taken with the Dark Energy Camera (DECam) we present Toyz, an open source Python package for viewing and analyzing images and data stored on a remote server or cluster. Users connect to the Toyz web application via a web browser, making it an convenient tool for students to visualize and interact with astronomical data without having to install any software on their local machines. In addition it provides researchers with an easy-to-use tool that allows them to browse the files on a server and quickly view very large images ($>$ 2 Gb) taken with DECam and other cameras with a large FOV and create their own visualization tools that can be added on as extensions to the default Toyz framework.

We investigate the nature of the unusual eclipsing star OGLE LMC-ECL-11893 (OGLE J05172127-6900558) in the Large Magellanic Cloud recently reported by Dong et al. 2014. The eclipse period for this star is 468 days, and the eclipses exhibit a minimum of ~1.4 mag, preceded by a plateau of ~0.8 mag. Spectra and optical/IR photometry are consistent with the eclipsed star being a lightly reddened B9III star of inferred age ~150 Myr and mass of ~4 solar masses. The disk appears to have an outer radius of ~0.2 AU with predicted temperatures of ~1100-1400 K. We model the eclipses as being due to either a transiting geometrically thin dust disk or gaseous accretion disk around a secondary object; the debris disk produces a better fit. We speculate on the origin of such a dense circumstellar dust disk structure orbiting a relatively old low-mass companion, and on the similarities of this system to the previously discovered EE Cep.

<abridged> The large relative sizes of circumstellar and circumplanetary disks imply that they might be seen in eclipse in stellar light curves. We estimate that a survey of ~10^4 young (~10 Myr old) post-accretion pre-MS stars monitored for ~10 years should yield at least a few deep eclipses from circumplanetary disks and disks surrounding low mass companion stars. We present photometric and spectroscopic data for a pre-MS K5 star (1SWASP J140747.93-394542.6), a newly discovered ~0.9 Msun member of the ~16 Myr-old Upper Cen-Lup subgroup of Sco-Cen at a kinematic distance of 128 pc. SuperWASP and ASAS light curves for this star show a remarkably long, deep, and complex eclipse event centered on 29 April 2007. At least 5 multi-day dimming events of >0.5 mag are identified, with a >3.3 mag deep eclipse bracketed by two pairs of ~1 mag eclipses symmetrically occurring +-12 days and +-26 days before and after. Hence, significant dimming of the star was taking place on and off over at least a ~54 day period in 2007, and a strong >1 mag dimming event occurred over a ~12 day span. We place a firm lower limit on the period of 850 days (i.e. the orbital radius of the eclipser must be >1.7 AU and orbital velocity must be <22 km/s). The shape of the light curve is similar to the lop-sided eclipses of the Be star EE Cep. We suspect that this new star is being eclipsed by a low-mass object orbited by a dense inner disk, girded by at least 3 dusty rings of lower optical depth. Between these rings are at least two annuli of near-zero optical depth (i.e. gaps), possibly cleared out by planets or moons, depending on the nature of the secondary. For possible periods in the range 2.33-200 yr, the estimated total ring mass is ~8-0.4 Mmoon (if the rings have optical opacity similar to Saturn's rings), and the edge of the outermost detected ring has orbital radius ~0.4-0.09 AU.