Analytical Reference Functions <Emphasis Type="Italic">F</Emphasis>(λ) for the Sun's Limb Darkening and Its Absolute Continuum Intensities (λλ 300 to 1100 m)

It is shown (1) that the coefficients A_i of the limb darkening functions I(μ)/I_center = P₅ (μ) = ∑A_i μⁱ (i = 0... 5; μ = cos ϑ), which had been published by Neckel and Labs (Solar Phys. 153, 91, 1994), can well be approximated by analytical functions of wavelength λ, and (2) that at first sight purely formal extrapolation of the functions P₅(μ) to the very limb (μ = 0.0) is not meaningless: in combination with absolute intensities for the disk center these functions yield ‘limb intensities’ which all correspond to almost the same ‘limb temperature’, T_limb≈4746 K. Together these results lead to ‘reference functions’ which can quickly yield rather reliable values of the Sun's continuum intensities, for any values of μ and λ.     

Nonaxisymmetric instabilities of self-gravitating disks. I toroids

We perform an extensive linear investigation of the nonaxisymmetric disk modes referred to in the literature as P, I, and J modes in self-gravitating polytropic toroids with power law angular velocity distributions. For selected models, we also follow the development of instability from the linear regime through the quasi-linear regime to deep into the nonlinear regime. We consider modes with azimuthal dependence e ^imφ , where m is an integer and φ is the azimuthal angle. We find that instability sets in through m=2 barlike I modes at T/|W|∼0.16–0.18 depending upon the chosen angular velocity law where T is the rotational kinetic energy and W is the gravitational energy of the toroid. Instability in the barlike I mode peaks in strength around T/|W|=0.22–0.23 after which it weakens, eventually stabilizing around T/|W|∼0.25–0.26. One-armed modes (m=1 modes) become unstable just after instability in the m=2I modes sets in; instability in m=1 modes sets in at T/|W|∼0.19. They dominate the barlike I modes in toroids with T/|W|≳0.25. However, almost immediately after the m=1 mode overtakes the barlike I mode, higher-m J modes appear. J modes with m=2, 3, and 4 become unstable for T/|W|≳0.25–0.26, 0.23–0.25, and 0.25–0.26, respectively. m≥3J modes dominate the m=1 mode in toroids with T/|W|≳0.27. As T/|W| increases further, nonaxisymmetric instability sets in through higher and higher m modes. We find quantitative agreement between the early nonlinear behavior of the tested unstable toroids and our linear results. Quasi-linear modeling suggests that a gravitational self-interaction torque which arises early in the nonlinear regime saturates growth of the mode and leads to significant transport of mass and angular momentum. Neither I mode nor J mode instabilities produce prompt fission in toroids.     

Calculation of effective optical depth of absorption line formation in homogeneous semi-infinite planetary atmosphere during anisotropic scattering

The algorithm for determining effective optical thickness of absorption line formation in a plane-parallel homogeneous planetary atmosphere is presented. The case of anisotropic scattering is considered. The results of numerical calculations of τ _e (μ₀) at the scattering angle γ = π for some values of the single scattering albedo λ and the parameter of the Heyney-Greenstein scattering indicatrix g are given. The refined equation for the function T ^m (−μ, μ₀) is presented.     

Zonal structure and meridional drift of large-scale solar magnetic fields

Digitized synoptic charts of photospheric magnetic fields were analyzed for the past 4 incomplete solar activity cycles (1969–2000). The zonal structure and cyclic evolution of large-scale solar magnetic fields were investigated using the calculated values of the radial B _r, |B _r|, meridional B _θ, |B _θ|, and azimuthal B _φ, |B _φ| components of the solar magnetic field averaged over a Carrington rotation (CR). The time–latitude diagrams of all 6 parameters and their correlation analysis clearly reveal a zonal structure and two types of the meridional poleward drift of magnetic fields with the characteristic times of travel from the equator to the poles equal to ∼16–18 and ∼2–3 years. A conclusion is made that we observe two different processes of reorganization of magnetic fields in the Sun that are related to generation of magnetic fields and their subsequent redistribution in the process of emergence from the field generation region to the solar surface. Redistribution is supposed to be caused by some external forces (presumably, by sub-surface plasma flows in the convection zone).     

