Publications of the Astronomical Society of the Pacific

Old, digitized astronomical images taken before the human spacefaring age offer a rare glimpse of the sky before the era of artificial satellites. In this paper, we present the first optical searches for artificial objects with high specular reflections near the Earth. We follow the method proposed in Villarroel et al. and use a transient sample drawn from Solano et al. We use images from the First Palomar Sky Survey to search for multiple (within a plate exposure) transients that, in addition to being point-like, are aligned along a narrow band. We provide a shortlist of the most promising candidate alignments, including one with ∼3.9σ statistical significance. These aligned transients remain difficult to explain with known phenomena, even if rare optical ghosting producing point-like sources cannot be fully excluded at present. We explore remaining possibilities, including fast reflections from highly reflective objects in geosynchronous orbit, or emissions from artificial sources high above Earth’s atmosphere. We also find a highly significant (∼22σ) deficit of POSS-I transients within Earth's shadow when compared with the theoretical hemispheric shadow coverage at 42,164 km altitude. The deficit is still present though at reduced significance (∼7.6σ) when a more realistic plate-based coverage is considered. This study should be viewed as an initial exploration into the potential of archival photographic surveys to reveal transient phenomena, and we hope it motivates more systematic searches across historical data sets.

We introduce a stable, well tested Python implementation of the affine-invariant ensemble sampler for Markov chain Monte Carlo (MCMC) proposed by Goodman & Weare (2010). The code is open source and has already been used in several published projects in the astrophysics literature. The algorithm behind emcee has several advantages over traditional MCMC sampling methods and it has excellent performance as measured by the autocorrelation time (or function calls per independent sample). One major advantage of the algorithm is that it requires hand-tuning of only 1 or 2 parameters compared to ∼N² for a traditional algorithm in an N-dimensional parameter space. In this document, we describe the algorithm and the details of our implementation. Exploiting the parallelism of the ensemble method, emcee permits any user to take advantage of multiple CPU cores without extra effort. The code is available online at http://dan.iel.fm/emcee under the GNU General Public License v2.

Twenty-six years ago a small committee report, building on earlier studies, expounded a compelling and poetic vision for the future of astronomy, calling for an infrared-optimized space telescope with an aperture of at least 4 m. With the support of their governments in the US, Europe, and Canada, 20,000 people realized that vision as the 6.5 m James Webb Space Telescope. A generation of astronomers will celebrate their accomplishments for the life of the mission, potentially as long as 20 yr, and beyond. This report and the scientific discoveries that follow are extended thank-you notes to the 20,000 team members. The telescope is working perfectly, with much better image quality than expected. In this and accompanying papers, we give a brief history, describe the observatory, outline its objectives and current observing program, and discuss the inventions and people who made it possible. We cite detailed reports on the design and the measured performance on orbit.

The Milankovitch cycles of Earth result from gravitational interactions with other bodies in the solar system. These interactions lead to slow changes in the orbit and angular momentum vector of Earth, and correspondingly influence Earth’s climate evolution. Several studies have shown that Mars may play a significant role in these Milankovitch cycles, such as the 2.4 Myr eccentricity cycle related to perihelion precession dynamics. Here we provide the results of a detailed dynamical analysis that explores the Earth Milankovitch cycles as a function of the Martian mass to quantify the extent that Mars influences variations in Earth’s orbital eccentricity, the longitude of perihelion, the longitude of the ascending node, and obliquity (axial tilt). Our results show that, although the 405 kyr long-eccentricity metronome driven by g₂ (Venus) and g₅ (Jupiter) persists at all Mars masses, the ∼100 kyr short-eccentricity bands driven by g₄ (Mars) lengthen and gain power as Mars becomes more massive, consistent with enhanced coupling among inner-planet g-modes. The 2.4 Myr grand cycle is absent when Mars approaches zero mass, reflecting the movement of g₄ with the Martian mass. Meanwhile, Earth’s obliquity cycles driven by s₃ (Earth) and s₄ (Mars) lengthen from the canonical ∼41 kyr with increasing Mars mass, relocating to a dominant 45–55 kyr band when the mass of Mars is an order of magnitude larger than its present value. These results establish how Mars’ mass controls the architecture of Earth’s climate-forcing spectrum and that the Milankovitch spectrum of an Earth-like planet is a sensitive, interpretable probe of its planetary neighborhood.

The Zwicky Transient Facility (ZTF) is a new optical time-domain survey that uses the Palomar 48 inch Schmidt telescope. A custom-built wide-field camera provides a 47 deg² field of view and 8 s readout time, yielding more than an order of magnitude improvement in survey speed relative to its predecessor survey, the Palomar Transient Factory. We describe the design and implementation of the camera and observing system. The ZTF data system at the Infrared Processing and Analysis Center provides near-real-time reduction to identify moving and varying objects. We outline the analysis pipelines, data products, and associated archive. Finally, we present on-sky performance analysis and first scientific results from commissioning and the early survey. ZTF’s public alert stream will serve as a useful precursor for that of the Large Synoptic Survey Telescope.

This paper characterizes the actual science performance of the James Webb Space Telescope (JWST), as determined from the six month commissioning period. We summarize the performance of the spacecraft, telescope, science instruments, and ground system, with an emphasis on differences from pre-launch expectations. Commissioning has made clear that JWST is fully capable of achieving the discoveries for which it was built. Moreover, almost across the board, the science performance of JWST is better than expected; in most cases, JWST will go deeper faster than expected. The telescope and instrument suite have demonstrated the sensitivity, stability, image quality, and spectral range that are necessary to transform our understanding of the cosmos through observations spanning from near-earth asteroids to the most distant galaxies.

