Effects of metal abundances on the evolutionary period changes of classical Cepheids

Some peculiarities in the behaviour of a model self-gravitating system described by hydrodynamical equations and isothermal equation of state connected with the presence of thermodynamical fluctuations in real systems were investigated in numerical experiment. The values of density and velocity , , respectively, were computed by numerical code perturbed on each time-step and in each computational cell by random values , for modeling such fluctuations. Perturbed values _i = _i + _i ,v _i = _i + v _i were used to initiate the next step of computations. This procedure is equivalent to an introduction into original hydrodynamical equations of Langevin sources which are random functions. It is shown that these small fluctuations (= v =0,² =v ² = 10^–8) grow many times in marginally-stable state.     

О лагранжевых и зйлеровых решеиях обощенной задачи трех тел

(, 1969). ( ), ( ), , , . , (=), , , .. , . , , - ( ), ( ). , .

---

This paper is a continuation and a generalization of one published earlier (Duboshin, 1969): it discusses the problem whether there exist the Lagrangian and the Eulerian solutions of the generalized three-body (material points) problem. Every point in this generalized problem acts on another, one with a force (attractive or repulsive) directed along the straight line passing through these points, and in an arbitrary manner depending on time, mutual distance and its derivatives, the first and the second. Here, generally speaking, the third axiom of dynamics (law of action and reaction) is not presupposed as fulfilled, that is, it is supposed that every two material points interact in a different way.This most general assumption being made, we establish the conditions which must dictate the laws of the interactions, so that the three points can always remain at the apexes of the equilateral triangle (Langrangian solution), or remain always on a straight line (Eulerian solution).The author believes that such general treatment of the three-body problem can be useful for theoretical studies in celestial mechanics and also for practical applications in the study of isolated stellar systems.

     

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.     

The Harvard-Smithsonian reference atmosphere

We present a model of the solar atmosphere in the optical depth range from ₅₀₀₀ = 10^–8 to 25. It combines an improved model of the photosphere that incorporates recent EUV observations with a new model of the quiet lower chromosphere. The latter is based on OSO 4 observations of the Lyman continuum, on infrared observations, and on eclipse electron densities.Our model differs from the Bilderberg Continuum Atmosphere (BCA) in the low chromosphere (₅₀₀₀ < 10^–4), where deviations from local thermodynamic equilibrium in hydrogen and carbon have been taken into account. It also differs in the transition region between the chromosphere and the photosphere (10^–4 < ₅₀₀₀ < 10^–2), where the temperature is lower than in the BCA, and in the convective region (₅₀₀₀ 2), where the temperature is higher than in the BCA.     

High order symplectic integrators for perturbed Hamiltonian systems

A family of symplectic integrators adapted for the integration of perturbed Hamiltonian systems of the form H=A+B was given in (McLachlan, 1995). We give here a constructive proof that for all integer p, such integrator exists, with only positive steps, and with a remainder of order O(^p + ²²), where is the stepsize of the integrator. Moreover, we compute the analytical expressions of the leading terms of the remainders at all orders. We show also that for a large class of systems, a corrector step can be performed such that the remainder becomes O(^p +⁴²). The performances of these integrators are compared for the simple pendulum and the planetary three-body problem of Sun–Jupiter–Saturn.     

Оптимальное демпфирование нутационного движения спутников, стаБилизируемых вращением

, . , , . , , -, This article investigates the dynamics of a system for damping the nutating motion of a spin-stabilized satellite. The equations of motion of the satellite-damper system are derived omitting consideration of the influence of external torques. The conditions of stability of the stationary spinning are obtained and the optimal parameters of the satellite and the damper ensuring a maximal rate of damping of the nutation motion are determined.     

Structure and physics of solar faculae

