Self-similar regimes for the collapse of magnetic protostellar clouds

Self-similar solutions describing the homogeneous, free, isothermal collapse of protostellar clouds are considered. One such solution correponds to the critical case of the propagation of a rarefaction wave near the time when it is focused in the central region of the cloud. The speed of the rarefaction front is finite and equal to three times the isothermal sound speed. The asymptotic distributions of gas-dynamical quantities in the central part of the collapsing cloud and the surrounding envelopes are considered at both early and late stages of compression, after the formation of an opaque core (protostar). These solutions are used in a magneto-kinematic approximation to study the geometry and evolution of the large-scale magnetic fields of collapsing protostellar clouds. All the solutions are verified using direct numerical simulations. It is shown that an initially uniform magnetic field acquires an “hourglass” geometry with time. The characteristic opening angles in the self-similar solutions are in satisfactory agreement with observations.     

The nonthermal radio emission of WR 140 in a model with colliding two-phase winds

We present a model in which the nonthermal radio emission of binary systems containing Wolf-Rayet and O components is due to collisions between clouds belonging to dense phases of the wind of each star. The relativistic electrons are generated during the propagation of fast shock waves through the clouds and their subsequent de-excitation. The initial injection of superthermal particles is due to photoionization of the de-excited cold gas by hard radiation from the shock front. Therefore, the injection takes place in cloud regions fairly far from the front. Further, the superthermal electrons are accelerated by the betatron mechanism to relativistic energies during the isobaric compression of the cloud material, when most of the gas radiates its energy. Collisions between the clouds can occur far beyond the contact boundary between the rarefied wind components. Thus, the model avoids the problem of strong low-frequency absorption of the radiation.     

The dynamics of the shock disruption of interstellar clouds

The effect of the structural irregularity of an interstellar cloud on the dynamics of its disruption by a shock from a supernova is studied. Irregular clouds are disrupted twice as fast as spherical clouds. However, fragments of irregular clouds preserve their enhanced density for long times without being mixed with the intercloud gas. The fraction of shock energy that is converted to the kinetic energy of the fragments is 50% higher than in the disruption of a spherical cloud. Shocks are not able to trigger the gravitational compression of clouds.     

空气间隔装药爆破机理研究

利用爆轰波理论分析了空气间隔装药炮孔内一维不定常激波的相互作用及其在炮孔堵头、孔底的反射过程,同时分析了孔内各点的压力随时间的变化过程,介绍了空气间隔装药爆破的机理及设计参数。基于此,认为应充分利用空气间隔爆破结构的优势,并在梯段爆破中满足以下两个条件：（1）在设计过程中要尽量使稀疏波及从孔底反射的稀疏波传播过程能在整个孔内每一断面都作用到,即稀疏波到达孔底的时间要比从堵头反射的压力波到达孔底要早;（2）反射压力波应该到达空气与爆生气体接触面的时间比从孔底反射的稀疏波到达空气与爆生气体接触面的时间要早。由此通过计算得到了在梯段爆破工程中合理的空气层长度比例值约为30 %～42 %。计算结论与已有实测成果基本一致。     

Formation of low-mass X-ray novae with black holes from triple systems

We apply a population synthesis technique to study the formation and evolution of low-mass X-ray binaries with black holes, observed as X-ray novae, from hierarchical triple systems. A scenario is suggested in which an inner close binary system evolves into an X-ray system with a large mass ratio. The high rate of accretion onto the neutron star leads to a common envelope stage, which may result in the formation of a Thorne-Zytkow (TZ) object. During its evolution, the envelope of the TZ object expands, encompassing the third star. The recurrent common-envelope stage decreases the size of the orbit of the third star, leading to the formation of a lowmass X-ray nova with a black hole. The dynamical stability of triple systems automatically ensures that only lowmass X-ray novae form. We also consider the possible formation of an X-ray nova from a binary in the case of asymmetrical core collapse during a supernova explosion.     

Modeling of rapid variability in the spectral line profiles of Wolf-Rayet stars

Profiles of variable emission lines in the spectra of Wolf-Rayet stars are calculated using a stochastic cloud model for the inhomogeneous atmospheres of early-type stars. The model assumes that most of the line flux is formed in cold, dense condensations (clouds) that move through a rarified inter-cloud medium whose density monotonically decreases outwards. The formation of clouds is taken to be stochastic. Wavelet analysis is used to estimate the parameters of cloud ensembles. The model can reproduce the general pattern of line-profile variability observed in the spectra of Wolf-Rayet stars.     