The James Webb Space Telescope (JWST) is a large, infrared space telescope that has recently started its science program which will enable breakthroughs in astrophysics and planetary science. Notably, JWST will provide the very first observations of the earliest luminous objects in the universe and start a new era of exoplanet atmospheric characterization. This transformative science is enabled by a 6.6 m telescope that is passively cooled with a 5 layer sunshield. The primary mirror is comprised of 18 controllable, low areal density hexagonal segments, that were aligned and phased relative to each other in orbit using innovative image-based wave front sensing and control algorithms. This revolutionary telescope took more than two decades to develop with a widely distributed team across engineering disciplines. We present an overview of the telescope requirements, architecture, development, superb on-orbit performance, and lessons learned. JWST successfully demonstrates a segmented aperture space telescope and establishes a path to building even larger space telescopes.

We review recent determinations of the present‐day mass function (PDMF) and initial mass function (IMF) in various components of the Galaxy—disk, spheroid, young, and globular clusters—and in conditions characteristic of early star formation. As a general feature, the IMF is found to depend weakly on the environment and to be well described by a power‐law form for m≳1 M_⊙ and a lognormal form below, except possibly for early star formation conditions. The disk IMF for single objects has a characteristic mass around m_c ∼ 0.08 M_⊙ and a variance in logarithmic mass σ ∼ 0.7, whereas the IMF for multiple systems has m_c ∼ 0.2 M_⊙ and σ ∼ 0.6. The extension of the single MF into the brown dwarf regime is in good agreement with present estimates of L‐ and T‐dwarf densities and yields a disk brown dwarf number density comparable to the stellar one, n_BD ∼ n_* ∼ 0.1 pc⁻³. The IMF of young clusters is found to be consistent with the disk field IMF, providing the same correction for unresolved binaries, confirming the fact that young star clusters and disk field stars represent the same stellar population. Dynamical effects, yielding depletion of the lowest mass objects, are found to become consequential for ages ≳130 Myr. The spheroid IMF relies on much less robust grounds. The large metallicity spread in the local subdwarf photometric sample, in particular, remains puzzling. Recent observations suggest that there is a continuous kinematic shear between the thick‐disk population, present in local samples, and the genuine spheroid one. This enables us to derive only an upper limit for the spheroid mass density and IMF. Within all the uncertainties, the latter is found to be similar to the one derived for globular clusters and is well represented also by a lognormal form with a characteristic mass slightly larger than for the disk, m_c ∼ 0.2–0.3 M_⊙, excluding a significant population of brown dwarfs in globular clusters and in the spheroid. The IMF characteristic of early star formation at large redshift remains undetermined, but different observational constraints suggest that it does not extend below ∼1 M_⊙. These results suggest a characteristic mass for star formation that decreases with time, from conditions prevailing at large redshift to conditions characteristic of the spheroid (or thick disk) to present‐day conditions. These conclusions, however, remain speculative, given the large uncertainties in the spheroid and early star IMF determinations.

These IMFs allow a reasonably robust determination of the Galactic present‐day and initial stellar and brown dwarf contents. They also have important galactic implications beyond the Milky Way in yielding more accurate mass‐to‐light ratio determinations. The mass‐to‐light ratios obtained with the disk and the spheroid IMF yield values 1.8–1.4 times smaller than for a Salpeter IMF, respectively, in agreement with various recent dynamical determinations. This general IMF determination is examined in the context of star formation theory. None of the theories based on a Jeans‐type mechanism, where fragmentation is due only to gravity, can fulfill all the observational constraints on star formation and predict a large number of substellar objects. On the other hand, recent numerical simulations of compressible turbulence, in particular in super‐Alfvénic conditions, seem to reproduce both qualitatively and quantitatively the stellar and substellar IMF and thus provide an appealing theoretical foundation. In this picture, star formation is induced by the dissipation of large‐scale turbulence to smaller scales through radiative MHD shocks, producing filamentary structures. These shocks produce local nonequilibrium structures with large density contrasts, which collapse eventually in gravitationally bound objects under the combined influence of turbulence and gravity. The concept of a single Jeans mass is replaced by a distribution of local Jeans masses, representative of the lognormal probability density function of the turbulent gas. Objects below the mean thermal Jeans mass still have a possibility to collapse, although with a decreasing probability.

We present data from the Apache Point Observatory Lunar Laser-ranging Operation (APOLLO) covering the 15 yr span from 2006 April through the end of 2020. APOLLO measures the Earth–Moon separation by recording the round-trip travel time of photons from the Apache Point Observatory to five retro-reflector arrays on the Moon. The APOLLO data set, combined with the 50 yr archive of measurements from other lunar laser ranging (LLR) stations, can be used to probe fundamental physics such as gravity and Lorentz symmetry, as well as properties of the Moon itself. We show that range measurements performed by APOLLO since 2006 have a median nightly accuracy of 1.7 mm, which is significantly better than other LLR stations.

The James Webb Space Telescope (JWST) is NASA’s flagship mission successor to the highly successful Hubble Space Telescope. It is an infrared observatory featuring a cryogenic 6.6 m aperture, deployable Optical Telescope Element (OTE) with a payload of four science instruments (SIs) assembled into an Integrated Science Instrument Module (ISIM) that provide imagery and spectroscopy in the near-infrared band between 0.6 and 5 μm and in the mid-infrared band between 5 and 28.1 μm. JWST was successfully launched on 2021 December 25 aboard an Ariane 5 launch vehicle. All 50 major deployments were successfully completed on 2022 January 8. The observatory performed all midcourse correction maneuvers and achieved its operational mission orbit around the Sun–Earth second Lagrange point (L2). All commissioning and calibration activities have been completed, and JWST has begun its science mission. This paper will provide a description of the driving requirements and their technical challenges, the engineering processes involved in the design formulation, the resulting observatory design, the verification programs that proved it to be flightworthy, and the measured on-orbit performance of the observatory. Since companion papers will describe the details of the OTE and SIs, this paper will concentrate on describing the key features of the observatory architecture that accommodates these elements, particularly those features and capabilities associated with accommodating the radiometric and image-quality performance.

