Pak Shing Li

Simulating Radiating and Magnetized Flows in Multiple Dimensions with ZEUS-MP

This paper describes ZEUS-MP, a multi-physics, massively parallel, message-passing implementation of the ZEUS code. ZEUS-MP differs significantly from the thoroughly documented ZEUS-2D code, the completely undocumented (in peer-reviewed literature) ZEUS-3D code, and a marginally documented “version 1” of ZEUS-MP first distributed publicly in 1999. ZEUS-MP offers an MHD algorithm which is better suited for multidimensional flows than the ZEUS-2D module by virtue of modifications to the Method of Characteristics scheme first suggested by Hawley & Stone (1995). This MHD module is shown to compare quite favorably to the TVD scheme described by Ryu et al. (1998). ZEUS-MP is the first publicly-available ZEUS code to allow the advection of multiple chemical (or nuclear) species. Radiation hydrodynamic simulations are enabled via an implicit flux-limited radiation diffusion (FLD) module. The hydrodynamic, MHD, and FLD modules may be used, singly or in concert, in one, two, or three space dimensions. Additionally, so-called “1.5-D” and “2.5-D” grids, in which the “half-D” denotes a symmetry axis along which a constant but non-zero value of velocity or magnetic field is evolved, are supported. Self gravity may be included either through the assumption of a GM/r potential or a solution of Poisson’s equation using one of three linear solver packages (conjugategradient, multigrid, and FFT) provided for that purpose. Point-mass potentials are also supported. Because ZEUS-MP is designed for large simulations on parallel computing platforms, considerable attention is paid to the parallel performance characteristics of each module in the code. Strong-scaling tests involving pure hydrodynamics (with and without self-gravity), MHD, and RHD are performed in which large problems (256 3 zones) are distributed among as many as 1024 processors of an IBM SP3. Parallel efficiency is a strong function of the amount of communication required between processors in a given algorithm, but all modules are shown to scale well on up to 1024 processors for the chosen fixed problem size. Subject headings: hydrodynamics – methods:numerical – methods:parallel – MHD – radiative transfer

ASCA discovery of an x-ray pulsar in the error box of SGR 1900+14

We present a 2-10 keV ASCA observation of the field around the soft gamma repeater SGR 1900+14. One quiescent X-ray source was detected in this observation, and it was in the SGR error box. In 2-10 keV X-rays, its spectrum may be fitted by a power law with index -2.2, and its unabsorbed flux is 9.6 × 10-12 ergs cm-2 s-1. We also find a clear 5.16 s period. The properties of the three well-studied soft gamma repeaters are remarkably similar to one another, and they provide evidence that all of them are associated with young, strongly magnetized neutron stars in supernova remnants.

Two regimes of turbulent fragmentation and the stellar initial mass function from primordial to present-day star formation

The Padoan and Nordlund model of the stellar initial mass function (IMF) is derived from low-order statistics of supersonic turbulence, neglecting gravity (e.g., gravitational fragmentation, accretion, and merging). In this work, the predictions of that model are tested using the largest numerical experiments of supersonic hydrodynamic (HD) and magnetohydrodynamic (MHD) turbulence to date (~10003 computational zones) and three different codes (Enzo, Zeus, and the Stagger code). The model predicts a power-law distribution for large masses, related to the turbulence-energy power-spectrum slope and the shock-jump conditions. This power-law mass distribution is confirmed by the numerical experiments. The model also predicts a sharp difference between the HD and MHD regimes, which is recovered in the experiments as well, implying that the magnetic field, even below energy equipartition on the large scale, is a crucial component of the process of turbulent fragmentation. These results suggest that the stellar IMF of primordial stars may differ from that in later epochs of star formation, due to differences in both gas temperature and magnetic field strength. In particular, we find that the IMF of primordial stars born in turbulent clouds may be narrowly peaked around a mass of order 10 M☉, as long as the column density of such clouds is not much in excess of 1022 cm-2.

The Formation of Self-Gravitating Cores in Turbulent Magnetized Clouds

