Nonlinear Hydromagnetic Wave Support of a Stratified Molecular Cloud

We perform numerical simulations of nonlinear MHD waves in a gravitationally stratified molecular cloud that is bounded by a hot and tenuous external medium. We study the relation between the strength of the turbulence and various global properties of a molecular cloud, within a 1.5-dimensional approximation. Under the influence of a driving source of Alfvenic disturbances, the cloud is lifted up by the pressure of MHD waves and reaches a steady-state characterized by oscillations about a new time-averaged equilibrium state. The nonlinear effect results in the generation of longitudinal motions and many shock waves; however, the wave kinetic energy remains predominantly in transverse, rather than longitudinal, motions. There is an approximate equipartition of energy between the transverse velocity and fluctuating magnetic field (aspredicted by small-amplitude theory) in the region of the stratified cloud which contains most of the mass; however, this relation breaks down in the outer regions, particularly near the cloud surface, where the motions have a standing-wave character. This means that the Chandrasekhar-Fermi formula applied to molecular clouds must be significantly modified in such regions. Models of an ensemble of clouds show that, for various strengths of the input energy, the velocity dispersion in the cloud $\sigma \propto Z^{0.5}$, where $Z$ is a characteristic size of the cloud.Furthermore, $\sigma$ is always comparable to the mean Alfven velocity of the cloud, consistent with observational results.Comment: 16 pages, 15 figures, emulateapj, to appear in ApJ, 2003 Oct 1, higher resolution figures at http://www.astro.uwo.ca/~basu/pub.html or http://www.astro.uwo.ca/~kudoh/pub.htm

3-D General Relativistic MHD Simulations of Generating Jets

