A numerical solution of the dynamical instability of a thin accretion torus

We consider the non-axisymmetric, dynamical instability of a thin accretion torus with a non-zero shift of corotation radius. By numerical method we evaluated the wave number dependence of the linear rate of growth of instability and the co-rotation shift. The rate of growth is only slightly affected by the non-zero co-rotation shift, while the dispersion relation in the case of a shift is the same as that of the linear KdV equation. This shows that the “planet-like” solution found in numerical simulations of thin tori is very probably analogous to the soliton solution of the KdV equation.     

Galactic dynamo with helicity fluxes

We construct an approximation for the magnetic field of galaxies that takes into account the magnetic helicity conservation law. In our calculations, we use the fact that the galactic disk is fairly thin and, therefore, the magnetic field component perpendicular to the galactic disk can be neglected (the so-called no-z approximation). However, an averaging of the magnetic field over the entire galaxy, as was done in previous works, is not performed. Our results are compared both with the approximation that disregards the helicity flux and with the results obtained in models with helicity fluxes but without averaging. We show that, compared to the classical model, there are a number of new effects (for example, magnetic field oscillations) and, compared to the model with averaging, the behavior of the magnetic field “softens”: its stationary value is reached more slowly and the oscillation amplitude decreases. This is because the dissipative processes changing the magnetic field growth rate are taken into account in our model. In contrast to the model with averaging, here it becomes possible to construct the dependence of the magnetic field and helicity on the distance from the galactic center.     

The Effect of Radiative Efficiency on the Growth of the Black Hole Mass

The effect of variation of radiative efficiency on the growth of the black hole mass is discussed. In the process of accretion, the black hole's angular momentum varies, resulting in a variation of the radiative efficiency, and the growth of the black hole mass is thereby affected. For the exponential growth model of back holes and taking into account the effect of a varying radiative efficiency, the equation for the growth of the black hole mass is solved numerically. Compared to the model that assumes a constant radiative efficiency, the result indicates that variation of radiative efficiency has a marked, retarding effect on the growth of the black hole mass. This model can explain quantitatively some recent observational results.     

Ionization structure of the atmospheres and line profiles in the spectra of Wolf-Rayet stars

The ionization structure of the atmospheres of Wolf-Rayet (WR) and WC stars is studied. The stellar atmospheres were assumed to consist of helium, hydrogen, and carbon. Profiles of the C III l 5696 line are calculated, both for a spherically symmetric atmosphere with a density that decreases monotonically outward and for an atmosphere containing a dense condensation (inhomogeneity). The dependence of line profiles on the parameters of the inhomogeneity is investigated. It is shown that profiles of the C III λ 5696 line calculated assuming no inhomogeneities in the atmosphere are too weak, whereas assuming the existence of inhomogeneities enables one to reconcile the observed and calculated profiles. An equation is obtained relating the mass of an inhomogeneity to the flux in the detail of the total profile of the CIII λ 5696 line formed by that inhomogeneity. This equation is used to construct a stochastic cloud model of the atmosphere of a WR star, consisting of a large number of inhomogeneities in a homogeneous, spherically symmetric stellar wind. In the proposed model, the formation of inhomogeneities was treated as a random process. It is shown that in this model it is possible both to obtain an average line profile corresponding to the observed one and to reproduce the amplitude and overall pattern of variability of profiles in the spectra of Wolf-Rayet stars. Translated from Astrofizika, Vol. 42, No. 3, pp. 373–398, July–September, 1999.     

To bin or not to bin: Decorrelating the cosmic equation of state

The physics behind the acceleration of the cosmic expansion can be elucidated through comparison of the predictions of dark energy equations of state to observational data. In seeking to optimize this, we investigate the advantages and disadvantages of using principal component analysis, uncorrelated bandpowers and the equation of state within redshift bins. We demonstrate that no one technique is a panacea, with tension between clear physical interpretation from localization and from decorrelated errors, as well as model dependence and form dependence. Specific lessons include the critical role of proper treatment of the high redshift expansion history and the lack of a unique, well-defined signal to noise or figure of merit.     

