Astrophysical Radiation Hydrodynamics: The Prospects for Scaling

The general principles of scaling are discussed, followed by a survey of the important dimensionless parameters of fluid dynamics including radiation and magnetic fields, and of non-LTE spectroscopy. The values of the parameters are reviewed for a variety of astronomical and laboratory environments. It is found that parameters involving transport coefficients--the fluid and magnetic Reynolds numbers--have enormous values for the astronomical problems that are not reached in the lab. The parameters that measure the importance of radiation are also scarcely reached in the lab. This also means that the lab environments are much closer to LTE than the majority of astronomical examples. Some of the astronomical environments are more magnetically dominated than anything in the lab. The conclusion is that a good astronomical environment for simulation in a given lab experiment can be found, but that the reverse is much more difficult

Radiation Hydrodynamics

The discipline of radiation hydrodynamics is the branch of hydrodynamics in which the moving fluid absorbs and emits electromagnetic radiation, and in so doing modifies its dynamical behavior. That is, the net gain or loss of energy by parcels of the fluid material through absorption or emission of radiation are sufficient to change the pressure of the material, and therefore change its motion; alternatively, the net momentum exchange between radiation and matter may alter the motion of the matter directly. Ignoring the radiation contributions to energy and momentum will give a wrong prediction of the hydrodynamic motion when the correct description is radiation hydrodynamics. Of course, there are circumstances when a large quantity of radiation is present, yet can be ignored without causing the model to be in error. This happens when radiation from an exterior source streams through the problem, but the latter is so transparent that the energy and momentum coupling is negligible. Everything we say about radiation hydrodynamics applies equally well to neutrinos and photons (apart from the Einstein relations, specific to bosons), but in almost every area of astrophysics neutrino hydrodynamics is ignored, simply because the systems are exceedingly transparent to neutrinos, even though the energy flux in neutrinos may be substantial. Another place where we can do ''radiation hydrodynamics'' without using any sophisticated theory is deep within stars or other bodies, where the material is so opaque to the radiation that the mean free path of photons is entirely negligible compared with the size of the system, the distance over which any fluid quantity varies, and so on. In this case we can suppose that the radiation is in equilibrium with the matter locally, and its energy, pressure and momentum can be lumped in with those of the rest of the fluid. That is, it is no more necessary to distinguish photons from atoms, nuclei and electrons, than it is to distinguish hydrogen atoms from helium atoms, for instance. There are all just components of a mixed fluid in this case. So why do we have a special subject called ''radiation hydrodynamics'', when photons are just one of the many kinds of particles that comprise our fluid? The reason is that photons couple rather weakly to the atoms, ions and electrons, much more weakly than those particles couple with each other. Nor is the matter-radiation coupling negligible in many problems, since the star or nebula may be millions of mean free paths in extent. Radiation hydrodynamics exists as a discipline to treat those problems for which the energy and momentum coupling terms between matter and radiation are important, and for which, since the photon mean free path is neither extremely large nor extremely small compared with the size of the system, the radiation field is not very easy to calculate. In the theoretical development of this subject, many of the relations are presented in a form that is described as approximate, and perhaps accurate only to order of {nu}/c. This makes the discussion cumbersome. Why are we required to do this? It is because we are using Newtonian mechanics to treat our fluid, yet its photon component is intrinsically relativistic; the particles travel at the speed of light. There is a perfectly consistent relativistic kinetic theory, and a corresponding relativistic theory of fluid mechanics, which is perfectly suited to describing the photon gas. But it is cumbersome to use this for the fluid in general, and we prefer to avoid it for cases in which the flow velocity satisfies {nu} << c. The price we pay is to spend extra effort making sure that the source-sink terms relating to our relativistic gas component are included in the equations of motion in a form that preserves overall conservation of energy and momentum, something that would be automatic if the relativistic equations were used throughout

Aspherical Explosion Models for SN 1998bw/GRB 980425

