Source author record

Kedron Silsbee

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

9works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

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

9 published item(s)

preprint2022arXiv

Icy molecule desorption in interstellar grain collisions

Observations of gaseous complex organic molecules (COMs) in cold starless and prestellar cloud cores require efficient desorption of the COMs and their parent species from icy mantles on interstellar grains. With a simple astrochemical model, we investigate if mechanical removal of ice fragments in oblique collisions between grains in two size bins (0.01 and 0.1 micron) can substantially affect COM abundances. Two grain collision velocities were considered - 10 and 50 meters per second, corresponding to realistic grain relative speeds arising from ambipolar diffusion and turbulence, respectively. From the smaller grains, the collisions are assumed to remove a spherical cap with height equal to 1/3 and 1 ice mantle thickness, respectively. We find that the turbulence-induced desorption can elevate the gas-phase abundances of COMs by several orders of magnitude, reproducing observed COM abundances within an order of magnitude. Importantly, the high gaseous COM abundances are attained for long time-scales of up to 1 Myr and for a rather low methanol ice abundance, common for starless cores. The simple model, considering only two grain size bins and several assumptions, demonstrates a concept that must be tested with a more sophisticated approach.

preprint2019arXiv

O(N^2) fragmentation algorithm

Collisional fragmentation is a ubiquitous phenomenon arising in a variety of astrophysical systems, from asteroid belts to debris and protoplanetary disks. Numerical studies of fragmentation typically rely on discretizing the size distribution of colliding objects into a large number N of bins in mass space, usually logarithmically spaced. A standard approach for redistributing the debris produced in collisions into the corresponding mass bins results in O(N^3) calculation, which leads to significant computational overhead when N is large. Here we formulate a more efficient explicit O(N^2) fragmentation algorithm, which works when the size spectrum of fragments produced in an individual collision has a self-similar shape with only a single characteristic mass scale (which can have arbitrary dependence on the energy and masses of colliding objects). Fragment size spectra used in existing fragmentation codes typically possess this property. We also show that our O(N^2) approach can be easily extended to work with non-self-similar fragment size distributions, for which we provide a worked example. This algorithm offers a substantial speedup of fragmentation calculations for large N> 100, even over the implicit methods, making it an attractive tool for studying collisionally evolving systems.

preprint2016arXiv

Modeling the Nearly Isotropic Comet Population in Anticipation of LSST Observations

We run simulations to determine the expected distribution of orbital elements of nearly isotropic comets (NICs) in the outer solar system, assuming that these comets originate in the Oort Cloud at thousands of AU and are perturbed into the planetary region by the Galactic tide. We show that the Large Synoptic Survey Telescope (LSST) should detect and characterize the orbits of hundreds to thousands of NICs with perihelion distance outside 5 AU. Observing NICs in the outer solar system is our only way of directly detecting comets from the inner Oort Cloud, as these comets are dynamically excluded from the inner solar system by the giant planets. Thus the distribution of orbital elements constrains the spatial distribution of comets in the Oort cloud and the environment in which the solar system formed. Additionally, comet orbits can be characterized more precisely when they are seen far from the Sun as they have not been affected by non-gravitational forces.

preprint2016arXiv

Radiation pressure on fluffy submicron-sized grains

We investigate the claim that the ratio β of radiation pressure force to gravitational force on a dust grain in our solar system can substantially exceed unity for some grain sizes, provided that grain porosity is high enough. For model grains consisting of random aggregates of silicate spherules, we find that the maximum value of β is almost independent of grain porosity, but for small (<0.3 μm) grains, β actually decreases with increasing porosity. We also investigate the effect of metallic iron and amorphous carbon inclusions in the dust grains and find that while these inclusions do increase the radiation pressure cross-section, β remains below unity for grains with 3 pg of silicate material. These results affect the interpretation of the grain trajectories estimated from the {\it Stardust} mission, which were modeled assuming β values exceeding one. We find that radiation pressure effects are not large enough for particles Orion and Hylabrook captured by {\it Stardust} to be of interstellar origin given their reported impact velocities. We also consider the effects of solar radiation on transverse velocities and grain spin, and show that radiation pressure introduces both transverse velocities and equatorial spin velocities of several hundred meters per second for incoming interstellar grains at 2 AU. These transverse velocities are not important for modeling trajectories, but such spin rates may result in centrifugal disruption of aggregates.

preprint2015arXiv

Birth Locations of the Kepler Circumbinary Planets

The Kepler mission has discovered about a dozen circumbinary planetary systems, all containing planets on ~ 1 AU orbits. We place bounds on the locations in the circumbinary protoplanetary disk, where these planets could have formed through collisional agglomeration starting from small (km-sized or less) planetesimals. We first present a model of secular planetesimal dynamics that accounts for the (1) perturbation due to the eccentric precessing binary, as well as the (2) gravity and (3) gas drag from a precessing eccentric disk. Their simultaneous action leads to rich dynamics, with (multiple) secular resonances emerging in the disk. We derive analytic results for size-dependent planetesimal eccentricity, and demonstrate the key role of the disk gravity for circumbinary dynamics. We then combine these results with a simple model for collisional outcomes and find that in systems like Kepler 16, planetesimal growth starting with 10-100 m planetesimals is possible outside a few AU. The exact location exterior to which this happens is sensitive to disk eccentricity, density and precession rate, as well as to the size of the first generation of planetesimals. Strong perturbations from the binary in the inner part of the disk, combined with a secular resonance at a few AU inhibit the growth of km-sized planetesimals within 2 - 4 AU of the binary. In situ planetesimal growth in the Kepler circumbinary systems is possible only starting from large (few-km-sized) bodies in a low-mass disk with surface density less than 500 g/cm^2 at 1 AU.

