The 1.5 Seyfert galaxy NGC 3227 has been observed by several X-ray missions. We carried out combined analysis of the data obtained by more recent major observations of this source - two observations performed by XMM-Newton in 2000 and 2006 and six observations performed by Suzaku in 2008. A unified model was constructed which is consistent with all eight of the observations by the two satellites with large intensity and spectral changes. The model consists of a hard power law with the spectral index of 1.35-1.65 which is interpreted as the Comptonized emission from the corona above an accretion disk. In the high flux states an additional soft excess component dominates, which is consistent with a model with either a steeper power law with an index of 3.5-3.6 or the warm Comptonization component. These emissions from the central engine are absorbed by a neutral partial covering material and warm absorbers. A reflection component and several emission lines are also present. We examined the relationship between the intrinsic luminosity and the absorbers' physical parameters such as the column density, which suggests that the source expanded significantly during the bright states where the soft excess is greatly enhanced.

If a magnetically dominated magnetosphere is to extract a black hole's rotational energy and transmit it to distant regions, then the inner light surface of that magnetosphere must lie within the ergoregion. That inner light surface condition limits the angular velocity of magnetic field lines. We take the distribution of magnetic field line angular velocity on the horizon to be a useful proxy for inner light surface location and study how different distributions affect the structure of energy-extracting magnetospheres. Within magnetospheres that exhibit differential field line bending towards both the azimuthal axis and the equatorial plane, we find that the total Poynting flux energy directed outward along the azimuthal axis can vary by over a factor of 100 for a single value of black hole spin.

We study the effects of frame-dragging on the structure of force-free magnetospheres around rotating black holes. For slowly rotating black holes, analytic explorations often focus on transforming an exact force-free solution applicable to a static black hole into a solution approximately applicable to a slowly-rotating black hole via perturbations in black hole spin. We show that the single perturbed monopolar solution most commonly arrived at using such techniques is in fact a separatrix between two broad classes of solutions: those with poloidal magnetic fields that bend upwards towards the azimuthal axis and those that bend downwards towards the equatorial plane. This is because frame-dragging affects the toroidal magnetic field distributions, changing the force balance between poloidal magnetic field lines.

An energy extracting black hole magnetosphere can be defined by the location of its inner Alfv\'{e}n surface, which determines the rate of black hole energy extraction along a given magnetic field line. We study how the location of the inner Alfv\'{e}n surface can modify the structure of energy extracting black hole magnetospheres in the force-free limit. Hundreds of magnetospheres are numerically computed via a general relativistic extension of the Newtonian magnetofrictional method for a full range of black hole spins and flow parameters. We find that jet-like structures naturally form very close to the horizon for Alfv\'{e}n surfaces near the boundary of the ergosphere and that energy is extracted towards the equatorial plane for Alfv\'{e}n surfaces close to the horizon. This suggests that two broad classes of energy extracting black hole magnetospheres might exist; those that transmit extracted energy directly to distant observers, and those that transmit extracted energy to nearby accreting matter. Applied to transient high energy phenomena, we find that they might also differ in timescale by a factor of 20 or more.

Using the latest 70 month Swift-BAT catalog we examined hard X-ray selected Seyfert I galaxies which are relatively little known and little studied, and yet potentially promising to test the ionized relativistic reflection model. From this list we chose 13 sources which have been observed by XMM-Newton for less than 20 ks, in order to explore the broad band soft to hard X-ray properties with the analysis of combined XMM-Newton and Swift data. Out of these we found seven sources which exhibit potentially promising features of the relativistic disc reflection, such as a strong soft excess, a large Compton hump and/or a broadened Fe line. Longer observations of four of these sources with the currently operating satellite missions, such as Suzaku, XMM-Newton and NuStar and two others by such future missions as ASTRO-H, will be invaluable, in order to better understand the relativistic disc reflection closest to the central black hole and constrain such important effects of strong gravity as the black hole spin.