Cosmology without Big Bang within the framework of the Projective Unified Field Theory and the influence on the motion of cosmic objects

The Projective Unified Field Theory having been developed by the author since 1957 is applied to a closed homogeneous isotropic cosmological model. By postulating a (new) conservation law for the scalaric substrate mass (derived from the balance equation of matter and leading to a basically new view on the concept of mass in context with the Mach principle) we obtained a coupled system of differential equations for the world radius and the scalaric world function, which could be treated numerically by Maple. Detailed calculations show that the big bang singularity of the Einstein theory can be avoided. Since the mass density and the temperature of the background radiation exhibit maxima at the same time after the cosmological “big start”, this fact could be of interest for cosmogonic activities. Within the framework of this theory there are hints for an antigravitational world era with repulsive forces after the big start. Furthermore, the balance equations for the energy and the angular momentum of a body as well as the equation of motion of a body (with a series of consequences) are derived. In this context we found a formula for the time dependence of the “effective Newtonian gravitational constant”. Further results refer to certain aspects for understanding the observed rotational curves of cosmic objects within galaxies as well as to the conservation of the number of photons and baryons and their mutual ratio etc.     

Consequences of an isotropic static plasma sheet in models of the geomagnetic tail

The equation of momentum balance and magnetic flux conservation are given for a static tail model with an isotropic plasma sheet. The possibility of magnetic field leakage into the solar wind and across the neutral sheet is allowed. Numerical integrations for a wide variety of adjustable model parameters are presented that give the dependence on distance from Earth of all tail parameters (field strength inside and outside of the plasma sheet, plasma pressure, plasma sheet area, tail radius, and normal field component to the neutral sheet). The model gives good agreement with the observed distance dependence of the tail field strength, and accounts for the scatter in the data in terms of a mixture of the fields inside and outside the plasma sheet in the data averages. However, compared with the present interpretations of the observations the model gives a too large plasma pressure at large distances and a too small normal component to the neutral sheet. The discrepancies imply that plasma flow and/or pressure anisotropy are required for an adequate model.     

Formation of spherules in impact produced vapor plumes

We have constructed a numerical model of spherule formation in an impact produced vapor plume. This model tracks the expansion of the vapor plume using a one-dimensional Lagrangian hydrocode coupled with the ANEOS equation of state for silica. We then include the equations for nucleation and growth as described by homogeneous nucleation theory to describe the process of spherule formation. We use this model to determine the number and size of the spherules that an impact creates. We also explore when and where spherules are formed in the vapor plume, and how this affects the size of the spherules. In general we find that smaller spherules form in the outer, faster moving, portions of the vapor plume at earlier times. This work also explores the effect of impactor size and impact velocity on the resultant spherule size. We report a simple linear dependence on impactor size and a complex dependence on impact velocity. We find that a 10 km diameter asteroid impacting at a velocity of ～21 km/s creates spherules that are ～250 μm in diameter which is comparable to the spherules found in the K/Pg boundary layer.     

Inversion of crater morphometric data to gain insight on the cratering process

Abstract— In recent years, morphometric data for Venus and several outer planet satellites have been collected, so we now have observational data of complex craters formed in a large range of target properties. We present general inversion techniques that can utilize the morphometric data to quantitatively test various models of complex crater formation. The morphometric data we use in this paper are depth of a complex crater, the diameter at which the depth-diameter ratio changes, and onset diameters for central peaks, terraces, and peak rings. We tested the roles of impactor velocities and hydrostatic pressure vs. crustal strength, and we tested the specific models of acoustic fluidization (Melosh, 1982) and nonproportional growth (Schultz, 1988). Neither the acoustic fluidization model nor the nonproportional growth in their published formulations are able to successfully reproduce the data. No dependence on impactor velocity is evident from our inversions. Most of the morphometric data is consistent with a linear dependence on the ratio of crustal strength to hydrostatic pressure on a planet, or the factor c/p g.     

