We report studies on mitigation of optical effects of bright Low-Earth-Orbit (LEO) satellites on Vera C. Rubin Observatory and its Legacy Survey of Space and Time (LSST). These include options for pointing the telescope to avoid satellites, laboratory investigations of bright trails on the Rubin Observatory LSST Camera sensors, algorithms for correcting image artifacts caused by bright trails, experiments on darkening SpaceX Starlink satellites, and ground-based follow-up observations. Starlink satellites with no brightness mitigation are presently g ~ 5.1 mag, and an initial experiment "DarkSat" is g ~ 6.1 mag. Future Starlink darkening plans may reach g ~ 7 mag, a brightness level that enables non-linear image artifact correction to the same level as background noise. However, the satellite trails will still exist at S/N ~ 100, generating systematics that may impact data analysis and limiting some science. LEO satellite trails are wider than a point-spread function because satellites are slightly out of focus due to their finite distance; for Rubin Observatory's 8.4-m mirror and a satellite at 550 km, this widening effect is about 3 arcsec, which helps avoid saturation by decreasing the trail's peak surface brightness.

Spectroscopic Eclipsing Binaries (SEBs) are fundamental benchmarks in stellar astrophysics and today are observed in breathtaking detail by missions like TESS, Kepler, and APOGEE. We develop a methodology for simultaneous analysis of high precision Kepler light curves and high resolution near-IR spectra from APOGEE and present orbital solutions and evolutionary histories for a subset of SEBs within this overlap. Radial velocities extracted from APOGEE spectra using the Broadening Function technique are combined with Kepler light curves and to yield binary orbital solutions. The Broadening Function approach yields more precise radial velocities than the standard Cross-Correlation Function, which in turn yields more precise orbital parameters and enables the identification of tertiary stars. The orbital periods of these seven SEBs range from 4 to 40 days. Four of the systems (KIC 5285607, KIC 6864859, KIC 6778289, and KIC 4285087) are well-detached binaries. The remaining three systems have apparent tertiary companions, but each exhibits two eclipses along with at least one spectroscopically varying component (KIC 6449358, KIC 6131659, and KIC 6781535). Gaia distances are available for four targets which we use to estimate temperatures of both members of these SEBs. We explore evolutionary histories in H-R diagram space and estimate ages for this subset of our sample. Finally, we consider the implications for the formation pathways of close binary systems via interactions with tertiary companions. Our methodology combined with the era of big data and observation overlap opens up the possibility of discovering and analyzing large numbers of diverse SEBs, including those with high flux ratios and those in triple systems.

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.

