Abstract
We present the results of a large suite of three-dimensional models of the collapse of magnetic molecular cloud cores using the adaptive mesh refinement code Enzo2.2 in the ideal magnetohydrodynamics approximation. The cloud cores are initially either prolate or oblate, centrally condensed clouds with masses of 1.73 or 2.73 Mo, respectively. The radial density profiles are Gaussian, with central densities 20 times higher than boundary densities. A barotropic equation of state is used to represent the transition from low density isothermal phases, to high density optically thick phases. The initial magnetic field strength ranges from 6.3 to 100 p,G, corresponding to clouds that are strongly to marginally supercritical, respectively, in terms of the mass to magnetic flux ratio. The magnetic field is initially uniform and aligned with the clouds' rotation axes, with initial ratios of rotational to gravitational energy ranging from 10-4 to 0.1. Two significantly different outcomes for collapse result: (1) formation of single protostars with spiral arms, and (2) fragmentation into multiple protostar systems. The transition between these two outcomes depends primarily on the initial magnetic field strength, with fragmentation occurring for mass to flux ratios greater than about 14 times the critical ratio for prolate clouds. Oblate clouds typically fragment into several times more clumps than prolate clouds. Multiple, rather than binary, system formation is the general rule in either case, suggesting that binary stars are primarily the result of the orbital dissolution of multiple protostar systems.