Scientific CMOS (sCMOS) represents a novel type of X-ray focal plane detector developed over the past decade. Micrometeoroids and space/orbital debris (MMOD) pose significant hazards to onboard X-ray instruments and may degrade the performance of their focal plane detectors. In this work, we evaluate performance degradation in the sCMOS subjected to MMOD impacts by means of ground-based hypervelocity impact experiments. Projectile particles with diameters ranging from 1 to 50 μm, and average velocities between 1 and 6 km s⁻¹ were used. The impact target is a scaled dummy module of an X-ray telescope consisting of two micro-pore optics (MPO) mirrors in front as the focusing lens, and an aluminum-coated sCMOS serving as the focal plane detector. Because the sCMOS in this dummy module cannot be biased during the impact tests, its performance was characterized before and after the impact tests with a stand-alone sCMOS detector system. Scanning electron microscopy images show that projectile particles or secondary spallation products from the MPO reached the sCMOS and produced craters. Post-impact analyses reveal that significant increases in optical transmission is observed in some pixels, attributed to the localized peeling of the aluminum coating. For impacts involving projectile particles smaller than 10 μm, the dark current changes are minor and no dead pixels are observed. In contrast, impacts produced by 50 μm projectile particles caused 9 dead pixels. Excluding the damaged pixels, the sCMOS remains functional after the impact tests, exhibiting only minor performance degradation and thus demonstrating a degree of reliability against hypervelocity particle impacts.

Integral field spectroscopy has been added as a new observation mode to the Robert Stobie Spectrograph (RSS), the workhorse multi-mode instrument on the Southern African Large Telescope. RSS operates as an imaging spectrograph covering 320–900 nm with a spectral resolution–slit-width product up to 6600″. Using fiber optics and prismatic fold mirrors, we have been able to construct compact integral field units (IFUs) that fit within the same volume as the long-slit cassettes (134 mm ×130 mm × 8 mm). These “slit mask” IFUs (SMIs) direct the telescope beam into a two-dimensional sky-facing fiber array routed in the focal plane dimension into an 8$^{\prime} \rm{}\,\rm{}\,$one-dimensional pseudo-slit, with fiber output redirected back into the spectrograph collimator. The first completed unit, SMI-200, features 303 object fibers and 24 sky fibers, providing a spatial resolution of 0$\mathop{.}\limits{^{\prime\prime} }$8 (200 μm core diameter) over a field of view of 22$\mathop{.}\limits{^{\prime\prime} }$5 ×17$\mathop{.}\limits{^{\prime\prime} }$6. This paper describes the specific design considerations and design and fabrication strategies to maximize performance and minimize risk during construction given the demanding and highly constrained cassette geometry. We also detail mapping and laboratory characterization of the IFU. Laboratory measurements demonstrate a total throughput of 77%, but an effective throughput of only 55%–60% within the RSS collimator acceptance beam of f/4.2 due to losses dominated by focal ratio degradation induced by the detailed routing within the tight cassette volume.

Dust grains in the interstellar medium interact with photons across the electromagnetic spectrum. They are generally photon energy converters, absorbing short wavelength radiation and emitting long wavelength radiation. Sixty years ago in 1965, thermal emission from dust grains in the interstellar medium was discovered. This tutorial is a summary of the physics of thermal dust continuum emission and how to use observations of the intensity and flux density of dusty objects to calculate physical properties such as mass, column density, luminosity, dust temperature, and dust opacity spectral index. Equations are derived, when feasible, from first principles with all limits and assumptions explicitly stated. Properties of dust opacities appropriate for different astrophysical environments (e.g., diffuse ISM, dense cores, protoplanetary disks) are discussed and tabulated for the wavelengths of past, current, and future bolometer cameras. Corrections for observations at high redshift as well as the effects of telescope measurement limitations are derived. We also update the calculation of the mean molecular weight in different ISM environments and find that it is 1.404 per H atom, 2.809 per H₂ molecule, and 2.351 per gas particle assuming protosolar metallicity and the latest values of the ISM gas phase abundances of metals.

Path-length diversity methods may be used for adaptive optics systems to retrieve phase and amplitude information by measuring intensity across multiple planes. Observations that rely on free-space propagation, such as the nonlinear curvature wavefront sensor (WFS), have been shown to offer excellent sensitivity and robustness to scintillation. However, the default design results in a large opto-mechanical footprint due to unavoidable geometric-optics and wave-optics effects. Measurements recorded in a convergent beam would improve instrument compactness, while concentrating light into smaller detector regions of interest, improving signal-to-noise ratio and possibly wavefront reconstruction speed. In this paper, we study path-length diversity wavefront sensing using four planes of contemporaneous intensity measurements made in a convergent beam. We develop a physical optics propagation model and validate the model by performing wavefront reconstructions in both simulations and lab experiments. The manuscript’s core contribution is a practical, intensity-domain, Fourier-transform-based recipe to use a conventional multi-plane Gerchberg–Saxton (or comparable) reconstruction pipeline with convergent-beam measurements, enabling a compact optical layout. We find that this approach offers practical benefits over an equivalent free-space WFS, in particular reducing size, weight, complexity and cost.

BE Lyncis (BE Lyn) is a well-studied high-amplitude δ Scuti variable star. Recently, Niu et al. analyzed a 39 yr baseline of times of maximum light of BE Lyn, reporting that it is the most eccentric binary known (e ≈ 0.9989) and hosts the nearest black hole (BH) known to date. We analyze Hipparcos and Gaia astrometry of BE Lyn, predicting what the observed proper motion anomaly (PMA) over the 25 yr baseline between the two missions would be were the companion really a ≳17.5 M_⊙ BH. We find that the predicted PMA is at least an order of magnitude larger than the observed value of ≈1.7 ± 0.8 mas yr⁻¹, regardless of the assumed orientation of the orbit. We predict the expected Gaia DR3 RUWE for different orientations of the putative BH binary, finding that it ranges from ≈2.5 to 4.0, much larger than the reported value of 1.073. The observed value is instead consistent with a low-mass secondary or a single star. We find that BE Lyncis would have received a 7-parameter acceleration solution if it were a BH binary, in contradiction with its absence from the Gaia DR3 non-single star catalogs. Finally, we show that the reported orbit is impossible because the luminous star would overflow its Roche lobe at periastron, irrespective of inclination. We recommend caution in interpreting light-travel time effect models that require very high eccentricities, face-on inclinations, or large companion masses. The observed pulsation timing variations are most likely simply a result of red noise or pulsation phase evolution.

