We argue that impact velocities between dust grains with sizes less than $\sim 0.1$ $\mu m$ in molecular cloud cores are dominated by drift arising from ambipolar diffusion. This effect is due to the size dependence of the dust coupling to the magnetic field and the neutral gas. Assuming perfect sticking in collisions up to $\approx 50$ m/s, we show that this effect causes rapid depletion of small grains - consistent with starlight extinction and IR/microwave emission measurements, both in the core center ($n \sim 10^{6}$ cm$^{-3}$) and envelope ($n \sim 10^{4}$ cm$^{-3}$). The upper end of the size distribution does not change significantly if only velocities arising from this effect are considered. We consider the impact of an evolved dust size distribution on the gas temperature, and argue that if the depletion of small dust grains occurs as would be expected from our model, then the cosmic ray ionization rate must be well below $10^{-16}$ s$^{-1}$ at a number density of $10^{5}$ cm$^{-3}$.

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.

We develop a self-consistent model for the equilibrium gas temperature and size-dependent dust temperature in cold, dense pre-stellar cores, assuming an arbitrary power-law size distribution of dust grains. Compact analytical expressions applicable to a broad range of physical parameters are derived and compared with predictions of the commonly used standard model. It is suggested that combining the theoretical results with observations should allow us to constrain the degree of dust evolution and the cosmic-ray ionization rate in dense cores, and to help in discriminating between different regimes of cosmic-ray transport in molecular clouds. In particular, assuming a canonical MRN distribution of grain sizes, our theory demonstrates that the gas temperature measurements in the pre-stellar core L1544 are consistent with an ionization rate as high as $\sim 10^{-16}$ s$^{-1}$, an order of magnitude higher than previously thought.

Chondrules are silicate spheroids found in meteorites, serving as important fossil records of the early solar system. In order to form chondrules, chondrule precursors must be heated to temperatures much higher than the typical conditions in the current asteroid belt. One proposed mechanism for chondrule heating is the passage through bow shocks of highly eccentric planetesimals in the protoplanetary disk in the early solar system. However, it is difficult for planetesimals to gain and maintain such high eccentricities. In this paper, we present a new scenario in which planetesimals in the asteroid belt region are excited to high eccentricities by the Jovian sweeping secular resonance in a depleting disk, leading to efficient formation of chondrules. We study the orbital evolution of planetesimals in the disk using semi-analytic models and numerical simulations. We investigate the dependence of eccentricity excitation on the planetesimal's size as well as the physical environment, and calculate the probability for chondrule formation. We find that 50 - 2000 km planetesimals can obtain eccentricities larger than 0.6 and cause effective chondrule heating. Most chondrules form in high velocity shocks, in low density gas, and in the inner disk. The fraction of chondrule precursors which become chondrules is about 4 - 9 % between 1.5 - 3 AU. Our model implies that the disk depletion timescale is $\tau_\mathrm{dep}\approx 1~\mathrm{Myr}$, comparable to the age spread of chondrules; and that Jupiter formed before chondrules, no more than 0.7 Myr after the formation of the CAIs.

Understanding the cosmic ray (CR) ionization rate is crucial in order to simulate the dynamics of, and interpret the chemical species observed in molecular clouds. Calculating the CR ionization rate requires both accurate knowledge of the spectrum of MeV to GeV protons at the edge of the cloud as well as a model for the propagation of CRs into molecular clouds. Some models for the propagation of CRs in molecular clouds assume the CRs to stream freely along magnetic field lines, while in others they propagate diffusively due to resonant scattering off of magnetic disturbances excited by MHD turbulence present in the medium. We discuss the conditions under which CR diffusion can operate in a molecular cloud, calculate the local CR spectrum and ionization rate in both a free-streaming and diffusive propagation model, and highlight the different results from the two models. We also apply these two models to the propagation through the ISM to obtain the spectrum seen by Voyager 1, and show that such a spectrum favors a diffusive propagation model.

We study the combined impact of magnetic mirroring and focusing on the ionization by cosmic rays (CRs) in dense molecular clouds and circumstellar disks. We show that for effective column densities of up to $\sim10^{25}$ cm$^{-2}$ (where ionization is the main mechanism of energy losses by CRs) the two effects practically cancel each other out, provided the magnetic field strength has a single peak along field lines. In this case the ionization rate at a given location is controlled solely by attenuation of interstellar CRs due to energy losses. The situation is very different in the presence of magnetic pockets -- local minima of the field strength, where the CR density and thus ionization can be reduced drastically. We obtain simple analytical expressions allowing accurate calculation of the ionization rate in these regions.

It is likely that multiple bodies with masses between those of Mars and Earth ("planetary embryos") formed in the outer planetesimal disk of the solar system. Some of these were likely scattered by the giant planets into orbits with semi-major axes of hundreds of AU. Mutual torques between these embryos may lift the perihelia of some of them beyond the orbit of Neptune, where they are no longer perturbed by the giant planets so their semi-major axes are frozen in place. We conduct N-body simulations of this process, and its effect on smaller planetesimals in the region of the giant planets and the Kuiper belt. We find that (i) there is a significant possibility that one sub-Earth mass embryo, or possibly more, is still present in the outer solar system; (ii) the orbit of the surviving embryo(s) typically has perihelion of 40--70 AU, semi-major axis less than 200 AU, and inclination less than 30 degrees; (iii) it is likely that any surviving embryos could be detected by current or planned optical surveys or have a significant effect on solar-system ephemerides; (iv) whether or not an embryo has survived to the present day, their dynamical influence earlier in the history of the solar system can explain the properties of the detached disk (defined in this paper as containing objects with perihelia > 38 AU and semi-major axes between 80 and 500 AU).

We show that a black-hole binary with an external companion can undergo Lidov-Kozai cycles that cause a close pericenter passage, leading to a rapid merger due to gravitational-wave emission. This scenario occurs most often for systems in which the companion has mass comparable to the reduced mass of the binary and the companion orbit has semi-major axis within a factor of $\sim 10$ of the binary semi-major axis. Using a simple population-synthesis model and 3-body simulations, we estimate the rate of mergers in triple black hole systems in the field to be about six per Gpc$^3$ per year in the absence of natal kicks during black hole formation. This value is within the low end of the 90\% credible interval for the total black-hole black-hole merger rate inferred from the current LIGO results. There are many uncertainties in these calculations, the largest of which is the unknown distribution of natal kicks. Even modest natal kicks of $40\mbox{km s}^{-1}$ will reduce the merger rate by a factor of 40. A few percent of these systems will have eccentricity greater than 0.999 when they first enter the frequency band detectable by aLIGO (above 10 Hz).

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.

We investigate the claim that the ratio {\beta} 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 {\beta} is almost independent of grain porosity, but for small (<0.3 {\mu}m) grains, {\beta} 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, {\beta} 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 {\beta} 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.

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.

Detections of planets in eccentric, close (separations of ~20 AU) binary systems such as \alpha Cen or \gamma 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.

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 $\alpha$ Cen and $\gamma$ 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.

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.

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.