We have performed a first fully 3-D GRMHD simulation with Schwarzschild black hole with a free falling corona. The initial simulation results show that a jet is created as in previous axisymmetric simulations. However, the time to generate the jet is longer than in the 2-D simulations. We expect that due to the additional azimuthal dimension the dynamics of jet formation can be modified.Comment: 4 pages Proc. Oxford Radio Galaxy Workshop ed. R. Laing & K. Blundell (San Francisco: PASP) in press (revised

Thermodynamics of four-dimensional black objects in the warped compactification

We reinvestigate the thermodynamics of black objects (holes and strings) in four-dimensional braneworld models that are originally constructed by Emparan, Horowitz and Myers based on the anti-de Sitter (AdS) C-metric. After proving the uniqueness of slicing the AdS C-metric, we derive thermodynamic quantities of the black objects by means of the Euclidean formulation and find that we have no necessity of requiring any regularization to calculate their classical action. We show that there exist the Bekenstein-Hawking law and the thermodynamic first law. The thermodynamic mass of the localized black hole on a flat brane is negative, and it differs from the one previously derived. We discuss the thermodynamic stabilities and show that the BTZ black string is more stable than the localized black holes in a canonical ensemble, except for an extreme case. We also find a braneworld analogue of the Hawking-Page transition between the BTZ black string and thermal AdS branes. The localized black holes on a de Sitter brane is discussed by considering Nariai instanton, comparing the study of "black cigar" in the five-dimensional braneworld model.Comment: 15 pages, 4 figures, RevTex4, typos fixed, minor correction

Six-dimensional localized black holes: numerical solutions

To test the strong-gravity regime in Randall-Sundrum braneworlds, we consider black holes bound to a brane. In a previous paper, we studied numerical solutions of localized black holes whose horizon radii are smaller than the AdS curvature radius. In this paper, we improve the numerical method and discuss properties of the six dimensional (6D) localized black holes whose horizon radii are larger than the AdS curvature radius. At a horizon temperature $\mathcal{T} \approx 1/2\pi \ell$, the thermodynamics of the localized black hole undergo a transition with its character changing from a 6D Schwarzschild black hole type to a 6D black string type. The specific heat of the localized black holes is negative, and the entropy is greater than or nearly equal to that of the 6D black strings with the same thermodynamic mass. The large localized black holes show flattened horizon geometries, and the intrinsic curvature of the horizon four-geometry becomes negative near the brane. Our results indicate that the recovery mechanism of lower-dimensional Einstein gravity on the brane works even in the presence of the black holes.Comment: 17 pages, 9 figures, RevTeX4, typos correcte

Quadrupole formula for Kaluza-Klein modes in the braneworld

The quadrupole formula in four-dimensional Einstein gravity is a useful tool to describe gravitational wave radiation. We derive the quadrupole formula for the Kaluza-Klein (KK) modes in the Randall-Sundrum braneworld model. The quadrupole formula provides transparent representation of the exterior weak gravitational field induced by localized sources. We find that a general isolated dynamical source gives rise to the 1/r^2 correction to the leading 1/r gravitational field. We apply the formula to an evaluation of the effective energy carried by the KK modes from the viewpoint of an observer on the brane. Contrary to the ordinary gravitational waves (zero mode), the flux of the induced KK modes by the non-spherical part of the quadrupole moment vanishes at infinity and only the spherical part contributes to the flux. Since the effect of the KK modes appears in the linear order of the metric perturbations, the effective energy flux observed on the brane is not always positive, but can become negative depending on the motion of the localized sources.Comment: 9 pages, no figures, REVTeX 4; version accepted for publication in CQ

The Stability of Magnetized Rotating Plasmas with Superthermal Fields

During the last decade it has become evident that the magnetorotational instability is at the heart of the enhanced angular momentum transport in weakly magnetized accretion disks around neutron stars and black holes. In this paper, we investigate the local linear stability of differentially rotating, magnetized flows and the evolution of the magnetorotational instability beyond the weak-field limit. We show that, when superthermal toroidal fields are considered, the effects of both compressibility and magnetic tension forces, which are related to the curvature of toroidal field lines, should be taken fully into account. We demonstrate that the presence of a strong toroidal component in the magnetic field plays a non-trivial role. When strong fields are considered, the strength of the toroidal magnetic field not only modifies the growth rates of the unstable modes but also determines which modes are subject to instabilities. We find that, for rotating configurations with Keplerian laws, the magnetorotational instability is stabilized at low wavenumbers for toroidal Alfven speeds exceeding the geometric mean of the sound speed and the rotational speed. We discuss the significance of our findings for the stability of cold, magnetically dominated, rotating fluids and argue that, for these systems, the curvature of toroidal field lines cannot be neglected even when short wavelength perturbations are considered. We also comment on the implications of our results for the validity of shearing box simulations in which superthermal toroidal fields are generated.Comment: 24 pages, 12 figures. Accepted for publication in ApJ. Sections 2 and 5 substantially expanded, added Appendix A and 3 figures with respect to previous version. Animations are available at http://www.physics.arizona.edu/~mpessah/research

Multifluid magnetohydrodynamic turbulent decay

It is generally believed that turbulence has a significant impact on the dynamics and evolution of molecular clouds and the star formation which occurs within them. Non-ideal magnetohydrodynamic effects are known to influence the nature of this turbulence. We present the results of a suite of 512-cubed resolution simulations of the decay of initially super-Alfvenic and supersonic fully multifluid MHD turbulence. We find that ambipolar diffusion increases the rate of decay of the turbulence while the Hall effect has virtually no impact. The decay of the kinetic energy can be fitted as a power-law in time and the exponent is found to be -1.34 for fully multifluid MHD turbulence. The power spectra of density, velocity and magnetic field are all steepened significantly by the inclusion of non-ideal terms. The dominant reason for this steepening is ambipolar diffusion with the Hall effect again playing a minimal role except at short length scales where it creates extra structure in the magnetic field. Interestingly we find that, at least at these resolutions, the majority of the physics of multifluid turbulence can be captured by simply introducing fixed (in time and space) resistive terms into the induction equation without the need for a full multifluid MHD treatment. The velocity dispersion is also examined and, in common with previously published results, it is found not to be power-law in nature.Comment: 16 pages, 15 figures, Accepted for publication in Ap

Strong Brane Gravity and the Radion at Low Energies

For the 2-brane Randall-Sundrum model, we calculate the bulk geometry for strong gravity, in the low matter density regime, for slowly varying matter sources. This is relevant for astrophysical or cosmological applications. The warped compactification means the radion can not be written as a homogeneous mode in the orbifold coordinate, and we introduce it by extending the coordinate patch approach of the linear theory to the non-linear case. The negative tension brane is taken to be in vacuum. For conformally invariant matter on the positive tension brane, we solve the bulk geometry as a derivative expansion, formally summing the `Kaluza-Klein' contributions to all orders. For general matter we compute the Einstein equations to leading order, finding a scalar-tensor theory with $\omega(\Psi) \propto \Psi / (1 - \Psi)$, and geometrically interpret the radion. We comment that this radion scalar may become large in the context of strong gravity with low density matter. Equations of state allowing $(\rho - 3 P)$ to be negative, can exhibit behavior where the matter decreases the distance between the 2 branes, which we illustrate numerically for static star solutions using an incompressible fluid. For increasing stellar density, the branes become close before the upper mass limit, but after violation of the dominant energy condition. This raises the interesting question of whether astrophysically reasonable matter, and initial data, could cause branes to collide at low energy, such as in dynamical collapse.Comment: 24 pages, 3 figure