Aster G. Taylor

2026

1. The Astronomy of Halakhic Nightfall: Calculating Ts’eit HaKokhavim and Motsa’ei Shabbat

   Aster G. Taylor

   arXiv e-printsAug 2026

   Abs ADS arXiv

   0 citations

   Updated Aug 11, 2026

   In Jewish tradition, the boundary between days is not midnight, but nightfall. Nightfall is when one must pray ma’ariv, count the omer, and may engage in activity after Shabbat. While there are slightly different definitions of nightfall for Shabbat restrictions (Motsa’ei Shabbat) versus other purposes (Ts’eit HaKokhavim), most minhagim define nightfall to be when three stars of some specified size and proximity are visible in the sky. This definition presents some difficulty. Not only are these conditions difficult to define, they are further complicated by weather, cloud cover, and light pollution. Observant Jews therefore usually approximate these times by defining nightfall to be when the sun has reached a certain distance below the horizon. Although these approximations have been effective, modern astronomy enables us to calculate precisely when these conditions are met. This work applies these techniques to calculate the time of Ts’eit HaKokhavim and Motsa’ei Shabbat evaluates the accuracy of the approximations, and explores the effects of light pollution on these times. While Ts’eit HaKokhavim is reasonably well-approximated by standard methods, the conditions for Motsa’ei Shabbat generally occur after the time predicted by the approximations. Light pollution does not shift Ts’eit HaKokhavim by more than a minute but significantly changes Motsa’ei Shabbat. In fact, light pollution causes the conditions of Motsa’ei Shabbat to never be met on at least some nights in every population center. Finally, I provide a tool (halakhic-nightfall.streamlit.app) for calculating Ts’eit HaKokhavim and Motsa’ei Shabbat at an arbitrary time and location on Earth.

2. The Effects of Magnetic Accretion on the Spatial Extent of White Dwarf Pollution

   Dang Pham, Aster G. Taylor, and Tim Cunningham

   arXiv e-printsJul 2026

   Abs ADS arXiv

   0 citations

   Updated Aug 11, 2026

   Many white dwarfs are polluted by metals, which are generally understood to be the accreted remnants of a planetary system. Modeling these systems typically assumes that the metal concentration is homogeneous throughout the white dwarf’s atmosphere. However, the magnetic fields of a white dwarf may affect the accretion geometry of the white dwarf via magnetospheric accretion. Convection in the white dwarf’s photosphere will then transport the metals across the surface, with a structure set by the relative sinking versus spreading timescales. In this work, we construct models for the accretion geometry, subsequent spreading, and observed pollution of magnetic white dwarfs. We show that the magnetic fields will initially concentrate the pollution into a narrow region of the white dwarf’s surface. The relative spreading and sinking timescales determine whether the metals become uniformly distributed or remain confined to localized patches. If the magnetic field and spin poles are misaligned, then patchy white dwarfs exhibit periodically variable pollution signatures, which enable constraints on the patch area. We explore this model as a possible explanation for the recent detection of periodically variable pollution signatures in magnetic white dwarfs. Finally, we also demonstrate that the concentration of material due to the magnetic field may lead to systematic underestimates of the mass accretion rate onto these objects.

3. A retrieval framework for observationally constraining the parameters of circumplanetary disks

   Aster G. Taylor, and Fred C. Adams

   IcarusDec 2026

   Abs ADS arXiv

   0 citations

   Updated Aug 11, 2026

   As they form, giant planets are surrounded by disks of gas and dust sourced from the background circumstellar disk. Although there have been few detections to date, upcoming instruments are likely to discover many more of these systems in the coming decades. Accurate spectral modeling will enable these observations to constrain the properties of these forming systems. Towards this end, we have constructed a semianalytic model for the structure and radiative signatures of geometrically thick circumplanetary disks and their planet hosts. Fitting these radiative signatures to synthetic observations of a two-dimensional disk model then quantifies the parameter constraints that can be derived (subject to model assumptions). This machinery provides estimates of the values and uncertainties in system parameters, and some combinations of parameters have significantly smaller uncertainties than others. This model is then used to fit observations of real protoplanets, with good results. The derived parameters provide useful context about the local extinction, formation history, and initial entropy of these objects.

