Formation of a filamentary structure in planetary nebulae

Intricate filamentary structure and multiple shell-like appearance are very common phenomena in Planetary Nebulae.In addition, recent observations also indicate that the individual filaments present in these objects can have larger velocities than the adjacent smooth background (Pascoli, 1992 PASP 104, 350 and paper quoted therein).We have hypothesized that non linear hydrodynamical processes existing within the nebular gas are, possibly, responsible for these structures. As a matter of fact, it is argued that such a characteristic morphology, reinterpreted as a intermingled network of solitary waves or solitons, can be spontaneously generated in Planetary Nebulae as soon as one assumes that the nebular gas is permeated by a weak magnetic field whose strength is about 10^–5 to 10^–4 gauss.Main results of this work and further comments will be subsequently published in Ap&SS.     

Perturbation structures in magnetized relativistic winds

A perturbation formulation is presented for steady three-dimensional compressible structures embedded in a background relativistic magnetohydrodynamic (RMHD) radial outflow with spherical symmetry. The property of the two concurred RMHD fast and slow critical points is examined. Explicit perturbation solutions at large radii are derived analytically in a polytropic background wind. This perturbation approach for compressible structures together with our recent analyses on propagation of Alfvénic perturbations in a RMHD wind forms a useful starting basis for modeling a broad class of complex structures in magnetized relativistic astrophysical outflows.     

Shock modelling of Planetary Nebulae

The kinematics of Planetary Nebulae are analyzed in terms of the solutions to the equations of hydrodynamic equilibrium developed by J. Cantó. We apply our analysis to the Planetary Nebulae NGC 6905 and NGC 6537. A detailed spectroscopic study of these objects reveals the existence of high nuclear velocities, together with complex kinematic structures and unusual emission line intensities. Shock ionization clearly plays a key role in these nebulae. Remarkably good agreement is obtained when comparing the synthetic maps and spectra resulting from the shock solutions with the observational data.     

Bipolar outflows and the evolution of stars

Hypersonic bipolar outflows are a ubiquitous phenomena associated with both young and highly evolved stars. Observations of Planetary Nebulae, the nebulae surrounding Luminous Blue Variables such as η Carinae, Wolf Rayet bubbles, the circumstellar environment of SN 1987A and Young Stellar Objects all reveal high velocity outflows with a wide range of shapes. In this paper the current state of our theoretical understanding of these outflows is reviewed. Beginning with Planetary Nebulae considerable progress has been made in understanding bipolar outflows as the result of stellar winds interacting with the circumstellar environment. In what has been called the “Generalized Wind Blown Bubble” (GWBB) scenario, a fast tenuous wind from the central star expands into a ambient medium with an aspherical (toroidal) density distribution. Inertial gradients due to the gaseous torus quickly lead to an expanding prolate or bipolar shell of swept-up gas bounded by strong shock waves. Numerical simulations of the GWBB scenario show a surprisingly rich variety of gasdynamical behavior, allowing models to recover many of the observed properties of stellar bipolar outflows including the development of collimated supersonic jets. In this paper the physics behind the GWBB scenario is reviewed in detail and its strengths and weakness are considered. Alternative models involving MHD processes are also examined. Applications of these models to each of the principle classes of stellar bipolar outflow (YSO, PNe, LBV, SN87A) are then reviewed. Outstanding issues in the study of bipolar outflows are considered as are those questions which arise when the outflows are viewed as a single class of phenomena occurring across the HR diagram.     

GRMHD/RMHD simulations & stability of magnetized spine-sheath relativistic jets

A new general relativistic magnetohydrodynamics (GRMHD) code “RAISHIN” used to simulate jet generation by rotating and non-rotating black holes with a geometrically thin Keplarian accretion disk finds that the jet develops a spine-sheath structure in the rotating black hole case. Spine-sheath structure and strong magnetic fields significantly modify the Kelvin-Helmholtz (KH) velocity shear driven instability. The RAISHIN code has been used in its relativistic magnetohydrodynamic (RMHD) configuration to study the effects of strong magnetic fields and weakly relativistic sheath motion, c/2, on the KH instability associated with a relativistic, γ=2.5, jet spine-sheath interaction. In the simulations sound speeds up to and Alfvén wave speeds up to ∼0.56c are considered. Numerical simulation results are compared to theoretical predictions from a new normal mode analysis of the RMHD equations. Increased stability of a weakly magnetized system resulting from c/2 sheath speeds and stabilization of a strongly magnetized system resulting from c/2 sheath speeds is found.     