Dust grains in the interstellar medium interact with photons across the electromagnetic spectrum. They are generally photon energy converters, absorbing short wavelength radiation and emitting long wavelength radiation. Sixty years ago in 1965, thermal emission from dust grains in the interstellar medium was discovered. This tutorial is a summary of the physics of thermal dust continuum emission and how to use observations of the intensity and flux density of dusty objects to calculate physical properties such as mass, column density, luminosity, dust temperature, and dust opacity spectral index. Equations are derived, when feasible, from first principles with all limits and assumptions explicitly stated. Properties of dust opacities appropriate for different astrophysical environments (e.g., diffuse ISM, dense cores, protoplanetary disks) are discussed and tabulated for the wavelengths of past, current, and future bolometer cameras. Corrections for observations at high redshift as well as the effects of telescope measurement limitations are derived. We also update the calculation of the mean molecular weight in different ISM environments and find that it is 1.404 per H atom, 2.809 per H₂ molecule, and 2.351 per gas particle assuming protosolar metallicity and the latest values of the ISM gas phase abundances of metals.

Measurements of physical parameters for stars and (exo)planets are often quoted in units normalized to the Sun and/or Earth. The nominal total solar irradiance, ${S}_{\odot }^{{\rm{N}}}$, while based on a current best estimate with uncertainties, was adopted to be an exact reference value of 1361 W m⁻² by IAU 2015 Resolution B3, corresponding to “the mean total electromagnetic energy from the Sun, integrated over all wavelengths, incident per unit area per unit time at distance 1 au.” In the planetary and exoplanetary science literature, the units employed for “flux,” “insolation,” “instellation,” etc., are often cumbersome or inconsistent. To simplify the quoting of irradiance units for astronomical applications, we introduce the portmanteau solirad, short for solar irradiance, as an abbreviated version of the longer IAU term “nominal total solar irradiance.” The solirad (So) is a unit of irradiance, where 1 solirad = 1 So = 1361 W m⁻², equivalent to the IAU nominal total solar irradiance, and to an apparent bolometric magnitude of m_bol ≃ −26.832 mag (per IAU 2015 Resolution B2).

Searching for life elsewhere in the universe is one of the most highly prioritized pursuits in astronomy today. However, the ability to observe evidence of Earth-like life through biosignatures is limited by the number of planets in the solar neighborhood with conditions similar to Earth. The occurrence rate of Earth-like planets in the habitable zones of Sun-like stars, η_⊕, is therefore crucial for addressing the apparent lack of consensus on its value in the literature. Here we present a review of the current understanding of η_⊕. We first provide definitions for parameters that contribute to η_⊕. Then, we discuss the previous and current estimated parameter values and the context of the limitations on the analyses that produced these estimates. We compile an extensive list of the factors that go into any calculation of η_⊕, and how detection techniques and surveys differ in their sensitivity and ability to accurately constrain η_⊕. Understanding and refining the value of η_⊕ is crucial for upcoming missions and telescopes, such as the planned Habitable Worlds Observatory and the Large Interferometer for Exoplanets, which aim to search for biosignatures on exoplanets in the solar neighborhood.

Exozodiacal dust disks (exozodis) are populations of warm (∼300 K) or hot (∼1000 K) dust, located in or interior to a star’s habitable zone, detected around ∼25% of main-sequence stars as excess emission over the stellar photosphere at mid- or near-infrared wavelengths. Often too plentiful to be explained by an in situ planetesimal belt, exozodi dust is usually thought to be transported inwards from further out in the system. There is no consensus on which (if any) of various proposed dynamical models is correct, yet it is vital to understand exozodis given the risk they pose to direct imaging and characterisation of Earth-like planets. This article reviews current theoretical understanding of the origin and evolution of exozodi dust. It also identifies key questions pertinent to the potential for exozodis to impact exoplanet imaging and summarises current understanding of the answer to them informed by exozodi theory. These address how exozodi dust is delivered, its size and spatial distribution, and the effect of its composition on exozodi observability, as well as the connection between hot and warm exozodis. Also addressed are how common different exozodi levels are and how that level can be predicted from system properties, as well as the features that planets impart in dust distributions and how exozodis affect a planet’s physical properties and habitability. We conclude that exozodis present both a problem and an opportunity, e.g., by introducing noise that makes planets harder to detect, but also identifying systems in which ingredients conducive to life, like water and volatiles, are delivered to the habitable zone.

This tutorial is an introduction to observational studies of dust transport and evolution in protoplanetary disks. Spatially resolved observations of disks at multiple wavelengths can allow us to infer the distribution of various dust grains and gas species. Combining these observations offers a more complete understanding of dust structure and properties across different disk locations. For example, by better characterizing the disk's vertical structure, observations help to constrain the level of vertical settling and identify regions of high dust density, which are favorable for grain growth and planet formation. This tutorial describes various methodologies for inferring dust properties and vertical height of different tracers, as an introduction for beginners.

