BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Kaitlin M. Kratter

Kaitlin M. Kratter appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.SRastro-ph.EPastro-ph.GAastro-ph.HE

Catalog footprint

What is connected

35works
4topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

35 published item(s)

preprint2022arXiv

Multiplicity Statistics of Stars in the Sagittarius Dwarf Spheroidal Galaxy: Comparison to the Milky Way

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Orbital Architectures of Planet-Hosting Binaries II. Low Mutual Inclinations Between Planetary and Stellar Orbits

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Substructures in protoplanetary disks imprinted by compact planetary systems

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The 3D dust and opacity distribution of protoplanets in multi-fluid global simulations

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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The VLA/ALMA Nascent Disk and Multiplicity (VANDAM) Survey of Orion Protostars V. A Characterization of Protostellar Multiplicity

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

ARMADA I: Triple Companions Detected in B-Type Binaries alpha Del and nu Gem

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Dust Settling Instability in Protoplanetary Discs

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Dynamics of small grains in transitional discs

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Impact of Binary Stars on Planet Statistics -- I. Planet Occurrence Rates, Trends with Stellar Mass, and Wide Companions to Hot Jupiter Hosts

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Kinematic Analysis of a Protostellar Multiple System: Measuring the Protostar Masses and Assessing Gravitational Instability in the Disks of L1448 IRS3B and L1448 IRS3A

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The highly variable time evolution of star-forming cores identified with dendrograms

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The McDonald Accelerating Stars Survey (MASS): White Dwarf Companions Accelerating the Sun-like Stars 12 Psc and HD 159062

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Rich Lack Close Neighbours: The Dependence of Blue-Straggler Fraction on Metallicity

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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

The Formation and Evolution of Wide-Orbit Stellar Multiples In Magnetized Clouds

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Adaptive Optics imaging of VHS 1256-1257: A Low Mass Companion to a Brown Dwarf Binary System

Recently, Gauza et al. (2015) reported the discovery of a companion to the late M-dwarf, VHS J125601.92-125723.9 (VHS 1256-1257). The companion's absolute photometry suggests its mass and atmosphere are similar to the HR 8799 planets. However, as a wide companion to a late-type star, it is more accessible to spectroscopic characterization. We discovered that the primary of this system is an equal-magnitude binary. For an age $\sim300$ Myr the A and B components each have a mass of $64.6^{+0.8}_{-2.0}~M_{\mathrm{Jup}}$, and the b component has a mass of $11.2^{+9.7}_{-1.8}$, making VHS 1256-1257 only the third brown dwarf triple system. There exists some tension between the spectrophotometric distance of $17.2\pm2.6$ pc and the parallax distance of $12.7\pm1.0$ pc. At 12.7 pc VHS1256-1257 A and B would be the faintest known M7.5 objects, and are even faint outliers among M8 types. If the larger spectrophotmetric distance is more accurate than the parallax, then the mass of each component increases. In particular, the mass of the b component increases well above the deuterium burning limit to $\sim35~M_{\mathrm{Jup}}$ and the mass of each binary component increases to $73^{+20}_{-17}~M_{\mathrm{Jup}}$. At 17.1 pc, the UVW kinematics of the system are consistent with membership in the AB~Dor moving group. The architecture of the system resembles a hierarchical stellar multiple suggesting it formed via an extension of the star-formation process to low masses. Continued astrometric monitoring will resolve this distance uncertainty and will provide dynamical masses for a new benchmark system.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On the gravitational stability of gravito-turbulent accretion disks