Formation of close binary black holes merging due to gravitational-wave radiation

The conditions for the formation of close-binary black-hole systems merging over the Hubble time due to gravitational-wave radiation are considered in the framework of current ideas about the evolution of massive close-binary systems. The original systems whose mergers were detected by LIGO consisted of main-sequence stars with masses of 30–100M _⊙. The preservation of the compactness of a binary black hole during the evolution of its components requires either the formation of a common envelope, probably also with a low initial abundance of metals, or the presence of a “kick”—a velocity obtained during a supernova explosion accompanied by the formation of a black hole. In principle, such a kick can explain the relatively low frequency of mergers of the components of close-binary stellar black holes, if the characteristic speed of the kick exceeds the orbital velocities of the system components during the supernova explosion. Another opportunity for the components of close-binary systems to approach each other is related to their possible motion in a dense molecular cloud.     

Analysis of Cloud Scale Error of Low Resolution Satellite Based on GF-4 and Its Influence on Downward Radiation Calculation

In order to study the scale error of low resolution meteorological satellite cloud detection and its impact on the calculation of downlink radiation, cloud detection using high resolution stationary satellite GF-4 data and error analysis were carried out. Firstly, the cloud detection of GF-4 data is carried out by using visible channel threshold method and time series method, and the error of cloud detection results of Himawari-8 and FY-2 (FY-2G, FY-2E) is analyzed based on the results of GF-4 cloud detection.In the study area, FY-2G, FY-2E and Himawari-8 cloud images could distinguish the clouds and clear sky. The main reason for the error was the scale effect produced by different spatial resolution satellites(the differences caused by cloud detection algorithms are not discussed here).Most of the errors occurred in the areas of thin clouds and broken clouds.High resolution data could detect broken clouds, while low resolution data lead to false and missed detection. On this basis, the error of remote sensing calculation of short wave radiation was analyzed,and it was found that the error of the actual cloud amount in the pixel would bring significant error to the estimation of the downward radiation.The relative error of the instantaneous downward radiation in the selected test area was -173.52%, and the maximum relative error of shortwave radiation was -20.20%.The results show that the high resolution stationary satellite data can significantly improve the estimation accuracy of the downlink shortwave radiation in the regions with more broken clouds.     

Turning points in black holes astrophysics

Black hole astrophysics is expected to have a major breakthrough any day. Theoretical understanding is well advanced, so the breakthrough has to come from observers. Conclusive and direct evidence of event horizons, signatures of tidal deformation of gas clouds (and eventual disappearance of matter behind horizons) could be obtained in matter of months. At this juncture, we wish to summarize how the subject came to this stage and what is the state-of-the-art flow solution around a black hole. We touch upon the most crucial issue of (un)predictability of this subject. We believe that if viscosity is roughly what MRI simulations are pointing to us, we still require very detailed knowledge of properties of the companion, especially time dependent mass transfer from the companion to the black hole.     

Circumstellar nebulae of WR and OB stars with strong winds: Simultaneous excitation by stellar radiation and shocks

We use a two-phase model for the structure of the circumstellar nebulae of hot stars to analyze the radiative cooling of a dense, compact cloud behind the shock produced by the compression of the cloud by hot gas from the stellar wind, taking into account ionization and heating by radiation from the central star. We can distinguish three stages of the evolution of the cloud during its compression. In the first stage, relevant for the entire cloud before compression and the gas ahead of the shock front, the state of the gas is determined purely by ionization by the stellar radiation. The next stage is characterized by the simultaneous action of two gas excitation mechanisms—photoionization by the stellar radiation and shock heating. In this stage, the gas intensively radiates thermal energy received at the shock front. After radiative cooling, in the final stage, ionization and heating of the gas are again determined mainly by the star. To compute the spectrum of the cloud radiation, we solved for the propagation of a plane-parallel, homogeneous flux through the shock front in the radiation field of the hot star. The computations show that a combination of two excitation mechanisms considerably enriches the theoretical spectrum. The relative intensities of emission lines of a single cloud may resemble either those for an HII region or of a supernova remnant.     

Dust formation in binaries with OB and WR components in a two-phase stellar-wind model