Taking into account the effect of roughness (or local departures from sphericity) of the emitting layers in the chromosphere-corona transition zone (CCT) allows one to determine the optical depths of layers responsible for resolved structures in Cii, Ciii, Oiv, and Ovi lines. The result, at the top of the irregularities, is of the order of respectively ₁ 3.5,2.0, 1.6,0.5, and for the bottom of these irregularities, ₂ = 0.7, 0.4, 0.3, 0.25. The characteristic angle of these irregularities is, respectively, of the order of 35 °, 33 °, 35 °, and 41 °. For unresolved structures of Civ and Ovi (already analyzed in the spherical symmetry hypothesis in Paper III), one finds ₁ 0.6; 0.9 and ₂ 0.12; 0.2 in the case of quiet areas; in the case of active areas, the range is broader for Civ and Ovi, from 1.0 to 1.7 for ₁ and from 0.2 to 0.9 for ₂. The values obtained from Ovi are in reasonable agreement with each other for resolved and unresolved structures. And the obtained values of ₁ and ₂ correspond not too badly with the determinations made in Paper III, by methods not exceedingly influenced by the spherical symmetry hypothesis.     

Parameters of Dust Particles in the Martian Atmosphere

A critical analysis of the methods and results of estimating the optical thickness of the dust component in the Martian atmosphere ₀, the particle size r ₀, and the imaginary part of the refractive index n _ihas shown the following. (1) Observational data on the brightness distribution over the Martian disk as well as the phase dependences of diffusely reflected light and the azimuthal dependences of diffusely transmitted light are most appropriate to use only for verifying the reliability of the aerosol parameters determined by other methods. (2) If the morning and evening fogs in the atmosphere are disregarded, the Bouguer–Lambert–Beer method used to analyze the solar-brightness attenuation measured on the planetary surface yields overestimated extraatmospheric solar intensity I ₀and atmospheric optical depth ₀. At the Viking 1landing site, I ₀and ₀could be overestimated by a factor of 1.7 and by 0.35, respectively. (3) The aerosol size determined by analyzing measurements of the azimuthal dependences for the Martian sky brightness at low elevations of the Sun most likely corresponds to the fog particles. (4) If overestimated values of I ₀were used to standardize the observations of the solar radiation transmitted by the Martian atmosphere, then n _iwere also overestimated; using overestimated ₀also affected the reliability of the latter. (5) The problem of reliability of the available ₀and r ₀estimates for periods of high atmospheric transparency is yet to be solved. For the highest activity of the dust storm in 1971, it was found that 4.5 r ₀ 7.5 m for the lognormal particle size distribution with ²= 0.2 and the optical thickness of a dust cloud ₀ 15. (6) The spectral values of the apparent albedo of Mars measured in October 1971 at a phase angle of 42° in the spectral range 0.250 0.717 allowed the imaginary part of the refractive index to be estimated in terms of a model of a dust cloud composed of spherical particles with the lognormal size distribution with r ₀= 4.5 m and ²= 0.2.     

On the solar oblateness: The combined effect of a pole-equator difference in effective temperature and mechanical heating

With the help of a model atmosphere of the Sun we evaluate the pole-equator difference in flux (as measured by Dicke and Goldenberg) assuming the following type of pole-equator temperature difference (T=T _e –T _p ): (a) T 2K for > ₀ (₀ 0.05); (b) T 10K for < ₀.The small T at all optical depths given by (a) could, for example, be due to a pole-equator difference in effective temperatures. At small optical depths a difference in mechanical heating could give rise to the larger temperature difference given by (b). We compare the results of our calculations with Dicke and Goldenberg's observations.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.     

Radiative damping of trapped gravity waves in the solar atmosphere

The radiative damping of trapped gravity waves in an optically thin atmosphere is studied for a stratified Boussinesq fluid. The character of the atmospheric eigenmodes depends on the distribution of the Brunt-Väisälä frequency N and the radiative relaxation time . The calculations for simple layer models show that if N is large over some finite fraction of the trapping region, then modes of long lifetime can exist. In order to suppress gravity waves entirely, it is necessary that N < 1 over the entire trapping region. Qualitative application of the results to the solar atmosphere leads to the conclusion that gravity wave eigenmodes of the solar atmosphere, although damped, are by no means eliminated by radiative effects.     

