Astronomical catalogs derived from wide-field imaging surveys are an important tool for understanding the Universe. We construct an astronomical catalog from 55 TB of imaging data using Celeste, a Bayesian variational inference code written entirely in the high-productivity programming language Julia. Using over 1.3 million threads on 650,000 Intel Xeon Phi cores of the Cori Phase II supercomputer, Celeste achieves a peak rate of 1.54 DP PFLOP/s. Celeste is able to jointly optimize parameters for 188M stars and galaxies, loading and processing 178 TB across 8192 nodes in 14.6 minutes. To achieve this, Celeste exploits parallelism at multiple levels (cluster, node, and thread) and accelerates I/O through Cori's Burst Buffer. Julia's native performance enables Celeste to employ high-level constructs without resorting to hand-written or generated low-level code (C/C++/Fortran), and yet achieve petascale performance.

A low-background, high-purity germanium detector has been used to search for evidence of low-energy, bremsstrahlung-generated solar axions. An upper bound of $1.36\times 10^{-11}$ $(95% CL)$ is placed on the direct coupling of DFSZ model axions to electrons. The prospects for the sensitivity of the Majorana Demonstrator array of point-contact germanium detectors to solar axions are discussed in the context of the model-independent annual modulation due to the seasonal variation of the earth-sun distance.