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.     

Systematics of Dwarf Novae

From about 30 Dwarf Novae with the best determined distances the following relationships are found.

1. a tight correlation between absolute magnitude at maximum light, M_v(max), and orbital period, P.
2. a correlation between M_v(min) and P showing wide scatter.
3. a correlation between M_v(mean), the mean absolute magnitude averaged over normal outbursts, and P, again with wide scatter. The scatter is shown to correlate strongly with ouburst timescale T_n.
4. a strong correlation between range, M_v(min)-M_v(max), and T_n (the Kukarkin-Parenago relationship).
5. a strong correlation between range M_v(mean)-M_v(max), and both T_n and P.

This final correlation is interpreted in terms of the disc instability model of dwarf novae and successfully predicts the observed width of outburst versus P relationship.     

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.

     

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.     

Observation of the impulsive phase of a simple flare

We present a broad range of complementary observations of the onset and impulsive phase of a fairly large (1B, M1.2) but simple two-ribbon flare. The observations consist of hard X-ray flux measured by the SMM HXRBS, high-sensitivity measurements of microwave flux at 22 GHz from Itapetinga Radio Observatory, sequences of spectroheliograms in UV emission lines from Ov (T ≈ 2 × 10⁵ K) and Fexxi (T ≈ 1 × 10⁷ K) from the SMM UVSP, Hα and Hei D₃ cine-filtergrams from Big Bear Solar Observatory, and a magnetogram of the flare region from the MSFC Solar Observatory. From these data we conclude:

1. The overall magnetic field configuration in which the flare occurred was a fairly simple, closed arch containing nonpotential substructure.
2. The flare occurred spontaneously within the arch; it was not triggered by emerging magnetic flux.
3. The impulsive energy release occurred in two major spikes. The second spike took place within the flare arch heated in the first spike, but was concentrated on a different subset of field lines. The ratio of Ov emission to hard X-ray emission decreased by at least a factor of 2 from the first spike to the second, probably because the plasma density in the flare arch had increased by chromospheric evaporation.
4. The impulsive energy release most likely occurred in the upper part of the arch; it had three immediate products:

1. An increase in the plasma pressure throughout the flare arch of at least a factor of 10. This is required because the Fexxi emission was confined to the feet of the flare arch for at least the first minute of the impulsive phase.
2. Nonthermal energetic (～ 25 keV) electrons which impacted the feet of the arch to produce the hard X-ray burst and impulsive brightening in Ov and D₃. The evidence for this is the simultaneity, within ± 2 s, of the peak Ov and hard X-ray emissions.
3. Another population of high-energy (～100keV) electrons (decoupled from the population that produced the hard X-rays) that produced the impulsive microwave emission at 22 GHz. This conclusion is drawn because the microwave peak was 6 ± 3 s later than the hard X-ray peak.

     

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.

     

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.

     

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.

     

The observational discrimination of Friedmann-Lemaître models

It is the purpose of this paper to illustrate the interrelation between the problems of the ‘missing mass’, the galactic age and the cosmological constant A (or its equivalent quantum vacuum densityρ _v ). The inflationary picture of the early universe predicts that our present universe should have a very nearlyEuclidean metric. If we accept this concept, one would have to discriminate between two rather extreme Euclidean cosmological models:

1. The standard model with ∧=0 and a densityρ _c = 3H ₀ ² 8πG. There are difficulties ifH ₀≥5- km s^?1 Mpc^?1 and the galactic aget ₀≥14×10⁹ years.
2. The Euclidean Friedmann-Lemaître models with ∧>0, i.e., Δgt;0, i.e.,ρ _v =ρ _c ?ρ _o , whereρ _o is the present matter density, including the nonrelativistic dark matter. Hereρ _v ‘competes’ with the missing mass.

Measurements of apparent diameters of galaxies up to redshifts of 2 will permit one to discriminate between the models provided that size evolution of galaxies can be determined or neglected (see Figure 3).     

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.     

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.

     

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.

     

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 _⊙.

     

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.

     

Analysis of the properties of clusters of galaxies in the region of the Ursa Major supercluster

