Alfvén waves in the solar atmosphere

We examine the propagation of Alfvén waves in the solar atmosphere. The principal theoretical virtues of this work are: (i) The full wave equation is solved without recourse to the small-wavelength eikonal approximation (ii) The background solar atmosphere is realistic, consisting of an HSRA/VAL representation of the photosphere and chromosphere, a 200 km thick transition region, a model for the upper transition region below a coronal hole (provided by R. Munro), and the Munro-Jackson model of a polar coronal hole. The principal results are:

1. If the wave source is taken to be near the top of the convection zone, where n _H = 5.2 × 10¹⁶ cm^?3, and if B _⊙ = 10.5 G, then the wave Poynting flux exhibits a series of strong resonant peaks at periods downwards from 1.6 hr. The resonant frequencies are in the ratios of the zeroes of J ₀, but depend on B _⊙, and on the density and scale height at the wave source. The longest period peaks may be the most important, because they are nearest to the supergranular periods and to the observed periods near 1 AU, and because they are the broadest in frequency.
2. The Poynting flux in the resonant peaks can be large enough, i.e. P _⊙ ≈ 10⁴–10⁵ erg cm^?2s^?1, to strongly affect the solar wind.
3. ¦δv¦ and ¦δB¦ also display resonant peaks.
4. In the chromosphere and low corona, ¦δv ≈ 7–25 kms^?1 and ¦δB¦ ≈0.3–1.0 G if P _⊙≈10⁴-10⁵ erg cm^?2s^?1.
5. The dependences of ¦δv¦ and ¦δB¦ on height are reduced by finite wavelength effects, except near the wave source where they are enhanced.
6. Near the base, ¦δB¦ ≈ 350–1200 G if P _⊙ ～- 10⁴–10⁵. This means that nonlinear effects may be important, and that some density and vertical velocity fluctuations may be associated with the Alfvén waves.
7. Below the low corona most wave energy is kinetic, except near the base where it becomes mostly magnetic at the resonances.
8. ?₀ < δv ² > v _A or < δB ² > v _A/4π are not good estimators of the energy flux.
9. The Alfvén wave pressure tensor will be important in the transition region only if the magnetic field diverges rapidly. But the Alfvén wave pressure can be important in the coronal hole.

     

Energy changes of cosmic rays in the interplanetary region

A clarification and discussion of the energy changes experienced by cosmic rays in the interplanetary region is presented. It is shown that the mean time rate of change of momentum of cosmic rays reckoned for a fixed volume in a reference frame fixed in the solar system is 〈p〉 =p V·G/3 (p=momentum,V is the solar wind velocity andG=cosmic-ray density gradient). This result is obtained in three ways:

1. by a rearrangement and reinterpretation of the cosmic-ray continuity equation;
2. by using a scattering analysis based on that of Gleeson and Axford (1967);
3. by using a special scattering model in which cosmic-rays are trapped in ‘magnetic boxes’ moving with the solar wind.

The third method also gives the rate of change of momentum of particles within a moving ‘magnetic box’ as 〈p〉_ad = ?p ?·V/3, which is the adiabatic deceleration rate of Parker (1965). We conclude that ‘turnaround’ energy change effects previously considered separately are already included in the equation of transport for cosmic rays.     

Problems with non-thermal models for the narrow line gamma rays reported from SS 433

The jet/grain model proposed by Ramatyet al. (1984, hereafter abbreviated as RKL) for production of the narrow gamma-ray lines reported from SS433 is examined and shown to be untenable on numerous grounds. Most importantly:

1. The huge Coulomb collisional losses (W _c?2×10⁴¹ erg s^?1) from the jet, which would necessarily accompany non-thermal production of the gamma rays, demands a jet acceleration/collimation process acting over a very long range and with a power at least 10² times the Eddington limit for any stellar object.
2. There is a collisional thick target limit (irrespective of jet mass) to the gamma ray yield per interstellar proton. Consequently, the gamma-ray data demand an improbably high interstellar density (?10⁹ cm^?3).
3. For the grains to be kept cool enough (?3000 K) to survive the heating rateW _c either by radiation or jet expansion would demand a ‘jet’ wider than its length and so inconsistent with narrow lines. In the case of radiative cooling, the resultant IR flux would exceed the observed values by a factor ?10⁴.
4. Light scattered on the jet grain mass required would be highly polarized, contrary to observations, unless the jet was optically thick to grains, again precluding their radiative cooling.
5. To avoid unacceptable precessional broadening of the gamma-ray lines demands an emitting jet length ?0.5 days atv=0.26c. This increases the necessary mass loss rate by a factor ?10 over the values obtained by RKL who assumed a 4-day ‘flare’.
6. The model also predicts rest energy gamma-ray lines which are not observed.

     

