Radiation pressure effects in the oscillations of compressible rotating homogeneous spheroids

Earlier models of compressible, rotating, and homogeneous ellipsoids with gas pressure are generalized to include the presence of radiation pressure. Under the assumptions of a linear velocity field of the fluid and a bounded ellipsoidal surface, the dynamical behaviour of these models can be described by ordinary differential equations. These equations are used to study the finite oscillations of massive radiative models with masses 10M and 30M in which the effects of radiation pressure are expected to be important.Models with two different degrees of equilibrium are chosen: an equilibrium (i.e., dynamically stable) model with an initial asymmetric inward velocity, and a nonequilibrium model with a nonequilibrium central temperature and which falls inwards from rest. For each of these two degrees of equilibrium, two initial configurations are considered: rotating spheroidal and nonrotating spherical models.From the numerical integration of the differential equations for these models, we obtain the time evolution of their principal semi-diametersa ₁ anda ₃, and of their central temperatures, which are graphically displayed by making plots of the trajectories in the (a ₁,a ₃) phase space, and of botha ₁ and the total central pressureP _c against time.It is found that in all the equilibrium radiative models (in which radiation pressure is taken into account), the periods of the oscillations of botha ₁ andP _c are longer than those of the corresponding nonradiative models, while the reverse is true for the nonequilibrium radiative models. The envelopes of thea ₁ oscillations of the equilibrium radiative models also have much longer periods; this result also holds for the nonequilibrium models whenever the envelope is well defined. Further, as compared to the nonradiative models, almost all the radiative models collapse to smaller volumes before rebouncing, with the more massive model undergoing a larger collapse and attaining a correspondingly larger peakP _c.When the mass is increased, the dynamical behavior of the radiative model generally becomes more nonperiodic. The ratio of the central radiation pressure to the central gas pressure, which is small for low mass models, increases with mass, and at the center of the more massive model, the radiation pressure can be comparable in magnitude to the gas pressure. In all the radiative models, the average periods as well as the average amplitudes of both thea ₁ andP _c oscillations also increase with mass.When either rotation or radiation pressure effects or both are included in the equilibrium nonradiative model, the period of the envelope of thea ₁ oscillations is increased. The presence of rotation in the equilibrium radiative model, however, decreases this period.Some astrophysical implications of this work are briefly discussed.     

Hydrogen ionization and n=2 population for model spicules and prominences

Using slab model atmospheres that are irradiated from both sides by photospheric, chromospheric, and coronal radiation fields we have determined the ionization and excitation equilibrium for hydrogen.The model atom consists of two bound levels (n = 1 and n = 2) and a continuum. Ly- was assumed to be optically thick with the transition in detailed radiative balance. The Balmer continuum was assumed to be optically thin with the associated radiative ionization dominated by the photospheric radiation field (T _rad = 5940 K). The ionization equilibrium was determined from an exact treatment of the radiative transfer problem for the internally generated Ly-c field and the impressed chromospheric and coronal field (characterized by T _rad = 6500K).Our calculations corroborate the hypothesis that N₂, the n = 2 population density, is uniquely determined by the electron density N _e. We also present ionization curves for 6000K, 7500K, and 10000K models ranging in total hydrogen density from 1 × 10¹⁰/cm³ to 3 × 10¹²/cm³. Using these curves it is possible to obtain the total hydrogen density from the n = 2 population density in prominences and spicules.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

Clues to the mode of excitation of Fe x ions in the solar corona from the 1980 eclipse observations

The line and continuum intensities deduced from the multislit spectra of the (Fe X) coronal emission line taken at the 1980 eclipse are used to discuss the relative roles of radiative and collisional excitation mechanisms. It is shown that for R/R < 1.2, collisional excitation is the predominant mode. Collisional as well as radiative excitation is equally important for 1.2 < R/R < 1.4, whereas beyond 1.4 R radiative excitation becomes dominant. The line width measurements indicate that a large number of locations have half-widths around 1.3 Å. The maximum half-width is reached at 1.4 R with an average value of 1.6 Å.     

Satellite theory

