Analytical potential–density pairs for flat rings and toroidal structures

The Kuzmin–Toomre family of discs is used to construct potential–density pairs that represent flat ring structures in terms of elementary functions. Systems composed of two concentric flat rings, a central disc surrounded by one ring and a ring with a centre of attraction are also presented. The circular velocity of test particles and the epicyclic frequency of small oscillations about circular orbits are calculated for these structures. A few examples of three-dimensional potential–density pairs of 'inflated' flat rings (toroidal mass distributions) are presented.     

Exact potential–density pairs for flattened dark haloes

Cosmological simulations suggest that dark matter haloes are not spherical, but typically moderately to strongly triaxial systems. We investigate methods to convert spherical potential–density pairs into axisymmetric ones, in which the basic characteristics of the density profile (such as the slope at small and large radii) are retained. We achieve this goal by replacing the spherical radius r by an oblate radius m in the expression of the gravitational potential  Φ( r )  .
We extend and formalize the approach pioneered by Miyamoto & Nagai to be applicable to arbitrary potential–density pairs. Unfortunately, an asymptotic study demonstrates that, at large radii, such models always show a   R ⁻³  disc superposed on a smooth roughly spherical density distribution. As a result, this recipe cannot be used to construct simple flattened potential–density pairs for dynamical systems such as dark matter haloes. Therefore, we apply a modification of our original recipe that cures the problem of the discy behaviour. An asymptotic analysis now shows that the density distribution has the desired asymptotic behaviour at large radii (if the density falls less rapidly than   r ⁻⁴  ). We also show that the flattening procedure does not alter the shape of the density distribution at small radii: while the inner density contours are flattened, the slope of the density profile is unaltered.
We apply this recipe to construct a set of flattened dark matter haloes based on the realistic spherical halo models by Dehnen & McLaughlin. This example illustrates that the method works fine for modest flattening values, whereas stronger flattening values lead to peanut-shaped density distributions.     

Realistic triaxial density–potential–force profiles for stellar systems and dark matter haloes

Popular models for describing the luminosity–density profiles of dynamically hot stellar systems (e.g. Jaffe, Hernquist, Dehnen) were constructed to match the deprojected form of de Vaucouleurs' R ^1/4 light-profile. However, we now know that elliptical galaxies and bulges display a mass-dependent range of structural profiles. To compensate this, the model in Terzić & Graham was designed to closely match the deprojected form of Sérsic R ^{1/ n} light-profiles, including deprojected exponential light-profiles and galaxies with partially depleted cores. It is thus applicable for describing bulges in spiral galaxies, dwarf elliptical galaxies, both 'power-law' and 'core' elliptical galaxies, also dark matter haloes formed from Λ cold dark matter cosmological simulations. In this paper, we present a new family of triaxial density–potential–force triplets, which generalizes the spherical model reported in Terzić & Graham to three dimensions. If the (optional) power-law core is present, it is a five-parameter family, while in the absence of the core it reduces to three parameters. The isodensity contours in the new family are stratified on confocal ellipsoids and the potential and forces are expressed in terms of integrals which are easy to evaluate numerically. We provide the community with a suite of numerical routines for orbit integration, which feature: optimized computations of potential and forces for this family; the ability to run simulations on parallel platforms; and modular and easily editable design.     

New biorthogonal potential–density basis functions

We use the weighted integral form of spherical Bessel functions and introduce a new analytical set of complete and biorthogonal potential–density basis functions. The potential and density functions of the new set have finite central values and they fall off, respectively, similar to   r ^{−(1+ l )}  and   r ^{−(4+ l )}  at large radii, where l is the latitudinal quantum number of spherical harmonics. The lowest order term associated with   l = 0  is the perfect sphere of de Zeeuw. Our basis functions are intrinsically suitable for the modelling of three-dimensional, soft-centred stellar systems and they complement the basis sets of Clutton-Brock, Hernquist & Ostriker and Zhao. We test the performance of our functions by expanding the density and potential profiles of some spherical and oblate galaxy models.     

Triaxial stellar systems following the r1/n luminosity law: an analytical mass–density expression, gravitational torques and the bulge/disc interplay