Does the proton-to-electron mass ratio μ = m p /m e vary in the course of cosmological evolution?

The possible cosmological variation of the proton-to-electron mass ratio μ = m _p /m _e was estimated by measuring the H₂ wavelengths in the high-resolution spectrum of the quasar Q 0347-382. Our analysis yielded an estimate for the possible deviation ofμ value in the past, 10 Gyr ago: for the unweighted valueΔ μ / μ = (3.0±2.4)×10^-5; for the weightedvalueΔ μ / μ = (5.02±1.82)×10^-5.Since the significance of the both results does not exceed3σ, further observations are needed to increase the statistical significance. In any case, this result may be considered as the most stringent estimate on an upper limit of a possible variation of μ (95% C.L.):|Δ μ / μ| < 8× 10^-5 .This value serves as an effective tool for selection of models determining a relation between possible cosmological deviations of the fine-structure constant α and the elementary particle masses (m_p, m_e, etc.). This revised version was published online in July 2006 with corrections to the Cover Date.     

Transition probabilities and dissociation energy of astrophysical molecule CoH

The Franck-Condon (FC) factors (transition probabilities) and r-centroids have been evaluated by the reliable numerical integration procedure for the bands of the A ³ φ₄ →X ³ φ₄ system of astrophysical molecule CoH, using a suitable potential. The dissociation energy D ⁰ ₀ = 2.5 ± 0.05 eV for the electronic ground state of CoH has been estimated by fitting Hulburt-Hirschfelder function to the experimental potential energy curve, using the correlation coefficient. This revised version was published online in July 2006 with corrections to the Cover Date.     

The O III 52 ▪m/88 ▪m emission-line ratio in planetary nebulae

R-matrix calculations of electron impact excitation rates in 0 III are used to derive the electron-density-sensitive emission-line ratioR = 1 (2s ² 2p ²³ P ₂-2s ² 2P ²³ P ₁)/I (2s ²2p ²³ P ₁ - 2s ²2p ²³ P ₀) =I (52μm)/I (88μm) for a range of electron temperatures(Te = 5000-20000 K) and densities(N _e =10–10 ⁵ cm⁻³) applicable to planetary nebulae. Electron densities deduced from the observed values ofR in several planetary nebulae are in excellent agreement with those deduced from C1 I and Ar IV, which provides support for the accuracy of the atomic data adopted in the calculations.     

Millimeter and Submillimeter Emissions from Galactic and Extragalactic Photodissociation Regions

The galactic and Extragalactic photodissociation regions are primarily heated by photoelectrons ejected from the surface of interstellar dust grains by Far-ultraviolet (FUV) photons. But there is no direct mechanism to measure the photoelectric heating efficiency. To understand the role of dust grains in processing the Interstellar Radiation Field (ISRF) and heating the gas, we compare the intensities I _CII, I _CO and I _FIR for (² P _3/2→² P _1/2) & (J = 1→ 0) line emission of CII & CO at 158 μm & 2.6 mm and integrated far-infrared from number of photodissociation regions, HII regions, planetary nebulae, reflection nebulae and high latitude translucent clouds (HLCs). It is found that I _CII is linearly correlated with I _FIR. In the cold medium where cloud is exposed to weak radiations temperature is low and most of the cooling is due to [CII] emissions. As a result the ratio of I _CII/I _FIR provide indirect method to evaluate the photoelectric heating efficiency. For the neutral cold medium it is evaluated to be ∼0.028. The FUV radiation field G ₀ are estimated through the model calculation of I _CII and I _CO for different galactic and photodissociation regions. The intensity of FIR radiation I _FIR are well represented as 1.23×10⁻⁴ G ₀(ergs cm⁻² s⁻¹ sr⁻¹) almost same as estimated for HLCs by Ingalls et al (2002). Hydrogen density for each source has also been estimated.     

