Numerical simulation of unstable two-dimensional motions of a circumstellar shell

The growth of two-dimensional axisymmetric perturbations in the motion of a neutral shell formed in the interstellar medium when an ionization-shock front exits at the surface of a cloud is simulated numerically. The perturbations are assumed to emerge when the shock ahead of the ionization front reaches the cloud boundary. For long-wavelength perturbations, the accumulation of mass has been found to take place in radially oriented condensations in the shape of a rod pointed toward the star and widened at the opposite end as a result of instability. The shell fragmentation is accompanied by supersonic spouting of a hot plasma into a low-density medium. Flow nonstationarity is shown to affect significantly the gas density and velocity distributions both inside and in the immediate vicinity of the condensation. As one recedes from the ionization front, the density of charged particles changes only slightly, which is inconsistent with the power law of density decrease with increasing distance from the condensation center commonly used in interpreting observations.     

Instability of a thin photoevaporable circumstellar envelope

We investigate the stability of a dense neutral shell that is accelerated outward by the hot-gas pressure and that loses its mass through photoionization by radiation from the central star. We assume the H I shell to be thin and use the Lagrangian coordinates to describe its motion. We show that a flow accompanied by cumulative effects emerges during the nonlinear development of the instability. We estimate the influence of the radiative cooling rate on the motion and determine parameters of the gas in the cumulative region. The results obtained are compared with the observations of the nebulae NGC 7293 and NGC 2392.     

Evolution and Fragmentation of Wide-Angle Wind Driven Molecular Outflows

We present two dimensional cylindrically symmetric hydrodynamic simulations and synthetic emission maps of a stellar wind propagating into an infalling, rotating environment. The resulting outflow morphology, collimation and stability observed in these simulations have relevance to the study of young stellar objects, Herbig-Haro jets and molecular outflows. Our code follows hydrogen gas with molecular, atomic and ionic components tracking the associated time dependent molecular chemistry and ionization dynamics with radiative cooling appropriate for a dense molecular gas. We present tests of the code as well as new simulations which indicate the presence of instabilities in the wind-blown bubble’s swept-up shell.     

A WENO algorithm of the temperature and ionization profiles around a point source

We develop a numerical solver for radiative transfer problems based on the weighted essentially nonoscillatory (WENO) scheme modified with anti-diffusive flux corrections, in order to solve the temperature and ionization profiles around a point source of photons in the reionization epoch. Algorithms for such simulation must be able to handle the following two features: (1) the sharp profiles of ionization and temperature at the ionizing front (I-front) and the heating front (T-front), and (2) the fraction of neutral hydrogen within the ionized sphere is extremely small due to the stiffness of the rate equations of atom processes. The WENO scheme can properly handle these two features, as it has been shown to have high order of accuracy and good convergence in capturing discontinuities and complicated structures in fluid as well as to be significantly superior over piecewise smooth solutions containing discontinuities. With this algorithm, we show the time-dependence of the preheated shell around a UV photon source. In the first stage the I-front and T-front are coincident, and propagate with almost the speed of light. In later stage, when the frequency spectrum of UV photons is hardened, the speeds of propagation of the ionizing and heating fronts are both significantly less than the speed of light, and the heating front is always beyond the ionizing front. In the spherical shell between the I- and T-fronts, the IGM is heated, while atoms keep almost neutral. The time scale of the preheated shell evolution is dependent on the intensity of the photon source. We also find that the details of the pre-heated shell and the distribution of neutral hydrogen remained in the ionized sphere are actually sensitive to the parameters used. The WENO algorithm can provide stable and robust solutions to study these details.     

Matter infall in collapsing molecular cloud cores with an axial magnetic field

The magnetic fields affect collapse of molecular cloud cores. Here, we consider a collapsing core with an axial magnetic field and investigate its effect on infall of matter and formation of accretion disk. For this purpose, the equations of motion of ions and neutral infalling particles are numerically solved to obtain the streamlines of trajectories. The results show that in a non-steady state of ionization and ion–neutral coupling, which is not unexpected in the case of infall, the radius of accretion disk will be larger as a consequence of axial magnetic field.     

Cosmic ray ionization of lower Venus atmosphere

