Lithium abundances in atmospheres of carbon stars estimated with the help of modern lists of spectral lines

We determine lithium abundances in atmospheres of three carbon stars from synthetic spectrum fitting in the λλ 668–674 nm range using the Li I λ 670.8 nm resonance doublet. To produce synthetic spectra, we use a modified list of atomic lines from the VALD database and three alternative line lists of CN and C₂ molecules which are modifications of line lists from the Jorgensen’s website () and from the Kurucz database (1993, CD-ROM nos. 1–23). The spectral lines from these lists were tested by matching synthetic spectra to observed spectra of the sun, Arcturus, and early R star HD 100764. We perform analysis of the blends involving the Li I λ 670.8 nm doublet in spectra of N stars AW Cyg and UX Dra. The lithium abundances in HD 100764, AW Cyg, and UX Dra are estimated to be lgN(Li) ≈ 2, −1.4, and −0.9, respectively. Discrepancies of lithium abundances lgN(Li) obtained with the help of molecular line lists do not seem to exceed 0.2 dex.     

Pulsations of the sun and a beat period of 399 days

Measurements of the Doppler effect of the solar photosphere have been carried out at the Crimean Astrophysical Observatory for 37 years, beginning in 1974 (in total, 2188 days or 13 247 h). The measurements use the differential center-to-limb method of registration of line-of-sight velocity with a solar magnetograph (in the iron absorption line λ512.4 nm). As a result of the experiment, two global pulsations of the sun with periods P ₀ = 9600.606(12) and P ₁ = 9597.936(16) s have been discovered. The nature of the periods is unknown. The first pulsation was detected in 1974–1982; the second, during nearly all the 37 years. The 2008–2010 data confirm the stability of the initial phase of the P ₁ pulsation with a mean (differential) amplitude of 0.25 m/s. The fact that the beat period of the two pulsations coincides with the synodic period of Jupiter’s orbital revolution, i.e., 399 days, raises a new, complex problem for solar physics and cosmogony.     

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.     

Gamma-ray emission expected from Kepler’s SNR

Nonlinear kinetic theory of cosmic ray (CR) acceleration in supernova remnants (SNRs) is used to investigate the properties of Kepler’s SNR and, in particular, to predict the γ-eay spectrum expected from this SNR. Observations of the nonthermal radio and X-ray emission spectra as well as theoretical constraints for the total supernova (SN) explosion energy E _sn are used to constrain the astronomical and particle acceleration parameters of the system. Under the assumption that Kepler’s SN is a type Ia SN we determine for any given explosion energy E _sn and source distance d the mass density of the ambient interstellar medium (ISM) from a fit to the observed SNR size and expansion speed. This makes it possible to make predictions for the expected γ-eay flux. Exploring the expected distance range we find that for a typical explosion energy E _sn=10⁵¹ erg the expected energy flux of TeV γ-rays varies from 2×10⁻¹¹ to 10⁻¹³ erg/(cm² s) when the distance changes from d=3.4 kpc to 7 kpc. In all cases the γ-eay emission is dominated by π ⁰-decay γ-rays due to nuclear CRs. Therefore Kepler’s SNR represents a very promising target for instruments like H.E.S.S., CANGAROO and GLAST. A non-detection of γ-rays would mean that the actual source distance is larger than 7 kpc.     

Shear-free homogeneous cosmological model with heat flux

We present the condition of vanishing shear in a spatially homogeneous spacetime in terms of the Ricci rotation co-efficients corresponding to an orthonormal tetrad (ν ^α. _A η ^α) (whereν ^α is the unit vector along the time axis and _A η ^α are the three independent reciprocal group vectors). Assuming that the velocity vector can be expanded in the direction ofν ^α and any one of the _A η ^α’s it is shown that shear-free motion is possible only in case of some special Bianchi types, and these cases are studied assuming the velocity vector to be geodetic and that there may be a nonvanishing heat flux term.     

Effect of fluctuations in Doppler width on the center of strong absorption lines