We describe the formation of carbon dust in binary systems with hot components as a result of the collisions of clouds in a two-phase stellar-wind model. Calculations are made for the well studied system WR 140. The collisions lead to the formation of composite clouds and shock waves, with the temperature at the shock front equal to about 3×10⁸ K along both sides of the interface boundary. During isobaric deexcitation to (0.5–0.7)×10⁴ K, the cloud density increases by a factor of several thousand; its thickness in the direction of the shock decreases by the same factor. After deexcitation, the hydrogen inside the composite cloud is in its atomic state, while the carbon remains ionized. The deexcitation is followed by expansion of the cloud, which moves away from both stars. During the first 10⁶ s, its thickness remains relatively small, so that the expansion is one-dimensional. The radiation field inside the cloud decays, resulting in the recombination of the carbon. Further expansion of the cloud leads to adiabatic cooling, and the formation of dust particles becomes possible. After the dimensions of the cloud have become roughly the same in all directions, its expansion is isotropic, so that it becomes transparent within approximately 10⁶ s, and the dust is heated to (1.0–1.4)×10³ K, observed as an IR “lare.” The time required for the cloud to move from the exciting star and heat the dust is comparable to the observed delay in the increased IR emission relative to the time of periastron.     

A multicomponent model for computing the thermal structure of collapsing protostellar clouds

A model for simulating the thermal and dynamical evolution of protostellar clouds is presented. In the model, the dust and gas temperatures are treated separately, making it possible to more precisely describe the initial stages of the cloud’s gravitational contraction and collapse. The model is fast enough to be applied in hydrodynamical computations, and has a high enough accuracy for the results to be used to compute emission spectra and comparing them with observational data. Two problems are considered as test examples and simple applications: calculation of the structure of clouds in thermal and hydrostatic equilibrium, and modeling the evolution of a protostellar cloud in a spherically symmetric approximation, including the formation of the first hydrostatic core.     

VLBI of supernovae and gamma-ray bursts

Supernovae and gamma-ray bursts (GRBs) are among the brightest events in the universe. Excluding Type Ia supernovae and short GRBs, they are the result of the core collapse of a massive star with material being ejectedwith speeds of several 1000 km/s to nearly the speed of light, and with a neutron star or a black hole left over as the compact remnant of the explosion. Synchrotron radiation in the radio is generated in a shell when the ejecta interact with the surrounding medium and possibly also in the central region near the compact remnant itself. VLBI has allowed resolving some of these sources and monitoring their expansion in detail, thereby revealing characteristics of the dying star, the explosion, the expanding shock front, and the expected compact remnant. We report on updates of some of the most interesting results that have been obtained with VLBI so far. Movies of supernovae are available from our website. They show the evolution from shortly after the explosion to decades thereafter, in one case revealing an emerging compact central source, which may be associated with shock interaction near the explosion center or with the stellar corpse itself, a neutron star or a black hole.     

The mechanism of circumbinary envelope formation in close binaries

We consider the structure and formation of the circumbinary envelopes in semi-detached binary systems. Three-dimensional numerical simulations of the gas dynamics are used to study the flow pattern in a binary system after it has reached the steady-state accretion regime. The outer parts of the circumbinary envelope are replenished by periodic ejections from the accretion disk and circum-disk halo through the vicinity of the Lagrange point L₃. In this mechanism, the shape and position of a substantial part of the disk is specified by a precessional density wave. On timescales comparable to the orbital period, the precessional wave (and hence an appreciable fraction of the disk) will be virtually stationary in the observer’s frame, whereas the positions of other elements of the flow will vary due to the orbital rotation. The periodic variations of the positions of the disk and the bow shock formed when the inner parts of the circumbinary envelope flow around the disk result in variations in both the rate of angular-momentum transfer to the disk and the flow structure near L₃. All these factors lead to a periodic increase of the matter flow into the outer layers of the circumbinary envelope through the vicinity of L₃. The total duration of the ejection is approximately half the orbital period.     

Guest Editorial to the Special Issue: Lessons Learned from post-2008 Wenchuan Earthquake Community Recovery

