Interstellar Extinction and the Intrinsic Colors of Classical Cepheids in the Galaxy, the LMC, and the SMC

New methods are applied to samples of classical cepheids in the galaxy, the Large Magellanic Cloud, and the Small Magellanic Cloud to determine the interstellar extinction law for the classical cepheids, R _B:R _V:R _I:R _J:R _H:R _K= 4.190:3.190:1.884:0.851:0.501:0.303, the color excesses for classical cepheids in the galaxy, E(B-V)=-0.382-0.168logP+0.766(V-I), and the color excesses for classical cepheids in the LMC and SMC, E(B-V)=-0.374-0.166logP+0.766(V-I). The dependence of the intrinsic color (B-V)₀ on the metallicity of classical cepheids is discussed. The intrinsic color (V-I)₀ is found to be absolutely independent of the metallicity of classical cepheids. A high precision formula is obtained for calculating the intrinsic colors of classical cepheids in the galaxy: (<B>-<V>)₀=0.365(±0.011)+0.328(±0.012)logP.     

Stability of motion in the Sitnikov 3-body problem

We study the stability of motion in the 3-body Sitnikov problem, with the two equal mass primaries (m ₁ = m ₂ = 0.5) rotating in the x, y plane and vary the mass of the third particle, 0 ≤ m ₃ < 10⁻³, placed initially on the z-axis. We begin by finding for the restricted problem (with m ₃ = 0) an apparently infinite sequence of stability intervals on the z-axis, whose width grows and tends to a fixed non-zero value, as we move away from z = 0. We then estimate the extent of “islands” of bounded motion in x, y, z space about these intervals and show that it also increases as |z| grows. Turning to the so-called extended Sitnikov problem, where the third particle moves only along the z-axis, we find that, as m ₃ increases, the domain of allowed motion grows significantly and chaotic regions in phase space appear through a series of saddle-node bifurcations. Finally, we concentrate on the general 3-body problem and demonstrate that, for very small masses, m ₃ ≈ 10⁻⁶, the “islands” of bounded motion about the z-axis stability intervals are larger than the ones for m ₃ = 0. Furthermore, as m ₃ increases, it is the regions of bounded motion closest to z = 0 that disappear first, while the ones further away “disperse” at larger m ₃ values, thus providing further evidence of an increasing stability of the motion away from the plane of the two primaries, as observed in the m ₃ = 0 case.     

Integrability of the Yang-Mills Hamiltonian system

This paper considers the integrability of generalized Yang-Mills system with the HamiltonianH _a (p, q)=1/2(p ₁ ² +p ₂ ² +a ₁ q ₁ ² +a ₂ q ₂ ² )+1/4q ₁ ⁴ +1/4a ₃ q ₂ ⁴ + 1/2a ₄ q ₁ ² q ₂ ² . We prove that the system is integrable for the cases: (A)a ₁=a ₂,a ₃=a ₄=1; (b)a ₁=a ₂,a ₃=1,a ₄=3; (C)a ₁=a ₂/4,a ₃=16,a ₄=6. Our main result is the presentation of these integrals. Only for cases A and B does the Yang-Mills Hamiltonian possess the Painlevé property. Therefore the Painlevé test does not take account of the integrability for the case C.     

Improved theoretical line ratios for Ciii in the Sun

The recent twelve-state R-matrix calculations of electron excitation rates in Ciii by Berrington are used to derive level populations applicable to the solar transition region. Line ratios R = I(2p ^{2 3} P ^e - 2s2p ³ P ^°)/I(2s2p ¹ P ^° - 2s ^{2 1} S ^e ) and R ₂=I(2p ^{2 1} S ^e - 2s2p ¹ P ^°)/I(2p ^{2 3} P ^e - 2s2p ³ P ^°) deduced from these data in conjunction with the relevent transition probabilities are found to be in much better agreement with the observed quiet Sun values than those determined from the level population calculations of Keenan et al.     

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.     

The ratios I(K1)/I(H1) and I(K3)/I(H3) of ca II as diagnostics of the chromosphere above sunspots