4. The Two-Dimensional Structure of Circumplanetary Disks and their Radiative Signatures

   Aster G. Taylor, Fred C. Adams, and Nuria Calvet

   IcarusMar 2026

   Abs ADS arXiv

   1 citation

   Updated Aug 11, 2026

   During their formative stages, giant planets are fed by infalling material sourced from the background circumstellar disk. Due to conservation of angular momentum, the incoming gas and dust collects into a circumplanetary disk that processes the material before it reaches the central planet itself. This work investigates the complex vertical structure of these circumplanetary disks and calculates their radiative signatures. A self-consistent numerical model of the temperature and density structure of the circumplanetary environment reveals that circumplanetary disks are thick and hot, with aspect ratios H/R∼0.1−0.25 and temperatures approaching that of the central planet. The disk geometry has a significant impact on the radiative signatures, allowing future observations to determine critical system parameters. The resulting disks are gravitationally stable and viscosity is sufficient to drive the necessary disk accretion. However, sufficiently rapid mass accretion can trigger a thermal instability, which sets an upper limit on the mass accretion rate. This paper shows how the radiative signatures depend on the properties of the planetary system and discuss how the system parameters can be constrained by future observations.

2025

1. The Kinematic Age of 3I/ATLAS and Its Implications for Early Planet Formation

   Aster G. Taylor, and Darryl Z. Seligman

   The Astrophysical Journal LettersSep 2025

   Abs ADS arXiv

   45 citations

   Updated Aug 11, 2026

   The recent discovery of the third interstellar object (3I/ATLAS) expands the known census from two to three and significantly improves statistical inferences regarding the underlying galactic population. 3I/ATLAS exhibits detectable activity that may increase as the object approaches perihelion. In this paper, we argue that the cometary activity likely significantly contributes to 3I/ATLAS’s brightness, since the nuclear size of 3I/ATLAS assuming an asteroidal reflectance implies an untenable interstellar object mass per star. 3I/ATLAS exhibits an excess velocity of v∞=58 km/s relative to the Sun which is significantly larger than that of 1I/‘Oumuamua or 2I/Borisov. This velocity implies that 3I/ATLAS is relatively old in comparison to previous interstellar objects. Here, we calculate the posterior distribution of ages implied by the kinematics of the interstellar objects and find that 3I/ATLAS is likely ∼3−11 Gyr old, assuming that the interstellar object and stellar age-velocity relations are equivalent. We also calculate the distribution of host star metallicities and find that 3I/ATLAS likely originated in a lower-metallicity system than the previous interstellar objects. These results show that interstellar object formation is likely efficient at low metallicities and early in the history of the Galaxy. Finally, we estimate and approximate interstellar object formation rate throughout Galactic history implied by these three objects.

2024

1. Strong Nongravitational Accelerations and the Potential for Misidentification of Near-Earth Objects

   Aster G. Taylor, Darryl Z. Seligman, Matthew J. Holman, and 5 more authors

   The Astrophysical JournalNov 2024

   Abs ADS arXiv

   3 citations

   Updated Aug 11, 2026

   Nongravitational accelerations in the absence of observed activity have recently been identified on near-Earth objects (NEOs), opening the question of the prevalence of anisotropic mass loss in the near-Earth environment. Motivated by the necessity of nongravitational accelerations to identify 2010 VL65 and 2021 UA12 as a single object, we investigate the problem of linking separate apparitions in the presence of nongravitational perturbations. We find that nongravitational accelerations on the order of 1 × 10–9 au day−2 can lead to a change in plane-of-sky positions of ∼1 × 103 arcsec between apparitions. Moreover, we inject synthetic tracklets of hypothetical nongravitationally accelerating NEOs into the Minor Planet Center orbit identification algorithms. We find that at large nongravitational accelerations (∣A i ∣ ≥ 1 × 10−8 au day−2) these algorithms fail to link a significant fraction of these tracklets. We further show that if orbits can be determined for both apparitions, the tracklets will be linked regardless of nongravitational accelerations, although they may be linked to multiple objects. In order to aid in the identification and linkage of nongravitationally accelerating objects, we propose and test a new methodology to search for unlinked pairs. When applied to the current census of NEOs, we recover the previously identified case but identify no new linkages. We conclude that current linking algorithms are generally robust to nongravitational accelerations, but objects with large nongravitational accelerations may potentially be missed. While current algorithms are well-positioned for the anticipated increase in the census population from future survey missions, it may be possible to find objects with large nongravitational accelerations hidden in isolated tracklet pairs.

