Search

A key processing step in ground-based astronomy involves combining multiple noisy and blurry exposures to produce an image of the night sky with improved signal-to-noise ratio. Typically, this is achieved via image coaddition, and can be undertaken such that the resulting night sky image has enhanced spatial resolution. Yet, this task remains a formidable challenge despite decades of advancements. In this paper, we introduce ImageMM: a new framework based on the majorization-minimization algorithm for joint multi-frame astronomical image restoration and super-resolution. ImageMM uses multiple registered astronomical exposures to produce a non-parametric latent image of the night sky, prior to the atmosphere's impact on the observed exposures. Our framework also features a novel variational approach to compute refined point-spread functions of arbitrary resolution for the restoration and super-resolution procedure. Our algorithms, implemented in TensorFlow, leverage Graphics Processing Unit acceleration to produce latent images in near real-time, even when processing high-resolution exposures. We tested ImageMM on Hyper Suprime-Cam (HSC) exposures, which are a precursor for upcoming imaging data from the Rubin Observatory. The results are encouraging: ImageMM produces sharp latent images, in which spatial features of bright sources are revealed in unprecedented detail (e.g. the structure of spiral galaxies), and where faint sources that are usually indistinguishable from the noisy sky-background also become discernible, thus pushing the detection limits. Moreover, aperture photometry performed on the HSC pipeline coadd and ImageMM's latent images yielded consistent source detection and flux measurements, thereby demonstrating ImageMM's suitability for cutting-edge photometric studies with state-of-the-art astronomical imaging data.

The relationship between galaxy size and environment has remained enigmatic, with over a decade of conflicting results. We present one of the first comprehensive studies of the variation of galaxy radius with environment beyond the local Universe and demonstrate that large-scale environmental density is correlated with galaxy radius independent of stellar mass and galaxy morphology. We confirm with $>5\sigma$ confidence that galaxies in denser environments are up to $\sim25\%$ larger than their equally massive counterparts with similar morphology in less dense regions of the Universe. We achieve this result by correlating projected two-dimensional densities over $\sim360$ deg$^2$ with the structural parameters of $\sim3$ million Hyper Suprime-Cam galaxies at $0.3 \leq z < 0.7$ with $\log M/M_{\odot} \geq 8.9$. Compared to most previous studies, this sample is $\sim100-10,000$ times larger and goes $\sim1$ dex deeper in mass-completeness. We demonstrate that past conflicting results have been driven by small sample sizes and a lack of robust measurement uncertainties. We verify the presence of the above correlation separately for disk-dominated, bulge-dominated, star-forming, and quiescent subpopulations. We find the strength of the correlation to be dependent on redshift, stellar mass, and morphology. The correlation is strongest at lower redshifts and systematically weakens or disappears beyond $z \geq 0.5$. At $z\geq0.5$, more massive galaxies still display a statistically significant correlation. Although some existing theoretical frameworks can be selectively invoked to explain some of the observed correlations, our work demonstrates the need for more comprehensive theoretical investigations of the correlation between galaxy size and environment., Comment: 29 pages, 13 figures. Published in The Astrophysical Journal. We welcome comments and constructive criticism

Evaluating the accuracy and calibration of the redshift posteriors produced by photometric redshift (photo-z) estimators is vital for enabling precision cosmology and extragalactic astrophysics with modern wide-field photometric surveys. Evaluating photo-z posteriors on a per-galaxy basis is difficult, however, as real galaxies have a true redshift but not a true redshift posterior. We introduce PZFlow, a Python package for the probabilistic forward modeling of galaxy catalogs with normalizing flows. For catalogs simulated with PZFlow, there is a natural notion of "true" redshift posteriors that can be used for photo-z validation. We use PZFlow to simulate a photometric galaxy catalog where each galaxy has a redshift, noisy photometry, shape information, and a true redshift posterior. We also demonstrate the use of an ensemble of normalizing flows for photo-z estimation. We discuss how PZFlow will be used to validate the photo-z estimation pipeline of the Dark Energy Science Collaboration (DESC), and the wider applicability of PZFlow for statistical modeling of any tabular data., Comment: 19 pages, 13 figures, submitted to AJ