In this paper dynamical characteristics of satellites are outlined by classifying the satellites into three categories according to the values of the solar tidal factor (n/n)² which is the disturbing factor due to the sun and the oblateness factor of the primary planetJ ₂/a ². For inner satellites (n/n)² is much smaller thanJ ₂/a ² and there are several pairs among them, for which the mean motions are commensurable to each other, and for some of them secular accelerations in the mean longitudes have been detected. For outer satellites (n/n)² is much larger and the solar perturbations are dominant. For intermediary satellites the motion of the pole of the orbital plane is not so simple as those of the satellites of the other categories.     

On the evolution of obstars from the main sequence to the Helium exhaustion phase

The evolution of a 20M star of extreme Population I (X=0.602,Z=0.044) from the main sequence through the He-exhaustion has been calculated with particular attention to the treatment of the semiconvective zones. Two different relations as convective instability criteria have been adopted _{R a} +(/(4–3)) (d ln /d lnP) or _{R a} ). The He-burning phase is substantially different in the two cases. In the first evolutionary picture the star begins the He-burning phase as red supergiant, while in the second one as blue supergiant. The two evolutions are compared with the observational blue to red supergiants ration _B/n _R in connection with the problem of the existence of the neutrino processes and taking into account the possible presence of a high percentage of binary systems among supergiants.We draw the conclusion that the massive star evolution seems to be still completely uncertain and also not able to give a decisive test for or against the neutrino emission.     

Atmospheric Disturbances and Radiation Impulses Caused by Large-Meteoroid Impacts on the Surface of Mars. II. Radiation Impulse Characteristics and Parameters of the Warm Layer

This paper presents the results of the calculation of spectral and angular characteristics of radiation emitted by the disturbed region after the vertical impacts on the Martian surface of stony meteoroids with radii R ₀ from 1 to 100 m at speeds of 11–20 km/s. The time dependences are given for the density of the radiation flux incident on horizontal surface areas located at different distances from the impact point. For small impactors (R ₀= 1 m), the heating of the surface and surrounding gas by the radiation impulse from the hot region formed after the impact is insignificant even at the crater edge. However, the radiation impulse that heats up the surface is also emitted during the meteoroid flight through the atmosphere. Although this heating is insufficient to initiate evaporation, heat transfer by turbulent diffusion leads to the formation of a layer with temperatures that substantially exceed the initial temperature of the atmosphere. For large impactors (R ₀ = 100 m), the total specific impulse of the plume radiation gives rise to the emergence of the heated layer with a thickness on the order of several meters within several kilometers of the impact point. The formation of this warm layer may lead to the formation of a high-speed jet moving along the Martian surface as well as eddies at the front of the precursor with a subsequent intense rise of dust.     

Some requirements of a colliding comet source of gamma ray bursts

Colliding comets in the Solar System may be an important source of gamma ray bursts. The spherical gamma ray comet cloud required by the results of the Venera Satellites (Mazets and Golenetskii, 1987) and the BATSE detector on the Compton Satellite (Meeganet al., 1992a, b) is neither the Oort Cloud nor the Kuiper Belt. To satisfy observations ofN(>P _max) vsP _max for the maximum gamma ray fluxes,P _max > 10^–5 erg cm^–2 s^–1 (about 30 bursts yr^–1), the comet density,n, should increase asn a ¹ from about 40 to 100 AU wherea is the comet heliocentric distance. The turnover above 100 AU requiresn a ^–1/2 to 200 AU to fit the Venera results andn a ^1/4 to 400 AU to fit the BATSE data. Then the masses of comets in the 3 regions are from: 40–100 AU, about 9 earth masses,m _E; 100–200 AU about 25m _E; and 100–400 AU, about 900m _E. The flux of 10^–5 erg cm^–2 s^–1 corresponds to a luminosity at 100 AU of 3 × 10²⁶ erg s^–1. Two colliding spherical comets at a distance of 100 AU, each with nucleus of radiusR of 5 km, density of 0.5 g cm^–3 and Keplerian velocity 3 km s^–1 have a combined kinetic energy of 3 × 10²⁸ erg, a factor of about 100 greater than required by the burst maximum fluxes that last for one second. Betatron acceleration in the compressed magnetic fields between the colliding comets could accelerate electrons to energies sufficient to produce the observed high energy gamma rays. Many of the additional observed features of gamma ray bursts can be explained by the solar comet collision source.     

The solution of radiation transfer problem for a slab with space-dependent albedo and anisotropic scattering