2. Radiative Signatures of Circumplanetary Disks and Envelopes During the Late Stages of Giant Planet Formation

   Aster G. Taylor, and Fred C. Adams

   IcarusNov 2024

   Abs ADS arXiv

   6 citations

   Updated Aug 11, 2026

   During the late stages of giant planet formation, protoplanets are surrounded by a circumplanetary disk and an infalling envelope of gas and dust. For systems with sufficient cooling, material entering the sphere of influence of the planet falls inward and approaches ballistic conditions. Due to conservation of angular momentum, most of the incoming material falls onto the disk rather than directly onto the planet. This paper determines the spectral energy distributions of forming planets in this stage of evolution. Generalizing previous work, we consider a range of possible geometries for the boundary conditions of the infall and determine the two-dimensional structure of the envelope, as well as the surface density of the disk. After specifying the luminosity sources for the planet and disk, we calculate the corresponding radiative signatures for the system, including the emergent spectral energy distributions and emission maps. These results show how the observational appearance of forming planets depend on the input parameters, including the instantaneous mass, mass accretion rate, semimajor axis of the orbit, and the planetary magnetic field strength (which sets the inner boundary condition for the disk). We also consider different choices for the form of the opacity law and attenuation due to the background circumstellar disk. Although observing forming planets will be challenging, these results show how the observational signatures depend on the underlying properties of the planet/disk/envelope system.

3. The Dynamical Origins of the Dark Comets and a Proposed Evolutionary Track

   Aster G. Taylor, Jordan K. Steckloff, Darryl Z. Seligman, and 5 more authors

   IcarusJul 2024

   Abs ADS arXiv

   9 citations

   Updated Aug 11, 2026

   So-called ‘dark comets’ are small, morphologically inactive near-Earth objects (NEOs) that exhibit nongravitational accelerations inconsistent with radiative effects. These objects exhibit short rotational periods (minutes to hours), where measured. We find that the strengths required to prevent catastrophic disintegration are consistent with those measured in cometary nuclei and expected in rubble pile objects. We hypothesize that these dark comets are the end result of a rotational fragmentation cascade, which is consistent with their measured physical properties. We calculate the predicted size-frequency distribution for objects evolving under this model. Using dynamical simulations, we further demonstrate that the majority of these bodies originated from the ν6 resonance, implying the existence of volatiles in the current inner main belt. Moreover, one of the dark comets, (523599) 2003 RM, likely originated from the outer main belt, although a JFC origin is also plausible. These results provide strong evidence that volatiles from a reservoir in the inner main belt are present in the near-Earth environment.

4. Formation and Structure of Circumplanetary Disks and Envelopes During the Late Stages of Giant Planet Formation

   Aster G. Taylor, and Fred C. Adams

   IcarusJun 2024

   Abs ADS arXiv

   9 citations

   Updated Aug 11, 2026

   Giant planets are expected to form within circumplanetary disks, which shape their formation history and the local environment. Here, we consider the formation and structure of circumplanetary disks that arise during the late stages of giant planet formation. During this phase, when most of the final mass is accumulated, incoming material enters the Hill sphere and falls toward the planet. In the absence of torques, the falling parcels of gas conserve their specific angular momentum and collect into a circumplanetary disk. Generalizing previous work, we consider a range of possible geometries for the flow entering the sphere of influence of the planet. Specifically, we consider five geometric patterns for the inward flow, ranging from concentration toward the rotational poles of the system to isotropic flow to concentration along the equatorial plane. For each case, we derive analytic descriptions for the density field of the infall region, the disk surface density in the absence of viscosity, and steady-state solutions for viscous disks. These results, in turn, specify the luminosity contributions of the planet, the circumplanetary disk, and the envelope. These power sources, in conjunction with the surrounding material, collectively determine the observational appearance of the forming planet. We conclude with an approximate determination of these radiative signatures.

5. Seasonally Varying Outgassing as an Explanation for Dark Comet Accelerations

   Aster G. Taylor, Davide Farnocchia, David Vokrouhlický, and 3 more authors

   IcarusJan 2024

   Abs ADS arXiv

   5 citations

   Updated Aug 11, 2026

   Significant nonradial, nongravitational accelerations with magnitudes incompatible with radiation-driven effects have been reported in seven photometrically inactive near-Earth objects. Two of these objects exhibit large transverse accelerations (i.e., within the orbital plane but orthogonal to the radial direction), and six exhibit significant out-of-plane accelerations. Here, we find that anisotropic outgassing resulting from differential heating on a nucleus with nonzero spin-pole obliquity, averaged over an eccentric orbit, can explain these accelerations for most of the objects. This balanced outgassing model depends on three parameters — the spin pole orientation (R.A. and Dec.) and an acceleration magnitude. For these “dark comets” (excepting 2003 RM), we obtain parameter values that reproduce the observed nongravitational accelerations. We derive formulae for the component accelerations under certain assumptions for the acceleration scaling over heliocentric distance. Although we lack estimates of these objects’ spin axes to confirm our values, this mechanism is nevertheless a plausible explanation for the observed accelerations, and produces accurate perturbations to the heliocentric motions of most of these objects. This model may also be applied to active objects outside of the dark comets group.

2023