preprint2014arXiv

Planet formation in binaries: dynamics of planetesimals perturbed by the eccentric protoplanetary disk and the secondary

Detections of planets in eccentric, close (separations of ~20 AU) binary systems such as αCen or γCep provide an important test of planet formation theories. Gravitational perturbations from the companion are expected to excite high planetesimal eccentricities resulting in destruction, rather than growth, of objects with sizes of up to several hundred km in collisions of similar-size bodies. It was recently suggested that gravity of a massive axisymmetric gaseous disk in which planetesimals are embedded drives rapid precession of their orbits, suppressing eccentricity excitation. However, disks in binaries are themselves expected to be eccentric, leading to additional planetesimal excitation. Here we develop secular theory of eccentricity evolution for planetesimals perturbed by the gravity of an elliptical protoplanetary disk (neglecting gas drag) and the companion. For the first time we derive an expression for the disturbing function due to an eccentric disk, which can be used for a variety of other astrophysical problems. We obtain explicit analytical solutions for planetesimal eccentricity evolution and delineate four different regimes of dynamical excitation. We show that in systems with massive (>10^{-2}M_Sun) disks, planetesimal eccentricity is usually determined by the gravity of the eccentric disk alone, and is comparable to the disk eccentricity. As a result, the latter imposes a lower limit on collisional velocities of solids, making their growth problematic. This fragmentation barrier can be removed if the gaseous disk rapidly precesses or if its own self-gravity is efficient at lowering disk eccentricity.

preprint2014arXiv

Planet formation in stellar binaries I: planetesimal dynamics in massive protoplanetary disks

About $20\%$ of exoplanets discovered by radial velocity surveys reside in stellar binaries. To clarify their origin one has to understand the dynamics of planetesimals in protoplanetary disks within binaries. The standard description, accounting for only gas drag and gravity of the companion star has been challenged recently, as the gravity of the protoplanetary disk was shown to play a crucial role in planetesimal dynamics. An added complication is the tendency of protoplanetary disks in binaries to become eccentric, giving rise to additional excitation of planetesimal eccentricity. Here, for the first time, we analytically explore secular dynamics of planetesimals in binaries such as $α$ Cen and $γ$ Cep under the combined action of (1) gravity of the eccentric protoplanetary disk, (2) perturbations due to the (coplanar) eccentric companion, and (3) gas drag. We derive explicit solutions for the behavior of planetesimal eccentricity ${\bf e}_p$ in non-precessing disks (and in precessing disks in certain limits). We obtain the analytical form of the distribution of relative velocities of planetesimals, which is a key input for understanding their collisional evolution. Disk gravity strongly influences relative velocities and tends to push sizes of planetesimals colliding with comparable objects at the highest speed to small values, $\sim 1$ km. We also find that planetesimals in eccentric protoplanetary disks apsidally aligned with the binary orbit collide at lower relative velocities than in mis-aligned disks. Our results highlight a decisive role that disk gravity plays in planetesimal dynamics in binaries.

preprint2014arXiv

Planet formation in stellar binaries II: overcoming the fragmentation barrier in alpha Centauri and gamma Cephei-like systems

Planet formation in small-separation (~20 AU) eccentric binaries such as gamma Cephei or alpha Centauri is believed to be adversely affected by the presence of the stellar companion. Strong dynamical excitation of planetesimals by the eccentric companion can result in collisional destruction (rather than growth) of 1-100 km objects, giving rise to the "fragmentation barrier" for planet formation. We revise this issue using a novel description of secular dynamics of planetesimals in binaries, which accounts for the gravity of the eccentric, coplanar protoplanetary disk, as well as gas drag. By studying planetesimal collision outcomes we show, in contrast to many previous studies, that planetesimal growth and subsequent formation of planets (including gas giants) in AU-scale orbits within ~20 AU separation binaries may be possible, provided that the protoplanetary disks are massive (>10^{-2}M_\odot) and only weakly eccentric (disk eccentricity <0.01). These requirements are compatible with both the existence of massive (several M_J) planets in gamma Cep-like systems and the results of recent simulations of gaseous disks in eccentric binaries. Terrestrial and Neptune-like planets can also form in lower-mass disks at small (sub-AU) radii. We find that fragmentation barrier is less of a problem in eccentric disks which are apsidally aligned with the binary orbit. Alignment gives rise to special locations, where (1) relative planetesimal velocities are low and (2) the timescale of their drag-induced radial drift is long. This causes planetesimal pileup at such locations in the disk and promotes their growth.

preprint2010arXiv

Spinning dust emission: the effect of rotation around a non-principal axis

We investigate the rotational emission from dust grains that rotate around non- principal axes. We argue that in many phases of the interstellar medium, the smallest grains, which dominate spinning dust emission, are likely to have their nutation state (orientation of principal axes relative to the angular momentum vector) randomized during each thermal spike. We recompute the excitation and damping rates associated with rotational emission from the grain permanent dipole, grain-plasma interactions, infrared photon emission, and collisions. The resulting spinning dust spectra gener- ally show a shift toward higher emissivities and peak frequencies relative to previous calculations.