Horizontal velocities in the solar photosphere

Horizontal macroscopic velocities V _hor in the photosphere are studied. High-resolution spectrograms of quiet regions are analyzed for center-limb variation of rms Doppler shifts. The data are treated to assure that the observed velocities refer to constant size volumes on the Sun (800 × × 3000 × 250 km), independent of μ. Using known height variation of vertical velocities and calculated line formation heights, the height dependence of 〈V _hor〉 is obtained. From a value around 450 m s^?1 it decreases rapidly with increasing height. To study also small-scale velocities, the time evolution of subarcsecond size elements in the photospheric network (solar filigree) is studied on filtergrams. It is concluded that they show proper motions implying 〈V _hor〉 about 1 km s^?1.     

Internal structure and properties of Mars

Theoretical physical models of the Martia interior are presented in the light of new and revised data and constraints. These models include thermal evolution, densities, and seismic wave velocities. The interior of Mars appears to be Earth-like in many respects. Although thermal models indicate that Mars has passed its peak of evolution it may still have an asthenosphere and may be moderately active tectonically. Mars has an Fe-FeS core with a radius of and may be moderately active tectonically. Mars has an Fe-FeS core with a radius of 1500–2000 km. The mantle is enriched in FeO with an olivine composition of about Fo₇₅. Theoretically determined seismic wave velocities are relatively well constrained in the mantle with upper-mantle P_n velocities ranging from 7.64 to 7.80 km/sec. However, there are wide variations in V_P in the core dependent on composition. The shadow zone due to the core is larger than the Earth's.     

Korteweg-de Vries equation for ion acoustic soliton with negative ions in the presence of nonextensive electrons

Korteweg-de Vries (KdV) equation for electrostatic ion acoustic wave in a three component plasma containing positive and negative ions along with the nonextensive electrons is derived. Fast and slow ion acoustic modes which propagate with different velocities are excited. The effects of variation of quantities like q (nonextensive parameter), Q (mass ratio of positive to negative ion), μ (electron to positive ion number density ratio), θ _i (positive ion to electron temperature ratio) and θ _n (negative ion to electron temperature ratio) have been presented for fast and slow ion acoustic modes. Both compressive and rarefactive solitons are observed. It is found that the solitary excitations strongly depend on the mass and density ratios of the positive and negative ions as well as on nonextensive electron parameter.     

The 4-day rotation of the upper atmosphere of venus

Near u.v. photographs of Venus have been analysed. Evidence is presented that so-called Y- and ψ-shaped markings accelerate from local sunrise to some time during local afternoon; mean horizontal velocities are 83 msec^?1 (morning) and 122 msec^?1 (afternoon). Possibly these configurations are more stable than a normal cloud system. ‘Irregular’ features do not seem to share this acceleration; they circulate at higher altitudes, sometimes covering or distorting the stable configurations.     

The rate of flux pile-up magnetic reconnection in reduced MHD

The rate of two-dimensional flux pile-up magnetic reconnection is known to be severely limited by gas pressure in a low-beta plasma of the solar corona. As earlier perturbational calculations indicated, however, the pressure limitation should be less restrictive for three-dimensional flux pile-up. In this paper the maximum rate of reconnection is calculated in the approximation of reduced magnetohydrodynamics (RMHD), which is valid in the solar coronal loops. The rate is calculated for finite-magnitude reconnecting fields in the case of a strong axial field in the loop. Gas pressure effects are ignored in RMHD but a similar limitation on the rate of magnetic merging exists. Nevertheless, the magnetic energy dissipation rate and the reconnection electric field can increase by several orders of magnitude as compared with strictly two-dimensional pile-up. Though this is still not enough to explain the most powerful solar flares, slow coronal transients with energy release rates of order 10²⁵– 10²⁶ erg s^–1and heating of quiet coronal loops are within the compass of the model.     

Relative flow of H+ and O+ ions in the topside ionosphere at mid-latitudes

