Oscillations of optical emission from flare stars and coronal loop diagnostics

Based on an analogy between stellar and solar flares, we investigate the ten-second oscillations detected in the U and B bands on the star EV Lac. The emission pulsations are associated with fast magnetoacoustic oscillations in coronal loops. We have estimated the magnetic field, B ≈ 320 G; the temperature, T ≈ 3.7 × 10⁷ K; and the plasma density, n ≈ 1.6 × 10¹¹ cm^?3, in the region of energy release. We provide evidence suggesting that the optical emission source is localized at the loop footpoints.     

Variability of hot supergiant IRAS 19336-0400 in the early phase of planetary nebula ionization

We present photoelectric and spectral observations of a hot candidate proto-planetary nebula—early B-type supergiant with emission lines in spectrum—IRAS 19336-0400. The light and color curves display fast irregular brightness variations with maximum amplitudes \(\Delta V = 0_ \cdot ^m 30\), \(\Delta B = 0_ \cdot ^m 35\), \(\Delta U = 0_ \cdot ^m 40\) and color-brightness correlations. By the variability characteristics IRAS 19336-0400 appears similar to other hot proto-planetary nebulae. Based on low-resolution spectra in the range λ4000–7500 Å we have derived absolute intensities of the emission lines Hα, Hβ, Hγ, [S II], [N II], physical conditions in gaseous nebula: n _e = 10⁴ cm^?3, T _e = 7000 ± 1000 K. The emission line Hα, Hβ equivalent widths are found to be considerably variable and related to light changes. By UBV-photometry and spectroscopy the color excess has been estimated: E _B-V = 0.50–0.54. Joint photometric and spectral data analysis allows us to assume that the star variability is caused by stellar wind variations.     

3D MHD Coronal Oscillations about a Magnetic Null Point: Application of WKB Theory

This paper is a demonstration of how the WKB approximation can be used to help solve the linearised 3D MHD equations. Using Charpit’s method and a Runge?–?Kutta numerical scheme, we have demonstrated this technique for a potential 3D magnetic null point, B=[x,ε y,?(ε+1)z]. Under our cold-plasma assumption, we have considered two types of wave propagation: fast magnetoacoustic and Alfvén waves. We find that the fast magnetoacoustic wave experiences refraction towards the magnetic null point and that the effect of this refraction depends upon the Alfvén speed profile. The wave and thus the wave energy accumulate at the null point. We have found that current buildup is exponential and the exponent is dependent upon ε. Thus, for the fast wave there is preferential heating at the null point. For the Alfvén wave, we find that the wave propagates along the field lines. For an Alfvén wave generated along the fan plane, the wave accumulates along the spine. For an Alfvén wave generated across the spine, the value of ε determines where the wave accumulation will occur: fan plane (ε=1), along the x-axis (0<ε<1) or along the y-axis (ε>1). We have shown analytically that currents build up exponentially, leading to preferential heating in these areas. The work described here highlights the importance of understanding the magnetic topology of the coronal magnetic field for the location of wave heating.     

The hydrogen emission spectrum of a shock wave in a stellar atmosphere

Based on a self-consistent solution of the equations of gas dynamics, kinetics of hydrogen atomic level populations, and radiative transfer, we analyze the structure of a shock wave that propagates in a partially ionized hydrogen gas. We consider the radiative transfer at the frequencies of spectral lines by taking into account the effects of a moving medium in the observer's frame of reference. The flux in Balmer lines is shown to be formed behind the shock discontinuity at the initial hydrogen recombination stage. The Doppler shift of the emission-line profile is approximately one and a half times smaller than the gas flow velocity in the Balmer emission region, because the radiation field of the shock wave is anisotropic. At Mach numbers M₁?10 and unperturbed gas densities σ₁=10^?10 g cm^?3, the Doppler shift is approximately one third of the shock velocity U₁. The FWHM of the emission-line profile δ _? is related to the shock velocity by δ _? ≈ k _? U₁, where k _? = 1, 0.6, and 0.65 for the Hα, Hβ, and Hγ lines, respectively.     

Ion acceleration and Alfvén wave generation at the Earth’s bow shock