The transport of thermal radiation has been considered within a finite slab which absorb and scatter anisotropically. The problem involves the space-dependent single-scattering albedow(x). Two approximations are taken forw(x). In the first it is represented in exponential form asw(x)=w ₀ exp(–x/s), wherew ₀ ands are given constants andx is the optical variable. The second approximation assumes the formw(x) = _r=0 ^R d _r ^* p _r (x/a), whered _r ^* are known expansion coefficients anda is the half optical thickness of the slab. Analytic expressions for the forward, backward radiation intensities and fluxes are given in each approximation. The solution of the linear transport equation is performed on the basis of integral Fourier transforms.     

Bimodality and the Hale cycle

Because of the bimodal distribution of sunspot cycle periods, the Hale cycle (or double sunspot cycle) should show evidence of modulation between 20 and 24 yr, with the Hale cycle having an average length of about 22 yr. Indeed, such a modulation is observed. Comparison of consecutive pairs of cycles strongly suggests that even-numbered cycles are preferentially paired with odd-numbered following cycles. Systematic variations are hinted in both the Hale cycle period and R _sum (the sum of monthly mean sunspot numbers over consecutively paired sunspot cycles). The preferred even-odd cycle pairing suggests that cycles 22 and 23 form a new Hale cycle pair (Hale cycle 12), that cycle 23 will be larger than cycle 22 (in terms of R _M, the maximum smoothed sunspot number, and of the individual cycle value of R _sum), and that the length of Hale cycle 12 will be longer than 22 yr. Because of the strong correlation (r = 0.95) between individual sunspot cycle values of R _sum and R _M, having a good estimate of R _Mfor the present sunspot cycle (22) allows one to predict its R _sum, which further allows an estimation of both R _Mand R _sum for cycle 23 and an estimation of R _sum for Hale cycle 12. Based on Wilson's bivariate fit (r = 0.98), sunspot cycle 22 should have an R _Mequal to 144.4 ± 27.3 (at the 3- level), implying that its R _sum should be about 8600 ± 2200; such values imply that sunspot cycle 23 should have an R _sum of about 10500 ± 2000 and an R _Mof about 175 ± 40, and that Hale cycle 12 should have an R _sum of about 19100 ± 3000.     

Self-similar space-times

An exact solution of the transport equation in radiative transfer for an axially symmetric Rayleigh scattering problem in semi-infinite planetary atmosphere both for emergent intensity and intensity at any optical depth has been derived with the help of the Laplace transform and the Wiener-Hopf technique, and by use of the constancy of net flux. Chandrasekhar's results for emergent intensity have been verified. New expressions for theH _l andH _r functions have been obtained.     

Radiative heating and the buoyant rise of magnetic flux tubes in the solar interior

We study the effect of radiative heating on the evolution of thin magnetic flux tubes in the solar interior and on the eruption of magnetic flux loops to the surface. Magnetic flux tubes experience radiative heating because (1) the mean temperature gradient in the lower convection zone and the overshoot region deviates substantially from that of radiative equilibrium, and hence there is a non-zero divergence of radiative heat flux; and (2) the magnetic pressure of the flux tube causes a small change of the thermodynamic properties within the tube relative to the surrounding field-free fluid, resulting in an additional divergence of radiative heat flux. Our calculations show that the former constitutes the dominant source of radiative heating experienced by the flux tube.In the overshoot region, the radiative heating is found to cause a quasi-static rising of the toroidal flux tubes with an upward drift velocity 10^-3|| cm s^-1, where _e – _ad < 0 describes the subadiabaticity in the overshoot layer. The upward drift velocity does not depend sensitively on the field strength of the flux tubes. Thus in order to store toroidal flux tubes in the overshoot region for a period comparable to the length of the solar cycle, the magnitude of the subadiabaticity (< 0) in the overshoot region must be as large as 3 × 10^–4. We discuss the possibilities for increasing the magnitude of and for reducing the rate of radiative heating of the flux tubes in the overshoot region.Using numerical simulations we study the formation of -shaped emerging loops from toroidal flux tubes in the overshoot region as a result of radiative heating. The initial toroidal tube is assumed to be non-uniform in its thermodynamic properties along the tube and lies at varying depths beneath the base of the convection zone. The tube is initially in a state of neutral buoyancy with the internal density of the tube plasma equal to the local external density. We find from our numerical simulations that such a toroidal tube rises quasi-statically due to radiative heating. The top portion of the nonuniform tube first enters the convection zone and may be brought to an unstable configuration which eventually leads to the eruption of an anchored flux loop to the surface. Assuming reasonable initial parameters, our numerical calculations yield fairly short rise times (2–4 months) for the development of the emerging flux loops. This suggests that radiative heating is an effective way of causing the eruption of magnetic flux loops, leading to the formation of active regions at the surface.The National Solar Observatory is one of the National Optical Astronomy Observatories by the Association of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation.     