Stochastic temperatures and turbulence are characterized by average velocities u _th and < u _turb  > ≡ u ₀ and fluctuations u￠_th {u'_{th}} and u′ (<u′ > = 0). Thus, the Doppler width of a line also has a fluctuating component Dl￠_D \Delta {\lambda '_D} . Observed spectra correspond to the radiative flux averaged over time and over a star’s surface, <H_λ>. Usually, only the average velocities u _th and u ₀ are taken into account in photospheric models and these yield the Doppler width Dl_D⁽⁰⁾ \Delta \lambda_D^{(0)} of a line in the customary way. The fluctuations Dl￠_D \Delta {\lambda '_D} mean that near a line center the average absorption coefficient < α_λ > is larger than the usual α_λ, which depends only on the average velocities u _th and u ₀. This enhances the absorption line near the center and is not explained by the photospheric models. This new statistical effect depends on the wavelength of the line. A comparison of observed lines with model profiles yields an estimate for the average level of fluctuations in the Doppler width, h = á | Dl￠_D | ñ

Identification of absorption lines of rare-earth elements in the region of the spectrum λλ 610.25–610.57 nm roAp stars,HD 101065, HD 134214, HD 137949, and HD 24712

In the spectral region λλ = 610.25−610.57 nm of the roAp star HD 101065, we performed the identification of the rare-earth element absorption lines absent from the atomic spectral line databases VALD and DREAM. The identified lines were used for calculating the synthetic spectra of the roAp stars HD 137949, HD 134214, and HD 24712. The upper limit of the identified line estimates was determined. The rotation/pulsation parameter νsini and magnetic field modulus were determined using the lines Nd III 669.083 nm and Ca I 616.217 nm. The element abundance in the roAp stars HD 134214 and HD 24712 was determined for the first time using Nd III lines. Calculations were carried out with the help of the code SynthM developed by S.A. Khan.     

On accounting for temperature and pressure in determining the Rayleigh scattering in the earth’s atmosphere

We calculated the variations of Rayleigh optical depth with changes of pressure and temperature for three observation sites: Simferopol (φ = 44°57′N, λ = 34°8′E, h = 265 m above sea level), Nauchny (φ = 44°43′N, λ = 34°3′E, h = 583 m), and Ai-Petry meteorological station (φ = 44°24′N, λ = 34°6′E, h = 1180 m).     

Neutron stars in generalized f(R) gravity

A gravity theory is considered with the Einstein-Hilbert Lagrangean R+aR ²+bR _μν R ^μν , R _μν being Ricci’s tensor and R the curvature scalar. The parameters a and b are taken of order 1 km². Arguments are given which suggest that the effective theory so obtained might be a fair approximation of a viable theory. A numerical integration is performed of the field equations for a free neutron gas. The result is that the star mass increases with increasing central density until about 1 solar mass and then decreases. The baryon number increases monotonically, which suggests that the theory allows stars in equilibrium with arbitrary baryon number, no matter how large.     

Measurement of the orbital period of the X-ray burster GS 1826-238 based on observations of its optical brightness variations

The variability of the optical and X-ray fluxes from the binary GS 1826-238 is investigated. An epoch-folding analysis of the optical data obtained with the RTT-150 telescope in 2003–2004 has revealed periodic brightness variations in the source with a period P _orb = 2.24940 ± 0.00015 h with a high statistical significance. When estimating the detection significance of the periodic signal, we have specially taken into account the presence of a powerful aperiodic component (“red noise”) in the source’s brightness variability. The source’s power density spectra in the frequency range ∼10⁻⁵–0.01 Hz have been obtained. We have detected a statistically significant break in the power density spectrum of GS 1826-238 at a frequency ν _br ≈ (8.48 ± 0.14) × 10⁻⁵ Hz in both optical and X-ray energy bands. We have estimated the orbital period of the binary GS 1826-238 using the correlation between the break frequency in the power density spectrum and the orbital period of binaries, P _orb ∝ 1/ν _br, found by Gilfanov and Arefiev (2005): P _orb = 3.7 ± 0.8 h and P _orb = 11.3 ± 5.9 h when using Sco X-1 and 1H 16267-273, respectively, as reference sources. It seems to us that the method for estimating the orbital periods of low-mass X-ray binaries using the correlation P _orb ∝ 1/ν _br may turn out to be very promising, especially for persistent low-luminosity X-ray binaries.     

HI fluctuations at large redshifts: I-visibility correlation

We investigate the possibility of probing the large scale structure in the universe at large redshifts by studying fluctuations in the redshifted 1420 MHz emission from the neutral hydrogen (HI) at early epochs. The neutral hydrogen content of the universe is known from absorption studies forz ≲ 4.5. TheHI distribution is expected to be inhomogeneous in the gravitational instability picture and this inhomogeneity leads to anisotropy in the redshifted HI emission. The best hope of detecting this anisotropy is by using a large low-frequency interferometric instrument like the Giant Meter-Wave Radio Telescope (GMRT). We calculate the visibility correlation function 〈V_v(U) V_v′(U)〉 at two frequenciesi andv′ of the redshiftedHI emission for an interferometric observation. In particular we give numerical results for the two GMRT channels centered aroundν = 325 MHz andν = 610 MHz from density inhomogeneity and peculiar velocity of the HI distribution. The visibility correlation is- 10^-10-10^-9 Jy². We calculate the signal-to-noise for detecting the correlation signal in the presence of system noise and show that the GMRT might detect the signal for integration times - 100 hrs. We argue that the measurement of visibility correlation allows optimal use of the uncorrelated nature of the system noise across baselines and frequency channels. On leave from Harish-Chandra Research Institute, Chhatnag Road, Jhusi, Allahabad 211 019, India.     

Bragg diffraction and the iron crust of cold neutron stars

If cooled-down neutron stars have a thin atomic crystalline–iron crust, they must diffract X-rays of appropriate wavelength. If the diffracted beam is to be visible from Earth (an extremely rare but possible situation), the illuminating source must be very intense and near the reflecting star. An example is a binary system composed of two neutron stars in close orbit, one of them inert, the other an X-ray pulsar. (Perhaps an “anomalous” X-ray pulsar or magnetar, not powered by gas absorption from the companion or surrounding space, would be the cleanest example.) The observable to be searched for is a secondary peak added (quasi-) periodically to the main X-ray pulse. The distinguishing feature of this secondary peak is that it appears at wavelengths related by simple integer numbers, λ,λ/2,λ/3,…,λ/n because of Bragg’s diffraction law.     

Tidal torque constants as a critical test for an evolution of close eclipsing binaries

Individual tidal torque λ ₂,E ₂ and apsidal-motion k ₂ constants were calculated for 112 close eclipsing binaries (CEBs) with Detached components belonging to the Main Sequence (DMS-type) from the catalogue by Svechnikov and Perevozkina (Catalogue of orbital elements, masses and luminosities of variable stars of DMS-type and some results of its statistical treatment, Ural State University Press, Yekaterinburg, pp. 1–5, 1999) and for 95 detached binaries taken from the catalogue by Torres et al. (Astron. Astrophys. Rev. 18:67, 2010) on the base of theoretical evolutionary stellar models including tidal torque constants by Claret (Astron. Astrophys. 424:919, 2004). A method of the inversion of model track grid into isochrones was formulated as a complex interpolation procedure for DMS-binaries data. Sets of isochrones were computed in k ₂–M, k ₂–R, λ ₂–M, λ ₂–R, E ₂–M, and E ₂–R planes. Calculated tidal torque constants allow to test stellar structure theory by comparing observed and estimated values of apsidal motion period and analyzing the correlation between timescales of synchronization, circularization, magnetic braking, as well as nuclear burning of DMS-components.     

New evidence for strong nonthermal effects in Tycho’s supernova remnant

For the case of Tycho’s supernova remnant (SNR) we present the relation between the blast wave and contact discontinuity radii calculated within the nonlinear kinetic theory of cosmic ray (CR) acceleration in SNRs. It is demonstrated that these radii are confirmed by recently published Chandra measurements which show that the observed contact discontinuity radius is so close to the shock radius that it can only be explained by efficient CR acceleration which in turn makes the medium more compressible. Together with the recently determined new value E _sn=1.2×10⁵¹ erg of the SN explosion energy this also confirms our previous conclusion that a TeV γ-ray flux of (2–5)×10⁻¹³ erg/(cm² s) is to be expected from Tycho’s SNR. Chandra measurements and the HEGRA upper limit of the TeV γ-ray flux together limit the source distance d to 3.3≤d≤4 kpc.     

Spectroscopic study of the variable planetary nebula IC 4997

Based on our spectroscopic observations of the variable planetary nebula IC 4997 in 2003–2009, we have obtained the relative fluxes in optical emission lines. The interstellar extinction c = 0.35 has been found from the Balmer decrement by taking into account the effect of self-absorption in hydrogen lines in dense nebular regions. The variations in the Balmer decrement point to variability of the self-absorption and circumstellar extinction. We have investigated the variations in the relative intensities of some spectral lines and their ratios with time. The drop in the ratios F(λ4363)/F(Hγ) and F(λ363)/F(λ4959) that began back in 1990–1995 has continued, suggesting a decrease in the electron density and temperature in the central nebular region. The ratio F(λ6731)/F(λ6717) has remained constant. It gives an estimate for the electron density in the outer regions of IC 4997, N _e ∼ 10⁴ cm⁻³.     

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 λ.     

Exploring morphological correlations among H2CO, 12CO, MSX and continuum mappings

There are relatively few H₂CO mappings of large-area giant molecular cloud (GMCs). H₂CO absorption lines are good tracers for low-temperature molecular clouds towards star formation regions. Thus, the aim of the study was to identify H₂CO distributions in ambient molecular clouds. We investigated morphologic relations among 6-cm continuum brightness temperature (CBT) data and H₂CO (1₁₁−1₁₀; Nanshan 25-m radio telescope), ¹²CO (1–0; 1.2-m CfA telescope) and midcourse space experiment (MSX) data, and considered the impact of background components on foreground clouds. We report simultaneous 6-cm H₂CO absorption lines and H110α radio recombination line observations and give several large-area mappings at 4.8 GHz toward W49 (50′×50′), W3 (70′×90′), DR21/W75 (60′×90′) and NGC2024/NGC2023 (50′×100′) GMCs. By superimposing H₂CO and ¹²CO contours onto the MSX color map, we can compare correlations. The resolution for H₂CO, ¹²CO and MSX data was ∼10′, ∼8′ and ∼18.3″, respectively. Comparison of H₂CO and ¹²CO contours, 8.28-μm MSX colorscale and CBT data revealed great morphological correlation in the large area, although there are some discrepancies between ¹²CO and H₂CO peaks in small areas. The NGC2024/NGC2023 GMC is a large area of HII regions with a high CBT, but a H₂CO cloud to the north is possible against the cosmic microwave background. A statistical diagram shows that 85.21% of H₂CO absorption lines are distributed in the intensity range from −1.0 to 0 Jy and the ΔV range from 1.206 to 5 km s⁻¹.     

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.     

Galactic kinematics from OB3 stars with distances determined from interstellar Ca II lines

Based on data for 102 OB3 stars with known proper motions and radial velocities, we have tested the distances derived by Megier et al. from interstellar Ca II spectral lines. The internal reconciliation of the distance scales using the first derivative of the angular velocity of Galactic rotation Ω′₀ and the external reconciliation with Humphreys’s distance scale for OB associations refined by Mel’nik and Dambis show that the initial distances should be reduced by ≈20%. Given this correction, the heliocentric distances of these stars lie within the range 0.6–2.6 kpc. A kinematic analysis of these stars at a fixed Galactocentric distance of the Sun, R ₀ = 8 kpc, has allowed the following parameters to be determined: (1) the solar peculiar velocity components (u _⊙, v _⊙, ω _⊙) = (8.9, 10.3, 6.8) ± (0.6, 1.0, 0.4) km s⁻¹; (2) the Galactic rotation parameters Ω₀ = −31.5 ± 0.9 km s⁻¹ kpc⁻¹, Ω′₀ = +4.49 ± 0.12 km s⁻¹ kpc⁻², Ω″₀ = −1.05 ± 0.38 km s⁻¹ kpc⁻³ (the corresponding Oort constants are A = 17.9 ± 0.5 km s⁻¹ kpc⁻¹, B = −13.6 ± 1.0 km s⁻¹ kpc⁻¹ and the circular rotation velocity of the solar neighborhood is |V ₀| = 252 ± 14 km s⁻¹); (3) the spiral density wave parameters, namely: the perturbation amplitudes for the radial and azimuthal velocity components, respectively, f _R = −12.5±1.1 km s⁻¹ and f _ϑ = 2.0 ± 1.6 km s⁻¹; the pitch angle for the two-armed spiral pattern i = −5.3° ± 0.3°, with the wavelength of the spiral density wave at the solar distance being λ = 2.3 ± 0.2 kpc; the Sun’s phase in the spiral wave x _⊙ = −91° ± 4°.     

Studies of magnetic fields in the solar photosphere using the line-ratio method

The longitudinal magnetic field measured using the Fe I λ 525 and Fe I λ 524.7 nm lines and global magnetic field of the sun differ depending on the observatory. To study the cause of these discrepancies, we calculate the H _‖(525)/H _‖(524.7) ratios for various combinations of magnetic elements and compare them with the corresponding observed values. We use the standard quiet model of the solar photosphere suggesting that there are magnetic fields of different polarities in the range between zero and several kilogauss. The magnetic element distribution is found as a function of magnetic field strength and the parameters of this distribution are determined for which the calculated H _‖(525)/H _‖(524.7) ratio agrees with the observed one. The sigma-components are found to be shifted differently for various points of the Fe I λ 525 nm profile calculated for the inhomogeneous magnetic field. The farther the point is from the line center, the larger the sigma-components shift. Such a peculiarity of the profiles may be responsible for the discrepancies in the measured values of the global magnetic field obtained at different observatories. The increase in modulus of the global magnetic field during the maxima of solar activity can be due to a larger fraction of magnetic elements with kilogauss magnetic fields.