The processes of ion acceleration and Alfvén wave generation by accelerated particles at the Earth’s bow shock are studied within a quasi-linear approach. Steady-state ion and wave spectra are shown to be established in a time of 0.3–4 h, depending on the background level of Alfvénic turbulence in the solar wind. The Alfvén waves produced by accelerated ions are confined within the frequency range 10^?2–1 Hz and their spectral peak with a wave amplitude βB ～ B comparable to the interplanetary magnetic field strength B corresponds to the frequency v = (2–3) × 10^?2 Hz. The high-frequency part of the wave spectrum (v > 0.2 Hz) undergoes damping by thermal ions. The calculated spectra of the accelerated ions and the Alfvén waves generated by them reproduce the main features observed in experiments.     

Propagation and Damping of a Localized Impulsive Longitudinal Perturbation in Coronal Loops

We use linear analysis to simulate the evolution of a coronal loop in response to a localized impulsive event. The disturbance is modeled by injecting a narrow Gaussian velocity pulse near one footpoint of a loop in equilibrium. Three different damping mechanisms, namely viscosity, thermal conduction, and optically thin radiation, are included in the loop calculations. We consider homogeneous and gravitationally stratified, isothermal loops of varying length (50≤L≤400 Mm) and temperature (2≤T≤10 MK). We find that a localized pulse can effectively excite slow magnetoacoustic waves that propagate up along the loop. The amplitudes of the oscillations increase with decreasing loop temperature and increasing loop length and size of the pulse width. At T≥4 MK, the waves are dissipated by the combined effects of viscosity and thermal conduction, whereas at temperatures of 2 MK, or lower, wave dissipation is governed by radiative cooling. We predict periods in the range of 4.6?–?41.6 minutes. The wave periods remain unaltered by variations of the pulse size, decrease with the loop temperature, and increase almost linearly with the loop length. In addition, gravitational stratification results in a small reduction of the periods and amplification of the waves as they propagate up along the loop.     

Pulsations of microwave emission and flare plasma diagnostics

We consider the modulation of nonthermal gyrosynchrotron emission from solar flares by the ballooning and radial oscillations of coronal loops. The damping mechanisms for fast magnetoacoustic modes are analyzed. We suggest a method for diagnosing the plasma of flare loops that allows their main parameters to be estimated from peculiarities of the microwave pulsations. Based on observational data obtained with the Nobeyama Radioheliograph (17 GHz) and using a technique developed for the event of May 8, 1998, we determined the particle density n≈3.7×10¹⁰ cm^?3, the temperature T≈4×10⁷ K, and the magnetic field strength B≈220 G in the region of flare energy release. A wavelet analysis for the solar flare of August 28, 1999, has revealed two main types of microwave oscillations with periods P₁≈7, 14 s and P₂≈2.4 s, which we attribute to the ballooning and radial oscillations of compact and extended flare loops, respectively. An analysis of the time profile for microwave emission shows evidence of coronal loop interaction. We determined flare plasma parameters for the compact (T≈5.3×10⁷ K, n≈4.8≈10¹⁰ cm^?3, B≈280 G) and extended (T≈2.1≈10⁷ K, n≈1.2≈10¹⁰ cm^?3, B≈160 G) loops. The results of the soft X-ray observations are consistent with the adopted model.     

Spectral and photometric studies of polar CRTS CSS 130604 J 215427+155714

We present the results of spectroscopic and photometric studies of a new polar CRTS CSS130604 J 215427+155714, conducted at the telescopes of the SAO RAS. Analysis of the photometric series of observations allowed to clarify the orbital period of the system, P _o = 0_. ^d 0672879 (±0.0000003). We build radial velocity curves and trace the intensity variations in the Hβ and H_γ hydrogen lines and He II λ 4686 ?A ionized heliumline. Based on the Hβ and He II lines we build Doppler maps. It is shown that the line formation region is localized near the Lagrange point. The following parameter estimates of the system are obtained:M ₁ = 0.83 ± 0.10M _⊙, M ₂ = 0.15 ± 0.01M _⊙, q = M ₂/M ₁ = 0.18 ± 0.03, i = 53? ± 5?. Based on the results of spectral, photometric and previously published polarimetric observations the possible geometric model of the system is discussed.     

Kinematics of the Galaxy from OB Stars with Data from the Gaia DR2 Catalogue