We have investigated the structural and dynamical properties of triaxial stellar systems whose surface brightness profiles follow the   r ^{1/ n}   luminosity law – extending the analysis by Ciotti, who explored the properties of spherical   r ^{1/ n}   systems. A new analytical expression that accurately reproduces the spatial (i.e., deprojected) luminosity density profiles (error less than 0.1 per cent) is presented for detailed modelling of the Sérsic family of luminosity profiles. We evaluate both the symmetric and the non-axisymmetric components of the gravitational potential and force, and compute the torques as a function of position. For a given triaxiality, stellar systems with smaller values of n have a greater non-axisymmetric gravitational field component . We also explore the strength of the non-axisymmetric forces produced by bulges with differing n and triaxiality on systems having a range of bulge-to-disc ratios. The increasing disc-to-bulge ratio with increasing galaxy type (decreasing n ) is found to greatly reduce the amplitude of the non-axisymmetric terms, and therefore reduce the possibility that triaxial bulges in late-type systems may be the mechanism or perturbation for non-symmetric structures in the disc.
Using seeing-convolved   r ^{1/ n}   -bulge plus exponential-disc fits to the K -band data from a sample of 80 nearby disc galaxies, we probe the relations between galaxy type, Sérsic index n and the bulge-to-disc luminosity ratio. These relations are shown to be primarily a consequence of the relation between n and the total bulge luminosity. In the K band, the trend of decreasing bulge-to-disc luminosity ratio along the spiral Hubble sequence is predominantly, though not entirely, a consequence of the change in the total bulge luminosity; the trend between the total disc luminosity and Hubble type is much weaker.     

Simple potential–density pairs for flat rings

Potential–density pairs representing flat-ring structures in terms of elementary functions are presented. Structures representing one or several concentric flat rings, and discs surrounded by concentric flat rings are examined. The stability of concentric circular orbits of particles moving on a flat-ring structure is analyzed for radial perturbations.     

N-body simulations using customized potential–density pair basis sets

Potential–density pair basis sets can be used for highly efficient N -body simulation codes, but they suffer from a lack of versatility, i.e. a basis set has to be constructed for each different class of stellar system. We present numerical techniques for generating a biorthonormal potential–density pair basis set that has a general specified pair as its lowest-order member. We go on to demonstrate how the set can be used to construct N -body equilibria, which we then evolve using an N -body code that calculates forces using the basis set.     

Integrable models of galactic discs with double nuclei

We introduce a new class of 2D mass models, whose potentials are of Stäckel form in elliptic coordinates. Our model galaxies have two separate strong cusps that form double nuclei. The potential and surface density distributions are locally axisymmetric near the nuclei and become highly non-axisymmetric outside the nucleus. The surface density diverges toward the cuspy nuclei with the law     Our model is sustained by four general types of regular orbits: butterfly , nucleophilic banana , horseshoe and aligned loop orbits. Horseshoes and nucleophilic bananas support the existence of cuspy regions. Butterflies and aligned loops control the non-axisymmetric shape of outer regions. Without any need for central black holes, our distributed mass models resemble the nuclei of M31 and NGC 4486B. It is also shown that the self-gravity of the stellar disc can prevent the double nucleus to collapse.     

Density cusps: restrictions on non-axisymmetric models

Galactic nuclei are now generally thought to have density cusps in their centres, and the evidence is mounting that as a consequence they are unlikely to be triaxial. Self-consistent stellar dynamical models of non-axisymmetric cusps would be an interesting counter-argument to this conclusion. We consider 2D analogues of triaxial cusps: a sequence of non-axisymmetric, cuspy discs first described by Sridhar & Touma. Scale-free models with potential Φ ∝  r ^α are examined in detail. It is shown analytically for 0 < α ≲ 0.43 that self-consistent models with positive phase-space density do not exist. Numerical solutions of the combined Vlasov and Poisson equations suggest that the whole sequence of models with 0 < α < 1 are also unphysical. Together with existing work on cusps, we conclude on purely theoretical grounds that galactic nuclei are not expected to be triaxial.     

Non-Axisymmetric Oscillations of Thin Prominence Fibrils

