-
A thermodynamic criterion for the formation of Circumplanetary Disks
Authors:
Leonardo Krapp,
Kaitlin M. Kratter,
Andrew N. Youdin,
Pablo Benítez-Llambay,
Frédéric Masset,
Philip J. Armitage
Abstract:
The formation of circumplanetary disks is central to our understanding of giant planet formation, influencing their growth rate during the post-runaway phase and observability while embedded in protoplanetary disks. We use 3D global multifluid radiation hydrodynamics simulations with the FARGO3D code to define the thermodynamic conditions that enable circumplanetary disk formation around Jovian pl…
▽ More
The formation of circumplanetary disks is central to our understanding of giant planet formation, influencing their growth rate during the post-runaway phase and observability while embedded in protoplanetary disks. We use 3D global multifluid radiation hydrodynamics simulations with the FARGO3D code to define the thermodynamic conditions that enable circumplanetary disk formation around Jovian planets on wide orbits. Our simulations include stellar irradiation, viscous heating, static mesh refinement, and active calculation of opacity based on evolving dust fluids. We find a necessary condition for the formation of circumplanetary disks in terms of a mean cooling time: when the cooling time is at least one order of magnitude shorter than the orbital time scale, the specific angular momentum of the gas is nearly Keplerian at scales of $R_{\rm{Hill}}/3$. We show that the inclusion of multifluid dust dynamics favors rotational support because dust settling produces an anisotropic opacity distribution that favors rapid cooling. In all our models with radiation hydrodynamics, specific angular momentum decreases as time evolves in agreement with the formation of an inner isentropic envelope due to compressional heating. The isentropic envelope can extend up to $R_{\rm{Hill}}/3$ and shows negligible rotational support. Thus, our results imply that young gas giant planets may host spherical isentropic envelopes, rather than circumplanetary disks.
△ Less
Submitted 26 February, 2024; v1 submitted 22 February, 2024;
originally announced February 2024.
-
The Santa Barbara Binary-Disk Code Comparison
Authors:
Paul C. Duffell,
Alexander J. Dittmann,
Daniel J. D'Orazio,
Alessia Franchini,
Kaitlin M. Kratter,
Anna B. T. Penzlin,
Enrico Ragusa,
Magdalena Siwek,
Christopher Tiede,
Haiyang Wang,
Jonathan Zrake,
Adam M. Dempsey,
Zoltan Haiman,
Alessandro Lupi,
Michal Pirog,
Geoffrey Ryan
Abstract:
We have performed numerical calculations of a binary interacting with a gas disk, using eleven different numerical methods and a standard binary-disk setup. The goal of this study is to determine whether all codes agree on a numerically converged solution, and to determine the necessary resolution for convergence and the number of binary orbits that must be computed to reach an agreed-upon relaxed…
▽ More
We have performed numerical calculations of a binary interacting with a gas disk, using eleven different numerical methods and a standard binary-disk setup. The goal of this study is to determine whether all codes agree on a numerically converged solution, and to determine the necessary resolution for convergence and the number of binary orbits that must be computed to reach an agreed-upon relaxed state of the binary-disk system. We find that all codes can agree on a converged solution (depending on the diagnostic being measured). The zone spacing required for most codes to reach a converged measurement of the torques applied to the binary by the disk is roughly 1% of the binary separation in the vicinity of the binary components. For our disk model to reach a relaxed state, codes must be run for at least 200 binary orbits, corresponding to about a viscous time for our parameters, $0.2 (a^2 Ω_B /ν)$ binary orbits, where $ν$ is the kinematic viscosity. We did not investigate dependence on binary mass ratio, eccentricity, disk temperature, or disk viscosity; therefore, these benchmarks may act as guides towards expanding converged solutions to the wider parameter space but might need to be updated in a future study that investigates dependence on system parameters. We find the most major discrepancies between codes resulted from the dimensionality of the setup (3D vs 2D disks). Beyond this, we find good agreement in the total torque on the binary between codes, although the partition of this torque between the gravitational torque, orbital accretion torque, and spin accretion torque depends sensitively on the sink prescriptions employed. In agreement with previous studies, we find a modest difference in torques and accretion variability between 2D and 3D disk models. We find cavity precession rates to be appreciably faster in 3D than in 2D.
△ Less
Submitted 20 June, 2024; v1 submitted 20 February, 2024;
originally announced February 2024.
-
The Disk Orientations of Perseus Protostellar Multiples at 8 au Resolution
Authors:
Nickalas K. Reynolds,
John J. Tobin,
Patrick D. Sheehan,
Sarah I. Sadavoy,
Leslie W. Looney,
Kaitlin M. Kratter,
Zhi-Yun Li,
Dominique M. Segura-Cox,
Nathan A. Kaib
Abstract:
We present a statistical characterization of circumstellar disk orientations toward 12 protostellar multiple systems in the Perseus molecular cloud using the Atacama Large Millimeter/submillimeter Array at Band 6 (1.3 mm) with a resolution of 25 mas (8 au). This exquisite resolution enabled us to resolve the compact inner disk structures surrounding the components of each multiple system and to de…
▽ More
We present a statistical characterization of circumstellar disk orientations toward 12 protostellar multiple systems in the Perseus molecular cloud using the Atacama Large Millimeter/submillimeter Array at Band 6 (1.3 mm) with a resolution of 25 mas (8 au). This exquisite resolution enabled us to resolve the compact inner disk structures surrounding the components of each multiple system and to determine the projected 3-D orientation of the disks (position angle and inclination) to high precision. We performed a statistical analysis on the relative alignment of disk pairs to determine whether the disks are preferentially aligned or randomly distributed. We considered three subsamples of the observations selected by the companion separations, a <100 au, a >500 au, and a < 10,000 au. We found for the compact (< 100 au) subsample, the distribution of orientation angles is best described by an underlying distribution of preferentially aligned sources (within 30deg) but does not rule out distributions with 40% misaligned sources. The wide companion (>500 au) subsample appears to be consistent with a distribution of 40%-80% preferentially aligned sources. Similarly, the full sample of systems with companions (a< 10, 000 au) is most consistent with a fractional ratio of at most 80% preferentially aligned source and rules out purely randomly aligned distributions. Thus our results imply the compact sources (<100 au) and the wide companions (>500 au) are statistically different.
△ Less
Submitted 20 December, 2023;
originally announced December 2023.
-
A Fast second-order solver for stiff multifluid dust and gas hydrodynamics
Authors:
Leonardo Krapp,
Juan Garrido-Deutelmoser,
Pablo Benítez-Llambay,
Kaitlin M. Kratter
Abstract:
We present MDIRK: a Multifluid second-order Diagonally-Implicit Runge-Kutta method to study momentum transfer between gas and an arbitrary number ($N$) of dust species. The method integrates the equations of hydrodynamics with an Implicit Explicit (IMEX) scheme and solves the stiff source term in the momentum equation with a diagonally-implicit asymptotically stable Runge-Kutta method (DIRK). In p…
▽ More
We present MDIRK: a Multifluid second-order Diagonally-Implicit Runge-Kutta method to study momentum transfer between gas and an arbitrary number ($N$) of dust species. The method integrates the equations of hydrodynamics with an Implicit Explicit (IMEX) scheme and solves the stiff source term in the momentum equation with a diagonally-implicit asymptotically stable Runge-Kutta method (DIRK). In particular, DIRK admits a simple analytical solution that can be evaluated with $\mathcal{O}(N)$ operations, instead of standard matrix inversion, which is $\mathcal{O}(N)^3$. Therefore the analytical solution significantly reduces the computational cost of the multifluid method, making it suitable for studying the dynamics of systems with particle-size distributions. We demonstrate that the method conserves momentum to machine precision and converges to the correct equilibrium solution with constant external acceleration. To validate our numerical method we present a series of simple hydrodynamic tests, including damping of sound waves, dusty shocks, a multi-fluid dusty Jeans instability, and a steady-state gas-dust drift calculation. The simplicity of MDIRK lays the groundwork to build fast high-order asymptotically stable multifluid methods.
△ Less
Submitted 26 February, 2024; v1 submitted 29 September, 2023;
originally announced October 2023.
-
The multiplicity of massive stars in the Scorpius OB1 association through high-contrast imaging
Authors:
Tinne Pauwels,
Maddalena Reggiani,
Hugues Sana,
Alan Rainot,
Kaitlin Kratter
Abstract:
One of the most remarkable properties of massive stars is that almost all of them are found in binaries or higher-order multiple systems. Observations that cover the full companion mass ratio and separation regime are essential to constrain massive star and binary formation theories. We used VLT/SPHERE to characterise the multiplicity properties of 20 OB stars in the active star-forming region Sco…
▽ More
One of the most remarkable properties of massive stars is that almost all of them are found in binaries or higher-order multiple systems. Observations that cover the full companion mass ratio and separation regime are essential to constrain massive star and binary formation theories. We used VLT/SPHERE to characterise the multiplicity properties of 20 OB stars in the active star-forming region Sco OB1. We simultaneously observed with the IFS and IRDIS instruments, obtaining high-contrast imaging observations that cover a field of view of 1".73 x 1".73 in YJH bands and 11" x 12".5 in $K_1$ and $K_2$ bands, respectively, corresponding to a separation range between $\sim$200 and 9000 AU. The observations reach contrast magnitudes down to $ΔK_1 \sim 13$, allowing us to detect companions at the stellar-substellar boundary. In total, we detect 789 sources, most of which are likely background or foreground objects. We obtain SPHERE companion fractions of $2.3 \pm 0.4$ and $4.2 \pm 0.8$ for O- and B-type stars, respectively. Including all previously detected companions, we find a total multiplicity fraction of $0.89\pm0.07$ for our sample in the range of $\sim$0-12000 AU. In conclusion, SPHERE explores an as of yet uncharted territory of companions around massive stars, which is crucial to ultimately improve our understanding of massive star and binary formation.
△ Less
Submitted 25 July, 2023;
originally announced July 2023.
-
Direct images and spectroscopy of a giant protoplanet driving spiral arms in MWC 758
Authors:
Kevin Wagner,
Jordan Stone,
Andrew Skemer,
Steve Ertel,
Ruobing Dong,
Dániel Apai,
Eckhart Spalding,
Jarron Leisenring,
Michael Sitko,
Kaitlin Kratter,
Travis Barman,
Mark Marley,
Brittany Miles,
Anthony Boccaletti,
Korash Assani,
Ammar Bayyari,
Taichi Uyama,
Charles E. Woodward,
Phil Hinz,
Zackery Briesemeister,
Kellen Lawson,
François Ménard,
Eric Pantin,
Ray W. Russell,
Michael Skrutskie
, et al. (1 additional authors not shown)
Abstract:
Understanding the driving forces behind spiral arms in protoplanetary disks remains a challenge due to the faintness of young giant planets. MWC 758 hosts such a protoplanetary disk with a two-armed spiral pattern that is suggested to be driven by an external giant planet. We present new thermal infrared observations that are uniquely sensitive to redder (i.e., colder or more attenuated) planets t…
▽ More
Understanding the driving forces behind spiral arms in protoplanetary disks remains a challenge due to the faintness of young giant planets. MWC 758 hosts such a protoplanetary disk with a two-armed spiral pattern that is suggested to be driven by an external giant planet. We present new thermal infrared observations that are uniquely sensitive to redder (i.e., colder or more attenuated) planets than past observations at shorter wavelengths. We detect a giant protoplanet, MWC 758c, at a projected separation of ~100 au from the star. The spectrum of MWC 758c is distinct from the rest of the disk and consistent with emission from a planetary atmosphere with Teff = 500 +/- 100 K for a low level of extinction (AV<30), or a hotter object with a higher level of extinction. Both scenarios are commensurate with the predicted properties of the companion responsible for driving the spiral arms. MWC 758c provides evidence that spiral arms in protoplanetary disks can be caused by cold giant planets or by those whose optical emission is highly attenuated. MWC 758c stands out both as one of the youngest giant planets known, and also as one of the coldest and/or most attenuated. Furthermore, MWC 758c is among the first planets to be observed within a system hosting a protoplanetary disk.
△ Less
Submitted 8 July, 2023;
originally announced July 2023.
-
RAM: Rapid Advection Algorithm on Arbitrary Meshes
Authors:
Pablo Benítez-Llambay,
Leonardo Krapp,
Ximena S. Ramos,
Kaitlin M. Kratter
Abstract:
The study of many astrophysical flows requires computational algorithms that can capture high Mach number flows, while resolving a large dynamic range in spatial and density scales. In this paper we present a novel method, RAM: Rapid Advection Algorithm on Arbitrary Meshes. RAM is a time-explicit method to solve the advection equation in problems with large bulk velocity on arbitrary computational…
▽ More
The study of many astrophysical flows requires computational algorithms that can capture high Mach number flows, while resolving a large dynamic range in spatial and density scales. In this paper we present a novel method, RAM: Rapid Advection Algorithm on Arbitrary Meshes. RAM is a time-explicit method to solve the advection equation in problems with large bulk velocity on arbitrary computational grids. In comparison with standard up-wind algorithms, RAM enables advection with larger time steps and lower truncation errors. Our method is based on the operator splitting technique and conservative interpolation. Depending on the bulk velocity and resolution, RAM can decrease the numerical cost of hydrodynamics by more than one order of magnitude. To quantify the truncation errors and speed-up with RAM, we perform one and two-dimensional hydrodynamics tests. We find that the order of our method is given by the order of the conservative interpolation and that the effective speed up is in agreement with the relative increment in time step. RAM will be especially useful for numerical studies of disk-satellite interaction, characterized by high bulk orbital velocities, and non-trivial geometries. Our method dramatically lowers the computational cost of simulations that simultaneously resolve the global disk and well inside the Hill radius of the secondary companion.
△ Less
Submitted 9 May, 2023;
originally announced May 2023.
-
Binary Formation in a 100 $μ$m-dark Massive Core
Authors:
Shuo Kong,
Héctor G. Arce,
John J. Tobin,
Yichen Zhang,
María José Maureira,
Kaitlin M. Kratter,
Thushara G. S. Pillai
Abstract:
We report high-resolution ALMA observations toward a massive protostellar core C1-Sa ($\sim$30 M$_\odot$) in the Dragon Infrared Dark Cloud. At the resolution of 140 AU, the core fragments into two kernels (C1-Sa1 and C1-Sa2) with a projected separation of $\sim$1400 AU along the elongation of C1-Sa, consistent with a Jeans length scale of $\sim$1100 AU. Radiative transfer modeling using RADEX ind…
▽ More
We report high-resolution ALMA observations toward a massive protostellar core C1-Sa ($\sim$30 M$_\odot$) in the Dragon Infrared Dark Cloud. At the resolution of 140 AU, the core fragments into two kernels (C1-Sa1 and C1-Sa2) with a projected separation of $\sim$1400 AU along the elongation of C1-Sa, consistent with a Jeans length scale of $\sim$1100 AU. Radiative transfer modeling using RADEX indicates that the protostellar kernel C1-Sa1 has a temperature of $\sim$75 K and a mass of 0.55 M$_\odot$. C1-Sa1 also likely drives two bipolar outflows, one being parallel to the plane-of-the-sky. C1-Sa2 is not detected in line emission and does not show any outflow activity but exhibits ortho-H$_2$D$^+$ and N$_2$D$^+$ emission in its vicinity, thus it is likely still starless. Assuming a 20 K temperature, C1-Sa2 has a mass of 1.6 M$_\odot$. At a higher resolution of 96 AU, C1-Sa1 begins to show an irregular shape at the periphery, but no clear sign of multiple objects or disks. We suspect that C1-Sa1 hosts a tight binary with inclined disks and outflows. Currently, one member of the binary is actively accreting while the accretion in the other is significantly reduced. C1-Sa2 shows hints of fragmentation into two sub-kernels with similar masses, which requires further confirmation with higher sensitivity.
△ Less
Submitted 3 May, 2023;
originally announced May 2023.
-
The Small Separation A-Star Companion Population: First Results with CHARA/MIRC-X
Authors:
Matthew De Furio,
Tyler Gardner,
John Monnier,
Michael R. Meyer,
Kaitlin Kratter,
Gail Schaefer,
Narsireddy Anugu,
Claire L. Davies,
Stefan Kraus,
Cyprien Lanthermann,
Jean-Baptiste Le Bouquin,
Jacob Ennis
Abstract:
We present preliminary results from our long-baseline interferometry (LBI) survey to constrain the multiplicity properties of intermediate-mass A-type stars within 80pc. Previous multiplicity studies of nearby stars exhibit orbital separation distributions well-fitted with a log-normal with peaks > 15au, increasing with primary mass. The A-star multiplicity survey of De Rosa et al. (2014), sensiti…
▽ More
We present preliminary results from our long-baseline interferometry (LBI) survey to constrain the multiplicity properties of intermediate-mass A-type stars within 80pc. Previous multiplicity studies of nearby stars exhibit orbital separation distributions well-fitted with a log-normal with peaks > 15au, increasing with primary mass. The A-star multiplicity survey of De Rosa et al. (2014), sensitive beyond 30au but incomplete below 100 au, found a log-normal peak around 390au. Radial velocity surveys of slowly-rotating, chemically peculiar Am stars identified a significant number of very close companions with periods $\leq$ 5 days, ~ 0.1au, a result similar to surveys of O- and B-type primaries. With the improved performance of LBI techniques, we can probe these close separations for normal A-type stars where other surveys are incomplete. Our initial sample consists of 27 A-type primaries with estimated masses between 1.44-2.49M$_{\odot}$ and ages 10-790Myr, which we observed with the MIRC-X instrument at the CHARA Array. We use the open source software CANDID to detect five companions, three of which are new, and derive a companion frequency of 0.19$^{+0.11}_{-0.06}$ over mass ratios 0.25-1.0 and projected separations 0.288-5.481 au. We find a probability of 10$^{-6}$ that our results are consistent with extrapolations based on previous models of the A-star companion population, over mass ratios and separations sampled. Our results show the need to explore these very close separations to inform our understanding of stellar formation and evolution processes.
△ Less
Submitted 2 November, 2022;
originally announced November 2022.
-
The McDonald Accelerating Stars Survey (MASS): Architecture of the Ancient Five-Planet Host System Kepler-444
Authors:
Zhoujian Zhang,
Brendan P. Bowler,
Trent J. Dupuy,
Timothy D. Brandt,
G. Mirek Brandt,
William D. Cochran,
Michael Endl,
Phillip J. MacQueen,
Kaitlin M. Kratter,
Howard T. Isaacson,
Kyle Franson,
Adam L. Kraus,
Caroline V. Morley,
Yifan Zhou
Abstract:
We present the latest and most precise characterization of the architecture for the ancient ($\approx 11$ Gyr) Kepler-444 system, which is composed of a K0 primary star (Kepler-444 A) hosting five transiting planets, and a tight M-type spectroscopic binary (Kepler-444 BC) with an A-BC projected separation of 66 au. We have measured the system's relative astrometry using the adaptive optics imaging…
▽ More
We present the latest and most precise characterization of the architecture for the ancient ($\approx 11$ Gyr) Kepler-444 system, which is composed of a K0 primary star (Kepler-444 A) hosting five transiting planets, and a tight M-type spectroscopic binary (Kepler-444 BC) with an A-BC projected separation of 66 au. We have measured the system's relative astrometry using the adaptive optics imaging from Keck/NIRC2 and Kepler-444 A's radial velocities from the Hobby Eberly Telescope, and re-analyzed relative radial velocities between BC and A from Keck/HIRES. We also include the Hipparcos-Gaia astrometric acceleration and all published astrometry and radial velocities into an updated orbit analysis of BC's barycenter. These data greatly extend the time baseline of the monitoring and lead to significant updates to BC's barycentric orbit compared to previous work, including a larger semi-major axis ($a = 52.2^{+3.3}_{-2.7}$ au), a smaller eccentricity ($e = 0.55 \pm 0.05$), and a more precise inclination ($i =85.4^{+0.3}_{-0.4}$ degrees). We have also derived the first dynamical masses of B and C components. Our results suggest Kepler-444~A's protoplanetary disk was likely truncated by BC to a radius of $\approx 8$ au, which resolves the previously noticed tension between Kepler-444 A's disk mass and planet masses. Kepler-444 BC's barycentric orbit is likely aligned with those of A's five planets, which might be primordial or a consequence of dynamical evolution. The Kepler-444 system demonstrates that compact multi-planet systems residing in hierarchical stellar triples can form at early epochs of the Universe and survive their secular evolution throughout cosmic time.
△ Less
Submitted 13 October, 2022;
originally announced October 2022.
-
ARMADA II: Further Detections of Inner Companions to Intermediate Mass Binaries with Micro-Arcsecond Astrometry at CHARA and VLTI
Authors:
Tyler Gardner,
John D. Monnier,
Francis C. Fekel,
Jean-Baptiste Le Bouquin,
Adam Scovera,
Gail Schaefer,
Stefan Kraus,
Fred C. Adams,
Narsireddy Anugu,
Jean-Philippe Berger,
Theo Ten Brummelaar,
Claire L. Davies,
Jacob Ennis,
Douglas R. Gies,
Keith J. C. Johnson,
Pierre Kervella,
Kaitlin M. Kratter,
Aaron Labdon,
Cyprien Lanthermann,
Johannes Sahlmann,
Benjamin R. Setterholm
Abstract:
We started a survey with CHARA/MIRC-X and VLTI/GRAVITY to search for low mass companions orbiting individual components of intermediate mass binary systems. With the incredible precision of these instruments, we can detect astrometric "wobbles" from companions down to a few tens of micro-arcseconds. This allows us to detect any previously unseen triple systems in our list of binaries. We present t…
▽ More
We started a survey with CHARA/MIRC-X and VLTI/GRAVITY to search for low mass companions orbiting individual components of intermediate mass binary systems. With the incredible precision of these instruments, we can detect astrometric "wobbles" from companions down to a few tens of micro-arcseconds. This allows us to detect any previously unseen triple systems in our list of binaries. We present the orbits of 12 companions around early F to B-type binaries, 9 of which are new detections and 3 of which are first astrometric detections of known RV companions. The masses of these newly detected components range from 0.45-1.3 solar masses. Our orbits constrain these systems to a high astrometric precision, with median residuals to the orbital fit of 20-50 micro-arcseconds in most cases. For 7 of these systems we include newly obtained radial velocity data, which help us to identify the system configuration and to solve for masses of individual components in some cases. Although additional RV measurements are needed to break degeneracy in the mutual inclination, we find that the majority of these inner triples are not well-aligned with the wide binary orbit. This hints that higher mass triples are more misaligned compared to solar and lower mass triples, though a thorough study of survey biases is needed. We show that the ARMADA survey is extremely successful at uncovering previously unseen companions in binaries. This method will be used in upcoming papers to constrain companion demographics in intermediate mass binary systems down to the planetary mass regime.
△ Less
Submitted 1 September, 2022;
originally announced September 2022.
-
Statistics of Galactic-Scale Quasar Pairs at Cosmic Noon
Authors:
Yue Shen,
Hsiang-Chih Hwang,
Masamune Oguri,
Nianyi Chen,
Tiziana Di Matteo,
Yueying Ni,
Simeon Bird,
Nadia Zakamska,
Xin Liu,
Yu-Ching Chen,
Kaitlin M. Kratter
Abstract:
The statistics of galactic-scale quasar pairs can elucidate our understanding of the dynamical evolution of supermassive black hole (SMBH) pairs, the duty cycles of quasar activity in mergers, or even the nature of dark matter, but have been challenging to measure at cosmic noon, the prime epoch of massive galaxy and SMBH formation. Here we measure a double quasar fraction of…
▽ More
The statistics of galactic-scale quasar pairs can elucidate our understanding of the dynamical evolution of supermassive black hole (SMBH) pairs, the duty cycles of quasar activity in mergers, or even the nature of dark matter, but have been challenging to measure at cosmic noon, the prime epoch of massive galaxy and SMBH formation. Here we measure a double quasar fraction of $\sim 6.2\pm0.5\times 10^{-4}$ integrated over $\sim 0.3-3$ arcsec separations (projected physical separations of $\sim 3-30\,{\rm kpc}$ at $z\sim 2$) in luminous ($L_{\rm bol}>10^{45.8}\,{\rm erg\,s^{-1}}$) unobscured quasars at $1.5<z<3.5$, using Gaia EDR3-resolved pairs around SDSS DR16 quasars. The measurement was based on a sample of 60 Gaia-resolved double quasars (out of 487 Gaia pairs dominated by quasar+star superpositions) at these separations, corrected for pair completeness in Gaia, which we quantify as functions of pair separation, magnitude of the primary, and magnitude contrast. The double quasar fraction increases towards smaller separations by a factor of $\sim 5$ over these scales. The division between physical quasar pairs and lensed quasars in our sample is currently unknown, requiring dedicated follow-up observations (in particular, deep, sub-arcsec-resolution IR imaging for the closest pairs). Intriguingly, at this point the observed pair statistics are in rough agreement with theoretical predictions both for the lensed quasar population in mock catalogs and for dual quasars in cosmological hydrodynamic simulations. Upcoming wide-field imaging/spectroscopic space missions such as Euclid, CSST and Roman, combined with targeted follow-up observations, will conclusively measure the abundances and host galaxy properties of galactic-scale quasar pairs, offset AGNs, and sub-arcsec lensed quasars across cosmic time.
△ Less
Submitted 9 August, 2022;
originally announced August 2022.
-
Multiplicity Statistics of Stars in the Sagittarius Dwarf Spheroidal Galaxy: Comparison to the Milky Way
Authors:
Victoria Bonidie,
Travis Court,
Christine Mazzola Daher,
Catherine E. Fielder,
Carles Badenes,
Jeffrey Newman,
Maxwell Moe,
Kaitlin M. Kratter,
Matthew G. Walker,
Steven R. Majewski,
Christian R. Hayes,
Sten Hasselquist,
Keivan Stassun,
Marina Kounkel,
Don Dixon,
Guy S. Stringfellow,
Joleen Carlberg,
Borja Anguiano,
Nathan De Lee,
Nicholas Troup
Abstract:
We use time-resolved spectra from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) to examine the distribution of radial velocity (RV) variations in 249 stars identified as members of the Sagittarius (Sgr) dwarf spheroidal (dSph) galaxy by Hayes et al (2020). We select Milky Way (MW) stars that have stellar parameters ($log(g)$, $T_{eff}$, and $[Fe/H]$) similar to those of the S…
▽ More
We use time-resolved spectra from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) to examine the distribution of radial velocity (RV) variations in 249 stars identified as members of the Sagittarius (Sgr) dwarf spheroidal (dSph) galaxy by Hayes et al (2020). We select Milky Way (MW) stars that have stellar parameters ($log(g)$, $T_{eff}$, and $[Fe/H]$) similar to those of the Sagittarius members by means of a k-d tree of dimension 3. We find that the shape of the distribution of RV shifts in Sgr dSph stars is similar to that measured in their MW analogs, but the total fraction of RV variable stars in the Sgr dSph is larger by a factor of $\sim 2$. After ruling out other explanations for this difference, we conclude that the fraction of close binaries in the Sgr dSph is intrinsically higher than in the MW. We discuss the implications of this result for the physical processes leading to the formation of close binaries in dwarf spheroidal and spiral galaxies.
△ Less
Submitted 20 April, 2022;
originally announced April 2022.
-
Substructures in protoplanetary disks imprinted by compact planetary systems
Authors:
Juan Garrido-Deutelmoser,
Cristobal Petrovich,
Leonardo Krapp,
Kaitlin M. Kratter,
Ruobing Dong
Abstract:
The substructures observed in protoplanetary disks may be the signposts of embedded planets carving gaps or creating vortices. The inferred masses of these planets often fall in the Jovian regime despite their low abundance compared to lower-mass planets, partly because previous works often assume that a single substructure (a gap or vortex) is caused by a single planet. In this work, we study the…
▽ More
The substructures observed in protoplanetary disks may be the signposts of embedded planets carving gaps or creating vortices. The inferred masses of these planets often fall in the Jovian regime despite their low abundance compared to lower-mass planets, partly because previous works often assume that a single substructure (a gap or vortex) is caused by a single planet. In this work, we study the possible imprints of compact systems composed of Neptune-like planets ($\sim10-30\;M_\oplus$) and show that long-standing vortices are a prevalent outcome when their inter-planetary separation ($Δa$) falls below $\sim8$ times $H_{\rm p}$ -- the average disk's scale height at the planets locations. In simulations where a single planet is unable to produce long-lived vortices, two-planet systems can preserve them for at least $5,000$ orbits in two regimes: i) fully-shared density gaps with elongated vortices around the stable Lagrange points $L_4$ and $L_5$ for the most compact planet pairs ($Δa \lesssim 4.6\; H_{\rm p}$); ii) partially-shared gaps for more widely spaced planets ($Δa \sim 4.6 - 8\;H_{\rm p}$) forming vortices in a density ring between the planets through the Rossby wave instability. The latter case can produce vortices with a wide range of aspect ratios down to $\sim3$ and can occur for planets captured into the 3:2 (2:1) mean-motion resonances for disk's aspects ratios of $h\gtrsim 0.033$ ($h\gtrsim 0.057$). We suggest that their long lifetimes are sustained by the interaction of spiral density waves launched by the neighboring planets. Overall, our results show that distinguishing imprint of compact systems with Neptune-mass planets are long-lived vortices inside the density gaps, which in turn are shallower than single-planet gaps for a fixed gap width.
△ Less
Submitted 19 April, 2022;
originally announced April 2022.
-
The Origin and Evolution of Multiple Star Systems
Authors:
Stella S. R. Offner,
Maxwell Moe,
Kaitlin M. Kratter,
Sarah I. Sadavoy,
Eric L. N. Jensen,
John J. Tobin
Abstract:
Observational advances over the last decade have enabled high-resolution, interferometric studies of forming multiple systems, statistical surveys of multiplicity in star-forming regions, and new insights into disk evolution and planetary architectures in these systems. In this review, we compile the results of observational and theoretical studies of stellar multiplicity. We summarize the populat…
▽ More
Observational advances over the last decade have enabled high-resolution, interferometric studies of forming multiple systems, statistical surveys of multiplicity in star-forming regions, and new insights into disk evolution and planetary architectures in these systems. In this review, we compile the results of observational and theoretical studies of stellar multiplicity. We summarize the population statistics spanning system evolution from the protostellar phase through the main-sequence phase and evaluate the influence of the local environment. In short, most stars are born in multiple stellar systems, most main sequence stars are members of multiple systems, but most star systems are single. We describe current models for the origin of stellar multiplicity and review the landscape of numerical simulations and assess their consistency with observations. We review the properties of disks and discuss the impact of multiplicity on planet formation and system architectures. Finally, we summarize open questions and discuss the technical requirements for future observational and theoretical progress.
△ Less
Submitted 7 December, 2022; v1 submitted 18 March, 2022;
originally announced March 2022.
-
Orbital Architectures of Planet-Hosting Binaries II. Low Mutual Inclinations Between Planetary and Stellar Orbits
Authors:
Trent J. Dupuy,
Adam L. Kraus,
Kaitlin M. Kratter,
Aaron C. Rizzuto,
Andrew W. Mann,
Daniel Huber,
Michael J. Ireland
Abstract:
Planet formation is often considered in the context of one circumstellar disk around one star. Yet stellar binary systems are ubiquitous, and thus a substantial fraction of all potential planets must form and evolve in more complex, dynamical environments. We present the results of a five-year astrometric monitoring campaign studying 45 binary star systems that host Kepler planet candidates. The p…
▽ More
Planet formation is often considered in the context of one circumstellar disk around one star. Yet stellar binary systems are ubiquitous, and thus a substantial fraction of all potential planets must form and evolve in more complex, dynamical environments. We present the results of a five-year astrometric monitoring campaign studying 45 binary star systems that host Kepler planet candidates. The planet-forming environments in these systems would have literally been shaped by the binary orbits that persist to the present day. Crucially, the mutual inclinations of star-planet orbits can only be addressed by a statistical sample. We describe in detail our sample selection and Keck/NIRC2 laser guide star adaptive optics observations collected from 2012 to 2017. We measure orbital arcs, with a typical accuracy of ~0.1 mas/yr, that test whether the binary orbits tend to be aligned with the edge-on transiting planet orbits. We rule out randomly-distributed binary orbits at 4.7$σ$, and we show that low mutual inclinations are required to explain the observed orbital arcs. If the stellar orbits have a field binary-like eccentricity distribution, then the best match to our observed orbital arcs is a distribution of mutual inclinations ranging from 0-30 degrees. We discuss the implications of such widespread planet-binary alignment in the theoretical context of planet formation and circumstellar disk evolution.
△ Less
Submitted 31 January, 2022;
originally announced February 2022.
-
Probing the low-mass end of the companion mass function for O-type stars
Authors:
M. Reggiani,
A. Rainot,
H. Sana,
L. A. Almeida,
S. Caballero-Nieves,
K. Kratter,
S. Lacour,
J. -B. LeBouquin,
H. Zinnecker
Abstract:
Almost all massive stars are part of a binary system. Given the wide range of separations at which these companions are found, several observational techniques have been adopted to characterize them, but contrasts greater than 4 in the H-band have never been reached between 0".1 and 1". We used VLT/SPHERE to observe simultaneously with the IRDIS and IFS sub-systems 18 O-type stars within 6 kpc and…
▽ More
Almost all massive stars are part of a binary system. Given the wide range of separations at which these companions are found, several observational techniques have been adopted to characterize them, but contrasts greater than 4 in the H-band have never been reached between 0".1 and 1". We used VLT/SPHERE to observe simultaneously with the IRDIS and IFS sub-systems 18 O-type stars within 6 kpc and ages between 1-5 Myrs to probe the existence of stellar companions in the angular separation range from 0".15 to 6" down to very low mass ratios. The IFS YJH- band observations have allowed us to probe the presence of sub-solar companions in a 1".7x1".7 field-of-view down to magnitude limits of deltaH=10 at 0".4. In the wider 12"x12" IRDIS field-of-view, we reached contrasts of deltaK=12 at 1", enabling us to look for even fainter companions. This paper presents five newly discovered intermediate (<1") separation companions, three of which are smaller than 0.2M_sun. If confirmed by future analyses of proper motions, these new detections represent the lowest-mass companions ever found around O-type stars. Assuming that all sources detected within 1" are physically bound, the observed fraction of companions for O-type stars between 0".15 and 0".9 is 0.39+/-0.15, whereas it increases to 1.6+/-0.3 in the separation range from 0".9 to 6". These findings clearly support the notion that massive stars form almost exclusively in multiple systems, and that larger AO-assisted coronagraphic surveys are crucial in placing constraints on the multiplicity properties of massive star companions in regions of the parameter space that have previously gone unexplored, and demonstrate that the companion mass function is populated down to the lowest stellar masses.
△ Less
Submitted 20 December, 2021;
originally announced December 2021.
-
The VLA/ALMA Nascent Disk and Multiplicity (VANDAM) Survey of Orion Protostars V. A Characterization of Protostellar Multiplicity
Authors:
John J. Tobin,
Stella S. R. Offner,
Kaitlin M. Kratter,
S. Thomas Megeath,
Patrick D. Sheehan,
Leslie W. Looney,
Ana Karla Diaz-Rodriguez,
Mayra Osorio,
Guillem Anglada,
Sarah I. Sadavoy,
Elise Furlan,
Dominique Segura-Cox,
Nicole Karnath,
Merel L. R. van 't Hoff,
Ewine F. van Dishoeck,
Zhi-Yun Li,
Rajeeb Sharma,
Amelia M. Stutz,
Lukasz Tychoniec
Abstract:
We characterize protostellar multiplicity in the Orion molecular clouds using ALMA 0.87~mm and VLA 9~mm continuum surveys toward 328 protostars. These observations are sensitive to projected spatial separations as small as $\sim$20~au, and we consider source separations up to 10$^4$~au as potential companions. The overall multiplicity fraction (MF) and companion fraction (CF) for the Orion protost…
▽ More
We characterize protostellar multiplicity in the Orion molecular clouds using ALMA 0.87~mm and VLA 9~mm continuum surveys toward 328 protostars. These observations are sensitive to projected spatial separations as small as $\sim$20~au, and we consider source separations up to 10$^4$~au as potential companions. The overall multiplicity fraction (MF) and companion fraction (CF) for the Orion protostars are 0.30$\pm$0.03 and 0.44$\pm$0.03, respectively, considering separations from 20 to 10$^4$~au. The MFs and CFs are corrected for potential contamination by unassociated young stars using a probabilistic scheme based on the surface density of young stars around each protostar. The companion separation distribution as a whole is double peaked and inconsistent with the separation distribution of solar-type field stars, while the separation distribution of Flat Spectrum protostars is consistent solar-type field stars. The multiplicity statistics and companion separation distributions of the Perseus star-forming region are consistent with those of Orion. Based on the observed peaks in the Class 0 separations at $\sim$100~au and $\sim$10$^3$~au, we argue that multiples with separations $<$500~au are likely produced by both disk fragmentation and turbulent fragmentation with migration, and those at $\ga$10$^3$~au result primarily from turbulent fragmentation. We also find that MFs/CFs may rise from Class 0 to Flat Spectrum protostars between 100 and 10$^3$~au in regions of high YSO density. This finding may be evidence for migration of companions from $>$10$^3$~au to $<$10$^3$~au, and that some companions between 10$^3$ and 10$^4$~au must be (or become) unbound.
△ Less
Submitted 8 November, 2021;
originally announced November 2021.
-
A Sublime Opportunity: The Dynamics of Transitioning Cometary Bodies and the Feasibility of $\textit{In Situ}$ Observations of The Evolution of Their Activity
Authors:
Darryl Z. Seligman,
Kaitlin M. Kratter,
W. Garrett Levine,
Robert Jedicke
Abstract:
The compositional and morphological evolution of minor bodies in the Solar System is primarily driven by the evolution of their heliocentric distances, as the level of incident solar radiation regulates cometary activity. We investigate the dynamical transfer of Centaurs into the inner Solar System, facilitated by mean motion resonances with Jupiter and Saturn. The recently discovered object, P/20…
▽ More
The compositional and morphological evolution of minor bodies in the Solar System is primarily driven by the evolution of their heliocentric distances, as the level of incident solar radiation regulates cometary activity. We investigate the dynamical transfer of Centaurs into the inner Solar System, facilitated by mean motion resonances with Jupiter and Saturn. The recently discovered object, P/2019 LD2, will transition from the Centaur region to the inner Solar System in 2063. In order to contextualize LD2, we perform N-body simulations of a population of Centaurs and JFCs. Objects between Jupiter and Saturn with Tisserand parameter $T_J\sim$3 are transferred onto orbits with perihelia $q<4$au within the next 1000 years with notably high efficiency. Our simulations show that there may be additional LD2-like objects transitioning into the inner Solar System in the near-term future, all of which have low $Δ$V with respect to Jupiter. We calculate the distribution of orbital elements resulting from a single Jovian encounter and show that objects with initial perihelia close to Jupiter are efficiently scattered to $q<4$au. Moreover, approximately $55\%$ of the transitioning objects in our simulated population experience at least 1 Jovian encounter prior to reaching $q<4$au. We demonstrate that a spacecraft stationed near Jupiter would be well-positioned to rendezvous, orbit match, and accompany LD2 into the inner Solar System, providing an opportunity to observe the onset of intense activity in a pristine comet $\textit{in situ}$. Finally, we discuss the prospect of identifying additional targets for similar measurements with forthcoming observational facilities.
△ Less
Submitted 6 October, 2021;
originally announced October 2021.
-
The 3D dust and opacity distribution of protoplanets in multi-fluid global simulations
Authors:
Leonardo Krapp,
Kaitlin M. Kratter,
Andrew N. Youdin
Abstract:
The abundance and distribution of solids inside the Hill sphere are central to our understanding of the giant planet dichotomy. Here, we present a three-dimensional characterization of the dust density, mass flux, and mean opacities in the envelope of sub-thermal and super-thermal mass planets. We simulate the dynamics of multiple dust species in a global protoplanetary disk model accounting for d…
▽ More
The abundance and distribution of solids inside the Hill sphere are central to our understanding of the giant planet dichotomy. Here, we present a three-dimensional characterization of the dust density, mass flux, and mean opacities in the envelope of sub-thermal and super-thermal mass planets. We simulate the dynamics of multiple dust species in a global protoplanetary disk model accounting for dust feedback. We find that the meridional flows do not effectively stir dust grains at scales of the Bondi sphere. Thus the dust-settling driven by the stellar gravitational potential sets the latitudinal dust density gradient within the planet envelope. Not only does the planet's potential enhance this gradient, but also the spiral wakes serve as another source of asymmetry. These asymmetries substantially alter the inferred mean Rosseland and Planck opacities. In cases with the moderate-to-strong dust settling, the opacity gradient can range from a few percent to more than two orders of magnitude between the mid-plane and the polar regions of the Bondi sphere. Finally, we show that this strong latitudinal opacity gradient can introduce a transition between optically thick and thin regimes at the scales of the planet envelope. We suggest that this transition is likely to occur when the equilibrium scale height of hundred-micron-sized particles is smaller than the Hill radius of the forming planet. This work calls into question the adoption of a constant opacity derived from well-mixed distributions and demonstrates the need for global radiation hydrodynamics models of giant planet formation which account for dust dynamics.
△ Less
Submitted 9 March, 2022; v1 submitted 5 October, 2021;
originally announced October 2021.
-
Stellar multiplicity and stellar rotation: Insights from APOGEE
Authors:
Christine Mazzola Daher,
Carles Badenes,
Jamie Tayar,
Marc Pinsonneault,
Sergey E. Koposov,
Kaitlin Kratter,
Maxwell Moe,
Borja Anguiano,
Diego Godoy-Rivera,
Steven Majewski,
Joleen K. Carlberg,
Matthew G. Walker,
Rachel Buttry,
Don Dixon,
Javier Serna,
Keivan G. Stassun,
Nathan De Lee,
Jesús Hernández,
Christian Nitschelm,
Guy S. Stringfellow,
Nicholas W. Troup
Abstract:
We measure rotational broadening in spectra taken by the Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey to characterise the relationship between stellar multiplicity and rotation. We create a sample of 2786 giants and 24 496 dwarfs with stellar parameters and multiple radial velocities from the APOGEE pipeline, projected rotation speeds \vsini\ determined from our own pipel…
▽ More
We measure rotational broadening in spectra taken by the Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey to characterise the relationship between stellar multiplicity and rotation. We create a sample of 2786 giants and 24 496 dwarfs with stellar parameters and multiple radial velocities from the APOGEE pipeline, projected rotation speeds \vsini\ determined from our own pipeline, and distances, masses, and ages measured by Sanders \& Das. We use the statistical distribution of the maximum shift in the radial velocities, \drvm, as a proxy for the close binary fraction to explore the interplay between stellar evolution, rotation, and multiplicity. Assuming that the minimum orbital period allowed is the critical period for Roche Lobe overflow and rotational synchronization, we calculate theoretical upper limits on expected \vsini\ and \drvm\ values. These expectations agree with the positive correlation between the maximum \drvm\ and \vsini\ values observed in our sample as a function of \logg. We find that the fast rotators in our sample have a high occurrence of short-period ($\log(P/\text{d})\lesssim 4$) companions. We also find that old, rapidly-rotating main sequence stars have larger completeness-corrected close binary fractions than their younger peers. Furthermore, rapidly-rotating stars with large \drvm\ consistently show differences of 1-10 Gyr between the predicted gyrochronological and measured isochronal ages. These results point towards a link between rapid rotation and close binarity through tidal interactions. We conclude that stellar rotation is strongly correlated with stellar multiplicity in the field, and caution should be taken in the application of gyrochronology relations to cool stars.
△ Less
Submitted 1 March, 2022; v1 submitted 3 October, 2021;
originally announced October 2021.
-
Double-lined spectroscopic binaries in the APOGEE DR16 and DR17 data
Authors:
Marina Kounkel,
Kevin R. Covey,
Keivan G. Stassun,
Adrian M. Price-Whelan,
Jon Holtzman,
Drew Chojnowski,
Penélope Longa-Peña,
Carlos G. Román-Zúñiga,
Jesus Hernandez,
Javier Serna,
Carles Badenes,
Nathan De Lee,
Steven Majewski,
Guy S. Stringfellow,
Kaitlin M. Kratter,
Maxwell Moe,
Peter M. Frinchaboy,
Rachael L. Beaton,
José G. Fernández-Trincado,
Suvrath Mahadevan,
Dante Minniti,
Timothy C. Beers,
Donald P. Schneider,
Rodolfo H. Barbá,
Joel R. Brownstein
, et al. (3 additional authors not shown)
Abstract:
APOGEE spectra offer $\lesssim$1 km s$^{-1}$ precision in the measurement of stellar radial velocities (RVs). This holds even when multiple stars are captured in the same spectrum, as happens most commonly with double-lined spectroscopic binaries (SB2s), although random line of sight alignments of unrelated stars can also occur. We develop a code that autonomously identifies SB2s and higher order…
▽ More
APOGEE spectra offer $\lesssim$1 km s$^{-1}$ precision in the measurement of stellar radial velocities (RVs). This holds even when multiple stars are captured in the same spectrum, as happens most commonly with double-lined spectroscopic binaries (SB2s), although random line of sight alignments of unrelated stars can also occur. We develop a code that autonomously identifies SB2s and higher order multiples in the APOGEE spectra, resulting in 7273 candidate SB2s, 813 SB3s, and 19 SB4s. We estimate the mass ratios of binaries, and for a subset of these systems with sufficient number of measurements we perform a complete orbital fit, confirming that most systems with period $<$10 days have circularized. Overall, we find a SB2 fraction ($F_{SB2}$) $\sim$3\% among main sequence dwarfs, and that there is not a significant trend in $F_{SB2}$ with temperature of a star. We are also able to recover a higher $F_{SB2}$ in sources with lower metallicity, however there are some observational biases. We also examine light curves from TESS to determine which of these spectroscopic binaries are also eclipsing. Such systems, particularly those that are also pre- and post-main sequence, are good candidates for a follow-up analysis to determine their masses and temperatures.
△ Less
Submitted 22 July, 2021;
originally announced July 2021.
-
On the Origin of the Pluto System
Authors:
Robin M. Canup,
Kaitlin M. Kratter,
Marc Neveu
Abstract:
The goal of this chapter is to review hypotheses for the origin of the Pluto system in light of observational constraints that have been considerably refined over the 85-year interval between the discovery of Pluto and its exploration by spacecraft. We focus on the giant impact hypothesis currently understood as the likeliest origin for the Pluto-Charon binary, and devote particular attention to n…
▽ More
The goal of this chapter is to review hypotheses for the origin of the Pluto system in light of observational constraints that have been considerably refined over the 85-year interval between the discovery of Pluto and its exploration by spacecraft. We focus on the giant impact hypothesis currently understood as the likeliest origin for the Pluto-Charon binary, and devote particular attention to new models of planet formation and migration in the outer solar system. We discuss the origins conundrum posed by the system's four small moons. We also elaborate on the implications of these scenarios for the dynamical environment of the early transneptunian disk, the likelihood of finding a Pluto collisional family, and the origin of other binary systems in the Kuiper belt. Finally, we highlight outstanding open issues regarding the origins of the Pluto system and suggest areas of future progress.
△ Less
Submitted 16 July, 2021;
originally announced July 2021.
-
The McDonald Accelerating Stars Survey (MASS): Discovery of a Long-Period Substellar Companion Orbiting the Old Solar Analog HD 47127
Authors:
Brendan P. Bowler,
Michael Endl,
William D. Cochran,
Phillip J. MacQueen,
Justin R. Crepp,
Greg W. Doppmann,
Shannon Dulz,
Timothy D. Brandt,
G. Mirek Brandt,
Yiting Li,
Trent J. Dupuy,
Kyle Franson,
Kaitlin M. Kratter,
Caroline V. Morley,
Yifan Zhou
Abstract:
Brown dwarfs with well-determined ages, luminosities, and masses provide rare but valuable tests of low-temperature atmospheric and evolutionary models. We present the discovery and dynamical mass measurement of a substellar companion to HD 47127, an old ($\approx$7-10 Gyr) G5 main sequence star with a mass similar to the Sun. Radial velocities of the host star with the Harlan J. Smith Telescope u…
▽ More
Brown dwarfs with well-determined ages, luminosities, and masses provide rare but valuable tests of low-temperature atmospheric and evolutionary models. We present the discovery and dynamical mass measurement of a substellar companion to HD 47127, an old ($\approx$7-10 Gyr) G5 main sequence star with a mass similar to the Sun. Radial velocities of the host star with the Harlan J. Smith Telescope uncovered a low-amplitude acceleration of 1.93 $\pm$ 0.08 m s$^{-1}$ yr$^{-1}$ based on 20 years of monitoring. We subsequently recovered a faint ($ΔH$=13.14 $\pm$ 0.15 mag) co-moving companion at 1.95$''$ (52 AU) with follow-up Keck/NIRC2 adaptive optics imaging. The radial acceleration of HD 47127 together with its tangential acceleration from Hipparcos and Gaia EDR3 astrometry provide a direct measurement of the three-dimensional acceleration vector of the host star, enabling a dynamical mass constraint for HD 47127 B (67.5-177 $M_\mathrm{Jup}$ at 95% confidence) despite the small fractional orbital coverage of the observations. The absolute $H$-band magnitude of HD 47127 B is fainter than the benchmark T dwarfs HD 19467 B and Gl 229 B but brighter than Gl 758 B and HD 4113 C, suggesting a late-T spectral type. Altogether the mass limits for HD 47127 B from its dynamical mass and the substellar boundary imply a range of 67-78 $M_\mathrm{Jup}$ assuming it is single, although a preference for high masses of $\approx$100 $M_\mathrm{Jup}$ from dynamical constraints hints at the possibility that HD 47127 B could itself be a binary pair of brown dwarfs or that another massive companion resides closer in. Regardless, HD 47127 B will be an excellent target for more refined orbital and atmospheric characterization in the future.
△ Less
Submitted 3 May, 2021;
originally announced May 2021.
-
Which planets trigger longer-lived vortices: low-mass or high-mass?
Authors:
Michael Hammer,
Min-Kai Lin,
Kaitlin M. Kratter,
Paola Pinilla
Abstract:
Recent ALMA observations have found many protoplanetary discs with rings that can be explained by gap-opening planets less massive than Jupiter. Meanwhile, recent studies have suggested that protoplanetary discs should have low levels of turbulence. Past computational work on low-viscosity discs has hinted that these two developments might not be self-consistent because even low-mass planets can b…
▽ More
Recent ALMA observations have found many protoplanetary discs with rings that can be explained by gap-opening planets less massive than Jupiter. Meanwhile, recent studies have suggested that protoplanetary discs should have low levels of turbulence. Past computational work on low-viscosity discs has hinted that these two developments might not be self-consistent because even low-mass planets can be accompanied by vortices instead of conventional double rings. We investigate this potential discrepancy by conducting hydrodynamic simulations of growing planetary cores in discs with various aspect ratios ($H/r=0.04$, 0.06, 0.08) and viscosities ($1.5 \times 10^{-5} \lesssim α\lesssim 3 \times 10^{-4}$), having these cores accrete their gas mass directly from the disc. With $α< 10^{-4}$, we find that sub-Saturn-mass planets in discs with $H/r \le 0.06$ are more likely to be accompanied by dust asymmetries compared to Jupiter-mass planets because they can trigger several generations of vortices in succession. We also find that vortices with $H/r = 0.08$ survive $>6000$ planet orbits regardless of the planet mass or disc mass because they are less affected by the planet's spiral waves. We connect our results to observations and find that the outward migration of vortices with $H/r \ge 0.08$ may be able to explain the cavity in Oph IRS 48 or the two clumps in MWC 758. Lastly, we show that the lack of observed asymmetries in the disc population in Taurus is unexpected given the long asymmetry lifetimes in our low viscosity simulations ($α\sim 2 \times 10^{-5}$), a discrepancy we suggest is due to these discs having higher viscosities.
△ Less
Submitted 9 May, 2021; v1 submitted 6 April, 2021;
originally announced April 2021.
-
The McDonald Accelerating Stars Survey (MASS): White Dwarf Companions Accelerating the Sun-like Stars 12 Psc and HD 159062
Authors:
Brendan P. Bowler,
William D. Cochran,
Michael Endl,
Kyle Franson,
Timothy D. Brandt,
Trent J. Dupuy,
Phillip J. MacQueen,
Kaitlin M. Kratter,
Dimitri Mawet,
Garreth Ruane
Abstract:
We present the discovery of a white dwarf companion to the G1 V star 12 Psc found as part of a Keck adaptive optics imaging survey of long-term accelerating stars from the McDonald Observatory Planet Search Program. Twenty years of precise radial-velocity monitoring of 12 Psc with the Tull Spectrograph at the Harlan J. Smith telescope reveals a moderate radial acceleration ($\approx$10 m s$^{-1}$…
▽ More
We present the discovery of a white dwarf companion to the G1 V star 12 Psc found as part of a Keck adaptive optics imaging survey of long-term accelerating stars from the McDonald Observatory Planet Search Program. Twenty years of precise radial-velocity monitoring of 12 Psc with the Tull Spectrograph at the Harlan J. Smith telescope reveals a moderate radial acceleration ($\approx$10 m s$^{-1}$ yr $^{-1}$), which together with relative astrometry from Keck/NIRC2 and the astrometric acceleration between $Hipparcos$ and $Gaia$ DR2 yields a dynamical mass of $M_B$ = 0.605$^{+0.021}_{-0.022}$ $M_{\odot}$ for 12 Psc B, a semi-major axis of 40$^{+2}_{-4}$ AU, and an eccentricity of 0.84$\pm$0.08. We also report an updated orbit fit of the white dwarf companion to the metal-poor (but barium-rich) G9 V dwarf HD 159062 based on new radial velocity observations from the High-Resolution Spectrograph at the Hobby-Eberly Telescope and astrometry from Keck/NIRC2. A joint fit of the available relative astrometry, radial velocities, and tangential astrometric acceleration yields a dynamical mass of $M_B$ = 0.609$^{+0.010}_{-0.011}$ $M_{\odot}$ for HD 159062 B, a semi-major axis of 60$^{+5}_{-7}$ AU, and preference for circular orbits ($e$$<$0.42 at 95% confidence). 12 Psc B and HD 159062 B join a small list of resolved "Sirius-like" benchmark white dwarfs with precise dynamical mass measurements which serve as valuable tests of white dwarf mass-radius cooling models and probes of AGB wind accretion onto their main-sequence companions.
△ Less
Submitted 8 December, 2020;
originally announced December 2020.
-
ARMADA I: Triple Companions Detected in B-Type Binaries alpha Del and nu Gem
Authors:
Tyler Gardner,
John D. Monnier,
Francis C. Fekel,
Gail Schaefer,
Keith J. C. Johnson,
Jean-Baptiste Le Bouquin,
Stefan Kraus,
Narsireddy Anugu,
Benjamin R. Setterholm,
Aaron Labdon,
Claire L. Davies,
Cyprien Lanthermann,
Jacob Ennis,
Michael Ireland,
Kaitlin M. Kratter,
Theo Ten Brummelaar,
Judit Sturmann,
Laszlo Sturmann,
Chris Farrington,
Douglas R. Gies,
Robert Klement,
Fred C. Adams
Abstract:
Ground-based optical long-baseline interferometry has the power to measure the orbits of close binary systems at ~10 micro-arcsecond precision. This precision makes it possible to detect "wobbles" in the binary motion due to the gravitational pull from additional short period companions. We started the ARrangement for Micro-Arcsecond Differential Astrometry (ARMADA) survey with the MIRC-X instrume…
▽ More
Ground-based optical long-baseline interferometry has the power to measure the orbits of close binary systems at ~10 micro-arcsecond precision. This precision makes it possible to detect "wobbles" in the binary motion due to the gravitational pull from additional short period companions. We started the ARrangement for Micro-Arcsecond Differential Astrometry (ARMADA) survey with the MIRC-X instrument at the CHARA array for the purpose of detecting giant planets and stellar companions orbiting individual stars in binary systems. We describe our observations for the survey, and introduce the wavelength calibration scheme that delivers precision at the tens of micro-arcseconds level for <0.2 arcsecond binaries. We test our instrument performance on a known triple system kappa Peg, and show that our survey is delivering a factor of 10 better precision than previous similar surveys. We present astrometric detections of tertiary components to two B-type binaries: a 30-day companion to alpha Del, and a 50-day companion to nu Gem. We also collected radial velocity data for alpha Del with the Tennessee State University Automated Spectroscopic Telescope at Fairborn Observatory. We are able to measure the orbits and masses of all three components in these systems. We find that the previously published RV orbit for the inner pair of nu Gem is not consistent with our visual orbit. The precision achieved for these orbits suggests that our ARMADA survey will be successful at discovering new compact triple systems to A/B-type binary systems, leading to better statistics of hierarchical system architectures and formation history.
△ Less
Submitted 1 December, 2020;
originally announced December 2020.
-
Kinematic Analysis of a Protostellar Multiple System: Measuring the Protostar Masses and Assessing Gravitational Instability in the Disks of L1448 IRS3B and L1448 IRS3A
Authors:
Nickalas K. Reynolds,
John J. Tobin,
Patrick D. Sheehan,
Sarah I. Sadavoy,
Kaitlin M. Kratter,
Zhi-Yun Li,
Claire J. Chandler,
Dominique M. Segura-Cox,
Leslie W. Looney,
Michael M. Dunham
Abstract:
We present new Atacama Large Millimeter/submillimeter Array (ALMA) observations towards a compact (230~au separation) triple protostar system, L1448 IRS3B, at 879~\micron with \contbeam~resolution. Spiral arm structure within the circum-multiple disk is well resolved in dust continuum toward IRS3B, and we detect the known wide (2300~au) companion, IRS3A, also resolving possible spiral substructure…
▽ More
We present new Atacama Large Millimeter/submillimeter Array (ALMA) observations towards a compact (230~au separation) triple protostar system, L1448 IRS3B, at 879~\micron with \contbeam~resolution. Spiral arm structure within the circum-multiple disk is well resolved in dust continuum toward IRS3B, and we detect the known wide (2300~au) companion, IRS3A, also resolving possible spiral substructure. Using dense gas tracers, C17O, H13CO$+$, and H13CN, we resolve the Keplerian rotation for both the circum-triple disk in IRS3B and the disk around IRS3A. Furthermore, we use the molecular line kinematic data and radiative transfer modeling of the molecular line emission to confirm that the disks are in Keplerian rotation with fitted masses of $1.19^{+0.13}_{-0.07}$ for IRS3B-ab, $1.51^{+0.06}_{-0.07}$~Msun for IRS3A, and place an upper limit on the central protostar mass for the tertiary IRS3B-c of 0.2~Msun. We measure the mass of the fragmenting disk of IRS3B to be 0.29~Msun from the dust continuum emission of the circum-multiple disk and estimate the mass of the clump surrounding IRS3B-c to be 0.07~Msun. We also find that the disk around IRS3A has a mass of 0.04~Msun. By analyzing the Toomre~Q parameter, we find the IRS3A circumstellar disk is gravitationally stable (Q$>$5), while the IRS3B disk is consistent with a gravitationally unstable disk (Q$<$1) between the radii 200-500~au. This coincides with the location of the spiral arms and the tertiary companion IRS3B-c, supporting the hypothesis that IRS3B-c was formed in situ via fragmentation of a gravitationally unstable disk.
△ Less
Submitted 16 November, 2020;
originally announced November 2020.
-
The VLA/ALMA Nascent Disk and Multiplicity (VANDAM) Survey of Orion Protostars IV. Unveiling the Embedded Intermediate-Mass Protostar and Disk within OMC2-FIR3/HOPS-370
Authors:
John J. Tobin,
Patrick Sheehan,
Nickalas Reynolds,
S. Thomas Megeath,
Mayra Osorio,
Guillem Anglada,
Ana Karla Diaz-Rodriguez,
Elise Furlan,
Kaitlin Kratter,
Stella Offner,
Leslie Looney,
Mihkel Kama,
Zhi-Yun Li,
Merel van 't Hoff,
Sarah Sadavoy,
Nicole Karnath
Abstract:
We present ALMA (0.87~mm and 1.3~mm) and VLA (9~mm) observations toward the candidate intermediate-mass protostar OMC2-FIR3 (HOPS-370; L$_{bol}$~314~L$_{\odot}$) at $\sim$0.1" (40~au) resolution for the continuum emission and ~0.25" (100 au) resolution of nine molecular lines. The dust continuum observed with ALMA at 0.87~mm and 1.3~mm resolve a near edge-on disk toward HOPS-370 with an apparent r…
▽ More
We present ALMA (0.87~mm and 1.3~mm) and VLA (9~mm) observations toward the candidate intermediate-mass protostar OMC2-FIR3 (HOPS-370; L$_{bol}$~314~L$_{\odot}$) at $\sim$0.1" (40~au) resolution for the continuum emission and ~0.25" (100 au) resolution of nine molecular lines. The dust continuum observed with ALMA at 0.87~mm and 1.3~mm resolve a near edge-on disk toward HOPS-370 with an apparent radius of ~100 au. The VLA observations detect both the disk in dust continuum and free-free emission extended along the jet direction. The ALMA observations of molecular lines (H$_2$CO, SO, CH$_3$OH, $^{13}$CO, C$^{18}$O, NS, and H$^{13}$CN) reveal rotation of the apparent disk surrounding HOPS-370 orthogonal to the jet/outflow direction. We fit radiative transfer models to both the dust continuum structure of the disk and molecular line kinematics of the inner envelope and disk for the H$_2$CO, CH$_3$OH, NS, and SO lines. The central protostar mass is determined to be $\sim$2.5 M_sun with a disk radius of $\sim$94~au, when fit using combinations of the H$_2$CO, CH$_3$OH, NS, and SO lines, consistent with an intermediate-mass protostar. Modeling of the dust continuum and spectral energy distribution (SED) yields a disk mass of 0.035~M$_{\odot}$ (inferred dust+gas) and a dust disk radius of 62~au, thus the dust disk may have a smaller radius than the gas disk, similar to Class II disks. In order to explain the observed luminosity with the measured protostar mass, HOPS-370 must be accreting at a rate between 1.7 and 3.2$\times$10$^{-5}$~M$_{\odot}$~yr$^{-1}$.
△ Less
Submitted 10 November, 2020; v1 submitted 2 November, 2020;
originally announced November 2020.
-
A Search for In-Situ Field OB Star Formation in the Small Magellanic Cloud
Authors:
Irene Vargas-Salazar,
M. S. Oey,
Jesse R. Barnes,
Xinyi Chen,
N. Castro,
Kaitlin M. Kratter,
Timothy A. Faerber
Abstract:
Whether any OB stars form in isolation is a question central to theories of massive star formation. To address this, we search for tiny, sparse clusters around 210 field OB stars from the Runaways and Isolated O-Type Star Spectroscopic Survey of the SMC (RIOTS4), using friends-of-friends (FOF) and nearest neighbors (NN) algorithms. We also stack the target fields to evaluate the presence of an agg…
▽ More
Whether any OB stars form in isolation is a question central to theories of massive star formation. To address this, we search for tiny, sparse clusters around 210 field OB stars from the Runaways and Isolated O-Type Star Spectroscopic Survey of the SMC (RIOTS4), using friends-of-friends (FOF) and nearest neighbors (NN) algorithms. We also stack the target fields to evaluate the presence of an aggregate density enhancement. Using several statistical tests, we compare these observations with three random-field datasets, and we also compare the known runaways to non-runaways. We find that the local environments of non-runaways show higher aggregate central densities than for runaways, implying the presence of some "tips-of-iceberg" (TIB) clusters. We find that the frequency of these tiny clusters is low, $\sim 4-5\%$ of our sample. This fraction is much lower than some previous estimates, but is consistent with field OB stars being almost entirely runaway and walkaway stars. The lack of TIB clusters implies that such objects either evaporate on short timescales, or do not form, implying a higher cluster lower-mass limit and consistent with a relationship between maximum stellar mass ($m_{\rm max}$) and the mass of the cluster ($M_{\rm cl}$). On the other hand, we also cannot rule out that some OB stars may form in highly isolated conditions. Our results set strong constraints on the formation of massive stars in relative isolation.
△ Less
Submitted 25 September, 2020;
originally announced September 2020.
-
Open Cluster Chemical Homogeneity Throughout the Milky Way
Authors:
Vijith Jacob Poovelil,
G. Zasowski,
S. Hasselquist,
A. Seth,
John Donor,
Rachael L. Beaton,
K. Cunha,
Peter M. Frinchaboy,
D. A. García-Hernández,
K. Hawkins,
K. M. Kratter,
Richard R. Lane,
C. Nitschelm
Abstract:
The chemical homogeneity of surviving stellar clusters contains important clues about interstellar medium (ISM) mixing efficiency, star formation, and the enrichment history of the Galaxy. Existing measurements in a handful of open clusters suggest homogeneity in several elements at the 0.03 dex level. Here we present (i) a new cluster member catalog based only on APOGEE radial velocities and Gaia…
▽ More
The chemical homogeneity of surviving stellar clusters contains important clues about interstellar medium (ISM) mixing efficiency, star formation, and the enrichment history of the Galaxy. Existing measurements in a handful of open clusters suggest homogeneity in several elements at the 0.03 dex level. Here we present (i) a new cluster member catalog based only on APOGEE radial velocities and Gaia-DR2 proper motions, (ii) improved abundance uncertainties for APOGEE cluster members, and (iii) the dependence of cluster homogeneity on Galactic and cluster properties, using abundances of eight elements from the APOGEE survey for ten high-quality clusters. We find that cluster homogeneity is uncorrelated with Galactocentric distance, |Z|, age, and metallicity. However, velocity dispersion, which is a proxy for cluster mass, is positively correlated with intrinsic scatter at relatively high levels of significance for [Ca/Fe] and [Mg/Fe]. We also see a possible positive correlation at a low level of significance for [Ni/Fe], [Si/Fe], [Al/Fe], and [Fe/H], while [Cr/Fe] and [Mn/Fe] are uncorrelated. The elements that show a correlation with velocity dispersion are those that are predominantly produced by core-collapse supernovae (CCSNe). However, the small sample size and relatively low correlation significance highlight the need for follow-up studies. If borne out by future studies, these findings would suggest a quantitative difference between the correlation lengths of elements produced predominantly by Type~Ia SNe versus CCSNe, which would have implications for Galactic chemical evolution models and the feasibility of chemical tagging.
△ Less
Submitted 14 September, 2020;
originally announced September 2020.
-
A disk-driven resonance as the origin of high inclinations of close-in planets
Authors:
Cristobal Petrovich,
Diego J. Muñoz,
Kaitlin M. Kratter,
Renu Malhotra
Abstract:
The recent characterization of transiting close-in planets has revealed an intriguing population of sub-Neptunes with highly tilted and even polar orbits relative to their host star's equator. Any viable theory for the origin of these close-in, polar planets must explain (1) the observed stellar obliquities, (2) the substantial eccentricities, and (3) the existence of Jovian companions with large…
▽ More
The recent characterization of transiting close-in planets has revealed an intriguing population of sub-Neptunes with highly tilted and even polar orbits relative to their host star's equator. Any viable theory for the origin of these close-in, polar planets must explain (1) the observed stellar obliquities, (2) the substantial eccentricities, and (3) the existence of Jovian companions with large mutual inclinations. In this work, we propose a theoretical model that satisfies these requirements without invoking tidal dissipation or large primordial inclinations. Instead, tilting is facilitated by the protoplanetary disk dispersal during the late stage of planet formation, initiating a process of resonance sweeping and parametric instability. This mechanism consists of two steps. First, a nodal secular resonance excites the inclination to large values; then, once the inclination reaches a critical value, a linear eccentric instability is triggered, which detunes the resonance and ends inclination growth. The critical inclination is pushed to high values by general relativistic precession, making polar orbits an inherently post-Newtonian outcome. Our model predicts that polar, close-in sub-Neptunes coexist with cold Jupiters in low stellar obliquity orbits.
△ Less
Submitted 19 August, 2020;
originally announced August 2020.
-
The Close Binary Fraction as a Function of Stellar Parameters in APOGEE: A Strong Anti-Correlation With $α$ Abundances
Authors:
Christine N. Mazzola,
Carles Badenes,
Maxwell Moe,
Sergey E. Koposov,
Marina Kounkel,
Kaitlin Kratter,
Kevin Covey,
Matthew G. Walker,
Todd A. Thompson,
Brett Andrews,
Peter E. Freeman,
Borja Anguiano,
Joleen K. Carlberg,
Nathan M. De Lee,
Peter M. Frinchaboy,
Hannah M. Lewis,
Steven Majewski,
David Nidever,
Christian Nitschelm,
Adrian M. Price-Whelan,
Alexandre Roman-Lopes,
Keivan G. Stassun,
Nicholas W. Troup
Abstract:
We use observations from the APOGEE survey to explore the relationship between stellar parameters and multiplicity. We combine high-resolution repeat spectroscopy for 41,363 dwarf and subgiant stars with abundance measurements from the APOGEE pipeline and distances and stellar parameters derived using \textit{Gaia} DR2 parallaxes from \cite{Sanders2018} to identify and characterise stellar multipl…
▽ More
We use observations from the APOGEE survey to explore the relationship between stellar parameters and multiplicity. We combine high-resolution repeat spectroscopy for 41,363 dwarf and subgiant stars with abundance measurements from the APOGEE pipeline and distances and stellar parameters derived using \textit{Gaia} DR2 parallaxes from \cite{Sanders2018} to identify and characterise stellar multiples with periods below 30 years, corresponding to \drvm$\gtrsim$ 3 \kms, where \drvm\ is the maximum APOGEE-detected shift in the radial velocities. Chemical composition is responsible for most of the variation in the close binary fraction in our sample, with stellar parameters like mass and age playing a secondary role. In addition to the previously identified strong anti-correlation between the close binary fraction and \feh\, we find that high abundances of $α$ elements also suppress multiplicity at most values of \feh\ sampled by APOGEE. The anti-correlation between $α$ abundances and multiplicity is substantially steeper than that observed for Fe, suggesting C, O, and Si in the form of dust and ices dominate the opacity of primordial protostellar disks and their propensity for fragmentation via gravitational stability. Near \feh{} = 0 dex, the bias-corrected close binary fraction ($a<10$ au) decreases from $\approx$ 100 per cent at \alh{} = $-$0.2 dex to $\approx$ 15 per cent near \alh{} = 0.08 dex, with a suggestive turn-up to $\approx$20 per cent near \alh{} = 0.2. We conclude that the relationship between stellar multiplicity and chemical composition for sun-like dwarf stars in the field of the Milky Way is complex, and that this complexity should be accounted for in future studies of interacting binaries.
△ Less
Submitted 15 September, 2020; v1 submitted 17 July, 2020;
originally announced July 2020.
-
The Carina High-Contrast Imaging Project for massive Stars (CHIPS) I. Methodology and proof of concept on QZ Car ($\equiv$ HD93206)
Authors:
A. Rainot,
M. Reggiani,
H. Sana,
J. Bodensteiner,
C. A. Gomez-Gonzalez,
O. Absil,
V. Christiaens,
P. Delorme,
L. A. Almeida,
S. Caballero-Nieves,
J. De Ridder,
K. Kratter,
S. Lacour,
J. -B. Le Bouquin,
L. Pueyo,
H. Zinnecker
Abstract:
Massive stars like company. However, low-mass companions have remained extremely difficult to detect at angular separations ($ρ$) smaller than 1" (~1000-3000 au considering typical distance to nearby massive stars) given the large brightness contrast between the companion and the central star. Constraints on the low-mass end of the companions mass-function for massive stars are however needed, for…
▽ More
Massive stars like company. However, low-mass companions have remained extremely difficult to detect at angular separations ($ρ$) smaller than 1" (~1000-3000 au considering typical distance to nearby massive stars) given the large brightness contrast between the companion and the central star. Constraints on the low-mass end of the companions mass-function for massive stars are however needed, for example to help distinguishing between various scenarios for the formation of massive stars. To obtain statistically significant constraint on the presence of low-mass companions beyond the typical detection limit of current surveys ($Δ\mathrm{mag} \lesssim 5$ at $ρ\lesssim 1$"), we initiated a survey of O and Wolf-Rayet stars in the Carina region using the SPHERE coronagraphic instrument on the VLT. In this first paper, we aim to introduce the survey, to present the methodology and to demonstrate the capability of SPHERE for massive stars using the multiple system QZ~Car. High-contrast imaging techniques, such as angular- and spectral-differential imaging techniques as well as PSF-fitting, were applied to detect and measure the relative flux of companions in each spectral channel of the instrument. We detected 19 sources around the QZ~Car system with detection limits of 9~mag at $ρ> 200$~mas for IFS and as faint as 13~mag at $ρ$ > 1".8 for IRDIS (corresponding to sub-solar masses for potential companions). All but two are reported here for the first time. Based on this proof of concept, we showed that VLT/SPHERE allows us to reach the sub-solar mass regime of the companion mass function. This paves the way for this type of observation with a large sample of massive stars to provide novel constraints on the multiplicity of massive stars in a region of the parameter space that has remained inaccessible so far.
△ Less
Submitted 2 June, 2020;
originally announced June 2020.
-
Dust Settling Instability in Protoplanetary Discs
Authors:
Leonardo Krapp,
Andrew N. Youdin,
Kaitlin M. Kratter,
Pablo Benítez-Llambay
Abstract:
The streaming instability (SI) has been extensively studied in the linear and non-linear regimes as a mechanism to concentrate solids and trigger planetesimal formation in the midplane of protoplanetary discs. A related dust settling instability (DSI) applies to particles while settling towards the midplane. The DSI has previously been studied in the linear regime, with predictions that it could t…
▽ More
The streaming instability (SI) has been extensively studied in the linear and non-linear regimes as a mechanism to concentrate solids and trigger planetesimal formation in the midplane of protoplanetary discs. A related dust settling instability (DSI) applies to particles while settling towards the midplane. The DSI has previously been studied in the linear regime, with predictions that it could trigger particle clumping away from the midplane. This work presents a range of linear calculations and non-linear simulations, performed with FARGO3D, to assess conditions for DSI growth. We expand on previous linear analyses by including particle size distributions and performing a detailed study of the amount of background turbulence needed to stabilize the DSI. When including binned size distributions, the DSI often produces converged growth rates with fewer bins than the standard SI. With background turbulence, we find that the most favorable conditions for DSI growth are weak turbulence, characterized by $α\lesssim 10^{-6}$ with intermediate-sized grains that settle from one gas scale-height. These conditions could arise during a sudden decrease in disc turbulence following an accretion outburst. Ignoring background turbulence, we performed a parameter survey of local 2D DSI simulations. Particle clumping was either weak or occurred slower than particles settle. Clumping was reduced by a factor of two in a comparison 3D simulation. Overall, our results strongly disfavor the hypothesis that the DSI significantly promotes planetesimal formation. Non-linear simulations of the DSI with different numerical methods could support or challenge these findings.
△ Less
Submitted 15 September, 2020; v1 submitted 9 April, 2020;
originally announced April 2020.
-
The highly variable time evolution of star-forming cores identified with dendrograms
Authors:
Rachel A. Smullen,
Kaitlin M. Kratter,
Stella S. R. Offner,
Aaron T. Lee,
Hope How-Huan Chen
Abstract:
We investigate the time evolution of dense cores identified in molecular cloud simulations using dendrograms, which are a common tool to identify hierarchical structure in simulations and observations of star formation. We develop an algorithm to link dendrogram structures through time using the three-dimensional density field from magnetohydrodynamical simulations, thus creating histories for all…
▽ More
We investigate the time evolution of dense cores identified in molecular cloud simulations using dendrograms, which are a common tool to identify hierarchical structure in simulations and observations of star formation. We develop an algorithm to link dendrogram structures through time using the three-dimensional density field from magnetohydrodynamical simulations, thus creating histories for all dense cores in the domain. We find that the population-wide distributions of core properties are relatively invariant in time, and quantities like the core mass function match with observations. Despite this consistency, an individual core may undergo large (>40%), stochastic variations due to the redefinition of the dendrogram structure between timesteps. This variation occurs independent of environment and stellar content. We identify a population of short-lived (<200 kyr) overdensities masquerading as dense cores that may comprise ~20% of any time snapshot. Finally, we note the importance of considering the full history of cores when interpreting the origin of the initial mass function; we find that, especially for systems containing multiple stars, the core mass defined by a dendrogram leaf in a snapshot is typically less than the final system stellar mass. This work reinforces that there is no time-stable density contour that defines a star-forming core. The dendrogram itself can induce significant structure variation between timesteps due to small changes in the density field. Thus, one must use caution when comparing dendrograms of regions with different ages or environment properties because differences in dendrogram structure may not come solely from the physical evolution of dense cores.
△ Less
Submitted 2 April, 2020;
originally announced April 2020.
-
The Rich Lack Close Neighbours: The Dependence of Blue-Straggler Fraction on Metallicity
Authors:
Rosemary F. G. Wyse,
Maxwell Moe,
Kaitlin M. Kratter
Abstract:
Blue straggler stars (BSS) have been identified in star clusters and in field populations in our own Milky Way galaxy and in its satellite galaxies. They manifest as stars bluer and more luminous than the dominant old population, and usually have a spatial distribution that follows the old population. Their progenitors are likely to have been close binaries. We investigate trends of the BSS popula…
▽ More
Blue straggler stars (BSS) have been identified in star clusters and in field populations in our own Milky Way galaxy and in its satellite galaxies. They manifest as stars bluer and more luminous than the dominant old population, and usually have a spatial distribution that follows the old population. Their progenitors are likely to have been close binaries. We investigate trends of the BSS population in dwarf spheroidal galaxies (dSph) and in the bulge of the Milky Way and find an anti-correlation between the relative frequency of BSS and the metallicity of the parent population. The rate of occurrence of BSS in the metal-poor dwarf galaxies is approximately twice that found in the solar-metallicity bulge population. This trend of decreasing relative population of BSS with increasing metallicity mirrors that found for the close-binary fraction in the field population of the Milky Way. We argue that the dominant mode of BSS formation in low-density environments is likely to be mass transfer in close-binary systems. It then follows that the similarity between the trends for BSS in the dSph and field stars in our Galaxy supports the proposal that the small-scale fragmentation during star formation is driven by the same dominant physical process, despite the diversity in environments, plausibly gravitational instability of proto-stellar discs.
△ Less
Submitted 31 March, 2020;
originally announced April 2020.
-
The VLA/ALMA Nascent Disk and Multiplicity (VANDAM) Survey of Orion Protostars. A Statistical Characterization of Class 0 and I Protostellar Disks
Authors:
John J. Tobin,
Patrick Sheehan,
S. Thomas Megeath,
Ana Karla Diaz-Rodriguez,
Stella S. R. Offner,
Nadia M. Murillo,
Merel van 't Hoff,
Ewine F. van Dishoeck,
Mayra Osorio,
Guillem Anglada,
Elise Furlan,
Amelia M. Stutz,
Nickalas Reynolds,
Nicole Karnath,
William J. Fischer,
Magnus Persson,
Leslie W. Looney,
Zhi-Yun Li,
Ian Stephens,
Claire J. Chandler,
Erin Cox,
Michael M. Dunham,
Lukasz Tychoniec,
Mihkel Kama,
Kaitlin Kratter
, et al. (11 additional authors not shown)
Abstract:
We have conducted a survey of 328 protostars in the Orion molecular clouds with ALMA at 0.87 mm at a resolution of $\sim$0.1" (40 au), including observations with the VLA at 9 mm toward 148 protostars at a resolution of $\sim$0.08" (32 au). This is the largest multi-wavelength survey of protostars at this resolution by an order of magnitude. We use the dust continuum emission at 0.87 mm and 9 mm t…
▽ More
We have conducted a survey of 328 protostars in the Orion molecular clouds with ALMA at 0.87 mm at a resolution of $\sim$0.1" (40 au), including observations with the VLA at 9 mm toward 148 protostars at a resolution of $\sim$0.08" (32 au). This is the largest multi-wavelength survey of protostars at this resolution by an order of magnitude. We use the dust continuum emission at 0.87 mm and 9 mm to measure the dust disk radii and masses toward the Class 0, Class I, and Flat Spectrum protostars, characterizing the evolution of these disk properties in the protostellar phase. The mean dust disk radii for the Class 0, Class I, and Flat Spectrum protostars are 44.9$^{+5.8}_{-3.4}$, 37.0$^{+4.9}_{-3.0}$, and 28.5$^{+3.7}_{-2.3}$ au, respectively, and the mean protostellar dust disk masses are 25.9$^{+7.7}_{-4.0}$, 14.9$^{+3.8}_{-2.2}$, 11.6$^{+3.5}_{-1.9}$ Earth masses, respectively. The decrease in dust disk masses is expected from disk evolution and accretion, but the decrease in disk radii may point to the initial conditions of star formation not leading to the systematic growth of disk radii or that radial drift is keeping the dust disk sizes small. At least 146 protostellar disks (35% out of 379 detected 0.87 mm continuum sources plus 42 non-detections) have disk radii greater than 50 au in our sample. These properties are not found to vary significantly between different regions within Orion. The protostellar dust disk mass distributions are systematically larger than that of Class II disks by a factor of $>$4, providing evidence that the cores of giant planets may need to at least begin their formation during the protostellar phase.
△ Less
Submitted 13 January, 2020;
originally announced January 2020.
-
Astro2020 APC White Paper: Theoretical Astrophysics 2020-2030
Authors:
Juna A. Kollmeier,
Lauren Anderson,
Andrew Benson,
Tamara Bogdanovic,
Michael Boylan-Kolchin,
James S. Bullock,
Romeel Dave,
Federico Fraschetti,
Jim Fuller,
Philip F. Hopkins,
Manoj Kaplinghat,
Kaitlin Kratter,
Astrid Lamberts,
M. Coleman Miller,
James E. Owen,
E. Sterl Phinney,
Anthony L. Piro,
Hans-Walter Rix,
Brant Robertson,
Andrew Wetzel,
Coral Wheeler,
Andrew N. Youdin,
Matias Zaldarriaga
Abstract:
The past two decades have seen a tremendous investment in observational facilities that promise to reveal new and unprecedented discoveries about the universe. In comparison, the investment in theoretical work is completely dwarfed, even though theory plays a crucial role in the interpretation of these observations, predicting new types of phenomena, and informing observing strategies. In this whi…
▽ More
The past two decades have seen a tremendous investment in observational facilities that promise to reveal new and unprecedented discoveries about the universe. In comparison, the investment in theoretical work is completely dwarfed, even though theory plays a crucial role in the interpretation of these observations, predicting new types of phenomena, and informing observing strategies. In this white paper, we argue that in order to reach the promised critical breakthroughs in astrophysics over the next decade and well beyond, the national agencies must take a serious approach to investment in theoretical astrophysics research. We discuss the role of theory in shaping our understanding of the universe, and then we provide a multi-level strategy, from the grassroots to the national, to address the current underinvestment in theory relative to observational work.
△ Less
Submitted 18 December, 2019;
originally announced December 2019.
-
Dynamics of small grains in transitional discs
Authors:
Mark R. Krumholz,
Michael J. Ireland,
Kaitlin M. Kratter
Abstract:
Transitional discs have central regions characterised by significant depletion of both dust and gas compared to younger, optically-thick discs. However, gas and dust are not depleted by equal amounts: gas surface densities are typically reduced by factors of $\sim 100$, but small dust grains are sometimes depleted by far larger factors, to the point of being undetectable. While this extreme dust d…
▽ More
Transitional discs have central regions characterised by significant depletion of both dust and gas compared to younger, optically-thick discs. However, gas and dust are not depleted by equal amounts: gas surface densities are typically reduced by factors of $\sim 100$, but small dust grains are sometimes depleted by far larger factors, to the point of being undetectable. While this extreme dust depletion is often attributed to planet formation, in this paper we show that another physical mechanism is possible: expulsion of grains from the disc by radiation pressure. We explore this mechanism using 2D simulations of dust dynamics, simultaneously solving the equation of radiative transfer with the evolution equations for dust diffusion and advection under the combined effects of stellar radiation and hydrodynamic interaction with a turbulent, accreting background gas disc. We show that, in transition discs that are depleted in both gas and dust fraction by factors of $\sim 100-1000$ compared to minimum mass Solar nebular values, and where the ratio of accretion rate to stellar luminosity is low ($\dot{M}/L \lesssim 10^{-10}$ $M_\odot$ yr$^{-1}$ $L_\odot^{-1}$), radiative clearing of any remaining $\sim 0.5$ $μ$m and larger grains is both rapid and inevitable. The process is size-dependent, with smaller grains removed fastest and larger ones persisting for longer times. Our proposed mechanism thus naturally explains the extreme depletion of small grains commonly-found in transition discs. We further suggest that the dependence of this mechanism on grain size and optical properties may explain some of the unusual grain properties recently discovered in a number of transition discs. The simulation code we develop is freely available.
△ Less
Submitted 19 August, 2020; v1 submitted 14 December, 2019;
originally announced December 2019.
-
Impact of Binary Stars on Planet Statistics -- I. Planet Occurrence Rates, Trends with Stellar Mass, and Wide Companions to Hot Jupiter Hosts
Authors:
Maxwell Moe,
Kaitlin M. Kratter
Abstract:
Close binaries suppress the formation of circumstellar (S-type) planets and therefore significantly bias the inferred planet occurrence rates and statistical trends. After compiling various radial velocity and high-resolution imaging surveys, we determine that binaries with a < 1 au fully suppress S-type planets, binaries with a = 10 au host close planets at 15% the occurrence rate of single stars…
▽ More
Close binaries suppress the formation of circumstellar (S-type) planets and therefore significantly bias the inferred planet occurrence rates and statistical trends. After compiling various radial velocity and high-resolution imaging surveys, we determine that binaries with a < 1 au fully suppress S-type planets, binaries with a = 10 au host close planets at 15% the occurrence rate of single stars, and wide binaries with a > 200 au have a negligible effect on close planet formation. We show that F = 43% +/- 7% of solar-type primaries do not host close planets due to suppression by close stellar companions. By removing spectroscopic binaries from their samples, radial velocity surveys for giant planets boost their detection rates by a factor of 1/(1-F) = 1.8 +/- 0.2 compared to transiting surveys. This selection bias fully accounts for the discrepancy in hot Jupiter occurrence rates inferred from these two detection methods. Correcting for both planet suppression by close binaries and transit dilution by wide binaries, the occurrence rate of small planets orbiting single G-dwarfs is 2.1 +/- 0.3 times larger than the rate inferred from all G-dwarfs in the Kepler survey. Additionally, about half (but not all) of the observed increase in small, short-period planets toward low-mass hosts can be explained by the corresponding decrease in the binary fraction. Finally, we demonstrate that the apparent enhancement of wide stellar companions to hot Jupiter hosts is due to multiple selection effects. Very close binaries, brown dwarf companions, and massive planets with M_2 > 7 M_J within a < 0.2 au preferentially have metal-poor hosts and wide tertiary companions, but genuine hot Jupiters with M_p = 0.2-4 M_J that formed via core accretion instead favor metal-rich hosts and do not exhibit a statistically significant excess of wide stellar companions.
△ Less
Submitted 30 September, 2020; v1 submitted 3 December, 2019;
originally announced December 2019.
-
The Formation and Evolution of Wide-Orbit Stellar Multiples In Magnetized Clouds
Authors:
Aaron T. Lee,
Stella S. R. Offner,
Kaitlin M. Kratter,
Rachel A. Smullen,
Pak Shing Li
Abstract:
Stars rarely form in isolation. Nearly half of the stars in the Milky Way have a companion, and this fraction increases in star-forming regions. However, why some dense cores and filaments form bound pairs while others form single stars remains unclear. We present a set of three-dimensional, gravo-magnetohydrodynamic simulations of turbulent star-forming clouds, aimed at understanding the formatio…
▽ More
Stars rarely form in isolation. Nearly half of the stars in the Milky Way have a companion, and this fraction increases in star-forming regions. However, why some dense cores and filaments form bound pairs while others form single stars remains unclear. We present a set of three-dimensional, gravo-magnetohydrodynamic simulations of turbulent star-forming clouds, aimed at understanding the formation and evolution of multiple-star systems formed through large scale (>~$10^3$ AU) turbulent fragmentation. We investigate three global magnetic field strengths, with global mass-to-flux ratios of $μ_φ$=2, 8, and 32. The initial separations of protostars in multiples depends on the global magnetic field strength, with stronger magnetic fields (e.g., $μ_φ$=2) suppressing fragmentation on smaller scales. The overall multiplicity fraction (MF) is between 0.4-0.6 for our strong and intermediate magnetic field strengths, which is in agreement with observations. The weak field case has a lower fraction. The MF is relatively constant throughout the simulations, even though stellar densities increase as collapse continues. While the MF rarely exceeds 60% in all three simulations, over 80% of all protostars are part of a binary system at some point. We additionally find that the distribution of binary spin mis-alignment angles is consistent with a randomized distribution. In all three simulations, several binaries originate with wide separations and dynamically evolve to <~ $10^2$ AU separations. We show that a simple model of mass accretion and dynamical friction with the gas can explain this orbital evolution.
△ Less
Submitted 18 November, 2019;
originally announced November 2019.
-
Circumbinary Accretion from Finite and Infinite Disks
Authors:
Diego Muñoz,
Dong Lai,
Kaitlin Kratter,
Ryan Miranda
Abstract:
We carry out 2D viscous hydrodynamics simulations of circumbinary disk (CBD) accretion using {\footnotesize AREPO}. We resolve the accretion flow from a large-scale CBD down to the streamers and disks around individual binary components. Extending our recent studies \citep{mun19}, we consider circular binaries with various mass ratios ($0.1\leq q_{\rm{b}}\leq1$) and study accretion from ``infinite…
▽ More
We carry out 2D viscous hydrodynamics simulations of circumbinary disk (CBD) accretion using {\footnotesize AREPO}. We resolve the accretion flow from a large-scale CBD down to the streamers and disks around individual binary components. Extending our recent studies \citep{mun19}, we consider circular binaries with various mass ratios ($0.1\leq q_{\rm{b}}\leq1$) and study accretion from ``infinite'', steady-supply disks and from finite-sized, viscously spreading tori. For ``infinite'' disks, a global steady state can be reached, and the accretion variability has a dominant frequency ${\sim}0.2Ω_{\rm{b}}$ for $q_{\rm{b}}>0.5$ and $Ω_{\rm{b}}$ for $q_{\rm{b}}<0.5$, ($Ω_{\rm{b}}$ is the binary angular frequency). We find that the accretion ``eigenvalue'' $l_0$ -- the net angular momentum transfer from the disk to the binary per unit accreted mass -- is always positive and falls in the range ($0.65$-$0.85)a_{\rm b}^2Ω_{\rm{b}}$ (with $a_{\rm{b}}$ the binary separation), depending weakly on the mass ratio and viscosity. This leads to binary expansion when $q_{\rm{b}}\gtrsim0.3$. Accretion from a finite torus can be separated into two phases: an initial transient phase, corresponding to the filling of the binary cavity, followed by a viscous pseudo-stationary phase, during which the torus viscously spreads and accretes onto the binary. In the viscous phase, the net torque on the binary per unit accreted mass is close to $l_0$, the value derived for ``infinite'' disks. We conclude that similar-mass binaries accreting from CBDs gain angular momentum and expand over long time scales. This result significantly impacts the coalescence of supermassive binary black holes and newly formed binary stars. We offer a word of caution against conclusions drawn from simulations of transient accretion onto empty circumbinary cavities.
△ Less
Submitted 16 December, 2019; v1 submitted 10 October, 2019;
originally announced October 2019.
-
The VLA/ALMA Nascent Disk and Multiplicity (VANDAM) Survey of Orion Protostars I. Identifying and Characterizing the Protostellar Content of the OMC2-FIR4 and OMC2-FIR3 Regions
Authors:
John J. Tobin,
S. Thomas Megeath,
Merel van 't Hoff,
Ana Karla Diaz-Rodriguez,
Nickalas Reynolds,
Mayra Osorio,
Guillem Anglada,
Elise Furlan,
Nicole Karnath,
Stella S. R. Offner,
Patrick Sheehan,
Sarah I. Sadavoy,
Amelia M. Stutz,
William J. Fischer,
Mihkel Kama,
Magnus Persson,
James Di Francesco,
Leslie W. Looney,
Dan M. Watson,
Zhi-Yun Li,
Ian Stephens,
Claire J. Chandler,
Erin Cox,
Michael M. Dunham,
Kaitlin Kratter
, et al. (9 additional authors not shown)
Abstract:
We present ALMA (0.87~mm) and VLA (9~mm) observations toward OMC2-FIR4 and OMC2-FIR3 within the Orion integral-shaped filament that are thought to be the nearest regions of intermediate mass star formation. We characterize the continuum sources within these regions on $\sim$40~AU (0\farcs1) scales and associated molecular line emission at a factor of $\sim$30 better resolution than previous observ…
▽ More
We present ALMA (0.87~mm) and VLA (9~mm) observations toward OMC2-FIR4 and OMC2-FIR3 within the Orion integral-shaped filament that are thought to be the nearest regions of intermediate mass star formation. We characterize the continuum sources within these regions on $\sim$40~AU (0\farcs1) scales and associated molecular line emission at a factor of $\sim$30 better resolution than previous observations at similar wavelengths. We identify six compact continuum sources within OMC2-FIR4, four in OMC2-FIR3, and one additional source just outside OMC2-FIR4. This continuum emission is tracing the inner envelope and/or disk emission on less than 100~AU scales. HOPS-108 is the only protostar in OMC2-FIR4 that exhibits emission from high-excitation transitions of complex organic molecules (e.g., methanol and other lines) coincident with the continuum emission. HOPS-370 in OMC2-FIR3 with L~$\sim$~360~\lsun, also exhibits emission from high-excitation methanol and other lines. The methanol emission toward these two protostars is indicative of temperatures high enough to thermally evaporate methanol from icy dust grains; overall these protostars have characteristics similar to hot corinos. We do not identify a clear outflow from HOPS-108 in \twco, but find evidence of interaction between the outflow/jet from HOPS-370 and the OMC2-FIR4 region. The multitude of observational constraints indicate that HOPS-108 is likely a low to intermediate-mass protostar in its main mass accretion phase and it is the most luminous protostar in OMC2-FIR4. The high resolution data presented here are essential for disentangling the embedded protostars from their surrounding dusty environments and characterizing them.
△ Less
Submitted 1 October, 2019;
originally announced October 2019.
-
Instabilities in Multi-Planet Circumbinary Systems
Authors:
Adam P. Sutherland,
Kaitlin M. Kratter
Abstract:
The majority of the discovered transiting circumbinary planets are located very near the innermost stable orbits permitted, raising questions about the origins of planets in such perturbed environments. Most favored formation scenarios invoke formation at larger distances and subsequent migration to their current locations. Disk-driven planet migration in multi-planet systems is likely to trap pla…
▽ More
The majority of the discovered transiting circumbinary planets are located very near the innermost stable orbits permitted, raising questions about the origins of planets in such perturbed environments. Most favored formation scenarios invoke formation at larger distances and subsequent migration to their current locations. Disk-driven planet migration in multi-planet systems is likely to trap planets in mean motion resonances and drive planets inward into regions of larger dynamical perturbations from the binary. We demonstrate how planet-planet resonances can interact with the binary through secular forcing and mean-motion resonances, driving chaos in the system. We show how this chaos will shape the architecture of circumbinary systems, with specific applications to Kepler 47 and the Pluto-Charon system, limiting maximum possible stable eccentricities and indicating what resonances are likely to exist. We are also able to constrain the minimum migration rates of resonant circumbinary planets.
△ Less
Submitted 29 May, 2019;
originally announced May 2019.
-
Astro2020 Science White Paper: The Formation and Evolution of Multiple Star Systems
Authors:
John J. Tobin,
Marina Kounkel,
Stella Offner,
Patrick Sheehan,
Doug Johnstone,
S. Thomas Megeath,
Kaitlin M. Kratter,
Ian Stephens,
Zhi-Yun Li,
Sarah Sadavoy,
Leslie Looney,
Joel Green,
Rob Gutermuth,
Will Fischer,
Michael M. Dunham,
Yao-Lun Yang
Abstract:
Significant advances have been made over the past decade in the characterization of multiple protostar systems, enabled by the Karl G. Jansky Very Large Array (VLA), high-resolution infrared observations with the Hubble Space Telescope, and ground-based facilities. To further understand the mechanism(s) of multiple star formation, a combination of statistics, high-angular resolution radio/millimet…
▽ More
Significant advances have been made over the past decade in the characterization of multiple protostar systems, enabled by the Karl G. Jansky Very Large Array (VLA), high-resolution infrared observations with the Hubble Space Telescope, and ground-based facilities. To further understand the mechanism(s) of multiple star formation, a combination of statistics, high-angular resolution radio/millimeter continuum imaging, characterization of kinematic structure, magnetic fields via polarimetry, and comparison with numerical simulations are needed. Thus, understanding the origin of stellar multiplicity in different regimes of companion separation will soon be within reach. However, to overcome challenges that studies in this field are now confronted with, a range of new capabilities are required: a new millimeter/centimeter wave facility with 10 mas resolution at λ=1 cm, space-based near to far-infrared observatories, continued development of low to high resolution spectroscopy on 3m to 10m class telescopes, and an ELT-class telescope with near to mid-infrared imaging/spectroscopic capability.
△ Less
Submitted 17 April, 2019;
originally announced April 2019.
-
On the Mass Function, Multiplicity, and Origins of Wide-Orbit Giant Planets
Authors:
Kevin Wagner,
Dániel Apai,
Kaitlin M. Kratter
Abstract:
A major outstanding question regarding the formation of planetary systems is whether wide-orbit giant planets form differently than close-in giant planets. We aim to establish constraints on two key parameters that are relevant for understanding the formation of wide-orbit planets: 1) the relative mass function and 2) the fraction of systems hosting multiple companions. In this study, we focus on…
▽ More
A major outstanding question regarding the formation of planetary systems is whether wide-orbit giant planets form differently than close-in giant planets. We aim to establish constraints on two key parameters that are relevant for understanding the formation of wide-orbit planets: 1) the relative mass function and 2) the fraction of systems hosting multiple companions. In this study, we focus on systems with directly imaged substellar companions, and the detection limits on lower-mass bodies within these systems. First, we uniformly derive the mass probability distributions of known companions. We then combine the information contained within the detections and detection limits into a survival analysis statistical framework to estimate the underlying mass function of the parent distribution. Finally, we calculate the probability that each system may host multiple substellar companions. We find that 1) the companion mass distribution is rising steeply toward smaller masses, with a functional form of $N\propto M^{-1.3\pm0.3}$, and consequently, 2) many of these systems likely host additional undetected sub-stellar companions. Combined, these results strongly support the notion that wide-orbit giant planets are formed predominantly via core accretion, similar to the better studied close-in giant planets. Finally, given the steep rise in the relative mass function with decreasing mass, these results suggest that future deep observations should unveil a greater number of directly imaged planets.
△ Less
Submitted 26 July, 2019; v1 submitted 12 April, 2019;
originally announced April 2019.
-
Close companions around young stars
Authors:
Marina Kounkel,
Kevin Covey,
Maxwell Moe,
Kaitlin M. Kratter,
Genaro Suárez,
Keivan G. Stassun,
Carlos Román-Zúñiga,
Jesus Hernandez,
Jinyoung Serena Kim,
Karla Peña Ramírez,
Alexandre Roman-Lopes,
Guy S Stringfellow,
Karl O Jaehnig,
Jura Borissova,
Benjamin Tofflemire,
Daniel Krolikowski,
Aaron Rizzuto,
Adam Kraus,
Carles Badenes,
Penélope Longa-Peña,
Yilen Gómez Maqueo Chew,
Rodolfo Barba,
David L. Nidever,
Cody Brown,
Nathan De Lee
, et al. (4 additional authors not shown)
Abstract:
Multiplicity is a fundamental property that is set early during stellar lifetimes, and it is a stringent probe of the physics of star formation. The distribution of close companions around young stars is still poorly constrained by observations. We present an analysis of stellar multiplicity derived from APOGEE-2 spectra obtained in targeted observations of nearby star-forming regions. This is the…
▽ More
Multiplicity is a fundamental property that is set early during stellar lifetimes, and it is a stringent probe of the physics of star formation. The distribution of close companions around young stars is still poorly constrained by observations. We present an analysis of stellar multiplicity derived from APOGEE-2 spectra obtained in targeted observations of nearby star-forming regions. This is the largest homogeneously observed sample of high-resolution spectra of young stars. We developed an autonomous method to identify double lined spectroscopic binaries (SB2s). Out of 5007 sources spanning the mass range of $\sim$0.05--1.5 \msun, we find 399 binaries, including both RV variables and SB2s. The mass ratio distribution of SB2s is consistent with a uniform for $q<0.95$ with an excess of twins with $q>0.95$. The period distribution is consistent with what has been observed in close binaries ($<10$ AU) in the evolved populations. Three systems are found to have $q\sim$0.1, with a companion located within the brown dwarf desert. There are not any strong trends in the multiplicity fraction (MF) as a function of cluster age from 1 to 100 Myr. There is a weak dependence on stellar density, with companions being most numerous at $Σ_*\sim30$ stars/pc$^{-2}$, and decreasing in more diffuse regions. Finally, disk-bearing sources are deficient in SB2s (but not RV variables) by a factor of $\sim$2; this deficit is recovered by the systems without disks. This may indicate a quick dispersal of disk material in short-period equal mass systems that is less effective in binaries with lower $q$.
△ Less
Submitted 27 March, 2019; v1 submitted 25 March, 2019;
originally announced March 2019.
-
Dynamics of a Massive Binary at Birth
Authors:
Yichen Zhang,
Jonathan C. Tan,
Kei E. I. Tanaka,
James M. De Buizer,
Mengyao Liu,
Maria T. Beltran,
Kaitlin Kratter,
Diego Mardones,
Guido Garay
Abstract:
Almost all massive stars have bound stellar companions, existing in binaries or higher-order multiples. While binarity is theorized to be an essential feature of how massive stars form, essentially all information about such properties is derived from observations of already formed stars, whose orbital properties may have evolved since birth. Little is known about binarity during formation stages.…
▽ More
Almost all massive stars have bound stellar companions, existing in binaries or higher-order multiples. While binarity is theorized to be an essential feature of how massive stars form, essentially all information about such properties is derived from observations of already formed stars, whose orbital properties may have evolved since birth. Little is known about binarity during formation stages. Here we report high angular resolution observations of 1.3 mm continuum and H30alpha recombination line emission, which reveal a massive protobinary with apparent separation of 180 au at the center of the massive star-forming region IRAS07299-1651. From the line-of-sight velocity difference of 9.5 km/s of the two protostars, the binary is estimated to have a minimum total mass of 18 solar masses, consistent with several other metrics, and maximum period of 570 years, assuming a circular orbit. The H30alpha line from the primary protostar shows kinematics consistent with rotation along a ring of radius of 12 au. The observations indicate that disk fragmentation at several hundred au may have formed the binary, and much smaller disks are feeding the individual protostars.
△ Less
Submitted 18 March, 2019;
originally announced March 2019.
-
A Model-Independent Mass and Moderate Eccentricity for $β$ Pic b
Authors:
Trent J. Dupuy,
Timothy D. Brandt,
Kaitlin M. Kratter,
Brendan P. Bowler
Abstract:
We use a cross-calibration of Hipparcos and Gaia DR2 astrometry for $β$ Pic to measure the mass of the giant planet $β$ Pic b $(13\pm3$ $M_{\rm Jup})$ in a comprehensive joint orbit analysis that includes published relative astrometry and radial velocities. Our mass uncertainty is somewhat higher than previous work because our astrometry from the Hipparcos-Gaia Catalog of Accelerations accounts fo…
▽ More
We use a cross-calibration of Hipparcos and Gaia DR2 astrometry for $β$ Pic to measure the mass of the giant planet $β$ Pic b $(13\pm3$ $M_{\rm Jup})$ in a comprehensive joint orbit analysis that includes published relative astrometry and radial velocities. Our mass uncertainty is somewhat higher than previous work because our astrometry from the Hipparcos-Gaia Catalog of Accelerations accounts for the error inflation and systematic terms that are required to bring the two data sets onto a common astrometric reference frame, and because we fit freely for the host-star mass $(1.84\pm0.05$ $M_{\odot})$. This first model-independent mass for a directly imaged planet is inconsistent with cold-start models given the age of the $β$ Pic moving group $(22\pm6$ Myr) but consistent with hot- and warm-start models, concordant with past work. We find a higher eccentricity $(0.24\pm0.06)$ for $β$ Pic b compared to previous orbital fits. If confirmed by future observations, this eccentricity may help explain inner edge, scale height, and brightness asymmetry of $β$ Pic's disk. It could also potentially signal that $β$ Pic b has migrated inward to its current location, acquiring its eccentricity from interaction with the 3:1 outer Lindblad resonance in the disk.
△ Less
Submitted 30 December, 2018;
originally announced December 2018.