The Vera C. Rubin Observatory will, over a period of 10 years, repeatedly survey the southern sky. To ensure that images generated by Rubin meet the quality requirements for precision science, the observatory will use an Active Optics System (AOS) to correct for alignment and mirror surface perturbations introduced by gravity and temperature gradients in the optical system. To accomplish this Rubin will use out-of-focus images from sensors located at the edge of the focal plane to learn and correct for perturbations to the wavefront. We have designed and integrated a deep learning model for wavefront estimation into the AOS pipeline. In this paper, we compare the performance of this deep learning approach to Rubin's baseline algorithm when applied to images from two different simulations of the Rubin optical system. We show the deep learning approach is faster and more accurate, achieving the atmospheric error floor both for high-quality images, and low-quality images with heavy blending and vignetting. Compared to the baseline algorithm, the deep learning model is 40x faster, the median error 2x better under ideal conditions, 5x better in the presence of vignetting by the Rubin camera, and 14x better in the presence of blending in crowded fields. In addition, the deep learning model surpasses the required optical quality in simulations of the AOS closed loop. This system promises to increase the survey area useful for precision science by up to 8%. We discuss how this system might be deployed when commissioning and operating Rubin., Comment: 24 pages, 21 figures

We present the first set of trans-Neptunian objects (TNOs) observed on multiple nights in data taken from the DECam Ecliptic Exploration Project (DEEP). Of these 110 TNOs, 105 do not coincide with previously known TNOs and appear to be new discoveries. Each individual detection for our objects resulted from a digital tracking search at TNO rates of motion, using two to four hour exposure sets, and the detections were subsequently linked across multiple observing seasons. This procedure allows us to find objects with magnitudes $m_{VR} \approx 26$. The object discovery processing also included a comprehensive population of objects injected into the images, with a recovery and linking rate of at least $94\%$. The final orbits were obtained using a specialized orbit fitting procedure that accounts for the positional errors derived from the digital tracking procedure. Our results include robust orbits and magnitudes for classical TNOs with absolute magnitudes $H \sim 10$, as well as a dynamically detached object found at 76 au (semi-major axis $a\approx 77 \, \mathrm{au}$). We find a disagreement between our population of classical TNOs and the CFEPS-L7 three component model for the Kuiper belt., Comment: Accepted to AJ, companion paper do DEEP III. Objects will be released in the journal version (or contacting the authors)

We present a detailed study of the observational biases of the DECam Ecliptic Exploration Project's (DEEP) B1 data release and survey simulation software that enables direct statistical comparisons between models and our data. We inject a synthetic population of objects into the images, and then subsequently recover them in the same processing as our real detections. This enables us to characterize the survey's completeness as a function of apparent magnitudes and on-sky rates of motion. We study the statistically optimal functional form for the magnitude, and develop a methodology that can estimate the magnitude and rate efficiencies for all survey's pointing groups simultaneously. We have determined that our peak completeness is on average 80\% in each pointing group, and our magnitude drops to $25\%$ of this value at $m_{25} = 26.22$. We describe the freely available survey simulation software and its methodology. We conclude by using it to infer that our effective search area for objects at 40 au is $14.8\deg^2$, and that our lack of dynamically cold distant objects means that there at most $8\times 10^3$ objects with $60 < a < 80$ au and absolute magnitudes $H \leq 8$., Comment: Accepted to AJ, companion paper to DEEP VI

The DECam Ecliptic Exploration Project (DEEP) is a deep survey of the trans-Neptunian solar system being carried out on the 4-meter Blanco telescope at Cerro Tololo Inter-American Observatory in Chile using the Dark Energy Camera (DECam). By using a shift-and-stack technique to achieve a mean limiting magnitude of $r \sim 26.2$, DEEP achieves an unprecedented combination of survey area and depth, enabling quantitative leaps forward in our understanding of the Kuiper Belt populations. This work reports results from an analysis of twenty 3 sq.\ deg.\ DECam fields along the invariable plane. We characterize the efficiency and false-positive rates for our moving-object detection pipeline, and use this information to construct a Bayesian signal probability for each detected source. This procedure allows us to treat all of our Kuiper Belt Object (KBO) detections statistically, simultaneously accounting for efficiency and false positives. We detect approximately 2300 candidate sources with KBO-like motion at S/N $>6.5$. We use a subset of these objects to compute the luminosity function of the Kuiper Belt as a whole, as well as the Cold Classical (CC) population. We also investigate the absolute magnitude ($H$) distribution of the CCs, and find consistency with both an exponentially tapered power-law, which is predicted by streaming instability models of planetesimal formation, and a rolling power law. Finally, we provide an updated mass estimate for the Cold Classical Kuiper Belt of $M_{CC}(H_r < 12) = 0.0017^{+0.0010}_{-0.0004} M_{\oplus}$, assuming albedo $p = 0.15$ and density $\rho = 1$ g cm$^{-3}$., Comment: Accepted by PSJ

The Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) dataset will dramatically alter our understanding of the Universe, from the origins of the Solar System to the nature of dark matter and dark energy. Much of this research will depend on the existence of robust, tested, and scalable algorithms, software, and services. Identifying and developing such tools ahead of time has the potential to significantly accelerate the delivery of early science from LSST. Developing these collaboratively, and making them broadly available, can enable more inclusive and equitable collaboration on LSST science. To facilitate such opportunities, a community workshop entitled "From Data to Software to Science with the Rubin Observatory LSST" was organized by the LSST Interdisciplinary Network for Collaboration and Computing (LINCC) and partners, and held at the Flatiron Institute in New York, March 28-30th 2022. The workshop included over 50 in-person attendees invited from over 300 applications. It identified seven key software areas of need: (i) scalable cross-matching and distributed joining of catalogs, (ii) robust photometric redshift determination, (iii) software for determination of selection functions, (iv) frameworks for scalable time-series analyses, (v) services for image access and reprocessing at scale, (vi) object image access (cutouts) and analysis at scale, and (vii) scalable job execution systems. This white paper summarizes the discussions of this workshop. It considers the motivating science use cases, identified cross-cutting algorithms, software, and services, their high-level technical specifications, and the principles of inclusive collaborations needed to develop them. We provide it as a useful roadmap of needs, as well as to spur action and collaboration between groups and individuals looking to develop reusable software for early LSST science., Comment: White paper from "From Data to Software to Science with the Rubin Observatory LSST" workshop

In this Letter, we report the discovery of an ultraluminous fast-evolving transient in rest-frame UV wavelengths, MUSSES2020J, soon after its occurrence by using the Hyper Suprime-Cam (HSC) mounted on the 8.2 m Subaru telescope. The rise time of about 5 days with an extremely high UV peak luminosity shares similarities to a handful of fast blue optical transients whose peak luminosities are comparable with the most luminous supernovae while their timescales are significantly shorter (hereafter "fast blue ultraluminous transient," FBUT). In addition, MUSSES2020J is located near the center of a normal low-mass galaxy at a redshift of 1.063, suggesting a possible connection between the energy source of MUSSES2020J and the central part of the host galaxy. Possible physical mechanisms powering this extreme transient such as a wind-driven tidal disruption event and an interaction between supernova and circumstellar material are qualitatively discussed based on the first multiband early-phase light curve of FBUTs, although whether the scenarios can quantitatively explain the early photometric behavior of MUSSES2020J requires systematical theoretical investigations. Thanks to the ultrahigh luminosity in UV and blue optical wavelengths of these extreme transients, a promising number of FBUTs from the local to the high-z universe can be discovered through deep wide-field optical surveys in the near future., Comment: 14 pages, 5 figures, accepted for publication in The Astrophysical Journal Letters

Reinforcement learning (RL) tasks are typically framed as Markov Decision Processes (MDPs), assuming that decisions are made at fixed time intervals. However, many applications of great importance, including healthcare, do not satisfy this assumption, yet they are commonly modelled as MDPs after an artificial reshaping of the data. In addition, most healthcare (and similar) problems are offline by nature, allowing for only retrospective studies. To address both challenges, we begin by discussing the Semi-MDP (SMDP) framework, which formally handles actions of variable timings. We next present a formal way to apply SMDP modifications to nearly any given value-based offline RL method. We use this theory to introduce three SMDP-based offline RL algorithms, namely, SDQN, SDDQN, and SBCQ. We then experimentally demonstrate that only these SMDP-based algorithms learn the optimal policy in variable-time environments, whereas their MDP counterparts do not. Finally, we apply our new algorithms to a real-world offline dataset pertaining to warfarin dosing for stroke prevention and demonstrate similar results., Comment: Published at Conference on Health, Inference, and Learning (CHIL) 2022

Methods based on machine learning have recently made substantial inroads in many corners of cosmology. Through this process, new computational tools, new perspectives on data collection, model development, analysis, and discovery, as well as new communities and educational pathways have emerged. Despite rapid progress, substantial potential at the intersection of cosmology and machine learning remains untapped. In this white paper, we summarize current and ongoing developments relating to the application of machine learning within cosmology and provide a set of recommendations aimed at maximizing the scientific impact of these burgeoning tools over the coming decade through both technical development as well as the fostering of emerging communities., Comment: Contribution to Snowmass 2021. 32 pages