We have selected and analyzed a sample of OB stars with known line-of-sight velocities determined through ground-based observations and with trigonometric parallaxes and propermotions from the Gaia DR2 catalogue. Some of the stars in our sample have distance estimates made from calcium lines. A direct comparison with the trigonometric distance scale has shown that the calcium distance scale should be reduced by 13%. The following parameters of the Galactic rotation curve have been determined from 495 OB stars with relative parallax errors less than 30%: (U, V,W)_⊙ = (8.16, 11.19, 8.55)± (0.48, 0.56, 0.48) km s^?1, Ω₀ = 28.92 ± 0.39 km s^?1 kpc^?1, Ω'₀ = ?4.087 ± 0.083 km s^?1 kpc^?2, and Ω″ ₀ = 0.703 ± 0.067 km s?1 kpc^?3, where the circular velocity of the local standard of rest is V0 = 231 ± 5 km s^?1 (for the adopted R₀ = 8.0 ± 0.15 kpc). The parameters of the Galactic spiral density wave have been found from the series of radial, V_R, residual tangential, ΔV_circ, and vertical, W, velocities of OB stars by applying a periodogram analysis. The amplitudes of the radial, tangential, and vertical velocity perturbations are f_R = 7.1± 0.3 km s^?1, f_θ = 6.5 ± 0.4 km s^?1, and f_W = 4.8± 0.8 km s^?1, respectively; the perturbation wavelengths are λ_R = 3.3 ± 0.1 kpc, λ_θ = 2.3 ± 0.2 kpc, and λ_W = 2.6 ± 0.5 kpc; and the Sun’s radial phase in the spiral density wave is (χ_⊙)_R = ?135? ± 5?, (χ_⊙)θ = ?123? ± 8?, and (χ_⊙)_W = ?132? ± 21? for the adopted four-armed spiral pattern.     

Kinematics of B-F Stars as a Function of Their Dereddened Color from Gaia and PCRV Data

Parallaxes with an accuracy better than 10% and proper motions from the Gaia DR1 TGAS catalogue, radial velocities from the Pulkovo Compilation of Radial Velocities (PCRV), accurate Tycho-2 photometry, theoretical PARSEC, MIST, YaPSI, BaSTI isochrones, and the most accurate reddening and interstellar extinction estimates have been used to analyze the kinematics of 9543 thin-disk B-F stars as a function of their dereddened color. The stars under consideration are located on the Hertzsprung–Russell diagram relative to the isochrones with an accuracy of a few hundredths of a magnitude, i.e., at the level of uncertainty in the parallax, photometry, reddening, extinction, and the isochrones themselves. This has allowed us to choose the most plausible reddening and extinction estimates and to conclude that the reddening and extinction were significantly underestimated in some kinematic studies of other authors. Owing to the higher accuracy of TGAS parallaxes than that of Hipparcos ones, the median accuracy of the velocity components U, V, W in this study has improved to 1.7 km s^?1, although outside the range ?0.1 ^m < (B _T ? V _T )₀ < 0.5 ^m the kinematic characteristics are noticeably biased due to the incompleteness of the sample. We have confirmed the variations in the mean velocity of stars relative to the Sun and the stellar velocity dispersion as a function of their dereddened color known from the Hipparcos data. Given the age estimates for the stars under consideration from the TRILEGAL model and the Geneva–Copenhagen survey, these variations may be considered as variations as a function of the stellar age. A comparison of our results with the results of other studies of the stellar kinematics near the Sun has shown that selection and reddening underestimation explain almost completely the discrepancies between the results. The dispersions and mean velocities from the results of reliable studies fit into a ±2 km s^?1 corridor, while the ratios σ _V /σ _U and σ _W /σ _U fit into ±0.05. Based on all reliable studies in the range ?0.1 ^m < (B _T ? V _T )₀ < 0.5^m, i.e., for an age from 0.23 to 2.4 Gyr, we have found: W_⊙ = 7.15 km s^?1, \({\sigma _U} = 16.0{e^{1.29({B_T} - {V_T})o}}\), \({\sigma _V} = 10.9{e^{1.11({B_T} - {V_T})o}}\), \({\sigma _W} = 6.8{e^{1.46({B_T} - {V_T})o}}\), the stellar velocity dispersions in km s^?1 are proportional to the age in Gyr raised to the power β _U = 0.33, β _V = 0.285, and β _W = 0.37.     