The ratios I(K ₁)/I(H ₁) and I(K ₃)/I(H ₃) were calculated from four semi-empirical models of sunspot umbra. We determined the dependencies of both ratios of such parameters as temperature gradient and atmospheric opacity. A certain influence on the expected ratios I(K ₁)/I(H ₁) and I(K ₃)/I(H ₃) can also come from the FIP effect provided it exists in the chromosphere above sunspot umbra. Theoretical and observed values of I(K ₁)/I(H ₁) and I(K ₃)/I(H ₃) are compared. It is shown that for one of the sunspots we observed, the values obtained for the ratio I(K ₁)/I(H ₁) cannot be explained in terms of existing umbra models.     

Note on the Generalized Hansen and Laplace Coefficients

Recently, Breiter et al. [Celest. Mech. Dyn. Astron., 2004, 88, 153–161] reported the computation of Hansen coefficients X _k ^{γ ,m} for non-integer values of γ. In fact, the Hansen coefficients are closely related to the Laplace b _s ^(m), and generalized Laplace coefficients b _s,r ^(m) [Laskar and Robutel, 1995, Celest. Mech. Dyn. Astron., 62, 193–217] that do not require s,r to be integers. In particular, the coefficients X ₀ ^{γ ,m} have very simple expressions in terms of the usual Laplace coefficients b _{γ +2} ^(m), and all their properties derive easily from the known properties of the Laplace coefficients.     

Presence of LaO, ScO and VO Molecular Lines in Sunspot Umbral Spectra

High-resolution Fourier Transform Spectrometer sunspot umbral spectra of the National Solar Observatory/National Optical Astronomy Observatory at Kitt Peak were used to detect rotational lines from 19 electronic transition bands of the molecules LaO, ScO and VO, in the wavenumber range of 11 775 to 20 600 cm⁻¹. The presence of lines from the following transitions is confirmed: A ² Π _r1/2 – X ² Σ ⁺(0, 0; 0, 1), A ² Π _r3/2 – X ² Σ ⁺(1, 0), B ² Σ ⁺ – X ² Σ ⁺(0, 0; 0, 1; 1, 0) and C ² Π _r1/2 – A′²Δ _r3/2(0, 0; 1, 1) of LaO; A ² Π _r3/2 – X ² Σ ⁺(0, 0), A ² Π _r1/2 – X ² Σ ⁺(0, 0) and B ² Σ ⁺ – X ² Σ ⁺(0, 0) of ScO; and C ⁴ Σ ⁻ – X ⁴ Σ ⁻(0, 1; 1, 0; 0, 2) and (2, 0) of VO. However, the presence of A ² Π _r3/2 – X ² Σ ⁺(0, 0) and C ² Π _r3/2 – A′²Δ _r5/2(0, 0; 1, 1) of LaO and C ⁴ Σ ⁻ – X ⁴ Σ ⁻(0, 0) of VO are found to be doubtful because the lines are very weak, and detections are difficult owing to heavy blending by strong rotational lines of other molecules. Equivalent widths are measured for well-resolved lines and, thereby, the effective rotational temperatures are estimated for the systems for which the presence is confirmed.     

Meteoroid Stream Searching: The Use of the Vectorial Elements

In this initial study, we propose a new distance function D _V involving heliocentric vectorial orbital elements. The function measures differences between: the orbital energies, the angular momentums vectors and the Laplace vectors. In comparison with the widely used D _SH criterion of Southworth and Hawkins, D _D criterion of Drummond and their hybrid D _H by Jopek, the new function contains one invariant with respect to the principal secular perturbation: the orbital energy. The new function proved to be useful in the classification amongst the IAU2003 meteoroids which we searched for streams by D _V function and also using D _SH and D _N -function given by Valsecchi et al. For major streams, the results agree very well. For minor, and near-ecliptical streams the results sometimes differ markedly.     

Theoretical emission line ratios for Si XIII compared to solar observations