The Astropy Project (http://astropy.org) is, in its own words, "a community effort to develop a single core package for Astronomy in Python and foster interoperability between Python astronomy packages." For five years this project has been managed, written, and operated as a grassroots, self-organized, almost entirely volunteer effort while the software is used by the majority of the astronomical community. Despite this, the project has always been and remains to this day effectively unfunded. Further, contributors receive little or no formal recognition for creating and supporting what is now critical software. This paper explores the problem in detail, outlines possible solutions to correct this, and presents a few suggestions on how to address the sustainability of general purpose astronomical software.

Given the potential of ensemble asteroseismology for understanding fundamental properties of large numbers of stars, it is critical to determine the accuracy of the scaling relations on which these measurements are based. From several powerful validation techniques, all indications so far show that stellar radius estimates from the asteroseismic scaling relations are accurate to within a few percent. Eclipsing binary systems hosting at least one star with detectable solar-like oscillations constitute the ideal test objects for validating asteroseismic radius and mass inferences. By combining radial-velocity measurements and photometric time series of eclipses, it is possible to determine the masses and radii of each component of a double-lined spectroscopic binary. We report the results of a four-year radial-velocity survey performed with the \'echelle spectrometer of the Astrophysical Research Consortium's 3.5-m telescope and the APOGEE spectrometer at Apache Point Observatory. We compare the masses and radii of 10 red giants obtained by combining radial velocities and eclipse photometry with the estimates from the asteroseismic scaling relations. We find that the asteroseismic scaling relations overestimate red-giant radii by about 5% on average and masses by about 15% for stars at various stages of red-giant evolution. Systematic overestimation of mass leads to underestimation of stellar age, which can have important implications for ensemble asteroseismology used for Galactic studies. As part of a second objective, where asteroseismology is used for understanding binary systems, we confirm that oscillations of red giants in close binaries can be suppressed enough to be undetectable, an hypothesis that was proposed in a previous work.

We combine Kepler photometry with ground-based spectra to present a comprehensive dynamical model of the double red giant eclipsing binary KIC 9246715. While the two stars are very similar in mass (M1 = 2.171 [+0.006 / -0.008], M2 = 2.149 [+0.006 / -0.008] Msun) and radius (R1 = 8.37 [+0.03 / -0.07], R2 = 8.30 [+0.04 / -0.03] Rsun), an asteroseismic analysis finds one main set of solar-like oscillations with unusually low-amplitude, wide modes. A second set of oscillations from the other star may exist, but this marginal detection is extremely faint. Because the two stars are nearly twins, KIC 9246715 is a difficult target for a precise test of the asteroseismic scaling relations, which yield M = 2.17 +/- 0.14 Msun and R = 8.26 +/- 0.18 Rsun. Both stars are consistent with the inferred asteroseismic properties, but we suspect the main oscillator is Star 2 because it is less active than Star 1. We find evidence for stellar activity and modest tidal forces acting over the 171-day eccentric orbit, which are likely responsible for the essential lack of solar-like oscillations in one star and weak oscillations in the other. Mixed modes indicate the main oscillating star is on the secondary red clump (a core-He-burning star), and stellar evolution modeling supports this with a coeval history for a pair of red clump stars. This system is a useful case study and paves the way for a detailed analysis of more red giants in eclipsing binaries, an important benchmark for asteroseismology.

Red-giant stars are an incredible source of information for testing models of stellar evolution, as asteroseismology has opened up a window into their interiors. Such insights are a direct result of the unprecedented data from space missions Kepler and CoRoT as well as recent theoretical advances. Eclipsing binaries are also fundamental astrophysical objects, and when coupled with asteroseismology, binaries would provide two independent methods to obtain masses and radii and exciting opportunities to develop highly constrained stellar models. The possibility of discovering pulsating red giants in eclipsing binary systems is therefore an important goal that could potentially offer very robust characterization of these systems. Hitherto, only 1 case has been discovered with Kepler. We cross-correlated the detected red-giant and eclipsing-binary catalogs from Kepler data to find candidate systems. Light-curve modeling and mean asteroseismic properties are combined to yield measurements of periods, masses, radii, temperatures, eclipse timing variations, core rotation rates, and red-giant evolutionary state. After elimination of false positives, out of the 70 systems common to the red-giant and eclipsing-binary catalogs we find 13 new strong candidates that are pulsating red giants in eclipsing binaries, and 10 more that are pulsating red giants in hierarchical triple systems. The systems span a range of orbital eccentricities, periods, and spectral types for the red-giant companion. One case even suggests a system composed of two red-giants. The discovery of pulsating red giants in multiple systems provides exciting test beds for precise astrophysical modeling, and follow-up spectroscopic observations of the candidate systems are encouraged.

We present an improved method for determining the mass of neutron stars in eclipsing X-ray pulsar binaries and apply the method to six systems, namely Vela X-1, 4U 1538-52, SMC X-1, LMC X-4, Cen X-3, and Her X-1. In previous studies to determine neutron star mass, the X-ray eclipse duration has been approximated analytically by assuming the companion star is spherical with an effective Roche lobe radius. We use a numerical code based on Roche geometry with various optimizers to analyze the published data for these systems, which we supplement with new spectroscopic and photometric data for 4U 1538-52. This allows us to model the eclipse duration more accurately and thus calculate an improved value for the neutron star mass. The derived neutron star mass also depends on the assumed Roche lobe filling factor beta of the companion star, where beta = 1 indicates a completely filled Roche lobe. In previous work a range of beta between 0.9 and 1.0 was usually adopted. We use optical ellipsoidal lightcurve data to constrain beta. We find neutron star masses of 1.77 +/- 0.08 M_{sun} for Vela X-1, 0.87 +/- 0.07 M_{sun} for 4U 1538-52 (eccentric orbit), 1.00 +/- 0.10 M_{sun} for 4U 1538-52 (circular orbit), 1.04 +/- 0.09 M_{sun} for SMC X-1, 1.29 +/- 0.05 M_{sun} for LMC X-4, 1.49 +/- 0.08 M_{sun} for Cen X-3, and 1.07 +/- 0.36 M_{sun} for Her X-1. We discuss the limits of the approximations that were used to derive the earlier mass determinations, and we comment on the implications our new masses have for observationally refining the upper and lower bounds of the neutron star mass distribution.