Sur un formalisme explicitant la loi des distances des planetes

The well-known Titius-Bode law (T-B) giving distances of planets from the Sun was improved by Basano and Hughes (1979) who found: $$a_n = 0.285 \times 1.523^n ;$$ a _n being the semi-major axis expressed in astronomical units, of then-th planet. The integern is equal to 1 for Mercury, 2 for Venus etc. The new law (B-H) is more natural than the (T-B) one, because the valuen=?∞ for Mercury is avoided. Furthermore, it accounts for distances of all planets, including Neptune and Pluto. It is striking to note that this law:

1. does not depend on physical parameters of planets (mass, density, temperature, spin, number of satellites and their nature etc.).
2. shows integers suggesting an unknown, obscure wave process in the formation of the solar system.

In this paper, we try to find a formalism accounting for the B-H law. It is based on the turbulence, assumed to be responsible of accretion of matter within the primeval nebula. We consider the function $$\psi ^2 (r,t) = |u^2 (r,t) - u_0^2 |$$ , whereu ²(r, t) stands for the turbulence, i.e., the mean-square deviation velocities of particles at the pointr and the timet; andu ₀ ² is the value of turbulence for which the accretion process of matter is optimum. It is obvious that Ψ²(r _n,t₀) = 0 forr _n=0.285×1.523 ⁿ at the birth timet ₀ of proto-planets. Under these conditions, it is easily found that $$\psi ^2 (r,t_0 ) = \frac{{A^2 }}{r}\sin ^2 [\alpha log r - \Phi (t_0 )]$$ With α=7.47 and Φ(t ₀)=217.24 in the CGS system, the above function accounts for the B-H law. Another approach of the problem is made by considering fluctuations of the potentialU(r, t) and of the density of matter ρ(r, t). For very small fluctuations, it may be written down the Poisson equation $$\Delta \tilde U(r,t_0 ) + 4\pi G\tilde \rho (r,t_0 ) = 0$$ , withU(r, t)=U ₀(r)+?(r, t ₀ ) and \(\tilde \rho (r,t_0 )\) . It suffices to postulate \(\tilde \rho (r,t_0 ) = k[\tilde U(r,t_0 )/r^2 ](k = cte)\) for finding the solution $$\tilde U(r,t_0 ) = \frac{{cte}}{{r^{1/2} }}\cos [a\log r - \zeta (t_0 )]$$ . Fora=14.94 and ζ(t ₀)=434.48 in CGS system, the successive maxima of ?(r,t ₀) account again for the B-H law. In the last approach we try to write Ψ(r, t) under a wave function form $$\Psi ^2 (r,t) = \frac{{A^2 }}{r}\sin ^2 \left[ {\omega \log \left( {\frac{r}{v} - t} \right)} \right].$$ It is emphasized that all calculations are made under mathematical considerations.     

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.     

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.

     

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.

     

Large scale circulation in the convection zone and solar differential rotation

In this paper we study the dependence on depth and latitude of the solar angular velocity produced by a meridian circulation in the convection zone, assuming that the main mechanism responsible for setting up and driving the circulation is the interaction of rotation with convection. We solve the first order equations (perturbation of the spherically symmetric state) in the Boussinesq approximation and in the steady state for the axissymmetric case. The interaction of convection with rotation is modelled by a convective transport coefficient k _c = k _co + ?k _c2 P ₂(cos θ) where ? is the expansion parameter, P ₂ is the 2nd Legendre polynomial and k _c2 is taken proportional to the local Taylor number and the ratio of the convective to the total fluxes. We obtain the following results for a Rayleigh number 10³ and for a Prandtl number 1:

1. A single cell circulation extending from poles to the equator and with circulation directed toward the equator at the surface. Radial velocities are of the order of 10 cm s^?1 and meridional ones of the order of 150 cm s^?1.
2. A flux difference between pole and equator at the surface of about 5 percent, the poles being hotter.
3. An angular velocity increasing inwards.
4. Angular velocity constant surfaces of spheroidal shape. The model is consistent with the fact that the interaction of convection with rotation sets up a circulation (driven by the temperature gradient) which carries angular momentum toward the equator against the viscous friction. Unfortunately also a large flux variation at the surface is obtained. Nevertheless it seems that the model has the basic requisites for correct dynamo action.

     