Cosmic ray particles passing through dense lower atmosphere of Venus decay giving rise to various charged and neutral particles. The flux and degradation of dominant cascade particles namely neutrinos and pions are computed and ionization contributions at lower altitudes are estimated. Using the height profile of pion flux, the muon flux is computed and used to estimate ionization at lower altitudes. It is shown that cosmic ray produced ionization descends to much lower altitudes intercepting the thickness of Venus cloud deck. The dynamical features of Venus cloud deck are used to allow the likely charging and charge separation processes resulting into cloud-to-cloud lightning discharges.     

Nonlinear Dynamics of Ionization Fronts in HII Regions

Hydrodynamic instability of an accelerating ionization front (IF) is investigated with 2D hydrodynamic simulations, including absorption of incident photoionizing photons, recombination in the HII region, and radiative molecular cooling. When the amplitude of the perturbation is large enough, nonlinear dynamics of the IF triggered by the separation of the IF from the cloud surface is observed. This causes the second harmonic of the imposed perturbation to appear on the cloud surfaces, whereas the perturbation in density of ablated gas in the HII region remains largely single mode. This mismatch of modes between the IF and the density perturbation in the HII region prevents the strong stabilization effect seen in the linear regime. Large growth of the perturbation caused by Rayleigh-Taylor-like instability is observed late in time.     

Oscillatory propagation of an ionization-shock front

Based on numerical simulations, we show that the oscillatory propagation of a plane-parallel ionization-shock front is possible for a typical dependence of the interstellar-medium cooling function on density and temperature. In this case, the oscillation amplitude of the shock position in the presence of an ionization front can be several times larger than its value for a single shock wave. The variations in neutral and ionized gas velocities attributable to oscillations are comparable in order of magnitude and agree with the random velocities observed in H II regions.     

Effects of cosmic rays and density inhomogeneity of the circumstellar medium on the formation of a supernova shell

We consider the self-similar problem of a supernova explosion in a radially inhomogeneous medium by taking into account the generation of accelerated relativistic particles. The initial density of the medium is assumed to decrease with distance from the explosion center as a power law, ρ ₀ = A/r ^θ. We use a two-fluid approach in which the total pressure in the medium is the sum of the circumstellar gas pressure and the relativistic particle pressure. The relativistic particle pressure at the shock front is specified as an external parameter. This approach is applicable in the case where the diffusion coefficient of accelerated particles is small and the thickness of the shock front is much smaller than its radius. We have numerically solved a system of ordinary differential equations for the dimensionless quantities that describe the velocity and density behind the shock front as well as the nonrelativistic gas and relativistic particle pressures for various parameters of the inhomogeneity of the medium and various compression ratios of the medium at the shock front. We have established that the shock acceleration of cosmic rays affects most strongly the formation of a supernova shell (making it thinner) in a homogeneous circumstellar medium. A decrease in the circumstellar matter density with distance from the explosion center causes the effect of shock-accelerated relativistic particles on the supernova shell formation to weaken considerably. Inhomogeneity of the medium makes the shell thicker and less dense, while an increase in the compression ratio of the medium at the shock front causes the shell to become thinner and denser. As the relativistic particle density increases, the effect of circumstellar matter inhomogeneity on the shell formation becomes weaker.     

Charge transfer reactions at interfaces between neutral gas and plasma: Dynamical effects and X‐ray emission

