Relativistic Hydrodynamic Flows Using Spatial and Temporal Adaptive Structured Mesh Refinement

Astrophysical relativistic flow problems require high resolution three-dimensional numerical simulations. In this paper, we describe a new parallel three-dimensional code for simulations of special relativistic hydrodynamics (SRHD) using both spatially and temporally structured adaptive mesh refinement (AMR). We used the method of lines to discretize the SRHD equations spatially and a total variation diminishing (TVD) Runge-Kutta scheme for time integration. For spatial reconstruction, we have implemented piecewise linear method (PLM), piecewise parabolic method (PPM), third order convex essentially non-oscillatory (CENO) and third and fifth order weighted essentially non-oscillatory (WENO) schemes. Flux is computed using either direct flux reconstruction or approximate Riemann solvers including HLL, modified Marquina flux, local Lax-Friedrichs flux formulas and HLLC. The AMR part of the code is built on top of the cosmological Eulerian AMR code {\sl enzo}. We discuss the coupling of the AMR framework with the relativistic solvers. Via various test problems, we emphasize the importance of resolution studies in relativistic flow simulations because extremely high resolution is required especially when shear flows are present in the problem. We also present the results of two 3d simulations of astrophysical jets: AGN jets and GRB jets. Resolution study of those two cases further highlights the need of high resolutions to calculate accurately relativistic flow problems.Comment: 14 pages, 23 figures. A section on 3D GRB jet simulation added. Accepted by ApJ

A Call for the Structured Physicist Report

Introduction: The field of diagnostic radiology continues to struggle with the clinical adoption of the structured interpretive report, with many radiologists preferring a semistructured, free-text dictation style to a more rigid, highly structured approach that some professional leaders have promoted [1]. Although structured reporting compliance in the radiologist community has been difficult to achieve, diagnostic radiologists have been thinking about and discussing this important issue for many years; it is also a part of the ACR’s Imaging 3.0_ campaign [2]. In the breast imaging community, the well-established BI-RADS_ recommendations produce a very structured report, with a discussion of interpretive findings culminating in a numeric BI-RADS score ranging from 0 to 6 [3]. Unlike some interpretive radiology reports, which can be ambiguous in terms of the next course of action, the BI-RADS scale is not only a diagnostic scale but also prescriptive of what the necessary follow-up should be

Numerical Models of Binary Neutron Star System Mergers. I.: Numerical Methods and Equilibrium Data for Newtonian Models

The numerical modeling of binary neutron star mergers has become a subject of much interest in recent years. While a full and accurate model of this phenomenon would require the evolution of the equations of relativistic hydrodynamics along with the Einstein field equations, a qualitative study of the early stages on inspiral can be accomplished by either Newtonian or post-Newtonian models, which are more tractable. In this paper we offer a comparison of results from both rotating and non-rotating (inertial) frame Newtonian calculations. We find that the rotating frame calculations offer significantly improved accuracy as compared with the inertial frame models. Furthermore, we show that inertial frame models exhibit significant and erroneous angular momentum loss during the simulations that leads to an unphysical inspiral of the two neutron stars. We also examine the dependence of the models on initial conditions by considering initial configurations that consist of spherical neutron stars as well as stars that are in equilibrium and which are tidally distorted. We compare our models those of Rasio & Shapiro (1992,1994a) and New & Tohline (1997). Finally, we investigate the use of the isolated star approximation for the construction of initial data.Comment: 32 pages, 19 gif figures, manuscript with postscript figures available at http://www.astro.sunysb.edu/dswesty/docs/nspap1.p

A Co-moving Coordinate System for Relativistic Hydrodynamics

The equations of relativistic hydrodynamics are transformed so that steps forward in time preserves local simultaneity. In these variables, the space-time coordinates of neighboring points on the mesh are simultaneous according to co-moving observers. Aside from the time step varying as a function of the location on the mesh, the local velocity gradient and the local density then evolve according to non-relativistic equations of motion. Analytic solutions are found for two one-dimensional cases with constant speed of sound. One solution has a Gaussian density profile when mapped into the new coordinates. That solution is analyzed for the effects of longitudinal acceleration in relativistic heavy ion collisions at RHIC, especially in regards to two-particle correlation measurements of the longitudinal size