We study non-axisymmetric oscillations of thin prominence fibrils. A fibril is modeled by a straight thin magnetic tube with the ends frozen in dense plasmas. The density inside and outside the tube varies only along the tube and it is discontinuous at the tube boundary. Making a viable assumption that the tube radius is much smaller than its length, we show that the squares of the frequencies of non-axisymmetric tube oscillations are given by the eigenvalues of the Sturm–Liouville problem for a second-order ordinary differential equation on a finite interval with the zero boundary conditions. For an equilibrium density that is constant outside the tube and piecewise constant inside we derived a simple dispersion equation determining the frequencies of non-axisymmetric oscillations. We carry out a parametric study of this equation both analytically and numerically, restricting our analysis to the first even mode and the first odd mode. In particular, we obtained a criterion that allows to find out if each of these modes is a normal or leaky mode.     

A dynamo model for axisymmetric and non-axisymmetric solar magnetic fields

More and more observations are showing a relatively weak, but persistent, non-axisymmetric magnetic field co-existing with the dominant axisymmetric field on the Sun. Its existence indicates that the non-axisymmetric magnetic field plays an important role in the origin of solar activity. A linear non-axisymmetric  α²– Ω  dynamo model is derived to explore the characteristics of the axisymmetric  ( m = 0)  and the first non-axisymmetric  ( m = 1)  modes and to provide a theoretical basis with which to explain the 'active longitude', 'flip-flop' and other non-axisymmetric phenomena. The model consists of an updated solar internal differential rotation, a turbulent diffusivity varying with depth, and an α-effect working at the tachocline in a rotating spherical system. The difference between the  α²–Ω  and the  α–Ω  models and the conditions that favour the non-axisymmetric modes under solar-like parameters are also presented.     

Stability of power-law discs — I. The Fredholm integral equation

The power-law discs are a family of infinitesimally thin, axisymmetric stellar discs of infinite extent. The rotation curve can be rising, falling or flat. The self-consistent power-law discs are scale-free, so that all physical quantities vary as a power of radius. They possess simple equilibrium distribution functions depending on the two classical integrals, energy and angular momentum. While maintaining the scale-free equilibrium force law, the power-law discs can be transformed into cut-out discs by preventing stars close to the origin (and sometimes also at large radii) from participating in any disturbance. This paper derives the homogeneous Fredholm integral equation for the in-plane normal modes in the self-consistent and the cut-out power-law discs. This is done by linearizing the collisionless Boltzmann equation to find the response density corresponding to any imposed density and potential. The normal modes — that is, the self-consistent modes of oscillation — are found by requiring the imposed density to equal the response density. In practice, this scheme is implemented in Fourier space, by decomposing both imposed and response densities in logarithmic spirals. The Fredholm integral equation then relates the transform of the imposed density to the transform of the response density. Numerical strategies to solve the integral equation and to isolate the growth rates and the pattern speeds of the normal modes are discussed.     

Interface modes and their instabilities in accretion disc boundary layers

We study global non-axisymmetric oscillation modes trapped near the inner boundary of an accretion disc. Observations indicate that some of the quasi-periodic oscillations (QPOs) observed in the luminosities of accreting compact objects (neutron stars, black holes and white dwarfs) are produced in the innermost regions of accretion discs or boundary layers. Two simple models are considered in this paper. The magnetosphere–disc model consists of a thin Keplerian disc in contact with a uniformly rotating magnetosphere with and low plasma density, while the star–disc model involves a Keplerian disc terminated at the stellar atmosphere with high density and small density scaleheight. We find that the interface modes at the magnetosphere–disc boundary are generally unstable due to Rayleigh–Taylor and/or Kelvin–Helmholtz instabilities. However, differential rotation of the disc tends to suppress Rayleigh–Taylor instability, and a sufficiently high disc sound speed (or temperature) is needed to overcome this suppression and to attain net mode growth. On the other hand, Kelvin–Helmholtz instability may be active at low disc sound speeds. We also find that the interface modes trapped at the boundary between a thin disc and an unmagnetized star do not suffer Rayleigh–Taylor or Kelvin–Helmholtz instability, but can become unstable due to wave leakage to large disc radii and, for sufficiently steep disc density distributions, due to wave absorption at the corotation resonance in the disc. The non-axisymmetric interface modes studied in this paper may be relevant to the high-frequency QPOs observed in some X-ray binaries and in cataclysmic variables.     

