The initial stages of planet formation in circumstellar gas discs proceed via dust grains that collide and build up larger and larger bodies 1 . How this process continues from metre-sized boulders to kilometre-scale planetesimals is a major unsolved problem 2 : boulders stick together poorly 3 , and spiral into the protostar in a few hundred orbits due to a head wind from the slower rotating gas 4 . Gravitational collapse of the solid component has been suggested to overcome this barrier 1,5,6 . Even low levels of turbulence, however, inhibit sedimentation of solids to a sufficiently dense midplane layer 2, 7 , but turbulence must be present to explain observed gas accretion in protostellar discs 8 . Here we report the discovery of efficient gravitational collapse of boulders in locally overdense regions in the midplane. The boulders concentrate initially in transient high pressures in the turbulent gas 9 , and these concentra-1 arXiv:0708.3890v1 [astro-ph] 29 Aug 2007 tions are augmented a further order of magnitude by a streaming instability [10][11][12] driven by the relative flow of gas and solids. We find that gravitationally bound clusters form with masses comparable to dwarf planets and containing a distribution of boulder sizes. Gravitational collapse happens much faster than radial drift, offering a possible path to planetesimal formation in accreting circumstellar discs.Planet formation models typically treat turbulence as a diffusive process that opposes the gravitational sedimentation of solids to a high density midplane layer in circumstellar discs 7,13 . Recent models of solids moving in turbulent gas reveal that the turbulent motions not only mix them, but also concentrate metre-sized boulders in the transient gas overdensities 9 formed in magnetorotational turbulence 14 , in giant gaseous vortices 15,16 , and in spiral arms of self-gravitating discs 17 .Short-lived eddies at the dissipation scale of forced turbulence concentrate smaller millimetre-sized solids 18 . Some simulations mentioned above 9,11,12 were performed with the Pencil Code, which solves the magnetohydrodynamic (MHD) equations on a three-dimensional grid for a gas that interacts through drag forces with boulders. Boulders are represented as superparticles with independent positions and velocities, each having the mass of a huge number of boulders but the aerodynamic behaviour of a single boulder. We have now further developed the Pencil Code to include a fully parallel solver for the gravitational potential of the particles (see Supplementary Information). The particle density is mapped on the grid using the Triangular Shaped Cloud assignment scheme 19 and the gravitational potential of the solids is found using a Fast Fourier Transform method 20 . 2This allows us, for the first time, to simulate the dynamics of self-gravitating solid particles in magnetised, three-dimensional turbulence.We model a corotating, local box with linearised Keplerian shear that straddles the protoplanetary disc midplane and orbits the young star ...

The analysis of complex multiphysics astrophysical simulations presents a unique and rapidly growing set of challenges: reproducibility, parallelization, and vast increases in data size and complexity chief among them. In order to meet these challenges, and in order to open up new avenues for collaboration between users of multiple simulation platforms, we present yt a , an open source, communitydeveloped astrophysical analysis and visualization toolkit. Analysis and visualization with yt are oriented around physically relevant quantities rather than quantities native to astrophysical simulation codes. While originally designed for handling Enzo's structure adaptive mesh refinement (AMR) data, yt has been extended to work with several different simulation methods and simulation codes including Orion, RAMSES, and FLASH. We report on its methods for reading, handling, and visualizing data, including projections, multivariate volume rendering, multi-dimensional histograms, halo finding, light cone generation and topologically-connected isocontour identification. Furthermore, we discuss the underlying algorithms yt uses for processing and visualizing data, and its mechanisms for parallelization of analysis tasks.

This paper describes the open-source code Enzo, which uses block-structured adaptive mesh refinement to provide high spatial and temporal resolution for modeling astrophysical fluid flows. The code is Cartesian, can be run in 1, 2, and 3 dimensions, and supports a wide variety of physics including hydrodynamics, ideal and non-ideal magnetohydrodynamics, N-body dynamics (and, more broadly, self-gravity of fluids and particles), primordial gas chemistry, optically-thin radiative cooling of primordial and metal-enriched plasmas (as well as some optically-thick cooling models), radiation transport, cosmological expansion, and models for star formation and feedback in a cosmological context. In addition to explaining the algorithms implemented, we present solutions for a wide range of test problems, demonstrate the code's parallel performance, and discuss the Enzo collaboration's code development methodology.

We use three-dimensional hydrodynamical simulations to study the rapid infall phase of the common envelope interaction of a red giant branch star of mass equal to 0.88 M and a companion star of mass ranging from 0.9 down to 0.1 M . We first compare the results obtained using two different numerical techniques with different resolutions, and find overall very good agreement. We then compare the outcomes of those simulations with observed systems thought to have gone through a common envelope. The simulations fail to reproduce those systems in the sense that most of the envelope of the donor remains bound at the end of the simulations and the final orbital separations between the donor's remnant and the companion, ranging from 26.8 down to 5.9 R , are larger than the ones observed. We suggest that this discrepancy vouches for recombination playing an essential role in the ejection of the envelope and/or significant shrinkage of the orbit happening in the subsequent phase.