Physical Properties of Dust in the Martian Atmosphere: Analysis of Contradictions and Possible Ways of Their Resolution

Atmospheric aerosols play an important role in forming the Martian climate. However, the basic physical properties of the Martian aerosols are still poorly known; there are many contradictions in their estimates. We present an analytical overview of the published results and potentialities of various methods. We consider mineral dust. Zonally averaged data obtained from mapping IR instruments (TES and IRTM) give the optical thickness of mineral aerosols ₉ = 0.05–0.1 in the 9-m band for quite atmospheric conditions. There is a problem of comparing these estimates with those obtained in the visible spectral range. We suggest that the commonly used ratio _vis/₉ >2 depends on the interpretation and it may actually be smaller. The ratio _vis/₉ 1 is in better agreement with the IRIS data (materials like montmorillonite). If we assume that _vis/₉ = 1 and take into account the nonspherical particle shape, then the interpretation of ground-based integrated polarimetric observations ( < 0.04) can be reconciled with IR measurements from the orbit. However, for thin layers, the sensitivity of both methods to the optical thickness is poorly understood: on the one hand, polarimetry depends on the cloud cover and, on the other hand, the interpretation of IR measurements requires that the atmospheric temperature profile and the surface temperature and emissivity be precisely known. For quite atmospheric conditions, the local optical-thickness estimates obtained by the Bouguer–Lambert–Beer method and from the sky brightness measured from Viking 1 and 2 and Mars Pathfinder landers are much larger: = 0.3–0.6. Estimates of the contrasts in images from theVikingorbiters yield the same values. Thus, there is still a factor of 3 to 10 difference between different groups of optical-thickness estimates for the quiet atmosphere. This difference is probably explained by the contribution of condensation clouds and/or by local/time variations.     

Particle lifetimes in strong diffusion

The mean lifetime of a particle distribution, driven to isotropy by intense pitch-angle diffusion, is calculated by analytical means for conditions applicable to the Earth's magnetosphere. The resulting algebraic expressions reduce to [64La/35v _c ² (1–)] in the limit of a small equatorial loss cone (half-angle _c ), wherev is the particle speed,L is the magnetic shell parameter,a is the radius of the Earth, and is the particle albedo from the atmosphere at either foot of the field line. Distinction is made in the full expressions for between complete isotropy (caused by strong pitch-angle diffusion all along the field line) and incomplete isotropy (caused by strong diffusion that is localized at the magnetic equator) over the upward hemisphere in velocity space.     

Optically thick accretion onto black holes

Spherically symmetric, steady-state, optically thick accretion onto a nonrotating black hole with the mass of is studied. The gas accreting onto the black hole is assumed to be a fully ionized hydrogen plasma withn ₀=10⁸ cm^–3 andT ₀=10⁴ K far from the black hole, and a new approximate expression for the Eddington factor is introduced. The luminosity is estimated to beL=1.875×10³³ erg s^–1, which primarily arises from the optical surface (1) ofT10⁴ K. The accretion flow is characterized by 1 and (v/c)10. In the optically thin region, the flow remains isothermal, and the increase of temperature occurs at 1. The radiative equilibrium is strictly realized at (v/c)10.     

Les raies du triplet de l'helium de la Nova Delphini 1967