Contraction and Fragmentation of Magnetized Rotating Clouds and Formation of Binary Systems

Using three-dimensional (3D) magnetohydrodynamical (MHD) nested-grid simulations, the fragmentation of a rotating magnetized molecular cloud core is studied. An isothermal rotating magnetized cylindrical cloud in hydrostatic balance is considered. We studied non-axisymmetric evolution of the cloud. It is found that non-axisymmetry hardly evolves in the early phase, but it begins to grow after the gas contracts and forms a thin disk. The disk formation and thus growth of non-axisymmetric perturbation are strongly promoted by rotation and magnetic field strength. We found two types of fragmentations: fragmentation from a ring and that from a bar. These two types of fragmentations occur in thin adiabatic cores with the thickness being smaller than 1/4 of the radial size. For the fragments to survive, they should be formed in a heavily elongated barred core or a flat round disk. In the models showing fragmentation, outflows from respective fragments are found as well as that driven by the rotating bar or the disk.     

The non-axisymmetrical configuration of a large-scale solar and stellar magnetic field

The non-axisymmetric and nonlinear solutions of the magnetostatic equations are given in three-dimensional space of spherical coordinates (r, θ, ?). These solutions are applied to the large-scale solar magnetic field. Their basic features are similar to a dipole field near the polar regions and the polarity reverses near the equator. These features agree with observations for the large-scale solar magnetic field. The solutions can also be applied to investigating the connection between the structure of the magnetic field and the density distribution of the corona. It is shown that the tops of the closed magnetic field associate with density enhancements. Similar results may apply to the large-scale configuration of the stellar field.     

The Ω dependence in equations of motion

The equations of motion governing the evolution of a collisionless gravitating system of particles in an expanding universe can be cast in a form which is almost independent of the cosmological density parameter, Ω, and the cosmological constant, Λ. The new equations are expressed in terms of a time variable τ=ln D , where D is the linear rate of growth of density fluctuations. The dependence on the density parameter is proportional to ε=Ω^−0.2−1 times the difference between the peculiar velocity (with respect to τ) of particles and the gravity field (minus the gradient of the potential); or, before shell-crossing, times the sum of the density contrast and the velocity divergence. In a one-dimensional collapse or expansion, the equations are fully independent of Ω and Λ before shell crossing. In the general case, the effect of this weak Ω dependence is to enhance the rate of evolution of density perturbations in dense regions. In a flat universe with Λ7ne;0, this enhancement is less pronounced than in an open universe with Λ=0 and the same Ω. Using the spherical collapse model, we find that the increase of the rms density fluctuations in a low-Ω universe relative to that in a flat universe with the same linear normalization is ∼0.01ε(Ω)〈δ³〉, where δ is the density field in the flat universe. The equations predict that the smooth average velocity field scales like Ω^0.6, while the local velocity dispersion (rms value) scales, approximately, like Ω^0.5. High-resolution N -body simulations confirm these results and show that density fields, when smoothed on scales slightly larger than clusters, are insensitive to the cosmological model. Haloes in an open model simulation are more concentrated than haloes of the same M /Ω in a flat model simulation.     

On non-axisymmetric magnetic equilibria in stars

In previous work, stable approximately axisymmetric equilibrium configurations for magnetic stars were found by numerical simulation. Here, I investigate the conditions under which more complex, non-axisymmetric configurations can form. I present numerical simulations of the formation of stable equilibria from turbulent initial conditions and demonstrate the existence of non-axisymmetric equilibria consisting of twisted flux tubes lying horizontally below the surface of the star, meandering around the star in random patterns. Whether such a non-axisymmetric equilibrium or a simple axisymmetric equilibrium forms depends on the radial profile of the strength of the initial magnetic field. The results could explain observations of non-dipolar fields on stars such as the B0.2 main-sequence star τ Sco or the pulsar 1E 1207.4-5209. The secular evolution of these equilibria due to Ohmic and buoyancy processes is also examined.     