Tidal disruption of dark matter halos around proto-globular clusters

Tidal disruption of dark matter halos around proto-globular clusters in a halo of a small galaxy is studied in the context of the hierarchical clustering scenario by using semi-cosmological N-body/SPH simulations assuming the standard cold dark matter model ($\Omega_0 = 1$). Our analysis on formation and evolution of the galaxy and its substructures archives until $z = 2.0$. In such a high-redshift universe, the Einstein-de Sitter universe is still a good approximation for a recently favored $\Lambda$-dominated universe, and then our results does not depend on the choice of cosmology. In order to resolve small gravitationally-bound clumps around galaxies and consider radiative cooling below $T = 10^4 K$, we adopt a fine mass resolution (m_{\rm SPH} = 1.12 \times 10^3 \Msun). Because of the cooling, each clump immediately forms a `core-halo' structure which consists of a baryonic core and a dark matter halo. The tidal force from the host galaxy mainly strips the dark matter halo from clumps and, as a result, theses clumps get dominated by baryons. Once a clump is captured by the host halo, its mass drastically decreases each pericenter passage. At $z = 2$, more than half of the clumps become baryon dominated systems (baryon mass/total mass $> 0.5$). Our results support the tidal evolution scenario of the formation of globular clusters and baryon dominated dwarf galaxies in the context of the cold dark matter universe.Comment: 9page, 13 figures. Accepted for publication in ApJ. A high-resolution PDF of the paper can be obtained from http://th.nao.ac.jp/~takayuki/ApJ05

Cooling Radiation and the Lyman-alpha Luminosity of Forming Galaxies

We examine the cooling radiation from forming galaxies in hydrodynamic simulations of the LCDM model (cold dark matter with a cosmological constant), focusing on the Ly-alpha line luminosities of high-redshift systems. Primordial composition gas condenses within dark matter potential wells, forming objects with masses and sizes comparable to the luminous regions of observed galaxies. As expected, the energy radiated in this process is comparable to the gravitational binding energy of the baryons, and the total cooling luminosity of the galaxy population peaks at z ~= 2. However, in contrast to the classical picture of gas cooling from the \sim 10^6 K virial temperature of a typical dark matter halo, we find that most of the cooling radiation is emitted by gas with T < 20,000 K. As a consequence, roughly 50% of this cooling radiation emerges in the Ly-alpha line. While a galaxy's cooling luminosity is usually smaller than the ionizing continuum luminosity of its young stars, the two are comparable in the most massive systems, and the cooling radiation is produced at larger radii, where the Ly-alpha photons are less likely to be extinguished by dust. We suggest, in particular, that cooling radiation could explain the two large (\sim 100 kpc), luminous (L_{Ly-alpha} \sim 10^{44} erg s^{-1}) ``blobs'' of Ly-alpha emission found in Steidel et al.'s (1999) narrow band survey of a z = 3 proto-cluster. Our simulations predict objects of the observed luminosity at about the right space density, and radiative transfer effects can account for the observed sizes and line widths. We discuss observable tests of this hypothesis for the nature of the Ly-alpha blobs, and we present predictions for the contribution of cooling radiation to the Ly-alpha luminosity function of galaxies as a function of redshift.Comment: Submitted to ApJ. 28 pages including 9 PS figures. Version with color figures available at http://donald.astro.umass.edu/~fardal/papers/cooling/cooling.htm

COSMOS: A Hybrid N-Body/Hydrodynamics Code for Cosmological Problems

We describe a new hybrid N-body/hydrodynamical code based on the particle-mesh (PM) method and the piecewise-parabolic method (PPM) for use in solving problems related to the evolution of large-scale structure, galaxy clusters, and individual galaxies. The code, named COSMOS, possesses several new features which distinguish it from other PM-PPM codes. In particular, to solve the Poisson equation we have written a new multigrid solver which can determine the gravitational potential of isolated matter distributions and which properly takes into account the finite-volume discretization required by PPM. All components of the code are constructed to work with a nonuniform mesh, preserving second-order spatial differences. The PPM code uses vacuum boundary conditions for isolated problems, preventing inflows when appropriate. The PM code uses a second-order variable-timestep time integration scheme. Radiative cooling and cosmological expansion terms are included. COSMOS has been implemented for parallel computers using the Parallel Virtual Machine (PVM) library, and it features a modular design which simplifies the addition of new physics and the configuration of the code for different types of problems. We discuss the equations solved by COSMOS and describe the algorithms used, with emphasis on these features. We also discuss the results of tests we have performed to establish that COSMOS works and to determine its range of validity.Comment: 43 pages, 14 figures, submitted to ApJS and revised according to referee's comment

Tidal spin-up of stars in dense stellar cusps around massive black holes

We show that main-sequence stars in dense stellar cusps around massive black holes are likely to rotate at a significant fraction of the centrifugal breakup velocity due to spin-up by hyperbolic tidal encounters. We use realistic stellar structure models to calculate analytically the tidal spin-up in soft encounters, and extrapolate these results to close and penetrating collisions using smoothed particle hydrodynamics simulations. We find that the spin-up falls off only slowly with distance from the black hole because the increased tidal coupling in slower collisions at larger distances compensates for the decrease in the stellar density. We apply our results to the stars near the massive black hole in the Galactic Center. Over their lifetime, ~1 Msol main sequence stars in the inner 0.3 pc of the Galactic Center are spun-up on average to ~10%--30% of the centrifugal breakup limit. Such rotation is ~20--60 times higher than is usual for such stars and may affect their subsequent evolution and their observed properties.Comment: 25 pages, 7 figures. Submitted to Ap

Effects of a Supermassive Black Hole Binary on a Nuclear Gas Disk

We study influence of a galactic central supermassive black hole (SMBH) binary on gas dynamics and star formation activity in a nuclear gas disk by making three-dimensional Tree+SPH simulations. Due to orbital motions of SMBHs, there are various resonances between gas motion and the SMBH binary motion. We have shown that these resonances create some characteristic structures of gas in the nuclear gas disk, for examples, gas elongated or filament structures, formation of gaseous spiral arms, and small gas disks around SMBHs. In these gaseous dense regions, active star formations are induced. As the result, many star burst regions are formed in the nuclear region.Comment: 19 pages, 11 figures, accepted for publication in Ap

Cosmological Evolution of Supergiant Star-Forming Clouds

In an exploration of the birthplaces of globular clusters, we present a careful examination of the formation of self-gravitating gas clouds within assembling dark matter haloes in a hierarchical cosmological model. Our high-resolution smoothed particle hydrodynamical simulations are designed to determine whether or not hypothesized supergiant molecular clouds (SGMCs) form and, if they do, to determine their physical properties and mass spectra. It was suggested in earlier work that clouds with a median mass of several 10^8 M_sun are expected to assemble during the formation of a galaxy, and that globular clusters form within these SGMCs. Our simulations show that clouds with the predicted properties are indeed produced as smaller clouds collide and agglomerate within the merging dark matter haloes of our cosmological model. We find that the mass spectrum of these clouds obeys the same power-law form observed for globular clusters, molecular clouds, and their internal clumps in galaxies, and predicted for the supergiant clouds in which globular clusters may form. We follow the evolution and physical properties of gas clouds within small dark matter haloes up to z = 1, after which prolific star formation is expected to occur. Finally, we discuss how our results may lead to more physically motivated "rules" for star formation in cosmological simulations of galaxy formation.Comment: Accepted to The Astrophysical Journal; 17 pages, 8 figure