Saturation of the magnetorotational instability at large Elsasser number

The magnetorotational instability is investigated within the shearing box approximation in the large Elsasser number regime. In this regime, which is of fundamental importance to astrophysical accretion disk theory, shear is the dominant source of energy, but the instability itself requires the presence of a weaker vertical magnetic field. Dissipative effects are weaker still but not negligible. The regime explored retains the condition that (viscous and ohmic) dissipative forces do not play a role in the leading order linear instability mechanism. However, they are sufficiently large to permit a nonlinear feedback mechanism whereby the turbulent stresses generated by the MRI act on and modify the local background shear in the angular velocity profile. To date this response has been omitted in shearing box simulations and is captured by a reduced pde model derived here from the global MHD fluid equations using multiscale asymptotic perturbation theory. Results from numerical simulations of the reduced pde model indicate a linear phase of exponential growth followed by a nonlinear adjustment to algebraic growth and decay in the fluctuating quantities. Remarkably, the velocity and magnetic field correlations associated with these algebraic growth and decay laws conspire to achieve saturation of the angular momentum transport. The inclusion of subdominant ohmic dissipation arrests the algebraic growth of the fluctuations on a longer, dissipative time scale. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Diffusion coefficient of neutrinos in a pre-supernova star

The diffusion coefficient of neutrinos in the trapping regime in a presupernova star,where the scattering of neutrinos by nuclei is the dominant effect, is calculated within the framework of relativistic kinetic theory in the lowest order Chapman-Enskog approximation. The screening of the nuclei by ion-ion and electron-electron correlations is also taken into account.     

Local Stability Criterion for a Self-Gravitating Gaseous Disk of Spiral Galaxies

An improved linear stability theory of small-amplitude oscillations of a self-gravitating, infinitesimally thin gaseous disk of spiral galaxies has been developed by Bertin, Lau, Lin, Mark, Morozov, Polyachenko, and others in the approximation of moderately tightly wound gravity perturbations. In this regime, the generalized Lin–Shu type dispersion relation was also found by including higher order terms in the small parameter 1/kr for wavenumber k and radius r. It was shown that in the differentially rotating disks for nonaxisymmetric (spiral)perturbations Toomre's modified critical Q-parameter is larger than the standard one: the fact that the spiral perturbations in the nonuniformly rotating system are more unstable than the axisymmetric ones is taken into account in this modified local stability criterion. We use hydrodynamical simulations to test the validity of the modified local criterion. This revised version was published online in July 2006 with corrections to the Cover Date.     

Linear perturbations in a universe with a cosmological constant

There is now evidence that the cosmological constant Λ has a non-zero positive value. Alternative scenarios to a pure cosmological constant model are provided by quintessence, an effective negative pressure fluid permeating the Universe. Recent results indicate that the energy density ρ and the pressure p of this fluid are constrained by − ρ ≤ p ≲−0.6 ρ . As p =− ρ is equivalent to the pure cosmological constant model, it is appropriate to analyse this particular, but important, case further.
We study the linear theory of perturbations in a Friedmann–Robertson–Walker universe with a cosmological constant. We obtain the equations for the evolution of the perturbations in the fully relativistic case, for which we analyse the single-fluid and two-fluid cases. We obtain solutions to these equations in appropriate limits. We also study the Newtonian approximation. We find that for a positive cosmological constant universe (i) the perturbations will grow more slowly in the relativistic regime for a two-fluid composed of dark matter and radiation, and (ii) in the Newtonian regime the perturbations stop growing.     

Weak lensing predictions at intermediate scales

As pointed out in previous studies, the measurement of the skewness of the convergence field κ will be useful in breaking the degeneracy among the cosmological parameters constrained from weak lensing observations. The combination of shot noise and finite survey volume implies that such a measurement is likely to be performed in a range of intermediate scales (0.5 to 20 arcmin) where neither perturbation theory nor the hierarchical ansatz applies. Here we explore the behaviour of the skewness of κ at these intermediate scales, based on results for the non-linear evolution of the mass bispectrum. We combined different ray-tracing simulations to test our predictions, and we find that our calculations describe accurately the transition from the weakly non-linear to the strongly non-linear regime. We show that the single lens-plane approximation remains accurate even in the non-linear regime, and we explicitly calculate the corrections to this approximation. We also discuss the prospects of measuring the skewness in upcoming weak lensing surveys.     