Hydrodynamic and turbulent motions in the galactic disk,II

From a comparative study between stellar and gas data it is seen that turbulent and hydrodynamic motions in the Galaxy are common to both types of materials:

1. Galactic clusters have sizes and intrinsic dispersions compatible with the modified form of the Kolmogorov law seen in molecular clouds: undimensional velocities σ(km s^?1)=0.54d ^0.38 (pc). This indicates that ‘typic’ clusters were born from ‘typic’ dark clouds as these of the Lynds's catalogue (diametersd<10 pc, dispersions σ<1.5 km s^?1 hydrogen densitiesn _H>200 atom cm^?3). These clouds have mass enough to form galactic clusters (1000–3000M _⊙).
2. The cluster formation is related to the supersonic range of the Kolmogorov relationship σ(d>1 pc) while the AFGKM stars are related to the subsonic range of the same relationship σ(d<0.3 pc), the intermediate transition zone is probably related to OB stars and/or trapezia.
3. The effects of the magnetic fields in the clouds are also discussed. It seems to be that in the clouds the magnetic energy does not exceed the kinetic energy (proportional toσ ²(d)) and that this determinates the freezing criteria. The hypotheses introduced here can be checked with 21 cm Zeeman splitting.
4. Low-density globular clusters are also coherent with the Kolmogorov relationship. Some hypotheses about their origin and the type of clouds where they were born are discussed. This last part of the study lets open the possibility of further studies about evolution of globular clusters.

     

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.

     

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.

     

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

     

Recent progress in the theory of Trojan asteroids

In previous publications the author has constructed a long-periodic solution of the problem of the motion of the Trojan asteroids, treated as the case of 1:1 resonance in the restricted problem of three bodies. The recent progress reported here is summarized under three headings:

1. The nature on the long-periodic family of orbits is re-examined in the light of the results of the numerical integrations carried out by Deprit and Henrard (1970). In the vicinity of the critical divisor $$D_k \equiv \omega _1 - k\omega _2 ,$$ not accessible to our solution, the family is interrupted by bifurcations and shortperiodic bridges. Parametrized by the normalized Jacobi constant α², our family may, accordingly, be defined as the intersection of admissible intervals, in the form $$L = \mathop \cap \limits_j \left\{ {\left| {\alpha - \alpha _j } \right| > \varepsilon _j } \right\};j = k,k + 1, \ldots \infty .$$ Here, {α_j(m)} is the sequence of the critical α_j corresponding to the exactj: 1 commensurability between the characteristic frequencies ω₁ and ω₂ for a given value of the mass parameterm. Inasmuch as the ‘critical’ intervals |α?α_j|<ε_j can be shown to be disjoint, it follows that, despite the clustering of the sequence {α_j} at α=1, asj→∞, the family extends into the vicinity of the separatrix α=1, which terminates the ‘tadpole’ branch of the family.
2. Our analysis of the epicyclic terms of the solution, carrying the critical divisorD _k , supports the Deprit and Henrard refutation of the E. W. Brown conjecture (1911) regarding the termination of the tadpole branch at the Lagrangian pointL ₃. However, the conjecture may be revived in a refined form. “The separatrix α=1 of the tadpole branch spirals asymptotically toward a limit cycle centered onL ₃.”
3. The periodT(α,m) of the libration in the mean synodic longitude λ in the range $$\lambda _1 \leqslant \lambda \leqslant \lambda _2$$ is given by a hyperelliptic integral. This integral is formally expanded in a power series inm and α² or \(\beta \equiv \sqrt {1 - \alpha ^2 }\) .

The large amplitude of the libration, peculiar to our solution, is made possible by the mode of the expansion of the disturbing functionR. Rather than expanding about Lagrangian pointL ₄, with the coordinatesr=1, θ=π/3, we have expandedR about the circler=1. This procedure is equivalent to analytic continuation, for it replaces the circle of convergence centered atL ₄ by an annulus |r?1|<ε with 0≤θ<2π.     

Combined study of the solar neighbourhood kinematics: Spherical harmonics and Taylor expansions

In this paper we relate two methods of analyzing the kinematic parameters of the local macroscopic motions of the Galaxy:

1. The Ogorodnikov-Milne model (OM) that consists in the three-dimensional Taylor expansion of the mean velocity field.
2. The two-dimensional spherical harmonic development of the velocity components (SH). We present the theoretical relations between the SH coefficients and the second-order OM ones for the radial velocityv _r and the galactic heliocentric components of the velocityU, V, W. Only the hypothesis of separability of the stellar density function of the sample into angular and radial parts is needed. Also, we apply them to 4732 A-M stars included in the Figueras (1986) sample.

     

Small-scale fluctuations of relic radiation

Perturbations of the matter density in a homogeneous and isotropic cosmological model which leads to the formation of galaxies should, at later stages of evolution, cause spatial fluctuations of relic radiation. Silk assumed that an adiabatic connection existed between the density perturbations at the moment of recombination of the initial plasma and fluctuations of the observed temperature of radiation δT/T=δ _?m /3 _?m . It is shown in this article that such a simple connection is not applicable due to:

1. The long time of recombination;
2. The fact that when regions withM<10¹⁵ M _⊙ become transparent for radiation, the optical depth to the observer is still large due to Thompson scattering;
3. The spasmodic increase of δ _?m/?m in recombination.

As a result the expected temperature fluctuations of relic radiation should be smaller than adiabatic fluctuations. In this article the value of δT/T arising from scattering of radiation on moving electrons is calculated; the velocity field is generated by adiabatic or entropy density perturbations. Fluctuations of the relic radiation due to secondary heating of the intergalactic gas are also estimated. A detailed investigation of the spectrum of fluctuations may, in principle, lead to an understanding of the nature of initial density perturbations since a distinct periodic dependence of the spectral density of perturbations on wavelength (mass) is peculiar to adiabatic perturbations. Practical observations are quite difficult due to the smallness of the effects and the presence of fluctuations connected with discrete sources of radio emission.     

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.     

Dynamics of Hα spicules according to spectral observations at various heights of the solar chromosphere

Almost simultaneous height sequences of 69 spicules in the Hα line have been studied. The spectra are obtained at six heights during 6 s on the east side of the solar disk with the 53-cm Lyot coronagraph of Abastumani Astrophysical Observatory. Radial velocities V _r, total intensities or equivalent widths W, full widths at half maximum of intensity (FWHM) at all heights are determined (about 300 profiles of the Hα line). It is found that:

1. Absolute values of radial velocities increase linearly with the height (see Equation (1));
2. variation of the sign of the radial velocity along single spicules was never observed.

These results combined with the findings on the spicules radial velocities and shifts obtained earlier (Kulidzanishvili and Nikolsky, 1978; Nikolsky and Platova, 1970) led us to the conclusion that the 5-min tangential oscillations of spicules involve the entire spicule at once. The intensity height scales for single spicules and for the chromosphere ‘in toto’ are determined; they turned out to be 2.5 × 10³ km and 1.9 × 10³ km respectively (see Equations (2) and (3)). The dependence curve of the Hα line half-widths Δλ on the height h is drawn. The Hα line half-width for those spicule groups which are traced at all heights (10 spicules) decreases with the height (Figure 4); for the majority (～60 spicules) it remains essentially constant. Non-thermal ‘turbulent’ velocities V _t, in Hα spicules are defined. A mean value of the ‘turbulent’ velocity V _t at T = 6000° appeared to be 20–30km s^?1. The hydrogen concentration in the spicules at 5000 km is 6 × 10¹¹ cm^?3.     

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.

     

On the nature of some active regions in the microwave range

The radio emission of a selected number of solar active regions has been investigated with high angular resolution at two frequencies: 10 and 17 GHz. By comparing the results of the two observations the following conclusions can be drawn:

1. The brightness temperature distribution of an active region is often composed of very bright cores of small dimension (angular extent θ?20″) imbedded in extended halos of lower brightness.
2. The radio emission of such structures as well as the degree of polarization can be explained with a thermal process. The halos can originate by pure thermal bremsstrahlung while in the case of the very bright cores found at 10 GHz (brightness temperature T _b?1–9 × 10⁶K) the emission at the harmonics of the gyrofrequency is needed.

     

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.

     

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.