We analysed new simultaneous \emph{Suzaku} and \emph{Swift} data of the Broad Line Seyfert 1 (BLS1) galaxy CBS 126. A clear Fe emission line and a strong soft excess are present in the source spectra. We fit the spectra with a relativistic reflection model and find the model tends to fit the data with a high iron abundance possibly due to the large soft excess present. By checking the difference and the RMS spectra, we find there is likely an absorption edge at $\sim$ 0.89 keV, which might be caused by oxygen or neon. We also produced an analysis of the time-resolved spectra in order to examine the existence of the edge. Although high iron abundance is required more in the time-resolved spectra, a model of solar iron abundance together with an absorption edge gives a more physical explanation. Variation of the ionisation parameter is an alternative, plausible explanation for the excess seen in the difference spectra. It is difficult to know if there are warm absorbers in this source from the current data. To further investigate the presence of possible warm absorbers, higher signal-to-noise low-energy data are needed. The work presented here tentatively suggests that the spectra of a BLS1 can be explained by a relativistic reflection model similar to that often seen in their narrow-line cousins.

A number of scenarios have been proposed for the origin of the supermassive black holes (SMBHs) that are found in the centres of most galaxies. Many such scenarios predict a high-redshift population of massive black holes (MBHs), with masses in the range 100 to 100000 times that of the Sun. When the Laser Interferometer Space Antenna (LISA) is finally operational, it is likely that it will detect on the order of 100 of these MBH binaries as they merge. The differences between proposed population models produce appreciable effects in the portion of the population which is detectable by LISA, so it is likely that the LISA observations will allow us to place constraints on them. However, gravitational wave detectors such as LISA will not be able to detect all such mergers nor assign precise black hole parameters to the merger, due to weak gravitational wave signal strengths. This paper explores LISA's ability to distinguish between several MBH population models. In this way, we go beyond predicting a LISA observed population and consider the extent to which LISA observations could inform astrophysical modellers. The errors in LISA parameter estimation are applied with a direct method which generates random sample parameters for each source in a population realisation. We consider how the distinguishability varies depending on the choice of source parameters (1 or 2 parameters chosen from masses, redshift or spins) used to characterise the model distributions, with confidence levels determined by 1 and 2-dimensional tests based on the Kolmogorov-Smirnov test.