Galaxy clustering in the nonlinear regime

We calculate the r.m.s. density fluctuations using the ID Zel'dovich approximation implemented with a virialization scheme for the very dense fluctuations. The results are quite different from those found in the linear regime with strong implications for the normalization of the power spectrum of the initial density fluctuations. We also calculate the nonlinear two-point correlation function for matter with an initial spectrum (equivalent to white noise in 3D) of fluctuations; differences from the linear approximation appear at scales much larger than those considered previously of 10h ^–1 Mpc.Paper presented at the 11th European Regional Astronomical Meetings of the IAU on New Windows to the Universe, held 3–8 July, 1989, Tenerife, Canary Islands, Spain.     

Interference of spectral lines in thermal radiation from the lower atmosphere of Venus

The absorption spectrum and thermal radiation fluxes in the lower atmosphere of Venus are calculated using the theory of molecular state interference in the strong collision approximation. Comparison is made with the absorption and radiative transfer calculations in terms of the statistical theory of collisional line broadening and based on an empirical form factor. The calculations show that the line broadening mechanism does not affect the thermal regime of the atmosphere at heights above 60 km, but affects significantly the behavior of the greenhouse effect below the cloud layer.     

Predictions for the clustering properties of the Lyman-alpha forest – I. One-point statistics

We present predictions for the one-point probability distribution and cumulants of the transmitted QSO flux in the high redshift Lyman- α forest. We make use of the correlation between the Lyman- α optical depth and the underlying matter density predicted by gravitational instability theory and seen in numerical hydrodynamic simulations. We have modelled the growth of matter fluctuations using the non-linear shear‐free dynamics, an approximation which reproduces well the results of perturbation theory for the cumulants in the linear and weakly non-linear clustering regime. As high matter overdensities tend to saturate in spectra, the statistics of the flux distribution are dominated by weakly non-linear overdensities. As a result, our analytic approach can produce accurate predictions, when tested against N -body simulation results, even when the underlying matter field has root-mean-square fluctuations larger than unity. Our treatment can be applied to either Gaussian or non-Gaussian initial conditions. Here we concentrate on the former case, but also include a study of a specific non-Gaussian model. We discuss how the methods and predictions we present can be used as a tool to study the generic clustering properties of the Lyman- α forest at high redshift. With such an approach, rather than concentrating on simulating specific cosmological models, we may be in a position to directly test our assumptions for the Gaussian nature of the initial conditions, and the gravitational instability origin of structure itself. In a separate paper we present results for two-point statistics.     

On the post-NEWTONian approximation of theories of gravity

The physical meaning of the parameters of the post-NEWTON ian approximation is considered from a theoretical point of view. It is shown that the post-NEWTON ian approximation of any LORENTZ -covariant theory, which implies POISSON equations for every metric coefficient in the linear approximation, is the member of a three-parametric family of post-NEWTON ian metrics. Some tetrad theories are cited.     

Gravitational radiation from an isolated perfect fluid in higher multipole moments

The purpose of this paper is to give the unknown angular momentum loss of an isolated perfect fluid in any higher multipole moments in the linear approximation of general relativity theory. Also, we discuss the energy and linear momentum fluxes of the given source in higher multipole moments.     

Comment on “Neutrino oscillations in the early universe: how can large lepton asymmetry be generated?” [Astropart. Phys. 14 (2000) 79–90]