We use ZEUS-MP to perform high-resolution, three-dimensional, super-Alfvenic turbulent simulations in order to investigate the role of magnetic fields in self-gravitating core formation within turbulent molecular clouds. Statistical properties of our super-Alfvenic model without gravity agree with previous similar studies. Including self-gravity, our models give the following results. They are consistent with the turbulent frag- mentation prediction of the core mass distribution of Padoan & Nordlund. They also confirm that local gravitational collapse is not prevented by magnetohydrodynamic waves driven by turbulent flows, even when the turbulent Jeans mass exceeds the mass in the simulation volume. Comparison of results between 256 3 and 512 3 zone simulations reveals convergence in the collapse rate. Analysis of self-gravitating cores formed in the simulation shows the following: (1) All cores formed are magnetically supercritical by at least an order of magnitude. (2) A power-law relation between central magnetic field strength and density Bc / � 1=2 c is observed despite the cores being strongly supercritical. (3) Specific angular momentum j / R 3=2 for cores with radius R. (4) Most cores are prolate and triaxial in shape, in agreement with the results of Gammie and coworkers. We find a weak correlation between the minor axis of the core and the local magnetic field in our simulation at late times, different from the uncorrelated results reported by Gammie and coworkers. The core shape analysis and the power-law relationship between core mass and radius M / R 2:75 suggest the formation of some highly flattened cores. We identified 12 cloud cores with disklike appearance at a later stage of our high-resolution simulation. Instead of being tidally truncated or disrupted, the core disks survive and flourish while undergoing strong interactions. We discuss the physical properties of these disklike cores under the constraints of resolution limits. Subject headings: ISM: clouds — ISM: kinematics and dynamics — ISM: magnetic fields — methods: numerical — stars: formation — turbulence On-line material: mpeg animation

COMPARING NUMERICAL METHODS FOR ISOTHERMAL MAGNETIZED SUPERSONIC TURBULENCE

Many astrophysical applications involve magnetized turbulent flows with shock waves. Ab initio star formation simulations require a robust representation of supersonic turbulence in molecular clouds on a wide range of scales imposing stringent demands on the quality of numerical algorithms. We employ simulations of supersonic super-Alfvenic turbulence decay as a benchmark test problem to assess and compare the performance of nine popular astrophysical MHD methods actively used to model star formation. The set of nine codes includes: ENZO, FLASH, KT-MHD, LL-MHD, PLUTO, PPML, RAMSES, STAGGER, and ZEUS. These applications employ a variety of numerical approaches, including both split and unsplit, finite difference and finite volume, divergence preserving and divergence cleaning, a variety of Riemann solvers, and a range of spatial reconstruction and time integration techniques. We present a comprehensive set of statistical measures designed to quantify the effects of numerical dissipation in these MHD solvers. We compare power spectra for basic fields to determine the effective spectral bandwidth of the methods and rank them based on their relative effective Reynolds numbers. We also compare numerical dissipation for solenoidal and dilatational velocity components to check for possible impacts of the numerics on small-scale density statistics. Finally, we discuss the convergence of various characteristics for the turbulence decay test and the impact of various components of numerical schemes on the accuracy of solutions. The nine codes gave qualitatively the same results, implying that they are all performing reasonably well and are useful for scientific applications. We show that the best performing codes employ a consistently high order of accuracy for spatial reconstruction of the evolved fields, transverse gradient interpolation, conservation law update step, and Lorentz force computation. The best results are achieved with divergence-free evolution of the magnetic field using the constrained transport method and using little to no explicit artificial viscosity. Codes that fall short in one or more of these areas are still useful, but they must compensate for higher numerical dissipation with higher numerical resolution. This paper is the largest, most comprehensive MHD code comparison on an application-like test problem to date. We hope this work will help developers improve their numerical algorithms while helping users to make informed choices about choosing optimal applications for their specific astrophysical problems.

Magnetized interstellar molecular clouds – I. Comparison between simulations and Zeeman observations

The most accurate measurements of magnetic fields in star-forming gas are based on the Zeeman observations analyzed by Crutcher et al. (2010). We show that their finding that the 3D magnetic field scales approximately as density

A STABLE, ACCURATE METHODOLOGY FOR HIGH MACH NUMBER, STRONG MAGNETIC FIELD MHD TURBULENCE WITH ADAPTIVE MESH REFINEMENT: RESOLUTION AND REFINEMENT STUDIES

^{0.65}

AMBIPOLAR DIFFUSION HEATING IN TURBULENT SYSTEMS

can also be obtained from analysis of the observed line-of-sight fields. We present two large-scale AMR MHD simulations of several thousand

Photoionization Rates in Clumpy Molecular Clouds

M_\odot

ASCA Observation of the Quiescent X-ray Counterpart to SGR 1627-41

of turbulent, isothermal, self-gravitating gas, one with a strong initial magnetic field (Alfven Mach number