A dynamical model of type I supernova atmosphere with the velocity gradient

A compact structure of a low-mass Type I presupernovae is assumed to be an essential feature of the hydrodynamical problem dealing with the supernova Type I (SNI) envelope outbursts. This structure is characterized by a degenerate carbon-oxygen core, which suffers a thermonuclear explosion of carbon fuel (M ₀≃1.40M _⊙), and by a compact lowmass envelope (M _e ≲0.1M _⊙) with external radiusR _e≃10⁹ cm. The parameters, of this hydrostatic envelope are specified and then, for a relatively small explosion energy, ofW ₀≃(2–10)×10⁴⁹ erg, hydrodynamic problem of the envelope ejection is solved numerically. This energy comes from neutrino-induced detonative carbon burning. The resulting structure of the SNI atmosphere expanding with the velocity gradient can be employed for an interpretation of the observed SNI spectra. In accordance with our previous papers, the SNI light curves are considered to occur due to an additional slow (with time-scale 10⁶–10⁷ s) release of the bulk of the SNI energy,W≃10⁵¹, erg. The slow energy release does not, however, affect the structure of the outermost expanding layers of the envelope which are responsible for the SNI spectra. A short (Δt≃10⁻² s) burst of soft (2–10 keV) X-rays with total radiated energy of about 10⁴⁰ erg is found to appear 10–20 days before the SNI optical maximum.     

An extended investigation of Helios 1 and 2 observations: The interplanetary magnetic field between 0.3 and 1 AU

Helios-1 and 2 spacecraft allowed a detailed investigation of the radial dependence of the interplanetary magnetic field components between 0.3 and 1 AU. The behaviour of the radial component B _ris in a very good agreement with Parker's model (B _r r ^-2) and the azimuthal component B also shows a radial dependence which is close to theoretical predictions (B r ^-1). Experimental results for the normal component B and for the field magnitude B are consistent with those from previous investigations. The relative amplitude of the directional fluctuations with periods less than 12 hr is essentially independent of heliocentric distance, while their power decreases approximately as r ^–3 without any appreciable difference between higher and lower velocity regimes.Also at Laboratorio Plasma nello Spazio, CNR, Frascati.     

Robe’s problem: its extension to 2+2 bodies

In the problem of 2+2 bodies in the Robe’s setup, one of the primaries of mass m^*₁m^{*}_{1} is a rigid spherical shell filled with a homogeneous incompressible fluid of density ρ ₁. The second primary is a mass point m ₂ outside the shell. The third and the fourth bodies (of mass m ₃ and m ₄ respectively) are small solid spheres of density ρ ₃ and ρ ₄ respectively inside the shell, with the assumption that the mass and the radius of third and fourth body are infinitesimal. We assume m ₂ is describing a circle around m^*₁m^{*}_{1}. The masses m ₃ and m ₄ mutually attract each other, do not influence the motion of m^*₁m^{*}_{1} and m ₂ but are influenced by them. We also assume masses m ₃ and m ₄ are moving in the plane of motion of mass m ₂. In the paper, the equations of motion, equilibrium solutions, linear stability of m ₃ and m ₄ are analyzed. There are four collinear equilibrium solutions for the given system. The collinear equilibrium solutions are unstable for all values of the mass parameters μ,μ ₃,μ ₄. There exist an infinite number of non collinear equilibrium solutions each for m ₃ and m ₄, lying on circles of radii λ,λ′ respectively (if the densities of m ₃ and m ₄ are different) and the centre at the second primary. These solutions are also unstable for all values of the parameters μ,μ ₃,μ ₄, φ, φ′. Such a model may be useful to study the motion of submarines due to the attraction of earth and moon.     

Dynamical fine structure of a quiescent filament