New agegraphic dark energy in <Emphasis Type="Italic">f</Emphasis>(<Emphasis Type="Italic">R</Emphasis>) gravity

In this paper, we study a cosmological application of the new agegraphic dark energy density in the f(R) gravity framework. We employ the new agegraphic model of dark energy to obtain the equation of state for the new agegraphic energy density in a spatially flat universe. Our calculations show, taking n<0, that it is possible to have w _Λ crossing −1. This implies that one can generate a phantom-like equation of state from a new agegraphic dark energy model in a flat universe in the modified gravity cosmology framework. Also, we develop a reconstruction scheme for the modified gravity with f(R) action.     

On the electric field screening by electron–positron pairs in a pulsar magnetosphere

We present a steady one-dimensional model for a pulsar polar cap accelerator, where the field-aligned electric field and flow are solved self-consistently with a given current density. It is assumed that no particles return to the star. It is known that the space-charge-limited flow is accelerated to energies high enough to create electron–positron pairs if the assumed current density is high enough. We find that when pairs are created in such a space-charge-limited flow, the accelerating electric field is screened out within a short distance after pair creation, if the pair particle flux is larger than a critical value. We also find that a space charge density wave is excited in the screening region.
We find that a pair flux larger than the critical value M _c=10³–10⁵ must be reached in a layer with thickness equal to the braking distance for the decelerating component. Therefore, the required multiplicity – the number of pairs created by one primary particle – is too large to be realized in the actual pulsar magnetosphere. We suggest that in order to obtain a localized potential drop along the polar cap magnetic flux, one needs to take into account additional effects such as wave–particle interaction or quasi-periodic pair creation.     

Investigations of spherical kinematic mean-field dynamo models

The paper supplements an earlier one on the mean-field approach to spherical kinematic dynamo models (Rädler 1980a) by results of numerical investigations. A number of dynamo models working on the basis of the α²-mechanism are considered. Cases of pure α²-mechanism are studied, which includes only the simplest form of α-effect and no other induction effect, as well as cases with several additional effects due to fluctuating or mean motions. By the pure α²-mechanism axisymmetric and non-axisymmetric fields, can be excited and maintained with nearly equal ease. Part of the additional induction effects, however, clearly favour axisymmetric fields, and others non-axisymmetric fields. The non-axisymmetric fields are waves which travel in azimuthal direction, eastward or westward, depending on the models. For special dynamo models the transition from α² to αω-mechanism and properties of the latter are investigated. The results support the presumption that the αω-mechanism is able to maintain only axisymmetric but never non-axisymmetric fields. Conditions for the occurrence of non-oscillatory or oscillatory fields are discussed, and again the influence of additional induction effects is studied. There are further presented a model with βω-mechanism maintaining an axisymmetric non-oscillatory field, and models with two kinds of δω-mechanisms allowing axisymmetric non-oscillatory and oscillatory fields. Some ideas concerning dynamo models for the Earth, the Sun and magnetic stars are discussed. It seems possible to construct dynamo models for the Earth, on the basis of the α²-mechanism which explain not only the presence of a magnetic field with a strong dipole part but also the inclination of the dipole axis against the axis of rotation, the occurrence of higher multipoles and the westward drift of the non-axisymmetric parts. Models with αω, βω or δω-mechanism, which have to be considered in the case of a strong differential rotation inside the core, provide an explanation at first only of the axisymmetric parts of the field, and the non-axisymmetric parts have then to be interpreted, for example, as MAC-waves. As far as dynamo models for the Sun are concerned, in addition to the possibility of an αω-mechanism also that of a βω or δω-mechanism is discussed, which, however, does not look not very promising. In the models developed so far, which work with the αω-mechanism, only a few of the induction effects of fluctuating motions have been included; it seems necessary to investigate also influences of other effects. The sectorial structure of the solar magnetic field can hardly be understood in terms of the traditional mean-field concept. The magnetic stars possess fields which strongly deviate from symmetry with respect to the axis of rotation. The occurrence of such fields seems understandable only if there is no noticeable differential rotation. They can be maintained by an α²-mechanism but hardly by αω, βω or δω-mechanisms.