Probable causes of long-period variations in the Uranian geometrical albedo

We use the technique developed in our previous studies to determine the ratios of optical depth components τ _a /τ _R and τ_κ/τ _S on the basis of the observational data on the Uranian geometrical albedo for the years 1981, 1993, and 1995. Here τ _a and τ _R are the aerosol and gas scattering components, τ _S = τ _a + τ _R , and τ_κ is the absorption component of the effective optical depth at which the intensity of the diffusely reflected radiation is formed. The ratios turn out to be different for different years. This phenomenon is caused by the horizontal inhomogeneity of the aerosol component distribution over the Uranian disk.     

Linear waves in a radiating and scattering grey medium

The problem on linear waves in a radiating and scattering grey medium is studied using Whitham's method. Analysis of the basic equations distinguishes two limiting cases: the one is theequilibrium case in which the energy exchange between the gas and radiation plays an essential role, and the other is theScattering case in which the effect of energy exchange is negligible. A new type ofradiation acoustic wave with the speed is found in the scattering case. The governing equations for linearized one-dimensional flow are reduced to one equation of radiative acoustics valid to order 1/c, and the criterion for the two limiting cases is derived from studying this equation. The harmonic solution is analytically studied to show that theeffective optical depth corresponding to the wavelength of perturbation gives the measure of the interaction between the gas and radiation. When _eff1, the sound speeda _g ² =P _g / and the propagating speed of radiative disturbancea _f ² =fc ² appear as the modified classical and radiation-induced modes respectively, wheref is the Eddington factor. When _eff1, the isentropic sound speeda _s ² =(P _g +P _r / appears in the equilibrium case, and the radiation acoustic speeda _A ² appears in the scattering case. The dispersion relation of the harmonic solution is numerically calculated. The result shows that the wave pattern depends critically on the ratio=P _g /(P _g +P _r ). When , the speeda _S anda _A arise from the modified classical mode, and when , they originate from the radiation-induced mode.     

Polar coronal plumes

Polar coronal plumes are modeled using concentrations of magnetic flux at 1.01R , and assuming the field is current-free, or a potential field. Identifying the density enhancement of plumes with magnetic flux concentration produces good agreement between 1.01R and 1.10R , for model conditions of a large background magnetic field and a plume separation of 50 000 to 70 000 km at the base. Beyond 1.10R , both plumes and the potential field diverge very nearly as r ².Also Department of Astrogeophysics, University of Colorado, Boulder, Colo. 80309, U.S.A. Presently visiting Stanford University Institute for Plasma Research, Via Crespi, Stanford, Calif. 94303, U.S.A.     

Nonlinear pulsations of a very luminous star

The radial nonlinear pulsations of a model withM=0.8M ,M _bol=–6 mag andT _eff=5500 K have been studied. The pulsations are shown to exist in the form of the standing wave only in the innermost layers withR<0.2R _ph. In the outer layers, the standing wave transforms into running waves, the frequency of which decreases with an increasingR. the pulsation period at the photosphere is found to be twice as long as the pulsation period atR<0.2R _ph. The difference between the pulsation periods causes alternation of deep and shallow minima in the temporal dependences of the kinetic energy and radii and can be used to explain the nature of RV Tau variables. It is shown that at the distanceR>3R _ph, the time-independent mass flux caused by shocks takes place. The rate of mass loss is found to beM10^–4 M yr^–1.     

Perturbations of a close-earth satellite due to sunlight diffusely reflected from the Earth

A semianalytical method has been developed to calculate the radiation-pressure perturbations of a close-Earth satellite due to sunlight reflected from the Earth. It is assumed that the satellite is spherically symmetric and that the solar radiation is reflected from the Earth according to Lambert's Law. To account for the increasing reflectivity of the Earth toward the poles, its albedo is assumed to have a latitudinal dependence given bya=a ₀ +a ₂ sin². The effect of the terminator on the perturbations has been neglected. The perturbations within a particular revolution are given analytically, while the long-range perturbations are obtained by accumulation.     