The origin of scaling in the galaxy distribution

I present an analytic model for non‐linear clustering of the luminous (baryonic) material in a universe in which the gravitational field is dominated by dark matter. The model is based on a two-component generalization of the adhesion approximation in which the gravitational potential of the dark component is determined by the standard Zel'dovich approximation or one of its variants, or by an N ‐body simulation. The baryonic matter flow is dissipative and is driven by this dark matter gravitational potential. The velocity potential of the matter is described by a generalization of the Burgers equation: the random heat equation ('RH equation') with a spatially correlated Gaussian driving potential.
The properties of the RH equation are well understood: it is closely related to the equation for the Anderson model and to Brownian motion in a random potential: the solution can be expressed in terms of path integrals. Using this it is possible to derive the scaling properties of the solution and, in particular, those of the resultant velocity field. Even though the flow is non‐linear, the velocity field remains Gaussian and inherits its scaling properties from the gravitational potential. This provides an underlying dynamical reason why the density field in the baryonic component is lognormally distributed and manifests multifractal scaling.
By explicitly putting dark and luminous matter on different footings, the model provides an improved framework for considering the growth of large‐scale cosmic structure. It provides a solution for the velocity potential of the baryonic component in closed form (albeit a path integral) from which the statistical properties of the baryonic flow can be derived.     

Galaxies in the connection machine

Many problems in galactic dynamics require computer speeds orders of magnitude larger than what can be achieved on current single-processor machines. In the near future such speeds are likely to become available through computer architectures based on large-scale, fine-grained parallelism. An example of a highly parallel computer is the Connection Machine, with up to 65,636 processors. We have benchmarked gravitationalN-body algorithms on the Connection Machine, and compared those with similar benchmarks which we have obtained on more traditional vector supercomputers. Our conclusions are: (1) The direct summation algorithm, with of orderN ² interactions forN particles, can be made to run with high efficiency on either type of computer. As a result, the Connection Machine clearly wins in speed over all supercomuters tested, with the exception of an 8-processor ETA, which shows a comparable performance. (2) A more efficient tree algorithm reduces the growth of the number of interactions fromN ² toN logN. However, the greater complexity of this algorithm causes a considerable degradation of efficiency, by a factor which is larger on the Connection Machine than on vector supercomputers. As a result, out tree code runs at comparable speeds on both types of machines, with the notable exception of the 8-processor ETA, which has an extrapolated speed for running our tree code which is higher than any of the other machines we have tested.     

A comparison of the temperature and density structure in high and low speed thermal proton flows

The continuity, momentum and energy hydrodynamic equations for an H⁺-O⁺ topside ionosphere have been solved self-consistently for steady state conditions similar to those found outside the plasmasphere. Results are given for undisturbed and trough conditions with a range of H⁺ outflow velocities yielding subsonic and supersonic flow. In the formulation of the equations, account was taken of the velocity dependence of ion-neutral, ion-ion and ion-electron collision frequencies. In addition, parallel stress and the nonlinear acceleration term were retained in the H⁺ momentum equation. Results computed from this model show that, as a result of Joule (frictional) heating, the H⁺ temperature rises with increasing outflow velocity in the subsonic flow regime, reaching a maximum value of about 4000 K. For supersonic flow other terms in the H⁺ momentum equation become important and alter the H⁺ velocity profile such that convection becomes a heat sink in the 1000–1500 km altitude range. This, together with the reduced Joule heating resulting from the high-speed velocity dependence of the H⁺ collision frequencies, results in a decrease in the H⁺ temperature as the outflow velocity increases. However, for all outward flows the H⁺ temperature remains substantially greater than the O⁺ temperature. With identical upper boundary velocities, the H⁺ flow velocity is higher at low altitudes for trough conditions compared with non-trough conditions, but the H⁺ temperature in the trough is lower. The form of the H⁺ density profiles for supersonic flow does not in general differ greatly from those obtained with wholly subsonic flow conditions.     