We study the buildup of magnetic fields during the formation of Population III star-forming regions, by conducting cosmological simulations from realistic initial conditions and varying the Jeans resolution. To investigate this in detail, we start simulations from identical initial conditions, mandating 16, 32 and 64 zones per Jeans length, and studied the variation in their magnetic field amplification. We find that, while compression results in some amplification, turbulent velocity fluctuations driven by the collapse can further amplify an initially weak seed field via dynamo action, provided there is sufficient numerical resolution to capture vortical motions (we find this requirement to be 64 zones per Jeans length, slightly larger than, but consistent with previous work run with more idealized collapse scenarios). We explore saturation of amplification of the magnetic field, which could potentially become dynamically important in subsequent, fully-resolved calculations. We have also identified a relatively surprising phenomena that is purely hydrodynamic: the higher-resolved simulations possess substantially different characteristics, including higher infall-velocity, increased temperatures inside 1000 AU, and decreased molecular hydrogen content in the innermost region. Furthermore, we find that disk formation is suppressed in higher-resolution calculations, at least at the times that we can follow the calculation. We discuss the effect this may have on the buildup of disks over the accretion history of the first clump to form as well as the potential for gravitational instabilities to develop and induce fragmentation.

Numerical solutions of partial differential equations enable a broad range of scientific research. The Dedalus Project is a flexible, open-source, parallelized computational framework for solving general partial differential equations using spectral methods. Dedalus translates plain-text strings describing partial differential equations into efficient solvers. This paper details the numerical method that enables this translation, describes the design and implementation of the codebase, and illustrates its capabilities with a variety of example problems. The numerical method is a first-order generalized tau formulation that discretizes equations into banded matrices. This method is implemented with an object-oriented design. Classes for spectral bases and domains manage the discretization and automatic parallel distribution of variables. Discretized fields and mathematical operators are symbolically manipulated with a basic computer algebra system. Initial value, boundary value, and eigenvalue problems are efficiently solved using high-performance linear algebra, transform, and parallel communication libraries. Custom analysis outputs can also be specified in plain text and stored in self-describing portable formats. The performance of the code is evaluated with a parallel scaling benchmark and a comparison to a finite-volume code. The features and flexibility of the codebase are illustrated by solving several examples: the nonlinear Schrodinger equation on a graph, a supersonic magnetohydrodynamic vortex, quasigeostrophic flow, Stokes flow in a cylindrical annulus, normal modes of a radiative atmosphere, and diamagnetic levitation. The Dedalus code and the example problems are available online at http://dedalus-project.org/. CONTENTS

The observed shock wave positions and expansion in Cas A can be interpreted in a model of supernova interaction with a freely expanding stellar wind with a mass loss rate of ∼ 3 × 10 −5 M ⊙ yr −1 for a wind velocity of 10 km s −1 . The wind was probably still being lost at the time of the supernova, which may have been of Type IIn or IIb. The wind may play a role in the formation of very fast knots observed in Cas A. In this model, the quasi-stationary flocculi (QSFs) represent clumps in the wind, with a density contrast of several 10 3 compared to the smooth wind. The outer, unshocked clumpy wind is photoionized by radiation from the supernova, and is observed as a patchy HII region around Cas A. This gas has a lower density than the QSFs and is heated by nonradiative shocks driven by the blast wave. Denser clumps have recombined and are observed as HI compact absorption features towards Cas A.

The nonlinear evolution of the Kelvin-Helmholtz instability is a popular test for code verification. To date, most Kelvin-Helmholtz problems discussed in the literature are ill-posed: they do not converge to any single solution with increasing resolution. This precludes comparisons among different codes and severely limits the utility of the Kelvin-Helmholtz instability as a test problem. The lack of a reference solution has led various authors to assert the accuracy of their simulations based on ad-hoc proxies, e.g., the existence of small-scale structures. This paper proposes well-posed Kelvin-Helmholtz problems with smooth initial conditions and explicit diffusion. We show that in many cases numerical errors/noise can seed spurious small-scale structure in Kelvin-Helmholtz problems. We demonstrate convergence to a reference solution using both Athena, a Godunov code, and Dedalus, a pseudo-spectral code. Problems with constant initial density throughout the domain are relatively straightforward for both codes. However, problems with an initial density jump (which are the norm in astrophysical systems) exhibit rich behavior and are more computationally challenging. In the latter case, Athena simulations are prone to an instability of the inner rolled-up vortex; this instability is seeded by grid-scale errors introduced by the algorithm, and disappears as resolution increases. Both Athena and Dedalus exhibit late-time chaos. Inviscid simulations are riddled with extremely vigorous secondary instabilities which induce more mixing than simulations with explicit diffusion. Our results highlight the importance of running well-posed test problems with demonstrated convergence to a reference solution. To facilitate future comparisons, we include the resolved, converged solutions to the Kelvin-Helmholtz problems in this paper in machine-readable form.