Microwave radiation from magnetic stars

Cyclotron microwave emission from magnetic stars is considered, assuming that they have coronae with the temperatureT10⁷ K and the emission measureEM10⁵⁴ cm^–3. It has been shown that the cyclotron radiation from a star with a dipole magnetic field has a specific spectrum with a maximum in the frequency rangesv _o/2 >v >sv _o/2 (s being the number of cyclotron harmonic, andv _o the gyrofrequency corresponding to the polar magnetic field) and radiation flux decreasing towards lower frequencies asv _4/3. The frequency of the spectrum maximum depends on the angle between the line-of-sight and the magnetic axis of the star. The observed radiation from a rotating magnetic star can be modulated with a modulation depth of about 0.2 at frequencies near maximum. The radiation is partially circularly-polarized in the sense of an extraordinary mode. The degree of polarization is almost constant at frequenciesv >sv _o/2 and increases with frequency atv >sv _o/2. The estimation of cyclotron radio fluxes of the nearest magnetic stars shows that they are observable in microwaves by means of modern radio astronomy.     

Solar irradiance from Nimbus-7 compared with ground-based photometry

We have compared total solar irradiance from Nimbus-7 with ground-based photometry from the San Fernando Observatory (SFO) for 109 days between June 1 and December 31, 1988. We have also included in some analyses NOAA-9 SBUV2 data orF10.7 radio flux. The Nimbus-7 data are from orbital samples, averaged to the mean time of observation at SFO. Using the same parameters as in Chapmanet al. (1992), the multiple regression gives anR ² = 0.9131 and a solar minimum irradiance,S ₀, = 1371.76 ± 0.18 W m^–2 for the best fit.     

Further considerations on contracting solar nebula

Prentice (1978a) in his modern Laplacian theory of the origin of the solar system has established the scenario of the formation of the solar system on the basis of the usual laws of conservation of mass and angular momentum and the concept of supersonic turbulent convection that he has developed. In this, he finds the ratio of the orbital radii of successively disposed gaseous rings to be a constant - 1.69. This serves to provide a physical understanding of the Titius-Bode law of planetary distances. In an attempt to understand the law in an alternative way, Rawal (1984) starts with the concept of Roche limit. He assumes that during the collapse of the solar nebula, the halts at various radii are brought about by the supersonic turbulent convection developed by Prentice and arrives at the relation: R _p= Ra^p, where R _pare the radii of the solar nebula at various halts during the collapse, R the radius of the present Sun and a = 1.442. a is referred here as the Roche constant. In this context, it is shown here that Kepler's third law of planetary system assumes the form: T _p = T ₀(a^3/2)^p, where T _p are the orbital periods at the radii R _p, T ₀ - 0.1216d - 3 h, and a the Roche constant. We are inclined to interpret T ₀' to be the rotation period of the Sun at the time of its formation when it attained the present radius. It is also shown that the oribital periods T _pcorresponding to the radii R _psubmit themselves to the Laplace's resonance relation.     

The existence of long-lived rays of the coronal streamer belt – Radial density and velocity distributions of the solar wind flowing in them

A technique is proposed for separating the rays of the streamer belt with quasi-stationary and non-stationary solar wind (SW) flows. It is shown that the lifetime of rays with a quasi-stationary SW can exceed 20 days. A new method has been developed for measuring the relative density distribution of a quasi-stationary slow SW flowing along the streamer belt's ray of increased brightness, based on the LASCO/SOHO data. It is shown that the density n for such SW flows varies with the radius R according to the relationship nR ^–, where =₁3.3–3.9 within 4 R ₀ R 6 R ₀ (here R ₀ is the solar radius), and decreases gradually further away. It is also shown that the V(R)-profiles in some rays of the streamer belt differ little from each other, although the value of the mass flow density, j _E, at the Earth's orbit in them can vary more than by a factor of 4. This distinguishes in a crucial respect a slow SW in the streamer belt's rays from a fast SW originating in coronal holes, for which j _Econstant and the dependences V(R) in different fast flows can differ greatly.