We systematically analyzed the high-quality Suzaku data of 88 Seyfert galaxies. We obtained a clear relation between the absorption column density and the equivalent width of the 6.4 keV line above 10$^{23}$ cm$^{-2}$, suggesting a wide-ranging column density of $10^{23-24.5}$ cm$^{-2}$ with a similar solid and a Fe abundance of 0.7--1.3 solar for Seyfert 2 galaxies. The EW of the 6.4 keV line for Seyfert 1 galaxies are typically 40--120 eV, suggesting the existence of Compton-thick matter like the torus with a column density of $>10^{23}$ cm$^{-2}$ and a solid angle of $(0.15-0.4)*4pi$, and no difference of neutral matter is visible between Seyfert 1 and 2 galaxies. An absorber with a lower column density of $10^{21-23}$ cm$^{-2}$ for Compton-thin Seyfert 2 galaxies is suggested to be not a torus but an interstellar medium. These constraints can be understood by the fact that the 6.4 keV line intensity ratio against the 10--50 keV flux is almost identical within a range of 2--3 in many Seyfert galaxies. Interestingly, objects exist with a low EW, 10--30 eV, of the 6.4 keV line, suggesting that those torus subtends only a small solid angle of $<0.2*4pi$. Ionized Fe-K$\alpha$ emission or absorption lines are detected from several percents of AGNs. Considering the ionization state and equivalent width, emitters and absorbers of ionized Fe-K lines can be explained by the same origin, and highly ionized matter is located at the broad line region. The rapid increase in EW of the ionized Fe-K emission lines at $N_{H}>10^{23}$ cm$^{-2}$ is found, like that of the cold material. It is found that these features seem to change for brighter objects with more than several $10^{44}$ erg/s such that the Fe-K line features become weak. We discuss this feature, together with the torus structure.

Suzaku Hard X-ray Detector (HXD) achieved the lowest background level than any other previously or currently operational missions sensitive in the energy range of 10--600 keV, by utilizing PIN photodiodes and GSO scintillators mounted in the BGO active shields to reject particle background and Compton-scattered events as much as possible. Because it does not have imaging capability nor rocking mode for the background monitor, the sensitivity is limited by the reproducibility of the non X-ray background (NXB) model. We modeled the HXD NXB, which varies with time as well as other satellites with a low-earth orbit, by utilizing several parameters, including particle monitor counts and satellite orbital/attitude information. The model background is supplied as an event file in which the background events are generated by random numbers, and can be analyzed in the same way as the real data. The reproducibility of the NXB model depends on the event selection criteria (such as cut-off rigidity and energy band) and the integration time, and the 1sigma systematic error is estimated to be less than 3% (PIN 15--40 keV) and 1% (GSO 50--100 keV) for more than 10 ksec exposure.

A perturbative analysis is used to investigate the effect of rotation on the instability of a steady accretion shock (SASI) in a simple toy-model, in view of better understanding supernova explosions in which the collapsing core contains angular momentum. A cylindrical geometry is chosen for the sake of simplicity. Even when the centrifugal force is very small, rotation can have a strong effect on the non axisymmetric modes of SASI by increasing the growth rate of the spiral modes rotating in the same direction as the steady flow. Counter-rotating spiral modes are significantly damped, while axisymmetric modes are hardly affected by rotation. The growth rates of spiral modes have a nearly linear dependence on the specific angular momentum of the flow. The fundamental one-armed spiral mode (m=1) is favoured for small rotation rates, whereas stronger rotation rates favour the mode m=2. A WKB analysis of higher harmonics indicates that the efficiency of the advective-acoustic cycles associated to spiral modes is strongly affected by rotation in the same manner as low frequency modes, whereas the purely acoustic cycles are stable. These results suggest that the linear phase of SASI in rotating core-collapse supernovae naturally selects a spiral mode rotating in the same direction of the flow, as observed in the 3D numerical simulations of Blondin & Mezzacappa (2007). This emphasizes the need for a 3D approach of rotating core-collapse, before conclusions on the explosion mechanisms and pulsar kicks can be drawn.