A quiescent filament was observed near the center of the disk (N5, W5) with the MSDP spectrograph of the 50 cm refractor of the Pic-du-Midi Observatory on June 17, 1986. We focus our study on the statistical moments of the Dopplershift,V ₁, and the intensity,I ₁, at the center of a chord of the Hα profile (±0.256 Å), versus the minimum intensityI ₀. We use a statistical model simulating a numbern _max of threads (of optical thicknessτ ₀ and source functionS ₀), seen over the chromosphere. The threadsj along the same line-of-sighti are identical except for the velocityv _j (gaussian distributionv ₀,σ _v). We search for the best fit between the observed and simulated quantities:V ₁,σ (V ₁),I ₁,σ(I ₁), and the histogram of theI ₀ values over the field of view. A good fit is obtained with: (a) threads characterized byτ ₀ = 0.2,S ₀ = 0.06 (unit of the continuum at disk center), mean upward velocityv ₀ = 1.7 km s⁻¹ and gaussian-type velocity distributionσ _v = 3.5 km s⁻¹. Other possible values ofτ ₀ andσ _v are discussed; (b) underlying chromosphere deduced from observed quiet Sun (outside the filament) by modifying the chromospheric velocities: additional mean upward velocity 0.7 km s⁻¹, standard deviation reduced by a factorF _c ∼ 0.7. The results are discussed in connection with the values deduced from prominence observations.     

Static electrically charged fluids in terms of pressure: general property

A most general exact solution to the Einstein-Maxwell equations for static charged perfect fluid is sought in terms of pressure. Subsequently, metrics (e ^λ and e ^υ ), matter density and electric intensity are expressible in terms of pressure. Consequently, Pressure is found to be an invertible arbitrary function of ω(=c ₁+c ₂ r ²), where c ₁ and c ₂(≠0) are arbitrary constants, and r is the radius of star, i.e. p=p(ω). We present a general solution for charged pressure fluid in terms for ω. We list and discuss some old and new solutions which fall in this category.     

Collision induced rotational excitation of AlF (X 1Σ+) by para-H2 (j=0)

Rotational excitation cross sections and rate coefficients of AlF collisions with para-H₂ are computed at low temperature, i.e., for T≤70 K. Prior to collisional calculations, a four-dimensional (4D) potential energy surface (PES) for the AlF-H₂ system is calculated at the ab initio Coupled-Cluster level of the theory with an aug-cc-pVQZ Gaussian basis set. This 4D-PES is further reduced to a two-dimensional (2D) PES based on the considerations related to collisional studies with para-H₂. The [Al-F] and [H-H] bond lengths are frozen at their experimental equilibrium value r _e =1.654369 bohr and r _e =1.4011 bohr respectively. The interaction energy presents a global minimum located ∼63 cm⁻¹ below the AlF-H₂ dissociation limit. With this PES, cross sections are determined in the Close-Coupling (CC) approach and rate coefficients are inferred by averaging the cross sections over a Maxwell-Boltzman distribution of kinetic energies. These quantities are significantly magnified in comparison with their AlF-He counterparts. The already observed propensity towards ΔJ=1 transitions for AlF-He remains.     

The intensity,velocity and magnetic structure of a sunspot region

High resolution sunspot photographs in the blue, red and infrared continuum exposed on various days were used to derive the center-to-limb variation of the intensity ratio = I _sp / I _ph. Special care was taken to correct for image blurring, scattered light and the influence of line absorption.The observed increase of specific umbral intensities _u towards the limb leads to an extremely small temperature gradient in the umbra. From geometrical changes of the profiles (Wilson effect), we derived an umbral depression of about 650 km and a density scale height of about 450 km when H ^- is assumed to be the predominant source of absorption. The penumbral depression was found to be 50 km or less. The density scale height of the umbra as computed from the observed temperature distribution is 80 km in the case of hydrostatic equilibrium. We conclude that either magnetic pressure components produce deviations from hydrostatic equilibrium or that another source of absorption, dominating in the outer layers, has to be taken into account.     

On bridges B(pL, qS) around triangular libration points