Some factors affecting the growth and decay of plages

The Mount Wilson coarse array magnetograph data set is analyzed to examine the dependence of growth and decay rates on the tilt angles of the magnetic axes of the regions. It is found that there is a relationship between these quantities which is similar to that found earlier for sunspot groups. Regions near the average tilt angle show larger average (absolute) growth and decay rates. Thepercentage growth and decay rates show minima (in absolute values) at the average tilt angles because the average areas of regions are largest near this angle. This result is similar to that derived earlier for sunspot groups. As in the case of spot groups, this suggests that, for decay, the effect results from the fact that the average tilt angle may represent the simplest subsurface configuration of the flux loop or loops that make up the region. In the case of region growth, it was suggested that the more complicated loop configuration should result in increased magnetic tension in the flux loop, and thus in a slower ascent of the loop to the surface, and thus a slower growth rate. In order to examine this further, the growth and decay rates of plage regions were examined as functions of the magnetic complexity of the regions. In the case of decay, the result was as expected from the model suggested above - that is, the more complex regions decayed more slowly. But for growing regions the effect is the opposite to that expected (more complex regions grow faster, even in terms of percentage growth), so the explanation of the tilt angle effect for growing regions proposed earlier may not be valid.Operated by the Association of Universities for Research in Astronomy, Inc., under Cooperative Agreement with the National Science Foundation.     

Some doubts on the validity of the foreground Galactic contribution subtraction from microwave anisotropies

The Galactic foreground contamination in CMBR anisotropies, especially from the dust component, is not easily separable from the cosmological or extragalactic component. In this paper, some doubts will be raised concerning the validity of the methods used until now to remove Galactic dust emission and will show that none of them achieves its goal. First, I review the recent bibliography on the topic and discuss critically the methods of foreground subtraction: the cross-correlation with templates, analysis assuming the spectral shape of the Galactic components, the “maximum entropy method”, “internal linear combination”, and “wavelet-based high resolution fitting of internal templates”. Second, I analyse the Galactic latitude dependence from WMAP data. The frequency dependence is discussed with data in the available literature. The result is that all methods of subtracting the Galactic contamination are inaccurate. The Galactic latitude dependence analysis or the frequency dependence of the anisotropies in the range 50–250 GHz put a constraint on the maximum Galactic contribution in the power spectrum to be less than ∼ 10% (68% C. L.) for an ∼ 1 degree scale, and possibly higher for larger scales. The origin of most of the signals in the CMBR anisotropies is not Galactic. In any case, the subtraction of the galaxy is not accurate enough to allow a “precision Cosmology”; other sources of contamination (extragalactic, solar system) are also present.     

On the X-ray time-lags in the black hole candidates

It is shown that the energy dependence of the time-lags in Cygnus X-1 excludes any significant contribution of the standard reflected component to the observed lags. The conclusion is valid in the     frequency range where time-lags have been detected with sufficient significance. In fact, the data hint that the reflected component is working in the opposite direction, reducing the lags at energies where the contribution of the reflected component is significant.
We argue that the observed logarithmic dependence of time-lags on energy could be due to the small variations of the spectral index in the frame of a very simple phenomenological model. We assume that an optically thin flow/corona, emitting a power law like spectrum, is present at a range of distances from the compact object. The slope of the locally emitted spectrum is a function of distance, with the hardest spectrum emitted in the innermost region. If perturbations with different time-scales are introduced to the accretion flow at different radii, then X-ray lags naturally appear, caused by the inward propagation of perturbations on the diffusion time-scales.     

On the kinematic deconvolution of the local neighbourhood luminosity function