Integral field spectroscopy has been added as a new observation mode to the Robert Stobie Spectrograph (RSS), the workhorse multi-mode instrument on the Southern African Large Telescope. RSS operates as an imaging spectrograph covering 320–900 nm with a spectral resolution–slit-width product up to 6600″. Using fiber optics and prismatic fold mirrors, we have been able to construct compact integral field units (IFUs) that fit within the same volume as the long-slit cassettes (134 mm ×130 mm × 8 mm). These “slit mask” IFUs (SMIs) direct the telescope beam into a two-dimensional sky-facing fiber array routed in the focal plane dimension into an 8$^{\prime} \rm{}\,\rm{}\,$one-dimensional pseudo-slit, with fiber output redirected back into the spectrograph collimator. The first completed unit, SMI-200, features 303 object fibers and 24 sky fibers, providing a spatial resolution of 0$\mathop{.}\limits{^{\prime\prime} }$8 (200 μm core diameter) over a field of view of 22$\mathop{.}\limits{^{\prime\prime} }$5 ×17$\mathop{.}\limits{^{\prime\prime} }$6. This paper describes the specific design considerations and design and fabrication strategies to maximize performance and minimize risk during construction given the demanding and highly constrained cassette geometry. We also detail mapping and laboratory characterization of the IFU. Laboratory measurements demonstrate a total throughput of 77%, but an effective throughput of only 55%–60% within the RSS collimator acceptance beam of f/4.2 due to losses dominated by focal ratio degradation induced by the detailed routing within the tight cassette volume.

Dust grains in the interstellar medium interact with photons across the electromagnetic spectrum. They are generally photon energy converters, absorbing short wavelength radiation and emitting long wavelength radiation. Sixty years ago in 1965, thermal emission from dust grains in the interstellar medium was discovered. This tutorial is a summary of the physics of thermal dust continuum emission and how to use observations of the intensity and flux density of dusty objects to calculate physical properties such as mass, column density, luminosity, dust temperature, and dust opacity spectral index. Equations are derived, when feasible, from first principles with all limits and assumptions explicitly stated. Properties of dust opacities appropriate for different astrophysical environments (e.g., diffuse ISM, dense cores, protoplanetary disks) are discussed and tabulated for the wavelengths of past, current, and future bolometer cameras. Corrections for observations at high redshift as well as the effects of telescope measurement limitations are derived. We also update the calculation of the mean molecular weight in different ISM environments and find that it is 1.404 per H atom, 2.809 per H₂ molecule, and 2.351 per gas particle assuming protosolar metallicity and the latest values of the ISM gas phase abundances of metals.

Path-length diversity methods may be used for adaptive optics systems to retrieve phase and amplitude information by measuring intensity across multiple planes. Observations that rely on free-space propagation, such as the nonlinear curvature wavefront sensor (WFS), have been shown to offer excellent sensitivity and robustness to scintillation. However, the default design results in a large opto-mechanical footprint due to unavoidable geometric-optics and wave-optics effects. Measurements recorded in a convergent beam would improve instrument compactness, while concentrating light into smaller detector regions of interest, improving signal-to-noise ratio and possibly wavefront reconstruction speed. In this paper, we study path-length diversity wavefront sensing using four planes of contemporaneous intensity measurements made in a convergent beam. We develop a physical optics propagation model and validate the model by performing wavefront reconstructions in both simulations and lab experiments. The manuscript’s core contribution is a practical, intensity-domain, Fourier-transform-based recipe to use a conventional multi-plane Gerchberg–Saxton (or comparable) reconstruction pipeline with convergent-beam measurements, enabling a compact optical layout. We find that this approach offers practical benefits over an equivalent free-space WFS, in particular reducing size, weight, complexity and cost.

BE Lyncis (BE Lyn) is a well-studied high-amplitude δ Scuti variable star. Recently, Niu et al. analyzed a 39 yr baseline of times of maximum light of BE Lyn, reporting that it is the most eccentric binary known (e ≈ 0.9989) and hosts the nearest black hole (BH) known to date. We analyze Hipparcos and Gaia astrometry of BE Lyn, predicting what the observed proper motion anomaly (PMA) over the 25 yr baseline between the two missions would be were the companion really a ≳17.5 M_⊙ BH. We find that the predicted PMA is at least an order of magnitude larger than the observed value of ≈1.7 ± 0.8 mas yr⁻¹, regardless of the assumed orientation of the orbit. We predict the expected Gaia DR3 RUWE for different orientations of the putative BH binary, finding that it ranges from ≈2.5 to 4.0, much larger than the reported value of 1.073. The observed value is instead consistent with a low-mass secondary or a single star. We find that BE Lyncis would have received a 7-parameter acceleration solution if it were a BH binary, in contradiction with its absence from the Gaia DR3 non-single star catalogs. Finally, we show that the reported orbit is impossible because the luminous star would overflow its Roche lobe at periastron, irrespective of inclination. We recommend caution in interpreting light-travel time effect models that require very high eccentricities, face-on inclinations, or large companion masses. The observed pulsation timing variations are most likely simply a result of red noise or pulsation phase evolution.

ATUS, the Astronomical Telescope of the University of Stuttgart, is a fully remote-controlled 0.6 m f/8.17 Ritchey–Chrétien telescope optimized for high-cadence, high-fidelity photometry of transient sources. Observations are time-referenced with very high accuracy and precision, making it an ideal platform for time-domain astronomy and space situational awareness. Initially conceived to support instrument developments and operations of SOFIA, the Stratospheric Observatory for Infrared Astronomy, it evolved into a scientific instrument for various use cases in instrument development, astronomical research, and teaching. This paper presents an overview of its development and optimization to achieve diffraction-limited images and highly accurate pointing and tracking, even at high speeds. The findings and lessons learned are universally applicable to other telescopes that are currently at the planning stage, or where similar issues might be encountered.

Achieving good accuracy in the positions of the ALMA antennas is essential for high-fidelity interferometric observations. Scientific observations use a nearby unresolved phase-reference calibrator, but the finite separation angle from the target of interest, combined with antenna position errors, results in phase errors in the calibrated data of the target, which become worse at higher frequencies, larger calibrator separations, and longer antenna baselines. The ALMA antenna position accuracy goal is to achieve better than 0.5 mm for which a target-calibrator separation of 5^∘ produces a maximum phase noise of 0.5 rad at 869 GHz. In 2021, an antenna position monitoring campaign was made, consisting of 121 observations at 100 GHz for various array configurations including the most extended, with maximum baselines ranging between 3 and 16 km. The findings indicate that the antenna position uncertainty increases relative to the distance from the array center, following a power law with an index of 0.2 for short term observations, and 0.5 for long term monitoring. The uncertainty also grows over time, exhibiting a power law index of 0.3. As part of the error, the antenna positions show a drift over time whose amplitude is proportional to the distance, with a peak of the drift rate of 0.5 mm day ⁻¹ at 8 km. Antennas which are close to one another show a similar drift vector, which changes gradually with the distance, and forms two global drift patterns over the whole array, one parallel and one perpendicular to the ambient wind direction at the Earth’s surface. It is suggested that atmospheric turbulence, landscape features, and wind are the main drivers of the observed errors in antenna position. High-frequency (>300 GHz) interferometry can significantly reduce phase error by performing a short baseline calibration involving three point sources at lower frequencies (100 GHz) at the beginning of each scientific phase experiment.