The matter in the Big-Bang model is dust and not any arbitrary perfect fluid!

We consider a spherically symmetric general relativistic perfect fluid in its comoving frame. It is found that, by integrating the local energy momentum conservation equation, a general form of g ₀₀ can be obtained. During this study, we get a cue that an adiabatically evolving uniform density isolated sphere having ρ(r,t)=ρ ₀(t), should comprise “dust” having p ₀(t)=0; as recently suggested by Durgapal and Fuloria (J. Mod. Phys. 1:143, 2010) In fact, we offer here an independent proof to this effect. But much more importantly, we find that for the homogeneous and isotropic Friedmann-Robertson-Walker (FRW) metric having p(r,t)=p ₀(t) and ρ(r,t)=ρ ₀(t), \(g_{00} = e^{-2p_{0}/(p_{0} +\rho_{0})}\). But in general relativity (GR), one can choose an arbitrary t→t _?=f(t) without any loss of generality, and thus set g ₀₀(t _?)=1. And since pressure is a scalar, this implies that p ₀(t _?)=p ₀(t)=0 in the Big-Bang model based on the FRW metric. This result gets confirmed by the fact the homogeneous dust metric having p(r,t)=p ₀(t)=0 and ρ(r,t)=ρ ₀(t) and the FRW metric are exactly identical. In other words, both the cases correspond to the same Einstein tensor \(G^{a}_{b}\) because they intrinsically have the same energy momentum tensor \(T^{a}_{b}=\operatorname {diag}[\rho_{0}(t), 0,0, 0]\).     

Galactic kinematics from data on open star clusters from the MWSC catalogue

Open star clusters from the MWSC (Milky Way Star Clusters) catalogue have been used to determine the Galactic rotation parameters. The circular rotation velocity of the solar neighborhood around the Galactic center has been found from data on more than 2000 clusters of various ages to be V ₀ = 236 ± 6 km s^?1 for the adopted Galactocentric distance of the Sun R ₀ = 8.3 ± 0.2 kpc. The derived angular velocity parameters are Ω ₀ = 28.48 ± 0.36 km s^?1 kpc^?1, Ω’₀ = ?3.50 ± 0.08 km s^?1 kpc_?2, and Ω″₀ = 0.331 ± 0.037 km s^?1 kpc^?3. The influence of the spiral density wave has been detected only in the sample of clusters younger than 50 Myr. For these clusters the amplitudes of the tangential and radial velocity perturbations are f ^θ = 5.6 ± 1.6 km s^?1 and f _R = 7.7 ± 1.4 km s^?1, respectively; the perturbation wavelengths are λ _θ = 2.6 ± 0.5 kpc (i _θ = ^?11? ± 2?) and λ _R = 2.1 ± 0.5 kpc (i _R = ?9? ± 2?) for the adopted four-armed model (m = 4). The Sun’s phase in the spiral density wave is (χ_⊙)_θ = ?62? ± 9? and (χ_⊙)_R = ?85? ± 10? from the residual tangential and radial velocities, respectively.     

On a family of well behaved perfect fluid balls as astrophysical objects in general relativity

A family of well behaved perfect fluid balls has been derived starting with the metric potential g ₄₄=B(1+Cr ²) ⁿ for all positive integral values of n. For n≥4, the members of this family are seen to satisfy the various physical conditions e.g. c ² ρ≥p≥0,dp/dr<0,dρ/dr<0, along with the velocity of sound \((\sqrt{dp/c^{2}d\rho} )< 1\) and the adiabatic index ((p+c ² ρ)/p)(dp/(c ² dρ))>1. Also the pressure, energy density, velocity of sound and ratio of pressure and energy density are of monotonically decreasing towards the pressure free interface (r=a). The fluid balls join smoothly with the Schwarzschild exterior model at r=a. The well behaved perfect fluid balls so obtained are utilised to construct the superdense star models with their surface density 2×10¹⁴  gm/cm³. We have found that the maximum mass of the fluid balls corresponding to various values of n are decreasing with the increasing values of n. Over all maximum mass for the whole family turns out to be 4.1848M _Θ and the corresponding radius as 19.4144 km while the red shift at the centre and red shift at surface as Z ₀=1.6459 and Z _a =0.6538 respectively this all happens for n=4. It is interesting to note that for higher values of n viz n≥170, the physical data start merging with that of Kuchowicz superdense star models and hence the family of fluid models tends to the Kuchowicz fluid models as n→∞. Consequently the maximum mass of the family of solution can not be less than 1.6096 M _Θ which is the maximum mass occupied by the Kuchowicz superdense ball. Hence each member of the family for n≥4 provides the astrophysical objects like White dwarfs, Quark star, typical neutron star.     