When μ is smaller than Routh’s critical value μ ₁ = 0.03852 . . . , two planar Lyapunov families around triangular libration points exist, with the names of long and short period families. There are periodic families which we call bridges connecting these two Lyapunov families. With μ increasing from 0 to 1, how these bridges evolve was studied. The interval (0,1) was divided into six subintervals (0, μ ₅), (μ ₅, μ ₄), (μ ₄, μ ₃), (μ ₃, μ ₂), (μ ₂, μ ₁), (μ ₁, 1), and in each subinterval the families B(pL, qS) were studied, along with the families B(qS, qS′). Especially in the interval (μ ₂, μ ₁), the conclusion that the bridges B(qS, qS′) do not exist was obtained. Connections between the short period family and the bridges B(kS, (k + 1)S) were also studied. With these studies, the structure of the web of periodic families around triangular libration points was enriched.     

Om diagnostic of dilaton dark energy with two-parameter ω in potential [1+(σ−A)2]e (−Bσ)

Based on a new geometric diagnostic method-Om, we consider a new independent-model parametrization . When we work in potential W _σ [1+(σ−A)²]e ^(−Bσ), we investigate the evolutional behavior of Om with respect to red-shift z and the influence of coupling parameter α on the trajectory of Om with respect to z. We get that phantom model of Dilaton dark energy can avoid the future singularity “Big Rip”. The numerical results give current value of EOS which fits the latest observational data WMAP5+BAO+SNe very well.     

Variety of well behaved exact solutions of Einstein–Maxwell field equations: an application to Strange Quark stars, Neutron stars and Pulsars

We present a variety of well behaved classes of Charge Analogues of Tolman’s iv (1939). These solutions describe charged fluid balls with positively finite central pressure, positively finite central density; their ratio is less than one and causality condition is obeyed at the centre. The outmarch of pressure, density, pressure-density ratio and the adiabatic speed of sound is monotonically decreasing, however, the electric intensity is monotonically increasing in nature. These solutions give us wide range of parameter for every positive value of n for which the solution is well behaved hence, suitable for modeling of super dense stars. keeping in view of well behaved nature of these solutions, one new class of solutions is being studied extensively. Moreover, this class of solutions gives us wide range of constant K (0.3≤K≤0.91) for which the solution is well behaved hence, suitable for modeling of super dense stars like Strange Quark stars, Neutron stars and Pulsars. For this class of solutions the mass of a star is maximized with all degree of suitability, compatible with Quark stars, Neutron stars and Pulsars. By assuming the surface density ρ _b =2×10¹⁴ g/cm³ (like, Brecher and Caporaso in Nature 259:377, 1976), corresponding to K=0.30 with X=0.39, the resulting well behaved model has the mass M=2.12M _Θ, radius r _b ≈15.27 km and moment of inertia I=4.482×10⁴⁵ g cm²; for K=0.4 with X=0.31, the resulting well behaved model has the mass M=1.80M _Θ, radius r _b ≈14.65 km and moment of inertia I=3.454×10⁴⁵ g cm²; and corresponding to K=0.91 with X=0.135, the resulting well behaved model has the mass M=0.83M _Θ, radius r _b ≈11.84 km and moment of inertia I=0.991×10⁴⁵ g cm². For n=0 we rediscovered Pant et al. (in Astrophys. Space Sci. 333:161, 2011b) well behaved solution. These values of masses and moment of inertia are found to be consistent with other models of Neutron stars and Pulsars available in the literature and are applicable for the Crab and the Vela Pulsars.     

Accelerated expansion from?a?modified-quadratic gravity