Low mass, self-gravitating accretion disks admit quasi-steady,`gravito-turbulent' states in which cooling balances turbulent viscous heating. However, numerical simulations show that gravito-turbulence cannot be sustained beyond dynamical timescales when the cooling rate or corresponding turbulent viscosity is too large. The result is disk fragmentation. We motivate and quantify an interpretation of disk fragmentation as the inability to maintain gravito-turbulence due to formal secondary instabilities driven by: 1) cooling, which reduces pressure support; and/or 2) viscosity, which reduces rotational support. We analyze the axisymmetric gravitational stability of viscous, non-adiabatic accretion disks with internal heating, external irradiation, and cooling in the shearing box approximation. We consider parameterized cooling functions in 2D and 3D disks, as well as radiative diffusion in 3D. We show that generally there is no critical cooling rate/viscosity below which the disk is formally stable, although interesting limits appear for unstable modes with lengthscales on the order of the disk thickness. We apply this new linear theory to protoplanetary disks subject to gravito-turbulence modeled as an effective viscosity, and cooling regulated by dust opacity. We find that viscosity renders the disk beyond $\sim 60$AU dynamically unstable on radial lengthscales a few times the local disk thickness. This is coincident with the empirical condition for disk fragmentation based on a maximum sustainable stress. We suggest turbulent stresses can play an active role in realistic disk fragmentation by removing rotational stabilization against self-gravity, and that the observed transition in behavior from gravito-turbulent to fragmenting may reflect instability of the gravito-turbulent state itself.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Orbital Stability of Multi-Planet Systems: Behavior at High Masses

In the coming years, high contrast imaging surveys are expected to reveal the characteristics of the population of wide-orbit, massive, exoplanets. To date, a handful of wide planetary mass companions are known, but only one such multi-planet system has been discovered: HR8799. For low mass planetary systems, multi-planet interactions play an important role in setting system architecture. In this paper, we explore the stability of these high mass, multi-planet systems. While empirical relationships exist that predict how system stability scales with planet spacing at low masses, we show that extrapolating to super-Jupiter masses can lead to up to an order of magnitude overestimate of stability for massive, tightly packed systems. We show that at both low and high planet masses, overlapping mean motion resonances trigger chaotic orbital evolution, which leads to system instability. We attribute some of the difference in behavior as a function of mass to the increasing importance of second order resonances at high planet-star mass ratios. We use our tailored high mass planet results to estimate the maximum number of planets that might reside in double component debris disk systems, whose gaps may indicate the presence of massive bodies.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Planet Scattering Around Binaries: Ejections, Not Collisions

Transiting circumbinary planets discovered by Kepler provide unique insight into binary star and planet formation. Several features of this new found population, for example the apparent pile-up of planets near the innermost stable orbit, may distinguish between formation theories. In this work, we determine how planet-planet scattering shapes planetary systems around binaries as compared to single stars. In particular, we look for signatures that arise due to differences in dynamical evolution in binary systems. We carry out a parameter study of N-body scattering simulations for four distinct planet populations around both binary and single stars. While binarity has little influence on the final system multiplicity or orbital distribution, the presence of a binary dramatically effects the means by which planets are lost from the system. Most circumbinary planets are lost due to ejections rather than planet-planet or planet-star collisions. The most massive planet in the system tends to control the evolution. Systems similar to the only observed multi-planet circumbinary system, Kepler-47, can arise from much more tightly packed, unstable systems. Only extreme initial conditions introduce differences in the final planet populations. Thus, we suggest that any intrinsic differences in the populations are imprinted by formation.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Orbital Architectures of Planet-Hosting Binaries: I. Forming Five Small Planets in the Truncated Disk of Kepler-444A

We present the first results from our Keck program investigating the orbital architectures of planet-hosting multiple star systems. Kepler-444 is a metal-poor triple star system that hosts five sub-Earth-sized planets orbiting the primary star (Kepler-444A), as well as a spatially unresolved pair of M dwarfs (Kepler-444BC) at a projected distance of 1.8" (66 AU). We combine our Keck/NIRC2 adaptive optics astrometry with multi-epoch Keck/HIRES RVs of all three stars to determine a precise orbit for the BC pair around A, given their empirically constrained masses. We measure minimal astrometric motion ($1.0\pm0.6$ mas yr$^{-1}$, or $0.17\pm0.10$ km s$^{-1}$), but our RVs reveal significant orbital velocity ($1.7\pm0.2$ km s$^{-1}$) and acceleration ($7.8\pm0.5$ m s$^{-1}$ yr$^{-1}$). We determine a highly eccentric stellar orbit ($e=0.864\pm0.023$) that brings the tight M dwarf pair within $5.0^{+0.9}_{-1.0}$ AU of the planetary system. We validate that the system is dynamically stable in its present configuration via n-body simulations. We find that the A$-$BC orbit and planetary orbits are likely aligned (98%) given that they both have edge-on orbits and misalignment induces precession of the planets out of transit. We conclude that the stars were likely on their current orbits during the epoch of planet formation, truncating the protoplanetary disk at $\approx$2 AU. This truncated disk would have been severely depleted of solid material from which to form the total $\approx$1.5 $M_{\rm Earth}$ of planets. We thereby strongly constrain the efficiency of the conversion of dust into planets and suggest that the Kepler-444 system is consistent with models that explain the formation of more typical close-in Kepler planets in normal, not truncated, disks.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Planet-disc interaction on a freely moving mesh

General-purpose, moving-mesh schemes for hydrodynamics have opened the possibility of combining the accuracy of grid-based numerical methods with the flexibility and automatic resolution adaptivity of particle-based methods. Due to their supersonic nature, Keplerian accretion discs are in principle a very attractive system for applying such freely moving mesh techniques. However, the high degree of symmetry of simple accretion disc models can be difficult to capture accurately by these methods, due to the generation of geometric grid noise and associated numerical diffusion, which is absent in polar grids. To explore these and other issues, in this work we study the idealized problem of two-dimensional planet-disc interaction with the moving-mesh code AREPO. We explore the hydrodynamic evolution of discs with planets through a series of numerical experiments that vary the planet mass, the disc viscosity and the mesh resolution, and compare the resulting surface density, vortensity field and tidal torque with results from the literature. We find that the performance of the moving-mesh code in this problem is in accordance with published results, showing good consistency with grid codes written in polar coordinates. We also conclude that grid noise and mesh distortions do not introduce excessive numerical diffusion. Finally, we show how the moving-mesh approach can naturally increase resolution in regions of high densityaround planets and planetary wakes, while retaining the background flow at low resolution. This provides an alternative to the difficult task of implementing adaptive mesh refinement in conventional polar-coordinate codes.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Simulated Photoevaporative Mass Loss from Hot Jupiters in 3D

Ionizing stellar photons heat the upper regions of planetary atmospheres, driving atmospheric mass loss. Gas escaping from several hot, hydrogen-rich planets has been detected using UV and X-ray transmission spectroscopy. Because these planets are tidally locked, and thus asymmetrically irradiated, escaping gas is unlikely to be spherically symmetric. In this paper, we focus on the effects of asymmetric heating on local outflow structure. We use the Athena code for hydrodynamics to produce 3D simulations of hot Jupiter mass loss that jointly model wind launching and stellar heating via photoionization. Our fiducial planet is an inflated, hot Jupiter with radius $R_p=2.14 R_{\rm Jup}$ and mass $M_p = 0.53 M_{\rm Jup}$. We irradiate the initially neutral, atomic hydrogen atmosphere with 13.6 eV photons and compute the outflow's ionization structure. There are clear asymmetries in the atmospheric outflow, including a neutral shadow on the planet's nightside. Given an incident ionizing UV flux comparable to that of the Sun, we find a steady-state mass loss rate of ~$2\times10^{10}$ g s$^{-1}$. The total mass loss rate and the outflow substructure along the substellar ray show good agreement with earlier 1D models, for two different fluxes. Our 3D data cube can be used to generate the outflow's extinction spectrum during transit. As a proof of concept, we find absorption of stellar Lyman-alpha at Doppler-shifted velocities of up to $\pm 50$ km s$^{-1}$. Our work provides a starting point for further 3D models that can be used to predict observable signatures of hot Jupiter mass loss.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Stellar orbit evolution in close circumstellar disc encounters

The formation and early evolution of circumstellar discs often occurs within dense, newborn stellar clusters. For the first time, we apply the moving-mesh code AREPO, to circumstellar discs in 3-D, focusing on disc-disc interactions that result from stellar fly-bys. Although a small fraction of stars are expected to undergo close approaches, the outcomes of the most violent encounters might leave an imprint on the discs and host stars that will influence both their orbits and their ability to form planets. We first construct well-behaved 3-D models of self-gravitating discs, and then create a suite of numerical experiments of parabolic encounters, exploring the effects of pericenter separation r_p, disc orientation and disc-star mass ratio (M_d/M_*) on the orbital evolution of the host stars. Close encounters (2r_p<~ disc radius) can truncate discs on very short time scales. If discs are massive, close encounters facilitate enough orbital angular momentum extraction to induce stellar capture. We find that for realistic primordial disc masses M_d<~0.1M_*, non-colliding encounters induce minor orbital changes, which is consistent with analytic calculations of encounters in the linear regime. The same disc masses produce entirely different results for grazing/colliding encounters. In the latter case, rapidly cooling discs lose orbital energy by radiating away the energy excess of the shock-heated gas, thus causing capture of the host stars into a bound orbit. In rare cases, a tight binary with a circumbinary disc forms as a result of this encounter.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Effect of Magnetic Misalignment on Protobinary Evolution

The majority of solar-type stars reside in multiple systems, especially binaries. They form in dense cores of molecular clouds that are observed to be significantly magnetized. Our previous study shows that magnetic braking can tighten the binary separation during the protostellar mass accretion phase by removing the angular momentum of the accreting material. Recent numerical calculations of single star formation have shown that misalignment between the magnetic field and rotation axis may weaken both magnetic braking and the associated magnetically driven outflows. These two effects allow for disk formation even in strongly magnetized cores. Here we investigate the effects of magnetic field misalignment on the properties of protobinaries. Somewhat surprisingly, the misaligned magnetic field is more efficient at tightening the binary orbit compared to the aligned field. The main reason is that the misalignment weakens the magnetically-driven outflow, which allows more material to accrete onto the binary. Even though the specific angular momentum of this inner material is higher than in the aligned case, it is insufficient to compensate for the additional mass. A corollary of this result is that a weaker field is required to achieve the same degree of inward migration when the field is tilted relative to the rotation axis. Large field misalignment also helps to produce rotationally-supported circumbinary disks even for relatively strong magnetic fields, by weakening the magnetically-dominated structure close to the binary. Our result may provide an explanation for the circumbinary disks detected in recent SMA and ALMA observations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Planet Packing in Circumbinary Systems

The recent discovery of planets orbiting main sequence binaries will provide crucial constraints for theories of binary and planet formation. The formation pathway for these planets is complicated by uncertainties in the formation mechanism of the host stars. In this paper, we compare the dynamical states of single and binary star planetary systems. Specifically, we pose two questions: (1) What does it mean for a circumbinary system to be dynamically packed? (2) How many systems are required to differentiate between a population of packed or sparse planets? We determine when circumbinary systems become dynamically unstable as a function of the separation between the host-stars and the inner planet, and the first and second planets. We show that these represent unique stability constraints compared to single-star systems. We find that although the existing Kepler data is insufficient to distinguish between a population of packed or sparse circumbinary systems, a more thorough study of circumbinary TTVs combined with an order of magnitude increase in the number of systems may prove conclusive. Future space missions such as TESS provide the best opportunity for increasing the sample size.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Spin-orbit alignment in the very low mass binary regime: The L dwarf tight binary 2MASSW J0746425+200032AB

Studies of solar-type binaries have found coplanarity between the equatorial and orbital planes of systems with $<$40 AU separation. By comparison, the alignment of the equatorial and orbital axes in the substellar regime, and the associated implications for formation theory, are relatively poorly constrained. Here we present the discovery of the rotation period of 3.32 $\pm$ 0.15 hours from 2MASS J0746+20A - the primary component of a tight (2.7 AU) ultracool dwarf binary system (L0+L1.5). The newly discovered period, together with the established period via radio observations of the other component, and the well constrained orbital parameters and rotational velocity measurements, allow us to infer alignment of the equatorial planes of both components with the orbital plane of the system to within 10 degrees. This result suggests that solar-type binary formation mechanisms may extend down into the brown dwarf mass range, and we consider a number of formation theories that may be applicable in this case. This is the first such observational result in the very low mass binary regime. In addition, the detected period of 3.32 $\pm$ 0.15 hours implies that the reported radio period of 2.07 $\pm$ 0.002 hours is associated with the secondary star, not the primary, as was previously claimed. This in turn refutes the claimed radius of 0.78 $\pm$ 0.1 $R_{J}$ for 2MASS J0746+20A, which we demonstrate to be 0.99 $\pm$ 0.03 $R_{J}$.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

An Empirical Correction for Activity Effects on the Temperatures, Radii, and Estimated Masses of Low-Mass Stars and Brown Dwarfs

We present empirical relations for determining the amount by which the effective temperatures and radii -- and therefore the estimated masses -- of low-mass stars and brown dwarfs are altered due to chromospheric activity. We base our relations on a large set of low-mass stars in the field with Halpha activity measurements, and on a set of low-mass eclipsing binaries with X-ray activity measurements from which we indirectly infer the Halpha activity. Both samples yield consistent relations linking the amount by which an active object's temperature is suppressed, and its radius inflated, to the strength of its Halpha emission. These relations are found to approximately preserve bolometric luminosity. We apply these relations to the peculiar brown-dwarf eclipsing binary 2M0535-05, in which the active, higher-mass brown dwarf has a cooler temperature than its inactive, lower-mass companion. The relations correctly reproduce the observed temperatures and radii of 2M0535-05 after accounting for the Halpha emission; 2M0535-05 would be in precise agreement with theoretical isochrones were it inactive. The relations that we present are applicable to brown dwarfs and low-mass stars with masses below 0.8 Msun and for which the activity, as measured by the fractional Halpha luminosity, is in the range -4.6 < log LHa/Lbol < -3.3. We expect these relations to be most useful for correcting radius and mass estimates of low-mass stars and brown dwarfs over their active lifetimes (few Gyr) and when the ages or distances (and therefore luminosities) are unknown. We also discuss the implications of this work for improved determinations of young cluster initial mass functions.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Circumbinary Chaos: Using Pluto's Newest Moon to Constrain the Masses of Nix & Hydra

The Pluto system provides a unique local laboratory for the study of binaries with multiple low mass companions. In this paper, we study the orbital stability of P4, the most recently discovered moon in the Pluto system. This newfound companion orbits near the plane of the Pluto-Charon binary, roughly halfway between the two minor moons Nix and Hydra. We use a suite of few body integrations to constrain the masses of Nix and Hydra, and the orbital parameters of P4. For the system to remain stable over the age of the Solar System, the masses of Nix and Hydra likely do not exceed 5e16 kg and 9e16 kg, respectively. These upper limits assume a fixed mass ratio between Nix and Hydra at the value implied by their median optical brightness. Our study finds that stability is more sensitive to their total mass and that a downward revision of Charon's eccentricity (from our adopted value of 0.0035) is unlikely to significantly affect our conclusions. Our upper limits are an order of magnitude below existing astrometric limits on the masses of Nix and Hydra. For a density at least that of ice, the albedos of Nix and Hydra would exceed 0.3. This constraint implies they are icy, as predicted by giant impact models. Even with these low masses, P4 only remains stable if its eccentricity e < 0.02. The 5:1 commensurability with Charon is particularly unstable, Combining stability constraints with the observed mean motion places the preferred orbit for P4 just exterior to the 5:1 resonance. These predictions will be tested when the New Horizons satellite visits Pluto. Based on the results for the Pluto-Charon system, we expect that circumbinary, multi-planet systems will be more widely spaced than their singleton counterparts. Further, circumbinary exoplanets close to the three-body stability boundary, such as those found by Kepler, are less likely to have other companions nearby.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Correcting for Activity Effects on the Temperatures, Radii, and Estimated Masses of Low-Mass Stars and Brown Dwarfs

We present empirical relations for determining the amount by which the effective temperatures and radii---and therefore the estimated masses---of low-mass stars and brown dwarfs are altered due to chromospheric activity. Accurate estimates of stellar radii are especially important in the context of searches for transiting exoplanets, which rely upon the assumed stellar radius/density to infer the planet radius/density. Our relations are based on a large set of well studied low-mass stars in the field and on a set of benchmark low-mass eclipsing binaries. The relations link the amount by which an active object's temperature is suppressed, and its radius inflated, to the strength of its Halpha emission. These relations are found to approximately preserve bolometric luminosity. We apply these relations to the peculiar brown-dwarf eclipsing binary 2M0535-05, in which the active, higher-mass brown dwarf has a cooler temperature than its inactive, lower-mass companion. The relations correctly reproduce the observed temperatures and radii of 2M0535-05 after accounting for the Halpha emission; 2M0535-05 would be in precise agreement with theoretical isochrones were it inactive. The relations that we present are applicable to brown dwarfs and low-mass stars with masses below 0.8 Msun and for which the activity, as measured by Halpha, is in the range -4.6 < log Lha/Lbol < -3.3. We expect these relations to be most useful for correcting radius and mass estimates of low-mass stars and brown dwarfs over their active lifetimes (few Gyr). We also discuss the implications of this work for determinations of young cluster IMFs.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Star Hoppers: Planet Instability and Capture in Evolving Binary Systems

Many planets are observed in stellar binary systems, and their frequency may be comparable to that of planetary systems around single stars. Binary stellar evolution in such systems influences the dynamical evolution of the resident planets. Here we study the evolution of a single planet orbiting one star in an evolving binary system. We find that stellar evolution can trigger dynamical instabilities that drive planets into chaotic orbits. This instability leads to planet-star collisions, exchange of the planet between the binary stars ("star-hoppers"), and ejection of the planet from the system. The means by which planets can be recaptured is similar to the pull-down capture mechanism for irregular solar system satellites. Because planets often suffer close encounters with the primary on the asymptotic giant branch, captures during a collision with the stellar envelope are also possible. Such capture could populate the habitable zone around white dwarfs.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The triple evolution dynamical instability: Stellar collisions in the field and the formation of exotic binaries

Physical collisions and close approaches between stars play an important role in the formation of exotic stellar systems. Standard theories suggest that collisions are rare, occurring only via random encounters between stars in dense clusters. We present a different formation pathway, the triple evolution dynamical instability (TEDI), in which mass loss in an evolving triple star system causes orbital instability. The subsequent chaotic orbital evolution of the stars triggers close encounters, collisions, exchanges between the stellar components, and the dynamical formation of eccentric compact binaries (including Sirius like binaries). We demonstrate that the rate of stellar collisions due to the TEDI is approximately 10^{-4} yr^{-1} per Milky-Way Galaxy, which is nearly 30 times higher than the total collision rate due to random encounters in the Galactic globular clusters. Moreover, we find that the dominant type of stellar collisions is qualitatively different; most collisions involve asymptotic giant branch stars, rather than main sequence, or slightly evolved stars, which dominate collisions in globular clusters. The TEDI mechanism should lead us to revise our understanding of collisions and the formation of compact, eccentric binaries in the field.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Fragment Production and Survival in Irradiated Disks: A Comprehensive Cooling Criterion

Accretion disks that become gravitationally unstable can fragment into stellar or sub-stellar companions. The formation and survival of these fragments depends on the precarious balance between self-gravity, internal pressure, tidal shearing, and rotation. Disk fragmentation depends on two key factors (1) whether the disk can get to the fragmentation boundary of Q=1, and (2) whether fragments can survive for many orbital periods. Previous work suggests that to reach Q=1, and have fragments survive, a disk must cool on an orbital timescale. Here we show that disks heated primarily by external irradiation always satisfy the standard cooling time criterion. Thus even though irradiation heats disks, and makes them more stable in general, once they reach the fragmentation boundary, they fragment more easily. We derive a new cooling criterion that determines fragment survival, and calculate a pressure modified Hill radius, which sets the maximum size of pressure-supported objects in a Keplerian disk. We conclude that fragmentation in protostellar disks might occur at slightly smaller radii than previously thought, and recommend tests for future simulations that will better predict the outcome of fragmentation in real disks.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Spin down of protostars through gravitational torques

Young protostars embedded in circumstellar discs accrete from an angular momentum-rich mass reservoir. Without some braking mechanism, all stars should be spinning at or near break-up velocity. In this paper, we perform simulations of the self-gravitational collapse of an isothermal cloud using the ORION adaptive mesh refinement code and investigate the role that gravitational torques might play in the spin-down of the dense central object. While magnetic effects likely dominate for low mass stars, high mass and Population III stars might be less well magnetised. We find that gravitational torques alone prevent the central object from spinning up to more than half of its breakup velocity, because higher rotation rates lead to bar-like deformations that enable efficient angular momentum transfer to the surrounding medium. We also find that the long-term spin evolution of the central object is dictated by the properties of the surrounding disc. In particular, spiral modes with azimuthal wavenumber $m=2$ couple more effectively to its spin than the lopsided $m=1$ mode, which was found to inhibit spin evolution. We suggest that even in the absence of magnetic fields, gravitational torques may provide an upper limit on stellar spin, and that moderately massive circumstellar discs can cause long-term spin down.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Formation of Binaries

Binary stars produce an array of dramatic astrophysical phenomena. They allow us to probe stellar structure, nuclear physics, and gravitational wave physics. They also produce the powerful supernovae that allow us to measure the scale of the universe. Despite their importance and ubiquity, many questions remain unanswered as to how the star formation process produces the wide array of stellar multiples that we observe. A complete model for binary formation encompasses three main components. We must know: (1) the primordial population of systems, (2) the influence of the dynamical processes that reshape this distribution as stars form and natal clusters disperse, and (3) the role of binary stellar evolution. In this article I review the most prominent theories for binary formation: turbulent core fragmentation, disk fragmentation, and competitive accretion. I argue that turbulent core fragmentation at all masses, with disk fragmentation added in at the upper end of the mass spectrum, might explain the trend towards increasing multiplicity for higher mass stars. In addition, I provide a brief overview of the observational statistics and of some of the important processes that modify the primordial distribution of stellar orbits.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Runts of the Litter: Why planets formed through gravitational instability can only be failed binary stars

Recent direct imaging discoveries suggest a new class of massive, distant planets around A stars. These widely separated giants have been interpreted as signs of planet formation driven by gravitational instability, but the viability of this mechanism is not clear cut. In this paper, we first discuss the local requirements for fragmentation and the initial fragment mass scales. We then consider whether the fragment's subsequent growth can be terminated within the planetary mass regime. Finally, we place disks in the larger context of star formation and disk evolution models. We find that in order for gravitational instability to produce planets, disks must be atypically cold in order to reduce the initial fragment mass. In addition, fragmentation must occur during a narrow window of disk evolution, after infall has mostly ceased, but while the disk is still sufficiently massive to undergo gravitational instability. Under more typical conditions, disk-born objects will likely grow well above the deuterium burning planetary mass limit. We conclude that if planets are formed by gravitational instability, they must be the low mass tail of the distribution of disk-born companions. To validate this theory, on-going direct imaging surveys must find a greater abundance of brown dwarf and M-star companions to A-stars. Their absence would suggest planet formation by a different mechanism such as core accretion, which is consistent with the debris disks detected in these systems.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

On the Role of Disks in the Formation of Stellar Systems: A Numerical Parameter Study of Rapid Accretion

We study rapidly accreting, gravitationally unstable disks with a series of global, three dimensional, numerical experiments using the code ORION. In this paper we conduct a numerical parameter study focused on protostellar disks, and show that one can predict disk behavior and the multiplicity of the accreting star system as a function of two dimensionless parameters which compare the disk's accretion rate to its sound speed and orbital period. Although gravitational instabilities become strong, we find that fragmentation into binary or multiple systems occurs only when material falls in several times more rapidly than the canonical isothermal limit. The disk-to-star accretion rate is proportional to the infall rate, and governed by gravitational torques generated by low-m spiral modes. We also confirm the existence of a maximum stable disk mass: disks that exceed ~50% of the total system mass are subject to fragmentation and the subsequent formation of binary companions.

0 citations0 reviewsNo code linkedOpen access