Some properties of finite energy constant-α force-free magnetic fields in a half-space

Some useful properties of a finite energy, constant-α, force-free magnetic field B _α occupying a half-space D are presented. In particular:

1. Fourier and Green representations of B _α are obtained and used to derive conditions for the existence and uniqueness of a B _α having a given normal component B _z on the boundary ?D.
2. The asymptotic behaviour of B _α at infinity as well as stability results against changes in the boundary condition on ?D and in the value of α are established.
3. The energy of B _α is shown to be smaller than the energy of the open field having the same B _z on ?D, thus confirming an earlier conjecture (Aly, 1984).
4. B _α is proved to not be a Taylor-Heyvaerts-Priest state, in spite of the fact that its relative helicity H is finite and that it is the only solution of the Lagrange-Euler equation associated with the problem of minimizing the energy among all the fields having the same value of H and the same B _z on ?D.

     

Effect of perturbed potentials on the stability of libration points in the restricted problem

The location and the stability in the linear sense of the libration points in the restricted problem have been studied when there are perturbations in the potentials between the bodies. It is seen that if the perturbing functions satisfy certain conditions, there are five libration points, two triangular and three collinear. It is further observed that the collinear points are unstable and for the triangular points, the range of stability increases or decreases depending upon whetherP> or <0 wherep depends upon the perturbing functions. The theory is verified in the following four cases:

1. There are no perturbations in the potentials (classical problem).
2. Only the bigger primary is an oblate spheroid whose axis of symmetry is perpendicular to the plane of relative motion (circular) of the primaries.
3. Both the primaries are oblate spheroids whose axes of symmetry are perpendicular to the plane of relative motion (circular) of the primaries.
4. The primaries are spherical in shape and the bigger is a source of radiation.

     

On the velocity of jets powering hot spots similar to those in Cygnus A

Hot spots similar to those in the radio galaxy Cygnus A can be explained by the strong shock produced by a supersonic but classical jet \(\left( {u_{jet}< c/\sqrt 3 } \right)\) . The high integrated radio luminosity (L?2×10⁴⁴ erg s^?1) and the strength of mean magnetic field (B?2×10^?4 G) suggest the hot spots are the downstream flow of a very strong shock which generates the ultrarelativistic electrons of energy ?≥20 MeV. The fully-developed subsonic turbulence amplifies the magnetic field of the jet up to 1.6×10^?4 G by the dynamo effect. If we assume that the post-shock pressure is dominated by relativistic particles, the ratio between the magnetic energy density to the energy density in relativistic particles is found to be ?2×10^?2, showing that the generally accepted hypothesis of equipartition is not valid for hot spots. The current analysis allows the determination of physical parameters inside hot spots. It is found that:

1. The velocity of the upstream flow in the frame of reference of the shock isu ₁?0.2c. Radio observations indicate that the velocity of separation of hot spots isu _sep?0.05c, so that the velocity of the jet isu _jet=u ₁+u _sep?0.25c.
2. The density of the thermal electrons inside the hot spot isn ₂?5×10^?3 e ^? cm^?3 and the mass ejected per year to power the hot spot is ?4M ₀yr^?1.
3. The relativistic electron density is less than 20% of the thermal electron density inside the hot spot and the spectrum is a power law which continues to energies as low as 30 MeV.
4. The energy density of relativistic protons is lower than the energy density of relativistic electrons unlike the situation for cosmic rays in the Galaxy.

     

Study of physical properties of spectroscopic binary stars

The main results of a study of a catalogue of physical parameters of 1041 spectroscopic binaries are presented. The distribution of spectroscopic binaries over all main parametersM ₁, a, e, M₁/M₂, P, and certain dependencies between some of them have been found.

