Radiation of accretion disks in quasars and galactic nuclei

The continuum spectra of accretion disks around supermassive Schvarzshild black hole are calculated for a set of accretion disk parameters. In every point of the disk surface the atmosphere is considered to be similar to the atmosphere of the star with the same values of the surface gravity and effective temperature. All important effects of general relativity affected the spectrum shape are taken into account.     

2D and 3D MHD simulations of disk accretion by rotating magnetized stars: Search for variability

We performed 2D and full 3D magnetohydrodynamic simulations of disk accretion to a rotating star with an aligned or misaligned dipole magnetic field. We investigated the rotational equilibrium state and derived from simulations the ratio between two main frequencies: the spin frequency of the star and the orbital frequency at the inner radius of the disk. In 3D simulations we observed different features related to the non-axisymmetry of the magnetospheric flow. These features may be responsible for high-frequency quasi-periodic oscillations (QPOs). Variability at much lower frequencies may be connected with restructuring of the magnetic flux threading the inner regions of the disk. Such variability is specifically strong at the propeller stage of evolution.     

X-ray observations of black-hole accretion disks

Our current theoretical and observational understandings of the accretion disks around Galactic black-holes are reviewed. Historically, a simple phenomenological accretion disk model has been used to interpret X-ray observations. Although such a phenomenological interpretation is still useful, high quality X-ray data from contemporary instruments allow us to test more realistic accretion disk models. In a simple and ideal case, the standard optically thick accretion disk model is successful to explain observations, such that the inner disk radius is constant at three times the Schwarzschild radius over large luminosity variations. However, when disk luminosity is close to or exceeds the Eddington luminosity, the standard disk model breaks, and we have to consider the “slim disk” solution in which radial energy advection is dominant. Recent observations of Ultra-luminous X-ray sources (ULXs), which may not be explained by the standard disk model, strongly suggest the slim disk solution. We compare theoretical X-ray spectra from the slim disk with observed X-ray spectra of ULXs. We have found that the slim disk model is successful to explain ULX spectra, in terms of the massive stellar black-holes with several tens of solar mass and the super-Eddington mass accretion rates. In order to explain the large luminosities (>10⁴⁰ ergs s⁻¹) of ULXs, “intermediate black-holes” (>100M) are not required. Slim disks around massive stellar black-holes of up to several tens of solar mass would naturally explain the observed properties of ULXs.     

Galactic outflows and the pollution of the galactic environment by supernovae

We here explore the effects of the SN explosions into the halo of star-forming galaxies like the Milky Way. Successive randomly distributed and clustered SNe explosions cause the formation of hot superbubbles that drive either fountains or galactic winds above the galactic disk, depending on the amount and concentration of energy that is injected by the SNe. In a galactic fountain, the ejected gas is re-captured by the gravitational potential and falls back onto the disk. From 3D non-equilibrium radiative cooling hydrodynamical simulations of these fountains, we find that they may reach altitudes up to about 5 kpc in the halo and thus allow for the formation of the so called intermediate-velocity-clouds (IVCs) which are often observed in the halos of disk galaxies. The high-velocity-clouds that are also observed but at higher altitudes (of up to 12 kpc) require another mechanism to explain their production. We argue that they could be formed either by the capture of gas from the intergalactic medium and/or by the action of magnetic fields that are carried to the halo with the gas in the fountains. Due to angular momentum losses to the halo, we find that the fountain material falls back to smaller radii and is not largely spread over the galactic disk. Instead, the SNe ejecta fall nearby the region where the fountain was produced, a result which is consistent with recent chemical models of the galaxy. The fall back material leads to the formation of new generations of molecular clouds and to supersonic turbulence feedback in the disk.     

Clump-driven galaxy evolution and quasar activity

The collapse of a protogalaxy composed of dark matter and primordial gas has been investigated by numerical simulations and analytical multi-zone modelling in an attempt to examine the early evolution of disk galaxies. The importance of ample interstellar matter existing in young galactic disks has been highlighted. Confrontation of the theoretical results with the available observational data has led to a new picture for disk galaxy evolution, in which the bulge is the secondary object formed from disk matter. Occurrence of quasar activity is also discussed in relation to the dynamical evolution of the host galaxy.     