A Suzaku observation of the nucleus of the radio-loud AGN Centaurus A in 2005 has yielded a broadband spectrum spanning 0.3 to 250 keV. The net exposure times after screening were: 70 ks per X-ray Imaging Spectrometer (XIS) camera, 60.8 ks for the Hard X-ray Detector (HXD) PIN, and 17.1 ks for the HXD-GSO. The hard X-rays are fit by two power-laws of the same slope, absorbed by columns of 1.5 and 7 * 10^{23} cm^{-2} respectively. The spectrum is consistent with previous suggestions that the power-law components are X-ray emission from the sub-pc VLBI jet and from Bondi accretion at the core, but it is also consistent with a partial covering interpretation. The soft band is dominated by thermal emission from the diffuse plasma and is fit well by a two-temperature VAPEC model, plus a third power-law component to account for scattered nuclear emission, jet emission, and emission from X-ray Binaries and other point sources. Narrow fluorescent emission lines from Fe, Si, S, Ar, Ca and Ni are detected. The Fe K alpha line width yields a 200 light-day lower limit on the distance from the black hole to the line-emitting gas. Fe, Ca, and S K-shell absorption edges are detected. Elemental abundances are constrained via absorption edge depths and strengths of the fluorescent and diffuse plasma emission lines. The high metallicity ([Fe/H]=+0.1) of the circumnuclear material suggests that it could not have originated in the relatively metal-poor outer halo unless enrichment by local star formation has occurred. Relative abundances are consistent with enrichment from Type II and Ia supernovae.

The in-orbit performance and calibration of the Hard X-ray Detector (HXD) on board the X-ray astronomy satellite Suzaku are described. Its basic performances, including a wide energy bandpass of 10-600 keV, energy resolutions of ~4 keV (FWHM) at 40 keV and ~11% at 511 keV, and a high background rejection efficiency, have been confirmed by extensive in-orbit calibrations. The long-term gains of PIN-Si diodes have been stable within 1% for half a year, and those of scintillators have decreased by 5-20%. The residual non-X-ray background of the HXD is the lowest among past non-imaging hard X-ray instruments in energy ranges of 15-70 and 150-500 keV. We provide accurate calibrations of energy responses, angular responses, timing accuracy of the HXD, and relative normalizations to the X-ray CCD cameras using multiple observations of the Crab Nebula.

Bearing in mind the application to the theory of core-collapse supernovae, we performed a global linear analysis on the stability of spherically symmetric accretion flows through a standing shock wave onto a proto neutron star. As unperturbed flows, we adopted the spherically symmetric steady solutions to the Euler equations obtained with realistic equation of state and formulae for neutrino reaction rates taken into account. Then we solved the equations for linear perturbations numerically, and obtained the eigen frequencies and eigen functions. We found (1) the flows are stable for all modes if the neutrino luminosity is lower than $\sim 1\times 10^{52}$ ergs/s for $\dot{M}=1.0M_{\odot}/{\rm s}$. (2) For larger luminosities, the non-radial instabilities are induced, probably via the advection-acoustic cycles. Interestingly, the modes with $\ell=2$ and 3 become unstable at first for relatively low neutrino luminosities, e.g. $\gtrsim 2-3\times 10^{52}$ ergs/s for the same accretion rate, whereas the $\ell=1$ mode is the most unstable for higher luminosities, $\sim 3-7\times 10^{52}$ ergs/s. These are all oscillatory modes. (3) For still larger luminosities, $\sim 7\times 10^{52}$ ergs/s for $\dot{M}=1.0M_{\odot}/{\rm s}$, non-oscillatory modes, both radial and non-radial, become unstable. These non-radial modes were identified as convection. We confirmed the results obtained by numerical simulations that the instabilities induced by the advection-acoustic cycles are more important than the convection for lower neutrino luminosities.