We analyze the properties of the clusters of galaxies in the region of the Ursa Major (UMa) supercluster using observational data from SDSS and 2MASS catalogs. The region studied includes a supercluster (with a galaxy and cluster overdensity of 3 and 15, respectively) and field clusters inside the 150-Mpc diameter surrounding region. The total dynamical mass of 10 clusters of galaxies in UMa is equal to 2.25 × 10¹⁵ M _⊙, and the mass of 11 clusters of galaxies in the UMa neighborhood is equal to 1.70 × 10¹⁵ M _⊙. The fraction of early-type galaxies brighter than M _K ^* + 1 in the virialized regions of clusters is, on the average, equal to 70%, and it is virtually independent on the mass of the cluster. The fraction of these galaxies and their average photometric parameters are almost the same both for UMa clusters and for the clusters located in its surroundings. Parameters of the clusters of galaxies, such as infrared luminosities up to a fixed magnitude, the mass-to-luminosity ratio, and the number of galaxies have almost the same correlations with the cluster mass as in other samples of galaxies clusters. However, the scatter of these parameters for UMa member clusters is twice smaller than the corresponding scatter for field clusters, possibly, due to the common origin of UMa clusters and synchronized dynamical evolution of clusters in the supercluster.     

Flares and changing magnetic fields

An observational study of maps of the longitudinal component of the photospheric fields in flaring active regions leads to the following conclusions:

1. The broad-wing Hα kernels characteristic of the impulsive phase of flares occur within 10″ of neutral lines encircling features of isolated magnetic polarity (‘satellite sunspots’).
2. Photospheric field changes intimately associated with several importance 1 flares and one importance 2B flare are confined to satellite sunspots, which are small (10″ diam). They often correspond to spot pores in white-light photographs.
3. The field at these features appears to strengthen in the half hour just before the flares. During the flares the growth is reversed, the field drops and then recovers to its previous level.
4. The magnetic flux through flare-associated features changes by about 4 × 10¹⁹ Mx in a day. The features are the same as the ‘Structures Magnétiques Evolutives’ of Martres et al. (1968a).
5. An upper limit of 10²¹ Mx is set for the total flux change through McMath Regions 10381 and 10385 as the result of the 2B flare of 24 October, 1969.
6. Large spots in the regions investigated did not evince flux changes or large proper motions at flare time.
7. The results are taken to imply that the initial instability of a flare occurs at a neutral point, but the magnetic energy lost cannot yet be related to the total energy of the subsequent flare.
8. No unusual velocities are observed in the photosphere at flare time.

     

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.

     

On pulsar-driven isothermal blast-wave models of supernova remanants

We investigate the ‘equilibrium’ and stability of spherically-symmetric self-similar isothermal blast waves with a continuous post-shock flow velocity expanding into medium whose density varies asr ^?ω ahead of the blast wave, and which are powered by a central source (a pulsar) whose power output varies with time ast ^ω?3. We show that:

1. for ω<0, no physically acceptable self-similar solution exists;
2. for ω>3, no solution exists since the mass swept up by the blast wave is infinite;
3. ? must exceed zero in order that the blast wave expand with time, but ?<2 in order that the central source injects a finite total energy into the blast wave;
4. for 3>ω_min(?)>ω>ω_max(?)>0, where $$\begin{gathered} \omega _{\min } (\varphi ){\text{ }} = {\text{ }}2[5{\text{ }} - {\text{ }}\varphi {\text{ }} + {\text{ }}(10{\text{ }} + {\text{ 4}}\varphi {\text{ }} - {\text{ 2}}\varphi ^2 )^{1/2} ]^2 [2{\text{ }} + {\text{ (10 }} + {\text{ 4}}\varphi {\text{ }} - {\text{ 2}}\varphi ^2 {\text{)}}^{{\text{1/2}}} ]^{ - 2} , \hfill \\ \omega _{\max } (\varphi ){\text{ }} = {\text{ }}2[5{\text{ }} - {\text{ }}\varphi {\text{ }} - {\text{ }}(10{\text{ }} + {\text{ 4}}\varphi {\text{ }} - {\text{ 2}}\varphi ^2 )^{1/2} ]^2 [2{\text{ }} - {\text{ (10 }} + {\text{ 4}}\varphi {\text{ }} - {\text{ 2}}\varphi ^2 {\text{)}}^{{\text{1/2}}} ]^{ - 2} , \hfill \\ \end{gathered} $$ two critical points exist in the flow velocity versus position plane. The physically acceptable solution must pass through the origin with zero flow speed and through the blast wave. It must also pass throughboth critical points if \(\varphi > \tfrac{5}{3}\) , while if \(\varphi< \tfrac{5}{3}\) it must by-pass both critical points. It is shown that such a solution exists but a proper connection at the lower critical point (for ?>5/3) (through whichall solutions pass with thesame slope) has not been established;
5. for 3>ω>ω_min(?) it is shown that the two critical points of (iv) disappear. However a new pair of critical points form. The physically acceptable solution passing with zero flow velocity through the origin and also passing through the blast wave mustby-pass both of the new critical points. It is shown that the solution does indeed do so;
6. for 3>ω_min(?)>ω_max(?)>ω it is shown that the dependence of the self-similar solution on either ω or ? is non-analytic and therefore, inferences drawn from any solutions obtained in ω>ω_max(?) (where the dependence of the solutionis analytic on ω and ?) are not valid when carried over into the domain 3>ω_min(?)>ω_max(?)>ω;
7. all of the physically acceptable self-similar solutions obtained in 3>ω>0 are unstable to short wavelength, small amplitude but nonself-similar radial velocity perturbations near the origin, with a growth which is a power law in time;
8. the physical self-similar solutions are globally unstable in a fully nonlinear sense to radial time-dependent flow patterns. In the limit of long times, the nonlinear growth is a power law in time for 5<ω+2?, logarithmic in time for 5>ω+2?, and the square of the logarithm in time for 5=ω+2?.

The results of (vii) and (viii) imply that the memory of the system to initial and boundary values does not decay as time progresses and so the system does not tend to a self-similar form. These results strongly suggest that the evolution of supernova remnants is not according to the self-similar form.     

A statistical study of the spiral spots on the solar disc

After adding the data observed in the years from 1979 to 1982 to those obtained earlier (Ding et al., 1981), we re-examine the previous results and conclude:

1. The longitudinal distribution of spiral spots on the solar disc is generally the same as that of sunspot groups with areas of S _p ≥ 400, but their active longitudes seem to be more concentrated.
2. The distribution of spiral patterns in the southern and northern hemispheres shows that the differential rotation may be a fundamental solar dynamo for the formation of the spiral spots.
3. The statistical directions of the emerging twisted magnetic vectors in the active regions in the southern and northern hemispheres are synchronously inverse with a period of about two years. This period seems to be detected in other solar observations.

     

Source characteristics of main and post-burst-increase phases of solar bursts at 17 GHz

Twenty four solar bursts of peak fluxes above 50 sfu are analyzed which were observed with the 17 GHz interferometer at Nobeyama during the period from 1978 September to 1979 December. Source characteristics and their temporal evolutions are investigated on a statistical basis with high time resolutions up to 0.8 s. Use of a model-fitting technique recently developed by Kosugi (1982) is made to derive both the position of centroid and size (～ FWHM) of burst source with an uncertainty of a few arc sec. The results of this study are the following:

1. Two different phases in the burst, that is to say, the main phase and the post-burst-increase (PBI) phase, are distinguished clearly not only by the morphological difference of flux time profile, but also by the differences of brightness temperature (10⁷-?10⁹ K vs 10⁵–10⁷ K), circular polarization degree (0–50% vs 0–10%), and size (?5–25″ vs 10–70″). There is no definite correlation between the peak fluxes in the two phases.
2. The majority of the selected bursts (21 of 24) show in the main phase source characteristics of the impulsive burst. The total flux varies rapidly (characteristic time scale defined by FWHM ? 100 s), often associated with the rapid shift of position and the rapid change of polarization degree. The source height of the impulsive source is lower than that of the PBI source. On the other hand, the type IV_μ source, seen in three events, shows a gradual variation and the source ascends to a height of ～ 40 000 km above the photosphere.
3. In the PBI phase, the expansion and ascension of the source occur in general (21 of 23 for the former and 12 of 15 for the latter). The velocities of both the movements are of the order of 5 km s^?1.