Magnetically driven jets from accretion disks: The effect of magneto-rotational instability

We present the results of 2.5-dimensional MHD simulations for jet formation from accretion disks in a situation such that the magneto-rotational instability is occurring in the disk. When there is no initial perturbation in the disk, the surface layer falls faster than the equatorial part to make a surface avalanche and associated jets. However, if we input an initially large perturbation (> 10% of sound speed) in the disk, the dense region of the disk falls on an orbital time scale to make a more violent accretion and jet. In this case, the accretion rate and mass loss rates are an order of magnitude larger than those in the case without initial perturbation. The speed of the jet is of order the Keplerian velocity of the disk regardless of the instability.     

Spacecraft formation flying control around L2 sun-earth libration point using on–off SDRE approach

In this paper, on–off SDRE control approach is presented for spacecraft formation flying control around sun-earth L2 libration point. Orbits around libration points are significant targets for many space missions mainly because of efficient fuel consumption. Furthermore, less propellant usage can be achieved by considering optimal control approaches in spacecraft formation flying control design. Among various nonlinear and optimal control methods, SDRE has shown to be a popular controller in various missions due to the privileges including efficiency, accuracy and robustness. The spacecraft are assumed to have on–off thrusters as actuators. It requires them to be fed with a sequence of on–off pulses which is regarded as a challenge for spacecraft designers. Hence, the main contribution of this paper is designing an on–off SDRE approach for the formation flight around sun-earth L2 point with uncertainty with energy and accuracy considerations. Including on–off input as a constraint is not feasible for SDRE implementation because it makes the system non-affine. An alternative is utilizing an integral action technique and an auxiliary control to make the system affine which leads to on–off SDRE approach. It has also been shown that the proposed method is robust against parametric uncertainties of the states. Present study aims to design an energy-beneficial, simple and attractive controller for a complex nonlinear system with on–off inputs and uncertainty in CRTBP. Simulation results show that the on–off SDRE control could provide the formation flight around L2 point with high accuracy using less energy consumption.     

Time dependent accretion onto a black hole in a disk plus corona system

Hard X-ray lightcurves exhibit delays of 1 s with respect to the soft X-ray lightcurves when the microquasar GRS 1915+105 is in the state of frequent, regular outbursts. Such outbursts are supposed to be driven by the radiation pressure instability of the inner disk parts. The hard X-ray delays are then caused by the time needed for the adjustment of the corona to changing conditions in the underlying disk. We support this claim by the computation of the time evolution of the disk, including a non-stationary evaporation of the disk and mass exchange with the corona.     

Cold gas contents of galaxies and its relation to fueling: Gas dynamics and related phenomena in bar-driven evolution of disk galaxies

The distribution and the dynamics of the cold gas in the circumnuclear regions (r ≲ [1 − a few] kpc) of disk galaxies have been observed at high resolutions of a few arcsec in λ2.6 mm CO (J=1→0) line emission. Distinct CO features are commonly found in the observed disk galaxies and they are different from galaxy to galaxy. They are explained by means of orbit-resonance theories and dynamical evolution. The evolutionary sequence in time is constructed based on the studies of dynamics in theoretical works and numerical simulations in the literature. Not only the behavior of the cold gas but also the starburst, outbreaks of the active galactic nucleus, and the evolution of global structures in disk galaxies are discussed in the bar-driven evolution scenario.     

定常对称背涡流态下细长体绕流结构的演变

应用数值计算的方法模拟了细长体截面绕流结构的演变过程.指出随着细长体背涡的发展,导致截面流场的拓扑结构发生变化,会出现一种临界流动状态.并用微分方程的定性理论分析了此时流场中出现的一种高阶奇点.这种高阶奇点的指数为-3/2,它是结构不稳定的,稍有扰动就会产生分叉,使流场的拓扑结构发生变化.得出了定常对称背涡流态下细长体的空间绕流结构图.     

The search for extra-solar planetary systems.