Magnetic-field dependence of chemical anomalies in CP stars

The dependence of the degree of anomaly of parameter Z of Geneva photometry (Z0 = Z _CP ?Z _norm.) on the average surface magnetic field Bs is analyzed. The Z0 value is proportional to the degree of anomaly of chemical composition. It was found that Bs → 0 corresponds Z0 → ?0.010÷ ?0.015, i.e., part of CP stars are virtually devoid of magnetic field, but exhibit chemical anomalies. This effect may be due to selection whereby only objects with strong chemical anomalies are classified as CP stars, thereby producing a deficit of stars with relatively weak anomalies. Moreover, CP stars have other sources of stabilization of their atmospheres besides the magnetic field, e.g., slow rotation. Formulas relating Z0 to Bs are derived.     

Synchronous <Emphasis Type="Italic">UBVR</Emphasis> photometry and spectroscopy of DI Cephei

We present our synchronous spectroscopy and photometry of DI Cep, a classical T Tauri star. The equivalent widths and radial velocities of the individual components and Hα, Hβ, D1 and D2 Na I, and HeI λ5876 Å emission line profiles exhibit variability. We have found a clear positive correlation between the brightness and equivalent width for the Hα and Hβ emission lines. The photometric and spectroscopic data are satisfactorily described in phases of a 9-day period. The expected magnetic field of the star has been estimated using existing magnetospheric models to be 655–1000 G. The star is suspected to be a binary.     

Variability and rapid evolution of the protoplanetary object IRAS 18062+2410=V886 her

We present our long-term photometric and spectroscopic observations of a high-latitude B supergiant with an infrared excess—the protoplanetary nebula IRAS 18062+2410. OurU BV observations in 2000–2006 have confirmed the rapid irregular photometric variability of the star with a maximum amplitude as high as 0 _. ^m 4 in V that we found previously. The B—V and U—B color indices vary with amplitudes as high as 0 _. ^m 10 and 0 _. ^m 25, respectively, and show no clear correlation with the brightness. Our V-band CCD observations on 11 nights in 2006 have revealed brightness trends during the night. The variability of IRAS18062+2410 is similar in pattern to the light variations in other hot post-AGB objects and some of the nuclei of young planetary nebulae. We assume that pulsations and a variable stellar wind can be responsible for the variability of these stars. In addition to the rapid variability, our 12-year-long observations have revealed a systematic decline in the mean brightness of IRAS 18062+2410. This may be related to a rise in the temperature of the star at constant luminosity as a result of its evolution. Low-resolution spectroscopic observations have shown a systematic increase in the equivalent widths of the Hα, Hβ, [NII]λ6584 Å, OI λ8446 Å, and [OII] λ7320–7330 Å emission lines. The changes in the star’s emission line spectrum are probably caused by an increase in the degree of ionization of the gas shell due to a rise in the temperature of the ionizing star. Our photometric and spectroscopic observations of IRAS 18062+2410 confirm the previously made assumptions that the star evolves very rapidly to the region of planetary nebulae.     

3D stellar reddening map from 2MASS photometry: An improved version