Many scenarios have been proposed for the origin of the supermassive black holes (SMBHs) that are found in the centres of most galaxies. Many of these formation scenarios predict a high-redshift population of intermediate-mass black holes (IMBHs), with masses in the range 100 to 100000 times that of the Sun. A powerful way to observe these IMBHs is via gravitational waves the black holes emit as they merge. The statistics of the observed black hole population should, in principle, allow us to discriminate between competing astrophysical scenarios for the origin and formation of SMBHs. However, gravitational wave detectors such as LISA will not be able to detect all such mergers nor assign precise black hole parameters to the merger, due to weak gravitational wave signal strengths. In order to use LISA observations to infer the statistics of the underlying population, these errors must be taken into account. We describe here a method for folding the LISA gravitational wave parameter error estimates into an `error kernel' designed for use at the population model level. The effects of this error function are demonstrated by applying it to several recent models of black hole mergers, and some tentative conclusions are made about LISA's ability to test scenarios of the origin and formation of supermassive black holes.

We propose a model in which intermediate-mass black holes (IMBHs) with mass of ~10000 Msun are formed in early dark matter halos. We carry out detailed stellar evolution calculations for accreting primordial stars including annihilation energy of dark matter particles. We follow the stellar core evolution consistently up to gravitational collapse. We show that very massive stars, as massive as 10000 Msun, can be formed in an early dark matter halo. Such stars are extremely bright with Log L/Lsun > 8.2. They gravitationally collapse to form IMBHs. These black holes could have seeded the formation of early super-massive blackholes.

We calculate the evolution of zero-metallicity Population III (Pop III) stars whose mass grows from the initial mass of $\sim 1M_{\odot}$ by accreting the surrounding gases. Our calculations cover a whole evolutionary stages from the pre-main sequence, via various nuclear burning stages, through the final core collapse or pair-creation instability phases. We adopt the following stellar mass-dependent accretion rates which are derived from cosmological simulations of early structure formation based on the low mass dark matter halos at redshifts $z \sim 20$: (1) the accretion rates for the first generation (Pop III.1) stars and (2) the rates for zero-metallicity but the second generation (Pop III.2) stars which are affected by radiation from the Pop III.1 stars. For comparison, we also study the evolution with the mass-dependent accretion rates which are affected by radiatibe feedback. We show that the final mass of Pop III.1 stars can be as large as $\sim 1000M_{\odot}$, beyond the mass range ($140 - 300M_{\odot}$) for the pair-instability supernovae. Such massive stars undergo core-collapse to form intermediate-mass black holes, which may be the seeds for merger trees to supermassive black holes. On the other hand, Pop III.2 stars become less massive ($\lsim 40 - 60M_{\odot}$), being in the mass range of ordinary iron core-collapse stars. Such stars explode and eject heavy elements to contribute to chemical enrichment of the early universe as observed in the abundance patterns of extremely metal-poor stars in the Galactic halo.

We study magnetohydrodynamic (MHD) standing shocks in inflowing plasmas in a black hole magnetosphere. Fast and intermediate shock formation is explored in Schwarzschild and Kerr geometry to illustrate general relativistic effects. We find that non-equatorial standing MHD shocks are physically possible, creating a very hot plasma region close to the event horizon. Shocked downstream plasmas can be heated or magnetized depending on the values of various magnetic field-aligned parameters. Then we may expect high-energy thermal/nonthermal emissions from the shocked region. We present the properties of non-equatorial MHD shocks and discuss the shocked plasma region in the black hole magnetosphere. We also investigate the effects of the poloidal magnetic field and the black hole spin on the properties of shocks, and show that both effects can modify the distribution of the shock front and shock strength. We find for strong MHD shock formation that fast rotating magnetic fields are necessary. The physics of non-equatorial MHD shocks in the black hole magnetosphere could be very important when we are to construct the central engine model of various astrophysical phenomena.

We calculate evolution, collapse, explosion, and nucleosynthesis of Population III very-massive stars with 500$M_{\odot}$ and 1000$M_{\odot}$. Presupernova evolution is calculated in spherical symmetry. Collapse and explosion are calculated by a two-dimensional code, based on the bipolar jet models. We compare the results of nucleosynthesis with the abundance patterns of intracluster matter, hot gases in M82, and extremely metal-poor stars in the Galactic halo. It was found that both 500$M_{\odot}$ and 1000$M_{\odot}$ models enter the region of pair-instability but continue to undergo core collapse. In the presupernova stage, silicon burning regions occupy a large fraction, more than 20% of the total mass. For moderately aspherical explosions, the patterns of nucleosynthesis match the observational data of both intracluster medium and M82. Our results suggest that explosions of Population III core-collapse very-massive stars contribute significantly to the chemical evolution of gases in clusters of galaxies. For Galactic halo stars, our [O/Fe] ratios are smaller than the observational abundances. However, our proposed scenario is naturally consistent with this outcome. The final black hole masses are $\sim 230M_{\odot}$ and $\sim 500M_{\odot}$ for the $500M_{\odot}$ and 1000$M_{\odot}$ models, respectively. This result may support the view that Population III very massive stars are responsible for the origin of intermediate mass black holes which were recently reported to be discovered.

Recent years have seen some significant progress in theoretical studies of physics of dense matter. Combined with the observational data now available from the successful launch of Chandra and XMM/Newton X-ray space missions as well as various lower-energy band observations, these developments now offer the hope for distinguishing various competing neutron star thermal evolution models. For instance, the latest theoretical and observational developments may already exclude both nucleon and kaon direct Urca cooling. In this way we can now have a realistic hope for determining various important properties, such as the composition, superfluidity, the equation of state and stellar radius. These developments should help us obtain deeper insight into the properties of dense matter.

Fast and slow magnetosonic shock formation is presented for stationary and axisymmetric magnetohydrodynamical (MHD) accretion flows onto a black hole. The shocked black hole accretion solution must pass through magnetosonic points at some locations outside and inside the shock location. We analyze critical conditions at the magnetosonic points and the shock conditions. Then, we show the restrictions on the flow parameters for strong shocks. We also show that a very hot shocked plasma is obtained for a very high-energy inflow with small number density. Such a MHD shock can appear very close to the event horizon, and can be expected as a source of high-energy emissions. Examples of shocked MHD accretion flows are presented in the Schwarzschild case.

We present 6.4 keV iron K-alpha fluorescent line profiles predicted for a relativistic black hole accretion disk in the presence of a spiral motion in Kerr geometry, the work extended from an earlier literature motivated by recent magnetohydrodynamic (MHD) simulations. The velocity field of the spiral motion, superposed on the background Keplerian flow, results in a complicated redshift distribution in the accretion disk. An X-ray source attributed to a localized flaring region on the black hole symmetry axis illuminates the iron in the disk. The emissivity form becomes very steep because of the light bending effect from the primary X-ray source to the disk. The predicted line profile is calculated for various spiral waves, and we found, regardless of the source height, that: (i) a multiple-peak along with a classical double-peak structure generally appears, (ii) such a multiple-peak can be categorized into two types, sharp sub-peaks and periodic spiky peaks, (iii) a tightly-packed spiral wave tends to produce more spiky multiple peaks, whereas (iv) a spiral wave with a larger amplitude seems to generate more sharp sub-peaks, (v) the effect seems to be less significant when the spiral wave is centrally concentrated, (vi) the line shape may show a drastic change (forming a double-peak, triple-peak or multiple-peak feature) as the spiral wave rotates with the disk. Our results emphasize that around a rapidly-rotating black hole an extremely redshifted iron line profile with a noticeable spike-like feature can be realized in the presence of the spiral wave. Future X-ray observations, from {\it Astro-E2} for example, will have sufficient spectral resolution for testing our spiral wave model which exhibits unique spike-like features.

We explore isothermal shock formation in non-equatorial, adiabatic accretion flows onto a rotating black hole, with possible application to some active galactic nuclei (AGNs). The isothermal shock jump conditions as well as the regularity condition, previously developed for one-dimensional (1D) flows in the equatorial plane, are extended to two-dimensional (2D), non-equatorial flows, to explore possible geometrical effects. The basic hydrodynamic equations with these conditions are self-consistently solved in the context of general relativity to explore the formation of stable isothermal shocks. We find that strong shocks are formed in various locations above the equatorial plane, especially around a rapidly-rotating black hole with the prograde flows (rather than a Schwarzschild black hole). The retrograde flows are generally found to develop weaker shocks. The energy dissipation across the shock in the hot non-equatorial flows above the cooler accretion disk may offer an attractive illuminating source for the reprocessed features, such as the iron fluorescence lines, which are often observed in some AGNs.

With the successful launch of Chandra and XMM/Newton X-ray space missions combined with the lower-energy band observations, time has arrived when careful comparison of thermal evolution theories of isolated neutron stars with observations will offer a better hope for distinguishing among various competing neutron star cooling theories. For instance, the latest theoretical and observational developments may already exclude both nucleon and kaon direct Urca cooling. In this way we can now have a realistic hope for determining various important properties, such as the composition, superfluidity, the equation of state and stellar radius. These developments should help us obtain deeper insight into the properties of dense matter.

With the successful launch of Chandra and XMM/Newton X-ray space missions combined with the lower-energy band observations, we are in the position where careful comparison of neutron star cooling theories with observations will make it possible to distinguish among various competing theories. For instance, the latest theoretical and observational developments already exclude both nucleon and kaon direct URCA cooling. In this way we can now have realistic hope for determining various important properties, such as the composition, degree of superfluidity, the equation of state and steller radius. These developments should help us obtain better insight into the properties of dense matter.

We extend the work by Appl and Camenzind (1988) for special relativistic magnetohydrodynamic (MHD) jets, to fully general relativistic studies of the standing shock formation for accreting MHD plasma in a rotating, stationary and axisymmetric black hole magnetosphere. All the postshock physical quantities are expressed in terms of the relativistic compression ratio, which can be obtained in terms of preshock quantities. Then, the downstream state of a shocked plasma is determined by the upstream state of the accreting plasma. In this paper sample solutions are presented for slow magnetosonic shocks for accreting flows in the equatorial plane. We find that some properties of the slow magnetosonic shock for the rotating magnetosphere can behave like a fast magnetosonic shock. In fact, it is confirmed that in the limit of weak gravity for the upstream non-rotating accretion plasma where the magnetic field lines are leading and rotating, our results are very similar to the fast magnetosonic shock solution by Appl and Camenzind (1988). However, we find that the situation becomes far more complicated due to the effects of strong gravity and rotation, such as the frame dragging-effects. We show the tendency that the large spin of the black hole makes the slow magnetosonic shock strong for the accretion solutions with the same energy-flux.

Very recently the Chandra First Light Observation discovered a point-like source in Cassiopeia A (Cas A) supernova remnant. This detection was subsequently confirmed by the analyses of the archival data from both ROSAT and Einstein observations. Here we compare the results from these observations with the scenarios involving both black holes (BH) and neutron stars (NS). If this point source is a BH we offer as a promising model a disk-corona type model with low accretion rate where a soft photon source at ~ 0.1 keV is Comptonized by higher energy electrons in the corona. If it is a NS the dominant radiation observed by Chandra most likely originates from smaller, hotter regions of the stellar surface, but we argue that it is still worthwhile to compare the cooler component from the rest of the surface with cooling theories. We emphasize that the detection of this point source itself should potentially provide enormous impacts on the theories of supernova explosion, progenitor scenario, compact remnant formation, accretion to compact objects, and neutron star thermal evolution.

The thermal evolution of a neutron star is studied by including the energy loss due to axion emission. Two axion models and three types of neutron-star matter equation of state are used with the effects of nucleon superfluidity properly taken into account. In comparison with the observational data of PSR0656+14 from ROSAT, the upper limits on the axion mass are found to be $m_a < 0.06 - 0.3$ eV and 0.08$ -$ 0.8 eV for the KSVZ and DFSZ axion models, respectively, with the soft equation of state giving the most stringent limits.

In this paper we show that perturbations of the accretion flow within the central engines of active galctic nuclei (AGN) are likely to form shock waves in the accreting plasma. Such shocks, which may be either collisional or collisionless, can contribute to the observed high energy temporal and spectral variability. Our Rationale is the following: Observations show that the continuum emission probably originates in an optically thin, hot plasma in the AGN central engine. The flux and spectrum from this hot plasma varies significantly over light crossing timescales. several authors have suggested that macroscopic perturbations contained within this plasma are the sources of this variability. In order to produce the observed emission the perturbations must be radiatively coupled with the optically thin hot matter and must also move with high velocities. We suggest that shocks, which can be very effective in randomizing the bulk motion of the perturbations, are responsible for this coupling. Shocks should form in the central engine, because the temperatures and magnetic fields are probably reduced below their virial values by radiative dissipation. Perturbations moving at Keplerian speeds, or strong nonlinear excitations, result in supersonic and super-Alfvenic velocities leading to shock waves whithin the hot plasma We Show that even a perturbation significantly smaller than the emitting region can form a shock which significantly modifies the continuum emission in an AGN, and that the spectral and temporal variability from such a shock generally resemble those of radio quiet AGN. As an example, the shock inducing perturbation in our model is a small main sequence star, the capturing and eventual accretion of which are known to be a plausible process. We argue that