We comment on the recent paper by A.D. Dolgov, S.H. Hansen, S. Pastor and D.V. Semikoz (DHPS) [Astropart. Phys. 14 (2000) 79] on the generation of neutrino asymmetries from active–sterile neutrino oscillations. We demonstrate that the approximate asymmetry evolution equation obtained therein is an expansion, up to a minor discrepancy, of the well-established static approximation equation, valid only when the supposedly new higher order correction term is small. In the regime where this so-called “back-reaction” term is large and artificially terminates the asymmetry growth, their evolution equation ceases to be a faithful approximation to the quantum kinetic equations simply because pure Mikheyev–Smirnov–Wolfenstein (MSW) transitions have been neglected. At low temperatures the MSW effect is the dominant asymmetry amplifier. Neither the static nor the DHPS approach contains this important physics. Therefore we conclude that the DHPS results have sufficient veracity at the onset of explosive asymmetry generation, but are invalid in the ensuing low temperature epoch where MSW conversions are able to enhance the asymmetry to values of order 0.2–0.37. DHPS do claim to find a significant final asymmetry for very large δm² values. However, for this regime the effective potential they employed is not valid.     

Modelling the Propagation of a Weak Fast-Mode MHD Shock Wave near a 2D Magnetic Null Point Using Nonlinear Geometrical Acoustics

We present the results of analytical modelling of fast-mode magnetohydrodynamic wave propagation near a 2D magnetic null point. We consider both a linear wave and a weak shock and analyse their behaviour in cold and warm plasmas. We apply the nonlinear geometrical acoustics method based on the Wentzel–Kramers–Brillouin approximation. We calculate the wave amplitude, using the ray approximation and the laws of solitary shock wave damping. We find that a complex caustic is formed around the null point. Plasma heating is distributed in space and occurs at a caustic as well as near the null point due to substantial nonlinear damping of the shock wave. The shock wave passes through the null point even in a cold plasma. The complex shape of the wave front can be explained by the caustic pattern.     

Cosmological Perturbation Theory and the Spherical Collapse model — I. Gaussian initial conditions

We present a simple and intuitive approximation for solving the perturbation theory (PT) of small cosmic fluctuations. We consider only the spherically symmetric or monopole contribution to the PT integrals, which yields the exact result for tree-graphs (i.e. at leading order). We find that the non-linear evolution in Lagrangian space is then given by a simple local transformation over the initial conditions, although it is not local in Euler space. This transformation is found to be described by the spherical collapse (SC) dynamics, as it is the exact solution in the shearless (and therefore local) approximation in Lagrangian space. Taking advantage of this property, it is straightforward to derive the one-point cumulants, ξ_J, for both the unsmoothed and smoothed density fields to arbitrary order in the perturbative regime. To leading-order this reproduces, and provides us with a simple explanation for, the exact results obtained by Bernardeau. We then show that the SC model leads to accurate estimates for the next corrective terms when compared with the results derived in the exact perturbation theory making use of the loop calculations. The agreement is within a few per cent for the hierarchical ratios S _J  = ξ _J /ξ ^{J −1}₂. We compare our analytic results with N -body simulations, which turn out to be in very good agreement up to scales where σ ≈ 1. A similar treatment is presented to estimate higher order corrections in the Zel'dovich approximation. These results represent a powerful and readily usable tool to produce analytical predictions that describe the gravitational clustering of large-scale structure in the weakly non-linear regime.     

Resonant excitation of white dwarf oscillations in compact object binaries – I. The no back reaction approximation

We consider the evolution of white dwarfs with compact object companions (specifically black holes with masses up to  ∼10⁶  M_⊙  , neutron stars, and other white dwarfs). We suppose that the orbits are initially quite elliptical and then shrink and circularize under the action of gravitational radiation. During this evolution, the white dwarfs will pass through resonances when harmonics of the orbital frequency match the stellar oscillation eigenfrequencies. As a star passes through these resonances, the associated modes will be excited and can be driven to amplitudes that are so large that there is a back reaction on the orbit which, in turn, limits the growth of the modes. A formalism is presented for describing this dynamical interaction for a non-rotating star in the linear approximation when the orbit can be treated as non-relativistic. A semi-analytical expression is found for computing the resonant energy transfer as a function of stellar and orbital parameters for the regime where back reaction may be neglected. This is used to calculate the results of passage through a sequence of resonances for several hypothetical systems. It is found that the amplitude of the  ℓ= m = 2  f -mode can be driven into the non-linear regime for appropriate initial conditions. We also discuss where the no back reaction approximation is expected to fail, and the qualitative effects of back reaction.     