1. Fitting the Light Curve of 1I/’Oumuamua with a Nonprincipal Axis Rotational Model and Outgassing Torques

   Aster G. Taylor, Darryl Z. Seligman, Olivier R. Hainaut, and 1 more author

   Planetary Science JournalOct 2023

   Abs ADS arXiv

   4 citations

   Updated Aug 11, 2026

   In this paper, we investigate the nonprincipal axis (NPA) rotational state of 1I/‘Oumuamua—the first interstellar object discovered traversing the inner solar system—from its photometric light curve. Building upon Mashchenko, we develop a model which incorporates NPA rotation and Sun-induced, time-varying outgassing torques to generate synthetic light curves of the object. The model neglects tidal forces, which are negligible compared to outgassing torques over the distances at which ‘Oumuamua was observed. We implement an optimization scheme that incorporates the NPA rotation model to calculate the initial rotation state of the object. We find that an NPA rotation state with an average period of 〈P〉 ≃ 7.34 hr best reproduces the photometric data. The discrepancy between this period and previous estimates is due to continuous period modulation induced by outgassing torques in the rotational model, as well as different periods being used. The best fit to the 2017 October data does not reproduce the 2017 November data (although the later measurements are too sparse to fit). The light curve is consistent with there being no secular evolution of the angular momentum, which is somewhat in tension with the empirical correlations between nuclear spin-up and cometary outgassing. The complex rotation of ‘Oumuamua may be the result of primordial rotation about the smallest principal axis if (i) the object experienced hypervolatile outgassing and (ii) our idealized outgassing model is accurate.

2. Interstellar Comets from Post–Main Sequence Systems as Tracers of Extrasolar Oort Clouds

   W. Garrett Levine, Aster G. Taylor, Darryl Z. Seligman, and 4 more authors

   Planetary Science JournalJul 2023

   Abs ADS arXiv

   13 citations

   Updated Aug 11, 2026

   Interstellar small bodies are unique probes into the histories of exoplanetary systems. One hypothesized class of interlopers are "Jurads," exocomets released into the Milky Way during the post-main-sequence as the thermally pulsing asymptotic giant branch (AGB) host stars lose mass. In this study, we assess the prospects for the Legacy Survey of Space and Time (LSST) to detect a Jurad and examine whether such an interloper would be observationally distinguishable from exocomets ejected during the (pre-)main-sequence. Using analytic and numerical methods, we estimate the fraction of exo-Oort Cloud objects that are released from 1-8 M ⊙ stars during post-main-sequence evolution. We quantify the extent to which small bodies are altered by the increased luminosity and stellar outflows during the AGB, finding that some Jurads may lack hypervolatiles and that stellar winds could deposit dust that covers the entire exocomet surface. Next, we construct models of the interstellar small body reservoir for various size-frequency distributions and examine the LSST’s ability to detect members of those hypothesized populations. Combining these analyses, we highlight the joint constraints that the LSST will place on power-law size-frequency distribution slopes, characteristic sizes, and the total mass sequestered in the minor planets of exo-Oort Clouds. Even with the LSST’s increased search volume compared to contemporary surveys, we find that detecting a Jurad is unlikely but not infeasible given the current understanding of (exo)planet formation.

3. Dark Comets? Unexpectedly Large Nongravitational Accelerations on a Sample of Small Asteroids

   Darryl Z. Seligman, Davide Farnocchia, Marco Micheli, and 15 more authors

   Planetary Science JournalFeb 2023

   Abs ADS arXiv

   43 citations

   Updated Aug 11, 2026

   We report statistically significant detections of nonradial, nongravitational accelerations based on astrometric data in the photometrically inactive objects 1998 KY26, 2005 VL1, 2016 NJ33, 2010 VL65, 2016 RH120, and 2010 RF12. The magnitudes of the nongravitational accelerations are greater than those typically induced by the Yarkovsky effect, and there is no radiation-based, nonradial effect that can be so large. Therefore, we hypothesize that the accelerations are driven by outgassing and calculate implied H2O production rates for each object. We attempt to reconcile outgassing-induced acceleration with the lack of visible comae or photometric activity via the absence of surface dust and low levels of gas production. Although these objects are small, and some are rapidly rotating, the surface cohesive forces are stronger than the rotational forces, and rapid rotation alone cannot explain the lack of surface debris. It is possible that surface dust was removed previously, perhaps via outgassing activity that increased the rotation rates to their present-day value. We calculate dust production rates of order  10-4 g s-1 in each object, assuming that the nuclei are bare, within the upper limits of dust production from a sample stacked image of 1998 KY26 of M˙Dust<0.2 g s-1. This production corresponds to brightness variations of order  0.0025%, which are undetectable in extant photometric data. We assess the future observability of each of these targets and find that the orbit of 1998 KY26-which is also the target of the extended Hayabusa2 mission-exhibits favorable viewing geometry before 2025.