We report on a 100 ks Suzaku observation of the bright, nearby (z=0.008486) Seyfert 1.9 galaxy MCG -5-23-16. The broad-band (0.4-100 keV) X-ray spectrum allows us to determine the nature of the high energy emission with little ambiguity. The X-ray continuum consists of a cutoff power-law of photon index $\Gamma=1.9$, absorbed through Compton-thin matter of column density $N_{\rm H}=1.6\times10^{22}$ cm$^{-2}$. A soft excess is observed below 1 keV and is likely a combination of emission from scattered continuum photons and distant photoionized gas. The iron K line profile is complex, showing narrow neutral iron K$\alpha$ and K$\beta$ emission, as well as a broad line which can be modeled by a moderately inclined accretion disk. The line profile shows either the disk is truncated at a few tens of gravitational radii, or the disk emissivity profile is relatively flat. A strong Compton reflection component is detected above 10 keV, which is best modeled by a combination of reflection off distant matter and the accretion disk. The reflection component does not appear to vary. The overall picture is that this Seyfert 1.9 galaxy is viewed at moderate (50 degrees) inclination through Compton-thin matter at the edge of a Compton-thick torus covering $2\pi$ steradians, consistent with unified models.

We investigated the structure of the spherically symmetric accretion flows through the standing shock wave onto the proto-neutron star in the post-bounce phase of the collapse-driven supernova. We assume that the accretion flow is in a steady state controlled by the neutrino luminosity and mass accretion rate that are kept constant. We obtain solutions of the steady Euler equations for a wide range of neutrino luminosity and mass accretion rate. We employ a realistic EOS and neutrino-heating rates. More importantly, we take into account the effect of convection phenomenologically. For each mass accretion rate, we find the critical neutrino luminosity, above which there exists no steady solution. These critical points are supposed to mark the onset of the shock revival. As the neutrino luminosity increases for a given mass accretion rate, there appears a convectively unstable region at some point before the critical value is reached. We introduce a phenomenological energy flux by convection so that the negative entropy gradient should be canceled out. We find that the convection lowers the critical neutrino luminosity substantially. We also consider the effect of the self-gravity. It is found that the self-gravity is important only when the neutrino luminosity is high. The critical luminosity, however, is little affected if the energy transport by convection is taken into account.

We investigate the dynamics of an injected outflow propagating in a progenitor in the context of the collapsar model for gamma-ray bursts (GRBs) through two dimensional axisymmetric relativistic hydrodynamic simulations. Initially, we locally inject an outflow near the center of a progenitor. We calculate 25 models, in total, by fixing its total input energy to be 10^{51} ergs s^{-1} and radius of the injected outflow to be $7\times 10^7$ cm while varying its bulk Lorentz factor, $\Gamma_{0} = 1.05\sim 5$, and its specific internal energy, $\epsilon_0/c^2 = 0.1\sim 30$. The injected outflow propagates in the progenitor and drives a large-scale outflow or jet. We find a smooth but dramatic transition from a collimated jet to an expanding outflow among calculated models. The maximum Lorentz factor is, on the other hand, sensitive to both of $\Gamma_0$ and $\epsilon_0$; roughly $\Gamma_{\rm max} \sim \Gamma_0 (1+\epsilon_0/c^2)$. Our finding will explain a smooth transition between the GRBs, X-ray rich GRBs (XRRs) and X-ray Flashes (XRFs) by the same model but with different $\epsilon_0$ values.

In order to infer the effects of rotation on the revival of a stalled shock in supernova explosions, we investigated steady accretion flows with a standing shock. We first obtained a series of solutions for equations describing non-rotating spherically symmetric flows and confirmed the results of preceding papers that, for a given mass accretion rate, there is a critical luminosity of irradiating neutrinos, above which there exists no steady solution. Below the critical value, we found two branches of solutions; one is stable and the other is unstable against radial perturbations. With a simple argument based on the Riemann problem, we can identify the critical luminosity as the one, at which the stalled shock revives. We also obtained the condition satisfied by the flow velocity for the critical luminosity, which can be easily applied to the rotational case. If a collapsing star rotates, the accretion flow is non-spherical due to centrifugal forces. Flows are accelerated near the rotation axis whereas they are decelerated near the equatorial plane. As a result, the critical luminosity is lowered, that is, rotation assists the revival of a stalled shock. According to our calculations, the critical luminosity is $\sim25$% lower for the mass accretion rate of 1M$_{\odot}$/sec and the rotation frequency of 0.1 Hz at a radius of 1000 km than that of the spherically symmetric flow with the same mass accretion rate. We found that the condition of the flow velocity at the critical luminosity is first satisfied at the rotation axis. This suggests that the shock revival is triggered on the rotation axis and a jet-like explosion ensues.