The electron collision excitation rates recently calculated for transitions in Si xiii by Keenan et al. (1987) are used to derive the electron temperature sensitive ratio G(=(f + i)/r and the density sensitive ratio R(=f/i), where i, f, and r are the intercombination (1s ^{2 1} S – 1s2p ³ P _{1, 2}) forbidden (1s ^{2 1} S – 1s2s ³ S), and resonance (1s ^{2 1} S – 1s2p ¹ P), transitions respectively. Also estimated are the values of R in the low-density limit (R ₀) as a function of electron temperature. The theoretical G ratio at the temperature of maximum emissivity for Si xiii, G(T _m) = 0.70, is in much better agreement with the observed G for the 1985, May 5 flare determined by McKenzie et al. (G = 0.60 ± 0.07) than is the earlier calculation of Pradhan, who derived G(T _m) = 0.85. The error in the observed R ₀ ratio is so large that both our result and Pradhan's fall within the acceptable limits of uncertainty and hence one cannot estimate which of the two is the more accurate.     

A mass ratio limit for primordial black holes

We augment our scenario for the formation of astronomical objects from macroscopic superstrings by the assumption that the central matter keeps its identity in the fragmentation. From the condition that the angular momentum per mass squared of this matter should be less than the Kerr limit G/c, we obtain upper limits for the ratio of the mass of central black holes M(BH) to the mass M of the host object. This limit is M(BH)/M ≈ 0.001, and, expressed in observed quantities, approximately M(BH)/M ≈ σ²/(v · c) where σ is the r.m.s. velocity, v the rotational velocity and c the velocity of light. The valuesM(BH) agree with the observed behaviour both in order of magnitude and in the variation with velocity dispersion. (© 2003 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

MGS electron density profiles: Analysis of the peak magnitudes

We analyze here the behavior of the magnitudes of the F₁ and E peaks of the electron density profiles measured by the Radio Science Subsystem of the Mars Global Surveyor spacecraft, as a function of solar zenith angle χ and solar flux. For each of the 658 days of data in the six occultation seasons in the northern hemisphere, we choose one profile to analyze, which is that for which the F₁ peak is the median value. We assume that the variations of the measured peak densities can be represented as Aa(cosχ) and as Bb(F_10.7), where F_10.7 is the usual solar flux proxy, appropriately shifted to the orbital position of Mars. To minimize the effect of solar activity, we divide the data into 6 F_10.7 bins, fit the data in each bin, and derive the values of the exponent a and the coefficient AF10.7 for each bin. The median values that we derive for the exponent a is 0.46 for the F₁ peak, and 0.395 for the E peak. To minimize the effect of SZA, we divide the data into eight SZA bins, and derive the exponent b and the coefficient Bχ for each SZA bin. We argue that the last three SZA bins should be excluded because the fits were poor, due partly to the small number of data points in each of these bins. If we do so, the median values of b that we derive are 0.27 and 0.40 for the F₁ and E peaks, respectively. Finally we derive a 3-parameter fit to all the data, which expresses the variability of the peak densities as a function of a^′(cosχ) and b^′(F_10.7) simultaneously. The fitted values of the exponents a^′ and b^′ for the F₁ peak are 0.45 and 0.26, respectively; for the E peak, the values are 0.39 and 0.46, respectively. We compare our results to Chapman theory, and to those of other investigators.     

Families of orbits in planar anisotropic fields

The aim of the planar inverse problem of dynamics is: given a monoparametric family of curves f(x, y) = c, find the potential V (x, y) under whose action a material point of unit mass can describe the curves of the family. In this study we look for V in the class of the anisotropic potentials V(x, y) = v(a²x² + y²), (a=constant). These potentials have been used lately in the search of connections between classical, quantum, and relativistic mechanics. We establish a general condition which must be satisfied by all the families produced by an anisotropic potential. We treat special cases regarding the families (e. g. families traced isoenergetically) and we present certain pertinent examples of compatible pairs of families of curves and anisotropic potentials. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Density and temperature diagnostics of solar emission lines from NeV/MgV and SiVII/MgVII Ions

We present NeV/MgV and SiVII/MgVII theoretical line intensity ratios as a function of electron densityN _e and temperatureT _e . These are shown in the form of ratio-ratio diagrams, which should in principle allow bothN _e andT _e to be deduced for the emitting region of the solar plasma. We apply these diagnostics in the solar atmosphere, and discuss the available observations made from space. In most cases, however, we deduceN _e andT _e from the computed absolute line intensities in a spherically symmetric model atmosphere of the Sun. Possible future applications of this investigation to spectral data from the Coronal Diagnostic Spectrometer (CDS) on the Solar and Heliospheric Observatory (SOHO) are briefly discussed.     

Tilted Bianchi Type I Cosmological Models for Barotropic Perfect Fluid in General Relativity

In this paper, we have investigated that tilted Bianchi Type I cosmological models for stiff perfect fluid under a supplementary condition A = B ⁿ between metric potentials, is not possible. The tilted solution is also not possible when we assume A = t ^ℓ, B = t ^m , C = t ⁿ ; ℓ, m and n are constants for ε = p. Thus to preserve tilted nature of model, we assume p = γε, 0 ≤ γ ≤ 1 (barotropic equation of state) for the case A = t ^ℓ B = t ^m and C = t ⁿ . The physical and geometrical aspects of the models are also discussed. This revised version was published online in July 2006 with corrections to the Cover Date.     

Kukarkin-Parenago relationship,period gap,and cause of dwarf nova outbursts

It is shown that it can be useful to incorporate the orbital period P as an additional parameter into the Kukarkin-Parenago relationship. Then it turns out that the relationship splits into two, one holding for objects above and the other for objects below the well-known and much-discussed period gap between 2.2 and 2.8 hours A = -2.03 + 2.26 log C + 8.80 log P for P > 2.2 h A = 2.07 + 2.46 log C - 3.68 log P for P > 2.8 h (A amplitude (mag), C cycle length (d), P orbital period (h)). Future theories of the dwarf nova phenomenon and the origin of the period gap will have to deal with this.     

VV Orionis – improved elements

TheUBV light curves obtained by Duerbeck (1975) andHa (wide) and Ha (narrow) light curves obtained by Chambliss & Davan (1987) of the detached eclipsing binary VV Orionis (VV Ori) were analysed using the Wilson-Devinney method fixing the two parametersT _h (25,000 K) and q(0.4172), resulting in the following absolute elements:A = 13.605 ± 0.03 LR_⊙,R _h = 5.03 ±0.03R _⊙, R_c = 2.43 ±0.02R _⊙,M _bol,h = -5.18 ± 0.11,M _bol,c = -1.54 ± 0.06,m _h =10.81 + 0.42m _⊙ andm _c = 4.51 ± 0.41m _⊙. The de-reddened colours obtained from applying the reddening corrections ofE(B-V) = 0^m.05 andE(U-B) = O ^m .04, and the derived temperatures of the components, gave spectral types ofB 1.5V for the primary and 54-5V with anUV excess of 0 ^m ·3 for the secondary component. A comparison of the logL and logT _e of the components with the observed ZAMS shows the primary component to be a little above and the secondary component to be a little below/or on the ZAMS. A comparison of the properties of the components of VV Ori and a few other detached systems with the normal stars in the logL, logR and logT _e versus logm planes, indicated a need for either a readjustment of the scales of the above parameters or modifications in the theoretical models. From the position of the components on the evolutionary tracks of Pop I composition computed by Schaller et al. (1992) it is noticed that while the primary component of W Ori had slightly evolved along the main-sequence, its secondary is still unevolved. The age of VV Ori is found to be 10 ± 1 million years and it is at a distance of 368 ± 10 pc.     

Spin-forbidden resonance multiplets in light elements

We present new laboratory data on the multiplets 2s ^{2 1} S -2s2p ³ P, 2s ²2p ² P - 2s2p ^{2 4} P, and 2s ²2p ^{2 3} P - 2s2p ^{3 5} S in nitrogen, oxygen and fluorine, and discuss theZ-dependence of their wave-numbers. These multiplets are very faint in laboratory light sources, but can become prominent in astrophysical sources of low density. Our results confirm the solar identifications of the nitrogen and oxygen multiplets made by Burtonet al. Predicted positions of the corresponding multiplets in neon are given.     

The differential rotation of ε Eri from MOST data

From high-precision MOST photometry spanning 35 days the existence of two spots rotating with slightly differing periods is confirmed. From the marginal probability distribution of the derived differential rotation parameter k its expectation value as well as confidence limits are computed directly from the data. The result depends on the assumed range in inclination i, not on the shape of the prior distributions. Two cases have been considered: (a) The priors for angles, inclination i of the star and spot latitudes β 1,2, are assumed to be constant over i, β 1, and β 2; (b) the priors are assumed to be constant over cos i, sin β 1, and sin β 2. In both cases the full range of inclination is considered: 0° ≤ i ≤ 90°. Scalefree parameters, i. e. periods and spot areas (in case of small spots) are taken logarithmically. Irrespective of the shape of the prior, k is restricted to 0.03 ≤ k ≤ 0.10 (1 σ limits). The inclination i of the star is photometrically ill-defined. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)