The motion of a satellite in an axi-symmetric gravitational field

The equations for the variation of the osculating elements of a satellite moving in an axi-symmetric gravitational field are integrated to yield the complete first-order perturbations for the elements of the orbit. The expressions obtained include the effects produced by the second to eighth spherical harmonics. The orbital elements are presented in the most general form of summations by means of Hansen coefficients. Due to their general forms it is a simple matter to estimate the perturbations of any higher harmonic by simply increasing the index of summation. Finally, this paper gives the respective general expressions for the secular perturbations of the orbital elements. The formulae presented should be useful for the reductions of Earth-satellite observations and geopotential studies based on them.List of Symbols semi-major axis - C _jk ⁿ (, ) cosine functions of and - e eccentricity of the orbit - f acceleration vector of perturbing force - f sin²t - i inclination of the orbit - J _n coefficients in the potential expansion - M mean anomaly - n mean motion - p semi-latus rectum of the orbit - R, S, andW components of the perturbing acceleration - r radius-vector of satellite - r magnitude ofr - S _jk ⁿ (, ) sine functions of and - T time of perigee passage - u argument of latitude - U gravitational potential - true anomaly - V perturbing potential - G(M₊+m) (gravitational constant times the sum of the masses of Earth and satellite) - _n,k coefficients ofR component of disturbing acceleration (funtions off) - _n,k coefficients ofS andW components of disturbing acceleration (functions off) - mean anomaly at timet=0 - X ₀ ^n,m zero-order Hansen coefficients - argument of perigee - right ascension of the ascending node     

Thermal bifurcation in the upper photosphere inferred from heterodyne spectroscopy of OH rotational lines

We have obtained low noise (S/N > 10³), high spectral resolution (/ 10 ⁶) observations of two pure rotation transitions of OH from the solar photosphere. The observations were obtained using the technique of optically null-balanced infrared heterodyne spectroscopy, and consist of center-to-limb line profiles of a = 1 and a = 0 transition near 12 m. These lines should be formed in local thermodynamic equilibrium (LTE), and are diagnostics of the thermal structure of the upper photosphere. We find that the = 0R22(24.5)e line strengthens at the solar limb, in contradiction to the predictions of current one-dimensional photospheric models. Our data for this line support a two-dimensional model in which horizontal thermal fluctuations of order ±800 K occur in the region ₅₀₀₀ 10^–3–10^–2. This thermal bifurcation may be maintained by the presence of magnetic flux tubes, and may be related to the solar limb extensions observed in the 30–200 m region.Observations of the = 1R11(29.5)f line, at 885.643 cm^–1, show that it is anomalously weak in the photospheric spectrum. We argue that the source function in the core of this line has been substantially increased by interaction with the 9j-7i transition of Mgi at 885.524 cm^–1, which is itself too weak to appear in the disk center spectrum.Visiting Astronomer, Kitt Peak National Observatory, which is operated by the Association of Universities for Research in Astronomy, under contract with the National Science Foundation.     

Electrical conductivity of neutron star cores in the presence of a magnetic field

General theory of electrical conductivity of a multicomponent mixture of degenerate fermions in a magnetic fieldB, developed in the preceding article (this volume), is applied to a matter in neutron star interiors at densities ₀, where ₀ = 2.8×10¹⁴ g cm^–3 is the standard nuclear matter density. A model of free-particle mixture ofn, p, e is used, with account for appearance of ^–-hyperons at > _c , where _c 4₀. The electric resistivities along and acrossB, and , and the Hall resistivity _H are calculated and fitted by simple analytical formulae at _c and > _c for the cases of normal or superfluid neutrons provided other particles are normal. Charge transport alongB is produced by electrons, due to their Coulombic collisions with other charged particles; is independent ofB and almost independent of the neutron superfluidity. Charge transport acrossB at largeB may be essentially determined by other charged particles. If _c , one has = [1 + (B/B ₀)²] for the normal neutrons, and for the superfluid neutrons, while _H = B/B _e for both cases. HereB _e 10⁹ T ₈ ² G,B ₀10¹¹ T ₈ ² G, andT ₈ is temperature in units of 10⁸ K. Accordingly for the normal neutrons atBB ₀, the transverse resistivity suffers an enhancement, _{1/4 1}. When 5₀ andB varies from 0 toBB _p 10¹³ T ₈ ² G, increases by a factor of about 10³–10⁴ and _H changes sign. WhenBB _p , remains constant for the superfluid neutrons, and _H B ² for the normal neutrons, while _H B for any neutron state. Strong dependence of resistivity onB, T, and may affect evolution of magnetic fields in neutron star cores. In particular, the enhancement of at highB may noticeably speed up the Ohmic decay of those electric currents which are perpendicular toB.     