Theoretical results on the daily variation of O⁺ and H⁺ field-aligned velocities in the topside ionosphere are presented. The results are for an L = 3 magnetic field tube under sunspot minimum conditions at equinox. They come from calculations of time-dependent O⁺ and H⁺ continuity and momentum balance in a magnetic field tube which extends from the lower F2 region to the equatorial plane (Murphy et al., 1976).There are occasions when ion counterstreaming occurs, with the O⁺ velocity upward and H⁺ velocity downward. The conditions causing this counterstreaming are described: the H⁺ layer is descending whilst O⁺ is supplied from below either to increase the O⁺ concentration at fixed heights or to replace O⁺ ions lost by charge exchange with neutral H. It is suggested that the results of observations at Arecibo by Vickrey et al. (1976) of O⁺ and H⁺ concentrations and counterstreaming velocities are significantly affected by $\text{E}\text{×}\text{B}$ drift.     

Observations of photometric profiles of two planetary nebulae

It is noted that the Abel integral of a Gaussian function is a Gaussian function. This special property is used in order to develop an analytic method for studying photometric profile in Planetary Nebulae. Monochromatic observations have been made at the Haute-Provence Observatory using the RCA-CCD (320×512 pixels) mounted at theF/6 Newtonian focus of the 1.2 m telescope. Among the observed nebulae, two of them (NGC 6778 and NGC 7354) are regular enough to test the photometric profile method. The nebular emissivity per volume unit is then obtained for H⁺ and O⁺⁺ ions.     

Disappearence of equatorial Es, associated with magnetic field depressions

The phenomenon of the sudden disappearance of equatorial sporadic E at two stations, namely, Trivandrum and Kodaikanal is studied. It is established that whenever there is a sudden disappearance of E_s, there is a depression in the horizontal magnetic field (H) range. Electron velocities during the presence and absence of sporadic E have been estimated. These results show that the irregularities responsible for sporadic E are present even when the electron velocities are less than the ion thermal velocity.     

An HLLC Riemann solver for relativistic flows – II. Magnetohydrodynamics

An approximate Riemann solver for the equations of relativistic magnetohydrodynamics (RMHD) is derived. The Harten–Lax–van Leer contact wave (HLLC) solver, originally developed by Toro, Spruce and Spears, generalizes the algorithm described in a previous paper to the case where magnetic fields are present. The solution to the Riemann problem is approximated by two constant states bounded by two fast shocks and separated by a tangential wave. The scheme is Jacobian-free, in the sense that it avoids the expensive characteristic decomposition of the RMHD equations and it improves over the HLL scheme by restoring the missing contact wave.
Multidimensional integration proceeds via the single step, corner transport upwind (CTU) method of Colella, combined with the constrained transport (CT) algorithm to preserve divergence-free magnetic fields. The resulting numerical scheme is simple to implement, efficient and suitable for a general equation of state. The robustness of the new algorithm is validated against one- and two-dimensional numerical test problems.     

Dynamics of Hα spicules according to spectral observations at various heights of the solar chromosphere

Almost simultaneous height sequences of 69 spicules in the Hα line have been studied. The spectra are obtained at six heights during 6 s on the east side of the solar disk with the 53-cm Lyot coronagraph of Abastumani Astrophysical Observatory. Radial velocities V _r, total intensities or equivalent widths W, full widths at half maximum of intensity (FWHM) at all heights are determined (about 300 profiles of the Hα line). It is found that:

1. Absolute values of radial velocities increase linearly with the height (see Equation (1));
2. variation of the sign of the radial velocity along single spicules was never observed.

These results combined with the findings on the spicules radial velocities and shifts obtained earlier (Kulidzanishvili and Nikolsky, 1978; Nikolsky and Platova, 1970) led us to the conclusion that the 5-min tangential oscillations of spicules involve the entire spicule at once. The intensity height scales for single spicules and for the chromosphere ‘in toto’ are determined; they turned out to be 2.5 × 10³ km and 1.9 × 10³ km respectively (see Equations (2) and (3)). The dependence curve of the Hα line half-widths Δλ on the height h is drawn. The Hα line half-width for those spicule groups which are traced at all heights (10 spicules) decreases with the height (Figure 4); for the majority (～60 spicules) it remains essentially constant. Non-thermal ‘turbulent’ velocities V _t, in Hα spicules are defined. A mean value of the ‘turbulent’ velocity V _t at T = 6000° appeared to be 20–30km s^?1. The hydrogen concentration in the spicules at 5000 km is 6 × 10¹¹ cm^?3.     

On the Formation Height of the SDO/HMI Fe 6173 ? Doppler Signal

The Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO) is designed to study oscillations and the magnetic field in the solar photosphere. It observes the full solar disk in the Fe?i absorption line at 6173 Å. We use the output of a high-resolution, 3D, time-dependent, radiation-hydrodynamic simulation based on the CO ⁵ BOLD code to calculate profiles F(??,x,y,t) for the Fe?i 6173 Å line. The emerging profiles F(??,x,y,t) are multiplied by a representative set of HMI filter-transmission profiles R _i (??, 1??i??6) and filtergrams I _i (x,y,t; 1??i??6) are constructed for six wavelengths. Doppler velocities V _HMI(x,y,t) are determined from these filtergrams using a simplified version of the HMI pipeline. The Doppler velocities are correlated with the original velocities in the simulated atmosphere. The cross-correlation peaks near 100 km, suggesting that the HMI Doppler velocity signal is formed rather low in the solar atmosphere. The same analysis is performed for the SOHO/MDI Ni?i line at 6768 Å. The MDI Doppler signal is formed slightly higher at around 125 km. Taking into account the limited spatial resolution of the instruments, the apparent formation height of both the HMI and MDI Doppler signal increases by 40 to 50 km. We also study how uncertainties in the HMI filter-transmission profiles affect the calculated velocities.     

Correlations, spectra, and instability of phase-space density fluctuations in open-cluster models

The dynamical evolution of six open star cluster models is analyzed using the correlation and spectral analysis of phase-space density fluctuations. The two-time and mutual correlation functions are computed for the fluctuations of the phase-space density of cluster models. The data for two-time and two-particle correlations are used to determine the correlation time for phase-space density fluctuations ((0.1–1) τ _v.r., where τ _v.r. is the violent relaxation time of the model) and the average phase velocities of the propagation of such fluctuations in cluster models. These velocities are 2–20 times smaller than the root mean square velocities of the stars in the cluster core. The power spectra and dispersion curves of phase-space density fluctuations are computed using the Fourier transform of mutual correlation functions. The results confirm the presence of known unstable phase-space density fluctuations due to homologous fluctuations of the cluster cores. The models are found to exhibit a number of new unstable phase-space density fluctuations (up to 32–41 pairs of fluctuations with different complex conjugate frequencies in each model; the e-folding time of the amplitude growth of such fluctuations is (0.4–10) τ _v.r. and their phases are distributed rather uniformly). Astrophysical applications of the obtained results (irregular structure of open star clusters, formation and decay of quasi-stationary states in such clusters) are discussed.     

On the dynamics of the jovian ionosphere and thermosphere.: IV. Ion-neutral coupling

We use the fully coupled, three-dimensional, global circulation Jovian Ionospheric Model (JIM) to calculate the coupling between ions in the jovian auroral ovals and the co-existing neutral atmosphere. The model shows that ions subject to drift motion around the auroral oval, as a result of the E×B coupling between a meridional, equatorward electric field and the jovian magnetic field, generate neutral winds in the planetary frame of reference. Unconstrained by the magnetic field, these neutral winds have a greater latitudinal extent than the corresponding ion drifts. Values of the coupling coefficient, k(h), are presented as a function of altitude and cross-auroral electric field strength, for different incoming electron fluxes and energies. The results show that, with ion velocities of several hundred metres per second to over 1 km s⁻¹, k(h) can attain values greater than 0.5 at the ion production peak. This parameter is key to calculating the effective conductivities required to model magnetosphere-ionosphere coupling correctly. The extent to which angular momentum (and therefore energy) is transported vertically in JIM is much more limited than earlier, one-dimensional, studies have predicted.     

Plasma drift in the solar corona

Model calculations of plasma drifts in the solar corona were performed. We established that only drifts in crossed fields could result in velocities V of several hundred kilometers per second. Such velocities are typical of coronal mass ejections (CMEs). We derived an analytic expression for V where n, the expansion harmonic of the magnetic-field strength, varies with time. As follows from this expression, V is a power function of the distance with index (2?n) and the radial component changes sign (n?1) times in the latitude range from ?π/2 to +π/2. We found that if the magnetic dipole moment varies with time, the similarity between the spiral structures of coronal plasma is preserved when they displace within several solar radii and the density gradient at the conical boundaries increases (the apparent contrast is enhanced). There is a correspondence between the inferred model effects and the actually observed phenomena that accompany CMEs.