Trans-Neptunian Objects (TNOs) provide a window into the history of the Solar System, but they can be challenging to observe due to their distance from the Sun and relatively low brightness. Here we report the detection of 75 moving objects that we could not link to any other known objects, the faintest of which has a VR magnitude of $25.02 \pm 0.93$ using the KBMOD platform. We recover an additional 24 sources with previously-known orbits. We place constraints on the barycentric distance, inclination, and longitude of ascending node of these objects. The unidentified objects have a median barycentric distance of 41.28 au, placing them in the outer Solar System. The observed inclination and magnitude distribution of all detected objects is consistent with previously published KBO distributions. We describe extensions to KBMOD, including a robust percentile-based lightcurve filter, an in-line graphics processing unit (GPU) filter, new coadded stamp generation, and a convolutional neural network (CNN) stamp filter, which allow KBMOD to take advantage of difference images. These enchancements mark a significant improvement in the readiness of KBMOD for deployment on future big data surveys such as LSST., Comment: Accepted: Astronomical Journal

Ion irradiation provides a promising substitute to neutron tests for investigating the effects of radiation on materials for fission and fusion reactor plants. Here we show proton irradiation can quantitatively reproduce precipitation that leads to embrittlement in reactor pressure vessel steels, at dose rates 10,000 times greater than experienced in fission reactor operation. Small-angle neutron scattering (SANS) is used to characterize precipitate size distributions in copper-containing steels irradiated to average doses of approx. 7 mdpa with 5 MeV protons. Comparing our results with the literature on reactor pressure vessel steels containing at least 1 at.% nickel, we find a power-law scaling of dose with exponent 0.25--0.30 accounts for the effects of dose rate on precipitate volume fraction over 6 orders of magnitude in dose rate. In conjunction with dose rate, carbon is identified as performing a leading role in determining precipitate sizes, adding to the known effects of nickel, manganese and irradiation temperature. We discuss the composition of precipitates inferred from SANS, taking previous atom probe tomography studies into consideration., Comment: 18 pages, 9 figures

Vera C. Rubin Observatory is a ground-based astronomical facility under construction, a joint project of the National Science Foundation and the U.S. Department of Energy, designed to conduct a multi-purpose 10-year optical survey of the southern hemisphere sky: the Legacy Survey of Space and Time. Significant flexibility in survey strategy remains within the constraints imposed by the core science goals of probing dark energy and dark matter, cataloging the Solar System, exploring the transient optical sky, and mapping the Milky Way. The survey's massive data throughput will be transformational for many other astrophysics domains and Rubin's data access policy sets the stage for a huge potential users' community. To ensure that the survey science potential is maximized while serving as broad a community as possible, Rubin Observatory has involved the scientific community at large in the process of setting and refining the details of the observing strategy. The motivation, history, and decision-making process of this strategy optimization are detailed in this paper, giving context to the science-driven proposals and recommendations for the survey strategy included in this Focus Issue., Comment: Submitted as the opening paper of the Astrophysical Journal Focus Issue on Rubin LSST cadence and survey strategy

We present "Tracklet-less Heliocentric Orbit Recovery" (THOR), an algorithm for linking of observations of Solar System objects across multiple epochs that does not require intra-night tracklets or a predefined cadence of observations within a search window. By sparsely covering regions of interest in the phase space with "test orbits", transforming nearby observations over a few nights into the co-rotating frame of the test orbit at each epoch, and then performing a generalized Hough transform on the transformed detections followed by orbit determination (OD) filtering, candidate clusters of observations belonging to the same objects can be recovered at moderate computational cost and little to no constraints on cadence. We validate the effectiveness of this approach by running on simulations as well as on real data from the Zwicky Transient Facility (ZTF). Applied to a short, 2-week, slice of ZTF observations, we demonstrate THOR can recover 97.4% of all previously known and discoverable objects in the targeted ($a > 1.7$ au) population with 5 or more observations and with purity between 97.7% and 100%. This includes 10 likely new discoveries, and a recovery of an $e \sim 1$ comet C/2018 U1 (the comet would have been a ZTF discovery had THOR been running in 2018 when the data were taken). The THOR package and demo Jupyter notebooks are open source and available at https://github.com/moeyensj/thor., Comment: 22 pages, 11 figures