A method for inverting the statistical star counts equation, including proper motions, is presented; in order to break the degeneracy in that equation, it uses the supplementary constraints required by dynamical consistency. The inversion gives access to both the kinematics and the luminosity function of each population in three régimes: the singular ellipsoid, the constant ratio Schwarzschild ellipsoid plane-parallel models and the epicyclic model. This more realistic model is tailored to account for the local neighbourhood density and velocity distribution.
The first model is fully investigated, both analytically and by means of a non-parametric inversion technique, while the second model is shown to be formally its equivalent. The effect of noise and incompleteness in apparent magnitude is investigated. The third model is investigated by a     non-parametric inversion technique where positivity of the underlying luminosity function is explicitly accounted for.
It is argued that its future application to data such as the Tycho catalogue (and in the upcoming satellite GAIA ) could lead – provided that the vertical potential and or the asymmetric drift or w _⊙ are known – to a non-parametric determination of the local neighbourhood luminosity function without any reference to stellar evolution tracks. It should also yield the proportion of stars for each kinematic component and a kinematic diagnostic to split the thin disc from the thick disc or the halo.     

Environmental dependence of dark matter halo growth – I. Halo merger rates

In an earlier paper, we quantified the mean merger rate of dark matter haloes as a function of redshift z , descendant halo mass M ₀, and progenitor halo mass ratio ξ using the Millennium simulation of the Λ cold dark matter cosmology. Here, we broaden that study and investigate the dependence of the merger rate of haloes on their surrounding environment. A number of local mass overdensity variables, both including and excluding the halo mass itself, are tested as measures of a halo's environment. The simple functional dependence on   z , M ₀  , and ξ of the merger rate found in our earlier work, is largely preserved in different environments, but we find that the overall amplitude of the merger rate has a strong positive correlation with the environmental densities. For galaxy-mass haloes, we find mergers to occur ∼2.5 times more frequently in the densest regions than in voids at both   z = 0  and higher redshifts. Higher mass haloes show similar trends. We present a fitting form for this environmental dependence that is a function of both mass and local density and valid out to   z = 2  . The amplitude of the progenitor (or conditional) mass function shows a similar correlation with local overdensity, suggesting that the extended Press–Schechter model for halo growth needs to be modified to incorporate environmental effects.     

A revised flux vector for the collisional evolution of Keplerian orbits

In statistical Keplerian systems the disordered component of collisionally induced motion of matter introduces new terms into the flux vector. This contribution, which is calculated from a transport equation, tends to reduce the density gradient and causes the expansion which is observed in computer simulations of collisional systems. A quantitative comparison with Trulsen's (1972) simulations confirms the revised expression of the flux vector.     

Effect of nonextensive electron and ion on dust acoustic rogue waves in dusty plasma of opposite polarity

Propagation of nonlinear dust-acoustic waves in a magnetized collisionless plasma having positively, negatively charged dust grains and nonextensive distributed electrons and ions has been investigated. A reductive perturbation method is used to obtain a nonlinear Korteweg-de Vries (KdV) equation describing the model. The dynamics of the modulational instability gives rise to the formation of rogue waves that is described by a nonlinear Schrödinger equation. The dependence of rogue waves profiles on positive and negative charged dust cyclotron frequencies, nonextensive parameters of electrons and ions is investigated numerically. The result of the present investigation may be applicable to some plasma environments, such as cometary tails and upper mesosphere.     

Research on magnetic fields of galaxies using RZ-model

The generation of magnetic fields of galaxies is usually described by the dynamo mechanism.This process is characterized by the Steenbeck-Krause-Radler equation,which is the result of averaging the magnetohydrodynamics equations by distances which are associated with the size of turbulent cells in the interstellar medium.This equation is quite difficult to solve both from an analytical and numerical point of view.For galaxies,the no-z approximation is widely used.It describes the magnetic fields in thin discs.For such objects,where it is important to study the vertical structure of the field,it is not very applicable,so it is quite useful to adopt the RZ-model,which takes into account the dependence of the distance from the equatorial plane.During our research we have obtained the critical values of the dynamo number for galaxies with large half-thickness.We have also described typical z-structure for the magnetic field.Moreover,we have demonstrated that it is possible to generate dipolar magnetic fields.