The recent discovery of the unusual supernova SN1998bw and its apparent correlation with the gamma-ray burst GRB 980425 has raised new issues concerning both the GRB and supernovae. Although the spectra resemble those of TypeIc supernovae, there are distinct differences at early times and SN1998bw appeared to be unusually bright and red at maximum light. The apparent expansion velocities inferred by the Doppler shift of (unidentified) absorption features appeared to be high, making SN1998bw a possible candidate for a "hypernova" with explosion energies between 20 and 50E51 erg and ejecta masses in excess of 6 - 15 M_o. Based on light curve calculations for aspherical explosions and guided by the polarization observations of "normal" SNIc and related events, we present an alternative picture that allows SN1998bw to have an explosion energy and ejecta mass consistent with core collapse supernovae (although at the 'bright' end). We show that the LC of SN1998bw can be understood as result of an aspherical explosion along the rotational axis of a basically spherical, non-degenerate C/O core of massive star with an explosion energy of 2foe and a total ejecta mass of 2 M_o if it is seen from high inclinations with respect to the plane of symmetry. In this model, the high expansion velocities are a direct consequence of an aspherical explosion which, in turn, produces oblate iso-density contours. It suggests that the fundamental core-collapse explosion process itself is strongly asymmetric.Comment: 12 pages, 8 figures, latex, aas2pp4.sty, submitted to Ap

Three-Dimensional Simulations of Inflows Irradiated by a Precessing Accretion Disk in Active Galactic Nuclei: Formation of Outflows