Self-gravitating warped discs around supermassive black holes

We consider warped equilibrium configurations for stellar and gaseous discs in the Keplerian force field of a supermassive black hole, assuming that the self-gravity of the disc provides the only acting torques. Modelling the disc as a collection of concentric circular rings and computing the torques in the non-linear regime, we show that stable, strongly warped precessing equilibria are possible. These solutions exist for a wide range of disc-to-black-hole mass ratios   M _d/ M _bh  , can span large warp angles of up to  ±∼120°  , have inner and outer boundaries, and extend over a radial range of a factor of typically two to four. These equilibrium configurations obey a scaling relation such that in good approximation     where     is the (retrograde) precession frequency and Ω is a characteristic orbital frequency in the disc. Stability was determined using linear perturbation theory and, in a few cases, confirmed by numerical integration of the equations of motion. Most of the precessing equilibria are found to be stable, but some are unstable. The main result of this study is that highly warped discs near black holes can persist for long times without any persistent forcing other than by their self-gravity. The possible relevance of this to galactic nuclei is briefly discussed.     

The linear stability of dilute particulate rings

Irregular structure in planetary rings is often attributed to the intrinsic instabilities of a homogeneous state undergoing Keplerian shear. Previously these have been analysed with simple hydrodynamic models. We instead employ a kinetic theory, in which we solve the linearised moment equations derived in Shu and Stewart 1985 for a dilute ring. This facilitates an examination of velocity anisotropy and non-Newtonian stress, and their effects on the viscous and viscous/gravitational instabilities thought to occur in Saturn's rings. Because we adopt a dilute gas model, the applicability of our results to the actual dense rings of Saturn are significantly curtailled. Nevertheless this study is a necessary preliminary before an attack on the difficult problem of dense ring dynamics. We find the Shu and Stewart formalism admits analytic stability criteria for the viscous overstability, viscous instability, and thermal instability. These criteria are compared with those of a hydrodynamic model incorporating the effective viscosity and cooling function computed from the kinetic steady state. We find the two agree in the ‘hydrodynamic limit’ (i.e., many collisions per orbit) but disagree when collisions are less frequent, when we expect the viscous stress to be increasingly non-Newtonian and the velocity distribution increasingly anisotropic. In particular, hydrodynamics predicts viscous overstability for a larger portion of parameter space. We also numerically solve the linearised equations of the more accurate Goldreich and Tremaine 1978 kinetic model and discover its linear stability to be qualitatively the same as that of Shu and Stewart's. Thus the simple collision operator adopted in the latter would appear to be an adequate approximation for dilute rings, at least in the linear regime.     

Spiral mode of standing accretion shock instability in core-collapse supernovae

The study of standing accretion shock instability (SASI) in core-collapse supernova cores has been done with three-dimensional (3D) computer simulations. Rotations with various perturbations were introduced from outer boundary of an initial steady accreting flow. We found that one or two armed spiral accreting flow onto the proto-neutron star (PNS) is formed inside the shock wave depending on perturbations. The linear growth of spiral modes are clearly diagnosed by the mode analysis of the shock surface, and the lower m modes grow quickly in the linear regime.     

Post-Newtonian satellite orbits

The first post-Newtonian approximation of general relativity is used to account for the motion of solar system bodies and near-Earth objects which are slow moving and produce weak gravitational fields. The \(n\)-body relativistic equations of motion are given by the Einstein-Infeld-Hoffmann equations. For \(n=2\), we investigate the associated dynamics of two-body systems in the first post-Newtonian approximation. By direct integration of the associated planar equations of motion, we deduce a new expression that characterises the orbit of test particles in the first post-Newtonian regime generalising the well-known Binet equation for Newtonian mechanics. The expression so obtained does not appear to have been given in the literature and is consistent with classical orbiting theory in the Newtonian limit. Further, the accuracy of the post-Newtonian Binet equation is numerically verified by comparing secular variations of known expression with the full general relativistic orbit equation.