We present a multi-band study of three symbiotic binaries using combined ground- and space-based monitoring that spans up to 14 yr. These data sets enable a systematic investigation of variability on intermediate timescales (tens of days) and the detection of shorter-period signals. All systems display coherent photometric modulations that are distinct from the orbital cycles. In AX Per, a dominant 75 days signal and its 37 days harmonic are identified, which we interpret as pulsations of the cool giant. CI Cyg exhibits a stable modulation between 70 and 74 days, which likely arises from a combination of pulsation and circumstellar or disk-related variability. For Z And, we confirm a persistent modulation between 55 and 60 days, consistent with semiregular pulsations of the cool component. Additionally, space-based data reveal further short-period variability, including coherent signals at 26.7 and 66.6 minutes in Z And and CI Cyg, respectively, and a quasi-periodic modulation near 0.95 days in AX Per. These detections suggest the presence of rapid activity driven by accretion or rotation, superposed on the intermediate timescale behavior. Our results show that the observed variability in these symbiotic binaries reflects the combined effects of cool giant pulsation, circumstellar or disk activity, and possible rotation of the hot component. The multi-timescale behavior revealed here offers new constraints on mass transfer and activity cycles in interacting binaries.

We perform numerical simulations to investigate high-power wind accretion in massive binary systems undergoing enhanced mass-loss episodes. The primary star is taken in the mass range M₁ = 60–90 M_⊙, while the companion is a 30 M_⊙ hot star. We model binary orbits with eccentricities of e = 0–0.6 and orbital periods of P = 455–1155 days. We initiate strong eruptive events for the primary with mass-loss rates of ${\dot{M}}_{{\rm{w}}}=1{0}^{-2}$–10⁻¹M_⊙ yr⁻¹, lasting for 1.5 yr. A fraction of the ejected wind material is accreted by the companion, with the accretion efficiency determined by the orbital separation, eccentricity, and stellar mass ratio. We analyze the resulting accretion rates and provide an analytical relation describing their dependence on the stellar mass ratio, mass-loss rate, and orbital parameters. We find that although the accretion modifies the stellar parameters of the secondary, the companion remains in thermal equilibrium and does not undergo significant radial expansion. We further include wind mass loss from the companion during wind accretion and find a substantial reduction in accretion efficiency compared to no wind scenario. For longer orbital periods, the models yield negative accretion rates, implying that any captured material is expelled or prevented from settling onto the accretor. These results provide new insight into the role of eccentric orbits and extreme mass-loss events in shaping the mass-transfer processes in massive binaries.

Modern wide-field time-domain surveys facilitate the study of transient, variable and moving phenomena by conducting image differencing and relaying alerts to their communities. Machine learning tools have been used on data from these surveys and their precursors for more than a decade, and convolutional neural networks (CNNs), which make predictions directly from input images, saw particularly broad adoption through the 2010s. Since then, continually rapid advances in computer vision have transformed the standard practices around using such models. It is now commonplace to use standardized architectures pre-trained on large corpora of everyday images (e.g., ImageNet). In contrast, time-domain astronomy studies still typically design custom CNN architectures and train them from scratch. Here, we explore the effects of adopting various pre-training regimens and standardized model architectures on the performance of alert classification. We find that the resulting models match or outperform a custom, specialized CNN like what is typically used for filtering alerts. Moreover, our results show that pre-training on galaxy images from Galaxy Zoo tends to yield better performance than pre-training on ImageNet or training from scratch. We observe that the design of standardized architectures are much better optimized than the custom CNN baseline, requiring significantly less time and memory for inference despite having more trainable parameters. On the eve of the Legacy Survey of Space and Time and other image-differencing surveys, these findings advocate for a paradigm shift in the creation of vision models for alerts, demonstrating that greater performance and efficiency, in time and in data, can be achieved by adopting the latest practices from the computer vision field.

The Barbara A. Mikulski Archive for Space Telescopes (MAST) hosts science-ready data products from over 20 NASA missions plus community-contributed data collections and other select surveys. The data support forefront research in the ultraviolet, optical, and near-infrared wavelength bands. We have constructed bibliographies for each mission from publications in nearly 40 professional journals and identified more than 37,000 refereed articles where investigators made a science usage of data hosted in MAST. The publication rate over the last 50 yr shows that most MAST missions have had very high productivity during their in-service lifetimes and have remained so for years or decades afterward. Annual citations of these publications, a measure of impact on research, are robust for most missions, with citations that grow over more than a decade. Most of the citations come from about 10% of the articles within each mission. We examined the bibliographies of the active missions Hubble Space Telescope (HST) and James Webb Space Telescope (JWST) in greater detail. For HST, the rate of archival publications exceeded those authored by the original observing teams within a decade of launch and is now more than 3 times higher. Early indications hint that JWST archival articles could dominate the publication rate even sooner. The production of articles resulting from any given observing program can extend for decades. Programs with small and very large allocations of observing time tend to be particularly productive per unit of observing time. For HST in general, a first publication appears within 1.5 yr for 50% of observing programs and within 3.8 yr for 80% of programs. We discuss various external factors that affect publication metrics, their strengths and limitations for measuring scientific impact, and the challenges of making meaningful comparisons of publication metrics across missions.