We present three-dimensional (3-D) hydrodynamical simulations of gas flows in the vicinity of an active galactic nucleus (AGN) powered by a precessing accretion disk. We consider the effects of the radiation force from such a disk on its environment on a relatively large scale (up to ~10 pc. We implicitly include the precessing disk by forcing the disk radiation field to precess around a symmetry axis with a given period ($P$) and a tilt angle ($\Theta$). We study time evolution of the flows irradiated by the disk, and investigate basic dependencies of the flow morphology, mass flux, angular momentum on different combinations of $\Theta$ and $P$. We find the gas flow settles into a configuration with two components, (1) an equatorial inflow and (2) a bipolar inflow/outflow with the outflow leaving the system along the poles (the directions of disk normals). However, the flow does not always reach a steady state. We find that the maximum outflow velocity and the kinetic outflow power at the outer boundary can be reduced significantly with increasing $\Theta$. We also find that of the mass inflow rate across the inner boundary does not change significantly with increasing $\Theta$. (Abbreviated)Comment: Accepted for publication in ApJ. 15 pages, 7 figures. A version with full resolution figures can be downloaded from http://www.physics.unlv.edu/~rk/preprint/precess.pd

Neon Abundances from a Spitzer/IRS Survey of Wolf-Rayet Stars

We report on neon abundances derived from {\it Spitzer} high resolution spectral data of eight Wolf-Rayet (WR) stars using the forbidden line of [\ion{Ne}{3}] 15.56 microns. Our targets include four WN stars of subtypes 4--7, and four WC stars of subtypes 4--7. We derive ion fraction abundances $\gamma$ of Ne^{2+} for the winds of each star. The ion fraction abundance is a product of the ionization fraction $Q_{\rm i}$ in stage i and the abundance by number ${\cal A}_E$ of element E relative to all nuclei. Values generally consistent with solar are obtained for the WN stars, and values in excess of solar are obtained for the WC stars.Comment: to appear in Astrophysical Journa

Current status of NLTE analysis of stellar atmospheres

Various available codes for NLTE modeling and analysis of hot star spectra are reviewed. Generalizations of standard equations of kinetic equilibrium and their consequences are discussed.Comment: in Determination of Atmospheric Parameters of B-, A-, F- and G-Type Stars, E. Niemczura et al. eds., Springer, in pres

On the Mass of Population III Stars

Performing 1D hydrodynamical calculations coupled with non-equilibrium processes for H2 formation, we pursue the thermal and dynamical evolution of filamentary primordial clouds and attempt to make an estimate on the mass of population III stars. It is found that, almost independent of initial conditions, a filamentary cloud continues to collapse nearly isothermally due to H_2 cooling until the cloud becomes optically thick against the H_2 lines. During the collapse the cloud structure separates into two parts, i.e., a denser spindle and a diffuse envelope. The spindle contracts quasi-statically, and thus the line mass of the spindle keeps a characteristic value determined solely by the temperature ($\sim 800$ K). Applying a linear theory, we find that the spindle is unstable against fragmentation during the collapse. The wavelength of the fastest growing perturbation lessens as the collapse proceeds. Consequently, successive fragmentation could occur. When the central density exceeds $n_c \sim 10^{10-11} cm^{-3}$, the successive fragmentation may cease since the cloud becomes opaque against the H_2 lines and the collapse decelerates appreciably. The mass of the first star is then expected to be typically $\sim 3 M_\odot$, which may grow up to $\sim 16 M_\odot$ by accreting the diffuse envelope. Thus, the first-generation stars are anticipated to be massive but not supermassive.Comment: 23 pages, 6 figures, accepted by ApJ (April 10

The Evolution of Relativistic Binary Progenitor Systems

Relativistic binary pulsars, such as B1534+12 and B1913+16 are characterized by having close orbits with a binary separation of ~ 3 R_\sun. The progenitor of such a system is a neutron star, helium star binary. The helium star, with a strong stellar wind, is able to spin up its compact companion via accretion. The neutron star's magnetic field is then lowered to observed values of about 10^{10} Gauss. As the pulsar lifetime is inversely proportional to its magnetic field, the possibility of observing such a system is, thus, enhanced by this type of evolution. We will show that a nascent (Crab-like) pulsar in such a system can, through accretion-braking torques (i.e. the "propeller effect") and wind-induced spin-up rates, reach equilibrium periods that are close to observed values. Such processes occur within the relatively short helium star lifetimes. Additionally, we find that the final outcome of such evolutionary scenarios depends strongly on initial parameters, particularly the initial binary separation and helium star mass. It is, indeed, determined that the majority of such systems end up in the pulsar "graveyard", and only a small fraction are strongly recycled. This fact might help to reconcile theoretically expected birth rates with limited observations of relativistic binary pulsars.Comment: 24 pages, 10 Postscript figures, Submitted to The Astrophysical Journa

Physical State of Molecular Gas in High Galactic Latitude Translucent Clouds

The rotational transitions of carbon monoxide (CO) are the primary means of investigating the density and velocity structure of the molecular interstellar medium. Here we study the lowest four rotational transitions of CO towards high-latitude translucent molecular clouds (HLCs). We report new observations of the J = (4-3), (2-1), and (1-0) transitions of CO towards eight high-latitude clouds. The new observations are combined with data from the literature to show that the emission from all observed CO transitions is linearly correlated. This implies that the excitation conditions which lead to emission in these transitions are uniform throughout the clouds. Observed 13CO/12CO (1-0) integrated intensity ratios are generally much greater than the expected abundance ratio of the two species, indicating that the regions which emit 12CO (1-0) radiation are optically thick. We develop a statistical method to compare the observed line ratios with models of CO excitation and radiative transfer. This enables us to determine the most likely portion of the physical parameter space which is compatible with the observations. The model enables us to rule out CO gas temperatures greater than 30K since the most likely high-temperature configurations are 1 pc-sized structures aligned along the line of sight. The most probable solution is a high density and low temperature (HDLT) solution. The CO cell size is approximately 0.01 pc (2000 AU). These cells are thus tiny fragments within the 100 times larger CO-emitting extent of a typical high-latitude cloud. We discuss the physical implications of HDLT cells, and we suggest ways to test for their existence.Comment: 19 pages, 13 figures, 2 tables, emulateapj To be published in The Astrophysical Journa

Analysis of the Flux and Polarization Spectra of the Type Ia Supernova SN 2001el: Exploring the Geometry of the High-velocity Ejecta

SN 2001el is the first normal Type Ia supernova to show a strong, intrinsic polarization signal. In addition, during the epochs prior to maximum light, the CaII IR triplet absorption is seen distinctly and separately at both normal photospheric velocities and at very high velocities. The high-velocity triplet absorption is highly polarized, with a different polarization angle than the rest of the spectrum. The unique observation allows us to construct a relatively detailed picture of the layered geometrical structure of the supernova ejecta: in our interpretation, the ejecta layers near the photosphere (v \approx 10,000 km/s) obey a near axial symmetry, while a detached, high-velocity structure (v \approx 18,000-25,000 km/s) with high CaII line opacity deviates from the photospheric axisymmetry. By partially obscuring the underlying photosphere, the high-velocity structure causes a more incomplete cancellation of the polarization of the photospheric light, and so gives rise to the polarization peak and rotated polarization angle of the high-velocity IR triplet feature. In an effort to constrain the ejecta geometry, we develop a technique for calculating 3-D synthetic polarization spectra and use it to generate polarization profiles for several parameterized configurations. In particular, we examine the case where the inner ejecta layers are ellipsoidal and the outer, high-velocity structure is one of four possibilities: a spherical shell, an ellipsoidal shell, a clumped shell, or a toroid. The synthetic spectra rule out the spherical shell model, disfavor a toroid, and find a best fit with the clumped shell. We show further that different geometries can be more clearly discriminated if observations are obtained from several different lines of sight.Comment: 14 pages (emulateapj5) plus 18 figures, accepted by The Astrophysical Journa