1. It appears that among bright (m _v?3 ^m –5 ^m ) stars withM?1M _⊙, about 40% are apparently spectroscopic binaries with comparable masses of components.
2. The majority of spectroscopic binaries with the ratio of the large semiaxis of the orbit to the radius of the primarya/R ₁?20, have eccentricities close to zero. This is probably a consequence of the tidal circularization of orbits of close binaries by viscous friction.
3. The discovery of duplicity of double-line spectroscopic binaries is possible only if the semiamplitude of radial velocityK ₁ is almost 10 times higher than the semiamplitude of the radial velocity of a single-line spectroscopic binary of the same mass.
4. Double-line spectroscopic binaries witha/R _⊙?6(M ₁/M _⊙)^1/3,M ₁≈M ₂?1.5M _⊙ are almost almost absent, and the number of stars witha/R _⊙?6(M ₁/M _⊙)^1/3,M ₁≈1.5M _⊙ is relatively low.
5. The distribution of unevolved SB stars over the large semiaxis may be described by the expression d(N _d/N_t)≈0.2 d loga for 6(M ₁/M _⊙)^1/3?a/R _⊙?100.
6. The intial mass-function for primaries of spectroscopic binaries is the same Salpeter function dN _d≈M ₁ ^?2.35 dM ₁ for 1?M ₁/M _⊙?30.
7. It is possible to explain the observed ratio of the number of single-line spectroscopic binaries to the number of double-line binaries if one assumes that the average initial mass ratio is close to 1 and that the mass of the postmass-exchange remnant of the primary exceeds the theoretical one and/or that half of the angular momentum of the system is lost during mass-exchange.
8. The above-mentioned distributions ofM ₁ anda and assumptions on the mass of remnant and/or momentum loss also allow us to explain the observed shapes of dN/dM, dN/dq, and dN/da distributions after some selection effects are taken into account.

     

A regularization of the three-body problem

Letr ₁,r ₂,r ₃ be arbitrary coordinates of the non-zero interacting mass-pointsm ₁,m ₂,m ₃ and define the distancesR ₁=|r ₁?r ₃|,R ₂=|r ₂?r ₃|,R=|r ₁?r ₂|. An eight-dimensional regularization of the general three-body problem is given which is based on Kustaanheimo-Stiefel regularization of a single binary and possesses the properties:

1. The equations of motion are regular for the two-body collisionsR ₁→0 orR ₂→0.
2. Provided thatR?R ₁ orR?R ₂, the equations of motion are numerically well behaved for close triple encounters.

Although the requirementR? min (R ₁,R ₂) may involve occasional transformations to physical variables in order to re-label the particles, all integrations are performed in regularized variables. Numerical comparisons with the standard Kustaanheimo-Stiefel regularization show that the new method gives improved accuracy per integration step at no extra computing time for a variety of examples. In addition, time reversal tests indicate that critical triple encounters may now be studied with confidence. The Hamiltonian formulation has been generalized to include the case of perturbed three-body motions and it is anticipated that this procedure will lead to further improvements ofN-body calculations.     

Exotic ion H 3 ++ in strong magnetic fields

1. The exotic system H ₃ ⁺⁺ (which does not exist without magnetic field) exists in strong magnetic fields:
   1. In triangular configuration for B≈10⁸–10¹¹?G (under specific external conditions)
   2. In linear configuration for B>10¹⁰?G
2. In the linear configuration the positive z-parity states 1σ _g , 1π _u , 1δ _g are bound states
3. In the linear configuration the negative z-parity states 1σ _u , 1π _g , 1δ _u are repulsive states
4. The H ₃ ⁺⁺ molecular ion is the most bound one-electron system made from protons at B>3×10¹³?G

Possible application: The H ₃ ⁺⁺ molecular ion may appear as a component of a neutron star atmosphere under a strong surface magnetic field B=10¹²–10¹³?G.     

On the spin distributions of ΛCDM haloes

We used merger trees realizations, predicted by the extended Press-Schechter theory, in order to study the growth of angular momentum of dark matter haloes. Our results showed that:

1. The spin parameter λ′ resulting from the above method, is an increasing function of the present day mass of the halo. The mean value of λ′ varies from 0.0343 to 0.0484 for haloes with present day masses in the range of 10⁹h^?1 M _⊙ to 10¹⁴h^?1 M _⊙.
2. The distribution of λ′ is close to a log-normal, but, as it is already found in the results of N-body simulations, the match is not satisfactory at the tails of the distribution. A new analytical formula that approximates the results much more satisfactorily is presented.
3. The distribution of the values of λ′ depends only weakly on the redshift.
4. The spin parameter of an halo depends on the number of recent major mergers. Specifically the spin parameter is an increasing function of this number.

     

Spectral analysis of four quiescent prominences observed at the Peruvian eclipse

Two-dimensional distributions of kinetic temperature, density and turbulent velocity are obtained for four quiescent prominences observed at the Peruvian eclipse of 12 November, 1966.

1. The kinetic temperature derived from line widths is around 6000–7000 K in the central part of prominences and rises to 12000K in both edges and possibly in the top of prominences.
2. The turbulent velocity shows a similar tendency, being 7–9 km/sec in the central part and ≈ 20 km/sec in the outer part. The turbulent velocity also increases slowly towards higher heights in the prominence.
3. The electron density derived both from the Stark effect and the intensity ratio of the continuous spectra turns out to be about 10^10.2–10^10.6 cm^?3 in the central portion of two prominences.
4. From the width and the intensity, neutral helium lines are shown to originate in the same region as hydrogen and metallic lines where the kinetic temperature goes down to 6000 K. This indicates that neutral helium is emitted after the ionization due to UV radiation from the corona and the transition region.

     

A varying cosmological term as a link between cosmology and microphysics

The Weinberg relation (which connects the Hubble constantH to the mass of a typical elementary particle) is an empirical relation hitherto unexplained. I suggest an explanation based on the Zel'dovich energy tensor of vacuum in a Robertson-Walker universe with constant deceleration parameter,q = const. This model leads to

1. the Weinberg relation,
2. a varying cosmological term Λ scaling asH ²,
3. a varying gravitational constantG scaling asH,
4. a matter creation process throughout the universe at the rate 10^?47 g s^?1 cm³,
5. a deceleration parameter in the range -1 to 1/2, which allows a horizon-free universe and makes the lawG/H = constant, consistent with the Viking lander data on the orbit of planet Mars.

     

Drift theory of charged particles in electric and magnetic fields

In this paper we review the drift theory of charged particles in electric and magnetic fields. No new physical interpretations are added to this classical topic, but through an alternative, simplified derivation of the guiding centre velocity, several complexities are eliminated and possible misconceptions of the theory are clarified. It is shown that:

1. The curvature/gradient drift velocity in the magnetic field, averaged over a particle distribution function is to lowest order in the direction of?×B/B ², while the average particle velocity is in the direction ofB×? P withP the scalar particle pressure.
2. These drift directions are correct for first-order expansions of the particle distribution function, and only second-order or higher expansions change these directions.
3. The?×B/B ² drift, which is the standard gradient plus curvature drift, and which is usually considered as a ‘single particle’ drift, need not be ‘reconciled’ with theB×? P, or ‘macroscopic, collective’ drift, as is often asserted in the literature. They are in fact related per definition and we show how.
4. When viewed in fixed momentum intervals (p,p+dp), the so-called Compton-Getting factor enters into the electric field (E×B)/B ² drift term.
5. The results are independent of the scale length of variation ofE andB, in contrast to existing drift theory. We discuss the implications of this result for three important cases.

     

On the theoretical possibility of the libration cloud

We show within the framework of the restricted three body problem that

1. Only in the immediate neighbourhood of the Lagrangian pointsL ₄ andL ₅ the distribution of a cloud of particles tends to become uniform under the influence of random stochastic perturbations. No consequences can be derived from this fact for a tendency of dispersion of clouds librating at arbitrary distances around the Lagrangian points.
2. From an elementary inspection of the Jacobi integral we cannot conclude that mutual completely inelastic collisions tend to drive the particles away from the vicinity of the libration points.

Finally the motion of a particle within the libration cloud, approximated as a resisting medium, is briefly examined.     

Effects of temperature and velocity gradients on doppler widths

We have investigated how the gradients of temperature and expansion velocities will change the emergent profiles from an extended medium in spherical symmetry. Variation of the source function and expansion velocities are assumed. The following variations of temperature are employed:

1. T(r) ; T₀ (isothermal case)
2. T(r) ; T₀(r/r₀)^1/2
3. T(r) ; T₀(r/r₀)^-1
4. T(r) ; T₀(r/r₀)^-2
5. T(r) ; T₀(r/r₀)^-3

The profiles calculated present an interesting feature of broadening.     

Radio spectrum of two active regions

The results of the total solar eclipse of November 12, 1966, observed at 8 different wave-lengths between 3 and 21 cm, are studied and the spectrum of two active regions present on the disk is deduced. It is shown that the observed increase of the flux of the most intense source in the range 3–10 cm is due to geometrical effects. Neglecting the influence of the magnetic field, the following quantities are deduced.

1. the mean and central temperature of the coronal condensation.
2. the ∫ _corona N ²dh (N = electron density).

Both these quantities are in good agreement with optical observations.     

Orbital eccentricity study for the spherical component of our galaxy

An analysis of the data concerning high-velocity stars from Eggen's catalogue aimed at a determination of the approximate slope of the mass function for the spherical component of our Galaxy, and at estimating the local circular velocity, as well as the local rotation velocity, as by-products, has been performed. Our conclusions are that:

1. A linear dependence of the mass on the radius is very likely;
2. the value of the limiting radius is most likely equal to (40±10) kpc;
3. the two local velocities are approximately equal to each other, being both equal to (230±30) km s^?1;
4. the local escape velocity appears to be most likely equal to (520±30) km s^?1;
5. the total mass of a corona, obtained in this way, is (5±1)×10¹¹ M _⊙.

     

Problems of matter-antimatter boundary layers

This paper outlines the problems of the quasi-steady matter-antimatter boundary layers discussed in Klein-Alfvén's cosmological theory, and a crude model of the corresponding ambiplasma balance is presented:

1. At interstellar particle densities, no well-defined boundary layer can exist in presence of neutral gas, nor can such a layer be sustained in an unmagnetized fully ionized ambiplasma.
2. Within the limits of applicability of the present model, sharply defined boundary layers are under certain conditions found to exist in a magnetized ambiplasma. Thus, at beta values less than unity, a steep pressure drop of the low-energy components of matter and antimatter can be balanced by a magnetic field and the electric currents in the ambiplasma.
3. The boundary layer thickness is of the order of 2x ₀?10/BT ₀ ^1/4 metres, whereB is the magnetic field strength in MKS units andT ₀ the characteristic temperature of the low-energy components in the layer.

     

On the Hβ, [OIII] lines and continuum variability character in the nucleus of Seyfert galaxy NGC 1275

The radiation fluxes of the NGC 1275 galaxy central region are being observed on the 1.25-m telescope, using a scanning spectrophotometer with the entrance aperture 10″ in three Δλ=80 Å spectral regions: Hβ, 4959+5007 Å [OIII] and continuum. There were 35 nights of observations during 1982–1987. With the time resolution of half an hour 379 measurements were obtained in each spectral region. The analysis of these results shows:

1. The standard deviations of measurements in each spectral region 2–3 times exceed the errors of observations.
2. The radiation flux distribution resembles to normal one only for Hβ line.
3. Two-humps forms of continuum flux distribution curve is like that of radio emission in 8 mm and 2.6 cm wavelengths.
4. Various forms of fluxes distribution curves of Hβ and [OIII] lines permit us to suppose that the location of these lines emission regions near the sources of excitation are different.

     

Keynote address: Outstanding problems in solar physics

Celebrating the diamond jubilee of the Physics Research Laboratory (PRL) in Ahmedabad, India, we look back over the last six decades in solar physics and contemplate on the ten outstanding problems (or research foci) in solar physics:

1. The solar neutrino problem
2. Structure of the solar interior (helioseismology)
3. The solar magnetic field (dynamo, solar cycle, corona)
4. Hydrodynamics of coronal loops
5. MHD oscillations and waves (coronal seismology)
6. The coronal heating problem
7. Self-organized criticality (from nanoflares to giant flares)
8. Magnetic reconnection processes
9. Particle acceleration processes
10. Coronal mass ejections and coronal dimming

The first two problems have been largely solved recently, while the other eight selected problems are still pending a final solution, and thus remain persistent Challenges for Solar Cycle 24, the theme of this jubilee conference.