The aim of this paper has been to study the neutral helium triplet emission lines identified in the spectrum of the envelope of Nova Delphini. By comparing the observed flux of the neutral helium lines with that calculated theoretically by Robbins, we find that the optical thickness in the center of the line 3889 is of the order of 21.50 for summer 1969. The optical thickness obtained by this method is here denoted _tran(3889).On the other hand, we obtain the number of neutral helium atoms in the 2³S state [N(2³S)] by considering the equilibrium between the mechanisms that populate and depopulate this state. We then find that the depopulation by photoionization due to the radiation of Ly (Hi), transitions to the 2¹S, 2¹p and 2³p states by electron collision, photoionization due to the continuum radiation of the central star, are 82.70%, 13.20%, 2.40%, 0.90% and 0.80% respectively. We find that the mechanism of the photoionization by Ly is the dominant mechanism of depopulation of 2³S state. We calculated ( 3889) of the order of 82.37, fromN(2³S), obtained in the preceding paragraph. The optical thickness obtained by this method is here denoted _bal(3889).The difference between _tran(3889) and _bal(3889) is very large and it cannot be attributed to calculation errors. We have considered all the mechanisms that can depopulate the 2³S state, so we then conclude that the difference between _tran(3889) and _bal(3889) is due to the heterogeneity of the envelope of the Nova, already found by us in our previous study of the profiles of the permitted and forbidden lines. Finally, we find from this study a filling factor of the order of 0.30.     

Dimensional Analysis of the Flare Distribution Problem

Dimensional analysis is used to derive the distribution of solar flare energies,p() = A-^3/2, in accordance with recent observational and numerical results. Several other scalings, notably _fl ² , where _fl is the flare duration, are obtained as well.     

Stratification of the 5-min oscillation through the solar photosphere

The inversion of the radiative transfer equation has been applied to a time series of low spatial resolution Ki 7699 Å line spectra, allowing us to obtain a sequence of instantaneous models which are compared to the mean one to derive the 5-min oscillation stratification. The resulting models are reliable from log 0 to log -3.5. The increase of the velocity oscillation amplitude towards higher photospheric layers has been confirmed and quantified. Significant oscillation is observed in temperature.     

A simplified model solar atmosphere

A simplified representation of the temperature distribution in the solar photosphere is proposed: ( ₀) = ₀ - ₁ log ₀. An expression is derived for the emergent continuous spectrum from the simple model. The limitations and applications of the simple model are discussed.     

On the A-transform of the exponential integrals

In this paper, a technique of recursive analysis is developed for the integral transform A of the exponential integral functionsE _n which is denoted as _n (). The main result of this analysis enables us to establish a two-term recurrence formula for _n (0) and a three-term recurrence formula for _n (); 0. A computational algorithm based on these formulae is also constructed and its numerical results forn=2(1)25 are presented to 15-digit accuracy.     

The galactic diffuse component of soft X-rays and its source function

A new model for the source distribution of galactic soft X-ray (B and C band) emission is presented. From the mean dependence of count rates on galactic latitudeb (i.e., the brightness distribution), we derive the soft X-ray source functionQ as function of the optical depth by solving the equation of radiative transfer with the aid of a Laplace transform. Contrary to older Heaviside step models,Q is found to increase strongly, but not abruptly, in the range 1.5<<2.5, indicating a noticeable emission of X-rays from beyond theHi scale height. Using standard X-ray absorption cross-sections for the interstellar medium, we find that the B band X-ray emission coefficient is non-zero within theHi disk and has a maximum at az-value slightly above the Hi scale height. In the C band, the emission coefficient slightly decreases with increasingz, at least up to theHi scale height. A non-zero source function near the galactic plane implies that the interstellar medium (ISM) within theHi scale height is not only an absorbing layer but is mixed with X-ray emitting regions. The so-called local hot bubble is adopted as one of these regions. The maximum of the B band emission coefficient, together with the sharp increase ofQ, is strong evidence for the existence of a galactic soft X-ray halo, and, moreover, give rise to the assumption of a general intergalactic X-ray background. The effective absorption cross-sections given in the literature, based on an (pure) exponential dependence in the negative correlation between count rates andHi column densities, were biased to be too small, in particular in the B band. In replacing the Heaviside step (in the ISM) by a smoothed transition region, these inconsistencies become spurious.