I review the observational evidence for planetary systems around nearby stars and, using our own solar system as a guide, assess the stringent requirements that new searches need to meet in order to unambiguously establish the presence of another planetary system. Basically, these requirements are: 1 milliarcsecond or better positional accuracy for astrometric techniques, 9 orders of magnitude or better star to planet luminosity ratio discrimination at 0.5 to 1" separation in the optical for direct imaging techniques, 10 meters sec-1 or better radial velocity accuracy for reflex motion techniques and +/-1% or better brightness fluctuation accuracy for planet/star occultation measurements. The astrometric accuracy is in reach of HST, direct imaging will require much larger telescopes and/or a 50 times smoother mirror than HST while the reflex motion and occultation techniques best performed on the ground are just becoming viable and promise exciting new discoveries. On the other band, new indirect evidence on the existence of other planetary systems also comes from the observation of large dusty disks around nearby main sequence stars not too dissimilar from our sun. In one particular case, that of Beta Pictoris, a flattened disk seen nearly edge-on has been imaged in the optical and near IR down to almost 70 AU of the star. It probably represents a young planetary system in its clearing out phase as planetesimals collide, erode and are swept out of the inner system by radiation pressure. The hypothesized Kuiper belt around our solar system may be the analogous structure in a later evolutionary stage. Features of this type have been detected in the far IR and sub-millimeter wavelength regions around 50-100 nearby main sequence and pre-main sequence stars. I discuss a battery of new accurate observations planned in the near future of these objects some of which may actually harbour planets or planetesimals that will certainly dramatically improve our knowledge of planetary system formation processes and our peculiar position in this scheme.     

Contraction of magnetized gas to the central gravitator and the production of magnetically accelerated jets

A magnetodynamic model proposed and worked on by Uchida and collaborators for the star formation jets (bipolar flows), by taking a genetic point of view into account, is reviewed. A large scale magnetic field, which is week in the primordial gas but intensified in the gravitational contraction of the gas, assists the continuous accretion of the gas to the central gravitator by extracting angular momentum from the rotating disk, and this process creates a large amplitude torsional Alfven wave that swirls up the gas into the direction of the axis. This torsional Alfven wave, as it propagates, pinches the large scale field into a slender strong field structure which we identify with the collimated jet. The rationale for extending this mechanism to the AGN jet cases is given, and some results of application to the AGN jet case are presented, with interpretation of some characteristic features like the wiggling of the jets, extended radio lobes and the hotspots at the end of the jets.     

Composite profile of the Fe Kα spectral line emitted from a binary system of supermassive black holes

We used a model of a relativistic accretion disk around a supermassive black hole (SMBH), based on ray-tracing method in the Kerr metric, to study the variations of the composite Fe Kα line emitted from two accretion disks around SMBHs in a binary system. We assumed that the orbit of such a binary is approximately Keplerian, and simulated the composite line shapes for different orbital elements, accretion disk parameters and mass ratios of the components. The obtained results show that, if observed in the spectra of some SMBH binaries during their different orbital phases, such composite Fe Kα line profiles could be used to constrain the orbits and several properties of such SMBH binaries.     

The evolution of globular clusters

Taking a wide range of the initial conditions, including spatial density distribution and mass function etc, the dynamical evolution of globular clusters in the Milky Way is investigated in detail by means of Monte Carlo simulations. Four dynamic mechanisms are considered: stellar evaporation, stellar evolution, tidal shocks due to both the disk and bulge, and dynamical friction. It is found that stellar evaporation dominates the evolution of low-mass clusters and all four are important for massive ones. For both the power-law and lognormal initial clusters mass functions, we can find the best-fit models which can match the present-day observations with their main features of the mass function almost unchanged after evolution of several Gyr. This implies that it is not possible to determine the initial mass function only based on the observed ones today. Moreover, the dispersion of the modelled mass functions mainly depends on the potential wells of host galaxies with the almost constant peaks, which is consistent with current observations.     

A quantitative model for comet nucleus topography

Topography of the surface of a comet nucleus is likely rough at all scales smaller than the mean effective radius. We present a flexible and easily scalable model for quantitative calculations simulating the effects of comet nucleus topography on gas release and dust mantle evolution. The topographic features we describe must be large enough (typically> 10 m) so that they will not erode in one orbit of the nucleus around the Sun. The maximum effective size of a hill is about 1/√2 times the effective radius of the nucleus. If it is larger, then an ellipsoidal shape of the nucleus is more appropriate. The procedure described here also permits for inhomogeneous composition of the topographic features, leading to locally different rates of gas production (e.g., jet-like features and filaments) or different thicknesses of the dust mantle. It also can give rise to different temperature patches, locally varying albedos and emissivities, and may explain the formation of permanent dust mantles.     