Bidimensional spectroscopy of network bright points

We develop an automatic, computer controlled procedure to select and to analyze the Network Bright Points (NBPs) on solar images. These have been obtained at the Sac Peak Vacuum Tower Telescope by means of the Universal Birefringent Filter and Zeiss H filters, tuned, respectively, along the profiles of the H, Mg-b1, Na-D2, and H lines.A structure is identified as an NBP if at the wavelength H- 1.5 A its maximum intensity is greater than I + 3 and its area is greater than 1.5 arc sec² at I + 1.5, where I is the mean value and the standard deviation of the intensity distribution on the image. Each detected NBP is then searched and confirmed in all the remaining 31 images at different wavelengths.For each NBP several parameters are measured (position, area, mean and maximum contrast, Dopplergram velocity, compactness, and so on) and some identification constraints are applied.The statistical analysis of the various parameter distributions, for NBPs present within an active region and its surroundings, shows that two types of NBPs can be identified according to the value of their mean contrast C _min the H- 1.5 Å image (C _m 0.1 type I, C _m> 0.1 type II). The type I NBPs (all occurring on the boundaries of the supergranular network) appear to be much more frequent (180/26) than the type II ones.The size A of type I NBPs is less than 1.0 arc sec for H/H wings but of the order of 1.2 arc sec for Na-D2 and Mg-bl. The mean contrast C _m is around the value of 10% along the Na-D2 and Mg-bl profiles and of 20% along the H/H wings.The C _m - A scatter diagrams show, for the photospheric radiation (h < 100 km), a narrow range of variability for C _min correspondence with a wide range for A. For radiation orginated at higher levels (h > 200 km), the C _m- A scatter diagrams seem to indicate, even if with a large variance, that the highest C _m's tend to correspond to the highest A values.The mean Doppler shift is close to zero for Na-D2 and Mg-bl lines but negative (downward motion) for H and H lines.The type II NBPs tends to be preferentially located in the neighbourhood of small, compact sunspots and their detectability is almost constant through all the 4 studied line profiles. No conclusions can be derived on the mean size, contrast and Doppler shift values because their distributions are too dispersed. The only positive information is that its C _m- A scatter diagram, in H and H wings, indicates a wide range of variability for C _m in correspondence with very narrow range of variability for A.     

Pulsar radio emission from expanding charge sheets

We semi-quantitatively calculate the distribution of energy in frequency and angle emitted from a sheet of charges that are moving out relativistically along dipolar magnetic field lines originating near the magnetic polar caps of a rotating neutron star. The angular distribution is conical with the angle of maximum intensity varying with frequency as ^–1/4 for _c 2 _c /(R _M ²), whereR_M is the initial angular radius of the charge sheet at the surface of the star of radiusR. At higher frequencies the width of the angular cone remains constant. The radiation is linearly polarized with the polarization vector in the plane of the line of sight and the magnetic axis. A sheet of uniform charge density and finite thickness has a frequency spectrum that varies from ^–3/2 to ^–4 for _c and _c , respectively. These features are in good general agreement with the observed characteristics of the intensity, pulse shape, and frequency spectrum of the radio pulses from pulsars.Operated by Associated Universities, Inc., under contract with the National Science Foundation.     