An improved version of the 3D stellar reddening map in a space with a radius of 1200 pc around the Sun and within 600 pc of the Galactic midplane is presented. As in the previous 2010 and 2012 versions of the map, photometry with an accuracy better than 0.05 ^m in the J and Ks bands for more than 70 million stars from the 2MASS catalogue is used in the new version. However, the data reduction technique is considerably more complicated. As before, an analysis of the distribution of stars near the main-sequence turnoff on the (J ? Ks)?Ks diagram, where they form a distribution maximum, provides a basis for the method. The shift of this maximum, i.e., the mode (J ? Ks), along (J ? Ks) and Ks, given the spatial variations of the mean dereddened color (J ? Ks)₀ of these stars, is interpreted as a growth of the reddening with increasing distance. The main distinction of the new method is that instead of the fixed mean absolute magnitude, dereddened color, distance, and reddening for each cell, the individual values of these quantities are calculated for each star by iterations when solving the system of equations relating them. This has allowed one to increase the random accuracy of the map to 0.01 ^m and its spatial resolution to 20 pc in coordinates and distance and to 1° in longitude and latitude. Comparison with other reddening estimates for the same spatial cells and Gaia DR1 TGAS stars shows that the constructed map is one of the best maps for the space under consideration. Its systematic errors have been estimated to be σ(E(J ? Ks)) = 0.025 ^m , or σ(E(B ? V)) = 0.04 ^m . The main purpose of the map is to analyze the characteristics of Galactic structures, clouds, and cloud complexes. For this purpose, the reddening map within each spatial cell has also been computed by analyzing the reddening along each line of sight.     

Orbital stability of systems of closely-spaced planets,II: configurations with coorbital planets

We numerically investigate the stability of systems of 1 \({{\rm M}_{\oplus}}\) planets orbiting a solar-mass star. The systems studied have either 2 or 42 planets per occupied semimajor axis, for a total of 6, 10, 126, or 210 planets, and the planets were started on coplanar, circular orbits with the semimajor axes of the innermost planets at 1 AU. For systems with two planets per occupied orbit, the longitudinal initial locations of planets on a given orbit were separated by either 60° (Trojan planets) or 180°. With 42 planets per semimajor axis, initial longitudes were uniformly spaced. The ratio of the semimajor axes of consecutive coorbital groups in each system was approximately uniform. The instability time for a system was taken to be the first time at which the orbits of two planets with different initial orbital distances crossed. Simulations spanned virtual times of up to 1 × 10⁸, 5 × 10⁵, and 2 × 10⁵ years for the 6- and 10-planet, 126-planet, and 210-planet systems, respectively. Our results show that, for a given class of system (e.g., five pairs of Trojan planets orbiting in the same direction), the relationship between orbit crossing times and planetary spacing is well fit by the functional form log(t _c /t ₀) = b β + c, where t _c is the crossing time, t ₀ = 1 year, β is the separation in initial orbital semimajor axis (in terms of the mutual Hill radii of the planets), and b and c are fitting constants. The same functional form was observed in the previous studies of single planets on nested orbits (Smith and Lissauer 2009). Pairs of Trojan planets are more stable than pairs initially separated by 180°. Systems with retrograde planets (i.e., some planets orbiting in the opposite sense from others) can be packed substantially more closely than can systems with all planets orbiting in the same sense. To have the same characteristic lifetime, systems with 2 or 42 planets per orbit typically need to have about 1.5 or 2 times the orbital separation as orbits occupied by single planets, respectively.     

HI content in galactic disks: The role of gravitational instability

We examine the dependence of the total hydrogen mass M _HI in late-type star-forming galaxies on rotation velocity V _rot and optical size D ₂₅ or radial scale length R ₀ of the disk for two samples of galaxies: (i) isolated galaxies (AMIGA) and (ii) galaxies with edge-on disks (flat galaxies according to Karachentsev et al.). M _HI given in the HYPERLEDA database for flat galaxies have turned out to be, on average, overestimated by ~0.2 dex compared to isolated galaxies with similar V _rot or D ₂₅, which is apparently due to an overestimation of the self-absorption in the HI line. The hydrogen mass in the galaxies of both samples closely correlates with the total specific angular momentum of the galactic disk J, which is proportional to V _rot D ₂₅ or V _rot R ₀, with the low-surface-brightness galaxies lying along the common V _rot R ₀ sequence. We discuss the possibility of explaining the relationship between M _HI and V _rot D ₂₅ by assuming that the gas mass in the disk is regulated by the marginal gravitational stability condition for the gas layer. Comparison of the observed and theoretically expected dependences leads us to conclude that either the gravitational stability corresponds to higher values of the Toomre parameter than is usually assumed, or the threshold stability condition formost galaxies was fulfilled only in the past, when the gasmass in the disks was a factor of 2–4 higher than that at present (except for the galaxies with an anomalously high observed HI content). The latter condition requires that for most galaxies the conversion of gas into stars be not compensated by the external accretion of gas onto the disk.