Black hole spectral states and physical connections

The dramatic changes seen in the X-ray spectral and timing properties of accreting black hole candidates (BHCs) provide important clues about the accretion and jet formation processes that occur in these systems. Dividing the different source behaviors into spectral states provides a framework for studying BHCs. To date, there have been three main classification schemes with Luminosity-based, Component-based, or Transition-based criteria. The canonical, Luminosity-based criteria and physical models that are based on this concept do not provide clear explanations for several phenomena, including hysteresis of spectral states and the presence of jets. I discuss the re-definitions of states, focusing on an application of the Component-based states to more than 400 RXTE observations of the recurrent BHC 4U 1630–47. We compare the X-ray properties for the recent 2002–2004 outburst to those of an earlier (1998) outburst, during which radio jets were observed. The results suggest a connection between hysteresis of states and major jet ejections, and it is possible that both of these are related to the evolution of the inner radius of the optically thick accretion disk.     

Numerical Simulation for the 16 August 1999 EUV Brightenings

The 16 August 1999 EUV brightenings are numerically simulated by a third-order upwind compact scheme, basing on the TRACE observation. The present simulation can give a possible explanation to its formation and evolution. The numerical results show that the initial reconnection jets at around X-point are responsible for the occurrence of EUV brightening. The strong and superposed ejections caused by the first and second coalescence of magnetic islands are possibly related to the lifted material which initially appeared as absorption features and Later EUV-emitting structures respectively. The bi-directional reconnection jets may correspond to the lifted material that either continued to move upward along the apparently open field lines or fell down to the surface.      

The relevance of the Eddington limit to thick accretion disks

Thick accretion disks with narrow funnels around massive black holes are considered promising models for active galactic nuclei. These models assume a supercritical accretion rate and emit collimated beams with super-Eddington luminosities. We have made approximate calculations of the interaction between the emerging radiation and the walls for an optically thin funnel. The results are sensitive to the sound velocity and to the viscosity parameter α. They suggest that a significant particle luminosity can accompany the radiation in the super-Eddington case. By applying an Eddington type limit based on mechanical equilibrium to a suitably chosen interior surface in the disk, we find that energy transport, if radiative, can strongly limit the efficiency of these models.     

Reconstruction of the solar UV irradiance back to 1974

The variability of the solar UV irradiance has strong effects on the terrestrial atmosphere. In order to study the solar influence for times when no UV observations are available, it is necessary to reconstruct the variation of the UV irradiance with time on the basis of proxies. We present reconstructions of the solar UV irradiance based on the analysis of space-based and ground-based magnetograms of the solar disk going back to 1974. With COde for Solar Irradiance (COSI) we calculate solar intensity spectra for the quiet Sun and different active regions and combine them according to their fractional area on the solar disk, whereby their time-dependent contributions over the solar cycle lead to a variability in radiation. COSI calculates the continuum and line formation under conditions which are out of local thermodynamic equilibrium (non-LTE). The applied temperature and density structures include the chromosphere and transition region, which is particularly important for the UV. The reconstructions are compared with observations.     

METHIS: Mercury thermal infrared spectrometer

The BepiColombo mission to Mercury is devoted to the thorough exploration of Mercury and its environment, with the aim to understand the processes of planetary formation and evolution in the hottest part of the protoplanetary nebula. This mission represents an unique opportunity for the European community to extend the understanding of the Solar Nebula evolution from its outer edge – ideally represented by comets – to its inner and warmer edge. Obviously this exploration asks for a detailed knowledge of the main constituents of the matter present in the different Solar System areas. Spectroscopy is a powerful tool to acquire this knowledge. We have participated with a large consortium of European researchers to the development of the Rosetta imaging spectrometer. We propose here to use our experience to develop a newly designed spectrometer to investigate the mineralogical composition of the Mercurial surface. Given the particular thermodynamical situation of the Mercurial surface, we have developed a concept that combines a medium IR low spectral resolution imager with a moderate spectral resolution NIR point spectrometer. The main goal of METHIS is to provide the mineralogical characterisation of the surface with sufficient spectral resolution in a scientifically diagnostic spectral range.