A novel methodology for radiative transfer in a planetary atmosphere

A novel methodology for evaluating the field of anisotropically scattered radiation within a homogeneous slab atmosphere of arbitrary optical thickness is provided. It departs from the traditional radiative transfer approach in first considering that the atmosphere is illuminated by an isotropic light source. From the solution of this problem, it subsequently proceeds to that for the more conventional case of monodirectional illumination. The azimuthal dependence of the field is separated in the usual manner by an harmonic expansion, leaving a problem in four dimensions (=optical depth, ₀=thickness, , =directions of incidence and scattering) which, as is well known, is numerically extremely inconvenient. Two auxiliary radiative transfer formulations of increasing dimensionality are considered: (i) a transfer equation for the newly introduced functionb ^m(,,₀) with Sobolev's function ^m(,₀) playing the role of a source-function. Because the incident direction does not intervene, ^m is simply expressed as a single integral term involvingb ^m. For bottom illumination, an analogous equation holds for the other new functionh ^m(,,₀). However, simple reciprocity relations link the two functions so that it is only necessary to considerb ^m; (ii) a transfer equation for the other new functiona ^m(,,,₀) with a source-function provided by Sobolev's functionD ^m(,,₀). For bottom illumination, another functionf ^m(,,,₀) is introduced; by a similar argument using reciprocity relations,f ^m is reduced toa ^m rendering necessary only the consideration ofa ^m. However, a fundamental decomposition formula is obtained which shows thata ^m is expressible algebraically in terms of functions of a single angular variable. The functions ^m andD ^m are shown to be the values in the horizontal plane ofb ^m anda ^m, respectively. The other auxiliary functionsX ^m andY ^m are also expressed algebraically in terms ofb ^m. These results enable one to proceed to the final step of evaluating the radiation field for monodirectional illumination. The above reductions toalgebraic relations involving only the functionb ^m appear to be more advantageous than Sobolev's (1972) recent approach; they also circumvent some basic numerical difficulties in it. We believe the present approach may likewise prove to be superior to most (if not all) other methods of solution known heretofore.This paper presents the results of one phase of research carried out at the Jet Propulsion Laboratory under Contract No. NAS-7-100 sponsored by the National Aeronautics and Space Administration.     

Thermal and electrical conductivities of crystals in neutron stars and degenerate dwarfs

Thermal and electrical conductivities due to electron scattering on phonons are calculated for degenerate cores of white dwarfs and envelopes of neutron stars for wide ranges of density, temperature and ion charge. In the stellar zones, in which T _pi(Z^1/3e²/_F) (_piis the ion plasma frequency and _F the Fermi velocity of electrons), the main contribution into scattering comes from the Umklapp processes. In the zones with lowerT, the Umklapp processes are frozen out, that results in a sharp growth of electrical and thermal conductivities. This, for instance, should make nuclear burning more stable in such zones.     

Analysis of extreme-ultraviolet spectroheliograms of solar prominences

The optical depth at the head of the Lyman continuum, _H, is determined at a number of positions in three hedgerow prominences using spectroheliograms (5 × 5 resolution) of C III 977, LC 896, and O IV 554 observed with the Harvard experiment on Skylab. At heights greater than 10 above the limb the maximum value of _H is 30 to 50, which occurs at the central part of the prominences. For one of the prominences the determination of _H is found to be consistent with data from spectroheliograms of Mg X 625. The degree of ionization of hydrogen is estimated from the intensity of LC 896 at _H 1. In the central part of a model prominence N _P/N _HI1.9 for a reasonable range of the electron densities, where N _P and N _HI are the proton density and the neutral hydrogen density, respectively.     

Hα and radio observations of LSI+61o303

Quasi simultaneous H and radio observations of LSI+61^o303 during August-September 1993 are presented. The radio data show that during the studied epoch the outburst has peaked at radio phase 0.6, at level100 mJy. No significant variations in the H profile at phases 0.5 to 0.65 have been detected. A remarkable increase of the EW and FWHM of the H blue peak is observed at radio phase 0.23. Possible reasons are considered.     