Charge‐transfer is the main process linking neutrals and charged particles in the interaction regions of neutral (or partly ionized) gas with a plasma. In this paper we illustrate the importance of charge‐transfer with respect to the dynamics and the structure of neutral gas‐plasma interfaces. We consider the following phenomena: (1) the heliospheric interface ‐ region where the solar wind plasma interacts with the partly‐ionized local interstellar medium (LISM) and (2) neutral interstellar clouds embedded in a hot, tenuous plasma such as the million degree gas that fills the so‐called “Local Bubble”. In (1), we discuss several effects in the outer heliosphere caused by charge exchange of interstellar neutral atoms and plasma protons. In (2) we describe the role of charge exchange in the formation of a transition region between the cloud and the surrounding plasma based on a two‐component model of the cloud‐plasma interaction. In the model the cloud consists of relatively cold and dense atomic hydrogen gas, surrounded by hot, low density, fully ionized plasma. We discuss the structure of the cloud‐plasma interface and the effect of charge exchange on the lifetime of interstellar clouds. Charge transfer between neutral atoms and minor ions in the plasma produces X‐ray emission. Assuming standard abundances of minor ions in the hot gas surrounding the cold interstellar cloud, we estimate the X‐ray emissivity consecutive to the charge transfer reactions. Our model shows that the charge‐transfer X‐ray emission from the neutral cloud‐plasma interface may be comparable to the diffuse thermal X‐ray emission from the million degree gas cavity itself (© 2012 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Gravitational interaction between a planet and an optically thin disc

We investigate the gravitational interaction between a planet and an optically thin protoplanetary disc, performing local three-dimensional hydrodynamical simulations. In the present study, we take account of radiative energy transfer in optically thin discs. Before the stage of planetary accretion, dust opacity is expected to decrease significantly because of grain growth and planetesimal formation. Thus, it would be reasonable to consider optically thin discs in the disc–planet interaction. Furthermore, we focus on small planets that can neither capture disc gas nor open a disc gap. The one-sided torque exerted on a planet by an optically thin disc is examined for various values of the disc optical thickness (<1). In optically thin discs, the temperature behind the density waves is lower than the unperturbed value because of radiative cooling. Heating due to shock dissipation is less effective than radiative cooling. Because of radiative cooling, the density distribution around the planet is not axisymmetric, which exerts an additional torque on the planet. The torque enhancement becomes maximum when the cooling time is comparable with the Keplerian period. The enhancement is significant for low-mass planets. For planets with  3 M_⊕  , the additional one-sided torque can be 40 per cent of the torque in the isothermal case. The radiative cooling is expected to change the differential torque and the migration speed of planets, too.     

A possible EGRET source with a two-sided radio jet structure on a parsec scale – 3EG J0808+5114

Recent observations have revealed that damped Lyα clouds (DLAs) host star formation activity. In order to examine if such star formation activity can be triggered by ionization fronts, we perform high-resolution hydrodynamics and radiative transfer simulations of the effect of radiative feedback from propagating ionization fronts on high-density clumps. We examine two sources of ultraviolet (UV) radiation field to which high-redshift ( z ∼ 3) galaxies could be exposed: one corresponding to the UV radiation originating from stars within the DLA, itself, and the other corresponding to the UV background radiation. We find that, for larger clouds, the propagating I-fronts created by local stellar sources can trigger cooling instability and collapse of significant part, up to 85 per cent, of the cloud, creating conditions for star formation in a time-scale of a few Myr. The passage of the I-front also triggers collapse of smaller clumps (with radii below ∼4 pc), but in these cases the resulting cold and dense gas does not reach conditions conducive to star formation. Assuming that 85 per cent of the gas initially in the clump is converted into stars, we obtain a star formation rate of  ∼0.25 M_⊙ yr⁻¹ kpc⁻²  . This is somewhat higher than the value derived from recent observations. On the other hand, the background UV radiation which has harder spectrum fails to trigger cooling and collapse. Instead, the hard photons which have long mean free-path heat the dense clumps, which as a result expand and essentially dissolve in the ambient medium. Therefore, the star formation activity in DLAs is strongly regulated by the radiative feedback, both from the external UV background and internal stellar sources and we predict quiescent evolution of DLAs (not starburst-like evolution).     

Continuum shielding and flow dynamics in active galactic nuclei

We study the ionization, thermal structure and dynamics of active galactic nuclei (AGN) flows that are partially shielded from the central continuum. We utilize a detailed non-local thermodynamic equilibrium photoionization and radiative transfer code using exact (non-Sobolev) calculations. We find that shielding has a pronounced effect on the ionization, thermal structure and the dynamics of such flows. Moderate shielding is especially efficient in accelerating flows to high velocities because it suppresses the ionization level of the gas. The ionization structure of shielded gas tends to be distributed uniformly over a wide range of ionization levels. In such gas, radiation pressure due to trapped line photons can dominate over the thermal gas pressure and have a significant effect on the thermal stability of the flow. Heavily shielded flows are driven mainly by line radiation pressure, and so line locking has a large effect on the flow dynamics. We show that the observed 'L_α ghost' is a natural outcome in highly ionized flows that are shielded beyond the Lyman limit. We suggest that high-velocity AGN flows occupy only a small fraction of the volume and that their density depends only weakly on the velocity field.     

Charge state behaviour of projectiles under an acceleration mechanism in a hot plasma

We have studied the ionization states of most ions in solar flares when a stochastic acceleration mechanism is present. Calculations using a computer code (ESCAPE) designed to find the charge states of heavy and light ions under stochastic acceleration have shown an energy dependence on the ionization states, stronger for heavy ions. Moreover the charge states depend on source parameters as density or temperature, but if an acceleration mechanism is taken into account in impulsive solar energetic particles events, the same source temperature may not be inferred for all the ions.     

Revealing the “fingerprints” of the magnetic precursor of C-shocks

We present the first C-shock and radiative transfer model that calculates the evolution of the line profiles of neutral and ion species like SiO, H¹³CO⁺ and HN¹³C for different flow times along the propagation of the shock through the unperturbed gas. We find that the line profiles of SiO characteristic of the magnetic precursor stage have very narrow linewidths and are centered at velocities close to the ambient cloud velocity, as observed toward the young shocks in the L1448-mm outflow. Consistently with previous works, our model also reproduces the broad SiO emission detected in the high velocity gas in this outflow, for the downstream postshock gas in the shock. This implies that the different velocity components observed in L1448-mm are due to the coexistence of different shocks at different evolutionary stages.     

Expansion velocities of the planetary nebula NGC 7009—the high dispersion spectra line profiles

We investigated the line profiles of the high dispersion spectra of He?I, He?II, N?III, [N?II], [O?III], [Ar?III], [Ar?IV], [S?II], [S?III ], and [Cl?IV], secured at the center of the planetary nebula NGC 7009 with the fiber-fed Bohyunsan Echelle Spectrograph (BOES). The expansion velocities for the main shell and the faint outer thin shell were derived based on stronger double Gaussian profiles and the fast blue wing component, respectively. With the Keck 2D spectral images, we set the main shell boundaries as R～4″ and 6″ and the faint outer thin shell boundaries as R～11″ and 13″. The radial variation of expansion velocities of both shells does not follow a linear acceleration. The acceleration of the main shell gas seems to be in a retarded mode, while that of the outer faint shell gas seems to be affected by an additional force. The fast wing component appears also in the He?I, He?II, and N?III line profiles which are likely to have been formed in a compact region, located inside the main shell. The additional deceleration factor, such as the rotation of the central star or the dust gas in the gas shell might cause the non-linear acceleration. Based on the line profiles responsible for the three zones, we conclude that there were at least three major changes occurred in the central star temperature since the first outflow from the central star.     

Structure and Instabilities of Irradiated Protoplanetary Disks: Effects of Dust Particle Growth

We have constructed self-consistent temperature and density profiles of irradiated active protoplanetary disks, using a two-dimensional radiative transfer calculation. By means of these profiles we have studied the stabilization of the convective instability by radiative heating and the magnetorotational instability (MRI) via ohmic dissipation, taking into account the effect of dust particle growth. Simple chemistry such as ionization by cosmic rays and recombination on dust grains are used to calculate the ionization degree of gas in the disks. Our results show that the dust growth stabilizes the convective instability due to the 2D effect of radiative transfer, while it enhances the MRI through the decrease in the recombination of ions on the dust grains. In addition, the influences of the dust settling toward the midplane of the disks on the instabilities are discussed.     

A cloud collision model for water maser excitation

High-velocity collisions between small, dense, neutral clouds or between a dense cloud and a dense shell can provide the energy source required to excite H2O maser emission. The radiative precursor from the surface of the collisional shock front rapidly diffuses through the cloud, heating the dust grains but leaving the H2 molecules cool. Transient maser emission occurs as the conditions for the Goldreich and Kwan "hot-dust cold-gas" maser pump scheme are realized locally within the cloud. In time the local maser action quenches due to the heating of the H2 molecules by collisions against the grains. Although this model cannot explain the very long-lived steady maser features, it is quite successful in explaining a number of the observed properties of the high-velocity features in such sources as Orion, W51, and W49. In particular, it provides a natural explanation for the rapid time variations, the narrow line widths, juxtaposition of high- and low-velocity features, and the short lifetimes which are frequently observed for the so-called high-velocity maser "bullets" thought to be accelerated by strong stellar winds.