Lowell Observatory Main Page

Research involving Lowell Observatory staff 2021
(Articles and chapters)

This is a work ever in progress.

(Pulled from ADS* by sel on 2024-03-25)

*We are grateful for all the effort that went into making The SAO/NASA Astrophysics Data System (ADS) possible. The ADS is operated by the Smithsonian Astrophysical Observatory under NASA Cooperative Agreement NNX16AC86A and can be found at: https://ui.adsabs.harvard.edu/

If you notice publications that are missing, or ones that do not belong, please let us know (send email to sel .at. lowell .dot. edu).

For missing articles, please send either the ADS bibcode, or a standard short form journal citation.

Authors affiliated with "Lowell Obs" are highlighted.

Return to Lowell Annual Summaries

Years: 2021 Bottom

    2021

  1. Shipp, N., Erkal, D., Drlica-Wagner, A., Li, T., Pace, A., Koposov, S., Cullinane, L., Da Costa, G., Ji, A., Kuehn, K., Lewis, G., Mackey, D., Simpson, J., Wan, Z., Zucker, D., Bland-Hawthorn, J., Ferguson, P., Lilleengen, S., Lilleengen, S., 2021, ApJ, 923, 149, Measuring the Mass of the Large Magellanic Cloud with Stellar Streams Observed by S 5
    Stellar streams are excellent probes of the underlying gravitational potential in which they evolve. In this work, we fit dynamical models to five streams in the Southern Galactic hemisphere, combining observations from the Southern Stellar Stream Spectroscopic Survey (S 5), Gaia EDR3, and the Dark Energy Survey, to measure the mass of the Large Magellanic Cloud (LMC). With an ensemble of streams, we find a mass of the LMC ranging from ~14-19 1010 M , probed over a range of closest approach times and distances. With the most constraining stream (Orphan-Chenab), we measure an LMC mass of ${18.8}_{-4.0}^{+3.5}\times {10}^{10}\,{M}_{\odot }$ , probed at a closest approach time of 310 Myr and a closest approach distance of 25.4 kpc. This mass is compatible with previous measurements, showing that a consistent picture is emerging of the LMC's influence on structures in the Milky Way. Using this sample of streams, we find that the LMC's effect depends on the relative orientation of the stream and LMC at their point of closest approach. To better understand this, we present a simple model based on the impulse approximation and we show that the LMC's effect depends both on the magnitude of the velocity kick imparted to the stream and the direction of this kick.
  2. Kwiatkowski, T., Kolenczuk, P., Kryszczynska, A., Oszkiewicz, D., Kaminski, K., Kaminska, M., Troianskyi, V., Skiff, B., Moskowitz, N., Kashuba, V., Kim, M., Kim, T., Mottola, S., Santana-Ros, T., Kluwak, T., Buzzi, L., Bacci, P., Birtwhistle, P., Miles, R., Chatelain, J., 2021, A&A, 656, A126, Photometry and model of near-Earth asteroid 2021 DW1 from one apparition

    Aims: Very small asteroids (VSAs, objects with diameters smaller than about 150 m) can be spun up by the YORP effect to rotation periods as short as tens of seconds. This effect has been observed for many of them. It is also hypothesised, that in the same process their spin axes are asymptotically drawn to the position perpendicular to the orbital plane. So far this effect has been observed only for one VSA and needs further verification. For that, spin axes of several other VSAs should be determined by observing their brightness variations at many different positions on the sky.
    Methods: On 4 March 2021 at 9 UTC a 30-m in diameter near-Earth asteroid 2021 DW1 passed the Earth at a distance of 570 000 km, reaching the maximum brightness of V = 14.6 mag. We observed it photometrically from 2 March, when it was visible at V = 16.5 mag, until 7 March (V = 18.2 mag). During that time 2021 DW1 swept a 170 long arc in the northern sky, spanning solar phase angles in the range from 36 to 86. This made it an excellent target for physical characterisation, including spin axis and shape derivation.
    Results: Convex inversion of the asteroid lightcurves gives a sidereal period of rotation Psid = 0.013760 0.000001 h, and two solutions for the spin axis ecliptic coordinates: (A) 1 = 57 10, 1 = 29 10 and (B) 2 = 67 10, 2 = 40 10. The magnitude-phase curve can be fitted with a standard H, G function with H = 24.8 0.5 mag and an assumed G = 0.24. The asteroid colour indices are g i = 0.79 0.01 mag, and i z = 0.01 0.02 mag which indicates an S taxonomic class, with an average geometric albedo pV = 0.23 0.02. The asteroid effective diameter, derived from H and pV, is Deff = 30 10 m.
    Conclusions: It was found that the inclination of the spin axis of 2021 DW1 is not perpendicular to the orbital plane (obliquity = 54 10 or = 123 10). More spin axes of VSAs should be determined to check, if 2021 DW1 is an exception or a typical case.
  3. Avdellidou, C., Pajola, M., Lucchetti, A., Agostini, L., Delbo, M., Mazzotta Epifani, E., Bourdelle de Micas, J., Devogele, M., Fornasier, S., van Belle, G., Bruot, N., Dotto, E., Ieva, S., Cremonese, G., Palumbo, P., 2021, A&A, 656, L18, Characterisation of the main belt asteroid (223) Rosa. A proposed flyby target of ESA's JUICE mission
    Context. The ESA JUICE space mission, on its way to study Jupiter's environment and icy moons, will pass twice through the main asteroid belt. For this reason, the possibility to perform an asteroid flyby has been investigated.
    Aims: We aim to gain insight into the physical properties of the outer main belt asteroid (223) Rosa, which has been proposed as a potential JUICE flyby target.
    Methods: We report new visible and near-infrared spectroscopic observations at different rotation phases. Additionally, we perform a literature review of all the available physical properties, such as diameter, albedo, mass, and rotational period.
    Results: We find that asteroid Rosa is an X-type asteroid that shows no significant spectral variability combining the new and literature spectroscopic data. Its large size and orbital semimajor axis in the outer main belt indicate that Rosa does not belong to the Themis family, while its albedo is only marginally compatible with the family. Rosa's estimated density is in agreement with those of other low-albedo X-type asteroids. Hence, we propose that Rosa is a planetesimal that accreted in the protoplanetary disk beyond the snow line.
  4. Clark, C., van Belle, G., Horch, E., 2021, RNAAS, 5, 280, A New Stellar Companion to TYC 5493-889-1
    We present a serendipitous discovery of a new stellar companion to TYC 5493-889-1 detected with the Differential Speckle Survey Instrument at the 4.3 m Lowell Discovery Telescope. We also present photometric observations of TYC 5493-889-1, and determine a spectral type of F1V and a photometric distance of roughly 320 parsecs.
  5. Ferrero, I., Crocce, M., Tutusaus, I., Porredon, A., Blot, L., Fosalba, P., Carnero Rosell, A., Avila, S., Izard, A., Elvin-Poole, J., Chan, K., Camacho, H., Rosenfeld, R., Sanchez, E., Tallada-Crespi, P., Carretero, J., Sevilla-Noarbe, I., Gaztanaga, E., Andrade-Oliveira, F., De Vicente, J., Mena-Fernandez, J., Ross, A., Sanchez Cid, D., Ferte, A., Brandao-Souza, A., Fang, X., Krause, E., Gomes, D., Aguena, M., Allam, S., Annis, J., Bertin, E., Brooks, D., Carrasco Kind, M., Castander, F., Cawthon, R., Choi, A., Conselice, C., Costanzi, M., da Costa, L., Pereira, M., Diehl, H., Doel, P., Drlica-Wagner, A., Everett, S., Evrard, A., Flaugher, B., Frieman, J., Garcia-Bellido, J., Gerdes, D., Gruen, D., Gruendl, R., Gschwend, J., Gutierrez, G., Hinton, S., Hollowood, D., Honscheid, K., Hoyle, B., Huterer, D., James, D., Kuehn, K., Lima, M., Maia, M., Marshall, J., Menanteau, F., Miquel, R., Morgan, R., Muir, J., Ogando, R., Palmese, A., Paz-Chinchon, F., Percival, W., Plazas Malagon, A., Rodriguez-Monroy, M., Scarpine, V., Schubnell, M., Serrano, S., Smith, M., Soares-Santos, M., Suchyta, E., Swanson, M., Tarle, G., Thomas, D., To, C., Tucker, D., Varga, T., DES Collaboration, 2021, A&A, 656, A106, Dark Energy Survey Year 3 Results: Galaxy mock catalogs for BAO analysis
    The calibration and validation of scientific analysis in simulations is a fundamental tool to ensure unbiased and robust results in observational cosmology. In particular, mock galaxy catalogs are a crucial resource to achieve these goals in the measurement of baryon acoustic oscillation (BAO) in the clustering of galaxies. Here we present a set of 1952 galaxy mock catalogs designed to mimic the Dark Energy Survey Year 3 BAO sample over its full photometric redshift range 0.6 < zphoto < 1.1. The mocks are based upon 488 ICE-COLA fast N-body simulations of full-sky light cones and were created by populating halos with galaxies, using a hybrid halo occupation distribution - halo abundance matching model. This model has ten free parameters, which were determined, for the first time, using an automatic likelihood minimization procedure. We also introduced a novel technique to assign photometric redshift for simulated galaxies, following a two-dimensional probability distribution with VIMOS Public Extragalactic Redshift Survey data. The calibration was designed to match the observed abundance of galaxies as a function of photometric redshift, the distribution of photometric redshift errors, and the clustering amplitude on scales smaller than those used for BAO measurements. An exhaustive analysis was done to ensure that the mocks reproduce the input properties. Finally, mocks were tested by comparing the angular correlation function w(), angular power spectrum C, and projected clustering p(r) to theoretical predictions and data. The impact of volume replication in the estimate of the covariance is also investigated. The success in accurately reproducing the photometric redshift uncertainties and the galaxy clustering as a function of redshift render this mock creation pipeline as a benchmark for future analyses of photometric galaxy surveys.
  6. Pang, X., Yu, Z., Tang, S., Hong, J., Yuan, Z., Pasquato, M., Kouwenhoven, M., 2021, ApJ, 923, 20, Disruption of Hierarchical Clustering in the Vela OB2 Complex and the Cluster Pair Collinder 135 and UBC 7 with Gaia EDR3: Evidence of Supernova Quenching
    We identify hierarchical structures in the Vela OB2 complex and the cluster pair Collinder 135 and UBC 7 with Gaia EDR3 using the neural network machine-learning algorithm StarGO. Five second-level substructures are disentangled in Vela OB2, which are referred to as Huluwa 1 (Gamma Velorum), Huluwa 2, Huluwa 3, Huluwa 4, and Huluwa 5. For the first time, Collinder 135 and UBC 7 are simultaneously identified as constituent clusters of the pair with minimal manual intervention. We propose an alternative scenario in which Huluwa 1-5 have originated from sequential star formation. The older clusters Huluwa 1-3, with an age of 10-22 Myr, generated stellar feedback to cause turbulence that fostered the formation of the younger-generation Huluwa 4-5 (7-20 Myr). A supernova explosion located inside the Vela IRAS shell quenched star formation in Huluwa 4-5 and rapidly expelled the remaining gas from the clusters. This resulted in global mass stratification across the shell, which is confirmed by the regression discontinuity method. The stellar mass in the lower rim of the shell is 0.32 0.14 M higher than in the upper rim. Local, cluster-scale mass segregation is observed in the lowest-mass cluster Huluwa 5. Huluwa 1-5 (in Vela OB2) are experiencing significant expansion, while the cluster pair suffers from moderate expansion. The velocity dispersions suggest that all five groups (including Huluwa 1A and Huluwa 1B) in Vela OB2 and the cluster pair are supervirial and are undergoing disruption, and also that Huluwa 1A and Huluwa 1B may be a coeval young cluster pair. N-body simulations predict that Huluwa 1-5 in Vela OB2 and the cluster pair will continue to expand in the future 100 Myr and eventually dissolve.
  7. Adhikari, S., Shin, T., Jain, B., Hilton, M., Baxter, E., Chang, C., Wechsler, R., Battaglia, N., Bond, J., Bocquet, S., Choi, S., DeRose, J., Devlin, M., Dunkley, J., Evrard, A., Ferraro, S., Hill, J., Hughes, J., Gallardo, P., Lokken, M., MacInnis, A., Madhavacheril, M., McMahon, J., Nati, F., Newburgh, L., Niemack, M., Page, L., Palmese, A., Partridge, B., Rozo, E., Rykoff, E., Salatino, M., Schillaci, A., Sehgal, N., Sifon, C., To, C., Wollack, E., Wu, H., Xu, Z., Aguena, M., Allam, S., Amon, A., Annis, J., Avila, S., Bacon, D., Bertin, E., Bhargava, S., Brooks, D., Burke, D., Rosell, A., Kind, M., Carretero, J., Castander, F., Choi, A., Costanzi, M., da Costa, L., Vicente, J., Desai, S., Diehl, T., Doel, P., Everett, S., Ferrero, I., Ferte, A., Flaugher, B., Fosalba, P., Frieman, J., Garcia-Bellido, J., Gaztanaga, E., Gruen, D., Gruendl, R., Gschwend, J., Gutierrez, G., Hartley, W., Hinton, S., Hollowood, D., Honscheid, K., James, D., Jeltema, T., Kuehn, K., Kuropatkin, N., Lahav, O., Lima, M., Maia, M., Marshall, J., Martini, P., Melchior, P., Menanteau, F., Miquel, R., Morgan, R., L. C. Ogando, R., Paz-Chinchon, F., Malagon, A., Sanchez, E., Santiago, B., Scarpine, V., Serrano, S., Sevilla-Noarbe, I., Smith, M., Soares-Santos, M., Suchyta, E., E. C. Swanson, M., Varga, T., Wilkinson, R., Zhang, Y., Austermann, J., Beall, J., Becker, D., Denison, E., Duff, S., Hilton, G., Hubmayr, J., Ullom, J., Lanen, J., Vale, L., Vale, L., Vale, L., 2021, ApJ, 923, 37, Probing Galaxy Evolution in Massive Clusters Using ACT and DES: Splashback as a Cosmic Clock
    We measure the projected number density profiles of galaxies and the splashback feature in clusters selected by the Sunyaev-Zel'dovich effect from the Advanced Atacama Cosmology Telescope (AdvACT) survey using galaxies observed by the Dark Energy Survey (DES). The splashback radius is consistent with CDM-only simulations and is located at ${2.4}_{-0.4}^{+0.3}\,\mathrm{Mpc}\,{h}^{-1}$ . We split the galaxies on color and find significant differences in their profile shapes. Red and green-valley galaxies show a splashback-like minimum in their slope profile consistent with theory, while the bluest galaxies show a weak feature at a smaller radius. We develop a mapping of galaxies to subhalos in simulations and assign colors based on infall time onto their hosts. We find that the shift in location of the steepest slope and different profile shapes can be mapped to the average time of infall of galaxies of different colors. The steepest slope traces a discontinuity in the phase space of dark matter halos. By relating spatial profiles to infall time, we can use splashback as a clock to understand galaxy quenching. We find that red galaxies have on average been in clusters over 3.2 Gyr, green galaxies about 2.2 Gyr, while blue galaxies have been accreted most recently and have not reached apocenter. Using the full radial profiles, we fit a simple quenching model and find that the onset of galaxy quenching occurs after a delay of about a gigayear and that galaxies quench rapidly thereafter with an exponential timescale of 0.6 Gyr.
  8. van Belle, G., von Braun, K., Ciardi, D., Pilyavsky, G., Buckingham, R., Boden, A., Clark, C., Hartman, Z., van Belle, G., Bucknew, W., Cole, G., 2021, ApJ, 922, 163, Direct Measurements of Giant Star Effective Temperatures and Linear Radii: Calibration against Spectral Types and V - K Color
    We calculate directly determined values for effective temperature (T eff) and radius (R) for 191 giant stars based upon high-resolution angular size measurements from optical interferometry at the Palomar Testbed Interferometer. Narrow- to wideband photometry data for the giants are used to establish bolometric fluxes and luminosities through spectral energy distribution fitting, which allows for homogeneously establishing an assessment of spectral type and dereddened V 0 - K 0 color; these two parameters are used as calibration indices for establishing trends in T eff and R. Spectral types range from G0III to M7.75III, V 0 - K 0 from 1.9 to 8.5. For the V 0 - K 0 = {1.9, 6.5} range, median T eff uncertainties in the fit of effective temperature versus color are found to be less than 50 K; over this range, T eff drops from 5050 to 3225 K. Linear sizes are found to be largely constant at 11 R from G0III to K0III, increasing linearly with subtype to 50 R at K5III, and then further increasing linearly to 150 R by M8III. Three examples of the utility of this data set are presented: first, a fully empirical Hertzsprung-Russell diagram is constructed and examined against stellar evolution models; second, values for stellar mass are inferred based on measures of R and literature values for $\mathrm{log}g> ; finally, an improved calibration of an angular size prediction tool, based upon V and K values for a star, is presented.
  9. Massey, P., Neugent, K., Dorn-Wallenstein, T., Eldridge, J., Stanway, E., Levesque, E., 2021, ApJ, 922, 177, Testing Evolutionary Models with Red Supergiant and Wolf-Rayet Populations
    Despite the many successes that modern massive star evolutionary theory has enjoyed, reproducing the apparent trend in the relative number of red supergiants (RSGs) and Wolf-Rayet (WR) stars has remained elusive. Previous estimates show the RSG/WR ratio decreasing strongly with increasing metallicity. However, the evolutionary models have always predicted a relatively flat distribution for the RSG/WR ratio. In this paper we reexamine this issue, drawing on recent surveys for RSGs and WRs in the Magellanic Clouds, M31, and M33. The RSG surveys have used Gaia astrometry to eliminate foreground contamination and have separated RSGs from asymptotic giant branch stars using near-infrared colors. The surveys for WRs have utilized interference-filter imaging, photometry, and image subtraction techniques to identify candidates, which have then been confirmed spectroscopically. After carefully matching the observational criteria to the models, we now find good agreement in both the single-star Geneva and binary BPASS models with the new observations. The agreement is better when we shift the RSG effective temperatures derived from J - Ks photometry downwards by 200 K in order to agree with the Levesque TiO effective temperature scale. In an appendix we also present a source list of RSGs for the SMC which includes effective temperatures and luminosities derived from near-infrared 2MASS photometry, in the same manner as used for the other galaxies.
  10. Sharkey, B., Reddy, V., Malhotra, R., Thirouin, A., Kuhn, O., Conrad, A., Rothberg, B., Sanchez, J., Thompson, D., Veillet, C., 2021, ComEE, 2, 231, Lunar-like silicate material forms the Earth quasi-satellite (469219) 2016 HO3 Kamooalewa
    Little is known about Earth quasi-satellites, a class of near-Earth small solar system bodies that orbit the sun but remain close to the Earth, because they are faint and difficult to observe. Here we use the Large Binocular Telescope (LBT) and the Lowell Discovery Telescope (LDT) to conduct a comprehensive physical characterization of quasi-satellite (469219) Kamooalewa and assess its affinity with other groups of near-Earth objects. We find that (469219) Kamooalewa rotates with a period of 28.3 (+1.8/1.3) minutes and displays a reddened reflectance spectrum from 0.4-2.2 microns. This spectrum is indicative of a silicate-based composition, but with reddening beyond what is typically seen amongst asteroids in the inner solar system. We compare the spectrum to those of several material analogs and conclude that the best match is with lunar-like silicates. This interpretation implies extensive space weathering and raises the prospect that Kamo'oalewa could comprise lunar material.
  11. Friedrich, O., Andrade-Oliveira, F., Camacho, H., Alves, O., Rosenfeld, R., Sanchez, J., Fang, X., Eifler, T., Krause, E., Chang, C., Omori, Y., Amon, A., Baxter, E., Elvin-Poole, J., Huterer, D., Porredon, A., Prat, J., Terra, V., Troja, A., Alarcon, A., Bechtol, K., Bernstein, G., Buchs, R., Campos, A., Carnero Rosell, A., Carrasco Kind, M., Cawthon, R., Choi, A., Cordero, J., Crocce, M., Davis, C., DeRose, J., Diehl, H., Dodelson, S., Doux, C., Drlica-Wagner, A., Elsner, F., Everett, S., Fosalba, P., Gatti, M., Giannini, G., Gruen, D., Gruendl, R., Harrison, I., Hartley, W., Jain, B., Jarvis, M., MacCrann, N., McCullough, J., Muir, J., Myles, J., Pandey, S., Raveri, M., Roodman, A., Rodriguez-Monroy, M., Rykoff, E., Samuroff, S., Sanchez, C., Secco, L., Sevilla-Noarbe, I., Sheldon, E., Troxel, M., Weaverdyck, N., Yanny, B., Aguena, M., Avila, S., Bacon, D., Bertin, E., Bhargava, S., Brooks, D., Burke, D., Carretero, J., Costanzi, M., da Costa, L., Pereira, M., De Vicente, J., Desai, S., Evrard, A., Ferrero, I., Frieman, J., Garcia-Bellido, J., Gaztanaga, E., Gerdes, D., Giannantonio, T., Gschwend, J., Gutierrez, G., Hinton, S., Hollowood, D., Honscheid, K., James, D., Kuehn, K., Lahav, O., Lima, M., Maia, M., Menanteau, F., Miquel, R., Morgan, R., Palmese, A., Paz-Chinchon, F., Plazas, A., Sanchez, E., Scarpine, V., Serrano, S., Soares-Santos, M., Smith, M., Suchyta, E., Tarle, G., Thomas, D., To, C., Varga, T., Weller, J., Wilkinson, R., Wilkinson, R., DES Collaboration, 2021, MNRAS, 508, 3125, Dark Energy Survey year 3 results: covariance modelling and its impact on parameter estimation and quality of fit
    We describe and test the fiducial covariance matrix model for the combined two-point function analysis of the Dark Energy Survey Year 3 (DES-Y3) data set. Using a variety of new ansatzes for covariance modelling and testing, we validate the assumptions and approximations of this model. These include the assumption of Gaussian likelihood, the trispectrum contribution to the covariance, the impact of evaluating the model at a wrong set of parameters, the impact of masking and survey geometry, deviations from Poissonian shot noise, galaxy weighting schemes, and other sub-dominant effects. We find that our covariance model is robust and that its approximations have little impact on goodness of fit and parameter estimation. The largest impact on best-fitting figure-of-merit arises from the so-called fsky approximation for dealing with finite survey area, which on average increases the 2 between maximum posterior model and measurement by $3.7{{\ \rm per\ cent}}$ (2 18.9). Standard methods to go beyond this approximation fail for DES-Y3, but we derive an approximate scheme to deal with these features. For parameter estimation, our ignorance of the exact parameters at which to evaluate our covariance model causes the dominant effect. We find that it increases the scatter of maximum posterior values for m and 8 by about $3{{\ \rm per\ cent}}$ and for the dark energy equation-of-state parameter by about $5{{\ \rm per\ cent}}$.
  12. Rodriguez, C., Ferrero, G., Benvenuto, O., Gamen, R., Morrell, N., Barba, R., Arias, J., Massey, P., 2021, MNRAS, 508, 2179, Fundamental parameters of the massive eclipsing binary HM1 8
    We present a comprehensive study of the massive binary system HM1 8, based on multi-epoch high-resolution spectroscopy, V-band photometry, and archival X-ray data. Spectra from the OWN Survey, a high-resolution optical monitoring of Southern O and WN stars, are used to analyse the spectral morphology and perform quantitative spectroscopic analysis of both stellar components. The primary and secondary components are classified as O4.5 IV(f) and O9.7 V, respectively. From a radial velocity (RV) study, we derived a set of orbital parameters for the system. We found an eccentric orbit (e = 0.14 0.01) with a period of P = 5.87820 0.00008 d. Through the simultaneous analysis of the RVs and the V-band light curve, we derived an orbital inclination of 70.0 2.0 and stellar masses of $M_a=33.6^{+1.4}_{-1.2}~\text{M}_{\odot }$ for the primary, and $M_b=17.7^{+0.5}_{-0.7}~\text{M}_{\odot }$ for the secondary. The components show projected rotational velocities vasin i = 105 14 km s-1 and vbsin i = 82 15 km s-1, respectively. A tidal evolution analysis is also performed and found to be in agreement with the orbital characteristics. Finally, the available X-ray observations show no evidence of a colliding winds region; therefore, the X-ray emission is attributed to stellar winds.
  13. Tang, S., Robinson, T., Marley, M., Batalha, N., Lupu, R., Prato, L., 2021, ApJ, 922, 26, Impact of Water-latent Heat on the Thermal Structure of Ultra-cool Objects: Brown Dwarfs and Free-floating Planets
    Brown dwarfs are essential targets for understanding planetary and sub-stellar atmospheres across a wide range of thermal and chemical conditions. As surveys continue to probe ever deeper and as observing capabilities continue to improve, the number of known Y dwarfs-the coldest class of sub-stellar objects, with effective temperatures below about 600 K-is rapidly growing. Critically, this class of ultra-cool objects has atmospheric conditions that overlap with solar-system worlds and, as a result, tools and ideas developed from studying Earth, Jupiter, Saturn, and other nearby worlds are well suited for application to sub-stellar atmospheres. To that end, we developed a one-dimensional (vertical) atmospheric structure model for ultra-cool objects that includes moist adiabatic convection, as this is an important process for many solar-system planets. Application of this model across a range of effective temperatures (350, 300, 250, 200 K), metallicities ([M/H] of 0.0, 0.5, 0.7, 1.5), and gravities (log g of 4.0, 4.5, 4.7, 5.0) demonstrates strong impact of water-latent heat release on simulated temperature-pressure profiles. At the highest metallicities, water-vapor mixing ratios reach an Earth-like 3% with associated major alterations to the thermal structure in the atmospheric regions where water condenses. Spectroscopic and photometric signatures of metallicity and moist convection should be readily detectable at near- and mid-infrared wavelengths, especially with James Webb Space Telescope observations, and can help indicate the formation history of an object.
  14. Hsieh, H., Chandler, C., Denneau, L., Fitzsimmons, A., Erasmus, N., Kelley, M., Knight, M., Lister, T., Pittichova, J., Sheppard, S., Thirouin, A., Trujillo, C., Usher, H., Gomez, E., Chatelain, J., Greenstreet, S., Angel, T., Miles, R., Roche, P., Wooding, B., 2021, ApJL, 922, L9, Physical Characterization of Main-belt Comet (248370) 2005 QN173
    We report results from new and archival observations of the newly discovered active asteroid (248370) 2005 QN173 (also now designated Comet 433P), which has been determined to be a likely main-belt comet based on a subsequent discovery that it is recurrently active near perihelion. From archival data analysis, we estimate $g^{\prime} $ -, $r^{\prime} $ -, $i^{\prime} $ -, and $z^{\prime} $ -band absolute magnitudes for the nucleus of Hg = 16.62 0.13, Hr = 16.12 0.10, Hi = 16.05 0.11, and Hz = 15.93 0.08, corresponding to nucleus colors of $g^{\prime} -r^{\prime} =0.50\pm 0.16$ , $r^{\prime} -i^{\prime} =0.07\pm 0.15$ , and $i^{\prime} -z^{\prime} =0.12\pm 0.14$ ; an equivalent V-band absolute magnitude of HV = 16.32 0.08; and a nucleus radius of rn = 1.6 0.2 km (using a V-band albedo of pV = 0.054 0.012). Meanwhile, we find mean near-nucleus coma colors when 248370 is active of $g^{\prime} -r^{\prime} =0.47\pm 0.03$ , $r^{\prime} -i^{\prime} =0.10\pm 0.04$ , and $i^{\prime} -z^{\prime} =0.05\pm 0.05$ and similar mean dust tail colors, suggesting that no significant gas coma is present. We find approximate ratios between the scattering cross sections of near-nucleus dust (within 5000 km of the nucleus) and the nucleus of Ad/An = 0.7 0.3 on 2016 July 22 and 1.8 < Ad/An < 2.9 in 2021 July and August. During the 2021 observation period, the coma declined in intrinsic brightness by ~0.35 mag (or ~25%) in 37 days, while the surface brightness of the dust tail remained effectively constant over the same period. Constraints derived from the sunward extent of the coma and width of the tail as measured perpendicular to the orbit plane suggest that the terminal velocities of ejected dust grains are extremely slow (~1 m s-1 for 1 m particles), suggesting that the observed dust emission may be aided by rapid rotation of the nucleus lowering the effective escape velocity.
  15. Bernardinelli, P., Bernstein, G., Montet, B., Weryk, R., Wainscoat, R., Aguena, M., Allam, S., Andrade-Oliveira, F., Annis, J., Avila, S., Bertin, E., Brooks, D., Burke, D., Carnero Rosell, A., Carrasco Kind, M., Carretero, J., Cawthon, R., Conselice, C., Costanzi, M., da Costa, L., Pereira, M., De Vicente, J., Diehl, H., Everett, S., Ferrero, I., Flaugher, B., Frieman, J., Garcia-Bellido, J., Gaztanaga, E., Gerdes, D., Gruen, D., Gruendl, R., Gschwend, J., Gutierrez, G., Hinton, S., Hollowood, D., Honscheid, K., James, D., Kuehn, K., Kuropatkin, N., Lahav, O., Maia, M., Marshall, J., Menanteau, F., Miquel, R., Morgan, R., Ogando, R., Paz-Chinchon, F., Pieres, A., Malagon, A., Rodriguez-Monroy, M., Romer, A., Roodman, A., Sanchez, E., Schubnell, M., Serrano, S., Sevilla-Noarbe, I., Smith, M., Soares-Santos, M., Suchyta, E., Swanson, M., Tarle, G., To, C., Troxel, M., Varga, T., Walker, A., Zhang, Y., DES Collaboration, 2021, ApJL, 921, L37, C/2014 UN271 (Bernardinelli-Bernstein): The Nearly Spherical Cow of Comets
    Comet C/2014 UN271 (Bernardinelli-Bernstein), incoming from the Oort cloud, is remarkable in having the brightest (and presumably largest) nucleus of any well-measured comet and having been discovered at the heliocentric distance rh 29 au, farther than any Oort cloud comet. In this work, we describe the discovery process and observations and the properties that can be inferred from images recorded until the first reports of activity in 2021 June. The orbit has i = 95, with a perihelion of 10.97 au to be reached in 2031 and a previous aphelion at 40,400 260 au. Backward integration of the orbit under a standard Galactic tidal model and known stellar encounters suggests a perihelion of q 18 au on its previous perihelion passage 3.5 Myr ago; hence, the current data could be the first ever obtained of a comet that has not been inside Uranus's orbit in 4 Gyr. The photometric data show an unresolved nucleus with absolute magnitude Hr = 8.0, colors that are typical of comet nuclei or Damocloids, and no secular trend as it traversed the range 34-23 au. For the r-band geometric albedo pr, this implies a diameter of $150{({p}_{r}/0.04)}^{-0.5}$ km. There is strong evidence of brightness fluctuations at the 0.2 mag level, but no rotation period can be discerned. A coma, nominally consistent with a "stationary" 1/ surface brightness distribution, grew in scattering cross section at an exponential rate from Af 1 to 150 m as the comet approached from 28 to 20 au. The activity rate is consistent with a very simple model of sublimation of a surface species in radiative equilibrium with the Sun. The inferred enthalpy of sublimation matches those of CO2 and NH3. More volatile species, such as N2, CH4, and CO, must be far less abundant on the sublimating surfaces.
  16. Lopez-Valdivia, R., Sokal, K., Mace, G., Kidder, B., Hussaini, M., Nofi, L., Prato, L., Johns-Krull, C., Oh, H., Lee, J., Park, C., Oh, J., Kraus, A., Kaplan, K., Llama, J., Mann, A., Kim, H., Gully-Santiago, M., Lee, H., Pak, S., Hwang, N., Jaffe, D., 2021, ApJ, 921, 53, The IGRINS YSO Survey. I. Stellar Parameters of Pre-main-sequence Stars in Taurus-Auriga
    We present fundamental parameters for 110 canonical K- and M-type (1.3-0.13 M) Taurus-Auriga young stellar objects (YSOs). The analysis produces a simultaneous determination of effective temperature (Teff), surface gravity (log g), magnetic-field strength (B), and projected rotational velocity ( $v\sin i$ ). Our method employed synthetic spectra and high-resolution (R ~ 45,000) near-infrared spectra taken with the Immersion GRating INfrared Spectrometer (IGRINS) to fit specific K-band spectral regions most sensitive to those parameters. The use of these high-resolution spectra reduces the influence of distance uncertainties, reddening, and non-photospheric continuum emission on the parameter determinations. The median total (fit + systematic) uncertainties were 170 K, 0.28 dex, 0.60 kG, 2.5 km s-1 for Teff, log g, B, and $v\sin i$ , respectively. We determined B for 41 Taurus YSOs (upper limits for the remainder) and find systematic offsets (lower Teff, higher log g and $v\sin i$ ) in parameters when B is measurable but not considered in the fit. The average log g for the Class II and Class III objects differs by 0.23 0.05 dex, which is consistent with Class III objects being the more evolved members of the star-forming region. However, the dispersion in log g is greater than the uncertainties, which highlights how the YSO classification correlates with age (log g), yet there are exceptionally young (lower log g) Class III YSOs and relatively old (higher log g) Class II YSOs with unexplained evolutionary histories. The spectra from this work are provided in an online repository along with TW Hydrae Association comparison objects and the model grid used in our analysis.
  17. Casey, A., Ji, A., Hansen, T., Li, T., Koposov, S., Da Costa, G., Bland-Hawthorn, J., Cullinane, L., Erkal, D., Lewis, G., Kuehn, K., Mackey, D., Martell, S., Pace, A., Simpson, J., Zucker, D., 2021, ApJ, 921, 67, Signature of a Massive Rotating Metal-poor Star Imprinted in the Phoenix Stellar Stream
    The Phoenix stellar stream has a low intrinsic dispersion in velocity and metallicity that implies the progenitor was probably a low-mass globular cluster. In this work we use Magellan/Magellan Inamori Kyocera Echelle (MIKE) high-dispersion spectroscopy of eight Phoenix stream red giants to confirm this scenario. In particular, we find negligible intrinsic scatter in metallicity ( $\sigma ([\mathrm{Fe}\,{\rm\small{II}}/{\rm{H}}])={0.04}_{-0.03}^{+0.11}$ ) and a large peak-to-peak range in [Na/Fe] and [Al/Fe] abundance ratios, consistent with the light element abundance patterns seen in the most metal-poor globular clusters. However, unlike any other globular cluster, we also find an intrinsic spread in [Sr II/Fe] spanning ~1 dex, while [Ba II/Fe] shows nearly no intrinsic spread ( $\sigma ([\mathrm{Ba}\,{\rm\small{II}}/{\rm{H}}])={0.03}_{-0.02}^{+0.10}$ ). This abundance signature is best interpreted as slow-neutron-capture element production from a massive fast-rotating metal-poor star (15-20 M, vini/vcrit = 0.4, [Fe/H] = -3.8). The low inferred cluster mass suggests the system would have been unable to retain supernovae ejecta, implying that any massive fast-rotating metal-poor star that enriched the interstellar medium must have formed and evolved before the globular cluster formed. Neutron-capture element production from asymptotic giant branch stars or magneto-rotational instabilities in core-collapse supernovae provide poor fits to the observations. We also report one Phoenix stream star to be a lithium-rich giant (A(Li) = 3.1 0.1). At [Fe/H ] = -2.93; it is among the most metal-poor lithium-rich giants known. * This paper includes data gathered with the 6.5 m Magellan Telescopes located at Las Campanas Observatory, Chile.
  18. Ji, A., Koposov, S., Li, T., Erkal, D., Pace, A., Simon, J., Belokurov, V., Cullinane, L., Da Costa, G., Kuehn, K., Lewis, G., Mackey, D., Shipp, N., Simpson, J., Zucker, D., Hansen, T., Bland-Hawthorn, J., S5 Collaboration, 2021, ApJ, 921, 32, Kinematics of Antlia 2 and Crater 2 from the Southern Stellar Stream Spectroscopic Survey (S5)
    We present new spectroscopic observations of the diffuse Milky Way satellite galaxies Antlia 2 and Crater 2, taken as part of the Southern Stellar Stream Spectroscopic Survey (S5). The new observations approximately double the number of confirmed member stars in each galaxy and more than double the spatial extent of spectroscopic observations in Antlia 2. A full kinematic analysis, including Gaia EDR3 proper motions, detects a clear velocity gradient in Antlia 2 and a tentative velocity gradient in Crater 2. The velocity gradient magnitudes and directions are consistent with particle stream simulations of tidal disruption. Furthermore, the orbit and kinematics of Antlia 2 require a model that includes the reflex motion of the Milky Way induced by the Large Magellanic Cloud. We also find that Antlia 2's metallicity was previously overestimated, so it lies on the empirical luminosity-metallicity relation and is likely only now experiencing substantial stellar mass loss. Current dynamical models of Antlia 2 require it to have lost over 90% of its stars to tides, in tension with the low stellar mass loss implied by the updated metallicity. Overall, the new kinematic measurements support a tidal disruption scenario for the origin of these large and extended dwarf spheroidal galaxies.
  19. Baines, E., Thomas Armstrong, J., Clark, J., Gorney, J., Hutter, D., Jorgensen, A., Kyte, C., Mozurkewich, D., Nisley, I., Sanborn, J., Schmitt, H., van Belle, G., 2021, AJ, 162, 198, Angular Diameters and Fundamental Parameters of Forty-four Stars from the Navy Precision Optical Interferometer
    We measured the angular diameters of 44 stars with the Navy Precision Optical Interferometer, obtaining uncertainties on the limb-darkened diameter of 2% or less for all but four stars. We then used our diameters with Gaia or Hipparcos parallaxes to calculate each star's physical radius. We gathered information from the literature to determine bolometric flux and luminosity, and combined that with our diameters to produce an effective temperature. Our sample consists of mostly giant stars, and spans a wide range of spectral classes from B to M.
  20. Stern, S., Keeney, B., Hoover, R., Protopapa, S., White, O., Grundy, W., Cruikshank, D., New Horizons Team, 2021, AJ, 162, 207, New Investigations of Dark-floored Pits In the Volatile Ice of Sputnik Planitia on Pluto
    Sputnik Planitia, Pluto's gigantic, volatile ice glacier, hosts numerous scientific mysteries, including the presence of thousands of elongated pit structures. We examine various attributes of these pit structures in New Horizons data sets, revealing their length, aspect ratio, and orientation properties; we also study their reflectivities, colors, and compositions, and compare these attributes to some other relevant regions on Pluto. We then comment on origin mechanisms of the pits and also the fate of the missing volatiles represented by the pits on Sputnik Planitia. From a sample of 317 pits, we find typical length/width ratios of 2-4, with their major axis preferentially oriented approximately north-south. We also find that the floors of large pits in our sample have similar single-scattering albedos and colors to dark material on crater rims and floors (i.e., possible subsurface windows) in Burney basin. We also find that the base of the three pits in our sample, large enough to study with LEISA IR spectroscopy, display both CH4 and N2 absorption features, as do the dark regions in crater windows in Burney basin. Evidence for a sublimation erosion origin for the pits is supported over both the explosion/ejecta venting and structural collapse alternatives. Finally, we find that the mass lost by the pits on Sputnik Planitia most likely lies condensed elsewhere, on Pluto's surface, relocated there by volatile transport as opposed to removal by escape to space or photochemical conversion.
  21. Shin, T., Jain, B., Adhikari, S., Baxter, E., Chang, C., Pandey, S., Salcedo, A., Weinberg, D., Amsellem, A., Battaglia, N., Belyakov, M., Dacunha, T., Goldstein, S., Kravtsov, A., Varga, T., Abbott, T., Aguena, M., Alarcon, A., Allam, S., Amon, A., Andrade-Oliveira, F., Annis, J., Bacon, D., Bechtol, K., Becker, M., Bernstein, G., Bertin, E., Bocquet, S., Bond, J., Brooks, D., Buckley-Geer, E., Burke, D., Campos, A., Rosell, A., Kind, M., Carretero, J., Chen, R., Choi, A., Costanzi, M., da Costa, L., DeRose, J., Desai, S., De Vicente, J., Devlin, M., Diehl, H., Dietrich, J., Dodelson, S., Doel, P., Doux, C., Drlica-Wagner, A., Eckert, K., Elvin-Poole, J., Everett, S., Ferraro, S., Ferrero, I., Ferte, A., Flaugher, B., Frieman, J., Gallardo, P., Gatti, M., Gaztanaga, E., Gerdes, D., Gruen, D., Gruendl, R., Gutierrez, G., Harrison, I., Hartley, W., Hill, J., Hilton, M., Hinton, S., Hollowood, D., Hughes, J., James, D., Jarvis, M., Jeltema, T., Koopman, B., Krause, E., Kuehn, K., Kuropatkin, N., Lahav, O., Lima, M., Lokken, M., MacCrann, N., Madhavacheril, M., Maia, M., McCullough, J., McMahon, J., Melchior, P., Menanteau, F., Miquel, R., Mohr, J., Moodley, K., Morgan, R., Myles, J., Nati, F., Navarro-Alsina, A., Niemack, M., Ogando, R., Page, L., Palmese, A., Partridge, B., Paz-Chinchon, F., Pereira, M., Pieres, A., Malagon, A., Prat, J., Raveri, M., Rodriguez-Monroy, M., Rollins, R., Romer, A., Rykoff, E., Salatino, M., Sanchez, C., Sanchez, E., Santiago, B., Scarpine, V., Schillaci, A., Secco, L., Serrano, S., Sevilla-Noarbe, I., Sheldon, E., Sherwin, B., Sifon, C., Smith, M., Soares-Santos, M., Staggs, S., Suchyta, E., Swanson, M., Tarle, G., Thomas, D., To, C., Troxel, M., Tutusaus, I., Vavagiakis, E., Weller, J., Wollack, E., Yanny, B., Yin, B., Zhang, Y., 2021, MNRAS, 507, 5758, The mass and galaxy distribution around SZ-selected clusters
    We present measurements of the radial profiles of the mass and galaxy number density around Sunyaev-Zel'dovich (SZ)-selected clusters using both weak lensing and galaxy counts. The clusters are selected from the Atacama Cosmology Telescope Data Release 5 and the galaxies from the Dark Energy Survey Year 3 data set. With signal-to-noise ratio of 62 (45) for galaxy (weak lensing) profiles over scales of about 0.2-20 h-1 Mpc, these are the highest precision measurements for SZ-selected clusters to date. Because SZ selection closely approximates mass selection, these measurements enable several tests of theoretical models of the mass and light distribution around clusters. Our main findings are: (1) The splashback feature is detected at a consistent location in both the mass and galaxy profiles and its location is consistent with predictions of cold dark matter N-body simulations. (2) The full mass profile is also consistent with the simulations. (3) The shapes of the galaxy and lensing profiles are remarkably similar for our sample over the entire range of scales, from well inside the cluster halo to the quasilinear regime. We measure the dependence of the profile shapes on the galaxy sample, redshift, and cluster mass. We extend the Diemer & Kravtsov model for the cluster profiles to the linear regime using perturbation theory and show that it provides a good match to the measured profiles. We also compare the measured profiles to predictions of the standard halo model and simulations that include hydrodynamics. Applications of these results to cluster mass estimation, cosmology, and astrophysics are discussed.
  22. Cheng, T., Conselice, C., Aragon-Salamanca, A., Aguena, M., Allam, S., Andrade-Oliveira, F., Annis, J., Bluck, A., Brooks, D., Burke, D., Carrasco Kind, M., Carretero, J., Choi, A., Costanzi, M., da Costa, L., Pereira, M., De Vicente, J., Diehl, H., Drlica-Wagner, A., Eckert, K., Everett, S., Evrard, A., Ferrero, I., Fosalba, P., Frieman, J., Garcia-Bellido, J., Gerdes, D., Giannantonio, T., Gruen, D., Gruendl, R., Gschwend, J., Gutierrez, G., Hinton, S., Hollowood, D., Honscheid, K., James, D., Krause, E., Kuehn, K., Kuropatkin, N., Lahav, O., Maia, M., March, M., Menanteau, F., Miquel, R., Morgan, R., Paz-Chinchon, F., Pieres, A., Plazas Malagon, A., Roodman, A., Sanchez, E., Scarpine, V., Serrano, S., Sevilla-Noarbe, I., Smith, M., Soares-Santos, M., Suchyta, E., Swanson, M., Tarle, G., Thomas, D., To, C., 2021, MNRAS, 507, 4425, Galaxy morphological classification catalogue of the Dark Energy Survey Year 3 data with convolutional neural networks
    We present in this paper one of the largest galaxy morphological classification catalogues to date, including over 20 million galaxies, using the Dark Energy Survey (DES) Year 3 data based on convolutional neural networks (CNNs). Monochromatic i-band DES images with linear, logarithmic, and gradient scales, matched with debiased visual classifications from the Galaxy Zoo 1 (GZ1) catalogue, are used to train our CNN models. With a training set including bright galaxies (16 i < 18) at low redshift (z < 0.25), we furthermore investigate the limit of the accuracy of our predictions applied to galaxies at fainter magnitude and at higher redshifts. Our final catalogue covers magnitudes 16 i < 21, and redshifts z < 1.0, and provides predicted probabilities to two galaxy types - ellipticals and spirals (disc galaxies). Our CNN classifications reveal an accuracy of over 99 per cent for bright galaxies when comparing with the GZ1 classifications (i < 18). For fainter galaxies, the visual classification carried out by three of the co-authors shows that the CNN classifier correctly categorizes discy galaxies with rounder and blurred features, which humans often incorrectly visually classify as ellipticals. As a part of the validation, we carry out one of the largest examinations of non-parametric methods, including ~100 ,000 galaxies with the same coverage of magnitude and redshift as the training set from our catalogue. We find that the Gini coefficient is the best single parameter discriminator between ellipticals and spirals for this data set.
  23. Yu, Z., Martini, P., Penton, A., Davis, T., Malik, U., Lidman, C., Tucker, B., Sharp, R., Kochanek, C., Peterson, B., Aguena, M., Allam, S., Andrade-Oliveira, F., Annis, J., Asorey, J., Bertin, E., Brooks, D., Burke, D., Calcino, J., Carnero Rosell, A., Carollo, D., Carrasco Kind, M., Costanzi, M., da Costa, L., da Silva Pereira, M., Diehl, H., Everett, S., Ferrero, I., Flaugher, B., Frieman, J., Garcia-Bellido, J., Gaztanaga, E., Gerdes, D., Gruen, D., Gruendl, R., Gschwend, J., Gutierrez, G., Hinton, S., Hollowood, D., James, D., Kim, A., Kron, R., Kuehn, K., Kuropatkin, N., Lewis, G., Maia, M., March, M., Marshall, J., Menanteau, F., Miquel, R., Morgan, R., Moller, A., Palmese, A., Paz-Chinchon, F., Plazas, A., Sanchez, E., Scarpine, V., Serrano, S., Sevilla-Noarbe, I., Smith, M., Soares-Santos, M., Suchyta, E., Tarle, G., Thomas, D., To, C., Tucker, D., 2021, MNRAS, 507, 3771, OzDES Reverberation Mapping Programme: the first Mg II lags from 5 yr of monitoring
    Reverberation mapping is a robust method to measure the masses of supermassive black holes outside of the local Universe. Measurements of the radius-luminosity (R-L) relation using the Mg II emission line are critical for determining these masses near the peak of quasar activity at z 1-2, and for calibrating secondary mass estimators based on Mg II that can be applied to large samples with only single-epoch spectroscopy. We present the first nine Mg II lags from our 5-yr Australian Dark Energy Survey reverberation mapping programme, which substantially improves the number and quality of Mg II lag measurements. As the Mg II feature is somewhat blended with iron emission, we model and subtract both the continuum and iron contamination from the multiepoch spectra before analysing the Mg II line. We also develop a new method of quantifying correlated spectroscopic calibration errors based on our numerous, contemporaneous observations of F-stars. The lag measurements for seven of our nine sources are consistent with both the H and Mg II R-L relations reported by previous studies. Our simulations verify the lag reliability of our nine measurements, and we estimate that the median false positive rate of the lag measurements is $4{{\ \rm per\ cent}}$.
  24. Jacoby, G., Hillwig, T., Jones, D., Martin, K., De Marco, O., Kronberger, M., Hurowitz, J., Crocker, A., Dey, J., 2021, MNRAS, 506, 5223, Binary central stars of planetary nebulae identified with Kepler/K2
    We present the identification of 34 likely binary central stars (CSs) of planetary nebulae (PNe) from Kepler/K2 data, seven of which show eclipses. Of these, 29 are new discoveries. Two additional CSs with more complicated variability are also presented. We examined the light curves of all 'possible', 'likely', and 'true' PNe in every Kepler/K2 campaign (0 through 19) to identify CS variability that may indicate a binary CS. For Campaigns 0, 2, 7, 15, and 16, we find 6 likely or confirmed variables among 21 PNe. Our primary effort, though, was focused on Campaign 11 which targeted a Galactic bulge field containing approximately 183 PNe, in which we identified 30 candidate variable CSs. The periods of these variables range from 2.3 h to 30 d, and based on our analysis, most are likely to be close binary star systems. We present periods and preliminary classifications (eclipsing, double degenerate, or irradiated systems) for the likely binaries based on light-curve shape. From our total sample of 204 target PNe, with a correction for incompleteness due to magnitude limits, we calculate a binary fraction of PN central stars to be 20.7 per cent for all the observed PNe, or 23.5 per cent if we limit our sample only to 'true' PNe. However, these fractions are almost certainly lower limits due to the large angular size of the Kepler pixels, which leads to reduced sensitivity in detecting variability, primarily as a result of dilution and noise from the nebula and neighbouring stars. We discuss the binary population of CSs based on these results as part of the total known sample of close binary CSs.
  25. Marciniak, A., Durech, J., Ali-Lagoa, V., Ogoza, W., Szakats, R., Muller, T., Molnar, L., Pal, A., Monteiro, F., Arcoverde, P., Behrend, R., Benkhaldoun, Z., Bernasconi, L., Bosch, J., Brincat, S., Brunetto, L., Butkiewicz-Bak, M., Del Freo, F., Duffard, R., Evangelista-Santana, M., Farroni, G., Fauvaud, S., Fauvaud, M., Ferrais, M., Geier, S., Golonka, J., Grice, J., Hirsch, R., Horbowicz, J., Jehin, E., Julien, P., Kalup, C., Kaminski, K., Kaminska, M., Kankiewicz, P., Kecskemethy, V., Kim, D., Kim, M., Konstanciak, I., Krajewski, J., Kudak, V., Kulczak, P., Kundera, T., Lazzaro, D., Manzini, F., Medeiros, H., Michimani-Garcia, J., Morales, N., Nadolny, J., Oszkiewicz, D., Pakstiene, E., Pawowski, M., Perig, V., Pilcher, F., Pinel, P., Podlewska-Gaca, E., Polakis, T., Richard, F., Rodrigues, T., Rondon, E., Roy, R., Sanabria, J., Santana-Ros, T., Skiff, B., Skrzypek, J., Sobkowiak, K., Sonbas, E., Stachowski, G., Strajnic, J., Trela, P., Tychoniec, .., Urakawa, S., Verebelyi, E., Wagrez, K., Zejmo, M., Zukowski, K., 2021, A&A, 654, A87, Properties of slowly rotating asteroids from the Convex Inversion Thermophysical Model
    Context. Recent results for asteroid rotation periods from the TESS mission showed how strongly previous studies have underestimated the number of slow rotators, revealing the importance of studying those targets. For most slowly rotating asteroids (those with P > 12 h), no spin and shape model is available because of observation selection effects. This hampers determination of their thermal parameters and accurate sizes. Also, it is still unclear whether signatures of different surface material properties can be seen in thermal inertia determined from mid-infrared thermal flux fitting.
    Aims: We continue our campaign in minimising selection effects among main belt asteroids. Our targets are slow rotators with low light-curve amplitudes. Our goal is to provide their scaled spin and shape models together with thermal inertia, albedo, and surface roughness to complete the statistics.
    Methods: Rich multi-apparition datasets of dense light curves are supplemented with data from Kepler and TESS spacecrafts. In addition to data in the visible range, we also use thermal data from infrared space observatories (mainly IRAS, Akari and WISE) in a combined optimisation process using the Convex Inversion Thermophysical Model. This novel method has so far been applied to only a few targets, and therefore in this work we further validate the method itself.
    Results: We present the models of 16 slow rotators, including two updated models. All provide good fits to both thermal and visible data.The obtained sizes are on average accurate at the 5% precision level, with diameters found to be in the range from 25 to 145 km. The rotation periods of our targets range from 11 to 59 h, and the thermal inertia covers a wide range of values, from 2 to <400 J m2 s12 K1, not showing any correlation with the period.
    Conclusions: With this work we increase the sample of slow rotators with reliable spin and shape models and known thermal inertia by 40%. The thermal inertia values of our sample do not display a previously suggested increasing trend with rotation period, which mightbe due to their small skin depth.

    The photometric data with asteroid lightcurves are only available at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (ftp://130.79.128.5) or via http://cdsarc.u-strasbg.fr/viz-bin/cat/J/A+A/654/A87

  26. Lee, Y., Johnstone, D., Lee, J., Herczeg, G., Mairs, S., Contreras-Pena, C., Hatchell, J., Naylor, T., Bell, G., Bourke, T., Broughton, C., Francis, L., Gupta, A., Harsono, D., Liu, S., Park, G., Plovie, S., Moriarty-Schieven, G., Scholz, A., Sharma, T., Teixeira, P., Wang, Y., Aikawa, Y., Bower, G., Vivien Chen, H., Bae, J., Baek, G., Chapman, S., Ping Chen, W., Du, F., Dutta, S., Forbrich, J., Guo, Z., Inutsuka, S., Kang, M., Kirk, H., Kuan, Y., Kwon, W., Lai, S., Lalchand, B., Lane, J., Lee, C., Liu, T., Morata, O., Pearson, S., Pon, A., Sahu, D., Shang, H., Stamatellos, D., Tang, S., Xu, Z., Yoo, H., Rawlings, J., 2021, ApJ, 920, 119, The JCMT Transient Survey: Four-year Summary of Monitoring the Submillimeter Variability of Protostars
    We present the four-year survey results of monthly submillimeter monitoring of eight nearby (<500 pc) star-forming regions by the JCMT Transient Survey. We apply the Lomb-Scargle Periodogram technique to search for and characterize variability on 295 submillimeter peaks brighter than 0.14 Jy beam-1, including 22 disk sources (Class II), 83 protostars (Class 0/I), and 190 starless sources. We uncover 18 secular variables, all of them protostars. No single-epoch burst or drop events and no inherently stochastic sources are observed. We classify the secular variables by their timescales into three groups: Periodic, Curved, and Linear. For the Curved and Periodic cases, the detectable fractional amplitude, with respect to mean peak brightness, is ~4% for sources brighter than ~0.5 Jy beam-1. Limiting our sample to only these bright sources, the observed variable fraction is 37% (16 out of 43). Considering source evolution, we find a similar fraction of bright variables for both Class 0 and Class I. Using an empirically motivated conversion from submillimeter variability to variation in mass accretion rate, six sources (7% of our full sample) are predicted to have years-long accretion events during which the excess mass accreted reaches more than 40% above the total quiescently accreted mass: two previously known eruptive Class I sources, V1647 Ori and EC 53 (V371 Ser), and four Class 0 sources, HOPS 356, HOPS 373, HOPS 383, and West 40. Considering the full protostellar ensemble, the importance of episodic accretion on few years timescale is negligible-only a few percent of the assembled mass. However, given that this accretion is dominated by events on the order of the observing time window, it remains uncertain as to whether the importance of episodic events will continue to rise with decades-long monitoring.
  27. Buie, M., Keeney, B., Strauss, R., Blank, T., Moore, J., Porter, S., Wasserman, L., Weryk, R., Levison, H., Olkin, C., Leiva, R., Bardecker, J., Brown, M., Brown, L., Collins, M., Davidson, H., Dunham, D., Dunham, J., Eaccarino, J., Finley, T., Fuller, L., Garcia, M., George, T., Getrost, K., Gialluca, M., Givot, R., Gupton, D., Hanna, W., Hergenrother, C., Hernandez, Y., Hill, B., Hinton, P., Holt, T., Howell, R., Jewell, J., Kamin, R., Kammer, J., Kareta, T., Kayl, G., Keller, J., Kenyon, D., Kester, S., Kidd, J., Lauer, T., Leung, C., Lorusso, Z., Lundgren, C., Magana, L., Maley, P., Marchis, F., Marcialis, R., McCandless, A., McCrystal, D., McGraw, A., Miller, K., Mueller, B., Noonan, J., Olsen, A., Patton, A., Peluso, D., Person, M., Rigby, J., Rolfsmeier, A., Salmon, J., Samaniego, J., Sawyer, R., Schulz, D., Skrutskie, M., Smith, R., Spencer, J., Springmann, A., Stanbridge, D., Stoffel, T., Tamblyn, P., Tobias, B., Verbiscer, A., von Schalscha, M., Werts, H., Zhang, Q., 2021, PSJ, 2, 202, Size and Shape of (11351) Leucus from Five Occultations
    We present observations of five stellar occultations for (11351) Leucus and reports from two efforts on (21900) Orus. Both objects are prime mission candidate targets for the Lucy Discovery mission. Combined results for Leucus indicate a very dark surface with pV = 0.037 0.001, which is derived from the average of the multichord occultations. Our estimate of the triaxial ellipsoidal shape is for axial diameters of 63.8 36.6 29.6 km assuming that the spin pole is normal to the line of sight. The actual shape of the object is only roughly elliptical in profile at each epoch. Significant topography is seen with horizontal scales up to 30 km and vertical scales up to 5 km. The most significant feature is a large depression on the southern end of the object as seen from a terrestrial viewpoint. For this work we developed a method to correct for differential refraction, accounting for the difference in color between the target object and the reference stars for astrometry derived from ground-based images.
  28. Brown, M., Levison, H., Noll, K., Binzel, R., Buie, M., Grundy, W., Marchi, S., Olkin, C., Spencer, J., Statler, T., Weaver, H., 2021, PSJ, 2, 170, The Orbit and Density of the Jupiter Trojan Satellite System Eurybates-Queta
    We report observations of the Jupiter Trojan asteroid (3548) Eurybates and its satellite Queta with the Hubble Space Telescope and use these observations to perform an orbital fit to the system. Queta orbits Eurybates with a semimajor axis of 2350 11 km at a period of 82.46 0.06 days and an eccentricity of 0.125 0.009. From this orbit we derive a mass of Eurybates of 1.51 0.03 1017 kg, corresponding to an estimated density of 1.1 0.3 g cm-3, broadly consistent with densities measured for other Trojans, C-type asteroids in the outer main asteroid belt, and small icy objects from the Kuiper Belt. Eurybates is the parent body of the only major collisional family among the Jupiter Trojans; its low density suggests that it is a typical member of the Trojan population. Detailed study of this system in 2027 with the Lucy spacecraft flyby should allow significant insight into collisional processes among what appear to be the icy bodies of the Trojan belt.
  29. Levison, H., Olkin, C., Noll, K., Marchi, S., Bell, J., Bierhaus, E., Binzel, R., Bottke, W., Britt, D., Brown, M., Buie, M., Christensen, P., Emery, J., Grundy, W., Hamilton, V., Howett, C., Mottola, S., Patzold, M., Reuter, D., Spencer, J., Statler, T., Stern, S., Sunshine, J., Weaver, H., Wong, I., 2021, PSJ, 2, 171, Lucy Mission to the Trojan Asteroids: Science Goals
    The Lucy Mission is a NASA Discovery-class mission to send a highly capable and robust spacecraft to investigate seven primitive bodies near both the L4 and L5 Lagrange points with Jupiter: the Jupiter Trojan asteroids. These planetesimals from the outer planetary system have been preserved since early in solar system history. The Lucy mission will fly by and extensively study a diverse selection of Trojan asteroids, including all the recognized taxonomic classes, a collisional family member, and a near equal-mass binary. It will visit objects with diameters ranging from roughly 1 km to 100 km. The payload suite consists of a color camera and infrared imaging spectrometer, a high-resolution panchromatic imager, and a thermal infrared spectrometer. Additionally, two spacecraft subsystems will also contribute to the science investigations: the terminal tracking cameras will supplement imaging during closest approach and the telecommunication subsystem will be used to measure the mass of the Trojans. The science goals are derived from the 2013 Planetary Decadal Survey and include determining the surface composition, assessing the geology, determining the bulk properties, and searching for satellites and rings.
  30. Olkin, C., Levison, H., Vincent, M., Noll, K., Andrews, J., Gray, S., Good, P., Marchi, S., Christensen, P., Reuter, D., Weaver, H., Patzold, M., Bell, J., Hamilton, V., Dello Russo, N., Simon, A., Beasley, M., Grundy, W., Howett, C., Spencer, J., Ravine, M., Caplinger, M., 2021, PSJ, 2, 172, Lucy Mission to the Trojan Asteroids: Instrumentation and Encounter Concept of Operations
    The Lucy Mission accomplishes its science during a series of five flyby encounters with seven Trojan asteroid targets. This mission architecture drives a concept of operations design that maximizes science return, provides redundancy in observations where possible, features autonomous fault protection, and utilizes onboard target tracking near closest approach. These design considerations reduce risk during the relatively short time-critical periods when science data is collected. The payload suite consists of a color camera and infrared imaging spectrometer, a high-resolution panchromatic imager, and a thermal infrared spectrometer. The mission design allows for concurrent observations of all instruments. Additionally, two spacecraft subsystems will also contribute to the science investigations: the Terminal Tracking Cameras will obtain wide field-of-view imaging near closest approach to determine the shape of each of the Trojan targets and the telecommunication subsystem will carry out Doppler tracking of the spacecraft to determine the mass of each of the Trojan targets.
  31. Rivkin, A., Chabot, N., Stickle, A., Thomas, C., Richardson, D., Barnouin, O., Fahnestock, E., Ernst, C., Cheng, A., Chesley, S., Naidu, S., Statler, T., Barbee, B., Agrusa, H., Moskovitz, N., Terik Daly, R., Pravec, P., Scheirich, P., Dotto, E., Della Corte, V., Michel, P., Kuppers, M., Atchison, J., Hirabayashi, M., 2021, PSJ, 2, 173, The Double Asteroid Redirection Test (DART): Planetary Defense Investigations and Requirements
    The Double Asteroid Redirection Test (DART) is a Planetary Defense mission, designed to demonstrate the kinetic impactor technique on (65803) Didymos I Dimorphos, the secondary of the (65803) Didymos system. DART has four level 1 requirements to meet in order to declare mission success: (1) impact Dimorphos between 2022 September 25 and October 2, (2) cause at least a 73 s change in its binary orbit period via the impact, (3) measure the change in binary period to an uncertainty of 7.3 s or less, and (4) measure the momentum transfer efficiency () of the impact and characterize the resulting effects of the impact. The data necessary to achieve these requirements will be obtained and analyzed by the DART Investigation Team. We discuss the rationales for the data to be gathered, the analyses to be undertaken, and how mission success will be achieved.
  32. Domingue, D., Kitazato, K., Matsuoka, M., Yokota, Y., Tatsumi, E., Iwata, T., Abe, M., Ohtake, M., Matsuura, S., Schroder, S., Vilas, F., Barucci, A., Brunetto, R., Takir, D., Le Corre, L., Moskovitz, N., 2021, PSJ, 2, 178, Spectrophotometric Properties of 162173 Ryugu's Surface from the NIRS3 Opposition Observations
    Examination of the opposition geometry properties show that Ryugu's surface regolith is commensurate with laboratory studies of the photometric behavior of powdered carbonaceous chondrites. The regolith is consistent with a broad grain size distribution that contains a fine-grained component.
  33. Fayolle, M., Quirico, E., Schmitt, B., Jovanovic, L., Gautier, T., Carrasco, N., Grundy, W., Vuitton, V., Poch, O., Protopapa, S., Young, L., Cruikshank, D., Dalle Ore, C., Bertrand, T., Stern, A., New Horizons Surface Composition Science Theme Team, 2021, Icar, 367, 114574, Testing tholins as analogues of the dark reddish material covering Pluto's Cthulhu region
    Pluto's fly-by by the New Horizons spacecraft in July 2015 has revealed a dark reddish equatorial region, named Cthulhu, covered by a dark, non-icy material whose origin and composition have yet to be determined. It has been suggested that this material could form from the sedimentation of photochemical aerosols, originating from dissociation and ionisation processes in Pluto's high atmosphere (similarly to aerosols forming Titan's haze). This hypothesis is here further investigated by comparing New Horizons spectra collected both in the visible and the near-infrared to laboratory reflectance measurements of analogues of Pluto's aerosols (Pluto tholins). These aerosols were synthesised in conditions mimicking Pluto's atmosphere, and their optical and reflectance properties were determined, before being used in Hapke models. In particular, the single scattering albedo and phase function of Pluto tholins were retrieved through Hapke model inversion, performed from laboratory reflectance spectra collected under various geometries. From reconstructed reflectance spectra and direct comparison with New Horizons data, some of these tholins are shown to reproduce the photometric level (i.e. reflectance continuum) reasonably well in the near-infrared. Nevertheless, a misfit of the red visible slope still remains and tholins absorption bands present in the modelled spectra are absent in those collected by the New Horizons instruments. Several hypotheses are considered to explain the absence of these absorption features in LEISA data, namely high porosity effects or GCR irradiation. The formation of highly porous structures, which is currently our preferred scenario, could be promoted by either sublimation of ices initially mixed with the aerosols, or gentle deposition under Pluto's weak gravity.
  34. Erasmus, N., Kramer, D., McNeill, A., Trilling, D., Janse van Rensburg, P., van Belle, G., Tonry, J., Denneau, L., Heinze, A., Weiland, H., 2021, MNRAS, 506, 3872, Discovery of superslow rotating asteroids with ATLAS and ZTF photometry
    We present here the discovery of a new class of superslow rotating asteroids (Prot 1000 h) in data extracted from the Asteroid Terrestrial-impact Last Alert System (ATLAS) and Zwicky Transient Facility (ZTF) all-sky surveys. Of the 39 rotation periods we report here, 32 have periods longer than any previously reported unambiguous rotation periods currently in the Asteroid Light Curve Data base. In our sample, seven objects have a rotation period >4000 h and the longest period we report here is 4812 h (~200 d). We do not observe any correlation between taxonomy, albedo, or orbital properties with superslow rotating status. The most plausible mechanism for the creation of these very slow rotators is if their rotations were slowed by YORP spin-down. Superslow rotating asteroids may be common, with at least 0.4 per cent of the main-belt asteroid population with a size range between 2 and 20 km in diameter rotating with periods longer than 1000 h.
  35. Gregg, M., West, M., Lemaux, B., Kupper, A., 2021, MNRAS, 506, 2260, Discovery of tidal debris stars from G1/Mayall II in M31
    The object Mayall II or G1 is the brightest globular cluster belonging to M31. Because of its extreme properties for a globular cluster, it has been speculated that G1 is the remnant nucleus of a dwarf galaxy that has been stripped by the tidal field of M31. Using the Keck DEIMOS spectrograph, we have conducted a survey for tidally stripped stars from G1, obtaining a sample of 351 stellar velocities over ~320 sq. arcmin of sky centred on G1. 13 are within $25~{\, \rm km\, s^{-1}\, }$ of the systemic velocity of G1, and exhibit spatial and velocity correlations consistent with being dynamically associated with G1, and all 13 are well outside the tidal radius of the cluster. These 13 stars could be either (i) the remnants of an almost completely evaporated stellar envelope or (ii) G1 member stars lost through tidal interaction with M31. Estimates of the implied mass-loss rate based on our data suggest a short dissolution time-scale for G1, thus favouring the stellar envelope hypothesis for the origin of the tidal tail stars, or, at the very least, an advanced stage of cluster dissolution. In either case, G1 and by extension compact stellar systems in general have likely played a significant role in building the halo of M31.
  36. Vega-Ferrero, J., Dominguez Sanchez, H., Bernardi, M., Huertas-Company, M., Morgan, R., Margalef, B., Aguena, M., Allam, S., Annis, J., Avila, S., Bacon, D., Bertin, E., Brooks, D., Carnero Rosell, A., Carrasco Kind, M., Carretero, J., Choi, A., Conselice, C., Costanzi, M., da Costa, L., Pereira, M., De Vicente, J., Desai, S., Ferrero, I., Fosalba, P., Frieman, J., Garcia-Bellido, J., Gruen, D., Gruendl, R., Gschwend, J., Gutierrez, G., Hartley, W., Hinton, S., Hollowood, D., Honscheid, K., Hoyle, B., Jarvis, M., Kim, A., Kuehn, K., Kuropatkin, N., Lima, M., Maia, M., Menanteau, F., Miquel, R., Ogando, R., Palmese, A., Paz-Chinchon, F., Plazas, A., Romer, A., Sanchez, E., Scarpine, V., Schubnell, M., Serrano, S., Sevilla-Noarbe, I., Smith, M., Suchyta, E., Swanson, M., Tarle, G., Tarsitano, F., To, C., Tucker, D., Varga, T., Wilkinson, R., 2021, MNRAS, 506, 1927, Pushing automated morphological classifications to their limits with the Dark Energy Survey
    We present morphological classifications of ~27 million galaxies from the Dark Energy Survey (DES) Data Release 1 (DR1) using a supervised deep learning algorithm. The classification scheme separates: (a) early-type galaxies (ETGs) from late-type galaxies (LTGs); and (b) face-on galaxies from edge-on. Our convolutional neural networks (CNNs) are trained on a small subset of DES objects with previously known classifications. These typically have mr 17.7 mag; we model fainter objects to mr < 21.5 mag by simulating what the brighter objects with well-determined classifications would look like if they were at higher redshifts. The CNNs reach 97 per cent accuracy to mr < 21.5 on their training sets, suggesting that they are able to recover features more accurately than the human eye. We then used the trained CNNs to classify the vast majority of the other DES images. The final catalogue comprises five independent CNN predictions for each classification scheme, helping to determine if the CNN predictions are robust or not. We obtain secure classifications for ~87 per cent and 73 per cent of the catalogue for the ETG versus LTG and edge-on versus face-on models, respectively. Combining the two classifications (a) and (b) helps to increase the purity of the ETG sample and to identify edge-on lenticular galaxies (as ETGs with high ellipticity). Where a comparison is possible, our classifications correlate very well with Sersic index (n), ellipticity (), and spectral type, even for the fainter galaxies. This is the largest multiband catalogue of automated galaxy morphologies to date.
  37. Jenniskens, P., Lauretta, D., Towner, M., Heathcote, S., Jehin, E., Hanke, T., Cooper, T., Baggaley, J., Howell, J., Johannink, C., Breukers, M., Odeh, M., Moskovitz, N., Juneau, L., Beck, T., De Cicco, M., Samuels, D., Rau, S., Albers, J., Gural, P., 2021, Icar, 365, 114469, Meteor showers from known long-period comets
    What long-period comets with orbital periods >250 years cause detectable meteor showers on Earth? Low-light video cameras are used to track the motion of +4 to -5 magnitude meteors in our atmosphere by triangulation and calculate the meteoroid orbit in space. In recent years, the CAMS (Cameras for Allsky Meteor Surveillance) low-light video camera network was greatly expanded and, together with other video networks, now has increased the total video meteoroid orbit database to over 2.2 million orbits. Here, we searched this database for meteor showers associated with known long-period comets. Previously, five associations were known. Now, we find 14, as well as six uncertain but likely associations. These showers show a change of longitude of perihelion with node that is a strong function of inclination. Showers of longer duration show a steeper magnitude distribution index, presumably due to aging of the meteoroid population. Showers are generally detected only if the orbital period of the comet is less than 4000 years and the Earth-Comet orbital miss distance is 0.10 AU. The lack of an associated meteor shower sets lower limits on the orbital period of poorly observed comets.
  38. Fausnaugh, M., Morgan, E., Vanderspek, R., Pepper, J., Burke, C., Levine, A., Rudat, A., Villasenor, J., Vezie, M., Goeke, R., Ricker, G., Latham, D., Seager, S., Winn, J., Jenkins, J., Bakos, G., Barclay, T., Berta-Thompson, Z., Bouma, L., Boyd, P., Brasseur, C., Burt, J., Caldwell, D., Charbonneau, D., Christensen-Dalsgaard, J., Clampin, M., Collins, K., Colon, K., De Lee, N., Dunham, E., Fleming, S., Fong, W., Garcia Soto, A., Scott Gaudi, B., Guerrero, N., Hesse, K., Holman, M., Huang, C., Kaltenegger, L., Lissauer, J., McDermott, S., McLean, B., Mireles, I., Mullally, S., Oelkers, R., Paegert, M., Pal, A., Quintana, E., Rinehart, S., Rodriguez, D., Rose, M., Sasselov, D., Schlieder, J., Sha, L., Shporer, A., Smith, J., Stassun, K., Tenenbaum, P., Ting, E., Torres, G., Twicken, J., Vanderburg, A., Wohler, B., Yu, L., 2021, PASP, 133, 095002, The TESS Mission Target Selection Procedure
    We describe the target selection procedure by which stars are selected for 2 minute and 20 s observations by TESS. We first list the technical requirements of the TESS instrument and ground systems processing that limit the total number of target slots. We then describe algorithms used by the TESS Payload Operation Center (POC) to merge candidate targets requested by the various TESS mission elements (the Target Selection Working Group, TESS Asteroseismic Science Consortium, and Guest Investigator office). Lastly, we summarize the properties of the observed TESS targets over the two-year primary TESS mission. We find that the POC target selection algorithm results in 2.1-3.4 times as many observed targets as target slots allocated for each mission element. We also find that the sky distribution of observed targets is different from the sky distributions of candidate targets due to technical constraints that require a relatively even distribution of targets across the TESS fields of view. We caution researchers exploring statistical analyses of TESS planet-host stars that the population of observed targets cannot be characterized by any simple set of criteria applied to the properties of the input Candidate Target Lists.
  39. Ghaffari, Z., Haas, M., Chiaberge, M., Willner, S., Chini, R., Hildebrandt, H., de Propris, R., West, M., 2021, A&A, 653, A44, Clustering of red and blue galaxies around high-redshift 3C radio sources as seen by the Hubble Space Telescope
    To properly understand the evolution of high-redshift galaxy clusters, both passive and star-forming galaxies have to be considered. Here we study the clustering environment of 21 radio galaxies and quasars at 1 < z < 2.5 from the third Cambridge catalog (3C). We use optical and near-infrared Hubble Space Telescope images with a 2' field-of-view, where the filters encompass the rest-frame 4000 A break. Passive red and star-forming blue galaxies were separated in the color-magnitude diagram using a redshift-dependent cut derived from galaxy evolution models. We find that about 16 of 21 radio sources inhabit a galaxy overdensity on scales of 250 kpc (30) projected radius. The sample shows a diversity of red and blue overdensities and also sometimes a deficiency of blue galaxies in the center. The following tentative evolutionary trends are seen: extended proto-clusters with only weak overdensities at z > 1.6, red overdensities at 1.2 < z < 1.6, and red overdensities with an increased deficit of central blue galaxies at z < 1.2. Only a few 3C sources show a blue overdensity tracing active star-formation in the cluster centers; this rarity could indicate that the powerful quasar activity may quench star-formation in the vicinity of most radio sources. The derived number of central luminous red galaxies and the radial density profiles are comparable to those found in local clusters, indicating that some 3C clusters are already mass-rich and compact.
  40. Ellis, T., Boyajian, T., von Braun, K., Ligi, R., Mourard, D., Dragomir, D., Schaefer, G., Farrington, C., 2021, AJ, 162, 118, Directly Determined Properties of HD 97658 from Interferometric Observations
    We conducted interferometric observations with the CHARA Array of transiting super-Earth host HD 97658 and measured its limb-darkened angular diameter to be LD = 0.314 0.004 mas. The combination of the angular diameter with the Gaia EDR3 parallax value with zero-point correction ( = 46.412 0.022 mas, d =21.546 0.011 pc) yields a physical radius of R = 0.728 0.008 R. We also measured the bolometric flux of the star to be Fbol = 2.42 0.05 10-8 erg s-1 cm-2, which, together with angular size, allows a measurement of the effective temperature Teff = 5212 43 K. Our directly determined physical stellar properties are in good agreement with previous estimates derived from spectroscopy. We used our measurements in combination with stellar evolutionary models and properties of the transit of HD 97658 b to determine the mass and age of HD 97658 as well as constrain the properties of the planet. Our results and our analysis of the TESS light curve on the planet (TOI-1821) corroborate previous studies of this system with tighter uncertainties.
  41. Fortino, W., Bernstein, G., Bernardinelli, P., Aguena, M., Allam, S., Annis, J., Bacon, D., Bechtol, K., Bhargava, S., Brooks, D., Burke, D., Carretero, J., Choi, A., Costanzi, M., Costa, L., Pereira, M., De Vicente, J., Desai, S., Doel, P., Drlica-Wagner, A., Eckert, K., Eifler, T., Evrard, A., Ferrero, I., Frieman, J., Garcia-Bellido, J., Gaztanaga, E., Gerdes, D., Gruendl, R., Gschwend, J., Gutierrez, G., Hartley, W., Hinton, S., Hollowood, D., Honscheid, K., James, D., Jarvis, M., Kent, S., Kuehn, K., Kuropatkin, N., Maia, M., Marshall, J., Menanteau, F., Miquel, R., Morgan, R., Myles, J., Ogando, R., Palmese, A., Paz-Chinchon, F., Plazas, A., Roodman, A., Rykoff, E., Sanchez, E., Santiago, B., Scarpine, V., Schubnell, M., Serrano, S., Sevilla-Noarbe, I., Smith, M., Suchyta, E., Tarle, G., To, C., Tucker, D., Varga, T., Walker, A., Weller, J., Wester, W., DES Collaboration, 2021, AJ, 162, 106, Reducing Ground-based Astrometric Errors with Gaia and Gaussian Processes
    Stochastic field distortions caused by atmospheric turbulence are a fundamental limitation to the astrometric accuracy of ground-based imaging. This distortion field is measurable at the locations of stars with accurate positions provided by the Gaia DR2 catalog; we develop the use of Gaussian process regression (GPR) to interpolate the distortion field to arbitrary locations in each exposure. We introduce an extension to standard GPR techniques that exploits the knowledge that the 2D distortion field is curl-free. Applied to several hundred 90 s exposures from the Dark Energy Survey as a test bed, we find that the GPR correction reduces the variance of the turbulent astrometric distortions 12 , on average, with better performance in denser regions of the Gaia catalog. The rms per-coordinate distortion in the riz bands is typically 7 mas before any correction and 2 mas after application of the GPR model. The GPR astrometric corrections are validated by the observation that their use reduces, from 10 to 5 mas rms, the residuals to an orbit fit to riz-band observations over 5 yr of the r = 18.5 trans-Neptunian object Eris. We also propose a GPR method, not yet implemented, for simultaneously estimating the turbulence fields and the 5D stellar solutions in a stack of overlapping exposures, which should yield further turbulence reductions in future deep surveys.
  42. Cigan, P., Young, L., Gomez, H., Madden, S., De Vis, P., Hunter, D., Elmegreen, B., Brinks, E., LITTLE THINGS Team, 2021, AJ, 162, 83, Herschel Photometric Observations of LITTLE THINGS Dwarf Galaxies
    We present here far-infrared photometry of galaxies in a sample that is relatively unexplored at these wavelengths: low-metallicity dwarf galaxies with moderate star formation rates. Four dwarf irregular galaxies from the LITTLE THINGS survey are considered, with deep Herschel PACS and SPIRE observations at 100, 160, 250, 350, and 500 m. Results from modified blackbody fits indicate that these galaxies have low dust masses and cooler dust temperatures than more actively star-forming dwarfs, occupying the lowest LTIR and Mdust regimes seen among these samples. Dust-to-gas mass ratios of ~10-5 are lower, overall, than in more massive and active galaxies but are roughly consistent with the broken power-law relation between the dust-to-gas ratio and metallicity found for other low-metallicity systems. Chemical evolution modeling suggests that these dwarf galaxies are likely forming very little dust via stars or grain growth and have very high dust destruction rates.
  43. Wiseman, P., Sullivan, M., Smith, M., Frohmaier, C., Vincenzi, M., Graur, O., Popovic, B., Armstrong, P., Brout, D., Davis, T., Galbany, L., Hinton, S., Kelsey, L., Kessler, R., Lidman, C., Moller, A., Nichol, R., Rose, B., Scolnic, D., Toy, M., Zontou, Z., Asorey, J., Carollo, D., Glazebrook, K., Lewis, G., Tucker, B., Abbott, T., Aguena, M., Allam, S., Andrade-Oliveira, F., Annis, J., Bacon, D., Bertin, E., Brooks, D., Buckley-Geer, E., Burke, D., Carnero Rosell, A., Carrasco Kind, M., Carretero, J., Costanzi, M., da Costa, L., Pereira, M., Desai, S., Diehl, H., Doel, P., Everett, S., Ferrero, I., Flaugher, B., Fosalba, P., Frieman, J., Garcia-Bellido, J., Gaztanaga, E., Giannantonio, T., Gruen, D., Gruendl, R., Gschwend, J., Gutierrez, G., Hollowood, D., Honscheid, K., Hoyle, B., James, D., Krause, E., Kuehn, K., Kuropatkin, N., Maia, M., Marshall, J., Martini, P., Menanteau, F., Miquel, R., Morgan, R., Ogando, R., Palmese, A., Paz-Chinchon, F., Petravick, D., Pieres, A., Plazas Malagon, A., Romer, A., Sanchez, E., Scarpine, V., Schubnell, M., Serrano, S., Sevilla-Noarbe, I., Soares-Santos, M., Suchyta, E., Swanson, M., Tarle, G., Thomas, D., To, C., Varga, T., Walker, A., DES Collaboration, 2021, MNRAS, 506, 3330, Rates and delay times of Type Ia supernovae in the Dark Energy Survey
    We use a sample of 809 photometrically classified Type Ia supernovae (SNe Ia) discovered by the Dark Energy Survey (DES) along with 40 415 field galaxies to calculate the rate of SNe Ia per galaxy in the redshift range 0.2 < z < 0.6. We recover the known correlation between SN Ia rate and galaxy stellar mass across a broad range of scales 8.5 log (M*/M) 11.25. We find that the SN Ia rate increases with stellar mass as a power law with index 0.63 0.02, which is consistent with the previous work. We use an empirical model of stellar mass assembly to estimate the average star formation histories (SFHs) of galaxies across the stellar mass range of our measurement. Combining the modelled SFHs with the SN Ia rates to estimate constraints on the SN Ia delay time distribution (DTD), we find that the data are fit well by a power-law DTD with slope index = -1.13 0.05 and normalization A = 2.11 0.05 10-13 SNe M-1 yr-1, which corresponds to an overall SN Ia production efficiency $N_{\mathrm{Ia}}/M_* = 0.9~_{-0.7}^{+4.0} \times 10^{-3}~\mathrm{SNe}~\mathrm{M}_{\odot }^{-1}$. Upon splitting the SN sample by properties of the light curves, we find a strong dependence on DTD slope with the SN decline rate, with slower-declining SNe exhibiting a steeper DTD slope. We interpret this as a result of a relationship between intrinsic luminosity and progenitor age, and explore the implications of the result in the context of SN Ia progenitors.
  44. Myles, J., Alarcon, A., Amon, A., Sanchez, C., Everett, S., DeRose, J., McCullough, J., Gruen, D., Bernstein, G., Troxel, M., Dodelson, S., Campos, A., MacCrann, N., Yin, B., Raveri, M., Amara, A., Becker, M., Choi, A., Cordero, J., Eckert, K., Gatti, M., Giannini, G., Gschwend, J., Gruendl, R., Harrison, I., Hartley, W., Huff, E., Kuropatkin, N., Lin, H., Masters, D., Miquel, R., Prat, J., Roodman, A., Rykoff, E., Sevilla-Noarbe, I., Sheldon, E., Wechsler, R., Yanny, B., Abbott, T., Aguena, M., Allam, S., Annis, J., Bacon, D., Bertin, E., Bhargava, S., Bridle, S., Brooks, D., Burke, D., Carnero Rosell, A., Carrasco Kind, M., Carretero, J., Castander, F., Conselice, C., Costanzi, M., Crocce, M., da Costa, L., Pereira, M., Desai, S., Diehl, H., Eifler, T., Elvin-Poole, J., Evrard, A., Ferrero, I., Ferte, A., Flaugher, B., Fosalba, P., Frieman, J., Garcia-Bellido, J., Gaztanaga, E., Giannantonio, T., Hinton, S., Hollowood, D., Honscheid, K., Hoyle, B., Huterer, D., James, D., Krause, E., Kuehn, K., Lahav, O., Lima, M., Maia, M., Marshall, J., Martini, P., Melchior, P., Menanteau, F., Mohr, J., Morgan, R., Muir, J., Ogando, R., Palmese, A., Paz-Chinchon, F., Plazas, A., Rodriguez-Monroy, M., Samuroff, S., Sanchez, E., Scarpine, V., Secco, L., Serrano, S., Smith, M., Soares-Santos, M., Suchyta, E., Swanson, M., Tarle, G., Thomas, D., To, C., Varga, T., Weller, J., Wester, W., 2021, MNRAS, 505, 4249, Dark Energy Survey Year 3 results: redshift calibration of the weak lensing source galaxies
    Determining the distribution of redshifts of galaxies observed by wide-field photometric experiments like the Dark Energy Survey (DES) is an essential component to mapping the matter density field with gravitational lensing. In this work we describe the methods used to assign individual weak lensing source galaxies from the DES Year 3 Weak Lensing Source Catalogue to four tomographic bins and to estimate the redshift distributions in these bins. As the first application of these methods to data, we validate that the assumptions made apply to the DES Y3 weak lensing source galaxies and develop a full treatment of systematic uncertainties. Our method consists of combining information from three independent likelihood functions: self-organizing map p(z) (SOMPZ), a method for constraining redshifts from galaxy photometry; clustering redshifts (WZ), constraints on redshifts from cross-correlations of galaxy density functions; and shear ratios (SRs), which provide constraints on redshifts from the ratios of the galaxy-shear correlation functions at small scales. Finally, we describe how these independent probes are combined to yield an ensemble of redshift distributions encapsulating our full uncertainty. We calibrate redshifts with combined effective uncertainties of z ~ 0.01 on the mean redshift in each tomographic bin.
  45. Eisner, N., Nicholson, B., Barragan, O., Aigrain, S., Lintott, C., Kaye, L., Klein, B., Miller, G., Taylor, J., Zicher, N., Buchhave, L., Caldwell, D., Horner, J., Llama, J., Mortier, A., Rajpaul, V., Stassun, K., Sporer, A., Tkachenko, A., Jenkins, J., Latham, D., Ricker, G., Seager, S., Winn, J., Alhassan, S., Baeten, E., Bean, S., Bundy, D., Efremov, V., Ferstenou, R., Goodwin, B., Hof, M., Hoffman, T., Hubert, A., Lau, L., Lee, S., Maetschke, D., Peltsch, K., Rubio-Alfaro, C., Wilson, G., 2021, MNRAS, 505, 1827, Planet Hunters TESS III: two transiting planets around the bright G dwarf HD 152843
    We report on the discovery and validation of a two-planet system around a bright (V = 8.85 mag) early G dwarf (1.43 R, 1.15 M, TOI 2319) using data from NASA's Transiting Exoplanet Survey Satellite (TESS). Three transit events from two planets were detected by citizen scientists in the month-long TESS light curve (sector 25), as part of the Planet Hunters TESS project. Modelling of the transits yields an orbital period of $11.6264 _{ - 0.0025 } ^ { + 0.0022 }$ d and radius of $3.41 _{ - 0.12 } ^ { + 0.14 }$ R for the inner planet, and a period in the range 19.26-35 d and a radius of $5.83 _{ - 0.14 } ^ { + 0.14 }$ R for the outer planet, which was only seen to transit once. Each signal was independently statistically validated, taking into consideration the TESS light curve as well as the ground-based spectroscopic follow-up observations. Radial velocities from HARPS-N and EXPRES yield a tentative detection of planet b, whose mass we estimate to be $11.56 _{ - 6.14 } ^ { + 6.58 }$ M, and allow us to place an upper limit of 27.5 M (99 per cent confidence) on the mass of planet c. Due to the brightness of the host star and the strong likelihood of an extended H/He atmosphere on both planets, this system offers excellent prospects for atmospheric characterization and comparative planetology.
  46. Grayling, M., Gutierrez, C., Sullivan, M., Wiseman, P., Vincenzi, M., Gonzalez-Gaitan, S., Tucker, B., Galbany, L., Kelsey, L., Lidman, C., Swann, E., Smith, M., Frohmaier, C., Carollo, D., Glazebrook, K., Lewis, G., Moller, A., Hinton, S., Uddin, S., Abbott, T., Aguena, M., Avila, S., Bertin, E., Bhargava, S., Brooks, D., Carnero Rosell, A., Carrasco Kind, M., Carretero, J., Costanzi, M., da Costa, L., De Vicente, J., Desai, S., Diehl, H., Doel, P., Everett, S., Ferrero, I., Fosalba, P., Frieman, J., Garcia-Bellido, J., Gaztanaga, E., Gruen, D., Gruendl, R., Gschwend, J., Gutierrez, G., Hoyle, B., Kuehn, K., Kuropatkin, N., Lima, M., MacCrann, N., Marshall, J., Martini, P., Miquel, R., Morgan, R., Palmese, A., Paz-Chinchon, F., Plazas, A., Romer, A., Sanchez, C., Sanchez, E., Scarpine, V., Serrano, S., Sevilla-Noarbe, I., Soares-Santos, M., Suchyta, E., Tarle, G., Thomas, D., To, C., Varga, T., Walker, A., Wilkinson, R., DES Collaboration, 2021, MNRAS, 505, 3950, Understanding the extreme luminosity of DES14X2fna
    We present DES14X2fna, a high-luminosity, fast-declining Type IIb supernova (SN IIb) at redshift z = 0.0453, detected by the Dark Energy Survey (DES). DES14X2fna is an unusual member of its class, with a light curve showing a broad, luminous peak reaching Mr -19.3 mag 20 d after explosion. This object does not show a linear decline tail in the light curve until 60 d after explosion, after which it declines very rapidly (4.30 0.10 mag 100 d-1 in the r band). By fitting semi-analytic models to the photometry of DES14X2fna, we find that its light curve cannot be explained by a standard 56Ni decay model as this is unable to fit the peak and fast tail decline observed. Inclusion of either interaction with surrounding circumstellar material or a rapidly-rotating neutron star (magnetar) significantly increases the quality of the model fit. We also investigate the possibility for an object similar to DES14X2fna to act as a contaminant in photometric samples of SNe Ia for cosmology, finding that a similar simulated object is misclassified by a recurrent neural network (RNN)-based photometric classifier as an SN Ia in ~1.1-2.4 per cent of cases in DES, depending on the probability threshold used for a positive classification.
  47. Williams, P., Morrell, N., Boutsia, K., Massey, P., 2021, MNRAS, 505, 5029, The episodic dust-making Wolf-Rayet star HD 38030 in the Large Magellanic Cloud
    Mid-infrared photometry of the Wolf-Rayet star HD 38030 in the Large Magellanic Cloud from the NEOWISE-R mission show it to have undergone a dust-formation episode in 2018 and the dust to have cooled in 2019-20. New spectroscopy with the MagE spectrograph on the Magellan I Baade Telescope in 2019 and 2020 show absorption lines attributable to a companion of type near O9.7III-IV. We found a significant shift in the radial velocity of the C IV 5801-12 blend compared with the RVs measured in 1984 and 1993. The results combine to suggest that HD 38030 is a colliding-wind binary having short-lived dust formation episodes, like the Galactic systems WR 140 and WR 19, but at intervals in excess of 20 yr.
  48. Morgado, B., Sicardy, B., Braga-Ribas, F., Desmars, J., Gomes-Junior, A., Berard, D., Leiva, R., Ortiz, J., Vieira-Martins, R., Benedetti-Rossi, G., Santos-Sanz, P., Camargo, J., Duffard, R., Rommel, F., Assafin, M., Boufleur, R., Colas, F., Kretlow, M., Beisker, W., Sfair, R., Snodgrass, C., Morales, N., Fernandez-Valenzuela, E., Amaral, L., Amarante, A., Artola, R., Backes, M., Bath, K., Bouley, S., Buie, M., Cacella, P., Colazo, C., Colque, J., Dauvergne, J., Dominik, M., Emilio, M., Erickson, C., Evans, R., Fabrega-Polleri, J., Garcia-Lambas, D., Giacchini, B., Hanna, W., Herald, D., Hesler, G., Hinse, T., Jacques, C., Jehin, E., Jrgensen, U., Kerr, S., Kouprianov, V., Levine, S., Linder, T., Maley, P., Machado, D., Maquet, L., Maury, A., Melia, R., Meza, E., Mondon, B., Moura, T., Newman, J., Payet, T., Pereira, C., Pollock, J., Poltronieri, R., Quispe-Huaynasi, F., Reichart, D., de Santana, T., Schneiter, E., Sieyra, M., Skottfelt, J., Soulier, J., Starck, M., Thierry, P., Torres, P., Trabuco, L., Unda-Sanzana, E., Yamashita, T., Winter, O., Zapata, A., Zuluaga, C., 2021, A&A, 652, A141, Refined physical parameters for Chariklo's body and rings from stellar occultations observed between 2013 and 2020
    Context. The Centaur (10199) Chariklo has the first ring system discovered around a small object. It was first observed using stellar occultation in 2013. Stellar occultations allow sizes and shapes to be determined with kilometre accuracy, and provide the characteristics of the occulting object and its vicinity.
    Aims: Using stellar occultations observed between 2017 and 2020, our aim is to constrain the physical parameters of Chariklo and its rings. We also determine the structure of the rings, and obtain precise astrometrical positions of Chariklo.
    Methods: We predicted and organised several observational campaigns of stellar occultations by Chariklo. Occultation light curves were measured from the datasets, from which ingress and egress times, and the ring widths and opacity values were obtained. These measurements, combined with results from previous works, allow us to obtain significant constraints on Chariklo's shape and ring structure.
    Results: We characterise Chariklo's ring system (C1R and C2R), and obtain radii and pole orientations that are consistent with, but more accurate than, results from previous occultations. We confirm the detection of W-shaped structures within C1R and an evident variation in radial width. The observed width ranges between 4.8 and 9.1 km with a mean value of 6.5 km. One dual observation (visible and red) does not reveal any differences in the C1R opacity profiles, indicating a ring particle size larger than a few microns. The C1R ring eccentricity is found to be smaller than 0.022 (3), and its width variations may indicate an eccentricity higher than ~0.005. We fit a tri-axial shape to Chariklo's detections over 11 occultations, and determine that Chariklo is consistent with an ellipsoid with semi-axes of 143.81.5+1.4, 135.22.8+1.4, and 99.12.7+5.4 km. Ultimately, we provided seven astrometric positions at a milliarcsecond accuracy level, based on Gaia EDR3, and use it to improve Chariklo's ephemeris.

    Tables C.1, C.2 and lightcurves are only available at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (ftp://130.79.128.5) or via http://cdsarc.u-strasbg.fr/viz-bin/cat/J/A+A/652/A141.

  49. Lester, K., Matson, R., Howell, S., Furlan, E., Gnilka, C., Scott, N., Ciardi, D., Everett, M., Hartman, Z., Hirsch, L., 2021, AJ, 162, 75, Speckle Observations of TESS Exoplanet Host Stars. II. Stellar Companions at 1-1000 au and Implications for Small Planet Detection
    We present high-angular-resolution imaging observations of 517 host stars of TESS exoplanet candidates using the 'Alopeke and Zorro speckle cameras at Gemini North and South. The sample consists mainly of bright F, G, K stars at distances of less than 500 pc. Our speckle observations span angular resolutions of ~20 mas out to 1"2, yielding spatial resolutions of <10-500 au for most stars, and our contrast limits can detect companion stars 5-9 mag fainter than the primary at optical wavelengths. We detect 102 close stellar companions and determine the separation, magnitude difference, mass ratio, and estimated orbital period for each system. Our observations of exoplanet host star binaries reveal that they have wider separations than field binaries, with a mean orbital semimajor axis near 100 au. Other imaging studies have suggested this dearth of very closely separated binaries in systems which host exoplanets, but incompleteness at small separations makes it difficult to disentangle unobserved companions from a true lack of companions. With our improved angular resolution and sensitivity, we confirm that this lack of close exoplanet host binaries is indeed real. We also search for a correlation between planetary orbital radii versus binary star separation; but, given the very short orbital periods of the TESS planets, we do not find any clear trend. We do note that in exoplanet systems containing binary host stars, there is an observational bias against detecting Earth-size planet transits due to transit depth dilution caused by the companion star.
  50. Cantu, S., Pace, A., Marshall, J., Strigari, L., Crnojevic, D., Simon, J., Drlica-Wagner, A., Bechtol, K., Martinez-Vazquez, C., Santiago, B., Amara, A., Stringer, K., Diehl, H., Aguena, M., Allam, S., Avila, S., Brooks, D., Carnero Rosell, A., Carrasco Kind, M., Carretero, J., Costanzi, M., Da Costa, L., De Vicente, J., Desai, S., Doel, P., Eifler, T., Everett, S., Frieman, J., Garcia-Bellido, J., Gaztanaga, E., Gruen, D., Gruendl, R., Gschwend, J., Gutierrez, G., Hinton, S., Hollowood, D., Honscheid, K., James, D., Kuehn, K., Maia, M., Menanteau, F., Miquel, R., Palmese, A., Paz-Chinchon, F., Plazas, A., Sanchez, E., Scarpine, V., Schubnell, M., Serrano, S., Sevilla-Noarbe, I., Smith, M., Soares-Santos, M., Suchyta, E., Swanson, M., Tarle, G., Walker, A., Wilkinson, R., DES Collaboration, 2021, ApJ, 916, 81, A Deeper Look at DES Dwarf Galaxy Candidates: Grus I and Indus II
    We present deep g- and r-band Magellan/Megacam photometry of two dwarf galaxy candidates discovered in the Dark Energy Survey (DES), Grus I and Indus II (DES J2038-4609). For the case of Grus I, we resolved the main sequence turn-off (MSTO) and ~2 mags below it. The MSTO can be seen at g0 ~ 24 with a photometric uncertainty of 0.03 mag. We show Grus I to be consistent with an old, metal-poor (~13.3 Gyr, [Fe/H] ~ -1.9) dwarf galaxy. We derive updated distance and structural parameters for Grus I using this deep, uniform, wide-field data set. We find an azimuthally-averaged halflight radius more than two times larger (~151+21-31 pc; ~ $4\buildrel{\,\prime}\over{.} {16}_{-0.74}^{+0.54}$) and an absolute V-band magnitude ~-4.1 that is ~1 magnitude brighter than previous studies. We obtain updated distance, ellipticity, and centroid parameters that are in agreement with other studies within uncertainties. Although our photometry of Indus II is ~2-3 magnitudes deeper than the DES Y1 public release, we find no coherent stellar population at its reported location. The original detection was located in an incomplete region of sky in the DES Y2Q1 data set and was flagged due to potential blue horizontal branch member stars. The best-fit isochrone parameters are physically inconsistent with both dwarf galaxies and globular clusters. We conclude that Indus II is likely a false positive, flagged due to a chance alignment of stars along the line of sight.
  51. Vincenzi, M., Sullivan, M., Graur, O., Brout, D., Davis, T., Frohmaier, C., Galbany, L., Gutierrez, C., Hinton, S., Hounsell, R., Kelsey, L., Kessler, R., Kovacs, E., Kuhlmann, S., Lasker, J., Lidman, C., Moller, A., Nichol, R., Sako, M., Scolnic, D., Smith, M., Swann, E., Wiseman, P., Asorey, J., Lewis, G., Sharp, R., Tucker, B., Aguena, M., Allam, S., Avila, S., Bertin, E., Brooks, D., Burke, D., Carnero Rosell, A., Carrasco Kind, M., Carretero, J., Castander, F., Choi, A., Costanzi, M., da Costa, L., Pereira, M., De Vicente, J., Desai, S., Diehl, H., Doel, P., Everett, S., Ferrero, I., Fosalba, P., Frieman, J., Garcia-Bellido, J., Gaztanaga, E., Gerdes, D., Gruen, D., Gruendl, R., Gutierrez, G., Hollowood, D., Honscheid, K., Hoyle, B., James, D., Kuehn, K., Kuropatkin, N., Maia, M., Martini, P., Menanteau, F., Miquel, R., Morgan, R., Palmese, A., Paz-Chinchon, F., Plazas, A., Romer, A., Sanchez, E., Scarpine, V., Serrano, S., Sevilla-Noarbe, I., Soares-Santos, M., Suchyta, E., Tarle, G., Thomas, D., To, C., Varga, T., Walker, A., Wilkinson, R., DES Collaboration, 2021, MNRAS, 505, 2819, The Dark Energy Survey supernova programme: modelling selection efficiency and observed core-collapse supernova contamination
    The analysis of current and future cosmological surveys of Type Ia supernovae (SNe Ia) at high redshift depends on the accurate photometric classification of the SN events detected. Generating realistic simulations of photometric SN surveys constitutes an essential step for training and testing photometric classification algorithms, and for correcting biases introduced by selection effects and contamination arising from core-collapse SNe in the photometric SN Ia samples. We use published SN time-series spectrophotometric templates, rates, luminosity functions, and empirical relationships between SNe and their host galaxies to construct a framework for simulating photometric SN surveys. We present this framework in the context of the Dark Energy Survey (DES) 5-yr photometric SN sample, comparing our simulations of DES with the observed DES transient populations. We demonstrate excellent agreement in many distributions, including Hubble residuals, between our simulations and data. We estimate the core collapse fraction expected in the DES SN sample after selection requirements are applied and before photometric classification. After testing different modelling choices and astrophysical assumptions underlying our simulation, we find that the predicted contamination varies from 7.2 to 11.7 per cent, with an average of 8.8 per cent and an r.m.s. of 1.1 per cent. Our simulations are the first to reproduce the observed photometric SN and host galaxy properties in high-redshift surveys without fine-tuning the input parameters. The simulation methods presented here will be a critical component of the cosmology analysis of the DES photometric SN Ia sample: correcting for biases arising from contamination, and evaluating the associated systematic uncertainty.
  52. Abbott, T., Adamow, M., Aguena, M., Allam, S., Amon, A., Annis, J., Avila, S., Bacon, D., Banerji, M., Bechtol, K., Becker, M., Bernstein, G., Bertin, E., Bhargava, S., Bridle, S., Brooks, D., Burke, D., Carnero Rosell, A., Carrasco Kind, M., Carretero, J., Castander, F., Cawthon, R., Chang, C., Choi, A., Conselice, C., Costanzi, M., Crocce, M., da Costa, L., Davis, T., De Vicente, J., DeRose, J., Desai, S., Diehl, H., Dietrich, J., Drlica-Wagner, A., Eckert, K., Elvin-Poole, J., Everett, S., Evrard, A., Ferrero, I., Ferte, A., Flaugher, B., Fosalba, P., Friedel, D., Frieman, J., Garcia-Bellido, J., Gaztanaga, E., Gelman, L., Gerdes, D., Giannantonio, T., Gill, M., Gruen, D., Gruendl, R., Gschwend, J., Gutierrez, G., Hartley, W., Hinton, S., Hollowood, D., Honscheid, K., Huterer, D., James, D., Jeltema, T., Johnson, M., Kent, S., Kron, R., Kuehn, K., Kuropatkin, N., Lahav, O., Li, T., Lidman, C., Lin, H., MacCrann, N., Maia, M., Manning, T., Maloney, J., March, M., Marshall, J., Martini, P., Melchior, P., Menanteau, F., Miquel, R., Morgan, R., Myles, J., Neilsen, E., Ogando, R., Palmese, A., Paz-Chinchon, F., Petravick, D., Pieres, A., Plazas, A., Pond, C., Rodriguez-Monroy, M., Romer, A., Roodman, A., Rykoff, E., Sako, M., Sanchez, E., Santiago, B., Scarpine, V., Serrano, S., Sevilla-Noarbe, I., Smith, J., Smith, M., Soares-Santos, M., Suchyta, E., Swanson, M., Tarle, G., Thomas, D., To, C., Tremblay, P., Troxel, M., Tucker, D., Turner, D., Varga, T., Walker, A., Wechsler, R., Weller, J., Wester, W., Wilkinson, R., Yanny, B., Zhang, Y., Nikutta, R., Fitzpatrick, M., Jacques, A., Scott, A., Olsen, K., Huang, L., Herrera, D., Juneau, S., Nidever, D., Weaver, B., Adean, C., Correia, V., de Freitas, M., Freitas, F., Singulani, C., Vila-Verde, G., Linea Science Server, 2021, ApJS, 255, 20, The Dark Energy Survey Data Release 2
    We present the second public data release of the Dark Energy Survey, DES DR2, based on optical/near-infrared imaging by the Dark Energy Camera mounted on the 4 m Blanco telescope at Cerro Tololo Inter-American Observatory in Chile. DES DR2 consists of reduced single-epoch and coadded images, a source catalog derived from coadded images, and associated data products assembled from 6 yr of DES science operations. This release includes data from the DES wide-area survey covering ~5000 deg2 of the southern Galactic cap in five broad photometric bands, grizY. DES DR2 has a median delivered point-spread function FWHM of g = 1.11, r = 0.95, i = 0.88, z = 0.83, and Y = 0"90, photometric uniformity with a standard deviation of < 3 mmag with respect to Gaia DR2 G band, a photometric accuracy of ~11 mmag, and a median internal astrometric precision of ~27 mas. The median coadded catalog depth for a 1"95 diameter aperture at signal-to-noise ratio = 10 is g = 24.7, r = 24.4, i = 23.8, z = 23.1, and Y = 21.7 mag. DES DR2 includes ~691 million distinct astronomical objects detected in 10,169 coadded image tiles of size 0.534 deg2 produced from 76,217 single-epoch images. After a basic quality selection, benchmark galaxy and stellar samples contain 543 million and 145 million objects, respectively. These data are accessible through several interfaces, including interactive image visualization tools, web-based query clients, image cutout servers, and Jupyter notebooks. DES DR2 constitutes the largest photometric data set to date at the achieved depth and photometric precision.
  53. Hansen, C., Castillo-Rogez, J., Grundy, W., Hofgartner, J., Martin, E., Mitchell, K., Nimmo, F., Nordheim, T., Paty, C., Quick, L., Roberts, J., Runyon, K., Schenk, P., Stern, A., Umurhan, O., 2021, PSJ, 2, 137, Triton: Fascinating Moon, Likely Ocean World, Compelling Destination!
    Triton is an important signpost in understanding the diverse populations of both ocean worlds and Kuiper Belt objects (KBOs). As a likely ocean world, it is unique by virtue of its kidnapped history from the Kuiper Belt: its large orbital inclination makes it the only ocean world thought to be primarily heated by obliquity tides. It is volatile-rich due to its formation in the outer solar system and its unusual surface geology may be the product of cryovolcanism. Observations from New Horizons and Cassini motivate re-examination of Triton data sets and models, with value for comparative planetology of ocean worlds and KBOs, most notably with Europa, Enceladus, Titan, and Pluto. We re-explore old data sets with the new perspective of the importance of ocean worlds in our solar system and the search for life.
  54. Hirabayashi, M., Mimasu, Y., Sakatani, N., Watanabe, S., Tsuda, Y., Saiki, T., Kikuchi, S., Kouyama, T., Yoshikawa, M., Tanaka, S., Nakazawa, S., Takei, Y., Terui, F., Takeuchi, H., Fujii, A., Iwata, T., Tsumura, K., Matsuura, S., Shimaki, Y., Urakawa, S., Ishibashi, Y., Hasegawa, S., Ishiguro, M., Kuroda, D., Okumura, S., Sugita, S., Okada, T., Kameda, S., Kamata, S., Higuchi, A., Senshu, H., Noda, H., Matsumoto, K., Suetsugu, R., Hirai, T., Kitazato, K., Farnocchia, D., Naidu, S., Tholen, D., Hergenrother, C., Whiteley, R., Moskovitz, N., Abell, P., the Hayabusa2 extended mission study Group, 2021, AdSpR, 68, 1533, Hayabusa2 extended mission: New voyage to rendezvous with a small asteroid rotating with a short period
    Hayabusa2 is the Japanese Asteroid Return Mission and targeted the carbonaceous asteroid Ryugu, conducted by the Japan Aerospace Exploration Agency (JAXA). The goal of this mission was to conduct proximity operations including remote sensing observations, material sampling, and a Small Carry-On Impact experiment, as well as sample analyses. As of September 2020, the spacecraft is on the way back to Earth with samples from Ryugu with no critical issues after the successful departure in November 2019. Here, we propose an extended mission in which the spacecraft will rendezvous with a small asteroid with ~30 m - ~40 m in diameter that is rotating at a spin period of ~10 min after an additional ~10-year cruise phase. We introduce that two scenarios are suitable for the extended mission. In the first scenario, the spacecraft will perform swing-by maneuvers at Venus once and Earth twice to arrive at asteroid 2001 AV43. In the second scenario, it will perform swing-by maneuvers at Earth twice to reach asteroid 1998 KY26. In both scenarios, the mission will continue until the early 2030s. JAXA recently released the decision that the spacecraft will rendezvous with 1998 KY26. This paper focuses on our scientific assessments of the two scenarios but leaves the decision process to go to 1998 KY26 for future reports. Rendezvous operations will be planned to detail the physical properties and surrounding environments of the target, one of the smallest elements of small planetary bodies. By achieving the planned operations, the mission will provide critical hints on the violent histories of collisions and accumulations of small bodies in the solar system. Furthermore, the established scientific knowledge and techniques will advance key technologies for planetary defense.
  55. Devogele, M., Ferrais, M., Jehin, E., Moskovitz, N., Skiff, B., Levine, S., Gustafsson, A., Farnocchia, D., Micheli, M., Snodgrass, C., Borisov, G., Manfroid, J., Moulane, Y., Benkhaldoun, Z., Burdanov, A., Pozuelos, F., Gillon, M., de Wit, J., Green, S., Bendjoya, P., Rivet, J., Abe, L., Vernet, D., Chandler, C., Trujillo, C., 2021, MNRAS, 505, 245, (6478) Gault: physical characterization of an active main-belt asteroid
    In 2018 December, the main-belt asteroid (6478) Gault was reported to display activity. Gault is an asteroid belonging to the Phocaea dynamical family and was not previously known to be active, nor was any other member of the Phocaea family. In this work, we present the results of photometric and spectroscopic observations that commenced soon after the discovery of activity. We obtained observations over two apparitions to monitor its activity, rotation period, composition, and possible non-gravitational orbital evolution. We find that Gault has a rotation period of P = 2.4929 0.0003 h with a light-curve amplitude of 0.06 magnitude. This short rotation period close to the spin barrier limit is consistent with Gault having a density no smaller than = 1.85 g cm-3 and its activity being triggered by the YORP (Yarkovsky-O'Keefe-Radzievskii-Paddack) spin-up mechanism. Analysis of the Gault phase curve over phase angles ranging from 0.4 to 23.6 provides an absolute magnitude of H = 14.81 0.04, G1 = 0.25 0.07, and G2 = 0.38 0.04. Model fits to the phase curve find the surface regolith grain size constrained between 100 and 500 $\rm {\mu }$m. Using relations between the phase curve and albedo, we determine that the geometrical albedo of Gault is pv = 0.26 0.05 corresponding to an equivalent diameter of $D = 2.8^{+0.4}_{-0.2}$ km. Our spectroscopic observations are all consistent with an ordinary chondrite-like composition (S, or Q-type in the Bus-DeMeo taxonomic classification). A search through archival photographic plate surveys found previously unidentified detections of Gault dating back to 1957 and 1958. Only the latter had been digitized, which we measured to nearly double the observation arc of Gault. Finally, we did not find any signal of activity during the 2020 apparition or non-gravitational effects on its orbit.
  56. Clark, C., van Belle, G., Ciardi, D., Lund, M., Howell, S., Everett, M., 2021, tsc2, 112, A Dearth of Stellar Companions to M-dwarf TESS Objects of Interest
    The Transiting Exoplanet Survey Satellite (TESS) has proven to be a powerful resource for uncovering planets, and the M-dwarfs have been established as favorable planet hosts. It has also become apparent that stellar multiplicity has wide-ranging implications for exoplanet detection and characterization, and that speckle imaging is one of the most efficient tools for probing these multi-star systems. We therefore present high-resolution imaging observations of 63 M-dwarf TOIs using speckle imagers at the 3.5m WIYN telescope, the 4.3m Lowell Discovery Telescope (LDT), and the twin 8.1m Gemini North and South telescopes. However, only one companion was detected. This finding is in contrast to the established 46% binarity rate in exoplanet host stars and the established 27% stellar multiplicity rate for field M-dwarfs. These results indicate that M-dwarf TOIs have a much lower multiplicity rate than field M-dwarfs. Our observations also imply that planet signals detected from M-dwarf TOIs are more likely to be real than those from higher-mass stars. Finally, these data support the observation that exoplanet-hosting binary stars have, in general, wider separations than field binaries.
  57. Claytor, Z., van Saders, J., Llama, J., Sadowski, P., Quach, B., Avallone, E., 2021, tsc2, 120, Recovery of TESS Stellar Rotation Periods with Deep Learning
    TESS is poised to increase the number of stellar rotation period estimates by an order of magnitude, but the mission's systematics have complicated period searches. While several efforts attempt to solve this problem by removing systematics, standard methods of data reduction have shown limited success. Here I present a method to predict rotation periods from TESS full-frame image light curves using deep learning. This method relies on a training set of simulated light curves convolved with TESS galaxy light curves to emulate the instrumental noise and systematics observed in stellar signals. The simulations include surface differential rotation, spot evolution, and activity level to make the light curves as realistic as possible. Our approach allows the network to learn the difference between rotation signals and TESS systematics. With the added ability to predict uncertainty in the period, we can determine what regions of parameter space the predictions are most credible, producing a reliable set of rotation periods. I present the first set of rotation periods obtained with this method and explore TESS's insights to stellar structure and evolution through the lens of rotation.
  58. Li, Y., Pang, X., Tang, S., 2021, RNAAS, 5, 173, Evidence of Early-stage Tidal Structures of Open Clusters Revealed by Kinematics with Gaia EDR3
    Blanco 1, a 100 Myr open cluster in the solar neighborhood, is well known for its two 50 pc long tidal tails. Taking Blanco 1 as a reference, we find evidence of early-stage tidal disruption in two other open clusters of ~120 Myr: the Pleiades and NGC 2516, via Gaia EDR3 data. These two clusters have a total mass of 2-6 times that of Blanco 1. Despite having a similar age as Blanco 1, the Pleiades and NGC 2516 have a larger fraction of their members bound: 86% of their mass is inside the tidal radius, versus 63% for Blanco 1. However, a correlation between Blanco 1's 50 pc long tidal tails and the "kinematic tails" in velocity space is also found for the Pleiades and NGC 2516. This evidence supports the idea that the modest elongation seen in the spatial distribution for the Pleiades and NGC 2516 is a result of early-stage tidal disruption.
  59. Hansen, T., Ji, A., Da Costa, G., Li, T., Casey, A., Pace, A., Cullinane, L., Erkal, D., Koposov, S., Kuehn, K., Lewis, G., Mackey, D., Simpson, J., Shipp, N., Zucker, D., Bland-Hawthorn, J., S5 Collaboration, 2021, ApJ, 915, 103, S5: The Destruction of a Bright Dwarf Galaxy as Revealed by the Chemistry of the Indus Stellar Stream
    The recently discovered Indus stellar stream exhibits a diverse chemical signature compared to what is found for most other streams due to the abundances of two outlier stars, Indus_0 and Indus_13. Indus_13 exhibits an extreme enhancement in rapid neutron-capture (r-)process elements with [Eu/Fe] = + 1.81. It thus provides direct evidence of the accreted nature of r-process-enhanced stars. In this paper we present a detailed chemical analysis of the neutron-capture elements in Indus_13, revealing the star to be slightly actinide poor. The other outlier, Indus_0, displays a globular cluster-like signature with high N, Na, and Al abundances, while the rest of the Indus stars show abundances compatible with a dwarf galaxy origin. Hence, Indus_0 provides the first chemical evidence of a fully disrupted dwarf containing a globular cluster. We use the chemical signature of the Indus stars to discuss the nature of the stream progenitor which was likely a chemically evolved system, with a mass somewhere in the range from Ursa Minor to Fornax.
  60. Linscott, I., Bird, M., Hinson, D., Tyler, G., Stern, S., Vincent, M., DeBoy, C., Young, L., Beyer, R., Ennico, K., Moore, J., Gladstone, G., Olkin, C., Patzold, M., Schenk, P., Strobel, D., Summers, M., Weaver, H., Grundy, W., New Horizons Science Team, 2021, Icar, 363, 114430, High-resolution radiometry of Pluto at 4.2 cm with New Horizons
    The radio thermal emission from Pluto was observed from the New Horizons spacecraft at a wavelength of 4.2 cm along two scans across the planetary disk shortly after closest approach to Pluto on 14 July 2015. The measurements were performed as part of the New Horizons Radio Science Experiment (REX) using the 2.1 m High Gain Antenna (HGA) and the spacecraft's X-Band receiver. The HGA boresight first scanned along a diametric chord across the Pluto disk and then reversed direction to traverse a chord that crossed close to Pluto's winter pole. The diametric scan reveals a "hot spot" on the Pluto nightside associated with an optically bright region centered roughly at the planetocentric coordinates 280 E, 55 S, imaged in 2002-03 with the Hubble Space Telescope. The nightside was also found to be warmer than the dayside during the polar scan. The highest emission was not observed at the maximum southern latitude, however, but rather near the outbound Pluto limb at lower latitude. The REX emission profile from the polar scan is qualitatively consistent with a bright U-shaped polar cap observed on Pluto's Charon-facing hemisphere during the recurring Pluto/Charon mutual events in the late 1980's. The REX radiometer measurements show distinct variations in microwave brightness that constrain volatile transport models and provide unique information on the thermal structure and composition on the regions in winter night during the New Horizons encounter at Pluto.
  61. Grandis, S., Mohr, J., Costanzi, M., Saro, A., Bocquet, S., Klein, M., Aguena, M., Allam, S., Annis, J., Ansarinejad, B., Bacon, D., Bertin, E., Bleem, L., Brooks, D., Burke, D., Carnero Rosel, A., Carrasco Kind, M., Carretero, J., Castander, F., Choi, A., da Costa, L., De Vincente, J., Desai, S., Diehl, H., Dietrich, J., Doel, P., Eifler, T., Everett, S., Ferrero, I., Floyd, B., Fosalba, P., Frieman, J., Garcia-Bellido, J., Gaztanaga, E., Gruen, D., Gruendl, R., Gschwend, J., Gupta, N., Gutierrez, G., Hinton, S., Hollowood, D., Honscheid, K., James, D., Jeltema, T., Kuehn, K., Lahav, O., Lidman, C., Lima, M., Maia, M., March, M., Marshall, J., Melchior, P., Menanteau, F., Miquel, R., Morgan, R., Myles, J., Ogando, R., Palmese, A., Paz-Chinchon, F., Plazas, A., Reichardt, C., Romer, A., Sanchez, E., Scarpine, V., Serrano, S., Sevilla-Noarbe, I., Singh, P., Smith, M., Suchyta, E., Swanson, M., Tarle, G., Thomas, D., To, C., Weller, J., Wilkinson, R., Wu, H., 2021, MNRAS, 504, 1253, Exploring the contamination of the DES-Y1 cluster sample with SPT-SZ selected clusters
    We perform a cross validation of the cluster catalogue selected by the red-sequence Matched-filter Probabilistic Percolation algorithm (redMaPPer) in Dark Energy Survey year 1 (DES-Y1) data by matching it with the Sunyaev-Zel'dovich effect (SZE) selected cluster catalogue from the South Pole Telescope SPT-SZ survey. Of the 1005 redMaPPer selected clusters with measured richness $\hat{\lambda }\gt 40$ in the joint footprint, 207 are confirmed by SPT-SZ. Using the mass information from the SZE signal, we calibrate the richness-mass relation using a Bayesian cluster population model. We find a mass trend MB consistent with a linear relation (B ~ 1), no significant redshift evolution and an intrinsic scatter in richness of = 0.22 0.06. By considering two error models, we explore the impact of projection effects on the richness-mass modelling, confirming that such effects are not detectable at the current level of systematic uncertainties. At low richness SPT-SZ confirms fewer redMaPPer clusters than expected. We interpret this richness dependent deficit in confirmed systems as due to the increased presence at low richness of low-mass objects not correctly accounted for by our richness-mass scatter model, which we call contaminants. At a richness $\hat{\lambda }=40$, this population makes up ${\gt}12{{\ \rm per\ cent}}$ (97.5 percentile) of the total population. Extrapolating this to a measured richness $\hat{\lambda }=20$ yields ${\gt}22{{\ \rm per\ cent}}$ (97.5 percentile). With these contamination fractions, the predicted redMaPPer number counts in different plausible cosmologies are compatible with the measured abundance. The presence of such a population is also a plausible explanation for the different mass trends (B ~ 0.75) obtained from mass calibration using purely optically selected clusters. The mean mass from stacked weak lensing (WL) measurements suggests that these low-mass contaminants are galaxy groups with masses ~3-5 1013 M which are beyond the sensitivity of current SZE and X-ray surveys but a natural target for SPT-3G and eROSITA.
  62. Engle, A., Hanley, J., Dustrud, S., Thompson, G., Lindberg, G., Grundy, W., Tegler, S., 2021, PSJ, 2, 118, Phase Diagram for the Methane-Ethane System and Its Implications for Titan's Lakes
    On Titan, methane (CH4) and ethane (C2H6) are the dominant species found in the lakes and seas. In this study, we have combined laboratory work and modeling to refine the methane-ethane binary phase diagram at low temperatures and probe how the molecules interact at these conditions. We used visual inspection for the liquidus and Raman spectroscopy for the solidus. Through these methods, we determined a eutectic point of 71.15 0.5 K at a composition of 0.644 0.018 methane-0.356 0.018 ethane mole fraction from the liquidus data. Using the solidus data, we found a eutectic isotherm temperature of 72.2 K with a standard deviation of 0.4 K. In addition to mapping the binary system, we looked at the solid-solid transitions of pure ethane and found that, when cooling, the transition of solid I-III occurred at 89.45 0.2 K. The warming sequence showed transitions of solid III-II occurring at 89.85 0.2 K and solid II-I at 89.65 0.2 K. Ideal predictions were compared with molecular dynamics simulations to reveal that the methane-ethane system behaves almost ideally, and the largest deviations occur as the mixing ratio approaches the eutectic composition.
  63. MacKenzie, S., Birch, S., Horst, S., Sotin, C., Barth, E., Lora, J., Trainer, M., Corlies, P., Malaska, M., Sciamma-O'Brien, E., Thelen, A., Turtle, E., Radebaugh, J., Hanley, J., Solomonidou, A., Newman, C., Regoli, L., Rodriguez, S., Seignovert, B., Hayes, A., Journaux, B., Steckloff, J., Nna-Mvondo, D., Cornet, T., Palmer, M., Lopes, R., Vinatier, S., Lorenz, R., Nixon, C., Czaplinski, E., Barnes, J., Sittler, E., Coates, A., 2021, PSJ, 2, 112, Titan: Earth-like on the Outside, Ocean World on the Inside
    Thanks to the Cassini-Huygens mission, Titan, the pale orange dot of Pioneer and Voyager encounters, has been revealed to be a dynamic, hydrologically shaped, organic-rich ocean world offering unparalleled opportunities to explore prebiotic chemistry. And while Cassini-Huygens revolutionized our understanding of each of the three "layers" of Titanthe atmosphere, the surface, and the interiorwe are only beginning to hypothesize how these realms interact. In this paper, we summarize the current state of Titan knowledge and discuss how future exploration of Titan would address some of the next decade's most compelling planetary science questions. We also demonstrate why exploring Titan, both with and beyond the Dragonfly New Frontiers mission, is a necessary and complementary component of an Ocean Worlds Program that seeks to understand whether habitable environments exist elsewhere in our solar system.
  64. Chen, A., Huterer, D., Lee, S., Ferte, A., Weaverdyck, N., Alves, O., Leonard, C., MacCrann, N., Raveri, M., Porredon, A., Di Valentino, E., Muir, J., Lemos, P., Liddle, A., Blazek, J., Campos, A., Cawthon, R., Choi, A., Dodelson, S., Elvin-Poole, J., Gruen, D., Ross, A., Secco, L., Sevilla-Noarbe, I., Sheldon, E., Troxel, M., Zuntz, J., Abbott, T., Aguena, M., Allam, S., Annis, J., Avila, S., Bertin, E., Bhargava, S., Bridle, S., Brooks, D., Carnero Rosell, A., Carrasco Kind, M., Carretero, J., Costanzi, M., Crocce, M., da Costa, L., Pereira, M., Davis, T., Doel, P., Eifler, T., Ferrero, I., Fosalba, P., Frieman, J., Garcia-Bellido, J., Gaztanaga, E., Gerdes, D., Gruendl, R., Gschwend, J., Gutierrez, G., Hinton, S., Hollowood, D., Honscheid, K., Hoyle, B., James, D., Jarvis, M., Kuehn, K., Lahav, O., Maia, M., Marshall, J., Menanteau, F., Miquel, R., Morgan, R., Palmese, A., Paz-Chinchon, F., Plazas, A., Roodman, A., Sanchez, E., Scarpine, V., Schubnell, M., Serrano, S., Smith, M., Suchyta, E., Tarle, G., Thomas, D., To, C., Varga, T., Weller, J., Wilkinson, R., DES Collaboration, 2021, PhRvD, 103, 123528, Constraints on dark matter to dark radiation conversion in the late universe with DES-Y1 and external data
    We study a phenomenological class of models where dark matter converts to dark radiation in the low redshift epoch. This class of models, dubbed DMDR, characterizes the evolution of comoving dark-matter density with two extra parameters, and may be able to help alleviate the observed discrepancies between early and late-time probes of the Universe. We investigate how the conversion affects key cosmological observables such as the cosmic microwave background (CMB) temperature and matter power spectra. Combining 3x2pt data from Year 1 of the Dark Energy Survey, Planck-2018 CMB temperature and polarization data, supernovae (SN) Type Ia data from Pantheon, and baryon acoustic oscillation (BAO) data from BOSS DR12, MGS and 6dFGS, we place new constraints on the amount of dark matter that has converted to dark radiation and the rate of this conversion. The fraction of the dark matter that has converted since the beginning of the Universe in units of the current amount of dark matter, , is constrained at 68% confidence level to be <0.32 for DES-Y1 3x2pt data, <0.030 for CMB +SN +BAO data, and <0.037 for the combined dataset. The probability that the DES and CMB+SN+BAO datasets are concordant increases from 4% for the CDM model to 8% (less tension) for DMDR. The tension in S8=8{m/0.3 } between DES-Y1 3x2pt and CMB +SN +BAO is slightly reduced from 2.3 to 1.9 . We find no reduction in the Hubble tension when the combined data is compared to distance-ladder measurements in the DMDR model. The maximum-posterior goodness-of-fit statistics of DMDR and CDM model are comparable, indicating no preference for the DMDR cosmology over CDM .
  65. Tang, S., Stahl, A., Johns-Krull, C., Prato, L., Llama, J., 2021, JOSS, 6, 3095, IGRINS RV: A Python Package for Precision Radial Velocities with Near-Infrared Spectra
    No abstract found.
  66. Horch, E., Broderick, K., Casetti-Dinescu, D., Henry, T., Fekel, F., Muterspaugh, M., Willmarth, D., Winters, J., van Belle, G., Clark, C., Everett, M., 2021, AJ, 161, 295, Observations with the Differential Speckle Survey Instrument. X. Preliminary Orbits of K-dwarf Binaries and Other Stars
    This paper details speckle observations of binary stars taken at the Lowell Discovery Telescope, the WIYN telescope, and the Gemini telescopes between 2016 January and 2019 September. The observations taken at Gemini and Lowell were done with the Differential Speckle Survey Instrument (DSSI), and those done at WIYN were taken with the successor instrument to DSSI at that site, the NN-EXPLORE Exoplanet Star and Speckle Imager (NESSI). In total, we present 378 observations of 178 systems, and we show that the uncertainty in the measurement precision for the combined data set is 2 mas in separation, 1-2 in position angle depending on the separation, and 0.1 mag in magnitude difference. Together with data already in the literature, these new results permit 25 visual orbits and one spectroscopic-visual orbit to be calculated for the first time. In the case of the spectroscopic-visual analysis, which is done on the ternary star HD 173093, we calculate masses with a precision of better than 1% for all three stars in that system. Twenty-one of the visual orbits calculated have a K dwarf as the primary star; we add these to the known orbits of K-dwarf primary stars and discuss the basic orbital properties of these stars at this stage. Although incomplete, the data that exist so far indicate that binaries with K-dwarf primaries tend not to have low-eccentricity orbits at separations of one to a few tens of astronomical units, that is, on solar system scales.
  67. Stahl, A., Tang, S., Johns-Krull, C., Prato, L., Llama, J., Mace, G., Joon Lee, J., Oh, H., Luna, J., Jaffe, D., 2021, AJ, 161, 283, IGRINS RV: A Precision Radial Velocity Pipeline for IGRINS Using Modified Forward Modeling in the Near-infrared
    Application of the radial velocity (RV) technique in the near-infrared is valuable because of the diminished impact of stellar activity at longer wavelengths, making it particularly advantageous for the study of late-type stars but also for solar-type objects. In this paper, we present the IGRINS RV open-source python pipeline for computing infrared RV measurements from reduced spectra taken with IGRINS, an R / 45,000 spectrograph with simultaneous coverage of the H band (1.49-1.80 m) and K band (1.96-2.46 m). Using a modified forward-modeling technique, we construct high-resolution telluric templates from A0 standard observations on a nightly basis to provide a source of common-path wavelength calibration while mitigating the need to mask or correct for telluric absorption. Telluric standard observations are also used to model the variations in instrumental resolution across the detector, including a yearlong period when the K band was defocused. Without any additional instrument hardware, such as a gas cell or laser frequency comb, we are able to achieve precisions of 26.8 m s-1 in the K band and 31.1 m s-1 in the H band for narrow-line hosts. These precisions are empirically determined by a monitoring campaign of two RV standard stars, as well as the successful retrieval of planet-induced RV signals for both HD 189733 and Boo A; furthermore, our results affirm the presence of the Rossiter-McLaughlin effect for HD 189733. The IGRINS RV pipeline extends another important science capability to IGRINS, with publicly available software designed for widespread use. * https://github.com/shihyuntang/igrins_rv
  68. Benedict, G., Franz, O., Horch, E., Prato, L., Torres, G., McArthur, B., Wasserman, L., Latham, D., Stefanik, R., Latham, C., Skiff, B., 2021, AJ, 161, 285, Dissecting the Quadruple Binary Hyad vA 351 Masses for Three M Dwarfs and a White Dwarf
    We extend results first announced by Franz et al., that identified vA 351 = H346 in the Hyades as a multiple star system containing a white dwarf. With Hubble Space Telescope Fine Guidance Sensor fringe tracking and scanning, and more recent speckle observations, all spanning 20.7 years, we establish a parallax, relative orbit, and mass fraction for two components, with a period, $P=2.70\,\mathrm{yr}$ and total mass 2.1 ${{ \mathcal M }}_{\odot }$ . With ground-based radial velocities from the McDonald Observatory Otto Struve 2.1 m Telescope Sandiford Spectrograph, and Center for Astrophysics Digital Speedometers, spanning 37 years, we find that component B consists of BC, two M dwarf stars orbiting with a very short period ( ${P}_{\mathrm{BC}}=0.749$ days), having a mass ratio ${{ \mathcal M }}_{{\rm{C}}}$ / ${{ \mathcal M }}_{{\rm{B}}}$ = 0.95. We confirm that the total mass of the system can only be reconciled with the distance and component photometry by including a fainter, higher-mass component. The quadruple system consists of three M dwarfs (A, B, C) and one white dwarf (D). We determine individual M dwarf masses ${{ \mathcal M }}_{{\rm{A}}}$ = 0.53 0.10 ${{ \mathcal M }}_{\odot }$ , ${{ \mathcal M }}_{{\rm{B}}}$ = 0.43 0.04 ${{ \mathcal M }}_{\odot }$ , and ${{ \mathcal M }}_{{\rm{C}}}$ = 0.41 0.04 ${{ \mathcal M }}_{\odot }$ . The white dwarf mass, 0.54 0.04 ${{ \mathcal M }}_{\odot }$ , comes from cooling models, an assumed Hyades age of 670 Myr, and consistency with all previous and derived astrometric, photometric, and radial velocity results. Velocities from H and He I emission lines confirm the BC period derived from absorption lines, with similar (He I) and higher (H) velocity amplitudes. We ascribe the larger H amplitude to emission from a region each component shadows from the other, depending on the line of sight. * We dedicate this paper to John Stauffer, who died on 2021 January 29, in honor of his many contributions to the field.
  69. Knight, M., Schleicher, D., Farnham, T., 2021, PSJ, 2, 104, Narrowband Observations of Comet 46P/Wirtanen during Its Exceptional Apparition of 2018/19. II. Photometry, Jet Morphology, and Modeling Results
    We report on our extensive photometry and imaging of comet 46P/Wirtanen during its 2018/19 apparition and use these data to constrain the modeling of Wirtanen's activity. Narrowband photometry was obtained in 9 epochs from 2018 October through 2019 March as well as 10 epochs during the 1991, 1997, and 2008 apparitions. The ensemble photometry reveals a typical composition and a secular decrease in activity since 1991. Production rates were roughly symmetric around perihelion for the carbon-bearing species (CN, C3, and C2), but steeper for OH and NH outbound. Our imaging program emphasized CN, whose coma morphology and lightcurve yielded rotation periods reported in a companion paper (Farnham et al. 2021). Here, we compare the gas and dust morphology on the 18 nights for which observations of additional species were obtained. The carbon-bearing species exhibited similar morphology that varied with rotation. OH and NH had broad, hemispheric brightness enhancements in the tailward direction that did not change significantly with rotation, which we attribute to their originating from a substantial icy grain component. We constructed a Monte Carlo model that replicates the shape, motion, and brightness distribution of the CN coma throughout the apparition with a single, self-consistent solution in principal axis rotation. Our model yields a pole having (R.A., decl.) = 319, -5 (pole obliquity of 70) and two large sources (radii of 50 and 40) centered at near-equatorial latitudes and separated in longitude by 160. Applications of the model to explain observed behaviors are discussed.
  70. Sevilla-Noarbe, I., Bechtol, K., Carrasco Kind, M., Carnero Rosell, A., Becker, M., Drlica-Wagner, A., Gruendl, R., Rykoff, E., Sheldon, E., Yanny, B., Alarcon, A., Allam, S., Amon, A., Benoit-Levy, A., Bernstein, G., Bertin, E., Burke, D., Carretero, J., Choi, A., Diehl, H., Everett, S., Flaugher, B., Gaztanaga, E., Gschwend, J., Harrison, I., Hartley, W., Hoyle, B., Jarvis, M., Johnson, M., Kessler, R., Kron, R., Kuropatkin, N., Leistedt, B., Li, T., Menanteau, F., Morganson, E., Ogando, R., Palmese, A., Paz-Chinchon, F., Pieres, A., Pond, C., Rodriguez-Monroy, M., Smith, J., Stringer, K., Troxel, M., Tucker, D., de Vicente, J., Wester, W., Zhang, Y., Abbott, T., Aguena, M., Annis, J., Avila, S., Bhargava, S., Bridle, S., Brooks, D., Brout, D., Castander, F., Cawthon, R., Chang, C., Conselice, C., Costanzi, M., Crocce, M., da Costa, L., Pereira, M., Davis, T., Desai, S., Dietrich, J., Doel, P., Eckert, K., Evrard, A., Ferrero, I., Fosalba, P., Garcia-Bellido, J., Gerdes, D., Giannantonio, T., Gruen, D., Gutierrez, G., Hinton, S., Hollowood, D., Honscheid, K., Huff, E., Huterer, D., James, D., Jeltema, T., Kuehn, K., Lahav, O., Lidman, C., Lima, M., Lin, H., Maia, M., Marshall, J., Martini, P., Melchior, P., Miquel, R., Mohr, J., Morgan, R., Neilsen, E., Plazas, A., Romer, A., Roodman, A., Sanchez, E., Scarpine, V., Schubnell, M., Serrano, S., Smith, M., Suchyta, E., Tarle, G., Thomas, D., To, C., Varga, T., Wechsler, R., Weller, J., Wilkinson, R., DES Collaboration, 2021, ApJS, 254, 24, Dark Energy Survey Year 3 Results: Photometric Data Set for Cosmology
    We describe the Dark Energy Survey (DES) photometric data set assembled from the first three years of science operations to support DES Year 3 cosmologic analyses, and provide usage notes aimed at the broad astrophysics community. Y3 GOLD improves on previous releases from DES, Y1 GOLD, and Data Release 1 (DES DR1), presenting an expanded and curated data set that incorporates algorithmic developments in image detrending and processing, photometric calibration, and object classification. Y3 GOLD comprises nearly 5000 deg2 of grizY imaging in the south Galactic cap, including nearly 390 million objects, with depth reaching a signal-to-noise ratio 10 for extended objects up to iAB 23.0, and top-of-the-atmosphere photometric uniformity <3 mmag. Compared to DR1, photometric residuals with respect to Gaia are reduced by 50%, and per-object chromatic corrections are introduced. Y3 GOLD augments DES DR1 with simultaneous fits to multi-epoch photometry for more robust galactic color measurements and corresponding photometric redshift estimates. Y3 GOLD features improved morphological star-galaxy classification with efficiency >98% and purity >99% for galaxies with 19 < iAB < 22.5. Additionally, it includes per-object quality information, and accompanying maps of the footprint coverage, masked regions, imaging depth, survey conditions, and astrophysical foregrounds that are used to select the cosmologic analysis samples.
  71. Stern, S., Keeney, B., Singer, K., White, O., Hofgartner, J., Grundy, W., The New Horizons Team, 2021, PSJ, 2, 87, Some New Results and Perspectives Regarding the Kuiper Belt Object Arrokoth's Remarkable, Bright Neck
    One of the most striking and curious features of the small Kuiper Belt Object (KBO), Arrokoth, explored by New Horizons is the bright, annular neck it exhibits at the junction between its two lobes. Here we summarize past reported findings regarding the properties of this feature and then report new results regarding its dimensions, reflectivity and color, shape profile, and lack of identifiable craters. We conclude by enumerating possible origin scenarios for this unusual feature. New results include a new estimated measurement of the observed neck area of 8 1.5 km2, a total neck surface area of 32 km2, a 12.5:1 ratio of circumference to height, a normal reflectance histogram of the observed neck, and the fact that no significant (i.e., >2) neck color units were identified, meaning the neck's color is generally spatially uniform at the 1.5 km pixel-1 scale of the best color images. Although several origin hypotheses for the bright material in the neck are briefly discussed, none can be conclusively demonstrated to be the actual origin mechanism at this time; some future tests are identified.
  72. Doux, C., Chang, C., Jain, B., Blazek, J., Camacho, H., Fang, X., Gatti, M., Krause, E., MacCrann, N., Samuroff, S., Secco, L., Troxel, M., Zuntz, J., Aguena, M., Allam, S., Amon, A., Avila, S., Bacon, D., Bertin, E., Brooks, D., Burke, D., Carnero Rosell, A., Carrasco Kind, M., Carretero, J., Choi, A., Costanzi, M., Crocce, M., da Costa, L., Pereira, M., Davis, T., Dietrich, J., Doel, P., Ferrero, I., Ferte, A., Fosalba, P., Garcia-Bellido, J., Gaztanaga, E., Gerdes, D., Gruen, D., Gruendl, R., Gschwend, J., Gutierrez, G., Hartley, W., Hinton, S., Hollowood, D., Huterer, D., James, D., Kuehn, K., Kuropatkin, N., Maia, M., Marshall, J., Menanteau, F., Miquel, R., Morgan, R., Palmese, A., Paz-Chinchon, F., Plazas, A., Roodman, A., Sanchez, E., Schubnell, M., Serrano, S., Sevilla-Noarbe, I., Smith, M., Soares-Santos, M., Suchyta, E., Tarle, G., To, C., Varga, T., Weller, J., Wilkinson, R., DES Collaboration, 2021, MNRAS, 503, 3796, Consistency of cosmic shear analyses in harmonic and real space
    Recent cosmic shear studies have reported discrepancies of up to 1 on the parameter ${S_{8}=\sigma _{8}\sqrt{{\Omega _{\rm m}}/0.3}}$ between the analysis of shear power spectra and two-point correlation functions, derived from the same shear catalogues. It is not a priori clear whether the measured discrepancies are consistent with statistical fluctuations. In this paper, we investigate this issue in the context of the forthcoming analyses from the third year data of the Dark Energy Survey (DES Y3). We analyse DES Y3 mock catalogues from Gaussian simulations with a fast and accurate importance sampling pipeline. We show that the methodology for determining matching scale cuts in harmonic and real space is the key factor that contributes to the scatter between constraints derived from the two statistics. We compare the published scales cuts of the KiDS, Subaru-HSC, and DES surveys, and find that the correlation coefficients of posterior means range from over 80 per cent for our proposed cuts, down to 10 per cent for cuts used in the literature. We then study the interaction between scale cuts and systematic uncertainties arising from multiple sources: non-linear power spectrum, baryonic feedback, intrinsic alignments, uncertainties in the point spread function, and redshift distributions. We find that, given DES Y3 characteristics and proposed cuts, these uncertainties affect the two statistics similarly; the differential biases are below a third of the statistical uncertainty, with the largest biases arising from intrinsic alignment and baryonic feedback. While this work is aimed at DES Y3, the tools developed can be applied to Stage-IV surveys where statistical errors will be much smaller.
  73. Doux, C., Baxter, E., Lemos, P., Chang, C., Alarcon, A., Amon, A., Campos, A., Choi, A., Gatti, M., Gruen, D., Jarvis, M., MacCrann, N., Park, Y., Prat, J., Rau, M., Raveri, M., Samuroff, S., DeRose, J., Hartley, W., Hoyle, B., Troxel, M., Zuntz, J., Abbott, T., Aguena, M., Allam, S., Annis, J., Avila, S., Bacon, D., Bertin, E., Bhargava, S., Brooks, D., Burke, D., Carrasco Kind, M., Carretero, J., Cawthon, R., Costanzi, M., da Costa, L., Pereira, M., Desai, S., Diehl, H., Dietrich, J., Doel, P., Everett, S., Ferrero, I., Fosalba, P., Frieman, J., Garcia-Bellido, J., Gerdes, D., Giannantonio, T., Gruendl, R., Gschwend, J., Gutierrez, G., Hinton, S., Hollowood, D., Honscheid, K., Huff, E., Huterer, D., Jain, B., James, D., Krause, E., Kuehn, K., Kuropatkin, N., Lahav, O., Lidman, C., Lima, M., Maia, M., Menanteau, F., Miquel, R., Morgan, R., Muir, J., Ogando, R., Palmese, A., Paz-Chinchon, F., Plazas, A., Sanchez, E., Scarpine, V., Schubnell, M., Serrano, S., Sevilla-Noarbe, I., Smith, M., Suchyta, E., Swanson, M., Tarle, G., To, C., Tucker, D., Varga, T., Weller, J., Wilkinson, R., 2021, MNRAS, 503, 2688, Dark energy survey internal consistency tests of the joint cosmological probes analysis with posterior predictive distributions
    Beyond CDM, physics or systematic errors may cause subsets of a cosmological data set to appear inconsistent when analysed assuming CDM. We present an application of internal consistency tests to measurements from the Dark Energy Survey Year 1 (DES Y1) joint probes analysis. Our analysis relies on computing the posterior predictive distribution (PPD) for these data under the assumption of CDM. We find that the DES Y1 data have an acceptable goodness of fit to CDM, with a probability of finding a worse fit by random chance of p = 0.046. Using numerical PPD tests, supplemented by graphical checks, we show that most of the data vector appears completely consistent with expectations, although we observe a small tension between large- and small-scale measurements. A small part (roughly 1.5 per cent) of the data vector shows an unusually large departure from expectations; excluding this part of the data has negligible impact on cosmological constraints, but does significantly improve the p-value to 0.10. The methodology developed here will be applied to test the consistency of DES Year 3 joint probes data sets.
  74. Medan, I., Lepine, S., Hartman, Z., 2021, AJ, 161, 234, Bayesian Cross-matching of High Proper-motion Stars in Gaia DR2 and Photometric Metallicities for 1.7 million K and M Dwarfs
    We present a Bayesian method to cross-match 5,827,988 high proper-motion Gaia sources ( > 40 mas yr1) to various photometric surveys: Two Micron All Sky Survey, AllWISE data release from the Wide-field Infrared Explorer (WISE) mission, Galaxy Evolution Explorer, Radial Velocity Experiment, Sloan Digital Sky Survey, and Panoramic Survey Telescope and Rapid Response System (Pan-STARRS). To efficiently associate these objects across catalogs, we develop a technique that compares the multidimensional distribution of all sources in the vicinity of each Gaia star to a reference distribution of random field stars obtained by extracting all sources in a region on the sky displaced 2. This offset preserves the local field stellar density and magnitude distribution, allowing us to characterize the frequency of chance alignments. The resulting catalog with Bayesian probabilities >95 has a marginally higher match rate than current internal Gaia data release 2 (DR2) matches for most catalogs. However, a significant improvement is found with Pan-STARRS, where 99.8 of the sample within the Pan-STARRS footprint is recovered, as compared to a low 20.8 in Gaia DR2. Using these results, we train a Gaussian process regressor to calibrate two photometric metallicity relationships. For dwarfs of 3500 < Teff < 5280 K, we use metallicity values of 4378 stars from the Apache Point Observatory Galactic Evolution Experiment and Hejazi et al. to calibrate the relationship, producing results with a 1 precision of 0.12 dex and few systematic errors. We then indirectly infer the metallicity of 4018 stars with 2850 < Teff < 3500 K, which are wide companions of primaries whose metallicities are estimated with our first regressor, to produce a relationship with a 1 precision of 0.21 dex and significant systematic errors. Additional work is needed to better remove unresolved binaries from this second sample to reduce these systematic errors.
  75. Pang, X., Li, Y., Yu, Z., Tang, S., Dinnbier, F., Kroupa, P., Pasquato, M., Kouwenhoven, M., 2021, ApJ, 912, 162, 3D Morphology of Open Clusters in the Solar Neighborhood with Gaia EDR 3: Its Relation to Cluster Dynamics
    We analyze the 3D morphology and kinematics of 13 open clusters (OCs) located within 500 pc of the Sun, using Gaia EDR 3 and kinematic data from the literature. Members of OCs are identified using the unsupervised machine-learning method StarGO, using five parameters (X, Y, Z, ${\mu }_{\alpha }\cos \delta ,{\mu }_{\delta }$ ). The OC sample covers an age range of 25 Myr to 2.65 Gyr. We correct the asymmetric distance distribution that is due to parallax error using Bayesian inversion. The uncertainty in the corrected distance for a cluster at 500 pc is 3.0-6.3 pc, depending on the intrinsic spatial distribution of its members. We determine the 3D morphology of the OCs in our sample and fit the spatial distribution of stars within the tidal radius in each cluster with an ellipsoid model. The shapes of the OCs are well described with oblate spheroids (NGC 2547, NGC 2516, NGC 2451A, NGC 2451B, and NGC 2232), prolate spheroids (IC 2602, IC 4665, NGC 2422, Blanco 1, and Coma Berenices), or triaxial ellipsoids (IC 2391, NGC 6633, and NGC 6774). The semimajor axis of the fitted ellipsoid is parallel to the Galactic plane for most clusters. Elongated filament-like substructures are detected in three young clusters (NGC 2232, NGC 2547, and NGC 2451B), while tidal-tail-like substructures (tidal tails) are found in older clusters (NGC 2516, NGC 6633, NGC 6774, Blanco 1, and Coma Berenices). Most clusters may be supervirial and expanding. N-body models of rapid gas expulsion with a star formation efficiency of 1/3 are consistent with clusters more massive than 250 M, while clusters less massive than 250 M tend to agree with adiabatic gas expulsion models. Only five OCs (NGC 2422, NGC 6633, NGC 6774, Blanco 1, and Coma Berenices) show clear signs of mass segregation.
  76. Navarro-Meza, S., Aadland, E., Trilling, D., 2021, RNAAS, 5, 111, Asteroid Lightcurves and Detection, Shape, and Size Biases in Large-scale Surveys
    Most asteroids are somewhat elongated and have non-zero lightcurve amplitudes. Such asteroids can be detected in large-scale sky surveys even if their mean magnitudes are fainter than the stated sensitivity limits. We explore the detection of elongated asteroids under a set of idealized but useful approximations. We find that objects up to 1 mag fainter than a survey's sensitivity limit are likely to be detected, and that the effect is most pronounced for asteroids with lightcurve amplitudes 0.1-0.4 mag. This imposes a bias on the derived size and shape distributions of the population that must be properly accounted for.
  77. Levine, S., Zuluaga, C., Person, M., Sickafoose, A., Bosh, A., Collins, M., 2021, AJ, 161, 210, Occultation of a Large Star by the Large Plutino (28978) Ixion on 2020 October 13 UTC
    We observed the occultation of the star Gaia DR2 4056440205544338944 by (28978) Ixion. The event was observed from two Lowell Observatory sites, using the 4.3 m Lowell Discovery Telescope (LDT), near Happy Jack, AZ, USA, and a 0.32 m telescope co-mounted with the Titan Monitoring telescope on Lowell's Mars Hill campus in Flagstaff, AZ. The LDT chord, at 44.86 s, was roughly 30% longer than the longest predicted possible chord. Under the assumption of a spherical body, Ixion's fitted diameter D = 709.6 0.2 km. The LDT light-curve profile was used to place an upper limit on the surface pressure P < 2 bar on any possible atmosphere of Ixion. At the distance of Ixion, the occulted star had a fitted projected diameter of 19.25 0.3 km assuming uniform disk illumination, giving a stellar angular diameter of 0.675 0.010 mas. Using the Gaia EDR3 parallax of 0.565 mas, the stellar radius is ${130}_{-17}^{+20}\,{R}_{\odot }$ <!-- --> . The measured size is consistent with prior spectral classification of this star as a reddened mid-M giant. This is one of only a modest number of M5 III stars to have a directly measured diameter, and is more distant than most.
  78. Scipioni, F., White, O., Cook, J., Bertrand, T., Cruikshank, D., Grundy, W., Beddingfield-Cartwright, C., Binzel, R., Dalle Ore, C., Jennings, D., Moore, J., Olkin, C., Protopapa, S., Reuter, D., Schmitt, B., Singer, K., Spencer, J., Stern, S., Weaver, H., Verbiscer, A., Young, L., 2021, Icar, 359, 114303, Pluto's Sputnik Planitia: Composition of geological units from infrared spectroscopy
    We have compared spectroscopic data of Sputnik Planitia on Pluto, as acquired by New Horizons' Linear Etalon Imaging Spectral Array (LEISA) instrument, to the geomorphology as mapped by White et al. (2017) using visible and panchromatic imaging acquired by the LOng-Range Reconnaissance Imager (LORRI) and the Multi-spectral Visible Imaging Camera (MVIC). We have focused on 13 of the geologic units identified by White et al. (2017), which include the plains and mountain units contained within the Sputnik basin. We divided the map of Sputnik Planitia into 15 provinces, each containing one or more geologic units, and we use LEISA to calculate the average spectra of the units inside the 15 provinces. Hapke-based modeling was then applied to the average spectra of the units to infer their surface composition, and to determine if the composition resulting from the modeling of LEISA spectra reflects the geomorphologic analyses of LORRI data, and if areas classified as being the same geologically, but which are geographically separated, share a similar composition. We investigated the spatial distribution of the most abundant ices on Pluto's surface - CH4, N2, CO, H2O, and a non-ice component presumed to be a macromolecular carbon-rich material, termed a tholin, that imparts a positive spectral slope in the visible spectral region and a negative spectral slope longward of ~1.1 m. Because the exact nature of the non-ice component is still debated and because the negative spectral slope of the available tholins in the near infrared does not perfectly match the Pluto data, for spectral modeling purposes we reference it generically as the negative spectral slope endmember (NSS endmember). We created maps of variations in the integrated band depth (from LEISA data) and areal mass fraction (from the modeling) of the components. The analysis of correlations between the occurrences of the endmembers in the geologic units led to the observation of an anomalous suppression of the strong CH4 absorption bands in units with compositions that are dominated by H2O ice and the NSS endmember. Exploring the mutual variation of the CH4 and N2 integrated band depths with the abundance of crystalline H2O and NSS endmember revealed that the NSS endmember is primarily responsible for the suppression of CH4 absorptions in mountainous units located along the western edge of Sputnik Planitia. Our spectroscopic analyses have provided additional insight into the geological processes that have shaped Sputnik Planitia. A general increase in volatile abundance from the north to the south of Sputnik Planitia is observed. Such an increase first observed and interpreted by Protopapa et al., 2017 and later confirmed by climate modeling (Bertrand et al., 2018) is expressed geomorphologically in the form of preferential deposition of N2 ice in the upland and mountainous regions bordering the plains of southern Sputnik Planitia. Relatively high amounts of pure CH4 are seen at the southern Tenzing Montes, which are a natural site for CH4 deposition owing to their great elevation and the lower insolation they are presently receiving. The NSS endmember correlates the existence of tholins within certain units, mostly those coating the low-latitude mountain ranges that are co-latitudinal with the tholin-covered Cthulhu Macula. The spectral analysis has also revealed compositional differences between the handful of occurrences of northern non-cellular plains and the surrounding cellular plains, all of which are located within the portion of Sputnik Planitia that is presently experiencing net sublimation of volatiles, and which do not therefore exhibit a surface layer of bright, freshly-deposited N2 ice. The compositional differences between the cellular and non-cellular plains here hint at the effectiveness of convection in entraining and trapping tholins within the body of the cellular plains, while preventing the spread of such tholins to abutting non-cellular plains.
  79. Hillier, D., Aadland, E., Massey, P., Morrell, N., 2021, MNRAS, 503, 2726, BAT99-9 - a WC4 Wolf-Rayet star with nitrogen emission: evidence for binary evolution?
    An analysis of the Large Magellanic Cloud (LMC) WC4 star BAT99-9 (HD 32125, FD 4, Brey 7, WS 3) shows that the star still contains photospheric nitrogen. Three N emission features (N V 1238, 1242, N IV 1719, and N IV 3479-3485) are unambiguously identified in the spectrum. CMFGEN models of the star yield an N/C ratio of 0.004 0.002 (by number) and a C/He ratio of $0.15_{-0.05}^{+0.10}$ . Due to the similarity of BAT99-9 to other WC4 stars, and the good fit achieved by CMFGEN to both the classic WC4 spectrum and the N lines, we argue that the N lines are intrinsic to BAT99-9. An examination of a limited set of rotating models for single-star evolution at LMC and Galactic metallicities shows that a model with a Galactic metallicity (z = 0.014) and a progenitor mass of around 50 M can have an N/C ratio similar to, or larger than, what we observe for a significant fraction of its lifetime. However, the LMC models (z = 0.006) are inconsistent with the observations. Both the single and binary BPASS models predict that many WC stars can have an N/C ratio similar to, or larger than, what we observe for a significant fraction of their lifetime. While the binary models cover a wider range of luminosities and provide a somewhat better match to BAT99-9, it is not currently possible to rule out BAT99-9 being formed through single-star evolution, given the uncertainties in mass-loss rates, and the treatment of convection and mixing processes.
  80. Scheirich, P., Pravec, P., Kusnirak, P., Hornoch, K., McMahon, J., Scheeres, D., Capek, D., Pray, D., Kucakova, H., Galad, A., Vrastil, J., Krugly, Y., Moskovitz, N., Avner, L., Skiff, B., McMillan, R., Larsen, J., Brucker, M., Tubbiolo, A., Cooney, W., Gross, J., Terrell, D., Burkhonov, O., Ergashev, K., Ehgamberdiev, S., Fatka, P., Durkee, R., Schunova, E., Inasaridze, R., Ayvazian, V., Kapanadze, G., Gaftonyuk, N., Sanchez, J., Reddy, V., McGraw, L., Kelley, M., Molotov, I., 2021, Icar, 360, 114321, A satellite orbit drift in binary near-Earth asteroids (66391) 1999 KW4 and (88710) 2001 SL9 - Indication of the BYORP effect
    We obtained thorough photometric observations of two binary near-Earth asteroids (66391) Moshup = 1999 KW4 and (88710) 2001 SL9 taken from 2000 to 2019. We modeled the data and derived physical and dynamical properties of the binary systems. For (66391) 1999 KW4, we derived its mutual orbit's pole, semimajor axis and eccentricity that are in agreement with radar-derived values (Ostro et al., 2006. Science, 314, 1276-1280). However, we found that the data are inconsistent with a constant orbital period and we obtained unique solution with a quadratic drift of the mean anomaly of the satellite of -0.65 0.16 deg./yr2 (all quoted uncertainties correspond to 3). This means that the semimajor axis of the mutual orbit of the components of this binary system, determined a = 2.548 0.015 km by Ostro et al. (2006), increases in time with a mean rate of 1.2 0.3 cm/yr.

    For (88710) 2001 SL9, we determined that the mutual orbit has a pole within 10 of (L, B) = (302, -73) (ecliptic coordinates), and is close to circular (eccentricity < 0.07). The data for this system are also inconsistent with a constant orbital period and we obtained two solutions for the quadratic drift of the mean anomaly: 2.8 0.2 and 5.2 0.2 deg./yr2, implying that the semimajor axis of the mutual orbit of the components (estimated a ~ 1.6 km) decreases in time with a mean rate of -2.8 0.2 or -5.1 0.2 cm/yr for the two solutions, respectively.

    The expanding orbit of (66391) 1999 KW4 may be explained by mutual tides interplaying with binary YORP (BYORP) effect (McMahon and Scheeres, 2010a. Icarus 209, 494-509). However, a modeling of the BYORP drift using radar-derived shapes of the binary components predicted a much higher value of the orbital drift than the observed one. It suggests that either the radar-derived shape model of the secondary is inadequate for computing the BYORP effect, or the present theory of BYORP overestimates it. It is possible that the BYORP coefficient has instead an opposite sign than predicted; in that case, the system may be moving into an equilibrium between the BYORP and the tides.

    In the case of (88710) 2001 SL9, the BYORP effect is the only known physical mechanism that can cause the inward drift of its mutual orbit.

    Together with the binary (175706) 1996 FG3 which has a mean anomaly drift consistent with zero, implying a stable equilibrium between the BYORP effect and mutual body tides (Scheirich et al., 2015. Icarus 245, 56-63), we now have three distinct cases of well observed binary asteroid systems with their long-term dynamical models inferred. They indicate a presence of all the three states of the mutual orbit evolution - equilibrium, expanding and contracting - in the population of near-Earth binary asteroids.

  81. Drummond, J., Merline, W., Carry, B., Conrad, A., Tamblyn, P., Enke, B., Christou, J., Dumas, C., Chapman, C., Durda, D., Owen, W., Grundy, W., Reynolds, O., Buckman, M., 2021, Icar, 358, 114275, The orbit of asteroid (317) Roxane's satellite Olympias from Gemini, Keck, VLT and the SOR, and (22) Kalliope's Linus from the SOR
    A definitive orbit is derived for asteroid (317) Roxane's satellite Olympias [S/2009 (317)1] by combining the 2009 discovery images from Gemini North (Merline et al. 2009) with images from Keck and the VLT obtained in 2012, as well as images from its 2016-2017 apparition from the Starfire Optical Range. The orbit is retrograde with respect to the ecliptic but in the same sense as Roxane's spin. Olympias has a period of P=11.94400.0005 days, a semi-major axis of a=2453 km, and an orbital pole at RA=97, Dec=-71, or ecliptic coordinates =245, =-85, close to the south ecliptic pole. This satellite orbital pole is only 3 from Roxane's orbital pole (but in a retrograde sense) and restricts all observations of Olympias from Earth to within 4 of the satellite's orbital plane. By fitting the brightness ratios between Roxane (rotational period of 8.169610.00005 h) and Olympias as a Fourier series, we find a rotational period for Olympias of 8.25870.0001 h, making this an asynchronous wide binary. From the brightness ratios, and with the average infrared modeling diameter found in the literature of 19.160.39 km (error of the mean), we estimate triaxial ellipsoid radii of 14.58.57.2 km for Roxane and 3.62.52.0 km for Olympias. We can then apportion the mass between the two objects and find a density for both (assumed to be the same) of 2.160.18 g/cm3. There are only a few E-type binaries known and this is the first direct determination of E-type density from a binary. We suggest that the system was formed by the Escaping Ejecta Binary (EEB) mechanism of Durda et al. (2004a), probably forming closer together, and then undergoing the complex evolution steps described by Jacobson et al. (2014) involving synchronization, BYORP orbit expansion, loss of tidal lock, and then YORP spinup. Roxane and Olympias may be the only known EEB system to date.

    From the same 2016-2017 apparition the orbit of Linus around asteroid (22) Kalliope is derived from the SOR. This well-observed bright satellite is found to have a circular orbit with a period of P=3.59560.0004 days, in good agreement with the latest elements of Vachier et al. (2012) of P=3.59570.0001 days, and a semi-major axis of a=10996 km, somewhat greater than their a=108211 km for a slightly eccentric orbit (e=0.0070.010). With a diameter for Kalliope of 1616 km (Hanus et al. 2017), we derive a density for Kalliope of 3.720.25 g/cm3 from our one apparition study, the same as Hanus et al. (2017) but greater than the 3.240.16 of Vachier et al. (2012).

  82. Moskovitz, N., Devogele, M., Gustafsson, A., 2021, plde, 28, Unique Capabilities of the 4.3-m Lowell Discovery Telescope (LDT) for Planetary Defense
    The 4.3-m Lowell Discovery Telescope located in Happy Jack, Arizona provides a diverse instrument suite and a demonstrated ability to target faint, fast-moving NEOs. We will give recent examples of high priority NEOs observed with the LDT to serve as informative analogs to the 2021 PDC asteroid. These examples are connected to the following projects:

    (1) An ongoing astrometric program that extends the orbital arc of potentially hazardous asteroids and virtual impactors by achieving detections down to apparent magnitudes of V~25. (2) The Mission Accessible Near-Earth Object Survey (MANOS) collects comprehensive physical characterization data to constrain the compositional, rotational, morphological, and orbital characteristics of the lowest v objects in near-Earth space (e.g. [1-2]). (3) Rotational and orbital characterization of Earth's second known mini-moon, 2020 CD3, when it was at an apparent magnitude of V~23 [3]. (4) High precision rotational lightcurves of NASA DART (Double Asteroid Redirect Test) mission target 65803 Didymos have served as the primary input to mission-critical models that have determined the non-gravitational orbital evolution of Didymos' satellite Dimorphos.

    These examples of science and planetary defense use cases at the LDT have benefited from several unique features of the facility. The Cassegrain focus of LDT is populated by an instrument cube that allows for five instruments to be simultaneously mounted, with switching between instruments taking less than a minute. This versatility allows for the efficient collection of complex data products within a single night, e.g. visible spectra + near-infrared spectra + photometry + astrometry. The accessibility of the full instrument suite on any given night has allowed for flexible scheduling in which rapid response observations can be conducted on very short notice. This flexible schedule coupled with dedicated LDT access for Lowell scientists and partner institutions enables cadence monitoring programs that can be difficult to perform at other facilities. Astrometric, photometric, and spectroscopic tracking of the 2021 PDC asteroid from its discovery in April 2021 through October 2021 would be an ideal use case for the LDT. As a whole these features of LDT make it one of the most capable NEO characterization tools in the world.

  83. Devogele, M., Moskovitz, N., 2021, plde, 35, Polarimetry as a tool for physical characterization of potentially hazardous asteroids
    The degree of linear polarization of sunlight scattered by an asteroid contains valuable information for rapid characterization of the surface properties of Near-Earth objects. In the case of atmosphereless bodies the state of linear polarization varies as a function of the phase angle () and is described using the so-called Pr parameter (see left panel of Fig. 1).

    The properties of the phase-polarization curve (see right panel of Fig. 1) of an asteroid are mostly defined by its albedo (pV). Numerous calibrations between polarization and pV have been proposed for main-belt asteroids [1, 2]. However, main-belt asteroids rarely exceed phase angle > 30 while nearEarth object can be observed at phase angle as large as 100. These observation at higher phase angle to the NEOs allows for deeper characterization of the observed object, but there is currently a lack of observations of NEOs in polarimetry to accurately calibrate the albedo-polarization relationship at high phase angles.

    In this presentation I will discuss the current state of NEO observations in polarimetry and how polarimetry could be used to obtain reliable information on the geometric albedo of NEOs. With a proper calibration of the polarization-albedo relation, one could reduce the uncertainty on a newly discovered object by a factor of 10 with one single polarimetric observation obtained at a phase angle > 40.

    As an example, according to the current state of the albedo-polarization calibration at large phase angle and considering an H = 22.3 mag NEO (similar to 2021 PDC). Without any albedo information size estimation would be from 35 meters to 700 meters. However, one measurements of Pr at a phase angle of 40 would allow to obtain reliable information about the diameter of the object as follow:

    D < 80 meters if Pr < 1% 80 < D < 120 meters is 1 < Pr < 4% D > 120 meters if Pr > 4%

    while the typical uncertainties on Pr are generally lower than 0.1%.

    I will also discuss the specific case of 2021 PDC and what a few observations in polarimetry obtained with the FORS2 polarimeter on the Very Large Telescope could provide. FORS2 on the VLT is currently the only instrument capable of obtained reliable polarimetric observations on the V=21.5 mag objects.

  84. Medeiros, H., de Leon, J., Licandro, J., Serra-Ricart, M., Popescu, M., Pinilla-Alonso, N., Devogele, M., 2021, plde, 112, Photometric analysis of NEOs in Arecibo supporting from NEOROCKS project
    The near-Earth objects (NEOs) include asteroids and comets that are in orbit close to the Earth. In some cases, they can even cross the Earth's orbit. The NEOs are one interesting population of small bodies for several reasons. There are objects that constitute a potential impact hazard. Because of the close approach to the Earth, they are more accessible to spacecraft, to detailed observations from ground-based facilities, and also allow us to study asteroids two to three orders of magnitude smaller than those observable on the main belt asteroids.

    The NEO Rapid Observation, Characterization and Key Simulations (NEOROCKS) project is funded (2020-2022) through the H2020 European Commission programme to improve knowledge on NEOs by connecting expertise in performing small body astronomical observations and the related modelling needed to derive their dynamical and physical properties. The IAC, and in particular members of the Solar System Group, participates in the NEOROCKS project and is currently leading one specific task to do observations of NEOs in support of the Arecibo Planetary Radar. These supporting observations include times-series photometry (lightcurves), visible to near-infrared spectroscopy, and astrometry (orbit determination). To do that we are using the telescopes from Roque de Los Muchachos and Teide Observatories, located at La Palma and Tenerife islands, respectively (Spain).

    Considering the current situation of the Arecibo telescope, in this work we focus in the obtention and analysis of lightcurve data for those targets that were observed in the past by the Arecibo telescope and have good signal-to-noise (SNR) radar data. Our photometry is intented to support and complement radar data and to help deriving physical properties like rotational period and shape. We have observed 8 asteroids to obtain their lightcurve. Figure 1 shows one example of lightcurve from the asteroid (52768) 1998 OR2, that we are observing in studing with Arecibo team. These are preliminary results.

  85. Naidu, S., Chesley, S., Farnocchia, D., Pravec, P., Scheirich, P., Moskovitz, N., Thomas, C., Rivkin, A., 2021, plde, 198, Estimation of orbital parameters of Dimorphos from lightcurve mutual events
    Binary near-Earth asteroid (65803) Didymos is the target of NASA's Double Asteroid Redirection Test (DART) space mission and ESA's Hera mission. The DART mission is a planetary defense experiment with a planned launch in July 2021. The mission consists of a spacecraft that would impact Dimorphos, the satellite of Didymos, in October 2022 and change its orbital period around Didymos by > 73 seconds. The change would be measured by ground-based photometric and radar observations over several months after the impact. The Hera mission will characterize Didymos and its satellite starting in 2026.

    We estimated the pre-impact orbit of Dimorphos using the times of occultations and eclipses observed in lightcurves from 2003, 2015, 2017, and 2019 using a weighted least-squares approach. The data is consistent with three separate solutions corresponding to different values of Binary YORP, a radiative effect that could cause a drift in mean motion. The three nominal solutions are separated by about 50 in orbital phase during the planned DART mission impact in October 2022 and their formal 3-sigma uncertainties are about 30 for each. Data from the ongoing 2021 observation campaign (Thomas et al., this conference) will be used to refine the estimates further. Observations in December 2020 to March 2021 are expected to eliminate two of the three solutions and reduce the formal 3-sigma uncertainty in the orbital phase during the DART mission to be < 10.

  86. Thomas, C., Rivkin, A., Moskovitz, N., Pravec, P., Scheirich, P., de Leon, J., Mazzotta Epifani, E., Birlan, M., Carry, B., Chesley, S., Dotto, E., Farnham, T., Fohring, D., Granvik, M., Howell, E., Ieva, S., Knight, M., Lazzarin, M., La Forgia, F., May, E., Migliorini, A., Naidu, S., Popescu, M., 2021, plde, 201, Constraining the Orbital Parameters of the Didymos-Dimorphos System: Lightcurve Observations in Preparation for AIDA/DART
    The binary near-Earth asteroid (65803) Didymos is the target for the Asteroid Impact and Deflection Assessment (AIDA) mission, which is a concept with two primary spacecraft: NASA's DART (Double Asteroid Redirection Test) impactor and ESA's Hera orbiter (Cheng et al. 2018; Michel et al. 2018). DART is NASA's first planetary defense mission and will be the first demonstration of asteroid deflection by a kinetic impactor. The DART spacecraft is designed to impact Dimorphos, the secondary in the Didymos system, and modify its orbit through momentum transfer. DART will launch in mid-2021 and is scheduled to impact in fall 2022. The DART spacecraft will carry ASI's LICIACube (Light Italian Cubesat for Imaging of Asteroids, Dotto et al. 2021) to observe the DART impact event and the resulting impact ejecta.

    A key scientific goal of the DART and Hera missions is to measure and characterize the deflection caused by the DART impact. The impact will change the satellite orbit period, which will be measured by ground-based facilities in the post-impact period. In order to correctly interpret the data from the impact epoch, we need to understand the baseline, unperturbed dynamics of the system. The DART/Hera Observations Working Group is tasked with characterizing the Didymos-Dimorphos system properties with sufficient accuracy to measure the change in the binary orbital period to within 7.3 seconds. This measurement is a small, but observable fraction of the current orbital period of the satellite (Porb=11.92 hours). The observed period change is a critical input to the calculation of the momentum transfer enhancement parameter, "Beta" (). We are obtaining lightcurve observations during the current apparition (December 2020 to March 2021) to further characterize the system. Combining these observations with past data will provide us with the opportunity to establish the state of the system before impact to a high level of precision.

    We have two goals for our 2020-2021 Didymos observing effort: (1) Measure the amount of Binary YORP (BYORP) torque occurring in the system and (2) Establish whether or not the secondary is in synchronous rotation. An international group of observers affiliated with DART and Hera obtained time at 11 different facilities over the time period from December 2020 to March 2021. We will discuss our state of knowledge from previous observations (through 2019) and selected results from our 2020-2021 observations.

  87. Gustafsson, A., MacLennan, E., Moskovitz, N., Trilling, D., 2021, plde, 225, Regolith Constraints for Near-Earth Asteroid 2021 PDC
    Asteroid regolith grain size is diagnostic of several poorly understood geophysical processes and has direct planetary defense applications. For example, it is well known that a rubble pile size barrier exists around asteroid diameters of 200-300m [1]. Small bodies with diameters below this size can rotate quickly with periods less than 2.2 hours, however, large objects generally do not rotate that quickly. This rotation rate barrier is attributed to a transition from large, gravity dominated objects to small, strength dominated ones [2]. Yet, it is unclear how this transition affects surface properties. Thermal (e.g. [3]) and in-situ (e.g. [4]) data suggest that object diameter influences surface grain size -- large bodies have fine surface grains and small bodies have coarse grains. However, rotation rate [5] and size-independent processes like thermal fracturing [6] and tidal interactions with planets [7] may be equally important in affecting surface properties.

    For planetary defense risk assessment, quantifying the nature of the regolith and near-subsurface material of potential impactors can inform early predictions of the momentum transfer from a kinetic impactor [8].

    We will present an overview of Hapke radiative transfer and thermal modeling tools that we have adopted to constrain near-Earth asteroid surface regolith grain sizes. These can be used for hazard assessment of close-approaching near-Earth asteroids and for better understanding the surface characteristics of objects above and below the size barrier. In this work, we focus on the application of these tools with respect to characterizing surface properties of near-Earth asteroid 2021 PDC. Specifically, we will apply our modeling and data analysis tools to synthetic observations of 2021 PDC from the 4.3-m Lowell Discovery Telescope and NASA's Infrared Telescope Facility.

    We have adopted a probabilistic implementation of Hapke radiative transfer modeling [9] on unresolved asteroids to constrain surface regolith grain size. By collecting high signal-to-noise (S/N > 50) photometry of objects going through opposition, we can fit the data with our Hapke photometric model to constrain regolith grain size.

    Using multi-epoch thermal measurements that are collected at opposition and across a wide range of solar phase angles, we can estimate thermal inertia and, from it, the surface grain size following methods by [10] and [11]. A reasonable estimate for the rotation period is needed for this approach, but shape and spin properties are treated as unknown free-parameters when fitting to the observations. Detailed knowledge of the lightcurve period and amplitude can be used to constrain the range of these free parameters.

    We will compare the thermal-derived grain size estimates to those constrained by our implementation of the Hapke photometric model to better understand the surface properties of the body across multiple wavelengths. For example, a bare-rock thermal inertia value and coarse Hapke grain size may suggest an intact monolith, as opposed to a rubble pile comprised of various-sized boulders and fine surface grains.

  88. Huang, H., Eifler, T., Mandelbaum, R., Bernstein, G., Chen, A., Choi, A., Garcia-Bellido, J., Huterer, D., Krause, E., Rozo, E., Singh, S., Bridle, S., DeRose, J., Elvin-Poole, J., Fang, X., Friedrich, O., Gatti, M., Gaztanaga, E., Gruen, D., Hartley, W., Hoyle, B., Jarvis, M., MacCrann, N., Miranda, V., Rau, M., Prat, J., Sanchez, C., Samuroff, S., Troxel, M., Zuntz, J., Abbott, T., Aguena, M., Annis, J., Avila, S., Becker, M., Bertin, E., Brooks, D., Burke, D., Carnero Rosell, A., Carrasco Kind, M., Carretero, J., Castander, F., da Costa, L., De Vicente, J., Dietrich, J., Doel, P., Everett, S., Flaugher, B., Fosalba, P., Frieman, J., Gruendl, R., Gutierrez, G., Hinton, S., Honscheid, K., James, D., Kuehn, K., Lahav, O., Lima, M., Maia, M., Marshall, J., Menanteau, F., Miquel, R., Paz-Chinchon, F., Malagon, A., Romer, K., Roodman, A., Sanchez, E., Scarpine, V., Serrano, S., Sevilla, I., Smith, M., Soares-Santos, M., Suchyta, E., Swanson, M., Tarle, G., Thomas, D., Weller, J., DES Collaboration, 2021, MNRAS, 502, 6010, Dark energy survey year 1 results: Constraining baryonic physics in the Universe
    Measurements of large-scale structure are interpreted using theoretical predictions for the matter distribution, including potential impacts of baryonic physics. We constrain the feedback strength of baryons jointly with cosmology using weak lensing and galaxy clustering observables (3 2pt) of Dark Energy Survey (DES) Year 1 data in combination with external information from baryon acoustic oscillations (BAO) and Planck cosmic microwave background polarization. Our baryon modelling is informed by a set of hydrodynamical simulations that span a variety of baryon scenarios; we span this space via a Principal Component (PC) analysis of the summary statistics extracted from these simulations. We show that at the level of DES Y1 constraining power, one PC is sufficient to describe the variation of baryonic effects in the observables, and the first PC amplitude (Q1) generally reflects the strength of baryon feedback. With the upper limit of Q1 prior being bound by the Illustris feedback scenarios, we reach $\sim 20{{\ \rm per\ cent}}$ improvement in the constraint of $S_8=\sigma _8(\Omega _{\rm m}/0.3)^{0.5}=0.788^{+0.018}_{-0.021}$ compared to the original DES 3 2pt analysis. This gain is driven by the inclusion of small-scale cosmic shear information down to 2.5 arcmin, which was excluded in previous DES analyses that did not model baryonic physics. We obtain $S_8=0.781^{+0.014}_{-0.015}$ for the combined DES Y1+Planck EE+BAO analysis with a non-informative Q1 prior. In terms of the baryon constraints, we measure $Q_1=1.14^{+2.20}_{-2.80}$ for DES Y1 only and $Q_1=1.42^{+1.63}_{-1.48}$ for DESY1+Planck EE+BAO, allowing us to exclude one of the most extreme AGN feedback hydrodynamical scenario at more than 2.
  89. Mucesh, S., Hartley, W., Palmese, A., Lahav, O., Whiteway, L., Bluck, A., Alarcon, A., Amon, A., Bechtol, K., Bernstein, G., Carnero Rosell, A., Carrasco Kind, M., Choi, A., Eckert, K., Everett, S., Gruen, D., Gruendl, R., Harrison, I., Huff, E., Kuropatkin, N., Sevilla-Noarbe, I., Sheldon, E., Yanny, B., Aguena, M., Allam, S., Bacon, D., Bertin, E., Bhargava, S., Brooks, D., Carretero, J., Castander, F., Conselice, C., Costanzi, M., Crocce, M., da Costa, L., Pereira, M., De Vicente, J., Desai, S., Diehl, H., Drlica-Wagner, A., Evrard, A., Ferrero, I., Flaugher, B., Fosalba, P., Frieman, J., Garcia-Bellido, J., Gaztanaga, E., Gerdes, D., Gschwend, J., Gutierrez, G., Hinton, S., Hollowood, D., Honscheid, K., James, D., Kuehn, K., Lima, M., Lin, H., Maia, M., Melchior, P., Menanteau, F., Miquel, R., Morgan, R., Paz-Chinchon, F., Plazas, A., Sanchez, E., Scarpine, V., Schubnell, M., Serrano, S., Smith, M., Suchyta, E., Tarle, G., Thomas, D., To, C., Varga, T., Wilkinson, R., DES Collaboration, 2021, MNRAS, 502, 2770, A machine learning approach to galaxy properties: joint redshift-stellar mass probability distributions with Random Forest
    We demonstrate that highly accurate joint redshift-stellar mass probability distribution functions (PDFs) can be obtained using the Random Forest (RF) machine learning (ML) algorithm, even with few photometric bands available. As an example, we use the Dark Energy Survey (DES), combined with the COSMOS2015 catalogue for redshifts and stellar masses. We build two ML models: one containing deep photometry in the griz bands, and the second reflecting the photometric scatter present in the main DES survey, with carefully constructed representative training data in each case. We validate our joint PDFs for 10 699 test galaxies by utilizing the copula probability integral transform and the Kendall distribution function, and their univariate counterparts to validate the marginals. Benchmarked against a basic set-up of the template-fitting code BAGPIPES, our ML-based method outperforms template fitting on all of our predefined performance metrics. In addition to accuracy, the RF is extremely fast, able to compute joint PDFs for a million galaxies in just under 6 min with consumer computer hardware. Such speed enables PDFs to be derived in real time within analysis codes, solving potential storage issues. As part of this work we have developed GALPRO1, a highly intuitive and efficient PYTHON package to rapidly generate multivariate PDFs on-the-fly. GALPRO is documented and available for researchers to use in their cosmology and galaxy evolution studies.
  90. Jenniskens, P., Gabadirwe, M., Yin, Q., Proyer, A., Moses, O., Kohout, T., Franchi, F., Gibson, R., Kowalski, R., Christensen, E., Gibbs, A., Heinze, A., Denneau, L., Farnocchia, D., Chodas, P., Gray, W., Micheli, M., Moskovitz, N., Onken, C., Wolf, C., Devillepoix, H., Ye, Q., Robertson, D., Brown, P., Lyytinen, E., Moilanen, J., Albers, J., Cooper, T., Assink, J., Evers, L., Lahtinen, P., Seitshiro, L., Laubenstein, M., Wantlo, N., Moleje, P., Maritinkole, J., Suhonen, H., Zolensky, M., Ashwal, L., Hiroi, T., Sears, D., Sehlke, A., Maturilli, A., Sanborn, M., Huyskens, M., Dey, S., Ziegler, K., Busemann, H., Riebe, M., Meier, M., Welten, K., Caffee, M., Zhou, Q., Li, Q., Li, X., Liu, Y., Tang, G., McLain, H., Dworkin, J., Glavin, D., Schmitt-Kopplin, P., Sabbah, H., Joblin, C., Granvik, M., Mosarwa, B., Botepe, K., 2021, M&PS, 56, 844, The impact and recovery of asteroid 2018 LA
    The June 2, 2018 impact of asteroid 2018 LA over Botswana is only the second asteroid detected in space prior to impacting over land. Here, we report on the successful recovery of meteorites. Additional astrometric data refine the approach orbit and define the spin period and shape of the asteroid. Video observations of the fireball constrain the asteroid's position in its orbit and were used to triangulate the location of the fireball's main flare over the Central Kalahari Game Reserve. Twenty three meteorites were recovered. A consortium study of eight of these classifies Motopi Pan as an HED polymict breccia derived from howardite, cumulate and basaltic eucrite, and diogenite lithologies. Before impact, 2018 LA was a solid rock of ~156 cm diameter with high bulk density ~2.85 g cm3, a relatively low albedo pV ~ 0.25, no significant opposition effect on the asteroid brightness, and an impact kinetic energy of ~0.2 kt. The orbit of 2018 LA is consistent with an origin at Vesta (or its Vestoids) and delivery into an Earth impacting orbit via the 6 resonance. The impact that ejected 2018 LA in an orbit toward Earth occurred 22.8 3.8 Ma ago. Zircons record a concordant U Pb age of 4563 11 Ma and a consistent 207Pb/206Pb age of 4563 6 Ma. A much younger Pb Pb phosphate resetting age of 4234 41 Ma was found. From this impact chronology, we discuss what is the possible source crater of Motopi Pan and the age of Vesta's Veneneia impact basin.
  91. Nofi, L., Johns-Krull, C., Lopez-Valdivia, R., Biddle, L., Carvalho, A., Huber, D., Jaffe, D., Llama, J., Mace, G., Prato, L., Skiff, B., Sokal, K., Sullivan, K., Tayar, J., 2021, ApJ, 911, 138, Projected Rotational Velocities and Fundamental Properties of Low-mass Pre-main-sequence Stars in the Taurus-Auriga Star-forming Region
    The projected stellar rotational velocity ( $v\sin i$ ) is critical for our understanding of processes related to the evolution of angular momentum in pre-main-sequence stars. We present $v\sin i$ measurements of high-resolution infrared and optical spectroscopy for 70 pre-main-sequence stars in the Taurus-Auriga star-forming region, in addition to effective temperatures measured from line-depth ratios, as well as stellar rotation periods determined from optical photometry. From the literature, we identified the stars in our sample that show evidence of residing in circumstellar disks or multiple systems. The comparison of infrared $v\sin i$ measurements calculated using two techniques shows a residual scatter of 1.8 km s-1, defining a typical error floor for the $v\sin i$ of pre-main-sequence stars from infrared spectra. A comparison of the $v\sin i$ distributions of stars with and without companions shows that binaries/multiples typically have a higher measured $v\sin i$ , which may be caused by contamination by companion lines, shorter disk lifetimes in binary systems, or tidal interactions in hierarchical triples. A comparison of optical and infrared $v\sin i$ values shows no significant difference regardless of whether the star has a disk or not, indicating that CO contamination from the disk does not impact $v\sin i$ measurements above the typical 1.8 km s-1 error floor of our measurements. Finally, we observe a lack of a correlation between the $v\sin i$ , presence of a disk, and H-R diagram position, which indicates a complex interplay between stellar rotation and evolution of pre-main-sequence stars.
  92. Li, T., Koposov, S., Erkal, D., Ji, A., Shipp, N., Pace, A., Hilmi, T., Kuehn, K., Lewis, G., Mackey, D., Simpson, J., Wan, Z., Zucker, D., Bland-Hawthorn, J., Cullinane, L., Da Costa, G., Drlica-Wagner, A., Hattori, K., Martell, S., Sharma, S., S5 Collaboration, 2021, ApJ, 911, 149, Broken into Pieces: ATLAS and Aliqa Uma as One Single Stream
    We present the first spectroscopic measurements of the ATLAS and Aliqa Uma streams from the Southern Stellar Stream Spectroscopic Survey (S5), in combination with the photometric data from the Dark Energy Survey and astrometric data from Gaia. From the coherence of spectroscopic members in radial velocity and proper motion, we find that these two systems are extremely likely to be one stream with discontinuity in morphology and density on the sky (the "kink" feature). We refer to this entire stream as the ATLAS-Aliqa Uma stream, or the AAU stream. We perform a comprehensive exploration of the effect of baryonic substructures and find that only an encounter with the Sagittarius dwarf 0.5 Gyr ago can create a feature similar to the observed "kink." In addition, we also identify two gaps in the ATLAS component associated with the broadening in the stream width (the "broadening" feature). These gaps have likely been created by small mass perturbers, such as dark matter halos, as the AAU stream is the most distant cold stream known with severe variations in both the stream surface density and the stream track on the sky. With the stream track, stream distance, and kinematic information, we determine the orbit of the AAU stream and find that it has been affected by the Large Magellanic Cloud, resulting in a misalignment between the proper motion and stream track. Together with the Orphan-Chenab Stream, AAU is the second stream pair that has been found to be a single stream separated into two segments by external perturbation.
  93. Stringer, K., Drlica-Wagner, A., Macri, L., Martinez-Vazquez, C., Vivas, A., Ferguson, P., Pace, A., Walker, A., Neilsen, E., Tavangar, K., Wester, W., Abbott, T., Aguena, M., Allam, S., Bacon, D., Bechtol, K., Bertin, E., Brooks, D., Burke, D., Carnero Rosell, A., Carrasco Kind, M., Carretero, J., Costanzi, M., Crocce, M., da Costa, L., Pereira, M., De Vicente, J., Desai, S., Diehl, H., Doel, P., Ferrero, I., Garcia-Bellido, J., Gaztanaga, E., Gerdes, D., Gruen, D., Gruendl, R., Gschwend, J., Gutierrez, G., Hinton, S., Hollowood, D., Honscheid, K., Hoyle, B., James, D., Kuehn, K., Kuropatkin, N., Li, T., Maia, M., Marshall, J., Menanteau, F., Miquel, R., Morgan, R., Ogando, R., Palmese, A., Paz-Chinchon, F., Plazas, A., Roodman, A., Sanchez, E., Schubnell, M., Serrano, S., Sevilla-Noarbe, I., Smith, M., Soares-Santos, M., Suchyta, E., Tarle, G., Thomas, D., To, C., Varga, T., Wilkinson, R., Zhang, Y., DES Collaboration, 2021, ApJ, 911, 109, Identifying RR Lyrae Variable Stars in Six Years of the Dark Energy Survey
    We present a search for RR Lyrae stars using the full six-year data set from the Dark Energy Survey covering 5000 deg2 of the southern sky. Using a multistage multivariate classification and light-curve template-fitting scheme, we identify RR Lyrae candidates with a median of 35 observations per candidate. We detect 6971 RR Lyrae candidates out to 335 kpc, and we estimate that our sample is >70% complete at 150 kpc. We find excellent agreement with other wide-area RR Lyrae catalogs and RR Lyrae studies targeting the Magellanic Clouds and other Milky Way satellite galaxies. We fit the smooth stellar halo density profile using a broken-power-law model with fixed halo flattening (q = 0.7), and we find strong evidence for a break at ${R}_{0}={32.1}_{-0.9}^{+1.1}\,\mathrm{kpc}$ with an inner slope of ${n}_{1}=-{2.54}_{-0.09}^{+0.09}$ and an outer slope of ${n}_{2}=-{5.42}_{-0.14}^{+0.13}$ . We use our catalog to perform a search for Milky Way satellite galaxies with large sizes and low luminosities. Using a set of simulated satellite galaxies, we find that our RR Lyrae-based search is more sensitive than those using resolved stellar populations in the regime of large (rh 500 pc), low-surface-brightness dwarf galaxies. A blind search for large, diffuse satellites yields three candidate substructures. The first can be confidently associated with the dwarf galaxy Eridanus II. The second has a distance and proper motion similar to the ultrafaint dwarf galaxy Tucana II but is separated by 5 deg. The third is close in projection to the globular cluster NGC 1851 but is 10 kpc more distant and appears to differ in proper motion.
  94. To, C., Krause, E., Rozo, E., Wu, H., Gruen, D., Wechsler, R., Eifler, T., Rykoff, E., Costanzi, M., Becker, M., Bernstein, G., Blazek, J., Bocquet, S., Bridle, S., Cawthon, R., Choi, A., Crocce, M., Davis, C., DeRose, J., Drlica-Wagner, A., Elvin-Poole, J., Fang, X., Farahi, A., Friedrich, O., Gatti, M., Gaztanaga, E., Giannantonio, T., Hartley, W., Hoyle, B., Jarvis, M., MacCrann, N., McClintock, T., Miranda, V., Pereira, M., Park, Y., Porredon, A., Prat, J., Rau, M., Ross, A., Samuroff, S., Sanchez, C., Sevilla-Noarbe, I., Sheldon, E., Troxel, M., Varga, T., Vielzeuf, P., Zhang, Y., Zuntz, J., Abbott, T., Aguena, M., Amon, A., Annis, J., Avila, S., Bertin, E., Bhargava, S., Brooks, D., Burke, D., Carnero Rosell, A., Carrasco Kind, M., Carretero, J., Chang, C., Conselice, C., da Costa, L., Davis, T., Desai, S., Diehl, H., Dietrich, J., Everett, S., Evrard, A., Ferrero, I., Flaugher, B., Fosalba, P., Frieman, J., Garcia-Bellido, J., Gruendl, R., Gutierrez, G., Hinton, S., Hollowood, D., Honscheid, K., Huterer, D., James, D., Jeltema, T., Kron, R., Kuehn, K., Kuropatkin, N., Lima, M., Maia, M., Marshall, J., Menanteau, F., Miquel, R., Morgan, R., Muir, J., Myles, J., Palmese, A., Paz-Chinchon, F., Plazas, A., Romer, A., Roodman, A., Sanchez, E., Santiago, B., Scarpine, V., Serrano, S., Smith, M., Suchyta, E., Swanson, M., Tarle, G., Thomas, D., Tucker, D., Weller, J., Wester, W., Wilkinson, R., DES Collaboration, 2021, PhRvL, 126, 141301, Dark Energy Survey Year 1 Results: Cosmological Constraints from Cluster Abundances, Weak Lensing, and Galaxy Correlations
    We present the first joint analysis of cluster abundances and auto or cross-correlations of three cosmic tracer fields: galaxy density, weak gravitational lensing shear, and cluster density split by optical richness. From a joint analysis (4 2 pt +N ) of cluster abundances, three cluster cross-correlations, and the auto correlations of the galaxy density measured from the first year data of the Dark Energy Survey, we obtain m=0.30 5-0.038+0.055 and 8=0.78 3-0.054+0.064. This result is consistent with constraints from the DES-Y1 galaxy clustering and weak lensing two-point correlation functions for the flat CDM model. Consequently, we combine cluster abundances and all two-point correlations from across all three cosmic tracer fields (6 2 pt +N ) and find improved constraints on cosmological parameters as well as on the cluster observable-mass scaling relation. This analysis is an important advance in both optical cluster cosmology and multiprobe analyses of upcoming wide imaging surveys.
  95. Hsieh, H., Ishiguro, M., Knight, M., Moskovitz, N., Sheppard, S., Trujillo, C., 2021, PSJ, 2, 62, The Reactivation of Main-belt Comet 259P/Garradd (P/2008 R1)
    We present observations of main-belt comet (MBC) 259P/Garradd from 4 months prior to its 2017 perihelion passage to 5 months after perihelion using the Gemini North and South telescopes. The object was confirmed to be active during this period, placing it among seven MBCs confirmed to have recurrent activity. We find an average net pre-perihelion dust production rate for 259P in 2017 of ${\dot{M}}_{d}=(4.6\pm 0.2)$ kg s-1 (assuming grain densities of = 2500 kg m-3 and a mean effective particle size of ${\bar{a}}_{d}=2$ mm) and a best-fit start date of detectable activity of 2017 April 22 1, when the object was at a heliocentric distance of rh = 1.96 0.03 au and a true anomaly of = 3139 04. We estimate the effective active fraction of 259P's surface area to be from fact 7 10-3 to fact 6 10-2 (corresponding to effective active areas of Aact 8 103 m2 to Aact 7 104 m2) at the start of its 2017 active period. A comparison of estimated total dust masses measured for 259P in 2008 and 2017 shows no evidence of changes in activity strength between the two active apparitions. The heliocentric distance of 259P's activity onset point is much smaller than those of other MBCs, suggesting that its ice reservoirs may be located at greater depths than on MBCs farther from the Sun, increasing the time needed for a solar-irradiation-driven thermal wave to reach subsurface ice. We suggest that deeper ice on 259P could be a result of more rapid ice depletion caused by the object's closer proximity to the Sun compared to other MBCs.
  96. Hunter, D., Elmegreen, B., Archer, H., Simpson, C., Cigan, P., 2021, AJ, 161, 175, A Search for Correlations between Turbulence and Star Formation in LITTLE THINGS Dwarf Irregular Galaxies
    Turbulence has the potential for creating gas density enhancements that initiate cloud and star formation (SF), and it can be generated locally by SF. To study the connection between turbulence and SF, we looked for relationships between SF traced by FUV images, and gas turbulence traced by kinetic energy density (KED) and velocity dispersion (vdisp) in the LITTLE THINGS sample of nearby dIrr galaxies. We performed 2D cross-correlations between FUV and KED images, measured cross-correlations in annuli to produce correlation coefficients as a function of radius, and determined the cumulative distribution function of the cross-correlation value. We also plotted on a pixel-by-pixel basis the locally excess KED, vdisp, and H I mass surface density, HI, as determined from the respective values with the radial profiles subtracted, versus the excess SF rate density SFR, for all regions with positive excess SFR. We found that SFR and KED are poorly correlated. The excess KED associated with SF implies a 0.5% efficiency for supernova energy to pump local H I turbulence on the scale of the resolution here, which is a factor of 2 too small for all of the turbulence on a galactic scale. The excess vdisp in SF regions is also small, only 0.37 km s-1. The local excess in HI corresponding to an excess in SFR is consistent with a H I consumption time of 1.6 Gyr in the inner parts of the galaxies. The similarity between this timescale and the consumption time for CO implies that CO-dark molecular gas has comparable mass to H I in the inner disks.
  97. Watkins, J., Salzer, J., Van Sistine, A., Hayslip, A., Hoar, E., Rampalli, R., 2021, ApJS, 253, 39, The H Dots Survey. IV. A Fourth List of Faint Emission-line Objects
    We present the fourth catalog of serendipitously discovered compact extragalactic emission-line sourcesH Dots. A total of 454 newly discovered objects are included in the current survey list. These objects have been detected in searches of moderately deep narrowband images acquired for the ALFALFA H project. The catalog of H Dots presented in the current paper was derived from searches carried out using ALFALFA H images obtained with the KPNO 2.1 m telescope. This results in a substantially deeper sample of Dots compared to our previous lists, which were all discovered in images taken with the WIYN 0.9 m telescope. The median R-band magnitude of the current catalog is 21.59, more than 1.6 mag fainter than the median for the 0.9 m sample (a factor of 4.4 fainter). Likewise, the median emission-line flux of the detected sources is a factor of 4.3 fainter. The line flux completeness limit of the current sample is ~3 10-16 erg s-1 cm-2. We present accurate coordinates, apparent magnitudes, and narrowband line fluxes for each object in the sample. Unlike our previous lists of H Dots, the current sample does not include follow-up spectroscopy.
  98. Napier, K., Gerdes, D., Lin, H., Hamilton, S., Bernstein, G., Bernardinelli, P., Abbott, T., Aguena, M., Annis, J., Avila, S., Bacon, D., Bertin, E., Brooks, D., Burke, D., Carnero Rosell, A., Carrasco Kind, M., Carretero, J., Costanzi, M., da Costa, L., De Vicente, J., Diehl, H., Doel, P., Everett, S., Ferrero, I., Fosalba, P., Garcia-Bellido, J., Gruen, D., Gruendl, R., Gutierrez, G., Hollowood, D., Honscheid, K., Hoyle, B., James, D., Kent, S., Kuehn, K., Kuropatkin, N., Maia, M., Menanteau, F., Miquel, R., Morgan, R., Palmese, A., Paz-Chinchon, F., Plazas, A., Sanchez, E., Scarpine, V., Serrano, S., Sevilla-Noarbe, I., Smith, M., Suchyta, E., Swanson, M., To, C., Walker, A., Wilkinson, R., DES Collaboration, 2021, PSJ, 2, 59, No Evidence for Orbital Clustering in the Extreme Trans-Neptunian Objects
    The apparent clustering in longitude of perihelion and ascending node of extreme trans-Neptunian objects (ETNOs) has been attributed to the gravitational effects of an unseen 5-10 Earth-mass planet in the outer solar system. To investigate how selection bias may contribute to this clustering, we consider 14 ETNOs discovered by the Dark Energy Survey, the Outer Solar System Origins Survey, and the survey of Sheppard and Trujillo. Using each survey's published pointing history, depth, and TNO tracking selections, we calculate the joint probability that these objects are consistent with an underlying parent population with uniform distributions in and . We find that the mean scaled longitude of perihelion and orbital poles of the detected ETNOs are consistent with a uniform population at a level between 17% and 94% and thus conclude that this sample provides no evidence for angular clustering.
  99. Aguena, M., Benoist, C., da Costa, L., Ogando, R., Gschwend, J., Sampaio-Santos, H., Lima, M., Maia, M., Allam, S., Avila, S., Bacon, D., Bertin, E., Bhargava, S., Brooks, D., Carnero Rosell, A., Carrasco Kind, M., Carretero, J., Costanzi, M., De Vicente, J., Desai, S., Diehl, H., Doel, P., Everett, S., Evrard, A., Ferrero, I., Ferte, A., Flaugher, B., Fosalba, P., Frieman, J., Garcia-Bellido, J., Giles, P., Gruendl, R., Gutierrez, G., Hinton, S., Hollowood, D., Honscheid, K., James, D., Jeltema, T., Kuehn, K., Kuropatkin, N., Lahav, O., Melchior, P., Miquel, R., Morgan, R., Palmese, A., Paz-Chinchon, F., Plazas, A., Romer, A., Sanchez, E., Santiago, B., Schubnell, M., Serrano, S., Sevilla-Noarbe, I., Smith, M., Soares-Santos, M., Suchyta, E., Tarle, G., To, C., Tucker, D., Wilkinson, R., 2021, MNRAS, 502, 4435, The WaZP galaxy cluster sample of the dark energy survey year 1
    We present a new (2+1)D galaxy cluster finder based on photometric redshifts called Wavelet Z Photometric (WaZP) applied to DES first year (Y1A1) data. The results are compared to clusters detected by the South Pole Telescope (SPT) survey and the redMaPPer cluster finder, the latter based on the same photometric data. WaZP searches for clusters in wavelet-based density maps of galaxies selected in photometric redshift space without any assumption on the cluster galaxy populations. The comparison to other cluster samples was performed with a matching algorithm based on angular proximity and redshift difference of the clusters. It led to the development of a new approach to match two optical cluster samples, following an iterative approach to minimize incorrect associations. The WaZP cluster finder applied to DES Y1A1 galaxy survey (1511.13 deg2 up to mi = 23 mag) led to the detection of 60 547 galaxy clusters with redshifts 0.05 < z < 0.9 and richness Ngals 5. Considering the overlapping regions and redshift ranges between the DES Y1A1 and SPT cluster surveys, all SZ based SPT clusters are recovered by the WaZP sample. The comparison between WaZP and redMaPPer cluster samples showed an excellent overall agreement for clusters with richness Ngals ( for redMaPPer) greater than 25 (20), with 95 per cent recovery on both directions. Based on the cluster cross-match, we explore the relative fragmentation of the two cluster samples and investigate the possible signatures of unmatched clusters.
  100. Benedict, G., Franz, O., Horch, E., Torres, G., Prato, L., 2021, csss, 29, The Multiple-component Binary Hyad, vA 351
    We extend results first announced by Franz et al. (1998), that identified vA351 = H346 in the Hyades as a multiple star system containing a white dwarf. With Hubble Space Telescope Fine Guidance Sensor fringe tracking and scanning, and more recent speckle observations, all spanning 20.7 years, we establish a parallax, relative orbit, and mass fraction for two components, with a period, P = 2.70y and total mass 2.1M. With ground-based radial velocities, we find that component B consists of BC, two M dwarf stars orbiting with a very short period (PBC = 0.749 days), having a mass ratio MC/MB=0.95. We confirm that the total mass of the system can only be reconciled with the distance and component photometry by including a fainter, higher mass component. The quadruple system consists of three M dwarfs (A,B,C) and one white dwarf (D); MA=0.57M, MB=0.48M, and MC=0.45M. The WD mass, 0.53M, comes from cooling models, an assumed Hyades age of 670My, and consistency with all previous and derived astrometric, photometric, and RV results. Velocities from H and He I emission lines confirm the BC period derived from absorption lines, with similar (HeI) and higher (H) velocity amplitudes. We ascribe the larger H amplitude to emission from a region each component shadows from the other, depending on the line of sight.
  101. Coffaro, M., Stelzer, B., Orlando, S., Hall, J., Metcalfe, T., Wolter, U., Mittag, M., Sanz-Forcada, J., Ducci, L., 2021, csss, 38, The Sun as a young star: reproducing the X-ray cycle of Eridani with solar magnetic structures
    Epsilon Eri is a young solar-like star with a ~3 yr X-ray activity cycle, detected by us for the first time in a dedicated XMM-Newton long-term monitoring campaign. The magnetic structures on the Sun are intimately linked to the 11-yr activity cycle and they were spatially and temporally resolved throughout the solar cycle. However, for other stars these structures can not be spatially resolved with present-day X-ray instruments. We have, thus, developed a new technique which allows us to reproduce the stellar X-ray variability in terms of time-variations in the coverage of the corona with the same kind of magnetic structures observed on the Sun: active regions (ARs), cores of active regions (COs) and flares (FLs). This poster presents this new method and the results we obtained for the case of Epsilon Eri. Our approach is to simulate a grid of emission measure distributions (EMDs) derived from the analysis of regions observed in the solar corona to artificially reproduce a solar-like corona with the physical characteristics of Epsilon Eri. The three magnetic structures are allowed to contribute to the total coronal EMD with varying area coverage fraction. Thus, from a comparison between these pseudo-solar EMDs and the observations of Epsilon Eri, we are able to associate to each state of the X-ray activity cycle of Epsilon Eri the percentage of ARs, COs and FLs on the corona of the star. The observed amplitude of the X-ray luminosity in the cycle of Epsilon Eri is much smaller than on the Sun. Our analysis provides a physical explanation for this: the simulated EMDs indicate that in all phases of the X-ray cycle a large portion of the corona of Epsilon Eri is covered by active structures. Therefore, there is little space for adding further magnetic structures in the cycle maximum. In the future, this method will be applied to other stars providing an important contribution to better understand the solar-stellar corona connection.
  102. Stahl, A., Tang, S., Johns-Krull, C., Prato, L., Llama, J., Mace, G., Lee, J., Oh, H., Luna, J., Jaffe, D., 2021, csss, 66, IGRINS RV: A Precision RV Pipeline for IGRINS Using Modified Forward-Modeling in the Near-Infrared
    Applications of the radial velocity (RV) technique to the near infrared (NIR) are valuable for their diminished susceptibility to the impact of stellar activity and their suitability for studying late-type stars. In this paper, we present the \texttt{IGRINS RV} open source \texttt{python} pipeline for computing infrared RV measurements from reduced spectra taken with IGRINS, a R$\equiv \lambda/\Delta \lambda \sim$45,000 spectrograph with simultaneous coverage of the H band (1.49--1.80 $\mu$m) and K band (1.96--2.46 $\mu$m). Using a modified forward modeling technique, we construct high resolution telluric templates from A0 standard observations on a nightly basis to provide a source of common-path wavelength calibration while mitigating the need to mask or correct for telluric absorption. A0 standard observations are also used to model the variations in instrumental resolution across the detector, including a yearlong period when the K band was defocused. Without any additional instrument hardware, such as a gas cell or laser frequency comb, we are able to achieve precisions of 26.8 \ms in the K band and 31.1 \ms in the H band for narrow-line hosts. These precisions are validated by a monitoring campaign of two RV standard stars as well as the successful retrieval of planet-induced RV signals for both HD\,189733 and $\tau$\,Boo\,A; furthermore, our results affirm the presence of the Rossiter-McLaughlin effect for HD\,189733. The \texttt{IGRINS RV} pipeline extends another important science capability to IGRINS, with publicly available software designed for widespread use.
  103. Biddle, L., Llama, J., Cameron, A., Prato, L., Jardine, M., Johns-Krull, C., 2021, csss, 214, Amplitude Modulation of Short Timescale Variability
    Variability of Classical T Tauri Systems occurs over a range of timescales. A common source of this variability is caused by accretion shocks, which can hamper the detectability of young planets within these systems, further motivating the characterization of accretion signatures. We present analysis of small-amplitude photometric variability in the K2 lightcurve of CI Tau occurring on timescales of 1 d. Our findings reveal that the amplitude of this stochastic variability exhibits the same periodic signatures as detected in the large-amplitude variability, indicating that the physical mechanism modulating these brightness features is the same.
  104. Clark, C., van Belle, G., Horch, E., von Braun, K., 2021, csss, 236, First Results from the POKEMON Speckle Survey of Nearby M-dwarfs
    We present the first results from the POKEMON (Pervasive Overview of Kompanions of Every M-dwarf in Our Neighborhood) survey, the largest speckle survey of stellar multiplicity ever produced for the objects that comprise over 70% of the stars in our galaxy: the M-dwarfs. We have conducted a volume-limited survey through M9 that inspected, at diffraction-limited resolution, every M-dwarf out to 15pc, with additional brighter targets to 25pc. POKEMON utilized the Differential Speckle Survey Instrument (DSSI) at the 4.3m Lowell Discovery Telescope, along with the NN-Explore Exoplanet Stellar Speckle Imager (NESSI) on the 3.5-m WIYN telescope. We report the discovery of 30+ new companions to these nearby M-dwarfs. Given the priority these targets have for exoplanet studies with TESS, and in the future JWST - and the degree to which initially undetected multiplicity has skewed Kepler results - a comprehensive survey of our nearby low-mass neighbors provides a homogeneous, complete catalog of fundamental utility. Prior knowledge of secondary objects - or robust non-detections, as captured by this survey - immediately clarify the nature of exoplanet transit detections from these current and upcoming missions.
  105. Hartman, Z., Lepine, S., van Belle, G., 2021, csss, 245, Confirming the "Lobster" Diagram: Unresolved Companions in K+K Wide Binaries with TESS and Kepler
    We present an analysis of TESS and Kepler light curves for 4,947 K+K wide binaries from the SUPERWIDE Catalog. We use these systems to examine the usefulness of the "Lobster" diagram, a plot which allows for the identification of over-luminous components in wide binaries using Gaia data alone. These over-luminous components are believed to be unresolved binaries, making the wide binaries a triple (or even quadruple) system. To confirm this, we search through the high cadence light curves from TESS, Kepler and K2 and the data products produced from the MIT Quick Lookup Pipeline for signals from eclipses and from rotational modulation in spotted fast-rotating stars. We identify 64 eclipsing systems and 115 systems showing signs of fast rotation. We highlight the systems containing eclipsing binaries, fast rotators (P<5 days) and rotators (P>5 days) on the "Lobster" diagram. Eclipsing binaries are overwhelmingly found to be in areas that show a component is over-luminous. Fast rotators are also more likely to be found in these areas while stars showing rotation with periods > 5 days are more likely to be found where true wide binaries are believed to be located.
  106. Carvalho, A., Johns-Krull, C., Prato, L., Anderson, J., 2021, ApJ, 910, 33, Radial Velocity Monitoring of the Young Star Hubble 4: Disentangling Star-spot Lifetimes from Orbital Motion
    We studied the weak-lined T Tauri star Hubble 4, a known long-period binary, and its star-spot phenomena. We used optical radial velocity (RV) data taken over a span of 14 yr (2004-2010, 2017-2019) at the McDonald Observatory 2.7 m Harlan J. Smith Telescope and single epoch imaging from the Hubble Space Telescope (HST)/Wide Field Camera 3 instrument. The observed and apparent RV variations show contributions, respectively, from the binary motion as well as from a large spot group on one of the stars, presumed to be the primary. Fitting and removing the orbital signal from the RVs, we found the lower bound on the lifetime of a previously identified large spot group on the surface of the star to be at least 5.1 yr. An 5 yr lower limit is a long, but not unprecedented, duration for a single spot group. The later epoch data indicate significant spot evolution has occurred, placing an upper bound on the spot group lifetime at 12 yr. We find that pre-main-sequence evolutionary models for the age of Taurus (2 Myr), combined with component mass estimates from the literature, permit us to reproduce the HST relative photometry and the binary-induced contribution to the apparent RV variations. The long-lived star spot we find on Hubble 4 has significant implications for dynamo models in young stars, as it adds evidence for long lifetimes of magnetic field topologies. There are also significant implications for young star exoplanet searches, as long-lived coherent RV signals may be spot induced and not the result of planetary motion. * This paper includes data taken at The McDonald Observatory of The University of Texas at Austin.
  107. Gustafsson, A., Moskovitz, N., Cushing, M., Bida, T., Dunham, E., Roe, H., 2021, PASP, 133, 035001, Science Commissioning of NIHTS: The Near-infrared High Throughput Spectrograph on the Lowell Discovery Telescope
    The Near-Infrared High Throughput Spectrograph (NIHTS) is in operation on the 4.3 m Lowell Discovery Telescope (LDT) in Happy Jack, AZ. NIHTS is a low-resolution spectrograph (R 200) that operates from 0.86 to 2.45 microns. NIHTS is fed by a custom dichroic mirror which reflects near-infrared wavelengths to the spectrograph and transmits the visible to enable simultaneous imaging with the Large Monolithic Imager (LMI), an independent visible wavelength camera. The combination of premier tracking and acquisition capabilities of the LDT, a several arcminutes field of view on LMI, and high spectral throughput on NIHTS enables novel studies of a number of astrophysical and planetary objects including Kuiper Belt Objects, asteroids, comets, low mass stars, and exoplanet hosts stars. We present a summary of NIHTS operations, commissioning, data reduction procedures with two approaches for the correction of telluric absorption features, and an overview of select science cases that will be pursued by Lowell Observatory, Northern Arizona University, and LDT partners.
  108. Hilton, M., Sifon, C., Naess, S., Madhavacheril, M., Oguri, M., Rozo, E., Rykoff, E., Abbott, T., Adhikari, S., Aguena, M., Aiola, S., Allam, S., Amodeo, S., Amon, A., Annis, J., Ansarinejad, B., Aros-Bunster, C., Austermann, J., Avila, S., Bacon, D., Battaglia, N., Beall, J., Becker, D., Bernstein, G., Bertin, E., Bhandarkar, T., Bhargava, S., Bond, J., Brooks, D., Burke, D., Calabrese, E., Carrasco Kind, M., Carretero, J., Choi, S., Choi, A., Conselice, C., da Costa, L., Costanzi, M., Crichton, D., Crowley, K., Dunner, R., Denison, E., Devlin, M., Dicker, S., Diehl, H., Dietrich, J., Doel, P., Duff, S., Duivenvoorden, A., Dunkley, J., Everett, S., Ferraro, S., Ferrero, I., Ferte, A., Flaugher, B., Frieman, J., Gallardo, P., Garcia-Bellido, J., Gaztanaga, E., Gerdes, D., Giles, P., Golec, J., Gralla, M., Grandis, S., Gruen, D., Gruendl, R., Gschwend, J., Gutierrez, G., Han, D., Hartley, W., Hasselfield, M., Hill, J., Hilton, G., Hincks, A., Hinton, S., Ho, S., Honscheid, K., Hoyle, B., Hubmayr, J., Huffenberger, K., Hughes, J., Jaelani, A., Jain, B., James, D., Jeltema, T., Kent, S., Knowles, K., Koopman, B., Kuehn, K., Lahav, O., Lima, M., Lin, Y., Lokken, M., Loubser, S., MacCrann, N., Maia, M., Marriage, T., Martin, J., McMahon, J., Melchior, P., Menanteau, F., Miquel, R., Miyatake, H., Moodley, K., Morgan, R., Mroczkowski, T., Nati, F., Newburgh, L., Niemack, M., Nishizawa, A., Ogando, R., Orlowski-Scherer, J., Page, L., Palmese, A., Partridge, B., Paz-Chinchon, F., Phakathi, P., Plazas, A., Robertson, N., Romer, A., Carnero Rosell, A., Salatino, M., Sanchez, E., Schaan, E., Schillaci, A., Sehgal, N., Serrano, S., Shin, T., Simon, S., Smith, M., Soares-Santos, M., Spergel, D., Staggs, S., Storer, E., Suchyta, E., Swanson, M., Tarle, G., Thomas, D., To, C., Trac, H., Ullom, J., Vale, L., Van Lanen, J., Vavagiakis, E., De Vicente, J., Wilkinson, R., Wollack, E., Xu, Z., Zhang, Y., 2021, ApJS, 253, 3, The Atacama Cosmology Telescope: A Catalog of >4000 Sunyaev-Zeldovich Galaxy Clusters
    We present a catalog of 4195 optically confirmed Sunyaev-Zeldovich (SZ) selected galaxy clusters detected with signal-to-noise ratio >4 in 13,211 deg2 of sky surveyed by the Atacama Cosmology Telescope (ACT). Cluster candidates were selected by applying a multifrequency matched filter to 98 and 150 GHz maps constructed from ACT observations obtained from 2008 to 2018 and confirmed using deep, wide-area optical surveys. The clusters span the redshift range 0.04 < z < 1.91 (median z = 0.52). The catalog contains 222 z > 1 clusters, and a total of 868 systems are new discoveries. Assuming an SZ signal versus mass-scaling relation calibrated from X-ray observations, the sample has a 90% completeness mass limit of M 500c > 3.8 1014 M , evaluated at z = 0.5, for clusters detected at signal-to-noise ratio >5 in maps filtered at an angular scale of 2.4. The survey has a large overlap with deep optical weak-lensing surveys that are being used to calibrate the SZ signal mass-scaling relation, such as the Dark Energy Survey (4566 deg2), the Hyper Suprime-Cam Subaru Strategic Program (469 deg2), and the Kilo Degree Survey (825 deg2). We highlight some noteworthy objects in the sample, including potentially projected systems, clusters with strong lensing features, clusters with active central galaxies or star formation, and systems of multiple clusters that may be physically associated. The cluster catalog will be a useful resource for future cosmological analyses and studying the evolution of the intracluster medium and galaxies in massive clusters over the past 10 Gyr.
  109. Nadler, E., Drlica-Wagner, A., Bechtol, K., Mau, S., Wechsler, R., Gluscevic, V., Boddy, K., Pace, A., Li, T., McNanna, M., Riley, A., Garcia-Bellido, J., Mao, Y., Green, G., Burke, D., Peter, A., Jain, B., Abbott, T., Aguena, M., Allam, S., Annis, J., Avila, S., Brooks, D., Carrasco Kind, M., Carretero, J., Costanzi, M., da Costa, L., De Vicente, J., Desai, S., Diehl, H., Doel, P., Everett, S., Evrard, A., Flaugher, B., Frieman, J., Gerdes, D., Gruen, D., Gruendl, R., Gschwend, J., Gutierrez, G., Hinton, S., Honscheid, K., Huterer, D., James, D., Krause, E., Kuehn, K., Kuropatkin, N., Lahav, O., Maia, M., Marshall, J., Menanteau, F., Miquel, R., Palmese, A., Paz-Chinchon, F., Plazas, A., Romer, A., Sanchez, E., Scarpine, V., Serrano, S., Sevilla-Noarbe, I., Smith, M., Soares-Santos, M., Suchyta, E., Swanson, M., Tarle, G., Tucker, D., Walker, A., Wester, W., DES Collaboration, 2021, PhRvL, 126, 091101, Constraints on Dark Matter Properties from Observations of Milky Way Satellite Galaxies
    We perform a comprehensive study of Milky Way (MW) satellite galaxies to constrain the fundamental properties of dark matter (DM). This analysis fully incorporates inhomogeneities in the spatial distribution and detectability of MW satellites and marginalizes over uncertainties in the mapping between galaxies and DM halos, the properties of the MW system, and the disruption of subhalos by the MW disk. Our results are consistent with the cold, collisionless DM paradigm and yield the strongest cosmological constraints to date on particle models of warm, interacting, and fuzzy dark matter. At 95% confidence, we report limits on (i) the mass of thermal relic warm DM, mWDM>6.5 keV (free-streaming length, fs10 h-1 kpc ), (ii) the velocity-independent DM-proton scattering cross section, 0<8.8 10-29 cm2 for a 100 MeV DM particle mass [DM-proton coupling, cp(0.3 GeV )-2], and (iii) the mass of fuzzy DM, m>2.9 10-21 eV (de Broglie wavelength, dB0.5 kpc ). These constraints are complementary to other observational and laboratory constraints on DM properties.
  110. Oszkiewicz, D., Wilawer, E., Podlewska-Gaca, E., Kryszczynska, A., Kwiatkowski, T., Troianskyi, V., Kolenczuk, P., Fohring, D., Galad, A., Skiff, B., Geier, S., Borczyk, W., Moskovitz, N., Gajdos, S., Vilagi, J., Polcic, L., Kashuba, V., Benishek, V., Shevchenko, V., 2021, Icar, 357, 114158, First survey of phase curves of V-type asteroids
    The V-type asteroids are of major scientific interest as they may sample multiple differentiated planetesimals. Determination of their physical properties is crucial for understanding the diversity and multiplicity of planetesimals. Previous studies have suggested distinct polarimetric behaviours for the V-type asteroids. Similarly to phase-polarization curves, asteroid phase-magnitude curves (hereinafter called "phase curves") are also diagnostic of surface and regolith properties, and can be used to unveil a variety of distinct behaviours.

    We present well determined phase curves for ~20 V-type asteroids for the first time. Their phase curve parameters are consistent with those for moderate and high albedo asteroids. The computed median G12 parameter for the V-type asteroids is G12 = 0.14. We do not find substantial evidence for any clustering into distinct phase curve parameters groups. Only one asteroid (2763) Jeans shows exceptionally high G2 value. The derived median G12 may be used in single parameters fitting of V-type asteroids.

  111. Johnson, J., Grundy, W., Lemmon, M., Liang, W., Bell, J., Hayes, A., Deen, R., 2021, Icar, 357, 114261, Spectrophotometric properties of materials observed by Pancam on the Mars Exploration Rovers: 4. Final mission observations
    The last sets of Panoramic Camera (Pancam) visible/near-infrared (432-1009 nm) multispectral observations made under varying viewing and illumination geometries by the Mars Exploration Rovers Spirit and Opportunity were examined using radiative transfer models to study the surface scattering and microphysical nature of rock and soil units at both sites. Nearly 12,000 individual measurements were collected for this study of soil, dust, and rock units over phase angles of ~0 to ~150. Images were acquired on sols 1944-1946 (June 2009) at Troy, the final resting place of Spirit on the western side of Home Plate in Gusev crater, and by Opportunity at three locations on the western rim of Endeavour crater in Meridiani Planum between sols 2785 (November 2011) and 3867 (December 2014). Sky models were developed from observations of atmospheric opacity, which enabled corrections for diffuse skylight when combined with surface facet orientations determined from stereo images. Model results were improved by removing data affected by scattered light evident in some high phase angle images (resulting from minor dust contamination on the camera windows). At Troy, relatively dust-free "gray" rock units exhibited narrow, forward scattering behaviors akin to previous analyses of similar gray rock units at Gusev crater. Soils and "red" rocks coated with greater amounts of dust were more backscattering. Red rocks exhibited higher single scattering albedo (w), macroscopic roughness (), and opposition effect width (h) parameters, indicative of rough, low-porosity surfaces perhaps with more uniform grain size distributions. At Meridiani Planum, rubbly soils near Sao Gabriel crater and Cape Tribulation exhibited w values typical of previous soil analyses. However, the large drift "dust" deposits found in depressions on the northern tip of Cape York near Turkey Haven demonstrated elevated w values with a downturn toward 1009 nm, consistent with minor hydration of these materials. The dust deposits were modeled with the lowest values and highest h values of all soil units analyzed during the Opportunity mission, indicative of a smooth surface with homogeneous grain size distribution and/or lower porosity than other units. The dust unit scattering function was dissimilar to those for atmospheric and airfall-deposited dusts, however, suggesting that the originally deposited materials had been modified, perhaps by hydration and ongoing aeolian effects. Analyses of phase curve ratios among the units studied here and from laboratory data of analog soils suggested that surface scattering is a major control on the peak phase angle position of the "arch" in phase curve ratios, alongside the effects of particle-scale roughness.
  112. Kelsey, L., Sullivan, M., Smith, M., Wiseman, P., Brout, D., Davis, T., Frohmaier, C., Galbany, L., Grayling, M., Gutierrez, C., Hinton, S., Kessler, R., Lidman, C., Moller, A., Sako, M., Scolnic, D., Uddin, S., Vincenzi, M., Abbott, T., Aguena, M., Allam, S., Annis, J., Avila, S., Bacon, D., Bertin, E., Brooks, D., Burke, D., Carnero Rosell, A., Carrasco Kind, M., Carretero, J., Castander, F., Costanzi, M., da Costa, L., Desai, S., Diehl, H., Doel, P., Everett, S., Ferrero, I., Ferte, A., Flaugher, B., Fosalba, P., Garcia-Bellido, J., Gerdes, D., Gruen, D., Gruendl, R., Gschwend, J., Gutierrez, G., Hollowood, D., Honscheid, K., James, D., Kim, A., Kuehn, K., Kuropatkin, N., Lahav, O., Lima, M., Marshall, J., Martini, P., Menanteau, F., Miquel, R., Morgan, R., Ogando, R., Palmese, A., Paz-Chinchon, F., Plazas, A., Romer, A., Sanchez, C., Sanchez, E., Serrano, S., Sevilla-Noarbe, I., Suchyta, E., Tarle, G., Thomas, D., To, C., Varga, T., Walker, A., Wilkinson, R., DES Collaboration, 2021, MNRAS, 501, 4861, The effect of environment on Type Ia supernovae in the Dark Energy Survey three-year cosmological sample
    Analyses of Type Ia supernovae (SNe Ia) have found puzzling correlations between their standardized luminosities and host galaxy properties: SNe Ia in high-mass, passive hosts appear brighter than those in lower mass, star-forming hosts. We examine the host galaxies of SNe Ia in the Dark Energy Survey 3-yr spectroscopically confirmed cosmological sample, obtaining photometry in a series of 'local' apertures centred on the SN, and for the global host galaxy. We study the differences in these host galaxy properties, such as stellar mass and rest-frame U - R colours, and their correlations with SN Ia parameters including Hubble residuals. We find all Hubble residual steps to be >3 in significance, both for splitting at the traditional environmental property sample median and for the step of maximum significance. For stellar mass, we find a maximal local step of 0.098 0.018 mag; 0.03 mag greater than the largest global stellar mass step in our sample (0.070 0.017 mag). When splitting at the sample median, differences between local and global U - R steps are small, both 0.08 mag, but are more significant than the global stellar mass step (0.057 0.017 mag). We split the data into sub-samples based on SN Ia light-curve parameters: stretch (x1) and colour (c), finding that redder objects (c > 0) have larger Hubble residual steps, for both stellar mass and U - R, for both local and global measurements, of 0.14 mag. Additionally, the bluer (star-forming) local environments host a more homogeneous SN Ia sample, with local U - R rms scatter as low as 0.084 0.017 mag for blue (c < 0) SNe Ia in locally blue U - R environments.
  113. Person, M., Bosh, A., Zuluaga, C., Sickafoose, A., Levine, S., Pasachoff, J., Babcock, B., Dunham, E., McLean, I., Wolf, J., Abe, F., Becklin, E., Bida, T., Bright, L., Brothers, T., Christie, G., Durst, R., Gilmore, A., Hamilton, R., Harris, H., Johnson, C., Kilmartin, P., Kosiarek, M., Leppik, K., Logsdon, S., Lucas, R., Mathers, S., Morley, C., Nelson, P., Ngan, H., Pfuller, E., Natusch, T., Sallum, S., Savage, M., Seeger, C., Siu, H., Stockdale, C., Suzuki, D., Thanathibodee, T., Tilleman, T., Tristram, P., Vacca, W., Van Cleve, J., Varughese, C., Weisenbach, L., Widen, E., Wiedemann, M., 2021, Icar, 356, 113572, Haze in Pluto's atmosphere: Results from SOFIA and ground-based observations of the 2015 June 29 Pluto occultation
    On UT 29 June 2015, the occultation by Pluto of a bright star (r = 11.9) was observed from the Stratospheric Observatory for Infrared Astronomy (SOFIA) and several ground-based stations in New Zealand and Australia. Pre-event astrometry allowed for an in-flight update to the SOFIA team with the result that SOFIA was deep within the central flash zone (~22 km from center). Analysis of the combined data leads to the result that Pluto's middle atmosphere is essentially unchanged from 2011 and 2013 (Person et al. 2013; Bosh et al. 2015); there has been no significant expansion or contraction of the atmosphere. Additionally, our multi-wavelength observations allow us to conclude that a haze component in the atmosphere is required to reproduce the light curves obtained. This haze scenario has implications for understanding the photochemistry of Pluto's atmosphere.
  114. Lewis, B., Stansberry, J., Holler, B., Grundy, W., Schmitt, B., Protopapa, S., Lisse, C., Stern, S., Young, L., Weaver, H., Olkin, C., Ennico, K., New Horizons Science Team, 2021, Icar, 356, 113633, Distribution and energy balance of Pluto's nitrogen ice, as seen by New Horizons in 2015
    Pluto's surface is geologically complex because of volatile ices that are mobile on seasonal and longer time scales. Here we analyzed New Horizons LEISA spectral data to globally map the nitrogen ice, including nitrogen with methane diluted in it. Our goal was to learn about the seasonal processes influencing ice redistribution, to calculate the globally averaged energy balance, and to place a lower limit on Pluto's N2 inventory. We present the average latitudinal distribution of nitrogen and investigate the relationship between its distribution and topography on Pluto by using maps that include the shifted bands of methane in solid solution with nitrogen (which are much stronger than the 2.15- m nitrogen band) to more completely map the distribution of the nitrogen ice. We find that the global average bolometric albedo is 0 . 83 0 . 11 , similar to that inferred for Triton, and that a significant fraction of Pluto's N2 is stored in Sputnik Planitia. We also used the encounter-hemisphere distribution of nitrogen ice to infer the latitudinal distribution of nitrogen over the rest of Pluto, allowing us to calculate the global energy balance. Under the assumption that Pluto's nitrogen-dominated 11.5 bar atmosphere is in vapor pressure equilibrium with the nitrogen ice, the ice temperature is 36 . 93 0 . 10 K, as measured by New Horizons' REX instrument. Combined with our global energy balance calculation, this implies that the average bolometric emissivity of Pluto's nitrogen ice is probably in the range 0.47-0.72. This is consistent with the low emissivities estimated for Triton based on Voyager results, and may have implications for Pluto's atmospheric seasonal variations, as discussed below. The global pattern of volatile transport at the time of the encounter was from north to south, and the transition between condensation and sublimation within Sputnik Planitia is correlated with changes in the grain size and CH4 concentration derived from the spectral maps. The low emissivity of Pluto's N2 ice suggests that Pluto's atmosphere may undergo an extended period of constant pressure even as Pluto recedes from the Sun in its orbit.
  115. Cruikshank, D., Dalle Ore, C., Scipioni, F., Beyer, R., White, O., Moore, J., Grundy, W., Schmitt, B., Runyon, K., Keane, J., Robbins, S., Stern, S., Bertrand, T., Beddingfield, C., Olkin, C., Young, L., Weaver, H., Ennico, K., 2021, Icar, 356, 113786, Cryovolcanic flooding in Viking Terra on Pluto
    A prominent fossa trough (Uncama Fossa) and adjacent 28-km diameter impact crater (Hardie) in Pluto's Viking Terra, as seen in the high-resolution images from the New Horizons spacecraft, show morphological evidence of in-filling with a material of uniform texture and red-brown color. A linear fissure parallel to the trough may be the source of a fountaining event yielding a cryoclastic deposit having the same composition and color properties as is found in the trough and crater. Spectral maps of this region with the New Horizons LEISA instrument reveal the spectral signature of H2O ice in these structures and in distributed patches in the adjacent terrain in Viking Terra. A detailed statistical analysis of the spectral maps shows that the colored H2O ice filling material also carries the 2.2-m signature of an ammoniated component that may be an ammonia hydrate (NH3nH2O) or an ammoniated salt. This paper advances the view that the crater and fossa trough have been flooded by a cryolava debouched from Pluto's interior along fault lines in the trough and in the floor of the impact crater. The now frozen cryolava consisted of liquid H2O infused with the red-brown pigment presumed to be a tholin, and one or more ammoniated compounds. Although the abundances of the pigment and ammoniated compounds entrained in, or possibly covering, the H2O ice are unknown, the strong spectral bands of the H2O ice are clearly visible. In consideration of the factors in Pluto's space environment that are known to destroy ammonia and ammonia-water mixtures, the age of the exposure is of order 109 years. Ammoniated salts may be more robust, and laboratory investigations of these compounds are needed.
  116. Gabasova, L., Schmitt, B., Grundy, W., Bertrand, T., Olkin, C., Spencer, J., Young, L., Ennico, K., Weaver, H., Stern, S., New Horizons Composition Team, 2021, Icar, 356, 113833, Global compositional cartography of Pluto from intensity-based registration of LEISA data
    In 2015 the New Horizons spacecraft reached the Pluto system and returned unprecedentedly detailed measurements of its surface properties. These measurements have already been integrated into global reflectance, topography and narrow-band multispectral surface maps. However, analysis of the hyperspectral data from the Ralph/LEISA infrared spectrometer, which lets us analyse the surface composition, has until now been confined to the high-resolution encounter hemisphere of Pluto. We use an innovative technique - intensity-based registration - to co-register this high-resolution data with lower-resolution measurements taken during the spacecraft's approach, and present the first global qualitative composition maps for CH<!--l. 82 -->4, N<!--l. 82 -->2 and H<!--l. 82 -->2O ice, and a tholin-like red material. We compare these maps with the other maps produced for Pluto and study the global extent of the previously-described latitudinal distribution of the surface components, which is relatively longitudinally constant with the exception of Sputnik Planitia. We also correlate these compositional components with geological features and propose physical interpretations, which include: CH<!--l. 82 -->4-ice-rich dissected plateaus in high northern latitudes, CH<!--l. 82 -->4-rich eroded terrain with N<!--l. 82 -->2-rich infill in medium northern latitudes, CH<!--l. 82 -->4-rich bladed terrain in low northern latitudes, and a red material belt overlaying H<!--l. 82 -->2O ice in low southern latitudes.
  117. Lisse, C., Young, L., Cruikshank, D., Sandford, S., Schmitt, B., Stern, S., Weaver, H., Umurhan, O., Pendleton, Y., Keane, J., Gladstone, G., Parker, J., Binzel, R., Earle, A., Horanyi, M., El-Maarry, M., Cheng, A., Moore, J., McKinnon, W., Grundy, W., Kavelaars, J., Linscott, I., Lyra, W., Lewis, B., Britt, D., Spencer, J., Olkin, C., McNutt, R., Elliott, H., Dello-Russo, N., Steckloff, J., Neveu, M., Mousis, O., 2021, Icar, 356, 114072, On the origin & thermal stability of Arrokoth's and Pluto's ices
    In this paper we discuss in a thermodynamic, geologically empirical way the long-term nature of the stable majority ices that could be present in Kuiper Belt object (KBO) 2014 MU69 (also called Arrokoth; hereafter "MU69") after its 4.6 Gyr residence in the Edgeworth-Kuiper belt (EKB) as a cold classical object. We compare the upper bounds for the gas production rate (~1024 molecules/s) measured by the New Horizons (NH) spacecraft flyby on 01 Jan 2019 to estimates for the outgassing flux rates from a suite of common cometary and KBO ices at the average ~ 40 K sunlit surface temperature of MU69, but do not find the upper limit very constraining except for the most volatile of species (e.g. CO, N2, CH4). More constraining is the stability versus sublimation into vacuum requirement over Myr to Gyr, and from this we find only 3 common ices that are truly refractory: HCN, CH3OH, and H2O (in order of increasing stability), while NH3 and H2CO ices are marginally stable and may be removed by any positive temperature excursions in the EKB, as produced every 108-109 years by nearby supernovae and passing O/B stars. To date the NH team has reported the presence of abundant CH3OH and H2O on MU69's surface (Stern et al., 2019; Grundy et al., 2020). NH3 has been searched for, but not found. We predict that future absorption feature detections, if any are ever derived from higher signal-to-noise ratio spectra, will be due to an HCN or poly-H2CO based species. Consideration of the conditions present in the EKB region during the formation era of MU69 lead us to state that it is highly likely that it "formed in the dark", in an optically thick mid-plane, unable to see the nascent, variable, highly luminous Young Stellar Object (YSO)/TTauri Sun, and that KBOs contain HCN and CH3OH ice phases in addition to the H2O ice phases found in their short period (SP) comet descendants. Finally, when we apply our ice thermal stability analysis to bodies/populations related to MU69, we find that methanol ice is likely ubiquitous in the outer solar system; that if Pluto isn't a fully differentiated body, then it must have gained its hypervolatile ices from proto-planetary disk (PPD) sources in the first few Myr of the solar system's existence; and that hypervolatile rich, highly primordial comet C/2016 R2 was placed onto an Oort Cloud orbit on a similar few Myr timescale.
  118. Showalter, M., Benecchi, S., Buie, M., Grundy, W., Keane, J., Lisse, C., Olkin, C., Porter, S., Robbins, S., Singer, K., Verbiscer, A., Weaver, H., Zangari, A., Hamilton, D., Kaufmann, D., Lauer, T., Mehoke, D., Mehoke, T., Spencer, J., Throop, H., Parker, J., Stern, S., New Horizons Geology, G., 2021, Icar, 356, 114098, A statistical review of light curves and the prevalence of contact binaries in the Kuiper Belt
    We investigate what can be learned about a population of distant Kuiper Belt Objects (KBOs) by studying the statistical properties of their light curves. Whereas others have successfully inferred the properties of individual, highly variable KBOs, we show that the fraction of KBOs with low amplitudes also provides fundamental information about a population. Each light curve is primarily the result of two factors: shape and orientation. We consider contact binaries and ellipsoidal shapes, with and without flattening. After developing the mathematical framework, we apply it to the existing body of KBO light curve data. Principal conclusions are as follows. (1) When using absolute magnitude H as a proxy for the sizes of KBOs, it is more accurate to use the maximum of the light curve (minimum H) rather than the mean. (2) Previous investigators have noted that smaller KBOs tend to have higher-amplitude light curves, and have interpreted this as evidence that they are systematically more irregular in shape than larger KBOs; we show that a population of flattened bodies with uniform proportions, independent of size, could also explain this result. (3) Our method of analysis indicates that prior assessments of the fraction of contact binaries in the Kuiper Belt may be artificially low. (4) The pole orientations of some KBOs can be inferred from observed changes in their light curves over time scales of decades; however, we show that these KBOs constitute a biased sample, whose pole orientations are not representative of the population overall. (5) Although surface topography, albedo patterns, limb darkening, and other surface properties can affect individual light curves, they do not have a strong influence on the statistics overall. (6) Photometry from the Outer Solar System Origins Survey (OSSOS) survey is incompatible with previous results and its statistical properties defy easy interpretation. We also discuss the promise of this approach for the analysis of future, much larger data sets such as the one anticipated from the upcoming Vera C. Rubin Observatory.
  119. Singer, K., Grundy, W., White, O., Binzel, R., 2021, Icar, 356, 114269, Introduction to Icarus special issue "Pluto System, Kuiper Belt, and Kuiper Belt Objects"
    The Kuiper belt represents a vast treasure of scientific information. From evolved, planet-sized bodies to relatively primordial smaller bodies, the science that can be done with Kuiper belt data touches many aspects of solar system formation and evolution. These connections along with advancements in telescope and spacecraft technology have created the conditions for a golden age for Kuiper belt science, as exemplified by the breadth and depth of topics in this special issue.
  120. Jarvis, M., Bernstein, G., Amon, A., Davis, C., Leget, P., Bechtol, K., Harrison, I., Gatti, M., Roodman, A., Chang, C., Chen, R., Choi, A., Desai, S., Drlica-Wagner, A., Gruen, D., Gruendl, R., Hernandez, A., MacCrann, N., Meyers, J., Navarro-Alsina, A., Pandey, S., Plazas, A., Secco, L., Sheldon, E., Troxel, M., Vorperian, S., Wei, K., Zuntz, J., Abbott, T., Aguena, M., Allam, S., Avila, S., Bhargava, S., Bridle, S., Brooks, D., Carnero Rosell, A., Carrasco Kind, M., Carretero, J., Costanzi, M., da Costa, L., De Vicente, J., Diehl, H., Doel, P., Everett, S., Flaugher, B., Fosalba, P., Frieman, J., Garcia-Bellido, J., Gaztanaga, E., Gerdes, D., Gutierrez, G., Hinton, S., Hollowood, D., Honscheid, K., James, D., Kent, S., Kuehn, K., Kuropatkin, N., Lahav, O., Maia, M., March, M., Marshall, J., Melchior, P., Menanteau, F., Miquel, R., Ogando, R., Paz-Chinchon, F., Rykoff, E., Sanchez, E., Scarpine, V., Schubnell, M., Serrano, S., Sevilla-Noarbe, I., Smith, M., Suchyta, E., Swanson, M., Tarle, G., Varga, T., Walker, A., Wester, W., Wilkinson, R., Wilkinson, R., DES Collaboration, 2021, MNRAS, 501, 1282, Dark Energy Survey year 3 results: point spread function modelling
    We introduce a new software package for modelling the point spread function (PSF) of astronomical images, called PIFF (PSFs In the Full FOV), which we apply to the first three years (known as Y3) of the Dark Energy Survey (DES) data. We describe the relevant details about the algorithms used by PIFF to model the PSF, including how the PSF model varies across the field of view (FOV). Diagnostic results show that the systematic errors from the PSF modelling are very small over the range of scales that are important for the DES Y3 weak lensing analysis. In particular, the systematic errors from the PSF modelling are significantly smaller than the corresponding results from the DES year one (Y1) analysis. We also briefly describe some planned improvements to PIFF that we expect to further reduce the modelling errors in future analyses.
  121. Sampaio-Santos, H., Zhang, Y., Ogando, R., Shin, T., Golden-Marx, J., Yanny, B., Herner, K., Hilton, M., Choi, A., Gatti, M., Gruen, D., Hoyle, B., Rau, M., De Vicente, J., Zuntz, J., Abbott, T., Aguena, M., Allam, S., Annis, J., Avila, S., Bertin, E., Brooks, D., Burke, D., Carrasco Kind, M., Carretero, J., Chang, C., Costanzi, M., da Costa, L., Diehl, H., Doel, P., Everett, S., Evrard, A., Flaugher, B., Fosalba, P., Frieman, J., Garcia-Bellido, J., Gaztanaga, E., Gerdes, D., Gruendl, R., Gschwend, J., Gutierrez, G., Hinton, S., Hollowood, D., Honscheid, K., James, D., Jarvis, M., Jeltema, T., Kuehn, K., Kuropatkin, N., Lahav, O., Maia, M., March, M., Marshall, J., Miquel, R., Palmese, A., Paz-Chinchon, F., Plazas, A., Sanchez, E., Santiago, B., Scarpine, V., Schubnell, M., Smith, M., Suchyta, E., Tarle, G., Tucker, D., Varga, T., Wechsler, R., 2021, MNRAS, 501, 1300, Is diffuse intracluster light a good tracer of the galaxy cluster matter distribution?
    We explore the relation between diffuse intracluster light (central galaxy included) and the galaxy cluster (baryonic and dark) matter distribution using a sample of 528 clusters at 0.2 z 0.35 found in the Dark Energy Survey (DES) Year 1 data. The surface brightness of the diffuse light shows an increasing dependence on cluster total mass at larger radius, and appears to be self-similar with a universal radial dependence after scaling by cluster radius. We also compare the diffuse light radial profiles to the cluster (baryonic and dark) matter distribution measured through weak lensing and find them to be comparable. The IllustrisTNG galaxy formation simulation, TNG300, offers further insight into the connection between diffuse stellar mass and cluster matter distributions - the simulation radial profile of the diffuse stellar component does not have a similar slope with the total cluster matter content, although that of the cluster satellite galaxies does. Regardless of the radial trends, the amount of diffuse stellar mass has a low-scatter scaling relation with cluster's total mass in the simulation, out-performing the total stellar mass of cluster satellite galaxies. We conclude that there is no consistent evidence yet on whether or not diffuse light is a faithful radial tracer of the cluster matter distribution. Nevertheless, both observational and simulation results reveal that diffuse light is an excellent indicator of the cluster's total mass.
  122. Neugent, K., 2021, ApJ, 908, 87, The Red Supergiant Binary Fraction as a Function of Metallicity in M31 and M33
    Recent work measuring the binary fraction of evolved red supergiants (RSGs) in the Magellanic Clouds points to a value between 15% and 30%, with the majority of the companions being unevolved B-type stars as dictated by stellar evolution. Here I extend this research to the Local Group galaxies M31 and M33 and investigate the RSG binary fraction as a function of metallicity. Recent near-IR photometric surveys of M31 and M33 have led to the identification of a complete sample of RSGs down to a limiting $\mathrm{log}L/{L}_{\odot }\geqslant 4.2$ . To determine the binary fraction of these M31 and M33 RSGs, I used a combination of newly obtained spectroscopy to identify single RSGs and RSG+OB binaries, as well as archival UV, visible, and near-IR photometry to probabilistically classify RSGs as either single or binary based on their colors. I then adjusted the observed RSG+OB binary fraction to account for observational biases. The resulting RSG binary fraction in M33 shows a strong dependence on galactocentric distance, with the inner regions having a much higher binary fraction ( ${41.2}_{-7.3}^{+12.0} \% $ ) than the outer regions ( ${15.9}_{-1.9}^{+12.4} \% $ ). Such a trend is not seen in M31; instead, the binary fraction in lightly reddened regions remains constant at ${33.5}_{-5.0}^{+8.6} \% $ . I conclude that the changing RSG binary fraction in M33 is due to a metallicity dependence, with higher-metallicity environments having higher RSG binary fractions. This dependence most likely stems not from changes in the physical properties of RSGs due to metallicity but from changes in the parent distribution of OB binaries.
  123. Costanzi, M., Saro, A., Bocquet, S., Abbott, T., Aguena, M., Allam, S., Amara, A., Annis, J., Avila, S., Bacon, D., Benson, B., Bhargava, S., Brooks, D., Buckley-Geer, E., Burke, D., Carnero Rosell, A., Carrasco Kind, M., Carretero, J., Choi, A., da Costa, L., Pereira, M., De Vicente, J., Desai, S., Diehl, H., Dietrich, J., Doel, P., Eifler, T., Everett, S., Ferrero, I., Ferte, A., Flaugher, B., Fosalba, P., Frieman, J., Garcia-Bellido, J., Gaztanaga, E., Gerdes, D., Giannantonio, T., Giles, P., Grandis, S., Gruen, D., Gruendl, R., Gupta, N., Gutierrez, G., Hartley, W., Hinton, S., Hollowood, D., Honscheid, K., James, D., Jeltema, T., Krause, E., Kuehn, K., Kuropatkin, N., Lahav, O., Lima, M., MacCrann, N., Maia, M., Marshall, J., Menanteau, F., Miquel, R., Mohr, J., Morgan, R., Myles, J., Ogando, R., Palmese, A., Paz-Chinchon, F., Plazas, A., Rapetti, D., Reichardt, C., Romer, A., Roodman, A., Ruppin, F., Salvati, L., Samuroff, S., Sanchez, E., Scarpine, V., Serrano, S., Sevilla-Noarbe, I., Singh, P., Smith, M., Soares-Santos, M., Stark, A., Suchyta, E., Swanson, M., Tarle, G., Thomas, D., To, C., Tucker, D., Varga, T., Wechsler, R., Zhang, Z., DES, SPT Collaborations, 2021, PhRvD, 103, 043522, Cosmological constraints from DES Y1 cluster abundances and SPT multiwavelength data
    We perform a joint analysis of the counts of redMaPPer clusters selected from the Dark Energy Survey (DES) year 1 data and multiwavelength follow-up data collected within the 2500 deg2 South Pole Telescope (SPT) Sunyaev-Zel'dovich (SZ) survey. The SPT follow-up data, calibrating the richness-mass relation of the optically selected redMaPPer catalog, enable the cosmological exploitation of the DES cluster abundance data. To explore possible systematics related to the modeling of projection effects, we consider two calibrations of the observational scatter on richness estimates: a simple Gaussian model which account only for the background contamination (BKG), and a model which further includes contamination and incompleteness due to projection effects (PRJ). Assuming either a CDM +m or w CDM +m cosmology, and for both scatter models, we derive cosmological constraints consistent with multiple cosmological probes of the low and high redshift Universe, and in particular with the SPT cluster abundance data. This result demonstrates that the DES Y1 and SPT cluster counts provide consistent cosmological constraints, if the same mass calibration data set is adopted. It thus supports the conclusion of the DES Y1 cluster cosmology analysis which interprets the tension observed with other cosmological probes in terms of systematics affecting the stacked weak lensing analysis of optically selected low-richness clusters. Finally, we analyze the first combined optically SZ selected cluster catalog obtained by including the SPT sample above the maximum redshift probed by the DES Y1 redMaPPer sample (z =0.65 ). Besides providing a mild improvement of the cosmological constraints, this data combination serves as a stricter test of our scatter models: the PRJ model, providing scaling relations consistent between the two abundance and multiwavelength follow-up data, is favored over the BKG model.
  124. Strauss, R., Leiva, R., Keller, J., Wilde, E., Buie, M., Weryk, R., Kavelaars, J., Bridges, T., Wasserman, L., Trilling, D., Ainsworth, D., Anthony, S., Baker, R., Bardecker, J., Bean, J., Bock, S., Chase, S., Dean, B., Frei, C., George, T., Gill, H., Gimple, H., Givot, R., Hopfe, S., Cota, J., Kehrli, M., King, R., Haley, S., Lara, C., Lund, N., Mattes, M., McCandless, K., McCrystal, D., McRae, J., Rodriguez Melgarejo, L., Mendoza, P., Miller, A., Norfolk, I., Palmquist, B., Reaves, R., Rivard, M., von Schalscha, M., Schar, R., Stoffel, T., Swanson, D., Thompson, D., Wise, J., Woods, L., Yang, Y., 2021, PSJ, 2, 22, The Sizes and Albedos of Centaurs 2014 YY49 and 2013 NL24 from Stellar Occultation Measurements by RECON
    In 2019, the Research and Education Collaborative Occultation Network (RECON) obtained multiple-chord occultation measurements of two Centaur objects: 2014 YY49 on 2019 January 28 and 2013 NL24 on 2019 September 4. RECON is a citizen-science telescope network designed to observe high-uncertainty occultations by outer solar system objects. Adopting circular models for the object profiles, we derive a radius $r={16}_{-1}^{+2}$ km and a geometric albedo ${p}_{V}={0.13}_{-0.024}^{+0.015}$ for 2014 YY49 and a radius $r={66}_{-5}^{+5}$ km and a geometric albedo ${p}_{V}={0.045}_{-0.008}^{+0.006}$ for 2013 NL24. To the precision of these measurements, no atmosphere or rings are detected for either object. The two objects measured here are among the smallest distant objects measured with the stellar occultation technique. In addition to these geometric constraints, the occultation measurements provide astrometric constraints for these two Centaurs at a higher precision than has been feasible by direct imaging. To supplement the occultation results, we also present an analysis of color photometry from the Pan-STARRS surveys to constrain the rotational light curve amplitudes and spectral colors of these two Centaurs. We recommend that future work focus on photometry to more deliberately constrain the objects' colors and light curve amplitudes and on follow-on occultation efforts informed by this astrometry.
  125. Ye, Q., Knight, M., Kelley, M., Moskovitz, N., Gustafsson, A., Schleicher, D., 2021, PSJ, 2, 23, A Deep Search for Emission from "Rock Comet" (3200) Phaethon at 1 au
    We present a deep imaging and spectroscopic search for emission from (3200) Phaethon, a large near-Earth asteroid that appears to be the parent of the strong Geminid meteoroid stream, using the 4.3 m Lowell Discovery Telescope. Observations were conducted on 2017 December 14-18 when Phaethon passed only 0.07 au from the Earth. We determine the 3 upper level of dust and CN production rates to be 0.007-0.2 kg s-1 and 2.3 1022 molecules s-1 through narrowband imaging. A search in broadband images taken through the SDSS r' filter shows no 100 m class fragments in Phaethon's vicinity. A deeper but star-contaminated search also shows no sign of fragments down to 15 m. Optical spectroscopy of Phaethon and comet C/2017 O1 (ASASSN) as a comparison confirms the absence of cometary emission lines from Phaethon and yields 3 upper levels of CN, C2, and C3 of 1024-1025 molecules s-1, 2 orders of magnitude higher than the CN constraint placed by narrowband imaging, due to the much narrower on-sky aperture of the spectrographic slit. We show that narrowband imaging could provide an efficient way to look for weak gas emission from near-extinct bodies near the Earth, though these observations require careful interpretation. Assuming Phaethon's behavior is unchanged, our analysis shows that the DESTINY+ mission, currently planning to explore Phaethon in 2026, may not be able to directly detect a gas coma.
  126. Nesvorny, D., Li, R., Simon, J., Youdin, A., Richardson, D., Marschall, R., Grundy, W., 2021, PSJ, 2, 27, Binary Planetesimal Formation from Gravitationally Collapsing Pebble Clouds
    Planetesimals are compact astrophysical objects roughly 1-1000 km in size, massive enough to be held together by gravity. They can grow by accreting material to become full-size planets. Planetesimals themselves are thought to form by complex physical processes from small grains in protoplanetary disks. The streaming instability (SI) model states that millimeter/centimeter-sized particles (pebbles) are aerodynamically collected into self-gravitating clouds that then directly collapse into planetesimals. Here we analyze ATHENA simulations of the SI to characterize the initial properties (e.g., rotation) of pebble clouds. Their gravitational collapse is followed with the PKDGRAV N-body code, which has been modified to realistically account for pebble collisions. We find that pebble clouds rapidly collapse into short-lived disk structures from which planetesimals form. The planetesimal properties depend on the cloud's scaled angular momentum, $l=L/({{MR}}_{{\rm{H}}}^{2}{\rm{\Omega }})$ , where L and M are the angular momentum and mass, RH is the Hill radius, and is the orbital frequency. Low-l pebble clouds produce tight (or contact) binaries and single planetesimals. Compact high-l clouds give birth to binary planetesimals with attributes that closely resemble the equally sized binaries found in the Kuiper Belt. Significantly, the SI-triggered gravitational collapse can explain the angular momentum distribution of known equally sized binariesa result pending verification from studies with improved resolution. About 10% of collapse simulations produce hierarchical systems with two or more large moons. These systems should be found in the Kuiper Belt when observations reach the threshold sensitivity.
  127. Porredon, A., Crocce, M., Fosalba, P., Elvin-Poole, J., Carnero Rosell, A., Cawthon, R., Eifler, T., Fang, X., Ferrero, I., Krause, E., MacCrann, N., Weaverdyck, N., Abbott, T., Aguena, M., Allam, S., Amon, A., Avila, S., Bacon, D., Bertin, E., Bhargava, S., Bridle, S., Brooks, D., Carrasco Kind, M., Carretero, J., Castander, F., Choi, A., Costanzi, M., da Costa, L., Pereira, M., De Vicente, J., Desai, S., Diehl, H., Doel, P., Drlica-Wagner, A., Eckert, K., Ferte, A., Flaugher, B., Frieman, J., Garcia-Bellido, J., Gaztanaga, E., Gerdes, D., Giannantonio, T., Gruen, D., Gruendl, R., Gschwend, J., Gutierrez, G., Hartley, W., Hinton, S., Hollowood, D., Honscheid, K., Hoyle, B., James, D., Jarvis, M., Kuehn, K., Kuropatkin, N., Maia, M., Marshall, J., Menanteau, F., Miquel, R., Morgan, R., Palmese, A., Pandey, S., Paz-Chinchon, F., Plazas, A., Rodriguez-Monroy, M., Roodman, A., Samuroff, S., Sanchez, E., Scarpine, V., Serrano, S., Sevilla-Noarbe, I., Smith, M., Soares-Santos, M., Suchyta, E., Swanson, M., Tarle, G., To, C., Varga, T., Weller, J., Wilkinson, R., DES Collaboration, 2021, PhRvD, 103, 043503, Dark Energy Survey Year 3 results: Optimizing the lens sample in a combined galaxy clustering and galaxy-galaxy lensing analysis
    We investigate potential gains in cosmological constraints from the combination of galaxy clustering and galaxy-galaxy lensing by optimizing the lens galaxy sample selection using information from Dark Energy Survey (DES) Year 3 data and assuming the DES Year 1 METACALIBRATION sample for the sources. We explore easily reproducible selections based on magnitude cuts in i -band as a function of (photometric) redshift, zphot, and benchmark the potential gains against those using the well-established REDMAGIC [E. Rozo et al., Mon. Not. R. Astron. Soc. 461, 1431 (2016), 10.1093/mnras/stw1281] sample. We focus on the balance between density and photometric redshift accuracy, while marginalizing over a realistic set of cosmological and systematic parameters. Our optimal selection, the MAGLIM sample, satisfies i <4 zphot+18 and has 30 % wider redshift distributions but 3.5 times more galaxies than REDMAGIC. Assuming a w CDM model (i.e. with a free parameter for the dark energy equation of state) and equivalent scale cuts to mitigate nonlinear effects, this leads to 40% increase in the figure of merit for the pair combinations of m, w , and 8, and gains of 16% in 8, 10% in m, and 12% in w . Similarly, in CDM , we find an improvement of 19% and 27% on 8 and m, respectively. We also explore flux-limited samples with a flat magnitude cut finding that the optimal selection, i <22.2 , has 7 times more galaxies and 20 % wider redshift distributions compared to MAGLIM, but slightly worse constraints. We show that our results are robust with respect to the assumed galaxy bias and photometric redshift uncertainties with only moderate further gains from increased number of tomographic bins or the inclusion of bin cross-correlations, except in the case of the flux-limited sample, for which these gains are more significant.
  128. Hunter, D., Elmegreen, B., Goldberger, E., Taylor, H., Ermakov, A., Herrmann, K., Oh, S., Malko, B., Barandi, B., Jundt, R., 2021, AJ, 161, 71, Relationships between the Stellar, Gaseous, and Star Formation Disks in LITTLE THINGS Dwarf Irregular Galaxies: Indirect Evidence for Substantial Fractions of Dark Molecular Gas
    The stellar, gaseous and young stellar disks in the LITTLE THINGS sample of nearby dwarf irregular galaxies are fitted with functions to search for correlations between the parameters. We find that the H I radial profiles are generally flatter in the center and fall faster in the outer regions than the V-band profiles, while young stars are more centrally concentrated, especially if the H I is more centrally flat. This pattern suggests that the H I is turning into molecules in the center, and the molecular clouds are forming stars and FUV. A model that assumes the molecular surface density is proportional to the total gas surface density to a power of 1.5 or 2, in analogy with the Kennicutt-Schmidt relation, reproduces the relationship between the ratio of the visible to the H I scale length and the H I Sersic index. The molecular fraction is estimated as a function of radius for each galaxy by converting the FUV to a molecular surface density using conventional calibrations. The average molecular fraction inside 3RD is 23% 17%. However, the break in the stellar surface brightness profile has no unified tracer related to star formation.
  129. Massey, P., Neugent, K., Levesque, E., Drout, M., Courteau, S., 2021, AJ, 161, 79, The Red Supergiant Content of M31 and M33
    We identify red supergiants (RSGs) in our spiral neighbors M31 and M33 using near-IR (NIR) photometry complete to a luminosity limit of $\mathrm{log}L/{L}_{\odot }=4.0$ . Our archival survey data cover 5 deg2 of M31, and 3 deg2 for M33, and are likely spatially complete for these massive stars. Gaia is used to remove foreground stars, after which the RSGs can be separated from asymptotic giant branch (AGB) stars in the color-magnitude diagram. The photometry is used to derive effective temperatures and bolometric luminosities via MARCS stellar atmosphere models. The resulting H-R diagrams show superb agreement with the evolutionary tracks of the Geneva evolutionary group. Our census includes 6400 RSGs in M31 and 2850 RSGs in M33 within their Holmberg radii; by contrast, only a few hundred RSGs are known so far in the Milky Way. Our catalog serves as the basis for a study of the RSG binary frequency being published separately, as well as future studies relating to the evolution of massive stars. Here we use the matches between the NIR-selected RSGs and their optical counterparts to show that the apparent similarity in the reddening of OB stars in M31 and M33 is the result of Malmquist bias; the average extinction in M31 is likely higher than that of M33. As expected, the distribution of RSGs follows that of the spiral arms, while the much older AGB population is more uniformly spread across each galaxy's disk.
  130. Tanoglidis, D., Drlica-Wagner, A., Wei, K., Li, T., Sanchez, J., Zhang, Y., Peter, A., Feldmeier-Krause, A., Prat, J., Casey, K., Palmese, A., Sanchez, C., DeRose, J., Conselice, C., Gagnon, L., Abbott, T., Aguena, M., Allam, S., Avila, S., Bechtol, K., Bertin, E., Bhargava, S., Brooks, D., Burke, D., Rosell, A., Kind, M., Carretero, J., Chang, C., Costanzi, M., da Costa, L., De Vicente, J., Desai, S., Diehl, H., Doel, P., Eifler, T., Everett, S., Evrard, A., Flaugher, B., Frieman, J., Garcia-Bellido, J., Gerdes, D., Gruendl, R., Gschwend, J., Gutierrez, G., Hartley, W., Hollowood, D., Huterer, D., James, D., Krause, E., Kuehn, K., Kuropatkin, N., Maia, M., March, M., Marshall, J., Menanteau, F., Miquel, R., Ogando, R., Paz-Chinchon, F., Romer, A., Roodman, A., Sanchez, E., Scarpine, V., Serrano, S., Sevilla-Noarbe, I., Smith, M., Suchyta, E., Tarle, G., Thomas, D., Tucker, D., Walker, A., DES Collaboration, 2021, ApJS, 252, 18, Shadows in the Dark: Low-surface-brightness Galaxies Discovered in the Dark Energy Survey
    We present a catalog of 23,790 extended low-surface-brightness galaxies (LSBGs) identified in $\sim 5000\,{\deg }^{2}$ from the first three years of imaging data from the Dark Energy Survey (DES). Based on a single-component Sersic model fit, we define extended LSBGs as galaxies with g-band effective radii ${R}_{\mathrm{eff}}(g)\gt 2\buildrel{\prime\prime}\over{.} 5$ and mean surface brightness ${\bar{\mu }}_{\mathrm{eff}}(g)\gt 24.2\,\mathrm{mag}\,{\mathrm{arcsec}}^{-2}$ . We find that the distribution of LSBGs is strongly bimodal in (g - r) versus (g - i) color space. We divide our sample into red (g - i 0.60) and blue (g - i < 0.60) galaxies and study the properties of the two populations. Redder LSBGs are more clustered than their blue counterparts and are correlated with the distribution of nearby (z < 0.10) bright galaxies. Red LSBGs constitute 33% of our LSBG sample, and $\sim 30 \% $ of these are located within 1 of low-redshift galaxy groups and clusters (compared to 8% of the blue LSBGs). For nine of the most prominent galaxy groups and clusters, we calculate the physical properties of associated LSBGs assuming a redshift derived from the host system. In these systems, we identify 41 objects that can be classified as ultradiffuse galaxies, defined as LSBGs with projected physical effective radii ${R}_{\mathrm{eff}}\gt 1.5\,\mathrm{kpc}$ and central surface brightness ${\mu }_{0}(g)\gt 24.0\,\mathrm{mag}\,{\mathrm{arcsec}}^{-2}$ . The wide-area sample of LSBGs in DES can be used to test the role of environment on models of LSBG formation and evolution.
  131. Farnham, T., Knight, M., Schleicher, D., Feaga, L., Bodewits, D., Skiff, B., Schindler, J., 2021, PSJ, 2, 7, Narrowband Observations of Comet 46P/Wirtanen during Its Exceptional Apparition of 2018/19. I. Apparent Rotation Period and Outbursts
    We obtained broad- and narrowband images of the hyperactive comet 46P/Wirtanen on 33 nights during its 2018/2019 apparition, when the comet made a historic close approach to the Earth. With our extensive coverage, we investigated the temporal behavior of the comet on both seasonal and rotational timescales. We used CN observations to explore the coma morphology, revealing that there are two primary active areas that produce spiral structures. The direction of rotation of these structures changes from pre- to postperihelion, indicating that the Earth crossed the comet's equatorial plane sometime around perihelion. We also used the CN images to create photometric light curves that consistently show two peaks in the activity, confirming the two source regions. We measured the nucleus's apparent rotation period at a number of epochs using both the morphology and the light curves. These results all show that the rotation period is continuously changing throughout our observation window, increasing from 8.98 hr in early November to 9.14 hr around perihelion and then decreasing again to 8.94 hr in February. Although the geometry changes rapidly around perihelion, the period changes cannot be primarily due to synodic effects. The repetition of structures in the coma, both within a night and from night to night, strongly suggests that the nucleus is in a near-simple rotation state. We also detected two outbursts, one on December 12 and the other on January 28. Using the apparent velocities of the ejecta in these events, 68 5 and 162 15 m s-1, respectively, we derived start times of 2018 December 12 at 00:13 UT 7 minutes and 2019 January 27 at 20:01 UT 30 minutes.
  132. Fernandez-Valenzuela, E., Pinilla-Alonso, N., Stansberry, J., Emery, J., Perkins, W., Van Laerhoven, C., Gladman, B., Fraser, W., Cruikshank, D., Lellouch, E., Muller, T., Grundy, W., Trilling, D., Fernandez, Y., Dalle Ore, C., 2021, PSJ, 2, 10, Compositional Study of Trans-Neptunian Objects at > 2.2 m
    Using data from the Infrared Array Camera on the Spitzer Space Telescope, we present photometric observations of a sample of 100 trans-Neptunian objects (TNOs) beyond 2.2 m. These observations, collected with two broadband filters centered at 3.6 and 4.5 m, were done in order to study the surface composition of TNOs, which are too faint to obtain spectroscopic measurements. With this aim, we have developed a method for the identification of different materials that are found on the surfaces of TNOs. In our sample, we detected objects with colors that are consistent with the presence of small amounts of water, and we were able to distinguish between surfaces that are predominantly composed of complex organics and amorphous silicates. We found that 86% of our sample have characteristics that are consistent with a certain amount of water ice, and the most common composition (73% of the objects) is a mixture of water ice, amorphous silicates, and complex organics. Twenty-three percent of our sample may include other ices, such as carbon monoxide, carbon dioxide, methane, or methanol. Additionally, only small objects seem to have surfaces dominated by silicates. This method is a unique tool for the identification of complex organics and to obtain the surface composition of extremely faint objects. Furthermore, this method will be beneficial when using the James Webb Space Telescope for differentiating groups within the trans-Neptunian population.
  133. Lindor, B., Hartman, J., Bakos, G., Bhatti, W., Csubry, Z., Penev, K., Bieryla, A., Latham, D., Torres, G., Buchhave, L., de Val-Borro, M., Howard, A., Isaacson, H., Fulton, B., Boisse, I., Santerne, A., Hebrard, G., Kovacs, T., Huang, C., Dembicky, J., Falco, E., Everett, M., Horch, E., Lazar, J., Papp, I., Sari, P., 2021, AJ, 161, 64, HAT-P-68b: A Transiting Hot Jupiter around a K5 Dwarf Star
    We report the discovery by the ground-based Hungarian-made Automated Telescope Network (HATNet) survey of the transiting exoplanet HAT-P-68b, which has a mass of 0.724 0.043 MJ, and radius of 1.072 0.012 RJ. The planet is in a circular P = 2.2984 day orbit around a moderately bright V = 13.937 0.030 magnitude K-dwarf star of mass ${0.673}_{-0.014}^{+0.020}$ M, and radius 0.6726 0.0069 R. The planetary nature of this system is confirmed through follow-up transit photometry obtained with the Fred L. Whipple Observatory (FLWO) 1.2 m telescope, high-precision radial velocities measured using Keck I/High Resolution Echelle Spectrometer (HIRES), FLWO 1.5 m/Tillinghast Reflector Echelle Spectrograph (TRES), and Observatoire de Haute-Provence (OHP) 1.9 m/Sophie, and high-spatial-resolution speckle imaging from WIYN 3.5 m/DSSI. HAT-P-68 is at an ecliptic latitude of +3 and outside the field of view of both the NASA Transiting Exoplanet Survey Satellite primary mission and the K2 mission. The large transit depth of 0.036 mag (r band) makes HAT-P-68b a promising target for atmospheric characterization via transmission spectroscopy. * Based on observations obtained with the Hungarian-made Automated Telescope Network. Based in part on observations made with the Keck I telescope at Maunakea Observatory, Hawaii (Keck time awarded through NASA programs N133Hr and N169Hr). Based in part on observations obtained with the Tillinghast Reflector 1.5 m telescope and the 1.2 m telescope, both operated by the Smithsonian Astrophysical Observatory at the Fred Lawrence Whipple Observatory in Arizona. Based on radial velocities obtained with the Sophie spectrograph mounted on the 1.93 m telescope at Observatoire de Haute-Provence.
  134. Holler, B., Grundy, W., Buie, M., Noll, K., 2021, Icar, 355, 114130, The Eris/Dysnomia system I: The orbit of Dysnomia
    We present new results on the Eris/Dysnomia system including analysis of new images from the WFC3 instrument on the Hubble Space Telescope (HST). Seven HST orbits were awarded to program 15171 in January and February 2018, with the intervals between observations selected to sample Dysnomia over a full orbital period. Using relative astrometry of Eris and Dysnomia, we computed a best-fit Keplerian orbit for Dysnomia. Based on the Keplerian fit, we find an orbital period of 15.7858990.000050 days, which is in good agreement with recent work. We report a non-zero eccentricity of 0.0062 at the 6.2- level, despite an estimated eccentricity damping timescale of 17 Myr. Considering the volumes of both Eris and Dysnomia, the new system density was calculated to be 2.430.05 g cm-3, a decrease of ~4% from the previous value of 2.520.05 g cm-3. The new astrometric measurements were high enough precision to break the degeneracy of the orbit pole orientation, and indicate that Dysnomia orbits in a prograde manner. The obliquity of Dysnomia's orbit pole with respect to the plane of Eris' heliocentric orbit was calculated to be 78.290.65 and is in agreement with previous work; the next mutual events season will occur in 2239. The Keplerian orbit fit to all the data considered in this investigation can be excluded at the 6.3- level, but identifying the cause of the deviation was outside the scope of this work.
  135. Vielzeuf, P., Kovacs, A., Demirbozan, U., Fosalba, P., Baxter, E., Hamaus, N., Huterer, D., Miquel, R., Nadathur, S., Pollina, G., Sanchez, C., Whiteway, L., Abbott, T., Allam, S., Annis, J., Avila, S., Brooks, D., Burke, D., Carnero Rosell, A., Carrasco Kind, M., Carretero, J., Cawthon, R., Costanzi, M., da Costa, L., De Vicente, J., Desai, S., Diehl, H., Doel, P., Eifler, T., Everett, S., Flaugher, B., Frieman, J., Garcia-Bellido, J., Gaztanaga, E., Gerdes, D., Gruen, D., Gruendl, R., Gschwend, J., Gutierrez, G., Hartley, W., Hollowood, D., Honscheid, K., James, D., Kuehn, K., Kuropatkin, N., Lahav, O., Lima, M., Maia, M., March, M., Marshall, J., Melchior, P., Menanteau, F., Palmese, A., Paz-Chinchon, F., Plazas, A., Sanchez, E., Scarpine, V., Serrano, S., Sevilla-Noarbe, I., Smith, M., Suchyta, E., Tarle, G., Thomas, D., Weller, J., Zuntz, J., Zuntz, J., DES Collaboration, 2021, MNRAS, 500, 464, Dark Energy Survey Year 1 results: the lensing imprint of cosmic voids on the cosmic microwave background
    Cosmic voids gravitationally lens the cosmic microwave background (CMB) radiation, resulting in a distinct imprint on degree scales. We use the simulated CMB lensing convergence map from the Marenostrum Institut de Ciencias de l'Espai (MICE) N-body simulation to calibrate our detection strategy for a given void definition and galaxy tracer density. We then identify cosmic voids in Dark Energy Survey (DES) Year 1 data and stack the Planck 2015 lensing convergence map on their locations, probing the consistency of simulated and observed void lensing signals. When fixing the shape of the stacked convergence profile to that calibrated from simulations, we find imprints at the 3 significance level for various analysis choices. The best measurement strategies based on the MICE calibration process yield S/N 4 for DES Y1, and the best-fitting amplitude recovered from the data is consistent with expectations from MICE (A 1). Given these results as well as the agreement between them and N-body simulations, we conclude that the previously reported excess integrated Sachs-Wolfe (ISW) signal associated with cosmic voids in DES Y1 has no counterpart in the Planck CMB lensing map.
  136. Muir, J., Baxter, E., Miranda, V., Doux, C., Ferte, A., Leonard, C., Huterer, D., Jain, B., Lemos, P., Raveri, M., Nadathur, S., Campos, A., Chen, A., Dodelson, S., Elvin-Poole, J., Lee, S., Secco, L., Troxel, M., Weaverdyck, N., Zuntz, J., Brout, D., Choi, A., Crocce, M., Davis, T., Gruen, D., Krause, E., Lidman, C., MacCrann, N., Moller, A., Prat, J., Ross, A., Sako, M., Samuroff, S., Sanchez, C., Scolnic, D., Zhang, B., Abbott, T., Aguena, M., Allam, S., Annis, J., Avila, S., Bacon, D., Bertin, E., Bhargava, S., Bridle, S., Brooks, D., Burke, D., Carnero Rosell, A., Carrasco Kind, M., Carretero, J., Cawthon, R., Costanzi, M., da Costa, L., Pereira, M., Desai, S., Diehl, H., Dietrich, J., Doel, P., Estrada, J., Everett, S., Evrard, A., Ferrero, I., Flaugher, B., Frieman, J., Garcia-Bellido, J., Giannantonio, T., Gruendl, R., Gschwend, J., Gutierrez, G., Hinton, S., Hollowood, D., Honscheid, K., Hoyle, B., James, D., Jeltema, T., Kuehn, K., Kuropatkin, N., Lahav, O., Lima, M., Maia, M., Menanteau, F., Miquel, R., Morgan, R., Myles, J., Palmese, A., Paz-Chinchon, F., Plazas, A., Romer, A., Roodman, A., Sanchez, E., Scarpine, V., Serrano, S., Sevilla-Noarbe, I., Smith, M., Suchyta, E., Swanson, M., Tarle, G., Thomas, D., To, C., Tucker, D., Varga, T., Weller, J., Wilkinson, R., DES Collaboration, 2021, PhRvD, 103, 023528, DES Y1 results: Splitting growth and geometry to test CDM
    We analyze Dark Energy Survey (DES) data to constrain a cosmological model where a subset of parametersfocusing on mare split into versions associated with structure growth (e.g., mgrow) and expansion history (e.g., mgeo). Once the parameters have been specified for the CDM cosmological model, which includes general relativity as a theory of gravity, it uniquely predicts the evolution of both geometry (distances) and the growth of structure over cosmic time. Any inconsistency between measurements of geometry and growth could therefore indicate a breakdown of that model. Our growth-geometry split approach therefore serves both as a (largely) model-independent test for beyond- CDM physics, and as a means to characterize how DES observables provide cosmological information. We analyze the same multiprobe DES data as [Phys. Rev. Lett. 122, 171301 (2019), 10.1103/PhysRevLett.122.171301] : DES Year 1 (Y1) galaxy clustering and weak lensing, which are sensitive to both growth and geometry, as well as Y1 BAO and Y3 supernovae, which probe geometry. We additionally include external geometric information from BOSS DR12 BAO and a compressed Planck 2015 likelihood, and external growth information from BOSS DR12 RSD. We find no significant disagreement with mgrow=mgeo. When DES and external data are analyzed separately, degeneracies with neutrino mass and intrinsic alignments limit our ability to measure mgrow, but combining DES with external data allows us to constrain both growth and geometric quantities. We also consider a parametrization where we split both m and w , but find that even our most constraining data combination is unable to separately constrain mgrow and wgrow. Relative to CDM , splitting growth and geometry weakens bounds on 8 but does not alter constraints on h .
  137. Lauer, T., Postman, M., Weaver, H., Spencer, J., Stern, S., Buie, M., Durda, D., Lisse, C., Poppe, A., Binzel, R., Britt, D., Buratti, B., Cheng, A., Grundy, W., Horanyi, M., Kavelaars, J., Linscott, I., McKinnon, W., Moore, J., Nunez, J., Olkin, C., Parker, J., Porter, S., Reuter, D., Robbins, S., Schenk, P., Showalter, M., Singer, K., Verbiscer, A., Young, L., 2021, ApJ, 906, 77, New Horizons Observations of the Cosmic Optical Background
    We used existing data from the New Horizons Long-range Reconnaissance Imager (LORRI) to measure the optical-band (0.4 0.9 m) sky brightness within seven high-Galactic latitude fields. The average raw level measured while New Horizons was 42-45 au from the Sun is 33.2 0.5 nW m-2 sr-1. This is 10 as dark as the darkest sky accessible to the Hubble Space Telescope, highlighting the utility of New Horizons for detecting the cosmic optical background (COB). Isolating the COB contribution to the raw total required subtracting scattered light from bright stars and galaxies, faint stars below the photometric detection limit within the fields, and diffuse Milky Way light scattered by infrared cirrus. We removed newly identified residual zodiacal light from the IRIS 100 m all-sky maps to generate two different estimates for the diffuse Galactic light. Using these yielded a highly significant detection of the COB in the range 15.9 4.2 (1.8 stat., 3.7 sys.) nW m-2 sr-1 to 18.7 3.8 (1.8 stat., 3.3 sys.) nW m-2 sr-1 at the LORRI pivot wavelength of 0.608 m. Subtraction of the integrated light of galaxies fainter than the photometric detection limit from the total COB level left a diffuse flux component of unknown origin in the range 8.8 4.9 (1.8 stat., 4.5 sys.) nW m-2 sr-1 to 11.9 4.6 (1.8 stat., 4.2 sys.) nW m-2 sr-1. Explaining it with undetected galaxies requires the assumption that the galaxy count faint-end slope steepens markedly at V > 24 or that existing surveys are missing half the galaxies with V < 30.
  138. Biddle, L., Llama, J., Cameron, A., Prato, L., Jardine, M., Johns-Krull, C., 2021, ApJ, 906, 113, Amplitude Modulation of Short-timescale Hot Spot Variability
    Variability of Classical T Tauri stars (CTTS) occurs over a vast range of timescales. CTTS in particular are subject to variability caused by accretion shocks, which can occur stochastically, periodically, or quasi-periodically on timescales over a few days. The detectability of young planets within these systems is likely hampered by activity; therefore, it is essential that we understand the origin of young star variability over a range of timescales to help disentangle stellar activity from signatures of planetary origin. We present an analysis of the stochastic small-amplitude photometric variability in the K2 lightcurve of CI Tau occurring on timescales of 1 day. We find the amplitude of this variability exhibits the same periodic signatures as detected in the large-amplitude variability, indicating that the physical mechanism modulating these brightness features is the same. The periods detected are also in agreement with the rotation period of the star (6.6 days) and the orbital period of the planet (9.0 days) known to drive pulsed accretion onto the star.
  139. Christou, A., Borisov, G., Dell'Oro, A., Cellino, A., Devogele, M., 2021, Icar, 354, 113994, Composition and origin of L5 Trojan asteroids of Mars: Insights from spectroscopy
    We investigate the mineralogical makeup of L5 Martian Trojan asteroids via reflectance spectroscopy, paying special attention to (101429) 1998 VF31, the only L5 Trojan that does not belong to the Eureka family (Christou, 2013). We find that this asteroid most likely belongs to the Bus-Demeo S-complex, in agreement with Rivkin et al. (2007). We compare it with a variety of solar system bodies and obtain good spectral matches with Sq- or S-type asteroids, with spectra of the lunar surface and of Martian and lunar meteorites. Mixture fitting to spectral endmembers suggests a surface abundance of Mg-rich orthopyroxene and iron metal or, alternatively, a combination of plagioclase and metal with a small amount of Mg-poor orthopyroxene. The metallic component may be part of the intrinsic mineral makeup of the asteroid or an indication of extreme space weathering. In light of our findings, we discuss a number of origin scenarios for (101429). The asteroid could be genetically related to iron-rich primitive achondrite meteorites (Rivkin et al 2007), may have originated as impact ejecta from Mars - a scenario proposed recently for the Eureka family asteroids (Polishook et al., 2017) - or could represent a relic fragment of the Moon's original solid crust, a possibility raised by the asteroid's close spectral similarity to areas of the lunar surface. If, on the other hand, (101429) is a relatively recent addition to the Martian Trojan clouds (Christou et al., 2020), its origin is probably traced to high-inclination asteroid families in the Inner Main Belt. For the olivine-dominated Eureka family, we find that the two smaller asteroids in our sample are more spectrally similar to one another than to (5261) Eureka, the largest family member. Spectral profiles of these three asteroids are closely similar shortward of ~0.7 m but diverge at longer wavelengths. For the two smaller asteroids in particular, we find the spectra are virtually identical in the visible region and up to 0.8 m. We attribute spectral differences in the near-IR region to differences in either: degree of space weathering, olivine chemical composition and/or regolith grain size.
  140. Cabot, S., Roettenbacher, R., Henry, G., Zhao, L., Harmon, R., Fischer, D., Brewer, J., Llama, J., Petersburg, R., Szymkowiak, A., 2021, AJ, 161, 26, EXPRES. II. Searching for Planets around Active Stars: A Case Study of HD 101501
    By controlling instrumental errors to below 10 cm s-1, the EXtreme PREcision Spectrograph (EXPRES) allows for a more insightful study of photospheric velocities that can mask weak Keplerian signals. Gaussian processes (GP) have become a standard tool for modeling correlated noise in radial velocity data sets. While GPs are constrained and motivated by physical properties of the star, in some cases they are still flexible enough to absorb unresolved Keplerian signals. We apply GP regression to EXPRES radial velocity measurements of the 3.5 Gyr old chromospherically active Sun-like star, HD 101501. We obtain tight constraints on the stellar rotation period and the evolution of spot distributions using 28 seasons of ground-based photometry, as well as recent Transiting Exoplanet Survey Satellite data. Light-curve inversion was carried out on both photometry data sets to reveal the spot distribution and spot evolution timescales on the star. We find that the >5 m s-1 rms radial velocity variations in HD 101501 are well modeled with a GP stellar activity model without planets, yielding a residual rms scatter of 45 cm s-1. We carry out simulations, injecting and recovering signals with the GP framework, to demonstrate that high-cadence observations are required to use GPs most efficiently to detect low-mass planets around active stars like HD 101501. Sparse sampling prevents GPs from learning the correlated noise structure and can allow it to absorb prospective Keplerian signals. We quantify the moderate to high-cadence monitoring that provides the necessary information to disentangle photospheric features using GPs and to detect planets around active stars.
  141. Bryson, S., Kunimoto, M., Kopparapu, R., Coughlin, J., Borucki, W., Koch, D., Aguirre, V., Allen, C., Barentsen, G., Batalha, N., Berger, T., Boss, A., Buchhave, L., Burke, C., Caldwell, D., Campbell, J., Catanzarite, J., Chandrasekaran, H., Chaplin, W., Christiansen, J., Christensen-Dalsgaard, J., Ciardi, D., Clarke, B., Cochran, W., Dotson, J., Doyle, L., Duarte, E., Dunham, E., Dupree, A., Endl, M., Fanson, J., Ford, E., Fujieh, M., Gautier, T., Geary, J., Gilliland, R., Girouard, F., Gould, A., Haas, M., Henze, C., Holman, M., Howard, A., Howell, S., Huber, D., Hunter, R., Jenkins, J., Kjeldsen, H., Kolodziejczak, J., Larson, K., Latham, D., Li, J., Mathur, S., Meibom, S., Middour, C., Morris, R., Morton, T., Mullally, F., Mullally, S., Pletcher, D., Prsa, A., Quinn, S., Quintana, E., Ragozzine, D., Ramirez, S., Sanderfer, D., Sasselov, D., Seader, S., Shabram, M., Shporer, A., Smith, J., Steffen, J., Still, M., Torres, G., Troeltzsch, J., Twicken, J., Uddin, A., Van Cleve, J., Voss, J., Weiss, L., Welsh, W., Wohler, B., Zamudio, K., 2021, AJ, 161, 36, The Occurrence of Rocky Habitable-zone Planets around Solar-like Stars from Kepler Data
    We present the occurrence rates for rocky planets in the habitable zones (HZs) of main-sequence dwarf stars based on the Kepler DR25 planet candidate catalog and Gaia-based stellar properties. We provide the first analysis in terms of star-dependent instellation flux, which allows us to track HZ planets. We define as the HZ occurrence of planets with radii between 0.5 and 1.5 R orbiting stars with effective temperatures between 4800 and 6300 K. We find that for the conservative HZ is between ${0.37}_{-0.21}^{+0.48}$ (errors reflect 68% credible intervals) and ${0.60}_{-0.36}^{+0.90}$ planets per star, while the optimistic HZ occurrence is between ${0.58}_{-0.33}^{+0.73}$ and ${0.88}_{-0.51}^{+1.28}$ planets per star. These bounds reflect two extreme assumptions about the extrapolation of completeness beyond orbital periods where DR25 completeness data are available. The large uncertainties are due to the small number of detected small HZ planets. We find similar occurrence rates between using Poisson likelihood Bayesian analysis and using Approximate Bayesian Computation. Our results are corrected for catalog completeness and reliability. Both completeness and the planet occurrence rate are dependent on stellar effective temperature. We also present occurrence rates for various stellar populations and planet size ranges. We estimate with 95% confidence that, on average, the nearest HZ planet around G and K dwarfs is 6 pc away and there are 4 HZ rocky planets around G and K dwarfs within 10 pc of the Sun.
  142. Colton, N., Horch, E., Everett, M., Howell, S., Davidson, J., Baptista, B., Casetti-Dinescu, D., 2021, AJ, 161, 21, Identifying Bound Stellar Companions to Kepler Exoplanet Host Stars Using Speckle Imaging
    The Kepler mission and subsequent ground-based follow-up observations have revealed a number of exoplanet host stars with nearby stellar companions. This study presents speckle observations of 57 Kepler objects of interest (KOIs) that are also double stars, each observed over a 3-8 yr period, which has allowed us to track their relative motions with high precision. Measuring the position angle and separation of the companion with respect to the primary can help determine if the pair exhibits common proper motion, indicating it is likely to be a bound binary system. We report on the motions of 34 KOIs that have close stellar companions, three of which are triple stars, for a total of 37 companions studied. Eighteen of the 34 systems are confirmed exoplanet hosts, including one triple star, while four other systems have been subsequently judged to be false positives and twelve are yet to be confirmed as planet hosts. We find that 21 are most likely to be common proper motion pairs, 4 are line-of-sight companions, and 12 are of an uncertain disposition at present. The fraction of the confirmed exoplanet host systems that are common proper motion pairs is approximately 86% in this sample. In this subsample, the planets are exclusively found with periods of less than 110 days, so that in all cases the stellar companion is found at a much larger separation from the planet host star than the planet itself. A preliminary period-radius relation for the confirmed planets in our sample suggests no obvious differences at this stage with the full sample of known exoplanets.
  143. De Propris, R., West, M., Andrade-Santos, F., Ragone-Figueroa, C., Rasia, E., Forman, W., Jones, C., Kipper, R., Borgani, S., Lambas, D., Romashkova, E., Patra, K., 2021, MNRAS, 500, 310, Brightest cluster galaxies: the centre can(not?) hold
    We explore the persistence of the alignment of brightest cluster galaxies (BCGs) with their local environment. We find that a significant fraction of BCGs do not coincide with the centroid of the X-ray gas distribution and/or show peculiar velocities (they are not at rest with respect to the cluster mean). Despite this, we find that BCGs are generally aligned with the cluster mass distribution even when they have significant offsets from the X-ray centre and significant peculiar velocities. The large offsets are not consistent with simple theoretical models. To account for these observations BCGs must undergo mergers preferentially along their major axis, the main infall direction. Such BCGs may be oscillating within the cluster potential after having been displaced by mergers or collisions, or the dark matter halo itself may not yet be relaxed.
  144. Henghes, B., Lahav, O., Gerdes, D., Lin, H., Morgan, R., Abbott, T., Aguena, M., Allam, S., Annis, J., Avila, S., Bertin, E., Brooks, D., Burke, D., Carnero Rosell, A., Carrasco Kind, M., Carretero, J., Conselice, C., Costanzi, M., da Costa, L., De Vicente, J., Desai, S., Diehl, H., Doel, P., Everett, S., Ferrero, I., Frieman, J., Garcia-Bellido, J., Gaztanaga, E., Gruen, D., Gruendl, R., Gschwend, J., Gutierrez, G., Hartley, W., Hinton, S., Honscheid, K., Hoyle, B., James, D., Kuehn, K., Kuropatkin, N., Marshall, J., Melchior, P., Menanteau, F., Miquel, R., Ogando, R., Palmese, A., Paz-Chinchon, F., Plazas, A., Romer, A., Sanchez, C., Sanchez, E., Scarpine, V., Schubnell, M., Serrano, S., Smith, M., Soares-Santos, M., Suchyta, E., Tarle, G., To, C., Wilkinson, R., DES Collaboration, 2021, PASP, 133, 014501, Machine Learning for Searching the Dark Energy Survey for Trans-Neptunian Objects
    In this paper we investigate how implementing machine learning could improve the efficiency of the search for Trans-Neptunian Objects (TNOs) within Dark Energy Survey (DES) data when used alongside orbit fitting. The discovery of multiple TNOs that appear to show a similarity in their orbital parameters has led to the suggestion that one or more undetected planets, an as yet undiscovered "Planet 9", may be present in the outer solar system. DES is well placed to detect such a planet and has already been used to discover many other TNOs. Here, we perform tests on eight different supervised machine learning algorithms, using a data set consisting of simulated TNOs buried within real DES noise data. We found that the best performing classifier was the Random Forest which, when optimized, performed well at detecting the rare objects. We achieve an area under the receiver operating characteristic (ROC) curve, (AUC) = 0.996 0.001. After optimizing the decision threshold of the Random Forest, we achieve a recall of 0.96 while maintaining a precision of 0.80. Finally, by using the optimized classifier to pre-select objects, we are able to run the orbit-fitting stage of our detection pipeline five times faster.
  145. 144 publications and 3108 citations in 2021.

144 publications and 3108 citations total.

Top