The effect of chemical abundance oscillations on the pulsational properties of A 5M ⊙ hydrogen-helium star

A model of a first generation intermediate star of 5M , with Z=0 has been considered. The model is at an advanced stage of its evolution and has a double shell burning. It burns helium in the inner shell, and hydrogen, via CNO cycle, in the outer shell. =(log/log) _T and _T =(log/logT) were computed allowing for the oscillations of the relative mass abundance of the reagents in nuclear reactions. Including =(log/log) _T and =(log/logT) of mean molecular weight and the effect of the oscillations of abundances due to nuclear reactions, stability was studied. Contrary to the results of the static calculations, we found that instability due to the excitation mechanism provided by the high temperature sensitivity of energy generation rate propagates up to the surface. Thus the model in question was found to be unstable against radial adiabatic pulsations, in its fundamental mode.     

Gravitational radiation and spiralling time of close binary systems (V)

Forty-six binary systems with their primary component masses between 2M and 3M have been considered for gravitational radiation study. Power output by gravitational radiation (P _B ) and spiral time ₀ for all individual systems have been evaluated. A relation has been given betweenP _B and ₀. The rate of decrease of orbital period (P) has also been given for 10 eccentric orbit systems.     

Effect of cavitation on spherical blast waves

For spherical blast waves propagating through a self-gravitating gas with an energy inputE =E ₀ t , whereE is the energy released up to timet,E ₀ is a functional constant, and is a constant, kinetic, internal heat, and gravitational potential energies have been computed. Taking the parameterA ², which characterises the gravitational field, equal to 2, variations of the percentages of these energies for =0, 1/2, 4/3, and 3 with shock strength have been presented. For =3, the effect of cavitation on the percentages of kinetic energy and internal heat energies has been explored.     

Intensity fluctuations in Fraunhofer lines

An expression is derived for the fluctuation (t) in emergent intensity (observed at some wavelength in a Fraunhofer line or the continuum) caused by a perturbation in temperature (z, t) in the Sun's atmosphere. If the contribution function for the observed intensity is single-peaked near z and if (z) and p(z) are not too rapidly varying, then (t) m (z , t)+N p(z , t) where m and N depend on the structure of the atmosphere. We compute M, N, and contribution functions for several values of and in the inner wings of the K line (13933 Caii).Presently on leave of absence from the Institute for Astronomy, Honolulu, Hawaii.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

Planetary accretions in jet streams

The problem of planetary accretion in a jet stream is studied using the model developed by Alfvén and Arrhenius. We find that there are basically three types of planetary accretion: namely, fast process _c < _i , slow process _c ~ _i and delayed process _c > _i where _c is the characteristic time of the occurrence of catastrophic accretion and _i the time-scale of mass injection to the planetary system (3×10⁸ yr). These different time scales of accretion are found to be closely related to the primordial thermal profiles and equatorial inclinations of the planets. Finally, Venus' retrograde rotational spin is shown to be a possible result of accretion process in a jet stream.     

Expanding and superdense astrophysical systems in general relativistic hierarchical cosmology