We introduce a stable, well tested Python implementation of the affine-invariant ensemble sampler for Markov chain Monte Carlo (MCMC) proposed by Goodman & Weare (2010). The code is open source and has already been used in several published projects in the astrophysics literature. The algorithm behind emcee has several advantages over traditional MCMC sampling methods and it has excellent performance as measured by the autocorrelation time (or function calls per independent sample). One major advantage of the algorithm is that it requires hand-tuning of only 1 or 2 parameters compared to ∼N² for a traditional algorithm in an N-dimensional parameter space. In this document, we describe the algorithm and the details of our implementation. Exploiting the parallelism of the ensemble method, emcee permits any user to take advantage of multiple CPU cores without extra effort. The code is available online at http://dan.iel.fm/emcee under the GNU General Public License v2.

We review recent determinations of the present‐day mass function (PDMF) and initial mass function (IMF) in various components of the Galaxy—disk, spheroid, young, and globular clusters—and in conditions characteristic of early star formation. As a general feature, the IMF is found to depend weakly on the environment and to be well described by a power‐law form for m≳1 M_⊙ and a lognormal form below, except possibly for early star formation conditions. The disk IMF for single objects has a characteristic mass around m_c ∼ 0.08 M_⊙ and a variance in logarithmic mass σ ∼ 0.7, whereas the IMF for multiple systems has m_c ∼ 0.2 M_⊙ and σ ∼ 0.6. The extension of the single MF into the brown dwarf regime is in good agreement with present estimates of L‐ and T‐dwarf densities and yields a disk brown dwarf number density comparable to the stellar one, n_BD ∼ n_* ∼ 0.1 pc⁻³. The IMF of young clusters is found to be consistent with the disk field IMF, providing the same correction for unresolved binaries, confirming the fact that young star clusters and disk field stars represent the same stellar population. Dynamical effects, yielding depletion of the lowest mass objects, are found to become consequential for ages ≳130 Myr. The spheroid IMF relies on much less robust grounds. The large metallicity spread in the local subdwarf photometric sample, in particular, remains puzzling. Recent observations suggest that there is a continuous kinematic shear between the thick‐disk population, present in local samples, and the genuine spheroid one. This enables us to derive only an upper limit for the spheroid mass density and IMF. Within all the uncertainties, the latter is found to be similar to the one derived for globular clusters and is well represented also by a lognormal form with a characteristic mass slightly larger than for the disk, m_c ∼ 0.2–0.3 M_⊙, excluding a significant population of brown dwarfs in globular clusters and in the spheroid. The IMF characteristic of early star formation at large redshift remains undetermined, but different observational constraints suggest that it does not extend below ∼1 M_⊙. These results suggest a characteristic mass for star formation that decreases with time, from conditions prevailing at large redshift to conditions characteristic of the spheroid (or thick disk) to present‐day conditions. These conclusions, however, remain speculative, given the large uncertainties in the spheroid and early star IMF determinations.

These IMFs allow a reasonably robust determination of the Galactic present‐day and initial stellar and brown dwarf contents. They also have important galactic implications beyond the Milky Way in yielding more accurate mass‐to‐light ratio determinations. The mass‐to‐light ratios obtained with the disk and the spheroid IMF yield values 1.8–1.4 times smaller than for a Salpeter IMF, respectively, in agreement with various recent dynamical determinations. This general IMF determination is examined in the context of star formation theory. None of the theories based on a Jeans‐type mechanism, where fragmentation is due only to gravity, can fulfill all the observational constraints on star formation and predict a large number of substellar objects. On the other hand, recent numerical simulations of compressible turbulence, in particular in super‐Alfvénic conditions, seem to reproduce both qualitatively and quantitatively the stellar and substellar IMF and thus provide an appealing theoretical foundation. In this picture, star formation is induced by the dissipation of large‐scale turbulence to smaller scales through radiative MHD shocks, producing filamentary structures. These shocks produce local nonequilibrium structures with large density contrasts, which collapse eventually in gravitationally bound objects under the combined influence of turbulence and gravity. The concept of a single Jeans mass is replaced by a distribution of local Jeans masses, representative of the lognormal probability density function of the turbulent gas. Objects below the mean thermal Jeans mass still have a possibility to collapse, although with a decreasing probability.

CASA, the Common Astronomy Software Applications, is the primary data processing software for the Atacama Large Millimeter/submillimeter Array (ALMA) and the Karl G. Jansky Very Large Array (VLA), and is frequently used also for other radio telescopes. The CASA software can handle data from single-dish, aperture-synthesis, and Very Long Baseline Interferometery (VLBI) telescopes. One of its core functionalities is to support the calibration and imaging pipelines for ALMA, VLA, VLA Sky Survey, and the Nobeyama 45 m telescope. This paper presents a high-level overview of the basic structure of the CASA software, as well as procedures for calibrating and imaging astronomical radio data in CASA. CASA is being developed by an international consortium of scientists and software engineers based at the National Radio Astronomy Observatory (NRAO), the European Southern Observatory, the National Astronomical Observatory of Japan, and the Joint Institute for VLBI European Research Infrastructure Consortium (JIV-ERIC), under the guidance of NRAO.

We present a novel algorithm for scheduling the observations of time-domain imaging surveys. Our integer linear programming approach optimizes an observing plan for an entire night by assigning targets to temporal blocks, enabling strict control of the number of exposures obtained per field and minimizing filter changes. A subsequent optimization step minimizes slew times between each observation. Our optimization metric self-consistently weights contributions from time-varying airmass, seeing, and sky brightness to maximize the transient discovery rate. We describe the implementation of this algorithm on the surveys of the Zwicky Transient Facility and present its on-sky performance.