We investigate the late-time dynamics of a four-dimensional universe based on modified scalar field gravity in which the standard Einstein-Hilbert action R is replaced by f(φ)R+f(R) where f(φ)=φ ² and f(R)=AR ²+BR _μν R ^μν,(A,B)∈ℝ. We discussed two independent cases: in the first model, the scalar field potential is quartic and for this special form it was shown that the universe is dominated by dark energy with equation of state parameter w≈−0.2 and is accelerated in time with a scale factor evolving like a(t)∝t ^5/3 and B+3A≈0.036. When, B+3A→∞ which corresponds for the purely quadratic theory, the scale factor evolves like a(t)∝t ^1/2 whereas when B+3A→0 which corresponds for the purely scalar tensor theory we found when a(t)∝t ^1.98. In the second model, we choose an exponential potential and we conjecture that the scalar curvature and the Hubble parameter vary respectively like R=hH[(f)\dot]/f,h ? \mathbbRR=\eta H\dot{\phi}/\phi,\eta\in\mathbb{R} and H=g[(f)\dot]^c,(g,c) ? \mathbbRH=\gamma\dot{\phi}^{\chi},(\gamma,\chi)\in\mathbb{R}. It was shown that for some special values of  χ, the universe is free from the initial singularity, accelerated in time, dominated by dark or phantom energy whereas the model is independent of the quadratic gravity corrections. Additional consequences are discussed.     

Doubly averaged effect of the Moon and Sun on a high altitude Earth satellite orbit

Infinite series expansions are obtained for the doubly averaged effects of the Moon and Sun on a high altitude Earth satellite, and the results used to interpret numerically integrated examples. New in this paper are: (1) both sublunar and translunar satellites are considered; (2) analytic expansions include all powers in the satellite and perturbing body semi-major axes; (3) the fact that retrograde orbits have more benign eccentricity behavior than direct orbits should be exploited for high altitude satellite systems; and (4) near circular orbits can be maintained with small expenditures of fuel in the face of an exponential driving force one forI _ab, whereI _b=180°–I _a andI _a is somewhat less than 39.2° for sublunar orbits and somewhat greater than 39.2° for translunar orbits.Nomenclature a semi-major axis - A _lk coefficient defined in Equation (11) - B _lk coefficient defined in Equation (24) - C _km coefficient defined in Equation (25) - D, E, F coefficients in Equations (38), (39) - e eccentricity - H _k expression defined in Equation (34) - expression defined in Equation (35) - I inclination of satellite orbit on lunar (or solar) ring plane - J ₂ coefficient of second harmonic of Earth's gravitational potential (1082.637×10^–6 R _E ² ) - K _k, L_k, M_k expressions in Section 4 - expressions in Section 4 - p=a(1–e ²) semi-latus rectum - P _l Legendre polynomial of degreel - q argument of Legendre polynomial - radial distance of satellite - R _E Earth equatorial radius (6378.16 km) - R, S, W perturbing accelerations in the radial, tangential and orbit normal directions - syn synchronous orbit radius (42 164.2 km=6.6107R _E) - t time - T satellite orbital period - T orbital period of perturbing body (Moon) - T _e period of long periodic oscillations ine for |I|<I _a - T _s synodic period - U gravitational potential of lunar (or solar) ring - x, y, z Cartesian coordinates of a satellite with (x, y) being the ring plane - coefficient defined in Equation (20) - average change in orbital element over one orbit (=a, e, I, , ) - ₁,₂₃ unit vectors in thex, y, z coordinate directions - _r , _s , _w unit vectors in the radial, tangential and orbit normal directions - =+ angle along the orbital plane from the ascending node on the ring plane to the true position of the satellite - angle around the ring - gravitational constant times mass of Earth (3.986 013×10⁵ km s^–2) - gravitational constant times mass of Moon (or Sun) - _m gravitational constant times mass of Moon (/81.301) - _s gravitational constant time mass of Sun (332 946 ) - ratio of the circumference of a circle to its diameter - radius of lunar (or solar) ring - _m radius of lunar ring (60.2665R _E) - _s radius of solar ring (23455R _E) - true anomaly - argument of perigee - ₀ initial value of - _i critical value of in quadranti(i=1, 2, 3, 4) - longitude of ascending node on ring plane This work was sponsored by the Department of the Air Force.