A semi-continuous hierarchy, (i.e., one in which there are galaxies outside clusters, clusters outside superclusters etc.), is examined using an expression of the field equations of general relativity in a form due to Podurets, Misner and Sharp. It is shown (a) that for a sufficiently populous hierarchy, the thinning factor( _i+1/ _i [r _i /r _i+1] is approximately equal to the exponentN in a continuous density law (=aR ^–N) provided (r _i /r _i+1)^3-1; (b) that a hierarchical Universe will not look decidedly asymmetric to an observer like a human being because such salient observers live close to the densest elements of the hierarchy (viz stars), the probability of the Universe looking spherically symmetric (dipole anisotropy0.1 to such an observer being of order unity; (c) the existence of a semi-continuous or continuous hierarchy (Peebles) requires that 2 if galaxies, not presently bound to clusters were once members of such systems; (d) there are now in existence no less than ten arguments for believing 2, though recent number counts by Sandageet al. seem to be in contradiction to such a value; (e) Hubble's law, withH independent of distance, can be proved approximately in a relativistic hierarchy provided (i)N=2, (ii)2GM(R)/c ² R1; (iii)Rc (iv)M0 in a system of massM, sizeR (f) Hubble's law holds also in a hierarchy with density jumps; (g)H100 km s^–1 Mpc^–1; (h) objects forming the stellar level of the hierarchy (in a cosmology of the Wilson type) must once have had 2GM/c ² R1; (i) there is a finite pressurep=2Ga in all astrophysical systems (a=R ^N ,N2); (j) for the Galaxy, theory predictsp _G7×10^–12 dyn cm^–2, observation givesp _G5×10^–12 dyn cm^–2; (k) if the mass-defect (or excess binding energy) hypothesis is taken as a postulate, all non-collapsed astrophysical systems must be non-static, and any non-static, p0 systems must in any case be losing mass; (1) the predicted mass-loss rate from the Sun is 10¹² g s^–1, compared to 10¹¹ g s^–1 in the observed solar wind; (m) the mass-loss rates known by observation imply timescales of 5×10⁹ years for the Sun and 10¹⁰ years for other astrophysical systems; (n) degenerate superdense objects composed of fermions must haveN-2 if they were ever at their Schwarzschild radii and comprised a finite numberN _B of baryons; (o)N _B10⁵⁷N for degenerate fermion and boson systems; (p)285-4; (q) the metric coefficients for superdense bodies give equations of motion that imply equal maximum luminosities for all evolving superdense bodies (L _max10⁵⁹ erg s^–1); (r) larger bodies have longer time-scales of energy radiation atL _max (10^–5 s for stars,1 h for QSO's) (s) expansion velocities are c soon after the initial loss of equilibrium in a superdense object; (t) if the density parametera(t) in aR ^–N isa=a (non-atomic constants of physicsc, G, A), andA, thenN=2; (u) N2 is necessary to giveMM at the stellar level of the hierarchy;(v) systems larger than, and including, galaxies must have formed by clumping of smaller systems and not (as advocated by Wertz and others) in a multiple big bang.     

Theory of the Trojan asteroids,IV

In a previous publication (1977) the author has constructed a family () of long-periodic orbits in the Trojan case of the restricted problems of three bodies. Here he constructs the domain of the analytical solution of the problem of the motion, excluding the vicinity of thecritical divisor which vanishes at the exact commensurability of the natural frequencies ₁ and ₂. In terms of thecritical masses m_j(₂), or the associatedcritical energies _j ² (m), is the intersection of the intervals ofshallow resonance, of the form. Inasmuch as the intervals |₂– _j ² |<_j ofdeep resonance aredisjoint, it follows that (1) the disjointed family () embraces the tadpole branch, 0²1, lying in: and (2) despite the clustering of _j ² (m) atj=, the family () includes, for ²=1, an asymptoticseparatrix that terminates the branch in the vicinity of the Lagrangian pointL ₃.In a similar manner, the family () can be extended to the horseshoe branch 1<_{2 2} ² .     

Some applications of Jacobi dynamics to the stellar evolution

The stars in the Main Sequence are seen as a hierarchy of objects with different massesM and effective dynamical radiiR _eff=R/ given by the stellar radii and the coefficients for the inner structure of the stars.As seen in a previous work (Paper I), during the lifetime in the Main SequenceR _eff(t) remains a near invariant when compared to the variation in the time ofR(t) and (t).With such an effectiveR _eff one obtains the amounts of actionA _c(M), the effective densities _eff(M)=(M)³(M), the densities of action and of energy (or mean presures in the stellar interior)a _c(M),e _c(M), and the potential energiesE _p(M).The amounts of action areA _cM ^k withk1.87 for the M stars,k5/3 for the KGF stars, andk1.83 for the A and earlier stars, representing very simples conditions for the other dynamical parameters. For instancek5/3 means a near invariant effective density _eff for the KGF stars, while for such stars the mean densities and coefficients present the strongest variations with masses (M)M ^–1.81, (M)M^0.6.The cases for the M stars (e _c(M)M ^–1) and for the A and earlier stars (betweena _c(M)=constant and _eff(M)M ^–1) and also discussed. These conditions for the earlier stars also represent reasonable mean values for the whole stellar hierarchy in the range of masses 0.2M M25M .With all this, one can build dynamical HR diagrams withA _c(M), E_p(M), _eff M ^–p , etc., whose characteristics are analogous to these in the photometrical HR diagram. A comparison is made betweenA _c(M) from the models here and the HR diagram with the best known stars of luminosity classes IV, V, and white dwarfs.The comparison of the potential energiesE _p(M)M ^–p according to the stellar models used here and the observed frequency function (MM _–q (number of stars in a given interval of masses) from different authors suggests the possibility that the productE _p(M)(M) is a constant, but this must be confirmed with further studies of the function (M) and its fine structure.There are analogies between the formulation used here for the stellar hierarchy and other physical processes, for instance, in modified forms of the Kolmogorov law of turbulence and in the formulation used for the hierarchy of molecular clouds in gravitational equilibrium. Besides, the function of actionA _c(M) for the stars has analogous properties to the relations of angular momenta and massesJ(M) for different types of objects. The cosmological implications of all this are discussed.     

Infrared emission of dust grains in the clusters of galaxies

The observations of the reddening of the distant galaxies and the weak diffuse radiation in the clusters of galaxies can be interpreted as a consequence of the presence of dust grains in the intergalactic medium. When allowance is made for the destruction of the grains in collision with particles of the hot gas, its lifetime is about 10⁷–10⁸ yr at a gas concentrationn _g 10^–3 cm^–3. The detection of the infrared (IR) emission from the galaxy clusters might be the test for the proof of the presence of dust grains in the intergalactic medium. In this paper the estimates of the expected intensities and fluxes of IR emission for the spectral region 50–300 are presented for two galaxy clusters in Coma and Perseus. The parameters of the hot gas spatial distribution are chosen from X-ray observations. Having assumed that intergalactic dust can be ejected only from the galaxies, we used such a model for intergalactic dust grains which explains very well the interstellar dust effects. It is shown that the dust temperature, which is determined from the general energetic balance of the dust grains, can achieve some scores of degrees of Kelvin. Two models of the dust spatial distribution are considered. It is found that the maximum of IR flux for the Coma cluster lies near =100 and the same for the Perseus cluster near 50–70. The total fluxes of IR emission from these clusters are about 10⁵–10⁶ Jy and can be detected by modern observational methods.     

Mass and orbit estimation of Planet X via a family of comets

The crucial assumption of this paper is, that the observed clustering of aphelion distances of intermediate-period comets in the 70–90 AU range is due to the influence of a tenth planet, called Planet X. We contribute to the search for Planet X a new and extended evaluation of a family of comets assumed to be Planet X's family of comets.By averaging the aphelion distances of comets that belong to a transplutonic family of comets, we get Planet X's semi-major axis a _x = (83.0 ± 5.3) AU. The comets' orbits also yield the upper limit of the planet's orbital eccentricity e _x - 0.019. If this planet played an important part in sending quasi-periodic comet showers to the inner solar system, we can calculate its orbital inclination i _x = 46 .1 ± 3 .6. By distributing all planets' masses into the heliocentric, torus-like zones, in which they were formed, we get the density distribution of the primordial solar nebula. Extrapolating this distribution we find the mass of the planet M _x = (5.1_–2.4 ^+3.6 M _Earth. A few plausible assumptions (e.g. Uranus and Neptune perturbations being caused by Planet X) lead to Planet X's actual location with declination and eccliptic longitude being = 57 ± 17 and = 54 ± 34 , respectively (1989.5 position). In addition, we give Planet X's apparent brightness dependent on its unknown albedo. All those properties and predictions are more or less in agreement with earlier work on Planet X.