The Zwicky Transient Facility (ZTF) is a new robotic time-domain survey currently in progress using the Palomar 48-inch Schmidt Telescope. ZTF uses a 47 square degree field with a 600 megapixel camera to scan the entire northern visible sky at rates of ∼3760 square degrees/hour to median depths of g ∼ 20.8 and r ∼ 20.6 mag (AB, 5σ in 30 sec). We describe the Science Data System that is housed at IPAC, Caltech. This comprises the data-processing pipelines, alert production system, data archive, and user interfaces for accessing and analyzing the products. The real-time pipeline employs a novel image-differencing algorithm, optimized for the detection of point-source transient events. These events are vetted for reliability using a machine-learned classifier and combined with contextual information to generate data-rich alert packets. The packets become available for distribution typically within 13 minutes (95th percentile) of observation. Detected events are also linked to generate candidate moving-object tracks using a novel algorithm. Objects that move fast enough to streak in the individual exposures are also extracted and vetted. We present some preliminary results of the calibration performance delivered by the real-time pipeline. The reconstructed astrometric accuracy per science image with respect to Gaia DR1 is typically 45 to 85 milliarcsec. This is the RMS per-axis on the sky for sources extracted with photometric S/N ≥ 10 and hence corresponds to the typical astrometric uncertainty down to this limit. The derived photometric precision (repeatability) at bright unsaturated fluxes varies between 8 and 25 millimag. The high end of these ranges corresponds to an airmass approaching ∼2—the limit of the public survey. Photometric calibration accuracy with respect to Pan-STARRS1 is generally better than 2%. The products support a broad range of scientific applications: fast and young supernovae; rare flux transients; variable stars; eclipsing binaries; variability from active galactic nuclei; counterparts to gravitational wave sources; a more complete census of Type Ia supernovae; and solar-system objects.

Technology has advanced to the point that it is possible to image the entire sky every night and process the data in real time. The sky is hardly static: many interesting phenomena occur, including variable stationary objects such as stars or QSOs, transient stationary objects such as supernovae or M dwarf flares, and moving objects such as asteroids and the stars themselves. Funded by NASA, we have designed and built a sky survey system for the purpose of finding dangerous near-Earth asteroids (NEAs). This system, the “Asteroid Terrestrial-impact Last Alert System” (ATLAS), has been optimized to produce the best survey capability per unit cost, and therefore is an efficient and competitive system for finding potentially hazardous asteroids (PHAs) but also for tracking variables and finding transients. While carrying out its NASA mission, ATLAS now discovers more bright (m < 19) supernovae candidates than any ground based survey, frequently detecting very young explosions due to its 2 day cadence. ATLAS discovered the afterglow of a gamma-ray burst independent of the high energy trigger and has released a variable star catalog of 5 × 10⁶ sources. This is the first of a series of articles describing ATLAS, devoted to the design and performance of the ATLAS system. Subsequent articles will describe in more detail the software, the survey strategy, ATLAS-derived NEA population statistics, transient detections, and the first data release of variable stars and transient light curves.

I introduce batman, a Python package for modeling exoplanet transit and eclipse light curves. The batman package supports calculation of light curves for any radially symmetric stellar limb darkening law, using a new integration algorithm for models that cannot be quickly calculated analytically. The code uses C extension modules to speed up model calculation and is parallelized with OpenMP. For a typical light curve with 100 data points in transit, batman can calculate one million quadratic limb-darkened models in 30 s with a single 1.7 GHz Intel Core i5 processor. The same calculation takes seven minutes using the four-parameter nonlinear limb darkening model (computed to 1 ppm accuracy). Maximum truncation error for integrated models is an input parameter that can be set as low as 0.001 ppm, ensuring that the community is prepared for the precise transit light curves we anticipate measuring with upcoming facilities. The batman package is open source and publicly available at https://github.com/lkreidberg/batman.

The All-Sky Automated Survey for Supernovae (ASAS-SN) is working toward imaging the entire visible sky every night to a depth of $V\sim 17$ mag. The present data covers the sky and spans ∼2–5 years with ∼100–400 epochs of observation. The data should contain some ∼1 million variable sources, and the ultimate goal is to have a database of these observations publicly accessible. We describe here a first step, a simple but unprecedented web interface https://asas-sn.osu.edu/ that provides an up to date aperture photometry light curve for any user-selected sky coordinate. The V band photometry is obtained using a two-pixel (160) radius aperture and is calibrated against the APASS catalog. Because the light curves are produced in real time, this web tool is relatively slow and can only be used for small samples of objects. However, it also imposes no selection bias on the part of the ASAS-SN team, allowing the user to obtain a light curve for any point on the celestial sphere. We present the tool, describe its capabilities, limitations, and known issues, and provide a few illustrative examples.

This paper characterizes the actual science performance of the James Webb Space Telescope (JWST), as determined from the six month commissioning period. We summarize the performance of the spacecraft, telescope, science instruments, and ground system, with an emphasis on differences from pre-launch expectations. Commissioning has made clear that JWST is fully capable of achieving the discoveries for which it was built. Moreover, almost across the board, the science performance of JWST is better than expected; in most cases, JWST will go deeper faster than expected. The telescope and instrument suite have demonstrated the sensitivity, stability, image quality, and spectral range that are necessary to transform our understanding of the cosmos through observations spanning from near-earth asteroids to the most distant galaxies.

The K2 mission will make use of the Kepler spacecraft and its assets to expand upon Kepler’s groundbreaking discoveries in the fields of exoplanets and astrophysics through new and exciting observations. K2 will use an innovative way of operating the spacecraft to observe target fields along the ecliptic for the next 2–3 years. Early science commissioning observations have shown an estimated photometric precision near 400 ppm in a single 30 minute observation, and a 6-hr photometric precision of 80 ppm (both at V = 12). The K2 mission offers long-term, simultaneous optical observation of thousands of objects at a precision far better than is achievable from ground-based telescopes. Ecliptic fields will be observed for approximately 75 days enabling a unique exoplanet survey which fills the gaps in duration and sensitivity between the Kepler and TESS missions, and offers pre-launch exoplanet target identification for JWST transit spectroscopy. Astrophysics observations with K2 will include studies of young open clusters, bright stars, galaxies, supernovae, and asteroseismology.