Identifying the spatial extent of volcanic ash clouds in the atmosphere and forecasting their direction and speed of movement has important implications for the safety of the aviation industry, community preparedness and disaster response at ground level. Nine regional Volcanic Ash Advisory Centres were established worldwide to detect, track and forecast the movement of volcanic ash clouds and provide advice to en route aircraft and other aviation assets potentially exposed to the hazards of volcanic ash. In the absence of timely ground observations, an ability to promptly detect the presence and distribution of volcanic ash generated by an eruption and predict the spatial and temporal dispersion of the resulting volcanic cloud is critical. This process relies greatly on the heavily manual task of monitoring remotely sensed satellite imagery and estimating the eruption source parameters (e.g. mass loading and plume height) needed to run dispersion models. An approach for automating the quick and efficient processing of next generation satellite imagery (big data) as it is generated, for the presence of volcanic clouds, without any constraint on the meteorological conditions, (i.e. obscuration by meteorological cloud) would be an asset to efforts in this space. An automated statistics and physics-based algorithm, the Automated Probabilistic Eruption Surveillance algorithm is presented here for auto-detecting volcanic clouds in satellite imagery and distinguishing them from meteorological cloud in near real time. Coupled with a gravity current model of early cloud growth, which uses the area of the volcanic cloud as the basis for mass measurements, the mass flux of particles into the volcanic cloud is estimated as a function of time, thus quantitatively characterising the evolution of the eruption, and allowing for rapid estimation of source parameters used in volcanic ash transport and dispersion models.

     

The mass spectrum of black holes in close binary systems

We present the results of population syntheses obtained using our “scenario machine.” The mass spectra of black holes in X-ray binary systems before and after the stage of accretion from an optical companion are obtained for various evolutionary scenarios. The results of the model computations are compared to observational data. The observational data are used to estimate the fraction of a presupernova’s mass that collapses into a black hole. This model can explain the formation of low-mass (2–4M_⊙) black holes in binary systems with optical companions. We show that the number of low-mass black holes in the Galaxy is sufficiently high for them to be detected. The population-synthesis results suggest that the vast majority of low-mass black holes are formed via the accretion-induced collapse of neutron stars. The percentage of low-mass black holes in binary systems that form due to accretion-induced collapse is 2–15% of the total number of black holes in binaries, depending on the evolutionary scenario.     

The stability of an isothermal gaseous sphere with gravitation

A dispersion relation has been obtained for waves propagating in an inhomogeneous gas with viscosity and gravity in the direction of the density gradient. It has been shown that a stationary distribution of isothermal gas is unstable under gravitational—convective perturbations. It is suggested that the initialand boundary-value problems of the non-stationary dynamics of an isothermal gas with gravity but without viscosity are incorrect and require regularization, for example, by adding viscosity. The presence of viscosity does not prevent the development of the Jeans long-wave instability. Numerical experiments based on a fully three-dimensional, non-stationary numerical code have demonstrated the development of a collapse inside an isothermal cloud.     

不同应力路径下剪切带的数值模拟

采用回映应力更新算法,编写了基于伏斯列夫面的超固结黏土本构关系模型子程序,嵌入非线性有限元软件ABAQUS。通过对单元试验进行三轴压缩、三轴伸长及平面应变等问题的模型预测,再现了超固结黏土在不同初始超固结比和应力路径时的变形和强度特性,从而验证了子程序的正确性。借助该本构模型,对三轴压缩、三轴伸长及平面应变应力路径下超固结黏土体变形局部化问题,进行了三维数值模拟。分析结果表明：超固结黏土在三轴压缩及伸长状态时,土体变形局部化在应力-应变关系软化时出现,而平面应变状态时,在应力-应变关系硬化阶段出现,其超固结黏土的剪胀特性在剪切带的形成过程中起重要作用。     

The triggered origin of the solar system

The scenario of the triggered origin of the solar system suggests that the formation of our planetary system was initiated by the impact of an interstellar shock wave on a molecular cloud core. The strength of this scenario lies in its ability to explain the presence of short-lived radionuclides in the early solar system. According to the proposal, the radioactivities were produced in a stellar source, transported into the molecular cloud core by a shock wave and mixed into the collapsing system during the interaction between the shock wave and the core. We examine the viability of the scenario by presenting results from recent numerical simulations. The calculations show that molecular cloud cores can be triggered into collapse by the impact of a shock wave propagating at the velocity of 10–45 km s⁻¹. Some of the shock wave material incident on the core, typically 10–20%, can be injected into the collapsing system. The time scale of the process is ∼10⁴–10⁵ years, sufficiently short for the survival of the short-lived radioactivities. The simulations therefore confirm the viability of the scenario of the triggered origin of the solar system.