We apply XCLUMPY, an X-ray spectral model from a clumpy torus in an active galactic nucleus (AGN), to the broadband X-ray spectra of 10 obscured AGNs observed with both Suzaku and NuSTAR. The infrared spectra of these AGNs were analyzed with the CLUMPY code. Since XCLUMPY adopts the same clump distribution as that in the CLUMPY, we can directly compare the torus parameters obtained from the X-ray spectra and those from the infrared ones. The torus angular widths determined from the infrared spectra ($\sigma_{\mathrm{IR}}$) are systematically larger than those from the X-ray data ($\sigma_{\mathrm{X}}$); the difference ($\sigma_{\mathrm{IR}}-\sigma_{\mathrm{X}}$) correlates with the inclination angle determined from the X-ray spectrum. These results can be explained by the contribution from dusty polar outflows to the observed infrared flux, which becomes more significant at higher inclinations (more edge-on views). The ratio of the hydrogen column density and V-band extinction in the line of sight absorber shows large scatter ($\simeq$1 dex) around the Galactic value, suggesting that a significant fraction of AGNs have dust-rich circumnuclear environments.

We have analyzed the broadband X-ray spectra of active galactic nuclei (AGNs) in two "non-merging" luminous infrared galaxies (LIRGs) UGC 2608 and NGC 5135, utilizing the data of NuSTAR, Suzaku, XMM-Newton, and Chandra. Applying the X-ray clumpy-torus model (XCLUMPY: Tanimoto et al. 2019), we find that both sources have similar spectra characterized by Compton-thick (CT) absorption ($N_{\rm H} \sim$ 5-7 $\times$ $10^{24}$ cm$^{-2}$) and small torus angular width ($\sigma$ $<$ 20$^{\circ}$). The intrinsic 2-10 keV luminosities are $3.9^{+2.2}_{-1.7}$ $\times$ $10^{43}$ erg s$^{-1}$ (UGC 2608) and $2.0^{+3.3}_{-1.0}$ $\times$ $10^{43}$ erg s$^{-1}$ (NGC 5135). The [O IV]-to-nuclear-12 $\mu$m luminosity ratios are larger than those of typical Seyferts, which are consistent with the torus covering factors ($C_{\rm T} \lesssim$ 0.7) estimated from the torus angular widths and column densities by X-ray spectroscopy. The torus covering factors and Eddington ratios ($\lambda_{\rm Edd} \sim$ 0.1) follow the relation found by Ricci et al. (2017c) for local AGNs, implying that their tori become geometrically thin due to significant radiation pressure of the AGN that blows out some part of the tori. These results indicate that the CT AGNs in these "non-merger" LIRGs are just a normal AGN population seen edge-on through a large line-of-sight column density. They are in contrast to the buried CT AGNs in late-stage mergers that have large torus covering factors even at large Eddington ratios.

We report on the core-collapse supernova simulation we conducted for a 11.2M progenitor model in three-dimensional space up to 20 ms after bounce, using a radiation hydrodynamics code with full Boltzmann neutrino transport. We solve the six-dimensional Boltzmann equations for three neutrino species and the three-dimensional compressible Euler equations with Furusawa and Togashi's nuclear equation of state. We focus on the prompt convection at ~10 ms after bounce and investigate how neutrinos are transported in the convective matter. We apply a new analysis based on the eigenvalues and eigenvectors of the Eddington tensor and make a comparison between the Boltzmann transport results and the M1 closure approximation in the transition regime between the optically thick and thin limits. We visualize the eigenvalues and eigenvectors using an ellipsoid, in which each principal axis is parallel to one of the eigenvectors and has a length proportional to the corresponding eigenvalue. This approach enables us to understand the difference between the Eddington tensor derived directly from the Boltzmann simulation and the one given by the M1 prescription from a new perspective. We find that the longest principal axis of the ellipsoid is almost always nearly parallel to the energy flux in the M1 closure approximation whereas in the Boltzmann simulation it becomes perpendicular in some transition regions, where the mean free path is ~0.1 times the radius. In three spatial dimensions, the convective motions make it difficult to predict where this happens and possibly improves the closure relation there.

Using the Boltzmann-radiation-hydrodynamics code which solves the Boltzmann equation for the neutrino transport, we present the results of the simulations with the nuclear equations of state (EOSs) of Lattimer-Swesty (LS) and Furusawa-Shen (FS). We extend the simulation time of the LS model and conducted thorough investigations, though the previous paper reported some of the results briefly. Only the LS model shows the shock revival. This seems to originate from the nuclear composition: the different nuclear composition results in the different energy loss by the photodissociation and hence the different strength of the prompt convection and the later neutrino-driven convection. The proto-neutron star seen in the FS model is more compact than that in the LS model since the existence of multi-nuclear species softens the EOS. For the behavior of neutrinos, we examined the flux and the Eddington tensor of neutrinos. In the optically thick region, the diffusion of neutrinos and the dragging by the matter motion determine the flux. In the optically thin region, the free-streaming determines it. The Eddington tensor is compared with that obtained from the M1-closure relation. The M1-closure scheme overestimates the contribution from the velocity-dependent terms in the semi-transparent region.

The small energy exchange via nucleon recoils in neutrino-nucleon scattering is now supposed to be one of the important factors for successful explosion of core-collapse supernovae (CCSNe) as they can change neutrino spectra through accumulation of a large number of scatterings. In finite-difference methods employed for neutrino transport in CCSN simulations, we normally can not afford to deploy a large enough number of energy bins needed to resolve this small energy exchange and sub-grid techniques are employed one way or another. In this paper we study quantitatively with the Monte Carlo (MC) method how well such a treatment performs. We first investigate the effects of nucleon recoils on the neutrino spectra and confirm that the average energy is reduced by $\sim$15% for heavy-lepton neutrinos and by much smaller quantities for other types of neutrinos in a typical post-bounce situation. It is also observed that the nucleon scattering dominates the electron scattering in the thermalization of neutrino spectra in all flavors. We then study possible artifacts that the coarse energy grid may produce in the finite-difference methods. In order to mimic the latter calculation, we re-distribute MC particles in each energy bin after a certain interval in a couple of ways and study how the results are affected and depend on the energy-resolution. We also discuss possible implications of our results for the finite-difference methods.

Neutrinos are believed to have a key role in the explosion mechanism of core-collapse supernovae as they carry most of the energy released by the gravitational collapse of a massive star. If their flavor is converted fast inside the neutrino sphere, the supernova explosion may be influenced. This paper is reporting the results of the extended work of our previous paper. We perform a thorough survey of the ELN crossing in one of our self-consistent, realistic Boltzmann simulations in two spatial dimensions under axisymmetry for the existence of the crossings between $\nu_e$ and $\bar\nu_e$ angular distributions, or the electron lepton number (ELN) crossing. We report for the first time the positive detections deep inside the core of the massive star in the vicinity of neutrino sphere at $r$ $\approx$ 16 - 21 km. We find that low values of the electron fraction $Y_e$ produced by convective motions together with the appearance of light elements are critically important to give rise to the ELN crossing by enhancing the chemical potential difference between proton and neutron, and hence by mitigating the Fermi-degeneracy of $\nu_e$. Since the region of positive detection are sustained and, in fact, expanding with time, it may have an impact on the explosion of core-collapse supernovae, observational neutrino astronomy and nucleosynthesis of heavy nuclei.

Mixing and fallback models in faint supernova models are supposed to reproduce the abundance patterns of observed carbon-enhanced metal-poor (CEMP) stars in the Galactic halo. A fine tuning of the model parameters for individual stars is required to reproduce the observed ratios of carbon to iron. We focus on extremely metal-poor stars formed out of the ejecta from the mixing and fallback models using a chemical evolution model. Our chemical evolution models take into account the contribution of individual stars to chemical enrichment in host halos together with their evolution in the context of the hierarchical clustering. Parametrized models of mixing and fallback models for Pop. III faint supernovae are implemented in the chemical evolution models with merger trees to reproduce the observed CEMP stars. A variety of choices for model parameters on star formation and metal-pollution by faint supernovae is unable to reproduce the observed stars with [Fe/H] < -4 and [C/H] > -2, which are the majority of CEMP stars among the lowest metallicity stars. Only possible solution is to form stars from small ejecta mass, which produces an inconsistent metallicity distribution function. We conclude that not all the CEMP stars are explicable by the mixing and fallback models. We also tested the contribution of binary mass transfers from AGB stars that are also supposed to reproduce the abundances of known CEMP stars. This model reasonably reproduces the distribution of carbon and iron abundances simultaneously only if we assume that long-period binaries are favored at [Fe/H] < -3.5.

We report the discovery of a Compton-thick (CT) dust-obscured galaxy (DOG) at $z$ = 0.89, WISE J082501.48+300257.2 (WISE0825+3002), observed by Nuclear Spectroscopic Telescope Array (NuSTAR). X-ray analysis with the XCLUMPY model revealed that hard X-ray luminosity in the rest-frame 2-10 keV band of WISE0825+3002 is $L_{\rm X}$ (2-10 keV) = $4.2^{+2.8}_{-1.6} \times 10^{44}$ erg s$^{-1}$ while its hydrogen column density is $N_{\rm H}$ = $1.0^{+0.8}_{-0.4} \times 10^{24}$ cm$^{-2}$, indicating that WISE0825+3002 is a mildly CT active galactic nucleus (AGN). We performed the spectral energy distribution (SED) fitting with CIGALE to derive its stellar mass, star formation rate, and infrared luminosity. The estimated Eddington ratio based on stellar mass and integration of the best-fit SED of AGN component is $\lambda_{\rm Edd}$ = 0.70, which suggests that WISE0825+3002 harbors an actively growing black hole behind a large amount of gas and dust. We found that the relationship between luminosity ratio of X-ray and 6 $\mu$m, and Eddington ratio follows an empirical relation for AGNs reported by Toba et al. (2019a).

We present a linear stability analysis of the fast-pairwise neutrino flavor conversion based on a result of our latest axisymmetric core-collapse supernova (CCSN) simulation with full Boltzmann neutrino transport. In the CCSN simulation, coherent asymmetric neutrino emissions of electron-type neutrinos ($\nu_{\rm e}$) and their anti-particles ($\bar{\nu}_{\rm e}$), in which the asymmetry of $\nu_{\rm e}$ and $\bar{\nu}_{\rm e}$ is anti-correlated with each other, occur at almost the same time as the onset of aspherical shock expansion. We find that the asymmetric neutrino emissions play a crucial role on occurrences of fast flavor conversions. The linear analysis shows that unstable modes appear in both pre- and post-shock flows; for the latter they appear only in the hemisphere of higher $\bar{\nu}_{\rm e}$ emissions (the same hemisphere with stronger shock expansion). We analyze in depth the characteristics of electron-lepton-number (ELN) crossing by closely inspecting the angular distributions of neutrinos in momentum space. The ELN crossing happens in various ways, and the property depends on the radius: in the vicinity of neutron star, $\bar{\nu}_{\rm e}$ ($\nu_{\rm e}$) dominates over $\nu_{\rm e}$ ($\bar{\nu}_{\rm e}$) in the forward (backward) direction: at the larger radius the ELN crossing occurs in the opposite way. We also find that the non-radial ELN crossing occurs at the boundary between no ELN crossing and the radial one, which is an effect of genuine multi-D transport. Our findings indicate that the collective neutrino oscillation may occur more commonly in CCSNe and suggest that the CCSN community needs to accommodate these oscillations self-consistently in the modelling of CCSNe.

We make a strong case that the fast neutrino-flavor conversion, one of the collective flavor oscillation modes, commonly occurs in core-collapse supernovae (CCSNe). It is confirmed in the numerical data obtained in realistic simulations of CCSNe but the argument is much more generic and applicable universally: the coherent neutrino-nucleus scattering makes the electron lepton number (ELN) change signs at some inward direction and trigger the flavor conversion in the outward direction in the pre-shock region. Although the ELN crossing is tiny and that is why it has eluded recognition so far, it is still large enough to induce the flavor conversion. Our findings will have an important observational consequences for CCSNe neutrinos.

In this {\it Letter}, we present the result of an axisymmetric CCSN simulation conducted with appropriate treatments of neutrino transport and proper motions of proto-neutron star (PNS), in which a remarkable PNS acceleration is observed in association with asymmetric neutrino emissions by the time of $300$ms after bounce. We find that these asymmetric neutrino emissions play important roles in the acceleration of PNS in this phase. The correlation between the PNS proper motion and the asymmetric ejecta is similar to that in the neutron star (NS) kick of hydrodynamic origin. Both electron-type neutrinos ($\nu_{\rm e}$) and their anti-particles ($\bar{\nu}_{\rm e}$) have $\sim 10\%$ level of asymmetry between the northern and southern hemispheres, while other heavy-leptonic neutrinos ($\nu_x$) have much smaller asymmetry of $\sim 1 \%$. The emissions of $\bar{\nu}_{\rm e}$ and $\nu_x$ are higher in the hemisphere of stronger shock expansion whereas the $\nu_{\rm e}$ emission is enhanced in the opposite hemisphere: in total the neutrinos carry some linear momentum to the hemisphere of the stronger shock expansion. The asymmetry is attributed to the non-spherical distribution of electron-fraction ($Y_e$) in the envelope of PNS. Although it is similar to LESA (lepton-emission self-sustained asymmetry), the $Y_e$ asymmetry seems to be associated with the PNS motion: the latter triggers lateral circular motions in the envelope of PNS by breaking symmetry of the matter distribution there, which are then sustained by a combination of convection, lateral neutrino diffusion and matter-pressure gradient. Our findings may have an influence on the current theories on the NS kick mechanism although long-term simulations are required to assess their impact on the later evolution.

We present a new method for neutrino-matter coupling in multi-dimensional radiation-hydrodynamic simulations of core-collapse supernova (CCSN) with the full Boltzmann neutrino transport. The development is motivated by the fact that the accurate conservation of momentum is required for reliable numerical modelings of CCSN dynamics including a recoil of proto-neutron star (PNS). The new method is built on a hybrid approach, in which we use the energy-momentum tensor of neutrinos to compute the momentum feedback from neutrino to matter in the optically thick region while we employ the collision term in the optically thin region. In this method we utilize a general relativistic description of radiation-hydrodynamics with angular moments, which allows us to evaluate the momentum feedback from neutrino to matter without inconsistency with our Boltzmann solver. We demonstrate that the new method substantially improves the accuracy of linear momentum conservation in our CCSN simulations under reasonable angular resolutions in momentum space, alleviating the difficulty in giving the diffusion limit precisely with the discrete ordinate (Sn) method. It is the first-ever demonstration that the PNS kick can be handled directly and properly in multi-dimensional radiation-hydrodynamic simulations with the full Boltzmann neutrino transport.

We propose an anomaly detection technique for high-resolution X-ray spectroscopy. The method is based on the neural network architecture variational autoencoder, and requires only {\it normal} samples for training. We implement the network using Python taking account of the effect of Poisson statistics carefully, and deonstrate the concept with simulated high-resolution X-ray spectral datasets of one-temperature, two-temperature and non-equilibrium plasma. Our proposed technique would assist scientists in finding important information that would otherwise be missed due to the unmanageable amount of data taken with future X-ray observatories.

Neutrinos are densely populated deep inside the core of massive stars after their gravitational collapse to produce supernova explosions and form compact stars such as neutron stars (NS) and black holes (BH). It has been considered that they may change their flavor identities through so-called fast-pairwise conversions induced by mutual forward scatterings. If that is really the case, the dynamics of supernova explosion will be influenced, since the conversion may occur near the neutrino sphere, from which neutrinos are effectively emitted. In this paper, we conduct a pilot study of such possibilities based on the results of fully self-consistent, realistic simulations of a core-collapse supernova explosion in two spatial dimensions under axisymmetry. As we solved the Boltzmann equations for neutrino transfer in the simulation not as a post-process but in real time, the angular distributions of neutrinos in momentum space for all points in the core at all times are available, a distinct feature of our simulations. We employ some of these distributions extracted at a few selected points and times from the numerical data and apply linear analysis to assess the possibility of the conversion. We focus on the vicinity of the neutrino sphere, where different species of neutrinos move in different directions and have different angular distributions as a result. This is a pilot study for a more thorough survey that will follow soon. We find no positive sign of conversion unfortunately at least for the spatial points and times we studied in this particular model. We hence investigate rather in detail the condition for the conversion by modifying the neutrino distributions rather arbitrarily by hand.

We propose new diagnostics that utilize the [O IV] 25.89 $\mu$m and nuclear (subarcsecond scale) 12 $\mu$m luminosity ratio for identifying whether an AGN is deeply `buried' in their surrounding material. Utilizing a sample of 16 absorbed AGNs at redshifts lower than 0.03 in the Swift/BAT catalog observed with Suzaku, we find that AGNs with small scattering fractions ($<$0.5%) tend to show weaker [O IV]-to-12 $\mu$m luminosity ratios than the average of Seyfert 2 galaxies. This suggests that this ratio is a good indicator for identifying buried AGNs. Then, we apply this criterion to 23 local ultra/luminous infrared galaxies (U/LIRGs) in various merger stages hosting AGNs. We find that AGNs in most of mid- to late-stage mergers are buried, while those in earlier stage ones (including non-merger) are not. This result suggests that the fraction of buried AGNs in U/LIRGs increases as the galaxy-galaxy interaction becomes more significant.

We apply a new X-ray clumpy torus model called XCLUMPY (Tanimoto et al. 2019), where the clump distribution is assumed to be the same as in the infrared clumpy torus model (CLUMPY) by Nenkova et al. 2008, to the broadband X-ray spectra of type-1 active galactic nuclei (AGNs) for the first time. We analyze the archival data of IC 4329A and NGC 7469 observed with NuSTAR/Suzaku and NuSTAR/XMM-Newton, respectively, whose infrared spectra were studied with CLUMPY by Ichikawa et al. 2015 and optical extinctions ($A_{\rm V}$) of the tori were estimated. We consider two models, invoking (Model 1) a relativistic reflection component from the accretion disk and (Model 2) a partial absorber. Assuming that the narrow Fe K$\alpha$ emission line at 6.4 keV originates from the torus, we separate the contribution of the torus reflection components in the total spectra. Our models yield equatorial hydrogen column densities of the tori to be $N^{\rm Equ}_{\rm H} =$ (0.53--1.43) $\times 10^{23}~\rm cm^{-2}$ and $N^{\rm Equ}_{\rm H} =$ (0.84--1.43) $\times 10^{24}~\rm cm^{-2}$, for IC 4329A and NGC 7469, respectively. We find that the $N_{\rm H}/A_{\rm V}$ ratios in the tori are by factors of 25--68 (IC 4329A) and 2.4--3.9 (NGC 7469) smaller than that in the Galactic interstellar medium (ISM). These results suggest that a non-negligible fraction of AGNs are 'dust-rich' compared with the Galactic ISM, as opposite to the general trend previously reported in many obscured AGNs.

With the Boltzmann-radiation-hydrodynamics code, which we have developed to solve numerically the Boltzmann equations for neutrino transfer, the Newtonian hydrodynamics equations, and the Newtonian self-gravity simultaneously and consistently, we simulate the collapse of a rotating core of the progenitor with a zero-age-main-sequence mass of $11.2\,M_\odot$ and a shelluler rotation of $1\,{\rm rad\,s^{-1}}$ at the center. We pay particular attention in this paper to the neutrino distribution in phase space, which is affected by the rotation. By solving the Boltzmann equations directly, we can assess the rotation-induced distortion of the angular distribution in momentum space, which gives rise to the rotational component of the neutrino flux. We compare the Eddington tensors calculated both from the raw data and from the M1-closure approximation. We demonstrate that the Eddington tensor is determined by complicated interplays of the fluid velocity and the neutrino interactions and that the M1-closure, which assumes that the Eddington factor is determined by the flux factor, fails to fully capture this aspect, especially in the vicinity of the shock. We find that the error in the Eddington factor reaches $\sim 20\%$ in our simulation. This is due not to the resolution but to the different dependence of the Eddington and flux factors on the angular profile of the neutrino distribution function, and hence modification to the closure relation is needed.

This paper is a sequel to Takahashi et al. (2016), in which the authors investigated the influences of fluctuations in pre-shock accreting matter on the linear stability of the standing accretion shock in core-collapse supernovae (CCSNe). If one understands that this was concerning the effect of the outer boundary condition for the post-shock accretion flows, the present paper should be regarded as an investigation on the possible influences of the inner boundary conditions. More specifically, we impose a time-dependent, oscillating condition at the inner boundary, having in mind the injection of acoustic power by an oscillating proto-neutron star. We also consider possible correlations between the inner and outer boundary conditions as invoked in the argument for Lepton-number Emission Self-sustained Asymmetry, or LESA. In this paper, we conduct the linear stability analysis of the standing accretion shock commonly encountered in CCSNe based on Laplace transform. We find that the acoustic power enhances the standing accretion shock instability, or SASI, especially when the luminosity is low. On the other hand, the correlation between the fluctuations of neutrino luminosity at the neutrino sphere has little influences on the instability, changing the amplitudes of eigenmodes only slightly. We further investigate steady solution of perturbation equations, being motivated by LESA, and conclude that not the difference but the sum of the fluxes of electron-type neutrinos and anti-neutrino is the key ingredient to production of the self-sustained steady perturbed configuration.

Gamma-ray astronomy explores the most energetic photons in nature to address some of the most pressing puzzles in contemporary astrophysics. It encompasses a wide range of objects and phenomena: stars, supernovae, novae, neutron stars, stellar-mass black holes, nucleosynthesis, the interstellar medium, cosmic rays and relativistic-particle acceleration, and the evolution of galaxies. MeV gamma-rays provide a unique probe of nuclear processes in astronomy, directly measuring radioactive decay, nuclear de-excitation, and positron annihilation. The substantial information carried by gamma-ray photons allows us to see deeper into these objects, the bulk of the power is often emitted at gamma-ray energies, and radioactivity provides a natural physical clock that adds unique information. New science will be driven by time-domain population studies at gamma-ray energies. This science is enabled by next-generation gamma-ray instruments with one to two orders of magnitude better sensitivity, larger sky coverage, and faster cadence than all previous gamma-ray instruments. This transformative capability permits: (a) the accurate identification of the gamma-ray emitting objects and correlations with observations taken at other wavelengths and with other messengers; (b) construction of new gamma-ray maps of the Milky Way and other nearby galaxies where extended regions are distinguished from point sources; and (c) considerable serendipitous science of scarce events -- nearby neutron star mergers, for example. Advances in technology push the performance of new gamma-ray instruments to address a wide set of astrophysical questions.

We perform an extensive study of the influence of nuclear weak interactions on core-collapse supernovae (CCSNe), paying particular attention to consistency between nuclear abundances in the equation of state (EOS) and nuclear weak interactions. We compute properties of uniform matter based on the variational method. For inhomogeneous nuclear matter, we take a full ensemble of nuclei into account with various finite-density and thermal effects and directly use the nuclear abundances to compute nuclear weak interaction rates. To quantify the impact of a consistent treatment of nuclear abundances on CCSN dynamics, we carry out spherically symmetric CCSN simulations with full Boltzmann neutrino transport, systematically changing the treatment of weak interactions, EOSs, and progenitor models. We find that the inconsistent treatment of nuclear abundances between the EOS and weak interaction rates weakens the EOS dependence of both the dynamics and neutrino signals. We also test the validity of two artificial prescriptions for weak interactions of light nuclei and find that both prescriptions affect the dynamics. Furthermore, there are differences in neutrino luminosities by ~10% and in average neutrino energies by 0.25-1 MeV from those of the fiducial model. We also find that the neutronization burst neutrino signal depends on the progenitor more strongly than on the EOS, preventing a detection of this signal from constraining the EOS.

A critical mass ONe core with a high ignition density is considered to end in gravitational collapse leading to neutron star formation. Being distinct from a Fe core collapse, the final evolution involves combustion flame propagation, in which complex phase transition from ONe elements into the nuclear-statistical-equilibrium (NSE) state takes place. We simulate the core evolution from the O+Ne ignition until the bounce shock penetrates the whole core, using a state-of-the-art 1D Lagrangian neutrino-radiation-hydrodynamic code, in which important nuclear burning, electron capture, and neutrino reactions are taken into account. Special care is also taken in making a stable initial condition by importing the stellar EOS, which is used for the progenitor evolution calculation, and by improving the remapping process. We find that the central ignition leads to intense $\nu_e$ radiation with $L_{\nu_e} \gtrsim 10^{51}$ erg s$^{-1}$ powered by fast electron captures onto NSE isotopes. This pre-bounce $\nu_e$ radiation heats the surroundings by the neutrino-electron scattering, which acts as a new driving mechanism of the flame propagation together with the adiabatic contraction. The resulting flame velocity of $\sim10^8$ cm s$^{-1}$ will be more than one-order-of-magnitude faster than that of laminar flame driven by heat conduction. We also find that the duration of the pre-bounce $\nu_e$ radiation phase depends on the degree of the core hydrostatic/dynamical stability. Therefore, the future detection of the pre-bounce neutrino is important not only to discriminate the ONe core collapse from the Fe core collapse but also to constrain the progenitor hydrodynamical stability.

We report X-ray, optical, and near-infrared monitoring of the new X-ray transient MAXI J1820$+$070 discovered with MAXI on 2018 March 11. Its X-ray intensity reached $\sim 2$ Crab in 2--20 keV at the end of March, and then gradually decreased until the middle of June. In this period, the X-ray spectrum was described by Comptonization of the disk emission, with a photon index of $\sim$1.5 and an electron temperature of $\sim$50 keV, which is consistent with a black hole X-ray binary in the low/hard state. The electron temperature and the photon index were slightly decreased and increased with increasing flux, respectively. The source showed significant X-ray flux variation on timescales of seconds. This short-term variation was found to be associated with changes in the spectral shape, and the photon index became slightly harder at higher fluxes. This suggests that the variation was produced by a change in the properties of the hot electron cloud responsible for the strong Comptonization. Modeling a multi-wavelength SED around the X-ray flux peak at the end of March, covering the near-infrared to X-ray bands, we found that the optical and near-infrared fluxes were likely contributed substantially by the jet emission. Before this outburst, the source was never detected in the X-ray band with MAXI (with a 3$\sigma$ upper limit of $\sim$0.2 mCrab in 4--10 keV, obtained from the 7-year data in 2009--2016), whereas weak optical and infrared activity was found at their flux levels $\sim$3 orders of magnitude lower than the peak fluxes in the outburst.

We present the results from the Hitomi Soft Gamma-ray Detector (SGD) observation of the Crab nebula. The main part of SGD is a Compton camera, which in addition to being a spectrometer, is capable of measuring polarization of gamma-ray photons. The Crab nebula is one of the brightest X-ray / gamma-ray sources on the sky, and, the only source from which polarized X-ray photons have been detected. SGD observed the Crab nebula during the initial test observation phase of Hitomi. We performed the data analysis of the SGD observation, the SGD background estimation and the SGD Monte Carlo simulations, and, successfully detected polarized gamma-ray emission from the Crab nebula with only about 5 ks exposure time. The obtained polarization fraction of the phase-integrated Crab emission (sum of pulsar and nebula emissions) is (22.1 $\pm$ 10.6)% and, the polarization angle is 110.7$^o$ + 13.2 / $-$13.0$^o$ in the energy range of 60--160 keV (The errors correspond to the 1 sigma deviation). The confidence level of the polarization detection was 99.3%. The polarization angle measured by SGD is about one sigma deviation with the projected spin axis of the pulsar, 124.0$^o$ $\pm$0.1$^o$.

We investigate axisymmetric steady solutions of (magneto)hydrodynamics equations that describe approximately accretion flows through a standing shock wave and discuss the effects of rotation and magnetic field on the revival of the stalled shock wave in supernova explosions. We develop a new powerful numerical method to calculate the 2-dimensional (2D) steady accretion flows self-consistently. We first confirm the results of preceding papers that there is a critical luminosity of irradiating neutrinos, above which there exists no steady solution in spherical models. If a collapsing star has rotation and/or magnetic field, the accretion flows are no longer spherical owing to the centrifugal force and/or Lorentz force and the critical luminosity is modified.In fact we find that the critical luminosity is reduced by about 50% - 70% for rapid rotations and about 20% - 50% for strong toroidal magnetic fields, depending on the mass accretion rate. These results may be also interpreted as an existence of the critical specific angular momentum or critical magnetic field, above which there exists no steady solution and the standing shock wave will revive for a given combination of mass accretion rate and neutrino luminosity.

We propose a new class of method for solving nonlinear systems of equations, which, among other things,has four nice features: (i) it is inspired by the mathematical property of damped oscillators, (ii) it can be regarded as a simple extention to the Newton-Raphson(NR) method, (iii) it has the same local convergence as the NR method does, (iv) it has a significantly wider convergence region or the global convergence than that of the NR method. In this article, we present the evidence of these properties, applying our new method to some examples and comparing it with the NR method.

We carry out a comprehensive study of supernova ejecta-companion interaction in massive binary systems. We aim to physically understand the kinematics of the interaction and predict observational signatures. To do this we perform simulations over a vast parameter space of binary configurations, varying the masses of the progenitor and companion, structure of the companion, explosion energy, and orbital separation. Our results were not so consistent with classical models by Wheeler et al. 1975, sometimes deviating by an order of magnitude. We construct an alternative simple model which explains the simulated results reasonably well and can be used to estimate impact velocities for arbitrary explosion profiles and companion star structures. We then investigate the long term evolution after the supernova, where the companion can be inflated by the energy injected into the star. We find that the companion can become more than an order of magnitude overluminous straight after the supernova, but quickly fades away after $\sim10$ years and returns to its original luminosity in about a thermal timescale of the star. Finally, we also discuss the possible surface contamination of heavy elements from the slower ejecta.

We present the third MAXI/GSC catalog in the high Galactic-latitude sky ($|b| > 10^\circ$) based on the 7-year data from 2009 August 13 to 2016 July 31, complementary to that in the low Galactic-latitude sky ($|b| < 10^\circ$; Hori et al. 2018). We compile 682 sources detected at significances of $s_{\rm D,4-10~keV} \geq 6.5$ in the 4--10 keV band. A two-dimensional image fit based on the Poisson likelihood algorithm ($C$-statistics) is adopted for the detections and constraints on their fluxes and positions. The 4--10 keV sensitivity reaches $\approx 0.48$ mCrab, or $\approx 5.9 \times 10^{-12}$ erg cm$^{-2}$ s$^{-1}$, over the half of the survey area. Compared with the 37-month catalog (Hiroi et al. 2013), which adopted a threshold of $s_{\rm D,4-10~keV} \geq 7$, the source number increases by a factor of $\sim$1.4. The fluxes in the 3--4 keV and 10--20 keV bands are further estimated, and hardness ratios (HRs) are calculated using the 3--4 keV, 4--10 keV, 3--10 keV, and 10--20 keV band fluxes. We also make the 4--10 keV lightcurves in one year bins for all the sources and characterize their variabilities with an index based on a likelihood function and the excess variance. Possible counterparts are found from five major X-ray survey catalogs by Swift, Uhuru, RXTE, XMM-Newton, and ROSAT, and an X-ray galaxy-cluster catalog (MCXC). Our catalog provides the fluxes, positions, detection significances, HRs, one-year bin lightcurves, variability indices, and counterpart candidates.

We present a broadband (0.4-70 keV) X-ray spectral analysis of the luminous infrared galaxy (LIRG) system Mrk 463 observed with Nuclear Spectroscopic Telescope Array (NuSTAR), Chandra, and XMM-Newton, which contains double active galactic nuclei (AGNs) (Mrk 463E and Mrk 463W) with a separation of $\sim$ 3.8 kpc. Detecting their transmitted hard X-ray $>10$ keV continua with NuSTAR, we confirm that Mrk 463E and Mrk 463W have AGNs with intrinsic luminosities of (1.6-2.2) $\times 10^{43}$ and (0.5-0.6) $\times 10^{43}$ erg s$^{-1}$ (2-10 keV) obscured by hydrogen column densities of 8 $\times 10^{23}$ and 3 $\times 10^{23}$ cm$^{-2}$, respectively. Both nuclei show strong reflection components from cold matter. The luminosity ratio between X-ray (2-10 keV) and [O IV] 25.89 $\mu$m of Mrk 463E is $\sim$5 times smaller than those of normal Seyfert galaxies, suggesting that the intrinsic SED is X-ray weak relative to the UV luminosity. In fact, the bolometric AGN luminosity of Mrk 463E estimated from $L^{\prime}$-band (3.8 $\mu$m), [O IV] 25.89 $\mu$m, and [Ne V] 14.32 $\mu$m lines indicate a large bolometric-to-X-ray luminosity ratio, $\kappa_{\rm 2-10 keV} \approx$ 110-410, and a high Eddington ratio, $\lambda_{\rm Edd} \sim$ 0.4-0.8. We suggest that the merger triggered a rapid growth of the black hole in Mrk 463E, which is not yet deeply "buried" by circumnuclear dust. By contrast, the $L^{\prime}$-band luminosity of Mrk 463W is unusually small relative to the X-ray luminosity, suggesting that the Eddington ratio is low ($<10^{-3}$) and it might be still in an early phase of merger-driven AGN activity.

Collective neutrino oscillations are induced by the presence of neutrinos themselves. As such they are intrinsically nonlinear phenomena and are much more complex than linear counterparts such as the vacuum or MSW oscillations. They obey integro-differential equations, numerical solutions of which are also very challenging. If one focuses on the onset of the collective oscillations, on the other hand, the equations can be linearized and the technique of linear analysis can be employed. Unfortunately, however, it is well known that such an analysis, when applied with discretizations of continuous angular distributions, suffers from the appearance of so-called spurious modes, unphysical eigenmodes of the discretized linear equations. In this paper, we analyze in detail the origin of these unphysical modes and present a simple solution to this annoying problem. We have found that the spurious modes originate from the artificial production of pole singularities instead of a branch cut in the Riemann surface by the discretizations. The branching point singularities in the Riemann surface for the original undiscretized equations can be recovered by approximating the angular distributions with polynomials and then performing the integrals analytically. We demonstrate for some examples that this simple prescription removes the spurious modes indeed. We also propose an even simpler method: a piecewise linear approximation to the angular distribution. It is shown that the same methodology is applicable to the multi-energy case as well as to the dispersion relation approach that was proposed very recently.

We present results from the Hitomi X-ray observation of a young composite-type supernova remnant (SNR) G21.5$-$0.9, whose emission is dominated by the pulsar wind nebula (PWN) contribution. The X-ray spectra in the 0.8-80 keV range obtained with the Soft X-ray Spectrometer (SXS), Soft X-ray Imager (SXI) and Hard X-ray Imager (HXI) show a significant break in the continuum as previously found with the NuSTAR observation. After taking into account all known emissions from the SNR other than the PWN itself, we find that the Hitomi spectra can be fitted with a broken power law with photon indices of $\Gamma_1=1.74\pm0.02$ and $\Gamma_2=2.14\pm0.01$ below and above the break at $7.1\pm0.3$ keV, which is significantly lower than the NuSTAR result ($\sim9.0$ keV). The spectral break cannot be reproduced by time-dependent particle injection one-zone spectral energy distribution models, which strongly indicates that a more complex emission model is needed, as suggested by recent theoretical models. We also search for narrow emission or absorption lines with the SXS, and perform a timing analysis of PSR J1833$-$1034 with the HXI and SGD. No significant pulsation is found from the pulsar. However, unexpectedly, narrow absorption line features are detected in the SXS data at 4.2345 keV and 9.296 keV with a significance of 3.65 $\sigma$. While the origin of these features is not understood, their mere detection opens up a new field of research and was only possible with the high resolution, sensitivity and ability to measure extended sources provided by an X-ray microcalorimeter.

We present the first results of our spatially axisymmetric core-collapse supernova simulations with full Boltzmann neutrino transport, which amount to a time-dependent 5-dimensional (2 in space and 3 in momentum space) problem in fact. Special relativistic effects are fully taken into account with a two-energy-grid technique. We performed two simulations for a progenitor of 11.2M, employing different nuclear equations-of-state (EOS's): Lattimer and Swesty's EOS with the incompressibility of K = 220MeV (LS EOS) and Furusawa's EOS based on the relativistic mean field theory with the TM1 parameter set (FS EOS). In the LS EOS the shock wave reaches ~700km at 300ms after bounce and is still expanding whereas in the FS EOS it stalled at ~200km and has started to recede by the same time. This seems to be due to more vigorous turbulent motions in the former during the entire post-bounce phase, which leads to higher neutrino-heating efficiency in the neutrino-driven convection. We also look into the neutrino distributions in momentum space, which is the advantage of the Boltzmann transport over other approximate methods. We find non-axisymmetric angular distributions with respect to the local radial direction, which also generate off-diagonal components of the Eddington tensor. We find that the r {\theta}-component reaches ~10% of the dominant rr-component and, more importantly, it dictates the evolution of lateral neutrino fluxes, dominating over the {\theta}{\theta}-component, in the semi-transparent region. These data will be useful to further test and possibly improve the prescriptions used in the approximate methods.

We present data preprocessing based on an artificial neural network to estimate the parameters of the X-ray emission spectra of a single-temperature thermal plasma. The method finds appropriate parameters close to the global optimum. The neural network is designed to learn the parameters of the thermal plasma (temperature, abundance, normalisation, and redshift) of the input spectra. After training using 9000 simulated X-ray spectra, the network has grown to predict all the unknown parameters with uncertainties of about a few percent. The performance dependence on the network structure has been studied. We applied the neural network to an actual high-resolution spectrum obtained with {\it Hitomi}. The predicted plasma parameters agreed with the known best-fit parameters of the Perseus cluster within $\lesssim10$\% uncertainties. The result shows a possibility that neural networks trained by simulated data can be useful to extract a feature built in the data, which would reduce human-intensive preprocessing costs before detailed spectral analysis, and help us make the best use of large quantities of spectral data coming in the next decades.

The QED vacuum polarization in external monochromatic plane-wave electromagnetic fields is calculated with spatial and temporal variations of the external fields being taken into account. We develop a perturbation theory to calculate the induced electromagnetic current that appears in the Maxwell equations, based on Schwinger's proper-time method, and combine it with the so-called gradient expansion to handle the variation of external fields perturbatively. The crossed field, i.e., the long wavelength limit of the electromagnetic wave is first considered. The eigenmodes and the refractive indices as the eigenvalues associated with the eigenmodes are computed numerically for the probe photon propagating in some particular directions. In so doing, no limitation is imposed on the field strength and the photon energy unlike previous studies. It is shown that the real part of the refractive index becomes less than unity for strong fields, the phenomenon that has been known to occur for high-energy probe photons. We then evaluate numerically the lowest-order corrections to the crossed-field resulting from the field variations in space and time. It is demonstrated that the corrections occur mainly in the imaginary part of the refractive index.

The present paper investigates the temperature structure of the X-ray emitting plasma in the core of the Perseus cluster using the 1.8--20.0 keV data obtained with the Soft X-ray Spectrometer (SXS) onboard the Hitomi Observatory. A series of four observations were carried out, with a total effective exposure time of 338 ks and covering a central region $\sim7'$ in diameter. The SXS was operated with an energy resolution of $\sim$5 eV (full width at half maximum) at 5.9 keV. Not only fine structures of K-shell lines in He-like ions but also transitions from higher principal quantum numbers are clearly resolved from Si through Fe. This enables us to perform temperature diagnostics using the line ratios of Si, S, Ar, Ca, and Fe, and to provide the first direct measurement of the excitation temperature and ionization temperature in the Perseus cluster. The observed spectrum is roughly reproduced by a single temperature thermal plasma model in collisional ionization equilibrium, but detailed line ratio diagnostics reveal slight deviations from this approximation. In particular, the data exhibit an apparent trend of increasing ionization temperature with increasing atomic mass, as well as small differences between the ionization and excitation temperatures for Fe, the only element for which both temperatures can be measured. The best-fit two-temperature models suggest a combination of 3 and 5 keV gas, which is consistent with the idea that the observed small deviations from a single temperature approximation are due to the effects of projection of the known radial temperature gradient in the cluster core along the line of sight. Comparison with the Chandra/ACIS and the XMM-Newton/RGS results on the other hand suggests that additional lower-temperature components are present in the ICM but not detectable by Hitomi SXS given its 1.8--20 keV energy band.

The Hitomi SXS spectrum of the Perseus cluster, with $\sim$5 eV resolution in the 2-9 keV band, offers an unprecedented benchmark of the atomic modeling and database for hot collisional plasmas. It reveals both successes and challenges of the current atomic codes. The latest versions of AtomDB/APEC (3.0.8), SPEX (3.03.00), and CHIANTI (8.0) all provide reasonable fits to the broad-band spectrum, and are in close agreement on best-fit temperature, emission measure, and abundances of a few elements such as Ni. For the Fe abundance, the APEC and SPEX measurements differ by 16%, which is 17 times higher than the statistical uncertainty. This is mostly attributed to the differences in adopted collisional excitation and dielectronic recombination rates of the strongest emission lines. We further investigate and compare the sensitivity of the derived physical parameters to the astrophysical source modeling and instrumental effects. The Hitomi results show that an accurate atomic code is as important as the astrophysical modeling and instrumental calibration aspects. Substantial updates of atomic databases and targeted laboratory measurements are needed to get the current codes ready for the data from the next Hitomi-level mission.

In this paper, we systematically calculate the polarization in soft X-rays emitted from magnetized neutron stars, which are expected to be observed by the next-generation X-ray satellites. Magnetars are one of the targets for these observations. This is because thermal radiation is normally observed in the soft X-ray band, and it is thought to be linearly polarized because of different opacities for two polarization modes of photons in the magnetized atmosphere of neutron stars and the dielectric properties of the vacuum in strong magnetic fields. In their previous study, Taverna et al. illustrated how strong magnetic fields influence the behavior of the polarization observables for radiation propagating in vacuo without addressing a precise, physical emission model. In this paper, we pay attention to the conversion of photon polarization modes that can occur in the presence of an atmospheric layer above the neutron star surface, computing the polarization angle and fraction and systematically changing the magnetic field strength, radii of the emission region, temperature, mass, and radii of the neutron stars. We confirmed that if plasma is present, the effects of mode conversion cannot be neglected when the magnetic field is relatively weak, $B \sim 10^{13} \mathrm{G}$. Our results indicate that strongly magnetized ($B \gtrsim 10^{14} \mathrm{G}$) neutron stars are suitable to detect polarizations, but not-so-strongly magnetized ($B \sim 10^{13} \mathrm{G}$) neutron stars will be the ones to confirm the mode conversion.

We present Hitomi observations of N132D, a young, X-ray bright, O-rich core-collapse supernova remnant in the Large Magellanic Cloud (LMC). Despite a very short observation of only 3.7 ks, the Soft X-ray Spectrometer (SXS) easily detects the line complexes of highly ionized S K and Fe K with 16-17 counts in each. The Fe feature is measured for the first time at high spectral resolution. Based on the plausible assumption that the Fe K emission is dominated by He-like ions, we find that the material responsible for this Fe emission is highly redshifted at ~800 km/s compared to the local LMC interstellar medium (ISM), with a 90% credible interval of 50-1500 km/s if a weakly informative prior is placed on possible line broadening. This indicates (1) that the Fe emission arises from the supernova ejecta, and (2) that these ejecta are highly asymmetric, since no blue-shifted component is found. The S K velocity is consistent with the local LMC ISM, and is likely from swept-up ISM material. These results are consistent with spatial mapping that shows the He-like Fe concentrated in the interior of the remnant and the S tracing the outer shell. The results also show that even with a very small number of counts, direct velocity measurements from Doppler-shifted lines detected in extended objects like supernova remnants are now possible. Thanks to the very low SXS background of ~1 event per spectral resolution element per 100 ks, such results are obtainable during short pointed or slew observations with similar instruments. This highlights the power of high-spectral-resolution imaging observations, and demonstrates the new window that has been opened with Hitomi and will be greatly widened with future missions such as the X-ray Astronomy Recovery Mission (XARM) and Athena.

We report a Hitomi observation of IGR J16318-4848, a high-mass X-ray binary system with an extremely strong absorption of N_H~10^{24} cm^{-2}. Previous X-ray studies revealed that its spectrum is dominated by strong fluorescence lines of Fe as well as continuum emission. For physical and geometrical insight into the nature of the reprocessing material, we utilize the high spectroscopic resolving power of the X-ray microcalorimeter (the soft X-ray spectrometer; SXS) and the wide-band sensitivity by the soft and hard X-ray imager (SXI and HXI) aboard Hitomi. Even though photon counts are limited due to unintended off-axis pointing, the SXS spectrum resolves Fe K{\alpha_1} and K{\alpha_2} lines and puts strong constraints on the line centroid and width. The line width corresponds to the velocity of 160^{+300}_{-70} km s^{-1}. This represents the most accurate, and smallest, width measurement of this line made so far from any X-ray binary, much less than the Doppler broadening and shift expected from speeds which are characteristic of similar systems. Combined with the K-shell edge energy measured by the SXI and HXI spectra, the ionization state of Fe is estimated to be in the range of Fe I--IV. Considering the estimated ionization parameter and the distance between the X-ray source and the absorber, the density and thickness of the materials are estimated. The extraordinarily strong absorption and the absence of a Compton shoulder component is confirmed. These characteristics suggest reprocessing materials which are distributed in a narrow solid angle or scattering primarily with warm free electrons or neutral hydrogen.

The origin of the narrow Fe-K{\alpha} fluorescence line at 6.4 keV from active galactic nuclei has long been under debate; some of the possible sites are the outer accretion disk, the broad line region, a molecular torus, or interstellar/intracluster media. In February-March 2016, we performed the first X-ray microcalorimeter spectroscopy with the Soft X-ray Spectrometer (SXS) onboard the Hitomi satellite of the Fanaroff-Riley type I radio galaxy NGC 1275 at the center of the Perseus cluster of galaxies. With the high energy resolution of ~5 eV at 6 keV achieved by Hitomi/SXS, we detected the Fe-K{\alpha} line with ~5.4 {\sigma} significance. The velocity width is constrained to be 500-1600 km s$^{-1}$ (FWHM for Gaussian models) at 90% confidence. The SXS also constrains the continuum level from the NGC 1275 nucleus up to ~20 keV, giving an equivalent width ~20 eV of the 6.4 keV line. Because the velocity width is narrower than that of broad H{\alpha} line of ~2750 km s$^{-1}$, we can exclude a large contribution to the line flux from the accretion disk and the broad line region. Furthermore, we performed pixel map analyses on the Hitomi/SXS data and image analyses on the Chandra archival data, and revealed that the Fe-K{\alpha} line comes from a region within ~1.6 kpc from the NGC 1275 core, where an active galactic nucleus emission dominates, rather than that from intracluster media. Therefore, we suggest that the source of the Fe-K{\alpha} line from NGC 1275 is likely a low-covering fraction molecular torus or a rotating molecular disk which probably extends from a pc to hundreds pc scale in the active galactic nucleus system.

Extending the earlier measurements reported in Hitomi collaboration (2016, Nature, 535, 117), we examine the atmospheric gas motions within the central 100~kpc of the Perseus cluster using observations obtained with the Hitomi satellite. After correcting for the point spread function of the telescope and using optically thin emission lines, we find that the line-of-sight velocity dispersion of the hot gas is remarkably low and mostly uniform. The velocity dispersion reaches maxima of approximately 200~km~s$^{-1}$ toward the central active galactic nucleus (AGN) and toward the AGN inflated north-western `ghost' bubble. Elsewhere within the observed region, the velocity dispersion appears constant around 100~km~s$^{-1}$. We also detect a velocity gradient with a 100~km~s$^{-1}$ amplitude across the cluster core, consistent with large-scale sloshing of the core gas. If the observed gas motions are isotropic, the kinetic pressure support is less than 10\% of the thermal pressure support in the cluster core. The well-resolved optically thin emission lines have Gaussian shapes, indicating that the turbulent driving scale is likely below 100~kpc, which is consistent with the size of the AGN jet inflated bubbles. We also report the first measurement of the ion temperature in the intracluster medium, which we find to be consistent with the electron temperature. In addition, we present a new measurement of the redshift to the brightest cluster galaxy NGC~1275.

Thanks to its high spectral resolution (~5 eV at 6 keV), the Soft X-ray Spectrometer (SXS) on board Hitomi enables us to measure the detailed structure of spatially resolved emission lines from highly ionized ions in galaxy clusters for the first time. In this series of papers, using the SXS we have measured the velocities of gas motions, metallicities and the multi-temperature structure of the gas in the core of the Perseus cluster. Here, we show that when inferring physical properties from line emissivities in systems like Perseus, the resonant scattering (RS) effect should be taken into account. In the Hitomi waveband, RS mostly affects the FeXXV He$\alpha$ line ($w$) - the strongest line in the spectrum. The flux measured by Hitomi in this line is suppressed by a factor ~1.3 in the inner ~30 kpc, compared to predictions for an optically thin plasma; the suppression decreases with the distance from the center. The $w$ line also appears slightly broader than other lines from the same ion. The observed distortions of the $w$ line flux, shape and distance dependence are all consistent with the expected effect of the resonant scattering in the Perseus core. By measuring the ratio of fluxes in optically thick ($w$) and thin (FeXXV forbidden, He$\beta$, Ly$\alpha$) lines, and comparing these ratios with predictions from Monte Carlo radiative transfer simulations, the velocities of gas motions have been obtained. The results are consistent with the direct measurements of gas velocities from line broadening described elsewhere in this series, although the systematic and statistical uncertainties remain significant. Further improvements in the predictions of line emissivities in plasma models, and deeper observations with future X-ray missions will enable RS measurements to provide powerful constraints on the amplitude and anisotropy of clusters gas motions.

This paper is a sequel to our previous one (Kato et al.2015), which calculated the luminosities and spectra of electron-type anti-neutrinos ($\bar{\nu}_e$'s) from the progenitors of core-collapse supernovae. Expecting that a capability to detect electron-type neutrinos ($\nu_e$'s) will increase dramatically with the emergence of liquid-argon detectors such as DUNE, we broaden the scope in this study to include all-flavors of neutrinos emitted from the pre-bounce phase. We pick up three progenitor models of an electron capture supernova (ECSN) and iron-core collapse supernovae (FeCCSNe). We find that the number luminosities reach $\sim10^{57} \mathrm{s^{-1}}$ and $\sim10^{53} \mathrm{s^{-1}}$ at maximum for $\nu_e$ and $\bar{\nu}_e$, respectively. We also estimate the numbers of detection events at terrestrial neutrino detectors including DUNE, taking flavor oscillations into account and assuming the distance to the progenitors to be 200 pc. It is demonstrated that $\bar{\nu}_e$'s from the ECSN-progenitor will be undetected at almost all detectors, whereas we will be able to observe $\gtrsim$15900 $\nu_e$'s at DUNE for the inverted mass hierarchy. From the FeCCSN-progenitors, the number of $\bar{\nu}_e$ events will be largest for JUNO, 200-900 $\bar{\nu}_e$'s, depending on the mass hierarchy whereas the number of $\nu_e$ events at DUNE is $\gtrsim$2100 for the inverted mass hierarchy. These results imply that the detection of $\bar{\nu}_e$'s is useful to distinguish FeCCSN- from ECSN-progenitors, while $\nu_e$'s will provide us with detailed information on the collapse phase regardless of the type and mass of progenitor.

The mechanism driving core-collapse supernovae is sensitive to the interplay between matter and neutrino radiation. However, neutrino radiation transport is very difficult to simulate, and several radiation transport methods of varying levels of approximation are available. We carefully compare for the first time in multiple spatial dimensions the discrete ordinates (DO) code of Nagakura, Yamada, and Sumiyoshi and the Monte Carlo (MC) code Sedonu, under the assumptions of a static fluid background, flat spacetime, elastic scattering, and full special relativity. We find remarkably good agreement in all spectral, angular, and fluid interaction quantities, lending confidence to both methods. The DO method excels in determining the heating and cooling rates in the optically thick region. The MC method predicts sharper angular features due to the effectively infinite angular resolution, but struggles to drive down noise in quantities where subtractive cancellation is prevalent, such as the net gain in the protoneutron star and off-diagonal components of the Eddington tensor. We also find that errors in the angular moments of the distribution functions induced by neglecting velocity dependence are sub-dominant to those from limited momentum-space resolution. We briefly compare directly computed second angular moments to those predicted by popular algebraic two-moment closures, and find that the errors from the approximate closures are comparable to the difference between the DO and MC methods. Included in this work is an improved Sedonu code, which now implements a fully special relativistic, time-independent version of the grid-agnostic Monte Carlo random walk approximation.

We propose a mechanism of state transition in wind-fed black hole binaries (high-mass X-ray binaries) such as Cyg X-1 and LMC X-1. Modeling a line-driven stellar wind from the companion by two-dimensional hydrodynamical calculations, we investigate the processes of wind capture by and accretion onto the black hole. We assume that the wind acceleration is terminated at the He II ionization front because ions responsible for line-driven acceleration are ionized within the front, i.e. He III region. It is found that the mass accretion rate inferred from the luminosity is remarkably smaller than the capture rate. Considering the difference, we construct a model for the state transition based on the accretion flow being controlled by magneto-rotational instability. The outer flow is torus like, and plays an important role to trigger the transition. The model can explain why state transition does occur in Cyg X-1, while not in LMC X-1. Cyg X-1 exhibits a relatively low luminosity, and then the He II ionization front is located and can move between the companion and black hole, depending on its ionizing photon flux. On the other hand, LMC X-1 exhibits too high luminosity for the front to move considerably; the front is too close to the companion atmosphere. The model also predicts that each state of high-soft or low-hard would last fairly long because the luminosity depends weakly on the wind velocity. In the context of the model, the state transition is triggered by a fluctuation of the magnetic field when its amplitude becomes comparable to the field strength in the torus-like outer flow.

We developed a three-dimensional radiative transfer code for an ultra-relativistic background flow-field by using the Monte Carlo (MC) method in the context of gamma-ray burst (GRB) emission. For obtaining reliable simulation results in the coupled computation of MC radiation transport with relativistic hydrodynamics which can reproduce GRB emission, we validated radiative transfer computation in the ultra-relativistic regime and assessed the appropriate simulation conditions. The radiative transfer code was validated through two test calculations: (1) computing in different inertial frames and (2) computing in flow-fields with discontinuous and smeared shock fronts. The simulation results of the angular distribution and spectrum were compared among three different inertial frames and in good agreement with each other. If the time duration for updating the flow-field was sufficiently small to resolve a mean free path of a photon into ten steps, the results were thoroughly converged. The spectrum computed in the flow-field with a discontinuous shock front obeyed a power-law in frequency whose index was positive in the range from 1 to 10 MeV. The number of photons in the high-energy side decreased with the smeared shock front because the photons were less scattered immediately behind the shock wave due to the small electron number density. The large optical depth near the shock front was needed for obtaining high-energy photons through bulk Compton scattering. Even one-dimensional structure of the shock wave could affect the results of radiation transport computation.

To search for giant X-ray pulses correlated with the giant radio pulses (GRPs) from the Crab pulsar, we performed a simultaneous observation of the Crab pulsar with the X-ray satellite Hitomi in the 2 -- 300 keV band and the Kashima NICT radio observatory in the 1.4 -- 1.7 GHz band with a net exposure of about 2 ks on 25 March 2016, just before the loss of the Hitomi mission.The timing performance of the Hitomi instruments was confirmed to meet the timing requirement and about 1,000 and 100 GRPs were simultaneously observed at the main and inter-pulse phases, respectively, and we found no apparent correlation between the giant radio pulses and the X-ray emission in either the main or inter-pulse phases.All variations are within the 2 sigma fluctuations of the X-ray fluxes at the pulse peaks, and the 3 sigma upper limits of variations of main- or inter- pulse GRPs are 22\% or 80\% of the peak flux in a 0.20 phase width, respectively, in the 2 -- 300 keV band.The values become 25\% or 110\% for main or inter-pulse GRPs, respectively, when the phase width is restricted into the 0.03 phase.Among the upper limits from the Hitomi satellite, those in the 4.5-10 keV and the 70-300 keV are obtained for the first time, and those in other bands are consistent with previous reports.Numerically, the upper limits of main- and inter-pulse GRPs in the 0.20 phase width are about (2.4 and 9.3) $\times 10^{-11}$ erg cm$^{-2}$, respectively. No significant variability in pulse profiles implies that the GRPs originated from a local place within the magnetosphere and the number of photon-emitting particles temporally increases.However, the results do not statistically rule out variations correlated with the GRPs, because the possible X-ray enhancement may appear due to a $>0.02$\% brightening of the pulse-peak flux under such conditions.

We construct new equations of state for baryons at sub-nuclear densities for the use in core-collapse supernova simulations. The abundance of various nuclei is obtained together with thermodynamic quantities. The formulation is an extension of the previous model, in which we adopted the relativistic mean field theory with the TM1 parameter set for nucleons, the quantum approach for $d$, $t$, $h$ and $\alpha$ as well as the liquid drop model for the other nuclei under the nuclear statistical equilibrium. We reformulate the model of the light nuclei other than $d$, $t$, $h$ and $\alpha$ based on the quasi-particle description. Furthermore, we modify the model so that the temperature dependences of surface and shell energies of heavy nuclei could be taken into account. The pasta phases for heavy nuclei and the Pauli- and self-energy shifts for $d$, $t$, $h$ and $\alpha$ are taken into account in the same way as in the previous model. We find that nuclear composition is considerably affected by the modifications in this work, whereas thermodynamical quantities are not changed much. In particular, the washout of shell effect has a great impact on the mass distribution above $T \sim 1$ MeV. This improvement may have an important effect on the rates of electron captures and coherent neutrino scatterings on nuclei in supernova cores.

We have constructed a nuclear equation of state (EOS) that includes a full nuclear ensemble for use in core-collapse supernova simulations. It is based on the EOS for uniform nuclear matter that two of the authors derived recently, applying a variational method to realistic two- and there-body nuclear forces. We utilize an extended liquid drop model of heavy nuclei and a quantum-theoretical mass evaluation for light nuclei. In addition to realistic nuclear forces, the inclusion of in-medium effects on the full ensemble of nuclei makes the new EOS one of the most realistic EOS's for supernova simulations. We make comparisons with the FYSS EOS, which is based on the same formulation for the nuclear ensemble but adopts the relativistic mean field (RMF) theory with the TM1 parameter set for uniform nuclear matter. The new EOS is softer than the FYSS EOS around and above nuclear saturation densities. We find that neutron-rich nuclei with small mass numbers are more abundant in the new EOS than in the FYSS EOS because of the larger saturation densities and smaller symmetry energy of nuclei in the former. We apply the two EOS's to 1D supernova simulations and find that the new EOS gives lower electron fractions and higher temperatures in the collapse phase owing to the smaller symmetry energy. As a result, the inner core has smaller masses for the new EOS. It is more compact, on the other hand, due to the softness of the new EOS and bounces at higher densities. The ensuing outward propagations of the shock wave in the outer core are very similar in the two simulations, which may be an artifact, though, caused by the use of the same tabulated electron capture rates for heavy nuclei ignoring differences in the nuclear composition between the two EOS's in these computations.

The Crab nebula originated from a core-collapse supernova (SN) explosion observed in 1054 A.D. When viewed as a supernova remnant (SNR), it has an anomalously low observed ejecta mass and kinetic energy for an Fe-core collapse SN. Intensive searches were made for a massive shell that solves this discrepancy, but none has been detected. An alternative idea is that the SN1054 is an electron-capture (EC) explosion with a lower explosion energy by an order of magnitude than Fe-core collapse SNe. In the X-rays, imaging searches were performed for the plasma emission from the shell in the Crab outskirts to set a stringent upper limit to the X-ray emitting mass. However, the extreme brightness of the source hampers access to its vicinity. We thus employed spectroscopic technique using the X-ray micro-calorimeter onboard the Hitomi satellite. By exploiting its superb energy resolution, we set an upper limit for emission or absorption features from yet undetected thermal plasma in the 2-12 keV range. We also re-evaluated the existing Chandra and XMM-Newton data. By assembling these results, a new upper limit was obtained for the X-ray plasma mass of <~ 1Mo for a wide range of assumed shell radius, size, and plasma temperature both in and out of the collisional equilibrium. To compare with the observation, we further performed hydrodynamic simulations of the Crab SNR for two SN models (Fe-core versus EC) under two SN environments (uniform ISM versus progenitor wind). We found that the observed mass limit can be compatible with both SN models if the SN environment has a low density of <~ 0.03 cm-3 (Fe core) or <~ 0.1 cm-3 (EC) for the uniform density, or a progenitor wind density somewhat less than that provided by a mass loss rate of 10-5 Mo yr-1 at 20 km s-1 for the wind environment.

We have constructed the database of stars in the local group using the extended version of the SAGA (Stellar Abundances for Galactic Archaeology) database that contains stars in 24 dwarf spheroidal galaxies and ultra faint dwarfs. The new version of the database includes more than 4500 stars in the Milky Way, by removing the previous metallicity criterion of [Fe/H] <= -2.5, and more than 6000 stars in the local group galaxies. We examined a validity of using a combined data set for elemental abundances. We also checked a consistency between the derived distances to individual stars and those to galaxies in the literature values. Using the updated database, the characteristics of stars in dwarf galaxies are discussed. Our statistical analyses of alpha-element abundances show that the change of the slope of the [alpha/Fe] relative to [Fe/H] (so-called "knee") occurs at [Fe/H] = -1.0+-0.1 for the Milky Way. The knee positions for selected galaxies are derived by applying the same method. Star formation history of individual galaxies are explored using the slope of the cumulative metallicity distribution function. Radial gradients along the four directions are inspected in six galaxies where we find no direction dependence of metallicity gradients along the major and minor axes. The compilation of all the available data shows a lack of CEMP-s population in dwarf galaxies, while there may be some CEMP-no stars at [Fe/H] <~ -3 even in the very small sample. The inspection of the relationship between Eu and Ba abundances confirms an anomalously Ba-rich population in Fornax, which indicates a pre-enrichment of interstellar gas with r-process elements. We do not find any evidence of anti-correlations in O-Na and Mg-Al abundances, which characterises the abundance trends in the Galactic globular clusters.

We investigate the criterion for the acoustic mechanism to work successfully in core-collapse supernovae. The acoustic mechanism is an alternative to the neutrino-heating mechanism. It was proposed by Burrows et al., who claimed that acoustic waves emitted by $g$-mode oscillations in proto-neutron stars (PNS) energize a stalled shock wave and eventually induce an explosion. Previous works mainly studied to which extent the $g$-modes are excited in the PNS. In this paper, on the other hand, we investigate how strong the acoustic wave needs to be if it were to revive a stalled shock wave. By adding the acoustic power as a new axis, we draw a critical surface, an extension of the critical curve commonly employed in the context of neutrino heating. We perform both 1D and 2D parametrized simulations, in which we inject acoustic waves from the inner boundary. In order to quantify the power of acoustic waves, we use the extended Myers theory to take neutrino reactions into proper account. We find for the 1D simulations that rather large acoustic powers are required to relaunch the shock wave, since the additional heating provided by the secondary shocks developed from acoustic waves is partially canceled by the neutrino cooling that is also enhanced. In 2D, the required acoustic powers are consistent with those of Burrows et al. Our results seem to imply, however, that it is the sum of neutrino heating and acoustic powers that matters for shock revival.

We investigate the influences of the nuclear composition on the weak interaction rates of heavy nuclei during the core collapse of massive stars. The nuclear abundances in nuclear statistical equilibrium (NSE) are calculated by some equation of state (EOS) models including in-medium effects on nuclear masses. We systematically examine the sensitivities of electron capture and neutrino-nucleus scattering on heavy nuclei to the nuclear shell effects and the single nucleus approximation. We find that the washout of shell effects at high temperatures brings significant change to weak rates by smoothing the nuclear abundance distribution: the electron capture rate decreases by $\sim$20$\%$ in the early phase and increases by $\sim$40$\%$ in the late phase at most, while the cross section for neutrino-nucleus scattering is reduced by $\sim$15$\%$. This is because the open-shell nuclei become abundant instead of those with closed neutron shells as the shell effects disappear. We also find that the single-nucleus description based on the average values leads to underestimations of weak rates. Electron captures and neutrino coherent scattering on heavy nuclei are reduced by $\sim$80$\%$ in the early phase and by $\sim$5$\%$ in the late phase, respectively. These results indicate that NSE like EOS accounting for shell washout is indispensable for the reliable estimation of weak interaction rates in simulations of core-collapse supernovae.

High-resolution X-ray spectroscopy with Hitomi was expected to resolve the origin of the faint unidentified E=3.5 keV emission line reported in several low-resolution studies of various massive systems, such as galaxies and clusters, including the Perseus cluster. We have analyzed the Hitomi first-light observation of the Perseus cluster. The emission line expected for Perseus based on the XMM-Newton signal from the large cluster sample under the dark matter decay scenario is too faint to be detectable in the Hitomi data. However, the previously reported 3.5 keV flux from Perseus was anomalously high compared to the sample-based prediction. We find no unidentified line at the reported high flux level. Taking into account the XMM measurement uncertainties for this region, the inconsistency with Hitomi is at a 99% significance for a broad dark-matter line and at 99.7% for a narrow line from the gas. We do not find anomalously high fluxes of the nearby faint K line or the Ar satellite line that were proposed as explanations for the earlier 3.5 keV detections. We do find a hint of a broad excess near the energies of high-n transitions of Sxvi (E=3.44 keV rest-frame) -- a possible signature of charge exchange in the molecular nebula and another proposed explanation for the unidentified line. While its energy is consistent with XMM pn detections, it is unlikely to explain the MOS signal. A confirmation of this interesting feature has to wait for a more sensitive observation with a future calorimeter experiment.

We investigated r-process nucleosynthesis in magneto-rotational supernovae, based on a new explosion mechanism induced by the magneto-rotational instability. A series of axisymmetric magneto-hydrodynamical simulations with detailed microphysics including neutrino heating is performed, numerically resolving the magneto-rotational instability. Neutrino-heating dominated explosions, enhanced by magnetic fields, showed mildly neutron-rich ejecta producing nuclei up to $A \sim 130$ (i.e. the weak r-process), while explosion models with stronger magnetic fields reproduce a solar-like r-process pattern. More commonly seen abundance patterns in our models are in between the weak and regular r-process, producing lighter and intermediate mass nuclei. These {\it intermediate r-processes} exhibit a variety of abundance distributions, compatible with several abundance patterns in r-process-enhanced metal-poor stars. The amount of Eu ejecta $\sim 10^{-5} M_\odot$ in magnetically-driven jets agrees with predicted values in the chemical evolution of early galaxies. In contrast, neutrino-heating dominated explosions have a significant amount of Fe ($^{56}{\rm Ni}$) and Zn, comparable to regular supernovae and hypernovae, respectively. These results indicate magneto-rotational supernovae can produce a wide range of heavy nuclei from iron-group to r-process elements, depending on the explosion dynamics.

We show the softest ever spectrum from Cyg~X-1, detected in 2013 with Suzaku. This has the weakest high energy Compton tail ever seen from this object, so should give the cleanest view of the underlying disk spectrum, and hence the best determination of black hole spin from disk continuum fitting. Using the standard model of a disk with simple non-thermal Comptonisation to produce the weak high energy tail gives a high spin black hole. However, we get a significantly better fit by including an additional, low temperature thermal Comptonisation component, which allows a much lower black hole spin. Corroboration of the existence of an additional Compton component comes from the frequency dependent hard lags seen in the rapid variability in archival high/soft state data. These can not be explained if the continuum is a single non-thermal Comptonisation component, but are instead consistent with a radially stratified, multi zone Comptonisation spectrum, where the spectrum is softer further from the black hole. A complex multi-zone Comptonisation continuum is required to explain both spectra and timing together, and this has an impact on the derived black hole spin.

We analyze the ${\it Kepler}$ monitoring light curve of a blazar W2R 1926$+$42 to examine features of microvariability by means of the "shot analysis" technique. We select 195 intra-day, flare-like variations (shots) for the continuous light curve of Quarter 14 with a duration of 100 d. In the application of the shot analysis, an averaged profile of variations is assumed to converge with a universal profile which reflects a physical mechanism generating the microvariability in a blazar jet, although light-variation profiles of selected shots show a variety. A mean profile, which is obtained by aligning the peaks of the 195 shots, is composed of a spiky-shape shot component at $\pm$0.1 d (with respect to the time of the peak), and two slow varying components ranging from $-$0.50 d to $-$0.15 d and from 0.10 d to 0.45 d of the peak time. The former spiky feature is well represented by an exponential rise of 0.043$\pm$0.001 d and an exponential decay of 0.061$\pm$0.002 d. These timescales are consistent with that corresponding to a break frequency of a power spectrum density calculated from the obtained light curve. After verification with the Monte-Carlo method, the exponential shape, but not the observed asymmetry, of the shot component can be explained by noise variation. The asymmetry is difficult to explain through a geometrical effect (i.e. changes of the geometry of the emitting region), but is more likely to be caused by the production and dissipation of high-energy accelerated particles in the jet. Additionally, durations of the detected shots show a systematic variation with a dispersion caused by a statistical randomness. A comparison with the variability of Cygnus X-1 is also briefly discussed.

The origin of carbon-enhanced metal-poor (CEMP) stars plays a key role in characterising the formation and evolution of the first stars and the Galaxy since the extremely-poor (EMP) stars with [Fe/H] \leq -2.5 share the common features of carbon enhancement in their surface chemical compositions. The origin of these stars is not yet established due to the controversy of the origin of CEMP stars without the enhancement of s-process element abundances, i.e., so called CEMP-no stars. In this paper, we elaborate the s-process nucleosynthesis in the EMP AGB stars and explore the origin of CEMP stars. We find that the efficiency of the s-process is controlled by O rather than Fe at [Fe/H] \lesssim -2. We demonstrate that the relative abundances of Sr, Ba, Pb to C are explained in terms of the wind accretion from AGB stars in binary systems.

A detailed search for emission and absorption lines and assessing their upper limits are performed for Suzaku data. The method utilizes a matched-filtering approach to maximize the signal-to-noise ratio for a given energy resolution, which could be applicable to many types of line search. We first applied it to well-known AGN spectra that have been reported to have ultra-fast outflows, and find that our results are consistent with previous findings at the ~3{\sigma} level. We proceeded to search for emission and absorption features in the two bright magnetars 4U 0142+61 and 1RXS J1708-4009, applying the filtering method to Suzaku data. We found that neither source showed any significant indication of line features, even using long Suzaku observations and dividing their spectra into spin phases. The upper limits on the equivalent width of emission/absorption lines are constrained to be a few eV at ~1 keV, and a few hundreds of eV at ~10 keV. This strengthens previous reports that persistently bright magnetars do not show proton cyclotron absorption features in soft X-rays and, even if they exist, they would be broadened or much weaker than below the detection limit of X-ray CCD.

We have developed a new formulation to obtain self-gravitating, axisymmetric configurations in permanent rotation. The formulation is based on the Lagrangian variational principle with a triangulated mesh. It treats not only barotropic but also baroclinic equations of state. We compare the various stellar equilibria obtained by our new scheme with those by Hachisu's self-consistent field scheme for the barotropic case, and those by Fujisawa's self-consistent field scheme for the baroclinic case. Included in these rotational configurations are those with shellular-type rotations, which are commonly assumed in the evolution calculation of rotating stars. Although radiation processes, convections and meridional flows have not been taken into account in this study, we have in mind the application of this method to the two-dimensional evolution calculations of rotating stars, for which the Lagrangian formulation is best suited.

Clusters of galaxies are the most massive gravitationally-bound objects in the Universe and are still forming. They are thus important probes of cosmological parameters and a host of astrophysical processes. Knowledge of the dynamics of the pervasive hot gas, which dominates in mass over stars in a cluster, is a crucial missing ingredient. It can enable new insights into mechanical energy injection by the central supermassive black hole and the use of hydrostatic equilibrium for the determination of cluster masses. X-rays from the core of the Perseus cluster are emitted by the 50 million K diffuse hot plasma filling its gravitational potential well. The Active Galactic Nucleus of the central galaxy NGC1275 is pumping jetted energy into the surrounding intracluster medium, creating buoyant bubbles filled with relativistic plasma. These likely induce motions in the intracluster medium and heat the inner gas preventing runaway radiative cooling; a process known as Active Galactic Nucleus Feedback. Here we report on Hitomi X-ray observations of the Perseus cluster core, which reveal a remarkably quiescent atmosphere where the gas has a line-of-sight velocity dispersion of 164+/-10 km/s in a region 30-60 kpc from the central nucleus. A gradient in the line-of-sight velocity of 150+/-70 km/s is found across the 60 kpc image of the cluster core. Turbulent pressure support in the gas is 4% or less of the thermodynamic pressure, with large scale shear at most doubling that estimate. We infer that total cluster masses determined from hydrostatic equilibrium in the central regions need little correction for turbulent pressure.

The explosion mechanism of core-collapse supernovae has not been fully understood yet but multi-dimensional fluid instabilities such as standing accretion shock instability (SASI) and convection are now believed to be crucial for shock revival. Another multi-dimensional effect that has been recently argued is the asymmetric structures in progenitors, which are induced by violent convections in silicon/oxygen layers that occur before the onset of collapse, as revealed by recent numerical simulations of the last stage of massive star evolutions. Furthermore, it has been also demonstrated numerically that accretions of such non-spherical envelopes could facilitate shock revival. These two multi-dimensional may hence hold a key to successful explosions. In this paper, we performed a linear stability analysis of the standing accretion shock in core-collapse supernovae, taking into account non-spherical, unsteady accretion flows onto the shock to clarify the possible links between the two effects. We found that such pre-shock perturbations can excite the fluid instabilities efficiently and hence help the shock revive in core-collapse supernovae.

We propose to employ the circular polarization of gravitational waves emitted by core-collapse supernovae as an unequivocal indication of rapid rotation deep in their cores. It has been demonstrated by three dimensional simulations that non-axisymmetric accretion flows may develop spontaneously via hydrodynamical instabilities in the post-bounce cores. It is not surprising then that the gravitational waves emitted by such fluid motions are circularly polarized. We show in this letter that a network of the second generation detectors of gravitational waves worldwide may be able to detect such polarizations up to the opposite side of Galaxy as long as the rotation period is shorter than a few seconds prior to collapse.

From 2013 April to 2014 April, we performed an X-ray and optical simultaneous monitoring of the type 1.5 Seyfert galaxy NGC 3516. It employed Suzaku, and 5 Japanese ground-based telescopes, the Pirka, Kiso Schmidt, Nayuta, MITSuME, and the Kanata telescopes. The Suzaku observations were conducted seven times with various intervals ranging from days, weeks, to months, with an exposure of $\sim50$ ksec each. The optical $B$-band observations not only covered those of Suzaku almost simultaneously, but also followed the source as frequently as possible. As a result, NGC 3516 was found in its faint phase with the 2-10 keV flux of $0.21-2.70 \times 10^{-11}$ erg s$^{-1}$ cm$^{-2}$. The 2-45 keV X-ray spectra were composed of a dominant variable hard power-law continuum with a photon index of $\sim1.7$, and a non-relativistic reflection component with a prominent Fe-K$\alpha$ emission line. Producing the $B$-band light curve by differential image photometry, we found that the $B$-band flux changed by $\sim2.7 \times 10^{-11}$ erg s$^{-1}$ cm$^{-2}$, which is comparable to the X-ray variation, and detected a significant flux correlation between the hard power-law component in X-rays and the $B$-band radiation, for the first time in NGC 3516. By examining their correlation, we found that the X-ray flux preceded that of $B$ band by $2.0^{+0.7}_{-0.6}$ days ($1\sigma$ error). Although this result supports the X-ray reprocessing model, the derived lag is too large to be explained by the standard view which assumes a "lamppost"-type X-ray illuminator located near a standard accretion disk. Our results are better explained by assuming a hot accretion flow and a truncated disk.

We present a newly developed moving-mesh technique for the multi-dimensional Boltzmann-Hydro code for the simulation of core-collapse supernovae (CCSNe). What makes this technique different from others is the fact that it treats not only hydrodynamics but also neutrino transfer in the language of the 3+1 formalism of general relativity (GR), making use of the shift vector to specify the time evolution of the coordinate system. This means that the transport part of our code is essentially general relativistic although in this paper it is applied only to the moving curvilinear coordinates in the flat Minknowski spacetime, since the gravity part is still Newtonian. The numerical aspect of the implementation is also described in detail. Employing the axisymmetric two-dimensional version of the code, we conduct two test computations: oscillations and runaways of proto-neutron star (PNS). We show that our new method works fine, tracking the motions of PNS correctly. We believe that this is a major advancement toward the realistic simulation of CCSNe.

The bright and highly-variable Seyfert I Active Galactic Nucleus, IC4329A, was observed with Suzaku 5 times in 2007 August with intervals of ~ 5 days, and a net exposure of 24-31 ks each. Another longer observation was carried out in 2012 August with a net exposure of 118 ks. In the 6 observations, the source was detected in 2-45 keV with average 2-10 keV fluxes of $(0.67-1.2)\times10^{-10}$ erg cm$^{-2}$ s$^{-1}$. The intensity changed by a factor of 2 among the 5 observations in 2007, and 1.5 within the 2012 observation. Difference spectra among these observations revealed that the variability of IC4329A was carried mainly by a power-law component with a photon index $\Gamma\sim2.0$. However, in addition to this primary component and the associated reflection, the broad-band Suzaku data required another, harder, and less variable component with $\Gamma\sim1.4$. The presence of this new continuum was also confirmed by analyzing the same 6 data sets through the spectral decomposition technique developed by Noda et al. 2013a. This $\Gamma\sim1.4$ continuum is considered to be a new primary component that has not been recognized in the spectra of IC4329A so far, although it was recently identified in those of several other Seyfert I galaxies (Noda et al. 2013a, 2014).

Utilizing an archived Suzaku data acquired on 2007 December 25 for 46 ks, X-ray spectroscopic properties of the dipping and eclipsing low-mass X-ray binary EXO 0748$-$676 were studied. At an assumed distance of 7.1 kpc, the data gave a persistent unabsorbed luminosity of $3.4\times10^{36}$ erg cm$^{-2}$ s$^{-1}$ in 0.6 $-$ 55 keV. The source was in a relatively bright low/hard state, wherein the 0.6 $-$ 55 keV spectrum can be successfully explained by a "double-seed" Comptonization model, incorporating a common corona with an electron temperature of $\sim13$ keV. The seed photons are thought to be supplied from both the neutron star surface, and a cooler truncated disk. Compared to a sample of non-dipping low-mass X-ray binaries in the low/hard state, the spectrum is subject to stronger Comptonization, with a relatively larger Comptonizing $y$-parameter of $\sim1.4$ and a larger coronal optical depth of $\sim5$. This result, when attributed to the high inclination of EXO 0748$-$676, suggests that the Comptonizing corona may elongate along the disk plane, and give a longer path for the seed photons when viewed from edge-on inclinations.

We visually inspected the light curves of 7557 Kepler Objects of Interest (KOIs) to search for single transit events (STEs) possibly due to long-period giant planets. We identified 28 STEs in 24 KOIs, among which 14 events are newly reported in this paper. We estimate the radius and orbital period of the objects causing STEs by fitting the STE light curves simultaneously with the transits of the other planets in the system or with the prior information on the host star density. As a result, we found that STEs in seven of those systems are consistent with Neptune- to Jupiter-sized objects of orbital periods ranging from a few to $\sim$ $20\,\mathrm{yr}$. We also estimate that $\gtrsim20\%$ of the compact multi-transiting systems host cool giant planets with periods $\gtrsim 3\,\mathrm{yr}$ on the basis of their occurrence in the KOIs with multiple candidates, assuming the small mutual inclination between inner and outer planetary orbits.

We study transitions of hadronic matter (HM) to 3-flavor quark matter (3QM), regarding the conversion processes as combustion and describing them hydrodynamically. Under the assumption that HM is metastable with their free energies being larger than those of 3QM but smaller than those of 2-flavor quark matter (2QM), we consider in this paper the conversion induced by diffusions of seed 3QM. This is a sequel to our previous paper, in which the shock-induced conversion was studied in the same frame work. We not only pay attention to the jump condition on both sides of the conversion front but the structures inside the front are also considered by taking into account what happens during the conversion processes on the time scale of weak interactions. We employ for HM the Shen's EOS, which is based on the relativistic mean field theory, and the bag model-based EOS for QM just as in the previous paper. We demonstrated in that paper that in this combination of EOS's the combustion will occur for a wide range of the bag constant and strong coupling constant in the so-called endothermic regime, in which the Hugoniot curve for combustion runs below the initial state. We find that weak deflagration nearly always occurs and that weak detonation is possible only when the diffusion constant is (unrealistically) large and the critical strange fraction is small. The velocities of the conversion front are ~ $10^3-10^7$ cm/s depending on the initial temperature and density as well as the parameters in the QM EOS and become particularly small when the final state is in the mixed phase. Finally we study linear stability of the laminar weak-deflagration front and find that it is unstable in the exothermic regime (Darrius-Landau instability) but stable in the endothermic regime, which is quite contrary to the ordinary combustions.

We carried out 2D-axisymmetric MHD simulations of core-collapse supernovae for rapidly-rotating magnetized progenitors. By changing both the strength of the magnetic field and the spatial resolution, the evolution of the magnetorotational instability (MRI) and its impacts upon the dynamics are investigated. We found that the MRI greatly amplifies the seed magnetic fields in the regime where not the Alfv\'en mode but the buoyant mode plays a primary role in the exponential growth phase. The MRI indeed has a powerful impact on the supernova dynamics. It makes the shock expansion faster and the explosion more energetic, with some models being accompanied by the collimated-jet formations. These effects, however, are not made by the magnetic pressure except for the collimated-jet formations. The angular momentum transfer induced by the MRI causes the expansion of the heating region, by which the accreting matter gain an additional time to be heated by neutrinos. The MRI also drifts low-$Y_p$ matter from the deep inside of the core to the heating region, which makes the net neutrino heating rate larger by the reduction of the cooling due to the electron capture. These two effects enhance the efficiency of the neutrino heating, which is found to be the key to boost the explosion. Indeed we found that our models explode far more weakly when the net neutrino heating is switched off. The contribution of the neutrino heating to the explosion energy could reach 60\% even in the case of strongest magnetic field in the current simulations.

We study transitions of hadronic matter (HM) to 3-flavor quark matter (3QM) locally, regarding the conversion processes as combustion and describing them hydrodynamically. Not only the jump condition on both sides of the conversion front but the structures inside the front are also considered by taking into account what happens during the conversion processes on the time scale of weak interactions as well as equations of state (EOS's) in the mixed phase. Under the assumption that HM is metastable with their free energies being larger than those of 3QM but smaller than those of 2-flavor quark matter (2QM), we consider the transition via 2QM triggered by a rapid density rise in a shock wave. Based on the results, we discuss which combustion modes (strong/weak detonation) may be realized. HM is described by an EOS based on the relativistic mean field theory and 2, 3QM's are approximated by the MIT bag model. We demonstrate for a wide range of bag constant and strong coupling constant in this combination of EOS's that the combustion may occur in the so-called endothermic regime, in which the Hugoniot curve for combustion runs below the one for the shock wave in P-V plane, and which has no terrestrial counter part. We find that strong detonation always occurs. Depending on the EOS of quark matter (QM) as well as the density of HM and the Mach number of the detonation front, deconfinement from HM to 2QM is either completed or not completed in the shock wave. In the latter case, which is more likely if the EOS of QM ensures that deconfinement occurs above the nuclear saturation density and that the maximum mass of cold quark stars is larger than two solar mass, the conversion continues further via the mixing state of HM and 3QM on the time scale of weak interactions.

By performing neutrino-radiation hydrodynamic simulations in spherical symmetry (1D) and axial symmetry (2D) with different progenitor models by Woosley & Heger (2007) from 12 $M_{\odot}$ to 100 $M_{\odot}$, we find that all 1D runs fail to produce an explosion and several 2D runs succeed. The difference in the shock evolutions for different progenitors can be interpreted by the difference in their mass accretion histories, which are in turn determined by the density structures of progenitors. The mass accretion history has two phases in the majority of the models: the earlier phase in which the mass accretion rate is high and rapidly decreasing and the later phase with a low and almost constant accretion rate. They are separated by the so-called turning point, the origin of which is a change of the accreting layer. We argue that shock revival will most likely occur around the turning point and hence that its location in the $\dot M$-$L_\nu$ plane will be a good measure for the possibility of shock revival: if the turning point lies above the critical curve and the system stays there for a long time, shock revival will obtain. In addition, we develop a phenomenological model to approximately evaluate the trajectories in the $\dot M$-$L_\nu$ plane, which, after calibrating free parameters by a small number of 1D simulations, reproduces the location of the turning point reasonably well by using the initial density structure of progenitor alone. We suggest the application of the phenomenological model to a large collection of progenitors in order to infer without simulations which ones are more likely to explode.

Energetics of nuclear reaction is fundamentally important to understand the mechanism of pair instability supernovae (PISNe). Based on the hydrodynamic equations and thermodynamic relations, we derive exact expressions for energy conservation suitable to be solved in simulation. We also show that some formulae commonly used in the literature are obtained as approximations of the exact expressions. We simulate the evolution of very massive stars of ~100-320 Msun with zero- and 1/10 Zsun, and calculate further explosions as PISNe, applying each of the exact and approximate formulae. The calculations demonstrate that the explosion properties of PISN, such as the mass range, the 56Ni yield, and the explosion energy, are significantly affected by applying the different energy generation rates. We discuss how these results affect the estimate of the PISN detection rate, which depends on the theoretical predictions of such explosion properties.

We report a measurement of the neutrino-electron elastic scattering rate of 862 keV 7Be solar neutrinos based on a 165.4 kton-day exposure of KamLAND. The observed rate is 582 +/- 90 (kton-day)^-1, which corresponds to a 862 keV 7Be solar neutrino flux of (3.26 +/- 0.50) x 10^9 cm^-2s^-1, assuming a pure electron flavor flux. Comparing this flux with the standard solar model prediction and further assuming three flavor mixing, a nu_e survival probability of 0.66 +/- 0.14 is determined from the KamLAND data. Utilizing a global three flavor oscillation analysis, we obtain a total 7Be solar neutrino flux of (5.82 +/- 0.98) x 10^9 cm^-2s^-1, which is consistent with the standard solar model predictions.

Aiming to distinguish two types of progenitors of core collapse supernovae, i.e., one with a core composed mainly of oxygen and neon (abbreviated as ONe core) and the other with an iron core (or Fe core), we calculated the luminosities and spectra of neutrinos emitted from these cores prior to gravitational collapse, taking neutrino oscillation into account. We found that the total energies emitted as $\bar{\nu}_e$ from the ONe core are $\lesssim 10^{46}\ {\rm erg}$, which is much smaller than $\sim 10^{47}\ {\rm erg}$ for Fe cores. The average energy, on the other hand, is twice as large for the ONe core as those for the Fe cores. The neutrinos produced by the plasmon decays in the ONe core are more numerous than those from the electron-positron annihilation in both cores but they have much lower average energies $\lesssim 1\ {\rm MeV}$. Although it is difficult to detect the pre-supernova neutrinos from the ONe core even if it is located within 200$\ $pc from the earth, we expect $ \sim 9 - 43$ and $\sim 7 - 61$ events for Fe cores at KamLAND and Super-Kamiokande, respectively, depending on the progenitor mass and neutrino-mass hierarchy. These numbers might be increased by an order of magnitude if we envisage next-generation detectors such as JUNO. We will hence be able to distinguish the two types of progenitors by the detection or non-detection of the pre-supernova neutrinos if they are close enough ($\lesssim1\ {\rm kpc}$).

We investigate the possible effects of the supernova ejecta hitting the companion star in iPTF 13bvn, focusing on the observable features when it becomes visible. iPTF 13bvn is a type Ib supernova that may become the first case that its progenitor is identified as a binary by near future observations. According to calculations by Bersten et al. (2014), the progenitor should have a mass $\approx3.5M_\odot$ to reproduce the supernova light curve, and such compact stars could only be produced via binary evolution. This is one of the reasons that we expect the progenitor to be a binary, but it should be confirmed by observing the remaining companion after the supernova. Their evolutionary calculations suggest that the companion star will be an overluminous OB star at the moment of supernova. With a combination of hydrodynamical and evolutionary simulations, we find that the secondary star will be heated by the supernova ejecta and expand to have larger luminosities and lower surface effective temperatures. The star will look rather like a red super giant, and this should be taken into account when searching for the companion star in the supernova ejecta in future observations.

DIOS (Diffuse Intergalactic Oxygen Surveyor) is a small satellite aiming for a launch around 2020 with JAXA's Epsilon rocket. Its main aim is a search for warm-hot intergalactic medium with high-resolution X-ray spectroscopy of redshifted emission lines from OVII and OVIII ions. The superior energy resolution of TES microcalorimeters combined with a very wide field of view (30--50 arcmin diameter) will enable us to look into gas dynamics of cosmic plasmas in a wide range of spatial scales from Earth's magnetosphere to unvirialized regions of clusters of galaxies. Mechanical and thermal design of the spacecraft and development of the TES calorimeter system are described. We also consider revising the payload design to optimize the scientific capability allowed by the boundary conditions of the small mission.

This paper reports the results of Suzaku observation of the spectral variation of the black hole binary LMCX-1 in the soft state. The observationwas carried out in 2009 from July 21 to 24. the obtained net count rate was $\sim$30 counts s$^{-1}$ in the 0.5--50 keV band with $\sim$10% peak-to-peak flux variation. The time-averaged X-ray spectrum cannot be described by a multi-color disk and single Compton component with its reflection, but requires additional Comptonized emissions. This double Compton component model allows a slightly larger inner radius of the multi-color disk, implying a lower spin parameter. Significant spectral evolution was observed above 8 keV along with a flux decrease on a timescale of $\sim$10$^4$--10$^5$ s. By spectral fitting, we show that this behavior is well explained by changes in the hard Comptonized emission component in contrast to the maintained disk and soft Comptonized emission.

We report MAXI and Swift observations of short-term spectral softenings of the galactic black-hole X-ray binary Swift J1753.5-0127 in the low/hard state. These softening events are characterized by a simultaneous increase of soft X-rays (2-4 keV) and a decrease of hard X-rays (15-50 keV) lasting for a few tens of days. The X-ray energy spectra during the softening periods can be reproduced with a model consisting of a multi-color disk blackbody and its Comptonized component. The fraction of the Comptonized component decreased from 0.30 to 0.15 when the spectrum became softer; meanwhile the inner disk temperature (Tin) increased from 0.2 to 0.45 keV. These results imply that the softening events are triggered by a short-term increase of the mass accretion rate. During the observed spectral softening events, the disk flux (F) and Tin did not obey the relation: F is proportional to Tin^4, suggesting that the inner disk radius does not reach the innermost stable circular orbit.

Blazars show rapid and violent variabilities, which timescale are often less than a day. We studied intraday variations by applying a shot analysis technique to Kepler monitoring of blazar W2R 1926+42 in Quarter 14. We obtained a mean profile calculated from 195 rapid variations. The mean profile shows three components; one is a sharp structure distributed within $\pm$0.1 day of the peak, and two slow-varying components. This spiky-peak component reflects features of rapid variations directly. The profile of peak component shows an exponential rise and decay of which timescales are different, 0.0416 and 0.0588 day respectively. This component is too sharp to represent a standard function which is often used to express blazar variations. This asymmetric profile at the peak is difficult to be explained by a simple variation of the Doppler factor by changing a geometry of the emitting region. This result indicates that intraday variations arise from a production of high-energy accelerated particles in the jet.

The joint JAXA/NASA ASTRO-H mission is the sixth in a series of highly successful X-ray missions developed by the Institute of Space and Astronautical Science (ISAS), with a planned launch in 2015. The ASTRO-H mission is equipped with a suite of sensitive instruments with the highest energy resolution ever achieved at E > 3 keV and a wide energy range spanning four decades in energy from soft X-rays to gamma-rays. The simultaneous broad band pass, coupled with the high spectral resolution of Delta E < 7 eV of the micro-calorimeter, will enable a wide variety of important science themes to be pursued. ASTRO-H is expected to provide breakthrough results in scientific areas as diverse as the large-scale structure of the Universe and its evolution, the behavior of matter in the gravitational strong field regime, the physical conditions in sites of cosmic-ray acceleration, and the distribution of dark matter in galaxy clusters at different redshifts.

Thanks to extensive observations with X-ray missions and facilities working in other wavelengths, as well as rapidly--advancing numerical simulations of accretion flows, our knowledge of astrophysical black holes has been remarkably enriched. Rapid progress has opened new areas of enquiry, including measurements of black hole spin, the properties and driving mechanisms of jets and disk winds, the impact of feedback into local environments, the origin of periodic and aperiodic X-ray variations, and the nature of super-Eddington accretion flows, among others. The goal of this White Paper is to illustrate how ASTRO-H can make dramatic progress in the study of astrophysical black holes, particularly the study of black hole X-ray binaries.

We have developed a new formulation to obtain self-gravitating, axisymmetric configurations in permanent rotation. The formulation is based on the Lagrangian variational principle, and treats not only barotropic but also baroclinic equations of state, for which angular momentum distributions are not necessarily cylindrical. We adopt a Monte Carlo technique, which is analogous to those employed in other fields, e.g. nuclear physics, in minimizing the energy functional, which is evaluated on a triangulated mesh. This letter is a proof of principle and detailed comparisons with existing results will be reported in the sequel, but some test calculations are presented, in which we have achieved an error of $O(10^{-4})$ in the Virial relation. We have in mind the application of this method to two-dimensional calculations of the evolutions of rotating stars, for which the Lagrangian formulation is best suited.

We study the multi-dimensional properties of neutrino transfer inside supernova cores by solving the Boltzmann equations for neutrino distribution functions in genuinely six dimensional (6D) phase space. Adopting representative snapshots of the post-bounce core from other supernova simulations in three dimensions, we solve the temporal evolutions to stationary states of neutrino distribution functions by our Boltzmann solver. Taking advantage of the multi-angle and multi-energy feature realized by the S$_n$ method in our code, we reveal the genuine characteristics of spatially three dimensional (3D) neutrino transfer such as non-radial fluxes and non-diagonal Eddington tensors. In addition, we assess the ray-by-ray approximation, turning off the lateral-transport terms in our code. We demonstrate that the ray-by-ray approximation tends to propagate fluctuations in thermodynamical states around the neutrino-sphere along each radial ray and overestimate the variations between the neutrino distributions on different radial rays. We find that the difference in the densities and fluxes of neutrinos between the ray-by-ray approximation and the full Boltzmann transport becomes ~20%, which is also the case for the local heating rate, whereas the volume-integrated heating rate in the Boltzmann transport is found to be only slightly larger (~2%) than the counterpart in the ray-by-ray approximation due to cancellation among different rays. These results suggest that we had better assess carefully the possible influences of various approximations in the neutrino transfer employed in the current simulations on supernova dynamics. Detailed information on the angle and energy moments of neutrino distribution functions will be profitable for the future development of numerical methods in neutrino-radiation hydrodynamics.

We present X-ray observations of the active galactic nucleus (AGN) in NGC 4785. The source is a local Seyfert 2 which has not been studied so far in much detail. It was recently detected with high significance in the 15-60 keV band in the 66 month Swift/BAT all sky survey, but there have been no prior pointed X-ray observations of this object. With Suzaku, we clearly detect the source below 10 keV, and find it to have a flat continuum and prominent neutral iron fluorescence line with equivalent width >~1 keV. Fitting the broadband spectra with physical reflection models shows the source to be a bona fide Compton thick AGN with Nh of at least 2x10^{24} cm^{-2} and absorption-corrected 2-10 keV X-ray power L(2-10) ~ few times 10^{42} erg s^{-1}. Realistic uncertainties on L(2-10) computed from the joint confidence interval on the intrinsic power law continuum photon index and normalization are at least a factor of 10. The local bona fide Compton thick AGN population is highly heterogeneous in terms of WISE mid-infrared source colours, and the nucleus of NGC 4785 appears especially sub-dominant in the mid-infrared when comparing to other Compton thick AGN. Such sources would not be easily found using mid-infrared selection alone. The extent of host galaxy extinction to the nucleus is not clear, though NGC 4785 shows a complex core with a double bar and inner disk, adding to the list of known Compton thick AGN in barred host galaxies.

We analyzed the growth of non-spherical perturbations in supersonic accretion flows. We have in mind the application to the post-bounce phase of core-collapse supernovae (CCSNe). Such non-spherical perturbations have been suggested by a series of papers by Arnett, who has numerically investigated violent convections in the outer layers of pre-collapse stars. Moreover, Couch & Ott (2013) demonstrated in their numerical simulations that such perturbations may lead to a successful supernova even for a progenitor that fails to explode without the fluctuations. This study investigated the linear growth of perturbations during the infall onto a stalled shock wave. The linearized equations are solved as an initial and boundary value problem with the use of Laplace transform. The background is a Bondi accretion flow whose parameters are chosen to mimic the 15 $\mathrm{M_\odot}$ progenitor model by Woosley & Heger (2007), which is supposed to be a typical progenitor of CCSNe. We found that the perturbations that are given at a large radius grow as they flow down to the shock radius; the density perturbations can be amplified by a factor of 30, for example. We analytically showed that the growth rate is proportional to $l$, the index of the spherical harmonics. We also found that the perturbations oscillate in time with frequencies that are similar to those of the standing accretion shock instability. This may have an implication for shock revival in CCSNe, which will be investigated in our forthcoming paper in more detail.

The majority of massive stars are formed in binary systems. It is hence reasonable to expect that most core-collapse supernovae (CCSNe) take place in binaries and the existence of a companion star may leave some imprints in observed features. Having this in mind, we have conducted two-dimensional hydrodynamical simulations of the collisions of CCSNe ejecta with the companion star in an almost-equal-mass ($\sim 10M_\odot$) binary to find out possible consequences of such events. In particular we pay attention to the amount of mass removed and its dependence on the binary separation. In contrast to the previous surmise, we find that the companion mass is stripped not by momentum transfer but by shock heating. Up to $25\%$ of the original mass can be removed for the closest separations and the removed mass decreases as $M_{ub} \propto a^{-4.3}$ with the binary separation $a$. By performing some experimental computations with artificially-modified densities of incident ejecta, we show that if the velocity of ejecta is fixed, the density of incident ejecta is the single important parameter that actually determines the removed mass as $M_{ub} \propto \rho_{ej} ^{1.4} $. On the other hand, another set of simulations with modified velocities of incident ejecta demonstrate that the strength of the forward shock, which heats up the stellar material and causes the mass loss of the companion star, is actually the key parameter for the removed mass.

We propose a novel numerical method for solving multi-dimensional, special relativistic Boltzmann equations for neutrinos coupled to hydrodynamics equations. It is meant to be applied to simulations of core-collapse supernovae. We handle special relativity in a non-conventional way, taking account of all orders of v/c. Consistent treatment of advection and collision terms in the Boltzmann equations is the source of difficulties, which we overcome by employing two different energy grids: Lagrangian remapped and laboratory fixed grids. We conduct a series of basic tests and perform a one-dimensional simulation of core-collapse, bounce and shock-stall for a 15M_{sun} progenitor model with a minimum but essential set of microphysics. We demonstrate in the latter simulation that our new code is capable of handling all phases in core-collapse supernova. For comparison, a non-relativistic simulation is also conducted with the same code, and we show that they produce qualitatively wrong results in neutrino transfer. Finally, we discuss a possible incorporation of general relativistic effects in our method.

The type I Seyfert galaxy NGC 3227 was observed by Suzaku six times in 2008, with intervals of $\sim1$ week and net exposures of $\sim50$ ksec each. Among the six observations, the source varied by nearly an order of magnitude, being brightest in the 1st observation with a 2-10 keV luminosity of $1.2\times10^{42}$~erg~s$^{-1}$, while faintest in the 4th with $2.9\times10^{41}$~erg~s$^{-1}$. As it became fainter, the continuum in a 2-45 keV band became harder, while a narrow Fe-K$\alpha$ emission line, detected on all occasions at 6.4 keV of the source rest frame, remained approximately constant in the photon flux. Through a method of variability-assisted broad-band spectroscopy (e.g., Noda et al. 2013), the 2-45 keV spectrum of NGC 3227 was decomposed into three distinct components. One is a relatively soft power-law continuum with a photon index of $\sim 2.3$, weakly absorbed and highly variable on time scales of $\sim5$ ksec; it was observed only when the source was above a threshold luminosity of $\sim6.6 \times10^{41}$ erg s$^{-1}$ (in 2-10 keV), and was responsible for further source brightening beyond. Another is a harder and more absorbed continuum with a photon index of $\sim 1.6$, which persisted through the six observations and varied slowly on time scales of a few weeks by a factor of $\sim2$. This component, carrying a major fraction of the broad-band emission when the source is below the threshold luminosity, is considered as an additional primary emission. The last one is a reflection component with the narrow iron line, produced at large distances from the central black hole.

The effect of rotation on the explosion of core-collapse supernovae is investigated systematically in three-dimensional simulations. In order to obtain the critical conditions for explosion as a function of mass accretion rate, neutrino luminosity, and specific angular momentum, rigidly rotating matter was injected from the outer boundary with an angular momentum, which is increased every 500 ms. It is found that there is a critical value of the specific angular momentum, above which the standing shock wave revives, for a given combination of mass accretion rate and neutrino luminosity, i.e. an explosion can occur by rotation even if the neutrino luminosity is lower than the critical value for a given mass accretion rate in non-rotational models. The coupling of rotation and hydrodynamical instabilities plays an important role to characterize the dynamics of shock revival for the range of specific angular momentum that are supposed to be realistic. Contrary to expectations from past studies, the most rapidly expanding direction of the shock wave is not aligned with the rotation axis. Being perpendicular to the rotation axis on average, it can be oriented in various directions. Its dispersion is small when the spiral mode of the standing accretion shock instability (SASI) governs the dynamics, while it is large when neutrino-driven convection is dominant. As a result of the comparison between 2D and 3D rotational models, it is found that m=!0 modes of neutrino-driven convection or SASI are important for shock revival around the critical surface.

The bright type I Seyfert galaxy NGC 3516 was observed by {\it Suzaku} twice, in 2005 October 12--15 and 2009 October 28--November 2, for a gross time coverage of 242 and 544 ksec and a net exposure of 134 and 255 ksec, respectively. The 2--10 keV luminosity was $2.8 \times 10^{41}$ erg s$^{-1}$ in 2005, and $1.6 \times 10^{41}$ erg s$^{-1}$ in 2009. The 1.4--1.7 keV and 2--10 keV count rates both exhibited peak-to-peak variations by a factor of $\sim2$ in 2005, while $\sim 4$ in 2009. In either observation, the 15--45 keV count rate was less variable. The 2--10 keV spectrum in 2005 was significantly more convex than that in 2009. Through a count-count-plot technique, the 2--45 keV signals in both data were successfully decomposed in a model-independent way into two distinct broadband components. One is a variable emission with a featureless spectral shape, and the other is a non-varying hard component accompanied by a prominent Fe-K emission line at 6.33 keV (6.40 keV in the rest frame). The former was fitted successfully by an absorbed power-law model, while the latter requires a new hard continuum in addition to a reflection component from distant materials. The spectral and variability differences between the two observations are mainly attributed to long-term changes of this new hard continuum, which was stable on time scales of several hundreds ksec.

We report on the results from Suzaku observations of the Galactic black hole X-ray binary H 1743-322 in the low/hard state during its outburst in 2012 October. We appropriately take into account the effects of dust-scattering to accurately analyze the X-ray spectra. The time-averaged spectra in the 1-200 keV band are dominated by a hard power-law component of a photon index of \approx 1.6 with a high-energy cutoff at \approx 60 keV, which is well described with Comptonization of the disk emission by hot corona. We estimate the inner disk radius from the multi-color disk component, and find that it is 1.3-2.3 times larger than the radius in the high/soft state. This suggests that the standard disk was not extended to the innermost stable circular orbit. A reflection component from the disk is detected with R = \Omega/2\pi \approx 0.6 (\Omega is the solid angle). We also successfully estimate the stable disk component in a way independent of the time-averaged spectral modeling, by analyzing short-term spectral variability on the \sim 1-sec timescale. A weak low-frequency quasi-periodic oscillation (LF QPO) at 0.1-0.2 Hz is detected, whose frequency is found to correlate with the X-ray luminosity and photon index. This result may be explained by the evolution of the disk truncation radius.

The symbiotic X-ray binary 4U 1954+319 is a rare system hosting a peculiar neutron star (NS) and an M-type optical companion. Its ~5.4h NS spin period is the longest among all known accretion-powered pulsars and exhibited large (~7%) fluctuations over 8 years. A spin trend transition was detected with Swift/BAT around an X-ray brightening in 2012. The source was in quiescent and bright states before and after this outburst based on 60 ks Suzaku observations in 2011 and 2012. The observed continuum is well described by a Comptonized model with the addition of a narrow 6.4 keV Fe Kalpha line during the outburst. Spectral similarities to slowly rotating pulsars in high-mass X-ray binaries, its high pulsed fraction (~60-80%), and the location in the Corbet diagram favor high B-field (>~1e+12 G) over a weak field as in low-mass X-ray binaries. The observed low X-ray luminosity (1e+33-1e+35 erg/s), probable wide orbit, and a slow stellar wind of this SyXB make quasi-spherical accretion in the subsonic settling regime a plausible model. Assuming a ~1e+13 G NS, this scheme can explain the ~5.4 h equilibrium rotation without employing the magnetar-like field (~1e+16 G) required in the disk accretion case. The time-scales of multiple irregular flares (~50 s) can also be attributed to the free-fall time from the Alfven shell for a ~1e+13 G field. A physical interpretation of SyXBs beyond the canonical binary classifications is discussed.

In this study, we conduct three-dimensional hydrodynamic simulations systematically to investigate the flow patterns behind the accretion shock waves that are commonly formed in the post-bounce phase of core-collapse supernovae. Adding small perturbations to spherically symmetric, steady, shocked accretion flows, we compute the subsequent evolutions to find what flow pattern emerges as a consequence of hydrodynamical instabilities such as convection and standing accretion shock instability (SASI) for different neutrino luminosities and mass accretion rates. Depending on these two controlling parameters, various flow patterns are indeed realized. We classify them into three basic patterns and two intermediate ones; the former includes sloshing motion (SL), spiral motion (SP) and multiple buoyant bubble formation (BB); the latter consists of spiral motion with buoyant-bubble formation (SPB) and spiral motion with pulsationally changing rotational velocities (SPP). Although the post-shock flow is highly chaotic, there is a clear trend in the pattern realization. The sloshing and spiral motions tend to be dominant for high accretion rates and low neutrino luminosities, and multiple buoyant bubbles prevail for low accretion rates and high neutrino luminosities. It is interesting that the dominant pattern is not always identical between the semi-nonlinear and nonlinear phases near the critical luminosity; the intermediate cases are realized in the latter case. Running several simulations with different random perturbations, we confirm that the realization of flow pattern is robust in most cases.

We investigated the impacts of magnetorotational instability (MRI) on the dynamics of weakly magnetized, rapidly rotating core-collapse by conducting high resolution MHD simulations in axisymmetry with simplified neutrino transfer. We found that an initially sub-magnetar class magnetic field is drastically amplified by MRI and substantially affects the dynamics thereafter. Although the magnetic pressure is not strong enough to eject matter, the amplified magnetic field efficiently transfers angular momentum from higher to lower latitudes, which causes the expansion of the heating region at low latitudes due to the extra centrifugal force. This then enhance the efficiency of neutrino heating and eventually leads to neutrino-driven explosion. This is a new scenario of core-collapse supernovae that has never been demonstrated by numerical simulations so far.

We have fabricated multilayer readout wiring Transition Edge Sensors (TES), which enable us to realize both large effective area and high-energy resolution for future X-ray astrophysical missions, such as DIOS. By sandwiching a SiO2 insulation layer between Al superconducting signal and return lines, self/mutual inductances and self fielding of bias leads are expected to be reduced. We fabricated 4x4 and 20x20 TES array on the multilayer wiring and tested their performance. Under the low temperature condition, several pixels in the TES array showed sharp superconducting transitions at around ~300 mK. We also succeeded in detecting X-ray signals from the 4x4 TES, contrary to the previous results of 20x20 TES. We further investigated the reasons for the differences between the 4x4 TES and the 20x20 TES, and present future plans for improving the multilayer TES array fabrication.

We investigate the chemical enrichment of r-process elements in the early evolutionary stages of the Milky Way halo within the framework of hierarchical galaxy formation using a semi-analytic merger tree. In this paper, we focus on heavy r-process elements, Ba and Eu, of extremely metal-poor (EMP) stars and give constraints on their astronomical sites. Our models take into account changes of the surface abundances of EMP stars by the accretion of interstellar matter (ISM). We also consider metal-enrichment of intergalactic medium (IGM) by galactic winds and the resultant pre-enrichment of proto-galaxies. The trend and scatter of the observed r-process abundances are well reproduced by our hierarchical model with $\sim 10\%$ of core-collapse supernovae in low-mass end ($\sim 10M_{\odot}$) as a dominant r-process source and the star formation efficiency of $\sim 10^{-10} \hbox{yr}^{-1}$. For neutron star mergers as an r-process source, their coalescence timescale has to be $ \sim 10^7$yrs, and the event rates $\sim 100$ times larger than currently observed in the Galaxy. We find that the accretion of ISM is a dominant source of r-process elements for stars with [Ba/H] < -3.5. In this model, a majority of stars at [Fe/H] < -3 are formed without r-process elements but their surfaces are polluted by the ISM accretion. The pre-enrichment affects $\sim 4\%$ of proto-galaxies, and yet, is surpassed by the ISM accretion in the surface of EMP stars.

We have built a code to obtain the exact solutions of Riemann problems in ideal magnetohydrodynamics (MHD) for an arbitrary initial condition. The code can handle not only regular waves but also switch-on/off rarefactions and all types of non-regular shocks: intermediate shocks and switch-on/off shocks. Furthermore, the initial conditions with vanishing normal or transverse magnetic fields can be handled although the code is partly based on the algorithm proposed by Torrilhon (2002) (Torrilhon 2002 Exact Solver and Uniqueness Condition for Riemann problems of Ideal Magnetohydrodynamics. Research report 2002-06, Seminar for Applied Mathematics, ETH, Zurich), which cannot deal with all types of non-regular waves nor the initial conditions without normal or transverse magnetic fields. Our solver can find all the solutions for a given Riemann problem and hence, as demonstrated in this paper, one can investigate the structure of the solution space in detail. Therefore the solver is a powerful instrument to solve the outstanding problem of the existence and uniqueness of solutions of MHD Riemann problems. Moreover, the solver may be applied to numerical MHD schemes like the Godunov scheme in the future.

We study the enrichment histories for nine elements, C, four alpha-elements of Mg, Si, Ca, and Ti, Sc, and three iron-peak elements of Co, Ni, and Zn, by using a large number of stellar data, collected by the Stellar Abundances for Galactic Archaeology (SAGA) database. We find statistically significant changes, or breaks, of the mean abundance ratios to iron at three metallicities of [Fe/H]-1.8, -2.2, and -3.3. Across the first one, the mean abundance ratios decrease with the metallicity by similar extents for all the elements with the sufficient data. Across the latter two, downward trends with the metallicity are also detected but for limited elements, C, Co, Zn, and possibly Sc, and for two of Co and Zn, respectively. The breaks define four stellar populations with the different abundance patters which are dominant in each metallicity range divided by the breaks, Pop IIa, IIb, IIc, and IId in order of increasing metallicity. We also explore their spatial distributions with the spectroscopic distances to demonstrate that Pops IIa and IIb spread over the Galactic halo while Pops IIc and IId are observed near the Galactic plane. In particular, Pop IIc stars emerge around [Fe/H] -2.6 and coexist with Pop IIb stars, segregated by the spatial distributions. Our results reveal two distinct modes of star formation during the early stages of Galaxy formation, which are associated with the variations of IMF and the spatial distribution of remnant low-mass stars. For the two lower-metallicity populations, the enhancements of Zn and Co indicate a high-mass and top-heavy IMF together with the statistics on the carbon-enhanced stars. We discuss the relevance to the kinematically resolved structures of the Galactic halo and the possible sites of these populations within the framework of hierarchical structure formation scenario.

Hadronic matter undergoes a deconfinement transition to quark matter at high temperature and/or high density. It would be realized in collapsing cores of massive stars. In the framework of MIT bag model, the ambiguities of the interaction are encapsulated in the bag constant. Some progenitor stars that invoke the core collapses explode as supernovae, and other ones become black holes. The fates of core collapses are investigated for various cases. Equations of state including the hadron-quark phase transition are constructed for the cases of the bag constant B=90, 150 and 250 MeV fm^{-3}. To describe the mixed phase, the Gibbs condition is used. Adopting the equations of state with different bag constants, the core collapse simulations are performed for the progenitor models with 15 and 40Msolar. If the bag constant is small as B=90 MeV fm^{-3}, an interval between the bounce and black hole formation is shortened drastically for the model with 40Msolar and the second bounce revives the shock wave leading to explosion for the model with 15Msolar.

The neutron-star Low-Mass X-ray Binary Aquila X-1 was observed by Suzaku for seven times, from 2007 September 28 to October 30. The observations successfully traced an outburst decay in which the source luminosity decreased almost monotonically from $\sim 10^{37}$ erg s$^{-1}$ to $\sim 10^{34}$ erg s$^{-1}$, by $\sim 3$ orders of magnitude. To investigate luminosity-dependent changes in the accretion geometry, five of the seven data sets with a typical exposure of $\sim 18$ ks each were analyzed; the other two were utilized in a previous work \citep{Sakurai2012}. The source was detected up to 100 keV in the 2nd to the 4th observations, to 40 keV in the 5th, and to 10 keV on the last two occasions. All spectra were reproduced successfully by Comptonized blackbody model with relatively high ($\gtrsim 2.0$) optical depths, plus an additional softer optically-thick component. The faintest three spectra were reproduced alternatively by a single Comptonized blackbody model with a relatively low ($\lesssim 0.8$) optical depth. The estimated radius of the blackbody emission, including seed photons for the Comptonization, was $10 \pm 2$ km at a 0.8--100 keV luminosity of $2.4\times 10^{36}$ erg s$^{-1}$ (the 2nd to the 4th observations). In contrast, it decreased to $7 \pm 1$ km and further to $3 \pm 1$ km, at a luminosity of $(4.8-5.2)\times 10^{35}$ erg s$^{-1}$ (the 5th observation) and $\sim 2\times 10^{34}$ erg s$^{-1}$ (the 6th and 7th), respectively, regardless of the above model ambiguity. This can be taken as evidence for the emergence of a weak magnetosphere of the neutron star.

We explore photospheric emissions from stratified two-component jets, wherein a highly relativistic spine outflow is surrounded by a wider and less relativistic sheath outflow. Thermal photons are injected in regions of high optical depth and propagated until they escape at the photosphere. Due to the presence of shear in velocity (Lorentz factor) at the boundary of the spine and sheath region, a fraction of the injected photons are accelerated via a Fermi-like acceleration mechanism such that a high energy power-law tail is formed in the resultant spectrum. We show, in particular, that if a velocity shear with a considerable variance in the bulk Lorentz factor is present, the high energy part of observed Gamma-ray Bursts (GRBs) photon spectrum can be explained by this photon acceleration mechanism. We also show that the accelerated photons may also account for the origin of the extra hard power-law component above the bump of the thermal-like peak seen in some peculiar bursts (e.g., GRB 090510, 090902B, 090926A). It is demonstrated that time-integrated spectra can also reproduce the low energy spectrum of GRBs consistently due to a multi-temperature effect when time evolution of the outflow is considered. Finally, we show that the empirical Ep-Lp relation can be explained by differences in the outflow properties of individual sources.

We perform numerical experiments to investigate the influence of inelastic neutrino reactions with light nuclei on the standing accretion shock instability (SASI). The time evolution of shock waves is calculated with a simple light-bulb approximation for the neutrino transport and a multi-nuclei equation of state. The neutrino absorptions and inelastic interactions with deuterons, tritons, helions and alpha particles are taken into account in the hydrodynamical simulations in addition to the ordinary charged-current interactions with nucleons. Axial symmetry is assumed but no equatorial symmetry is imposed. We show that the heating rates of deuterons reach as high as 10% of those of nucleons around the bottom of the gain region. On the other hand, alpha particles are heated near the shock wave, which is important when the shock wave expands and the density and temperature of matter become low. It is also found that the models with heating by light nuclei have different evolutions from those without it in the non-linear phase of SASI. This results is because matter in the gain region has a varying density and temperature and there appear sub-regions that are locally rich in deuterons and alpha particles. Although the light nuclei are never dominant heating sources and they work favorably for shock revival in some cases and unfavorably in other cases, they are non-negligible and warrant further investigation.

We construct new equations of state for baryons at sub-nuclear densities for the use in core-collapse simulations of massive stars. The abundance of various nuclei is obtained together with thermodynamic quantities. A model free energy is constructed, based on the relativistic mean field theory for nucleons and the mass formula for nuclei with the proton number up to ~ 1000. The formulation is an extension of the previous model, in which we adopted the liquid drop model to all nuclei under the nuclear statistical equilibrium. We reformulate the new liquid drop model so that the temperature dependences of bulk energies could be taken into account. Furthermore, we extend the region in the nuclear chart, in which shell affects are included, by using theoretical mass data in addition to experimental ones. We also adopt a quantum theoretical mass evaluation of light nuclei, which incorporates the Pauli- and self-energy shifts that are not included in the ordinary liquid drop model. The pasta phases for heavy nuclei are taken into account in the same way as in the previous model. We find that the abundances of heavy nuclei are modified by the shell effects of nuclei and temperature dependence of bulk energies. These changes may have an important effect on the rates of electron captures and coherent neutrino scatterings on nuclei in supernova cores. The abundances of light nuclei are also modified by the new mass evaluation, which may affect the heating and cooling rates of supernova cores and shocked envelopes.

We report the first observation of 3$-$30 MeV prolonged gamma-ray emission that was abruptly terminated by lightning. The gamma-ray detection was made during winter thunderstorms on December 30, 2010 by the Gamma-Ray Observation of Winter THunderclouds (GROWTH) experiment carried out in a coastal area along the Sea of Japan. The gamma-ray flux lasted for less than 3 min, continuously hardening closer to the lightning occurrence. The hardening at energies of 3$-$10 MeV energies was most prominent. The gamma-ray flux abruptly ceased less than 800 ms before the lightning flash that occurred over 5 km away from the experimental site. In addition, we observed a clear difference in the duration of the 3$-$10 MeV gamma rays and those $>$10 MeV, suggesting that the area of $>$10 MeV gamma-ray emission is considerably smaller than that of the lower-energy gamma rays. This work may give a manifestation that a local region emitting prolonged gamma rays connects with a distant region to initiate lightning.

We present the discovery of a reflection-dominated active galactic nucleus (AGN) in the early-type radio-quiet galaxy ESO 565--G019 with Suzaku and Swift/BAT. The source X-ray spectrum below 10 keV is characteristic of other Compton-thick (CT) AGN, clearly showing an inverted continuum and prodigious fluorescence iron emission above ~3 keV. A Compton shoulder to the neutral Fe Kalpha line also appears to be present. There is evidence for long-term hard X-ray flux variability which we associate with changes in the intrinsic AGN power-law. The increasing sensitivity of ongoing and new hard X-ray surveys means that more such reflection-dominated AGN ought to be uncovered in the near future. ESO 565--G019 is hosted in an early-type galaxy whose morphology has been variously classified as either type E or type S0. Only about 20 bona fide CT-AGN have been identified in the local universe so far, and all exist in host galaxies with late Hubble types (S0 or later). CT columns of nuclear obscuring gas are uncommon in early-type galaxies in the local universe, so confirmation of the exact morphological class of ESO 565--G019 is important. Infrared photometry also shows the presence of large quantities of cool dust in the host, indicative of significant ongoing star-formation. ESO 565--G019 may be the first identified local example of minor-merger driven CT-AGN growth in an early-type host, or may be the result of interaction with its neighboring galaxy ESO 565--G018 in a wide pair.

We perform some experimental simulations in spherical symmetry and axisymmetry to understand the post-shock-revival evolution of core-collapse supernovae. Assuming that the stalled shock wave is relaunched by neutrino heating and employing the so-called light bulb approximation, we induce shock revival by raising the neutrino luminosity by hand up to the critical value, which is also de- termined by dynamical simulations. A 15M_{sun} progenitor model is employed. We incorporate nuclear network calculations with a consistent equation of state in the simulations to account for the energy release by nuclear reactions and their feedback to hydrodynamics. Varying the shock-relaunch time rather arbitrarily, we investigate the ensuing long-term evolutions systematically, paying particular attention to the explosion energy and nucleosynthetic yields as a function of this relaunch time, or equivalently the accretion rate at shock revival. We study in detail how the diagnostic explosion energy approaches the asymptotic value and which physical processes contribute to the explosion energy in what proportions as well as their dependence on the relaunch time and the dimension of dynamics. We find that the contribution of nuclear reactions to the explosion energy is comparable to or greater than that of neutrino heating. In particular, recombinations are dominant over burnings in the contributions of nuclear reactions. Interestingly 1D models studied in this paper cannot produce the appropriate explosion energy and nickel mass simultaneously, overproducing nickels, whereas this problem is resolved in 2D models if the shock is relaunched at 300-400ms after bounce.

Rapid spectral changes in the hard X-ray on a time scale down to ~0.1 s are studied by applying "shot analysis" technique to the Suzaku observations of the black hole binary Cygnus X-1, performed on 2008 April 18 during the low/hard state. We successfully obtained the shot profiles covering 10--200 keV with the Suzaku HXD-PIN and HXD-GSO detector. It is notable that the 100-200 keV shot profile is acquired for the first time owing to the HXD-GSO detector. The intensity changes in a time-symmetric way, though the hardness does in a time-asymmetric way. When the shot-phase-resolved spectra are quantified with the Compton model, the Compton y-parameter and the electron temperature are found to decrease gradually through the rising phase of the shot, while the optical depth appears to increase. All the parameters return to their time-averaged values immediately within 0.1 s past the shot peak. We have not only confirmed this feature previously found in energies below ~60 keV, but also found that the spectral change is more prominent in energies above ~100 keV, implying the existence of some instant mechanism for direct entropy production. We discuss possible interpretations on the rapid spectral changes in the hard X-ray band.

Unified X-ray spectral and timing studies of Cygnus X-1 in the low/hard and hard intermediate state were conducted in a model-independent manner, using broadband Suzaku data acquired on 25 occasions from 2005 to 2009 with a total exposure of ~ 450 ks. The unabsorbed 0.1--500 keV source luminosity changed over 0.8--2.8% of the Eddington limit for 14.8 solar masses. Variations on short (1--2 seconds) and long (days to months) time scales require at least three separate components: a constant component localized below ~2 keV, a broad soft one dominating in the 2--10 keV range, and a hard one mostly seen in 10--300 keV range. In view of the truncated disk/hot inner flow picture, these are respectively interpreted as emission from the truncated cool disk, a soft Compton component, and a hard Compton component. Long-term spectral evolution can be produced by the constant disk increasing in temperature and luminosity as the truncation radius decreases. The soft Compton component likewise increases, but the hard Compton does not, so that the spectrum in the hard intermediate state is dominated by the soft Compton component; on the other hand, the hard Compton component dominates the spectrum in the dim low/hard state, probably associated with a variable soft emission providing seed photons for the Comptonization.

We present observations of a transient He-like Fe K alpha absorption line in Suzaku observations of the black hole binary Cygnus X-1 on 2011 October 5 near superior conjunction during the high/soft state, which enable us to map the full evolution from the start and the end of the episodic accretion phenomena or dips for the first time. We model the X-ray spectra during the event and trace their evolution. The absorption line is rather weak in the first half of the observation, but instantly deepens for ~10 ks, and weakens thereafter. The overall change in equivalent width is a factor of ~3, peaking at an orbital phase of ~0.08. This is evidence that the companion stellar wind feeding the black hole is clumpy. By analyzing the line with a Voigt profile, it is found to be consistent with a slightly redshifted Fe XXV transition, or possibly a mixture of several species less ionized than Fe XXV. The data may be explained by a clump located at a distance of ~10^(10-12) cm with a density of ~10^((-13)-(-11)) g cm^-3, which accretes onto and/or transits the line-of-sight to the black hole, causing an instant decrease in the observed degree of the ionization and/or an increase in density of the accreting matter. Continued monitoring for individual events with future X-ray calorimeter missions such as ASTRO-H and AXSIO will allow us to map out the accretion environment in detail and how it changes between the various accretion states.

The Crab Nebula is one of the brightest and most stable sources in the X-ray sky. Year-scale flux variation from the object was recently revealed in the hard X-ray band by four satellites. This marked the first detection of year-scale variability from pulsar wind nebulae in the hard X-ray band. The Crab Nebula has been observed at least once a year for calibration purposes with the Suzaku Hard X-ray Detector (HXD) since its launch in 2005. In order to investigate possible spectral changes as well as flux variation, the archival data of the HXD were analyzed. The flux variation reported by other instruments was confirmed in the 25 -- 100 keV band by the HXD in a few percent level, but flux above 100 keV did not follow the trend in variation below 100 keV. The hardness ratios produced utilizing the PIN and GSO sensors installed in the HXD exhibit significant scattering, thereby indicating spectral variations in the hard X-ray. The spectral changes are quantified by spectral fitting with a broken power-law model. The difference between the two photon indexes of the broken power-law model in harder and softer energy bands is in the range of < 2.54. Taking into account flux variation of 6.3% and spectral variation time-scale of a few days, multi components of the broken power-law-shaped synchrotron emission with different cooling times are suggested.

We compare three different statistical models for the equation of state (EOS) of stellar matter at subnuclear densities and temperatures (0.5-10 MeV) expected to occur during the collapse of massive stars and supernova explosions. The models introduce the distributions of various nuclear species in nuclear statistical equilibrium, but use somewhat different nuclear physics inputs. It is demonstrated that the basic thermodynamical quantities of stellar matter under these conditions are similar, except in the region of high densities and low temperatures. We demonstrate that mass and isotopic distributions have considerable differences related to the different assumptions of the models on properties of nuclei at these stellar conditions. Overall, the three models give similar trends, but the details reflect the uncertainties related to the modeling of medium effects, such as the temperature and density dependence of surface and bulk energies of heavy nuclei, and the nuclear shell structure effects. We discuss importance of new physics inputs for astrophysical calculations from experimental data obtained in intermediate energy heavy-ion collisions, in particular, the similarities of the conditions reached during supernova explosions and multifragmentation reactions.

We performed the first numerical simulations of magnetorotational instability from a sub-magnetar-class seed magnetic field in core collapse supernovae. As a result of axisymmetric ideal MHD simulations, we found that the magnetic field is greatly amplified to magnetar-class strength. In saturation phase, a substantial part of the core is dominated by turbulence, and the magnetic field possesses dominant large scale components, comparable to the size of the proto-neutron star. A pattern of coherent chanel flows, which generally appears during exponential growth phase in previous local simulations, is not observed in our global simulations. While the approximate convergence in the exponential growth rate is attained by increasing spatial resolution, that of the saturation magnetic field is not achieved due to still large numerical diffusion. Although the effect of magnetic field on the dynamics is found to be mild, a simulation with a high-enough resolution might result in a larger impact.

The stellar initial mass function (IMF) plays a crucial role in determining the number of surviving stars in galaxies, the chemical composition of the interstellar medium, and the distribution of light in galaxies. A key unsolved question is whether the IMF is universal in time and space. Here we use state-of-the-art results of stellar evolution to show that the IMF of our Galaxy made a transition from an IMF dominated by massive stars to the present-day IMF at an early phase of the Galaxy formation. Updated results from stellar evolution in a wide range of metallicities have been implemented in a binary population synthesis code, and compared with the observations of carbon-enhanced metal-poor (CEMP) stars in our Galaxy. We find that applying the present-day IMF to Galactic halo stars causes serious contradictions with four observable quantities connected with the evolution of AGB stars. Furthermore, a comparison between our calculations and the observations of CEMP stars may help us to constrain the transition metallicity for the IMF which we tentatively set at [Fe/H] = -2. A novelty of the current study is the inclusion of mass loss suppression in intermediate-mass AGB stars at low-metallicity. This significantly reduces the overproduction of nitrogen-enhanced stars that was a major problem in using the high-mass star dominated IMF in previous studies. Our results also demonstrate that the use of the present day IMF for all time in chemical evolution models results in the overproduction of Type I.5 supernovae. More data on stellar abundances will help to understand how the IMF has changed and what caused such a transition.

We studied roles of a turbulent resistivity in the core-collapse of a strongly magnetized massive star, carrying out 2D-resistive-MHD simulations. The three cases with different initial strengths of magnetic field and rotation are investigated; 1. strongly magnetized rotating core; 2.moderately magnetized rotating core; 3. very strongly magnetized non-rotating core. In each case, both an ideal-MHD model and resistive-MHD models are computed. As a result of computations, each model shows a matter eruption helped by a magnetic acceleration (and also by a centrifugal acceleration in the rotating cases). We found that a resistivity attenuates the explosion in case~1 and 2, while it enhances the explosion in case~3. We also found that in the rotating cases, main mechanisms for the amplification of a magnetic field in the post-bounce phase are an outward advection of magnetic field and a winding of poloidal magnetic field-lines by differential rotation, which are somewhat dampened down with the presence of a resistivity. Although the magnetorotational instability seems to occur in the rotating models, it will play only a minor role in a magnetic field amplification. Another impact of resistivity is that on the aspect ratio. In the rotating cases, a large aspect ratio of the ejected matters, $> 2.5$, attained in a ideal-MHD model is reduced to some extent in a resistive model. These results indicate that a resistivity possibly plays an important role in the dynamics of strongly magnetized supernovae.

We develop a new semi-dynamical method to study shock revival by neutrino heating in core- collapse supernovae. Our new approach is an extension of the previous studies that employ spherically symmetric, steady, shocked accretion flows together with the light bulb approximation. The latter has been widely used in the supernova community for the phenomenological investigation of the criteria for successful supernova explosions. In the present approach, on the other hand, we get rid of the steady-state condition and take into account shock wave motions instead. We have in mind the scenario that not the critical luminosity but the critical fluctuation generated by hydrodynamical instabilities such as SASI and neutrino-driven convection in the post-shock region determines the onset of shock revival. After confirming that the new approach indeed captures the dynamics of revived shock wave qualitatively, we then apply the method to various initial conditions and find that there is a critical fluctuation for shock revival, which can be well fit by the following formula: f_crit ~ 0.8 * (M_in/1.4M_sun) * {1- (rsh/10^8cm)}, in which fcrit denotes the critical pressure fluctuation normalized by the unperturbed post-shock value. Min and rsh stand for the mass of the central compact object and the shock radius, respectively. The critical fluctuation decreases with the shock radius, whereas it increases with the mass of the central object. We discuss the possible implications of our results for 3D effects on shock revival, which is currently controversial in the supernova community.

We present a new series of supernova neutrino light curves and spectra calculated by numerical simulations for a variety of progenitor stellar masses (13-50Msolar) and metallicities (Z = 0.02 and 0.004), which would be useful for a broad range of supernova neutrino studies, e.g., simulations of future neutrino burst detection by underground detectors, or theoretical predictions for the relic supernova neutrino background. To follow the evolution from the onset of collapse to 20 s after the core bounce, we combine the results of neutrino-radiation hydrodynamic simulations for the early phase and quasi-static evolutionary calculations of neutrino diffusion for the late phase, with different values of shock revival time as a parameter that should depend on the still unknown explosion mechanism. We here describe the calculation methods and basic results including the dependence on progenitor models and the shock revival time. The neutrino data are publicly available electronically.

We have built a code to numerically solve the Riemann problem in ideal magnetohydrodynamics (MHD) for an arbitrary initial condition to investigate a variety of solutions more thoroughly. The code can handle not only regular solutions, in which no intermediate shocks are involved, but also all types of non-regular solutions if any. As a first application, we explored the neighborhood of the initial condition that was first picked up by Brio & Wu (1988) and has been frequently employed in the literature as a standard problem to validate numerical codes. Contrary to the conventional wisdom that there will always be a regular solution, we found an initial condition, for which there is no regular solution but a non-regular one. The latter solution has only regular solutions in its neighborhood and actually sits on the boundary of regular solutions. This implies that the regular solutions are not sufficient to solve the ideal MHD Riemann problem and suggests that at least some types of non-regular solutions are physical. We also demonstrate that the non-regular solutions are not unique. In fact, we found for the Brio & Wu initial condition that there are uncountably many non-regular solutions. This poses an intriguing question: why a particular non-regular solution is always obtained in numerical simulations? This has important ramifications to the discussion of which intermediate shocks are really admissible.

The joint JAXA/NASA ASTRO-H mission is the sixth in a series of highly successful X-ray missions initiated by the Institute of Space and Astronautical Science (ISAS). ASTRO-H will investigate the physics of the high-energy universe via a suite of four instruments, covering a very wide energy range, from 0.3 keV to 600 keV. These instruments include a high-resolution, high-throughput spectrometer sensitive over 0.3-2 keV with high spectral resolution of Delta E < 7 eV, enabled by a micro-calorimeter array located in the focal plane of thin-foil X-ray optics; hard X-ray imaging spectrometers covering 5-80 keV, located in the focal plane of multilayer-coated, focusing hard X-ray mirrors; a wide-field imaging spectrometer sensitive over 0.4-12 keV, with an X-ray CCD camera in the focal plane of a soft X-ray telescope; and a non-focusing Compton-camera type soft gamma-ray detector, sensitive in the 40-600 keV band. The simultaneous broad bandpass, coupled with high spectral resolution, will enable the pursuit of a wide variety of important science themes.

We investigate the nucleosynthesis in a massive star of 70 M_solar with solar metallicity in the main sequence stage. The helium core mass after hydrogen burning corresponds to 32 M_solar. Nucleosynthesis calculations have been performed during the stellar evolution and the jetlike supernova explosion of a collapsar model, where the weak s-, p-, and r-processes are taken into account. We confirm that s-elements of 60 < A < 90 are highly overproduced relative to the solar abundances in the hydrostatic nucleosynthesis. During oxygen burning, p-elements of A > 90 are produced via photodisintegrations of seed s-elements. However, the produced p-elements are disintegrated in later stages except for ^{180}Ta. In the explosive nucleosynthesis, elements of 90 < A < 160 are significantly overproduced relative to the solar values owing to the r-process. Only heavy p-elements (N > 50) are overproduced via the p-process. Compared with the previous study of r-process nucleosynthesis calculations in the collapsar model of 40 M_solar by Fujimoto et al. 2007, 2008, our jet model cannot contribute to the third peak of the solar r-elements and intermediate p-elements. Averaging the overproduction factors over the progenitor masses with the use of Salpeter's IMF, we suggest that the 70 M_solar star could contribute to the solar weak s-elements of 60 < A < 90 and neutron-rich elements of 90 < A < 160. We confirm the primary synthesis of light p-elements in the ejected matter of high peak temperature. The ejected matter has [Sr/Eu] \sim -0.4, which is different from that of a typical r-process-enriched star CS22892-052 ([Sr/Eu] \sim -1). We find that Sr-Y-Zr isotopes are primarily synthesized in the explosive nucleosynthesis in a similar process of the primary production of light p-elements, which has been considered as one of the sites of a lighter element primary process (LEPP).

We present high-resolution X-ray spectra of cloud-shock interaction regions in the eastern and northern rims of the Galactic supernova remnant Puppis A, using the Reflection Grating Spectrometer onboard the XMM-Newton satellite. A number of emission lines including K alpha triplets of He-like N, O, and Ne are clearly resolved for the first time. Intensity ratios of forbidden to resonance lines in the triplets are found to be higher than predictions by thermal emission models having plausible plasma parameters. The anomalous line ratios cannot be reproduced by effects of resonance scattering, recombination, or inner-shell ionization processes, but could be explained by charge-exchange emission that should arise at interfaces between the cold/warm clouds and the hot plasma. Our observations thus provide observational support for charge-exchange X-ray emission in supernova remnants.

To constrain the origin of the soft X-ray excess phenomenon seen in many active galactic nuclei, the intensity-correlated spectral analysis, developed by Noda et al. (2011b) for Markarian 509, was applied to wide-band (0.5-45 keV) Suzaku data of five representative objects with relatively weak reflection signature. They are the typical bare-nucleus type 1 Seyfert Fairall 9, the bright and typical type 1.5 Seyfert MCG-2-58-22, 3C382 which is one of the X-ray brightest broad line radio galaxies, the typical Seyfert-like radio loud quasar 4C+74.26, and the X-ray brightest radio quiet quasar MR2251-178. In all of them, soft X-ray intensities in energies below 3 keV were tightly correlated with that in 3-10 keV, but with significant positive offsets. These offsets, when calculated in finer energy bands, define a stable soft component in 0.5-3 keV. In each object, this component successfully explained the soft excess above a power-law fit. These components were interpreted in several alternative ways, including a thermal Comptonization component which is independent of the dominant power-law emission. This interpretation, considered physically most reasonable, is discussed from a viewpoint of Multi-Zone Comptonization, which was proposed for the black hole binary Cygnus X-1 (Makishima et al. 2008).

We present the result of an indirect search for high energy neutrinos from WIMP annihilation in the Sun using upward-going muon (upmu) events at Super-Kamiokande. Datasets from SKI-SKIII (3109.6 days) were used for the analysis. We looked for an excess of neutrino signal from the Sun as compared with the expected atmospheric neutrino background in three upmu categories: stopping, non-showering, and showering. No significant excess was observed. The 90% C.L. upper limits of upward-going muon flux induced by WIMPs of 100 GeV/c$^2$ were 6.4$\times10^{-15}$ cm$^{-2}$ sec$^{-1}$ and 4.0$\times10^{-15}$ cm$^{-2}$ sec$^{-1}$ for the soft and hard annihilation channels, respectively. These limits correspond to upper limits of 4.5$\times10^{-39}$ cm$^{-2}$ and 2.7$\times10^{-40}$ cm$^{-2}$ for spin-dependent WIMP-nucleon scattering cross sections in the soft and hard annihilation channels, respectively.

This is a status report on our endeavor to reveal the mechanism of core-collapse supernovae (CCSNe) by large-scale numerical simulations. Multi-dimensionality of the supernova engine, general relativistic magnetohydrodynamics, energy and lepton number transport by neutrinos emitted from the forming neutron star as well as nuclear interactions there, are all believed to play crucial roles in repelling infalling matter and producing energetic explosions. These ingredients are nonlinearly coupled with one another in the dynamics of core-collapse, bounce, and shock expansion. Serious quantitative studies of CCSNe hence make extensive numerical computations mandatory. Since neutrinos are neither in thermal nor in chemical equilibrium in general, their distributions in the phase space should be computed. This is a six dimensional (6D) neutrino transport problem and quite a challenge even for those with an access to the most advanced numerical resources such as the "K computer". To tackle this problem, we have embarked on multi-front efforts. In particular we report in this paper our recent progresses in the treatments of multi-dimensional (multi-D) radiation-hydrodynamics. We are currently proceeding on two different paths to the ultimate goal; in one approach we employ an approximate but highly efficient scheme for neutrino transport and treat 3D hydrodynamics and/or general relativity rigorously; some neutrino-driven explosions will be presented and comparisons will be made between 2D and 3D models quantitatively; in the second approach, on the other hand, exact but so far Newtonian Boltzmann equations are solved in two and three spatial dimensions; we will show some demonstrative test simulations. We will also address the perspectives of exa-scale computations on the next generation supercomputers.

We report on a discovery of a diffuse nebula containing a pointlike source in the southern blowout region of the Cygnus Loop supernova remnant, based on Suzaku and XMM-Newton observations. The X-ray spectra from the nebula and the pointlike source are well represented by an absorbed power-law model with photon indices of 2.2+/-0.1 and 1.6+/-0.2, respectively. The photon indices as well as the flux ratio of F_nebula/F_pointlike ~ 4 lead us to propose that the system is a pulsar wind nebula, although pulsations have not yet been detected. If we attribute its origin to the Cygnus Loop supernova, then the 0.5-8 keV luminosity of the nebula is computed to be 2.1e31 (d/540pc)^2 ergs/s, where d is the distance to the Loop. This implies a spin-down loss-energy E_dot ~ 2.6e35 (d/540pc)^2 ergs/s. The location of the neutron star candidate, ~2 degrees away from the geometric center of the Loop, implies a high transverse velocity of ~1850 (d/540pc)(t/10kyr)^{-1} km/s, assuming the currently accepted age of the Cygnus Loop.

We present analysis of 4U 1626-67, a 7.7 s pulsar in a low-mass X-ray binary system, observed with the hard X-ray detector of the Japanese X-ray satellite Suzaku in March 2006 for a net exposure of \sim88 ks. The source was detected at an average 10-60 keV flux of \sim4 x10^-10 erg cm^-2 s^-1. The phase-averaged spectrum is reproduced well by combining a negative and positive power-law times exponential cutoff (NPEX) model modified at \sim 37 keV by a cyclotron resonance scattering feature (CRSF). The phase-resolved analysis shows that the spectra at the bright phases are well fit by the NPEX with CRSF model. On the other hand, the spectrum in the dim phase lacks the NPEX high-energy cutoff component, and the CRSF can be reproduced by either an emission or an absorption profile. When fitting the dim phase spectrum with the NPEX plus Gaussian model, we find that the feature is better described in terms of an emission rather than an absorption profile. The statistical significance of this result, evaluated by means of an F-test, is between 2.91 x 10^-3 and 1.53 x 10^-5, taking into account the systematic errors in the background evaluation of HXD-PIN. We find that, the emission profile is more feasible than the absorption one for comparing the physical parameters in other phases. Therefore, we have possibly detected an emission line at the cyclotron resonance energy in the dim phase.

We report on the geometry of accretion disk and high energy coronae in the strong Comptonization state (the very high/steep power law/hard intermediate state) based on a Suzaku observation of the famous Galactic black hole GX 339-4. These data were taken just before the peak of the 2006-2007 outburst, and the average X-ray luminosity in the 0.7-200 keV band is estimated to be 2.9E38 erg/s for a distance of 8 kpc. We fit the spectrum with both simple (independent disk and corona) and sophisticated (energetically coupled disk and corona) models, but all fits imply that the underlying optically thick disk is truncated significantly before the innermost stable circular orbit around the black hole. We show this directly by a comparison with similarly broadband data from a disk dominated spectrum at almost the same luminosity observed by XMM-Newton and RXTE 3 days after the Suzaku observation. During the Suzaku observation, the QPO frequency changes from 4.3 Hz to 5.5 Hz, while the spectrum softens. The energetically coupled model gives a corresponding 5+/- 8 % decrease in derived inner radius of the disk. While this is not significant, it is consistent with the predicted change in QPO frequency from Lense-Thirring precession of the hot flow interior to the disk and/or a deformation mode of this flow, as a higher QPO frequency implies a smaller size scale for the corona. This is consistent with the truncated disk extending further inwards towards the black hole.

We investigate the possibility of the r-process during the magnetohydrohynamical explosion of supernova in a massive star of 13 solar mass with the effects of neutrinos induced. We adopt five kinds of initial models which include properties of rotation and the toroidal component of the magnetic field . The simulations which succeed the explosions are limitted to a concentrated magnetic field and strong differential rotation. Low $Y_{e}$ ejecta produce heavy elements and the third peak can be reprocuced. However, the second peak is low because $Y_{e}$ distribution as a function of radius is steep and ejecta corresponding to middle $Y_{e}$ is very few.

GUT monopoles captured by the Sun's gravitation are expected to catalyze proton decays via the Callan-Rubakov process. In this scenario, protons, which initially decay into pions, will ultimately produce \nu_{e}, \nu_{\mu} and \bar{\nu}_{\mu}. After undergoing neutrino oscillation, all neutrino species appear when they arrive at the Earth, and can be detected by a 50,000 metric ton water Cherenkov detector, Super-Kamiokande (SK). A search for low energy neutrinos in the electron total energy range from 19 to 55 MeV was carried out with SK and gives a monopole flux limit of F_M(\sigma_0/1 mb) < 6.3 \times 10^{-24} (\beta_M/10^{-3})^2 cm^{-2} s^{-1} sr^{-1} at 90% C.L., where \beta_M is the monopole velocity in units of the speed of light and \sigma_0 is the catalysis cross section at \beta_M=1. The obtained limit is more than eight orders of magnitude more stringent than the current best cosmic-ray supermassive monopole flux limit, F_M < 1 \times 10^{-15} cm^{-2} s^{-1} sr^{-1} for \beta_M < 10^{-3} and also two orders of magnitude lower than the result of the Kamiokande experiment, which used a similar detection method.

The neutron-star Low-Mass X-ray Binary Aquila X-1 was observed seven times in total with the Suzaku X-ray observatory from September 28 to October 30 in 2007, in the decaying phase of an outburst. In order to constrain the flux-dependent accretion geometry of this source over wider energy bands than employed in most of previous works, the present study utilized two out of the seven data sets. The 0.8-31 keV spectrum on September 28, taken with the XIS and HXD-PIN for an exposure of 13.8 ks, shows an absorbed 0.8-31 keV flux of $3.6\times 10^{-9}$ erg s$^{-1}$ cm$^{-2}$, together with typical characteristics of the soft state of this type of objects. The spectrum was successfully explained by an optically-thick disk emission plus a Comptonized blackbody component. Although these results are in general agreement with previous studies, the significance of a hard tail recently reported using the same data was inconclusive in our analysis. The spectrum acquired on October 9 for an exposure of 19.7 ks was detected over a 0.8-100 keV band with the XIS, HXD-PIN, and HXD-GSO, at an absorbed flux of $8.5\times 10^{-10}$ erg s$^{-1}$ cm$^{-2}$ (in 0.8-100 keV). It shows characteristics of the hard state, and was successfully explained by the same two continuum components but with rather different parameters including much stronger thermal Comptonization, of which the seed photon source was identified with blackbody emission from the neutron-star surface. As a result, the accretion flow in the hard state is inferred to take a form of an optically-thick and geometrically-thin disk down to a radius of $21\pm 4$ km from the neutron star, and then turn into an optically-thin nearly-spherical hot flow.

We develop a numerical code to calculate the neutrino transfer with multi-energy and multi-angle in three dimensions (3D) for the study of core-collapse supernovae. The numerical code solves the Boltzmann equations for neutrino distributions by the discrete-ordinate (S_n) method with a fully implicit differencing for time advance. The Boltzmann equations are formulated in the inertial frame with collision terms being evaluated to the zeroth order of v/c. A basic set of neutrino reactions for three neutrino species is implemented together with a realistic equation of state of dense matter. The pair process is included approximately in order to keep the system linear. We present numerical results for a set of test problems to demonstrate the ability of the code. The numerical treatments of advection and collision terms are validated first in the diffusion and free streaming limits. Then we compute steady neutrino distributions for a background extracted from a spherically symmetric, general relativistic simulation of 15Msun star and compare them with the results in the latter computation. We also demonstrate multi-D capabilities of the 3D code solving neutrino transfers for artificially deformed supernova cores in 2D and 3D. Formal solutions along neutrino paths are utilized as exact solutions. We plan to apply this code to the 3D neutrino-radiation hydrodynamics simulations of supernovae. This is the first article in a series of reports on the development.

We present the result of a systematic study of pileup phenomena seen in the X-ray Imaging Spectrometer, an X-ray CCD instrument, onboard the Suzaku observatory. Using a data set of observed sources in a wide range of brightness and spectral hardness, we characterized the pileup fraction, spectral hardening, and grade migration as a function of observed count rate in a frame per pixel. Using the pileup fraction as a measure of the degree of pileup, we found that the relative spectral hardening (the hardness ratio normalized to the intrinsic spectral hardness), branching ratio of split events, and that of detached events increase monotonically as the pileup fraction increases, despite the variety of brightness and hardness of the sample sources. We derived the pileup fraction as a function of radius used for event extraction. Upon practical considerations, we found that events outside of the radius with a pileup fraction of 1% or 3% are useful for spectral analysis. We present relevant figures, tables, and software for the convenience of users who wish to apply our method for their data reduction of piled-up sources.

We investigate the emergence of hyperons in black-hole-forming failed supernovae, which are caused by the dynamical collapse of nonrotating massive stars. We perform neutrino-radiation hydrodynamical simulations in general relativity adopting realistic hyperonic equation-of-state (EOS). Attractive and repulsive cases are examined for the potential of \Sigma\ hyperons. Since hyperons soften the EOS, they shorten the time interval from the bounce to black hole formation, which corresponds to the duration of neutrino emission. This effect is larger for the attractive case than the repulsive case because \Sigma\ hyperons appear more easily. In addition, we investigate the impacts of pions to find that they also promotes the recollapse towards the black hole formation.

The Type I Seyfert galaxy Markarian 509 was observed with Suzaku in 2010 November, for a gross time span of 2.2 days. Timing and spectral properties of the 0.5-45 keV X-rays, detected with the XIS and HXD, consistently revealed the presence of a soft spectral component that remained constant while the total X- ray intensity varied by \pm10%. This stable soft component, found in the 0.5-3.0 keV range, was interpreted as a result of thermal Comptonization in a corona with a temperature of ~ 0.5 keV and an optical depth of ~ 18. The time-avearged 0.5-45 keV Suzaku spectrum was reproduced successfully, as a combination of this thermal Comptonization component, a harder power-law of photon index ~ 1.8, moderate reflection, and an iron K-emission line. By analyzing four archival Suzaku datasets of the same object obtained in 2006, the thermal Comptonization component, which was stable during the 2.2 day pointing in 2010, was found to vary on time scales of a few weeks, independently of the power-law component. Implications of these results are discussed in terms of the "multi-zone Comptonization" view, obtained with Suzaku from the black hole binary Cygnus X-1.

We report on the X-ray spectral analysis of the black hole candidate XTE\ J1752--223 in the 2009--2010 outburst, utilizing data obtained with the MAXI/Gas Slit Camera (GSC), the Swift/XRT, and Suzaku, which work complementarily. As already reported by Nakahira et al. (2010) MAXI monitored the source continuously throughout the entire outburst for about eight months. All the MAXI/GSC energy spectra in the high/soft state lasting for 2 months are well represented by a multi-color disk plus power-law model. The innermost disk temperature changed from $\sim$0.7 keV to $\sim$0.4 keV and the disk flux decreased by an order of magnitude. Nevertheless, the innermost radius is constant at $\sim$41 $D_{3.5}(\cos{\it i})^{-1/2}$ km, where $D_{3.5}$ is the source distance in units of 3.5 kpc and $i$ the inclination. The multi-color disk parameters obtained with the MAXI/GSC are consistent with those with the Swift/XRT and Suzaku. The Suzaku data also suggests a possibility that the disk emission is slightly Comptonized, which could account for broad iron-K features reported previously. Assuming that the obtained innermost radius represents the innermost stable circular orbit for a non-rotating black hole, we estimate the mass of the black hole to be 5.51$\pm$0.28 $M_{\odot}$ $D_{3.5}(\cos{\it i})^{-1/2}$, where the correction for the stress-free inner boundary condition and color hardening factor of 1.7 are taken into account. If the inclination is less than 49$^{\circ}$ as suggested from the radio monitoring of transient jets and the soft-to-hard transition in 2010 April occurred at 1--4% of Eddignton luminosity, the fitting of the Suzaku spectra with a relativistic accretion-disk model derives constraints on the mass and the distance to be 3.1--55 $M_{\odot}$ and 2.3--22 {\rm kpc}, respectively. This confirms that the compact object in XTE J1752--223 is a black hole.

We calculate a new equation of state for baryons at sub-nuclear densities meant for the use in core-collapse simulations of massive stars. The abundance of various nuclei is obtained together with the thermodynamic quantities. The formulation is the NSE description and the liquid drop approximation of nuclei. The model free energy to minimize is calculated by relativistic mean field theory for nucleons and the mass formula for nuclei with the atomic number up to ~ 1000. We have also taken into account the pasta phase, thanks to which the transition to uniform nuclear matter in our EOS occurs in the conventional manner: nuclei are not dissociated to nucleons but survive right up to the transition to uniform nuclear matter. We find that the free energy and other thermodynamical quantities are not very different from those given in the Shen's EOS, one of the standard EOS's that adopt the single nucleus approximation. The average mass is systematically different, on the other hand, which may have an important ramification to the rates of electron captures and coherent neutrino scatterings on nuclei in supernova cores. It is also interesting that the root mean square of the mass number is not very different from the average mass number, since the former is important for the evaluation of coherent scattering rates on nuclei but has been unavailable so far. The EOS table is currently under construction, which will include the weak interaction rates.

Improvements of in-orbit calibration of GSO scintillators in the Hard X-ray Detector on board Suzaku are reported. To resolve an apparent change of the energy scale of GSO which appeared across the launch for unknown reasons, consistent and thorough re-analyses of both pre-launch and in-orbit data have been performed. With laboratory experiments using spare hardware, the pulse height offset, corresponding to zero energy input, was found to change by ~0.5 of the full analog voltage scale, depending on the power supply. Furthermore, by carefully calculating all the light outputs of secondaries from activation lines used in the in-orbit gain determination, their energy deposits in GSO were found to be effectively lower, by several percent, than their nominal energies. Taking both these effects into account, the in-orbit data agrees with the on-ground measurements within ~5%, without employing the artificial correction introduced in the previous work (Kokubun et al. 2007). With this knowledge, we updated the data processing, the response, and the auxiliary files of GSO, and reproduced the HXD-PIN and HXD-GSO spectra of the Crab Nebula over 12-300 keV by a broken powerlaw with a break energy of ~110 keV.

From 2005 to 2009, 25 observations of Cyg X-1 were performed with Suzaku, achieving a total exposure of 446 ks. In all observations, the source was found in the low/hard state, while the 1.5--12.0 keV count rate of the All-Sky Monitor onboard RXTE varied by a factor of ~ 3. In each observation, the 10--60 keV HXD-PIN spectrum and the 60--400 keV HXD-GSO spectrum were fitted successfully by a thermal Comptonization model plus reflection by a thick neutral material. As the soft X-ray intensity increased, the Compton y-parameter was found to decrease from 1.0 to 0.6, while the solid angle of reflection to increase by ~ 30%. Also conducted was timing analysis over a frequency range of 0.001--10 Hz. As the source became brighter in soft X-rays, the characteristic frequency of hard X-ray variation increased from 0.03 to 0.3 Hz, while the fractional hard X-ray variation integrated over 0.001--0.01 Hz decreased by a factor of ~ 5. The signals in the 60--200 keV band were generally found to vary on shorter time scales than those in the 10--60 keV band. These spectral and timing results can be consistently interpreted by presuming that increases in the mass accretion rate cause the Comptonizing hot corona to shrink, while the optically-thick disk to intrude deeper therein.

Focusing on hard X-ray variability, we reanalyzed Suzaku data of Type I Seyfert galaxy MCG-6-30-15 obtained in 2006. Intensity-sorted spectroscopy and a principal component analysis consistently revealed a very hard component that varies independently of the dominant power-law emission. Although the exact nature of this hard component is not yet identified, it can be modeled as a power-law with a photon index ~2 affected by a partial covering absorption, or as a thermal Comptonization emission with a relatively large optical depth. When this component is included in the fitting model, the time-averaged 2.5-55 keV spectrum of MCG-6-30-15 can be reproduced successfully by invoking a mildly broadened iron line with its emission region located at > 8 times the gravitational radii from the central black hole, and a moderate reflection with a covering fraction of ~3.4. This result implies that the solution of a highly spinning black hole in MCG-6-30-15, obtained by Miniutti et al. (2007, PASJ, 59, S315) using the same Suzaku data, is a model dependent result.

The Gamma-Ray Observation of Winter THunderclouds (GROWTH) experiment, consisting of two radiation-detection subsystems, has been operating since 2006 on the premises of Kashiwazaki-Kariwa nuclear power plant located at the coastal area of Japan Sea. By 2010 February, GROWTH detected 7 long-duration $\gamma$-rays emissions associated with winter thunderstorms. Of them, two events, obtained on 2007 December 13 and 2008 December 25, are reported.On both occasions, all inorganic scintillators (NaI, CsI, and BGO) of the two subsystems detected significant gamma-ray signals lasting for >1 minute. Neither of these two events were associated with any lightning. In both cases, the gamma-ray energy spectra extend to 10 MeV, suggesting that the detected gamma-rays are produced by relativistic electrons via bremsstrahlung. Assuming that the initial photon spectrum at the source is expressed by a power-law function,the observed photons can be interpreted as being radiated from a source located at a distance of 290-560 m for the 2007 event and 110-690 m for the 2008 one, both at 90% confidence level.Employing these photon spectra, the number of relativistic electrons is estimated as 10^9 - 10^{11}. The estimation generally agrees with those calculated based on the relativistic runaway electron avalanche model. A GROWTH photon spectrum, summed over 3 individual events including the present two events and another reported previously, has similar features including a cut-off energy, to an averaged spectrum of terrestrial gamma-ray flashes.

We investigate by two-dimensional axisymmetric relativistic hydrodynamical simulations (1) jet propagations through an envelope of a rapidly rotating and collapsing massive star, which is supposed to be a progenitor of long duration gamma ray bursts (GRBs), (2) breakouts and subsequent expansions into stellar winds and (3) accompanying photospheric emissions. We find that if the envelope rotates uniformly almost at the mass shedding limit, its outer part stops contracting eventually when the centrifugal force becomes large enough. Then another shock wave is formed, propagates outwards and breaks out of the envelope into the stellar wind. Which breaks out earlier, the jet or the centrifugal bounce-induced shock, depends on the timing of jet injection. If the shock breakout occurs earlier owing to a later injection, the jet propagation and subsequent photospheric emissions are affected substantially. We pay particular attention to observational consequences of the difference in the timing of jet injection. We calculate optical depths to find the location of photospheres, extracting densities and temperatures at appropriate retarded times from the hydrodynamical data. We show that the luminosity and observed temperature of the photospheric emissions are both much lower than those reported in previous studies. Although luminosities are still high enough for GRBs, the observed temperature are lower than the energy at the spectral peak expected by the Yonetoku-relation. This may imply that energy exchanges between photons and matter are terminated deeper inside or some non-thermal processes are operating to boost photon energies.

We explore the evolution of the emissions by accelerated electrons in shocked shells driven by jets in active galactic nuclei (AGNs). Focusing on powerful sources which host luminous quasars, we evaluated the broadband emission spectra by properly taking into account adiabatic and radiative cooling effects on the electron distribution. The synchrotron radiation and inverse Compton (IC) scattering of various photons that are mainly produced in the accretion disc and dusty torus are considered as radiation processes. We show that the resultant radiation is dominated by the IC emission for compact sources (< 10kpc), whereas the synchrotron radiation is more important for larger sources. We also compare the shell emissions with those expected from the lobe under the assumption that a fractions of the energy deposited in the shell and lobe carried by the non-thermal electrons are $\epsilon_e \sim 0.01$ and $\epsilon_{e, lobe} \sim 1$, respectively. Then, we find that the shell emissions are brighter than the lobe ones at infra-red and optical bands when the source size is > 10kpc, and the IC emissions from the shell at > 10 GeV can be observed with the absence of contamination from the lobe irrespective of the source size. In particular, it is predicted that, for most powerful nearby sources ($L_j \sim 10^{47} ergs s^{-1}$), TeV gamma-rays produced via the IC emissions can be detected by the modern Cherenkov telescopes such as MAGIC, HESS and VERITAS.

We discuss the characteristics of known extremely metal-poor (EMP) stars in the Galaxy using the Stellar Abundances for Galactic Archaeology (SAGA) database (Suda et al. 2008, PASJ, 60, 1159).The analyses of carbon-enhanced stars in our sample suggest that the nucleosynthesis in AGB stars can contribute to the carbon enrichment in a different way depending on whether the metallicity is above or below [Fe/H] ~ -2.5, which is consistent with the current models of stellar evolution at low metallicity. We find the transition of the initial mass function at [Fe/H] ~ -2 in the viewpoint of the distribution of carbon abundance and the frequency of carbon-enhanced stars. For observed EMP stars, we confirmed that some, not all, of observed stars might have undergone at least two types of extra mixing to change their surface abundances. One is to deplete the lithium abundance during the early phase of red giant branch. Another is to decrease the C/N ratio by one order of magnitude during the red giant branch phase. Observed small scatters of abundances for alpha-elements and iron-group elements suggest that the chemical enrichment of our Galaxy takes place in a well-mixed interstellar medium. We find that the abundance trends of alpha-elements are highly correlated with each other, while the abundances of iron-group elements are subject to different slopes relative to the iron abundance. This implies that the supernova yields of alpha-elements are almost independent of mass and metallicity, while those of iron-group elements have a metallicity dependence or mass dependence with the variable initial mass function.The occurrence of the hot bottom burning in the mass range of 5 <~ M / Msun <~ 6 is consistent with the initial mass function of the Galaxy peaked at ~ 10 - 12 Msun to be compatible with the statistics of carbon enhanced stars with and without s-process element (truncated)

In the formation process of black holes, the density and temperature of matter become sufficiently high for quarks and pions to appear. In this study we numerically investigate stellar core collapse and black hole formation taking into account the equations of state involving quarks and/or pions. In our simulations, we utilize a code that solves the general relativistic hydrodynamics and neutrino transfer equations simultaneously, treating neutrino reactions in detail under spherical symmetry. Initial models with three different masses, namely, 40, 100 and 375Msolar, are adopted. Our results show that quarks and pions shorten the duration of neutrino emission if the collapse bounces before black hole formation. In addition, pions increase the luminosity and average energy of neutrinos before black hole formation. We also find that the hadron-quark phase transition leads to an interesting evolution of temperature. Moreover, the neutrino event number is evaluated for the currently operating neutrino detector, SuperKamiokande, to confirm that it is not only detectable but also affected by the emergence of quarks and pions for Galactic events. While there are some issues, such as hyperons, beyond the scope of this study, this is the first serious attempt to assess the impact of quarks and pions in dynamical simulations of black hole formation and will serve as an important foundation for future studies.

Context. White dwarfs (WDs) in cataclysmic variables (CVs) are important experimental laboratories where the electron degeneracy is taking place on a macroscopic scale. Magnetic CVs increase in number especially in the hard X-ray band (>10 keV) thanks to sensitive hard X-ray missions. Aims. From X-ray spectroscopy, we estimate the masses of nearby WDs in moderately-magnetized CVs, or Intermediate Polars (IPs). Methods. Using the Suzaku satellite, we aquired wide-band spectra of 17 IPs, covering 3-50 keV. An accretion column model of Suleimanov et al. (2005) and an optically-thin thermal emission code were used to construct a spectral emission model of IPs with resolved Fe emission lines. By simultaneously fitting the Fe line complex and the hard X-ray continuum of individual spectra, the shock temperature and the WD mass were determined with a better accuracy than in previous studies. Results. We determined the WD masses of the 17 IPs with statistical fitting errors of ~0.1-0.2 Msun in many cases. The WD mass of a recently-found IP, IGR J17195-4100, was also estimated for the first time (1.03+0.24-0.22 Msun). The average WD mass of the sample is 0.88 \pm 0.25 Msun. When our results were compared with previous X-ray mass determinations, we found significant deviation in a few systems although the reason of this is unclear. The iron abundance of the accreting gas was also estimated, and confirmed the previously reported sub-solar tendency in all sources with better accuracy.

The black-hole binary system LMC X-3 has been observed by virtually every X-ray mission since the inception of X-ray astronomy. Among the persistent sources, LMC X-3 is uniquely both habitually soft and highly variable. Using a fully relativistic accretion-disk model, we analyze hundreds of spectra collected during eight X-ray missions that span 26 years. For a selected sample of 391 RXTE spectra we find that to within ~2 percent the inner radius of the accretion disk is constant over time and unaffected by source variability. Even considering an ensemble of eight X-ray missions, we find consistent values of the radius to within ~4-6 percent. Our results provide strong evidence for the existence of a fixed inner-disk radius. The only reasonable inference is that this radius is closely associated with the general relativistic innermost stable circular orbit (ISCO). Our findings establish a firm foundation for the measurement of black hole spin.

We present a new formula to numerically construct configurations in rotational equilibrium, which consist of multiple layers. Each layer rotates uniformly or differentially according to cylindrical rotation-laws that are different from layer to layer. Assuming a different barotropic equation of state (EOS) for each layer, we solve the Bernoulli equation in each layer separately and combine the solutions by imposing continuity of the pressure at each boundary of the layers. It is confirmed that a single continuous barotropic EOS is incompatible with the junction condition. Identifying appropriate variables to be solved, we construct a convergent iteration scheme. For demonstration, we obtain two-layered configurations, each layer of which rotates rapidly with either an "$\Omega$-constant law" or a "$j$-constant law" or a "$v$-constant law". Other rotation laws and/or a larger number of layers can be treated similarly. We hope that this formula will be useful in studying the stellar evolution in multi-dimension with the non-spherical configuration induced by rotation being fully taken into account.

We explore the emissions by accelerated electrons in shocked shells driven by jets in active galactic nuclei (AGNs). Focusing on powerful sources which host luminous quasars, the synchrotron radiation and inverse Compton (IC) scattering of various photons that are mainly produced in the core are considered as radiation processes. We show that the radiative output is dominated by the IC emission for compact sources (< 30kpc), whereas the synchrotron radiation is more important for larger sources. It is predicted that, for powerful sources ($L_j \sim 10^{47} ergs s^{-1}$), GeV-TeV gamma-rays produced via the IC emissions can be detected by the Fermi satellite and modern Cherenkov telescopes such as MAGIC, HESS and VERITAS if the source is compact.

The detection of neutrinos from massive stellar collapses can teach us a lot not only about source objects but also about microphysics working deep inside them. In this study we discuss quantitatively the possibility to extract information on the properties of dense and hot hadronic matter from neutrino signals coming out of black-hole-forming collapses of non-rotational massive stars. Based on our detailed numerical simulations we evaluate the event numbers for SuperKamiokande with neutrino oscillations being fully taken into account. We demonstrate that the event numbers from a Galactic event are large enough not only to detect it but also to distinguish one hadronic equation of state from another by our statistical method assuming the same progenitor model and non-rotation. This means that the massive stellar collapse can be a unique probe into hadron physics and will be a promising target of the nascent neutrino astronomy.

We study in detail the motions of three planets interacting with each other under the influence of a central star. It is known that the system with more than two planets becomes unstable after remaining quasi-stable for long times, leading to highly eccentric orbital motions or ejections of some of the planets. In this paper, we are concerned with the underlying physics for this quasi-stability as well as the subsequent instability and advocate the so-called "stagnant motion" in the phase space, which has been explored in the field of dynamical system. We employ the Lyapunov exponent, the power spectra of orbital elements and the distribution of the durations of quasi-stable motions to analyze the phase space structure of the three-planet system, the simplest and hopefully representative one that shows the instability. We find from the Lyapunov exponent that the system is almost non-chaotic in the initial quasi-stable state whereas it becomes intermittently chaotic thereafter. The non-chaotic motions produce the horizontal dense band in the action-angle plot whereas the voids correspond to the chaotic motions. We obtain power laws for the power spectra of orbital eccentricities. Power-law distributions are also found for the durations of quasi-stable states. All these results combined together, we may reach the following picture: the phase space consists of the so-called KAM tori surrounded by satellite tori and imbedded in the chaotic sea. The satellite tori have a self-similar distribution and are responsible for the scale-free power-law distributions of the duration times. The system is trapped around one of the KAM torus and the satellites for a long time (the stagnant motion) and moves to another KAM torus with its own satellites from time to time, corresponding to the intermittent chaotic behaviors.

We study how the influence of the shock wave appears in neutrino oscillations and the neutrino spectrum using density profile of adiabatic explosion model of a core-collapse supernova which is calculated in an implicit Lagrangian code for general relativistic spherical hydrodynamics. We calculate expected event rates of neutrino detection at SK and SNO for various theta_{13} values and both normal and inverted hierarchies. The predicted event rates of bar{nu}_e and nu_e depend on the mixing angle theta_{13} for the inverted and normal hierarchies, respectively, and the influence of the shock appears for about 2 - 8 s when sin^2 2 theta_{13} is larger than 10^{-3}. These neutrino signals for the shock propagation is decreased by < 30 % for bar{nu}_e in inverted (SK) or by < 15 % for nu_e in normal hierarchy (SNO) compared with the case without shock. The obtained ratio of the total event for high-energy neutrinos (20 MeV < E_{nu} < 60 MeV) to low-energy neutrinos (5 MeV < E_{nu} < 20 MeV) is consistent with the previous studies in schematic semi-analytic or other hydrodynamic models of the shock propagation. The time dependence of the calculated ratio of the event rates of high-energy to low-energy neutrinos is a very useful observable which is sensitive to theta_{13} and hierarchies. Namely, time-dependent ratio shows clearer signal of the shock propagation that exhibits remarkable decrease by at most factor \sim 2 for bar{nu}_e in inverted (SK), whereas it exhibits smaller change by \sim 10 % for nu_e in normal hierarchy (SNO). Observing time-dependent high-energy to low-energy ratio of the neutrino events thus would provide a piece of very useful information to constrain theta_{13} and mass hierarchy, and eventually help understanding the propagation how the shock wave propagates inside the star.

We report results of the 2006 April multi-wavelengths campaign of SS 433, focusing on X-ray data observed with Suzaku at two orbital phases (in- and out-of- eclipse) and simultaneous optical spectroscopic observations. By analyzing the Fe25 K_alpha lines originating from the jets, we detect rapid variability of the Doppler shifts, dz/dt ~ 0.019/0.33 day^-1, which is larger than those expected from the precession and/or nodding motion. This phenomenon probably corresponding to "jitter" motions observed for the first time in X-rays, for which significant variability both in the jet angle and intrinsic speed is required. From the time lag of optical Doppler curves from those of X-rays, we estimate the distance of the optical jets from the base to be ~(3-4) \times 10^14 cm. Based on the radiatively cooling jet model, we determine the innermost temperature of the jets to be T_0 = 13 +/- 2 keV and 16 +/- 3 keV (the average of the blue and red jets) for the out-of-eclipse and in-eclipse phase, respectively, from the line intensity ratio of Fe25 K_alpha and Fe26 K_alpha. While the broad band continuum spectra over the 5--40 keV band in eclipse is consistent with a multi-temperature bremsstrahlung emission expected from the jets, and its reflection component from cold matter, the out-of-eclipse spectrum is harder than the jet emission with the base temperature determined above, implying the presence of an additional hard component.

The wide-band Suzaku spectra of the black hole binary GX 339-4, acquired in 2007 February during the Very High state, were reanalyzed. Effects of event pileup (significant within ~ 3' of the image center) and telemetry saturation of the XIS data were carefully considered. The source was detected up to ~ 300$ keV, with an unabsorbed 0.5--200 keV luminosity of ~3.8 10^{38} erg/s at 8 kpc. The spectrum can be approximated by a power-law of photon index 2.7, with a mild soft excess and a hard X-ray hump. When using the XIS data outside 2' of the image center, the Fe-K line appeared extremely broad, suggesting a high black hole spin as already reported by Miller et al. (2008) based on the Suzaku data and other CCD data. When the XIS data accumulation is further limited to >3' to avoid event pileup, the Fe-K profile becomes narrower, and there appears a marginally better solution that suggests the inner disk radius to be 5-14 times the gravitational radius (1-sigma), though a maximally spinning black hole is still allowed by the data at the 90% confidence level. Consistently, the optically-thick accretion disk is inferred to be truncated at a radius 5-32 times the gravitational radius. Thus, the Suzaku data allow an alternative explanation without invoking a rapidly spinning black hole. This inference is further supported by the disk radius measured previously in the High/Soft state.

We propose a ray-tracing method to estimate gravitational waves (GWs) generated by anisotropic neutrino emission in supernova cores. To calculate the gravitational waveforms, we derive analytic formulae in a useful form, which are applicable also for three-dimensional computations. Pushed by evidence of slow rotation prior to core-collapse, we focus on asphericities in neutrino emission and matter motions outside the protoneutron star. Based on the two-dimensional (2D) models, which mimic SASI-aided neutrino heating explosions, we compute the neutrino anisotropies via the ray-tracing method in a post-processing manner and calculate the resulting waveforms. With these computations, it is found that the waveforms exhibit more variety in contrast to the ones previously estimated by the ray-by-ray analysis (e.g., Kotake et al. (2007)). In addition to a positively growing feature, which was predicted to determine the total wave amplitudes predominantly, the waveforms are shown to exhibit large negative growth for some epochs during the growth of SASI. These features are found to stem from the excess of neutrino emission in lateral directions, which can be precisely captured by the ray-tracing calculation. Due to the negative contributions and the neutrino absorptions appropriately taken into account by the ray-tracing method, the wave amplitudes become more than one-order-of magnitude smaller than the previous estimation, thus making their detections very hard for a galactic source.On the other hand, it is pointed out that the GW spectrum from matter motions have its peak near $\sim 100$ Hz, which could be characteristic for the SASI-induced supernova explosions.(abridged)

Nuclear matter is considered to be inhomogeneous at subnuclear densities that are realized in supernova cores and neutron star crusts, and the structures of nuclear matter change from spheres to cylinders, slabs, cylindrical holes and spherical holes as the density increases. In this letter, we discuss other possible structures, that is, gyroid and double-diamond morphologies, which are periodic bicontinuous structures discovered in a block copolymer. Utilizing the compressible liquid drop model, we evaluate their surface and Coulomb energies and show that there is a chance of gyroid appearance near the transition point from a cylinder to a slab. This interesting analogy between nuclear and polymer systems is not merely qualitative. The volume fraction at the phase transition is also similar for the two systems. Although the five shapes listed initially have been long thought to be the only major constituents of so-called nuclear pasta at subnuclear densities, our findings imply that this may not be the case and suggest that more detailed studies on nuclear pasta including the gyroid phase are needed.

We reinterpret the proof of the Riemannian Penrose inequality by H. Bray. The modified argument turns out to have a nice feature so that the flow of Riemannian metrics appearing Bray's proof gives a Lorentzian metric of a spacetime. We also discuss a possible extension of our approach to charged black holes.

The Suzaku X-ray satellite observed the nearby spiral galaxy NGC 4258 for a total good exposure time of 100 ks. We present an analysis of the Suzaku XIS data, in which we confirm that the 0.5--2 keV spectra of the interstellar medium (ISM) are well-represented by a two-temperature model. The cool and hot ISM temperatures are 0.23+-0.02 and 0.59 +-0.01 keV, respectively. Suzaku's excellent spectral sensitivity enables us to measure the metal abundances of O, Ne, Mg, Si and Fe of the ISM for the first time. The resultant abundance pattern of O, Mg, Si, and Fe is consistent with that of the new solar abundance table of Lodders (2003), rather than Anders & Grevesse (1989). This suggests that the metal enrichment processes of NGC 4258 and of our Galaxy are similar.

We study properties of gravitational waves based on the three-dimensional simulations, which demonstrate the neutrino-driven explosions aided by the standing accretion shock instability (SASI). Pushed by evidence supporting slow rotation prior to core-collapse, we focus on the asphericities in neutrino emissions and matter motions outside the protoneutron star. By performing a ray-tracing calculation in 3D, we estimate accurately the gravitational waveforms from anisotropic neutrino emissions. In contrast to the previous work assuming axisymmetry, we find that the gravitational waveforms vary much more stochastically because the explosion anisotropies depend sensitively on the growth of the SASI which develops chaotically in all directions. Our results show that the gravitational-wave spectrum has its peak near $\sim 100$ Hz, reflecting the SASI-induced matter overturns of $\sim O(10)$ ms. We point out that the detection of such signals, possibly visible to the LIGO-class detectors for a Galactic supernova, could be an important probe into the long-veiled explosion mechanism.

This paper is a sequel to our previous work for accretion onto a Schwarzschild black hole and the so-called standing accretion shock instability (SASI), in this paper we investigate non-axisymmetric perturbations for a Kerr black hole. The linear and non-linear phases for the shock evolution are analyzed in detail by both 2D general relativistic hydrodynamical simulations and linear analysis. Since the structure of steady axisymmetric accretion flows with a standing shock wave is very sensitive to the inner transonic flow, their properties such as Mach numbers, which are important for the stability, depend on the Kerr parameter very much. Although the essential features of the instability do not differ from the previous results for the Schwarzschild black hole, the frame dragging effects specific to the Kerr black hole is also evident. Interestingly, the oscillation periods of the fundamental unstable modes are dependent only on the shock radius irrespective of the injection parameters.

We present the first Suzaku observation of the new Soft Gamma Repeater SGR 0501+4516, performed on 2008 August 26, four days after the onset of bursting activity of this new member of the magnetar family. The soft X-ray persistent emission was detected with the X-ray Imaging Spectrometer (XIS) at a 0.5-10 keV flux of 3.8E-11 erg/s/cm2, with a spectrum well fitted by an absorbed blackbody plus power-law model. The source pulsation was confirmed at a period of 5.762072+/-0.000002 s, and 32 X-ray bursts were detected by the XIS, four of which were also detected at higher energies by the Hard X-ray Detector (HXD). The strongest burst, which occurred at 03:16:16.9 (UTC), was so bright that it caused instrumental saturation, but its precursor phase, lasting for about 200 ms, was detected successfully over the 0.5-200 keV range, with a fluence of ~2.1E-7 erg/cm2 and a peak intensity of about 89 Crab. The entire burst fluence is estimated to be ~50 times higher. The precursor spectrum was very hard, and well modeled by a combination of two blackbodies. We discuss the bursting activity and X/gamma-ray properties of this newly discovered Soft Gamma Repeater in comparison with other members of the class.

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.

In 2006 June, the obscured low luminosity active galactic nucleus in the nearby Seyfert 1.9 galaxy NGC 4258 was observed with Suzaku for ~ 100 ks. Utilizing the XIS and the HXD, the nucleus emission was detected over 2 to 40 keV range, with an unabsorbed 2--10 keV luminosity of 8 x 10 40 erg / s, and varied by a factor of ~ 2 during the observation. Its 2--40 keV spectrum is reproduced by a single power law with photon index of ~ 2.0, absorbed by an equivalent hydrogen column of ~ 1.0 x 10 23 cm2. The spectrum within 4' of the nucleus required also a softer thin-thermal emission, as well as an intermediate hardness component attributable to integrated point sources. A weak neutral Fe-Kalpha florescence line was detected at an equivalent width of ~ 40 eV. The cold reflection component was not required by the data, with the reflector solid angle Omega seen from the nucleus constrained as Omega / 2 pi < 0.3 assuming a general case of 60 deg inclination. The results suggest that the cold reflecting material around the nucleus is localized along our line of sight, rather than forming a thick torus.

We investigate the emergence of strange baryons in the dynamical collapse of a non-rotating massive star to a black hole by the neutrino-radiation hydrodynamical simulations in general relativity. By following the dynamical formation and collapse of nascent proto-neutron star from the gravitational collapse of a 40Msun star adopting a new hyperonic EOS table, we show that the hyperons do not appear at the core bounce but populate quickly at ~0.5-0.7 s after the bounce to trigger the re-collapse to a black hole. They start to show up off center owing to high temperatures and later prevail at center when the central density becomes high enough. The neutrino emission from the accreting proto-neutron star with the hyperonic EOS stops much earlier than the corresponding case with a nucleonic EOS while the average energies and luminosities are quite similar between them. These features of neutrino signal are a potential probe of the emergence of new degrees of freedom inside the black hole forming collapse.

The broad-band 1-300 keV Suzaku spectrum of IGR J17497-2821, the X-ray transient discovered by INTEGRAL in September 2006, is presented. Suzaku observed IGR J17497-2821 on September 25, eight days after its discovery, for a net exposure of about 53 ksec. During the Suzaku observation, IGR J17497-2821 is very bright, 2 x 10^37 erg/s at 8 kpc in the 1-300 keV range, and shows a hard spectrum, typical of black hole candidates in the low-hard state. Despite the multi-mission X-ray monitoring of the source, only with Suzaku is it possible to obtain a broad-band spectrum in the 1-300 keV range with a very high signal to noise ratio. A sum of a multi-color disc (DISKBB) and a thermal Comptonization component (COMPPS) with mild reflection is a good representation of our IGR J17497-2821 Suzaku spectrum. The spectral properties of the accretion disc as well as the cut-off energy in the spectrum at about 150 keV are clearly detected and constrained. We discuss the implications on the physical model used to interpret the data and the comparison with previous results.

We studied the effects of rotation on standing accretion shock instability (SASI) by performing three-dimensional hydrodynamics simulations. Taking into account a realistic equation of state and neutrino heating/cooling, we prepared a spherically symmmetric and steady accretion flow through a standing shock wave onto a proto-neutron star (PNS). When the SASI entered the nonlinear phase, we imposed uniform rotation on the flow advecting from the outer boundary of the iron core, whose specific angular momentum was assumed to agree with recent stellar evolution models. Using spherical harmonics in space and Fourier decompositions in time, we performed mode analysis of the nonspherical deformed shock wave to observe rotational effects on the SASI in the nonlinear phase. We found that rotation imposed on the axisymmetric SASI did not make any spiral modes and hardly affected sloshing modes, except for steady l=2, m=0 modes. In contrast, rotation imposed on the non-axisymmetric flow increased the amplitude of spiral modes so that some spiral flows accreting on the PNS were more clearly formed inside the shock wave than without rotation. The amplitudes of spiral modes increased significantly with rotation in the progressive direction.

Recently, stellar collapse involving black hole formation from massive stars is suggested to emit enormous fluxes of neutrinos on par with ordinary core-collapse supernovae. We investigate their detectability for the currently operating neutrino detector, SuperKamiokande. Neutrino oscillation is also taken into account for the evaluation. We find that the event number is larger than or comparable to that of supernova neutrinos and, hence, black hole formation is also a candidate for neutrino astronomy. Moreover, we find that the event number depends dominantly on the equation of state used in the computations of the black hole formation. This fact implies that the detection of neutrinos emitted from the black hole progenitors is very valuable to probe the properties of the equation of state for hot and/or dense matter.

A report is presented on Suzaku observations of the ultra-luminous X-ray source X-1 in the starburst galaxy M82, made three time in 2005 October for an exposure of ~ 30 ks each. The XIS signals from a region of radius 3 around the nucleus defined a 2-10 keV flux of 2.1 x 10^-11 erg s-1 cm-2 attributable to point sources. The 3.2-10 keV spectrum was slightly more convex than a power-law with a photon index of 1.7. In all observations, the HXD also detected signals from M82 up to ~ 20 keV, at a 12-20 keV flux of 4.4 x 10^-12 erg s-1 cm-2 . The HXD spectrum was steeper than that of the XIS. The XIS and HXD spectra can be jointly reproduced by a cutoff power-law model, or similar curved models. Of the detected wide-band signals, 1/3 to 2/3 are attributable to X-1, while the remainder to other discrete sources in M82. Regardless of the modeling of these contaminants, the spectrum attributable to X-1 is more curved than a power-law, with a bolometric luminosity of (1.5 -3) x 10 ^40 erg s-1. These results are interpreted as Comptonized emission from a black hole of 100-200 solar masses, radiating roughly at the Eddington luminosity.

We study the stability of standing shock waves in advection-dominated accretion flows into a Schwarzschild black hole by 2D general relativistic hydrodynamic simulations as well as linear analysis in the equatorial plane. We demonstrate that the accretion shock is stable against axisymmetric perturbations but becomes unstable to non-axisymmetric perturbations. The results of dynamical simulations show good agreement with linear analysis on the stability, oscillation and growing time scales. The comparison of different wave-travel times with the growth time scales of the instability suggests that the instability is likely to be of the Papaloizou-Pringle type, induced by the repeated propagations of acoustic waves. However, the wavelengths of perturbations are too long to clearly define the reflection point. By analyzing the non-linear phase in the dynamical simulations, it is shown that quadratic mode couplings precede the non-linear saturation. It is also found that not only short-term random fluctuations by turbulent motions but also quasi periodic oscillations take place on longer time scales in the non-linear phase. We give some possible implications of the instability for quasi periodic oscillations (QPOs) and the central engine for gamma ray bursts (GRBs).

We study the progenitor dependence of the black hole formation and its associated neutrino signals from the gravitational collapse of non-rotating massive stars, following the preceding study on the single progenitor model in Sumiyoshi et al. (2007). We aim to clarify whether the dynamical evolution toward the black hole formation occurs in the same manner for different progenitors and to examine whether the characteristic of neutrino bursts is general having the short duration and the rapidly increasing average energies. We perform the numerical simulations by general relativistic neutrino-radiation hydrodynamics to follow the dynamical evolution from the collapse of pre-supernova models of 40Msun and 50Msun toward the black hole formation via contracting proto-neutron stars. For the three progenitor models studied in this paper, we found that the black hole formation occurs in ~0.4-1.5 s after core bounce through the increase of proto-neutron star mass together with the short and energetic neutrino burst. We found that density profile of progenitor is important to determine the accretion rate onto the proto-neutron star and, therefore, the duration of neutrino burst. We compare the neutrino bursts of black hole forming events from different progenitors and discuss whether we can probe clearly the progenitor and/or the dense matter.

We investigate the total kinetic powers (L_j) and ages (t_age) of powerful jets in FR II radio galaxies by comparison of the dynamical model of expanding cocoons with observations. We select four FR II radio sources (Cygnus A, 3C 223, 3C 284, and 3C 219), for which the mass-density profiles of intracluster medium (ICM) are known in the literature. It is found that large fractions > 0.02 - 0.7 of the Eddington luminosity (L_Edd) are carried away as a kinetic power of jet. The upper limit of estimated 2 L_j / L_Edd are larger than unity (< 10) for some sources, suggesting a possibility of super-Eddington mass accretions. As a consequence of the large powers, we also find that the total energy stored in the cocoon (E_c) exceeds the energy derived from the minimum energy condition for the energy of radiating non-thermal electrons and magnetic fields (E_min): 4< E_c /E_min <310. This implies that most of the energy in cocoon is carried by invisible components such as thermal leptons (electron and positron) and/or protons.

We describe the construction of a database of extremely metal-poor (EMP) stars in the Galactic halo whose elemental abundances have been determined. Our database contains detailed elemental abundances, reported equivalent widths, atmospheric parameters, photometry, and binarity status, compiled from papers in the recent literature that report studies of EMP halo stars with [Fe/H] < -2.5. The compilation procedures for this database have been designed to assemble the data effectively from electronic tables available from online journals. We have also developed a data retrieval system that enables data searches by various criteria, and permits the user to explore relationships between the stored variables graphically. Currently, our sample includes 1212 unique stars (many of which are studied by more than one group) with more than 15000 individual reported elemental abundances, covering all of the relevant papers published by December 2007. We discuss the global characteristics of the present database, as revealed by the EMP stars observed to date. For stars with [Fe/H] < -2.5, the number of giants with reported abundances is larger than that of dwarfs by a factor of two. The fraction of carbon-rich stars (among the sample for which the carbon abundance is reported) amount to ~30 % for [Fe/H] < -2.5. We find that known binaries exhibit different distributions of orbital period, according to whether they are giants or dwarfs, and also as a function of metallicity, although the total sample of such stars is still quite small.

We construct an equation of state including the hadron-quark phase transition. The mixed phase is obtained by the Gibbs conditions for finite temperature. We adopt the equation of state based on the relativistic mean field theory for the hadronic phase taking into account pions. As for the quark phase, the MIT bag model of the deconfined 3-flavor strange quark matter is used. As a result, our equation of state is thermodynamically stable and exhibits qualitatively the desired properties of hadron-quark mixed matter, such as the temperature dependence of the transition density. The pions raise the transition density because they make the equation of state softer. Using the equation of state constructed here, we study its astrophysical implications. The maximum mass of compact stars is investigated, and our equation of state is consistent with recent observations. We also compute the collapse of a massive star with 100 solar masses using our equation of state and find that the interval time from the bounce to the black hole formation becomes shorter for the model with pions and quarks. The pions and quarks affect the total energy of the emitted neutrinos because the duration time of the neutrino emission becomes shorter. The neutrino luminosity rises under the effect of pions since the density of the proto-neutron star becomes high.

The black-hole binary Cygnus X-1 was observed for 17 ks with the Suzaku X-ray observatory in 2005 October, while it was in a low/hard state with a 0.7-300 keV luminosity of 4.6 x 10^37 erg/s. The XIS and HXD spectra, spanning 0.7-400 keV, were reproduced successfully incorporating a cool accretion disk and a hot Comptonizing corona. The corona is characterized by an electron temperature of ~100 keV, and two optical depths of ~0.4 and ~1.5 which account for the harder and softer continua, respectively. The disk has the innermost temperature of ~0.2 keV, and is thought to protrude half way into the corona. The disk not only provides seed photons to the Compton cloud, but also produces a soft spectral excess, a mild reflection hump, and a weakly broadened iron line. A comparison with the Suzaku data on GRO J1655-40 reveals several interesting spectral differences, which can mostly be attributed to inclination effects assuming that the disk has a flat geometry while the corona is grossly spherical. An intensity-sorted spectroscopy indicates that the continuum becomes less Comptonized when the source flares up on times scales of 1-200 s, while the underlying disk remains unchanged.

We investigate the tidal interactions of a red giant with a main sequence in the dense stellar core of globular clusters by Smoothed Particle Hydrodynamics method. Two models of $0.8 \msun$ red giant with the surface radii 20 and $85 R_\sun$ are used with 0.6 or $0.8M_\sun$ main sequence star treated as a point mass. We demonstrate that even for the wide encounters that two stars fly apart, the angular momentum of orbital motion can be deposited into the red giant envelope to such an extent as to trigger rotational mixing and to explain the fast rotation observed for the horizontal branch stars, and also that sufficient mass can be accreted on the main sequence stars to disguise their surface convective zone with the matter from the red giant envelope. On the basis of the present results, we discuss the parameter dependence of these transfer characteristics with non-linear effects taken into account, and derive fitting formulae to give the amounts of energy and angular momentum deposited into the red giant and of mass accreted onto the perturber as functions of stellar parameters and the impact parameter of encounter. These formulae are applicable to the encounters not only of the red giants but also of the main sequence stars, and useful in the study of the evolution of stellar systems with the star-star interactions taken into account.

We investigate the role of Alfven waves in the core-collapse supernova (SN) explosion. We assume that Alfven waves are generated by convections inside a proto-neutron star (PNS) and emitted from its surface. Then these waves propagate outwards, dissipate via nonlinear processes, and heat up matter around a stalled prompt shock. To quantitatively assess the importance of this process for the revival of the stalled shock, we perform 1D time-dependent hydrodynamical simulations, taking into account the heating via the dissipation of Alfven waves that propagate radially outwards along open flux tubes. We show that the shock revival occurs if the surface field strength is larger than ~2e15 G and if the amplitude of velocity fluctuation at the PNS surface is larger than 20% of the local sound speed. Interestingly, the Alfven wave mechanism is self-regulating in the sense that the explosion energy is not very sensitive to the surface field strength and initial amplitude of Alfven waves as long as they are larger than the threshold values given above.

We have studied non-axisymmetric standing accretion shock instability, or SASI, by 3D hydrodynamical simulations. This is an extention of our previous study on axisymmetric SASI. We have prepared a spherically symmetric and steady accretion flow through a standing shock wave onto a proto-neutron star, taking into account a realistic equation of state and neutrino heating and cooling. This unperturbed model is supposed to represent approximately the typical post-bounce phase of core-collapse supernovae. We then have added a small perturbation (~1%) to the radial velocity and computed the ensuing evolutions. Not only axisymmetric but non-axisymmetric perturbations have been also imposed. We have applied mode analysis to the non-spherical deformation of the shock surface, using the spherical harmonics. We have found that (1) the growth rates of SASI are degenerate with respect to the azimuthal index m of the spherical harmonics Y_l^m, just as expected for a spherically symmetric background, (2) nonlinear mode couplings produce only m=0 modes for the axisymmetric perturbations, whereas m=!0 modes are also generated in the non-axisymmetric cases according to the selection rule for the quadratic couplings, (3) the nonlinear saturation level of each mode is lower in general for 3D than for 2D because a larger number of modes are contributing to turbulence in 3D, (4) low l modes are dominant in the nonlinear phase, (5) the equi-partition is nearly established among different m modes in the nonlinear phase, (6) the spectra with respect to l obey power laws with a slope slightly steeper for 3D, and (7) although these features are common to the models with and without a shock revival at the end of simulation, the dominance of low l modes is more remarkable in the models with a shock revival.

A series of numerical simulations on magnetorotational core-collapse supernovae are carried out. Dipole-like configurations which are offset northward are assumed for the initially strong magnetic fields together with rapid differential rotations. Aims of our study are to investigate effects of the offset magnetic field on magnetar kicks and on supernova dynamics. Note that we study a regime where the proto-neutron star formed after collapse has a large magnetic field strength approaching that of a ``magnetar'', a highly magnetized slowly rotating neutron star. As a result, equatorially-asymmetric explosions occur with a formation of the bipolar jets. Resultant magnetar's kick velocities are $\sim 300-1000$ km s$^{-1}$. We find that the acceleration is mainly due to the magnetic pressure while the somewhat weaker magnetic tension works toward the opposite direction, which is due to stronger magnetic field in the northern hemisphere. Noted that observations of magnetar's proper motions are very scarce, our results supply a prediction for future observations. Namely, magnetars possibly have large kick velocities, several hundred km s$^{-1}$, as ordinary neutron stars do, and in an extreme case they could have those up to 1000 km s$^{-1}$.

A report is made on a comprehensive observation of a burst-like $\gamma$-ray emission from thunderclouds on the Sea of Japan, during strong thunderstorms on 2007 January 6. The detected emission, lasting for $\sim$40 seconds, preceded cloud-to-ground lightning discharges. The burst spectrum, extending to 10 MeV, can be interpreted as consisting of bremsstrahlung photons originating from relativistic electrons. This ground-based observation provides first clear evidence that strong electric fields in thunderclouds can continuously accelerate electrons beyond 10 MeV prior to lightning discharges.

We investigate the total kinetic powers (L_{j}) and ages (t_{age}) of powerful jets of four FR II radio sources (Cygnus A, 3C 223, 3C 284, and 3C 219) by the detail comparison of the dynamical model of expanding cocoons with observed ones. It is found that these sources have quite large kinetic powers with the ratio of L_{j} to the Eddington luminosity (L_{Edd}) resides in $0.02 <L_{j}/L_{Edd} <10$. Reflecting the large kinetic powers, we also find that the total energy stored in the cocoon (E_{c}) exceed the energy derived from the minimum energy condition (E_{min}): $2< E_{c}/E_{min} <160$. This implies that a large amount of kinetic power is carried by invisible components such as thermal leptons (electron and positron) and/or protons.

We study the black hole formation and the neutrino signal from the gravitational collapse of a non-rotating massive star of 40 Msun. Adopting two different sets of realistic equation of state (EOS) of dense matter, we perform the numerical simulations of general relativistic neutrino-radiation hydrodynamics under the spherical symmetry. We make comparisons of the core bounce, the shock propagation, the evolution of nascent proto-neutron star and the resulting re-collapse to black hole to reveal the influence of EOS. We also explore the influence of EOS on the neutrino emission during the evolution toward the black hole formation. We find that the speed of contraction of the nascent proto-neutron star, whose mass increases fast due to the intense accretion, is different depending on the EOS and the resulting profiles of density and temperature differ significantly. The black hole formation occurs at 0.6-1.3 sec after bounce when the proto-neutron star exceeds its maximum mass, which is crucially determined by the EOS. We find that the average energies of neutrinos increase after bounce because of rapid temperature increase, but at different speeds depending on the EOS. The duration of neutrino emission up to the black hole formation is found different according to the different timing of re-collapse. These characteristics of neutrino signatures are distinguishable from those for ordinary proto-neutron stars in successful core-collapse supernovae. We discuss that a future detection of neutrinos from black-hole-forming collapse will contribute to reveal the black hole formation and to constrain the EOS at high density and temperature.

The existence of various anomalous stars, such as the first stars in the universe or stars produced by stellar mergers, has been recently proposed. Some of these stars will result in black hole formation. In this study, we investigate iron core collapse and black hole formation systematically for the iron-core mass range of 3 - 30Msun, which has not been studied well so far. Models used here are mostly isentropic iron cores that may be produced in merged stars in the present universe but we also employ a model that is meant for a Population III star and is obtained by evolutionary calculation. We solve numerically the general relativistic hydrodynamics and neutrino transfer equations simultaneously, treating neutrino reactions in detail under spherical symmetry. As a result, we find that massive iron cores with ~10Msun unexpectedly produce a bounce owing to the thermal pressure of nucleons before black hole formation. The features of neutrino signals emitted from such massive iron cores differ in time evolution and spectrum from those of ordinary supernovae. Firstly, the neutronization burst is less remarkable or disappears completely for more massive models because the density is lower at the bounce. Secondly, the spectra of neutrinos, except the electron type, are softer owing to the electron-positron pair creation before the bounce. We also study the effects of the initial density profile, finding that the larger the initial density gradient is, the more steeply the neutronization burst declines. Further more, we suggest a way to probe into the black hole progenitors from the neutrino emission and estimate the event number for the currently operating neutrino detectors.

We perform two-dimensional, magnetohydrodynamical core-collapse simulations of massive stars accompanying the QCD phase transition. We study how the phase transition affects the gravitational waveforms near the epoch of core-bounce. As for initial models, we change the strength of rotation and magnetic fields. Particularly, the degree of differential rotation in the iron core (Fe-core) is changed parametrically. As for the microphysics, we adopt a phenomenological equation of state above the nuclear density, including two parameters to change the hardness before the transition. We assume the first order phase transition, where the conversion of bulk nuclear matter to a chirally symmetric quark-gluon phase is described by the MIT bag model. Based on these computations, we find that the phase transition can make the maximum amplitudes larger up to $\sim$ 10 percents than the ones without the phase transition. On the other hand, the maximum amplitudes become smaller up to $\sim$ 10 percents owing to the phase transition, when the degree of the differential rotation becomes larger. We find that even extremely strong magnetic fields $\sim 10^{17}$ G in the protoneutron star do not affect these results.

The so-called "acoustic revival mechanism" of core-collapse supernova proposed recently by the Arizona group is an interesting new possibility. Aiming to understand the elementary processes involved in the mechanism, we have calculated the eigen frequencies and eigen functions for the g-mode oscillations of a non-rotating proto neutron star. The possible excitation of these modes by the standing accretion shock instability, or SASI, is discussed based on these eigen functions. We have formulated the forced oscillations of $g$-modes by the external pressure perturbations exerted on the proto neutron star surface. The driving pressure fluctuations have been adopted from our previous computations of the axisymmetric SASI in the non-linear regime. We have paid particular attention to low l modes, since these are the modes that are dominant in SASI and that the Arizona group claimed played an important role in their acoustic revival scenario. Here l is the index of the spherical harmonic functions, $Y_l^m$. Although the frequency spectrum of the non-linear SASI is broadened substantially by non-linear couplings, the typical frequency is still much smaller than those of g-modes, the fact leading to a severe impedance mismatch. As a result, the excitations of various $g$-modes are rather inefficient and the energy of the saturated g-modes is $\sim 10^{50}$erg or smaller, with the g_2-mode being the largest in our model. Here the g_2-mode has two radial nodes and is confined to the interior of the convection region. The energy transfer rate from the g-modes to out-going sound waves is estimated from the growth of the g-modes and found to be $\sim 10^{51}$erg/s in the model studied in this paper.

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 present the results of six Suzaku observations of the recurrent black hole transient 4U1630-472 during its decline from outburst from February 8 to March 23 in 2006. All observations show the typical high/soft state spectral shape in the 2-50keV band, roughly described by an optically thick disk spectrum in the soft energy band plus a weak power-law tail that becomes dominant only above \~20keV. The disk temperature decreases from 1.4keV to 1.2keV as the flux decreases by a factor 2, consistent with a constant radius as expected for disk-dominated spectra. All the observations reveal significant absorption lines from highly ionized (H-like and He-like) iron Ka at 7.0keV and 6.7keV. The brightest datasets also show significant but weaker absorption structures between 7.8keV and 8.2keV, which we identify as a blend of iron Kb and nickel Ka absorption lines. The energies of these absorption lines suggest a blue shift with an outflow velocity of ~1000km/s.. The H-like iron Ka equivalent width remains approximately constant at ~30 eV over all the observations, while that of the He-like Ka line increases from 7 eV to 20 eV. Thus the ionization state of the material decreases, as expected from the decline in flux. We fit the profiles with Voigt functions (curve of growth) to derive absorbing columns, then use these together with detailed photo-ionization calculations to derive physical parameters of the absorbing material. The data then constrain the velocity dispersion of the absorber to 200-2000 km/s, and the size of the plasma as ~1E10 cm assuming a source distance of 10 kpc.

It is expected that matter composed of a perfect fluid cannot be at rest outside of a black hole if the spacetime is asymptotically flat and static (non-rotating). However, there has not been a rigorous proof for this expectation without assuming spheical symmetry. In this paper, we provide a proof of non-existence of matter composed of a perfect fluid in static black hole spacetimes under certain conditions, which can be interpreted as a relation between the stellar mass and the black hole mass.

We present the results of numerical experiments, in which we study how the asphericities induced by the growth of the standing accretion shock instability (SASI) produce the gravitational waveforms in the postbounce phase of core-collapse supernovae. To obtain the neutrino-driven explosions, we parameterize the neutrino fluxes emitted from the central protoneutron star and approximate the neutrino transfer by a light-bulb scheme. We find that the waveforms due to the anisotropic neutrino emissions show the monotonic increase with time, whose amplitudes are up to two order-of-magnitudes larger than the ones from the convective matter motions outside the protoneutron stars. We point out that the amplitudes begin to become larger when the growth of the SASI enters the nonlinear phase, in which the deformation of the shocks and the neutrino anisotropy become large. From the spectrum analysis of the waveforms, we find that the amplitudes from the neutrinos are dominant over the ones from the matter motions at the frequency below $\sim 100$ Hz, which are suggested to be within the detection limits of the detectors in the next generation such as LCGT and the advanced LIGO for a supernova at 10 kpc. As a contribution to the gravitational wave background, we show that the amplitudes from this source could be larger at the frequency above $\sim$ 1 Hz than the primordial gravitational wave backgrounds, but unfortunately, invisible to the proposed space-based detectors.

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.

The gravitational collapse of a non-rotating, black-hole-forming massive star is studied by neutrino-radiation-hydrodynamical simulations for two different sets of realistic equation of state of dense matter. We show that the event will produce as many neutrinos as the ordinary supernova, but with distinctive characteristics in luminosities and spectra that will be an unmistakable indication of black hole formation. More importantly, the neutrino signals are quite sensitive to the difference of equation of state and can be used as a useful probe into the properties of dense matter. The event will be unique in that they will be shining only by neutrinos (and, possibly, gravitational waves) but not by photons, and hence they should be an important target of neutrino astronomy.

We present the results of numerical experiments, in which we have investigated the influence of the inelastic neutrino-helium interactions on the standing accretion shock instability supposed to occur in the post-bounce supernova core. The axisymmetric hydrodynamical simulations of accretion flows through the standing accretion shock wave onto the protoneutron star show that the interactions are relatively minor and the linear growth of the shock instability is hardly affected. The extra heating given by the inelastic reactions becomes important for the shock revival after the instability enters the non-linear regime, but only when the neutrino luminosity is very close to the critical value, at which the shock would be revived without the interactions. We have also studied the dependence of the results on the initial amplitudes of perturbation and the temperatures of mu and tau neutrinos.

We investigate synthesis of heavy elements in a collapsar. We have calculated detailed composition of magnetically driven jets ejected from a collapsar, which is based on long-term, magneto-hydrodynamic simulations of a rapidly rotating massive star of 40Msun during core collapse. We follow evolution of abundances of about 4000 nuclides from the collapse phase to the ejection phase through the jet generation phase with use of two large nuclear reaction networks. We find that the r-process successfully operates in the jets, so that U and Th are synthesized abundantly when the progenitor has large magnetic field of 10^{12} G and rapidly rotating core. Abundance pattern inside the jets is similar compared to that of r-elements in the solar system. Heavy neutron-rich nuclei \sim 0.01Msun can be ejected from the collapsar. The detailed abundances depend on nuclear properties of mass model, beta-decay rate, and fission, for nuclei near the neutron drip line. Furthermore, we find that p-nuclei are produced without seed nuclei: not only light p-nuclei, such as Se74, Kr78, Sr84, and Mo92, but also heavy p-nuclei, In113, Sn115, and La138, can be abundantly synthesized in the jets. The amounts of p-nuclei in the ejecta are much greater than those in core-collapse supernovae (SNe). In particular, Mo92, In113, Sn115, and La138 deficient in the SNe, are significantly produced in the ejecta.

We perform two-dimensional, axisymmetric, magnetohydrodynamic simulations of the collapse of a rotating star of 40 Msun and in the light of the collapsar model of gamma-ray burst. Considering two distributions of angular momentum, up to \sim 10^{17} cm^2/s, and the uniform vertical magnetic field, we investigate the formation of an accretion disk around a black hole and the jet production near the hole. After material reaches to the black hole with the high angular momentum, the disk is formed inside a surface of weak shock. The disk becomes in a quasi-steady state for stars whose magnetic field is less than 10^{10} G before the collapse. We find that the jet can be driven by the magnetic fields even if the central core does not rotate as rapidly as previously assumed and outer layers of the star has sufficiently high angular momentum. The magnetic fields are chiefly amplified inside the disk due to the compression and the wrapping of the field. The fields inside the disk propagate to the polar region along the inner boundary near the black hole through the Alfv{\'e}n wave, and eventually drive the jet. The quasi-steady disk is not an advection-dominated disk but a neutrino cooling-dominated one. Mass accretion rates in the disks are greater than 0.01 Msun/sec with large fluctuations. The disk is transparent for neutrinos. The dense part of the disk, which locates near the hole, emits neutrino efficiently at a constant rate of < 8 \times 10^{51} erg/s. The neutrino luminosity is much smaller than those from supernovae after the neutrino burst.

We report a current status of our radiation-magnetohydrodynamic code for the study of core-collapse supernovae. In this contribution, we discuss the accuracy of our newly developed numerical code by presenting the test problem in a static background model. We also present the application to the spherically symmetric core-collapse simulations. Since close comparison with the previously published models is made, we are now applying it for the study of magnetorotational core-collapse supernovae.

We present results of our intermediate-band optical imaging survey for high-$z$ Ly$\alpha$ emitters (LAEs) using the prime focus camera, Suprime-Cam, on the 8.2m Subaru Telescope. In our survey, we use eleven filters; four broad-band filters ($B$, $R_{\rm c}$, $i^\prime$, and $z^\prime$) and seven intermediate-band filters covering from 500 nm to 720 nm; we call this imaging program as the Mahoroba-11. The seven intermediate-band filters are selected from the IA filter series that is the Suprime-Cam intermediate-band filter system whose spectral resolution is $R = 23$. Our survey has been made in a $34^\prime \times 27^\prime$ sky area in the Subaru XMM Newton Deep Survey field. We have found 409 IA-excess objects that provide us a large photometric sample of strong emission-line objects. Applying the photometric redshift method to this sample, we obtained a new sample of 198 LAE candidates at $3 < z < 5$. We found that there is no evidence for evolution of the number density and the star formation rate density for LAEs with $\log L({\rm Ly}\alpha) ({\rm erg s^{-1}}) > 42.67$ between $z \sim 3$ and 5.

Pop III stars are the first stars in the universe. They do not contain metals and their formation and evolution may be different from that of stars of later generations. In fact, according to the theory of star formation, Pop III stars might have very massive components ($\sim 100 - 10000M_\odot$). In this paper, we compute the spherically symmetric gravitational collapse of these Pop III massive stars. We solve the general relativistic hydrodynamics and neutrino transfer equations simultaneously, treating neutrino reactions in detail. Unlike supermassive stars ($\gtrsim 10^5 M_\odot$), the stars of concern in this paper become opaque to neutrinos. The collapse is simulated until after an apparent horizon is formed. We confirm that the neutrino transfer plays a crucial role in the dynamics of gravitational collapse, and find also that the $\beta$-equilibration leads to a somewhat unfamiliar evolution of electron fraction. Contrary to the naive expectation, the neutrino spectrum does not become harder for more massive stars. This is mainly because the neutrino cooling is more efficient and the outer core is more massive as the stellar mass increases. Here the outer core is the outer part of the iron core falling supersonically. We also evaluate the flux of relic neutrino from Pop III massive stars. As expected, the detection of these neutrinos is difficult for the currently operating detectors. However, if ever observed, the spectrum will enable us to obtain the information on the formation history of Pop III stars. We investigate 18 models covering the mass range of $300 - 10^4 M_\odot$, making this study the most detailed numerical exploration of spherical gravitational collapse of Pop III massive stars. This will also serve as an important foundation for multi-dimensional investigations.

We discuss effects of magnetic fields on proto-neutron star winds by performing numerical simulation. We assume that the atmosphere of proto-neutron star has a homogenous magnetic field (ranging from ~10^{12} G to ~10^{15} G) perpendicular to the radial direction and examine the dependence of the three key quantities (dynamical time scale, electron fraction, and entropy per baryon) for the successful r-process on the magnetic field strength. Our results show that even with a magneter-class field strength, ~10^{15} G, the feature of the wind dynamics varies only little from that of non-magnetic winds, and that the condition for successful r-process is not realized.

We have numerically studied the instability of the spherically symmetric standing accretion shock wave against non-spherical perturbations. We have in mind the application to the collapse-driven supernovae in the post bounce phase, where the prompt shock wave generated by core bounce is commonly stalled. We take an experimental stand point in this paper. Using spherically symmetric, completely steady, shocked accretion flows as unperturbed states, we have clearly observed both the linear growth and the subsequent nonlinear saturation of the instability. In so doing, we have employed a realistic equation of state together with heating and cooling via neutrino reactions with nucleons. We have done a mode analysis based on the spherical harmonics decomposition and found that the modes with l=1, 2 are dominant not only in the linear regime, but also after the nonlinear couplings generate various modes and the saturation occurs. Varying the neutrino luminosity, we have constructed the unperturbed states both with and without a negative entropy-gradient. We have found that in both cases the growth of the instability is similar, suggesting the convection does not play a dominant role, which also appears to be supported by the recent linear analysis of the convection in accretion flows by Foglizzo et al. The real part of the eigen frequency seems to be mainly determined by the advection time rather than by the sound-crossing time. Whatever the cause may be, the instability is favorable for the shock revival.

We study the evolution of supernova core from the beginning of gravitational collapse of a 15Msolar star up to 1 second after core bounce. We present results of spherically symmetric simulations of core-collapse supernovae by solving general relativistic neutrino-radiation-hydrodynamics in the implicit time-differencing. We aim to explore the evolution of shock wave in a long term and investigate the formation of protoneutron star together with supernova neutrino signatures. These studies are done to examine the influence of equation of state (EOS) on the postbounce evolution of shock wave in the late phase and the resulting thermal evolution of protoneutron star. We make a comparison of two sets of EOS, that is, by Lattimer and Swesty (LS-EOS) and by Shen et al.(SH-EOS). We found that, for both EOSs, the core does not explode and the shock wave stalls similarly in the first 100 milliseconds after bounce. The revival of shock wave does not occur even after a long period in either cases. However, the recession of shock wave appears different beyond 200 milliseconds after bounce, having different thermal evolution of central core. A more compact protoneutron star is found for LS-EOS than SH-EOS with a difference in the central density by a factor of ~2 and a difference of ~10 MeV in the peak temperature. Resulting spectra of supernova neutrinos are different to the extent that may be detectable by terrestrial neutrino detectors.

We study the explosion mechanism of collapse-driven supernovae by numerical simulations with a new nuclear EOS based on unstable nuclei. We report new results of simulations of general relativistic hydrodynamics together with the Boltzmann neutrino-transport in spherical symmetry. We adopt the new data set of relativistic EOS and the conventional set of EOS (Lattimer-Swesty EOS) to examine the influence on dynamics of core-collapse, bounce and shock propagation. We follow the behavior of stalled shock more than 500 ms after the bounce and compare the evolutions of supernova core.

We perform two-dimensional numerical simulations on the core-collapse of a massive star with strong magnetic fields and differential rotations using a numerical code ZEUS-2D. Changing field configurations and laws of differential rotation parametrically, we compute 14 models and investigate effects of these parameters on the dynamics. In our models, we do not solve the neutrino transport and instead employ a phenomenological parametric EOS that takes into account the neutrino emissions. As a result of the calculations, we find that the field configuration plays a significant role in the dynamics of the core if the initial magnetic field is large enough. Models with initially concentrated fields produce more energetic explosions and more prolate shock waves than the uniform field. Quadrapole-like fields produce remarkably collimated and fast jet, which might be important for gamma-ray bursts(GRB). The Lorentz forces exerted in the region where the plasma-beta is less than unity are responsible for these dynamics. The pure toroidal field, on the other hand, does not lead to any explosion or matter ejection. This suggests the presupernova models of Heger et al.(2003), in which toroidal fields are predominant, is disadvantageous for the magnetorotation-induced supernova considered here. Models with initially weak magnetic fields do not lead to explosion or matter ejection, either. In these models magnetic fields play no role as they do not grow on the timescale considered in this paper so that the magnetic pressure could be comparable to the matter pressure. This is because the exponential field growth as expected in MRI is not seen in our models. The magnetic field is amplified mainly by field-compression and field-wrapping in our simulations.

We investigate $r$-process nucleosynthesis during the magnetohydrodynamical (MHD) explosion of supernova in a massive star of 13 $M_{\odot}$. Contrary to the case of the spherical explosion, jet-like explosion due to the combined effects of the rotation and magnetic field lowers the electron fraction significantly inside the layers above the iron core. We find that the ejected material of low electron fraction responsible for the $r$-process comes out from the silicon rich layer of the presupernova model. This leads to the production up to the third peak in the solar $r$-process elements. We examine whether the fission affects the $r$-process paths by using the full nuclear reaction network with both the spontaneous and $\beta$-delayed fission included. Moreover, we pay particular attention how the mass formula affects the $r$-process peaks with use of two mass formulae. It is found that both formulae can reproduce the global abundance pattern up to the third peak though detailed distributions are rather different. We point out that there are variations in the $r$-process nucleosynthesis if the MHD effects play an important role in the supernova explosion.

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.

We study the hydrodynamical effects of two colliding shells, adopted to model internal shocks in various relativistic outflows such as gamma-ray bursts and blazars. We find that the density profiles are significantly affected by the propagation of rarefaction waves. A split-feature appears at the contact discontinuity of the two shells. The shell spreading with a few ten percent of the speed of light is also shown to be a notable aspect. The conversion efficiency of the bulk kinetic energy to internal one shows deviations from the widely-used inelastic two-point-mass-collision model. Observational implications are also shortly discussed.

We perform a series of two-dimensional magnetohydrodynamic simulations of supernova cores. Since the distributions of the angular momentum and the magnetic fields of strongly magnetized stars are quite uncertain, we systematically change the combinations of the strength of the angular momentum, the rotations law, the degree of differential rotation, and the profiles of the magnetic fields to construct the initial conditions. By so doing, we estimate how the rotation-induced anisotropic neutrino heating are affected by the strong magnetic fields through parity-violating effects and first investigate how the north-south asymmetry of the neutrino heating in a strongly magnetized supernova core could be. As for the microphysics, we employ a realistic equation of state based on the relativistic mean field theory and take into account electron captures and the neutrino transport via the neutrino leakage scheme. With these computations, we find that the parity-violating corrections reduce $ \lesssim 0.5 %$ of the neutrino heating rate than that without the magnetic fields in the vicinity of the north pole of a star, on the other hand, enhance about $ \lesssim 0.5 %$ in the vicinity of the south pole. If the global asymmetry of the neutrino heating in the both of the poles develops in the later phases, the newly born neutron star might be kicked toward the north pole in the subsequent time.

We present the results of indirect searches for Weakly Interacting Massive Particles (WIMPs) with 1679.6 live days of data from the Super-Kamiokande detector using neutrino-induced upward through-going muons. The search is performed by looking for an excess of high energy muon neutrinos from WIMP annihilations in the Sun, the core of the Earth, and the Galactic Center, as compared to the number expected from the atmospheric neutrino background. No statistically significant excess was seen. We calculate flux limits in various angular cones around each of the above celestial objects. We obtain conservative model-independent upper limits on WIMP-nucleon cross-section as a function of WIMP mass and compare these results with the corresponding results from direct dark matter detection experiments.

We present the results of a survey for Ly-alpha emitters at z ~ 5.8 using a new intermediate-band filter centered at lambda_c = 8275 AA with Delta-lambda_FWHM = 340 AA (i.e., the spectroscopic resolution is R ~ 23) with a combination with a traditional narrow-band centered at lambda_c = 8150 AA with Delta-lambda_FWHM = 120 AA (R ~ 68). Our observations were made with use of the Subaru Prime Focus Camera, Suprime-Cam, on the 8.2 m Subaru telescope in a sky area surrounding the high redshift quasar, SDSSp J104433.04-012522.2 at z=5.74, covering an effective sky area with ~ 720 arcmin^2. In this survey, we have found four Ly-alpha-emitter candidates from the intermediate-band image (z \~ 5.8 with Delta z ~ 0.3). Combined with our previous results based on the NB816 imaging, we discuss the star formation activity in galaxies between z ~ 5.7 and z ~ 5.9.

We study both analytically and numerically hydrodynamical effects of two colliding shells, the simplified models of the internal shock in various relativistic outflows such as gamma-ray bursts and blazars. We pay particular attention to three interesting cases: a pair of shells with the same rest mass density (``{\it equal rest mass density}''), a pair of shells with the same rest mass (``{\it equal mass}''), and a pair of shells with the same bulk kinetic energy (``{\it equal energy}'') measured in the intersteller medium (ISM) frame. We find that the density profiles are significantly affected by the propagation of rarefaction waves. A split-feature appears at the contact discontinuity of two shells for the ``equal mass'' case, while no significant split appears for the ``equal energy'' and ``equal rest mass density'' cases. The shell spreading with a few ten percent of the speed of light is also shown as a notable aspect caused by rarefaction waves. The conversion efficiency of the bulk kinetic energy to internal one is numerically evaluated. The time evolutions of the efficiency show deviations from the widely-used inellastic two-point-mass-collision model.

We have performed 2.5-dimensional general relativistic magnetohydrodynamic (MHD) simulations of collapsars including a rotating black hole. This paper is an extension of our previous paper (Mizuno et al. 2004). The current calculation focuses on the effect of black hole rotation using general relativistic MHD with simplified microphysics. Initially, we assume that the core collapse is failed in this star. A few solar mass rotating black hole is inserted by hand into the calculation. We consider two cases, the co-rotating case and counter-rotating case with respect to the black hole rotation. Although the counter-rotating case may be unrealistic for collapsar, we perform as a maximally dragging case of magnetic field. The simulation results show the formation of a disk-like structure and the generation of a jet-like outflow near the central black hole. The jet-like outflow propagates outwardly with the twisted magnetic field and becomes collimated. We have found that the jets are generated and accelerated mainly by the magnetic field. The total jet velocity in the rotating black hole case is comparable to that of the non-rotating black hole case (Mizuno et al. 2004), 0.3c. When the rotation of the black hole is faster, the magnetic field is twisted strongly owing to the frame-dragging effect. The magnetic energy stored by the twisting magnetic field is converted to kinetic energy of the jet directly rather than propagating as an Alfven wave. Thus, as the rotation of the black hole becomes faster, the poloidal velocity of the jet becomes faster. In the rapidly rotating black hole case the jet-like outflow can be produced by the frame dragging effect only through the twisting of magnetic field even if there is no stellar rotation.

We have performed 2.5-dimensional general relativistic magnetohydrodynamic (MHD) simulations of the gravitational collapse of a magnetized rotating massive star as a model of gamma ray bursts (GRBs). The current calculation focuses on general relativistic MHD with simplified microphysics (we ignore neutrino cooling, physical equation of state, and photodisintegration). Initially, we assume that the core collapse has failed in this star. A few $M_{\odot}$ black hole is inserted by hand into the calculation. The simulations presented in the paper follow the accretion of gas into a black hole that is assumed to have formed before the calculation begins.The simulation results show the formation of a disk-like structure and the generation of a jetlike outflow inside the shock wave launched at the core bounce. We have found that the jet is accelerated by the magnetic pressure and the centrifugal force and is collimated by the pinching force of the toroidal magnetic field amplified by the rotation and the effect of geometry of the poloidal magnetic field. The maximum velocity of the jet is mildly relativistic (0.3c). The velocity of the jet becomes larger as the initial rotational velocity of stellar matter gets faster. On the other hand, the dependence on the initial magnetic field strength is a bit more complicated: the velocity of the jet increases with the initial field strength in the weak field regime, then is saturated at some intermediate field strength, and decreases beyond the critical field strength. These results are related to the stored magnetic energy determined by the balance between the propagation time of the Alfven wave and the rotation time of the disk (or twisting time).

We investigate the dynamics and morphology of jets propagating into the interstellar medium using 2D relativistic hydrodynamic simulations. The calculations are performed assuming axisymmetric geometry and trace a long propagation of jets. The jets are assumed to be 'light' with the density ratio between the beam to the ambient gas much less than unity. We examine the mechanism for the appearance of vortices at the head of jets in the hot spot. Such vortices are known as a trigger of a deceleration phase which appears after a short phase in which one dimensional analysis is dominant. We find that an oblique shock at the boundary rim near the end of the beam strongly affects the flow structure in and around the hot spot. Weakly shocked gas passed through this oblique shock and becomes a trigger for the generation of vortices. We also find the parameter dependence of these effects for the propagation and dynamics of the jets. The jet with slower propagation velocity is weakly pinched, has large vortices and shows very complex structure at the head of the jets and extended synchrotron emissivity.

We perform a series of two-dimensional, axisymmetric, magnetohydrodynamic simulations of the rotational collapse of a supernova core. In order to calculate the waveforms of the gravitational wave, we derive the quadrupole formula including the contributions from the electromagnetic fields. Recent stellar evolution calculations imply that the magnetic fields of the toroidal components are much stronger than those of the poloidal ones at the presupernova stage. Thus, we systematically investigate the effects of the toroidal magnetic fields on the amplitudes and waveforms. Furthermore, we employ the two kinds of the realistic equation of states, which are often used in the supernova simulations. Then, we investigate the effects of the equation of states on the gravitational wave signals. With these computations, we find that the peak amplitudes are lowered by an order of 10% for the models with the strongest toroidal magnetic fields. However, the peak amplitudes are mostly within sensitivity range of laser interferometers such as TAMA and the first LIGO for a source at a distance of 10 kpc. Furthermore, we point out that the amplitudes of second peaks are still within the detection limit of the first LIGO for the source, although the characteristics of second peaks are reduced by the magnetic fields. We stress the importance of the detection, since it will give us information about the angular momentum distribution of massive evolved stars. When we compare the gravitational waves from the two realistic equation of states, significant differences are not found, except that the typical frequencies of the gravitational wave become slightly higher for the softer equation of state.

We have performed 2.5-dimensional general relativistic magnetohydrodynamic (MHD) simulations of the gravitational collapse of a magnetized rotating massive star as a model of gamma ray bursts (GRBs). This simulation showed the formation of a disk-like structure and the generation of a jet-like outflow inside the shock wave launched at the core bounce. We have found the jet is accelerated by the magnetic pressure and the centrifugal force and is collimated by the pinching force of the toroidal magnetic field amplified by the rotation and the effect of geometry of the poloidal magnetic field. The maximum velocity of the jet is mildly relativistic ($\sim$ 0.3 c).

We study characteristics of the relativistic equation of state (EOS) for collapse-driven supernovae, which is derived by relativistic nuclear many body theory. Recently the relativistic EOS table has become available as a new complete set of physical EOS for numerical simulations of supernova explosion. We examine this EOS table by using general relativistic hydrodynamics of the gravitational collapse and bounce of supernova cores. In order to study dense matter in dynamical situation, we perform simplified calculations of core collapse and bounce by following adiabatic collapse with the fixed electron fraction for a series of progenitor models. This is intended to give us ``approximate models'' of prompt explosion. We investigate the profiles of thermodynamical quantities and the compositions during collapse and bounce. We also perform the calculations with the Lattimer-Swesty EOS to compare the properties of dense matter. As a measure of the stiffness of the EOS, we examine the explosion energy of the prompt explosion with electron capture totally suppressed. We study the derivative of the thermodynamical quantities obtained by the relativistic EOS to discuss the convective condition in neutron-rich environment, which may be important in the delayed explosion.

We present an optical multicolor-imaging study of the galaxy cluster Abell 521 at $z = 0.25$, using Suprime-Cam on the Subaru Telescope, covering an area of $32 \times 20$ arcmin$^2$ ($9.4 \times 5.8 h_{50}^{-2}$ Mpc$^2$ at $z = 0.25$). Our imaging data taken with both a narrow-band filter, $NB816$ ($\lambda_0 = 8150$\AA and $\Delta \lambda = 120$\AA), and broad-band filters, $B,V,R_{\rm C}, i^\prime$, and $z^\prime$ allow us to find 165 H$\alpha$ emitters. We obtain the H$\alpha$ luminosity function (LF) for the cluster galaxies within 2 Mpc; the Schechter parameters are $\alpha = -0.75 \pm 0.23$, $\phi^\star = 10^{-0.25 \pm 0.20}$ Mpc$^{-3}$, and $L^\star = 10^{42.03 \pm 0.17}$ erg s$^{-1}$. Although the faint end slope, $\alpha$, is consistent with that of the local cluster H$\alpha$ LFs, the characteristic luminosity, $L^\star$, is about 6 times (or $\approx 2$ mag) brighter. This strong evolution implies that Abell 521 contains more active star-forming galaxies than the local clusters, being consistent with the observed Butcher-Oemler effect. However, the bright $L^\star$ of Abell 521 may be, at least in part, due to the dynamical condition of this cluster.

We present the results of a survey for Ly$\alpha$ emitters at $z\approx 5.7$ based on optical narrow-band ($\lambda_{\rm c} = 8150$ \AA ~ and $\Delta\lambda = 120$ \AA), and broad-band ($B$, $R_{\rm C}$, $I_{\rm C}$, and $z^\prime$) observations of the field surrounding the high redshift quasar, SDSSp J104433.04$-$012522.2, on the 8.2 m Subaru Telescope with the Subaru Prime Focus Camera, Suprime-Cam. This survey covers a sky area of $\approx 720$ arcmin$^2$ and a co-moving volume of $\simeq 2 \times 10^5 ~ h_{0.7}^{-3}$ Mpc$^3$. We have found 20 candidates of Ly$\alpha$ emitters at $z \approx$ 5.7 with $\Delta z \approx 0.1$. Two of them have been confirmed star-forming galaxies at $z=5.655$ and $z=5.687$ from our follow-up optical spectroscopy. We discuss star-formation properties of the 20 objects from a statistical point of view. Our survey leads to a new estimate of the star formation rate density at $z \approx 5.7$, $\sim 1.2 \times 10^{-3} h_{0.7} M_\odot$ yr$^{-1}$ Mpc$^{-3}$.

Several high-z (z > 5.7) quasars have been found in the course of Sloan Digital Sky Survey. The presence of such very high-z quasars is expected to give constraints on early structure formation. On one hand, it is suggested that these most luminous objects at high redshift are biased toward the highly magnified objects by gravitational lensing. To clarify the effect of gravitational lensing on the high-z quasars, we began the imaging survey of intervening lensing galaxies. Indeed our previous optical image showed that SDSSp J104433.04+012502.2 at z=5.74 is gravitationally magnified by a factor 2. In this paper, we report our new optical imaging of other two high-z quasars, SDSSp J103027.10+052455.0 at z=6.28 and SDSSp J130608.26+035626.3 at z=5.99. Since we find neither intervening galaxy nor counter image with i^{\prime} < 25.4-25.8 around each quasar, we conclude that they are not strongly magnified regardless that a lens galaxy is dusty.

We have done a series of two-dimensional hydrodynamic simulations of the rotational collapse of a supernova core in axisymmetry. We have employed a realistic equation of state (EOS) and taken into account electron captures and neutrino transport by the so-called leakage scheme. It is an important progress to apply the realistic EOS coupled with the microphysics to 2-D simulations for computing gravitational radiation in rotational core collapse. We have used the quadrupole formula to calculate the amplitudes and the waveforms of gravitational wave assuming the Newtonian gravity. From these computations, we have extended the conventional category of the gravitational waveforms. Our results have shown that the peak amplitudes of gravitational wave are mostly within the sensitivity range of the laser interferometers such as TAMA and first LIGO for a source at a distance of 10 kpc. Furthermore we have found that the amplitudes of the second peaks are within the detection limit of first LIGO for the source and first pointed out the importance of the detections, since they will give us the information as to the angular momentum distribution of evolved massive stars.

We perform 2-dimensional relativistic hydrodynamical simulations in the context of collapsar model. Calculations of explosive nucleosynthesis are also accomplished. We investigate the influence of the structure of the progenitor and energy deposition rate on the resulting explosive nucleosynthesis, assuming that 56Ni is mainly synthesized in the jet launched by the neutrino heating. We show the amount of 56Ni is very sensitive to the energy deposition rate. Thus we conclude that it is quite natural not to detect an underlying supernova in some X-ray afterglows as in GRB 010921. We also point out the possibility that the relative abundance of the elements with intermediate mass number such as Si and S in the X-ray afterglow of GRB 011211 may be naturally explained if the energy deposition rate at the central engine is relatively long because little amount of 56Ni should be synthesized under such an environment. If this discussion is true, there should be correlation between the line features in the X-ray afterglow and duration of the GRB. It should be noted that the duration of GRB 011211 is 270 seconds, making it the longest burst ever observed by Beppo-SAX although it suffers from the effect of red-shift (z_host=2.14), which supports our conclusion. Our results also suggest that the type I collapsar model in which the energy deposition rate is relatively low (\dot{E} \sim 10^{51} erg/s) might have difficulty in reproducing the observed amount of 56Ni in a hypernova such as SN 1998bw. This means that the mechanism of the central engine of a hypernova accompanying GRB may be constrained by the discussion of explosive nucleosynthesis.

We have carried out a deep imaging survey for H$\alpha$ emitting galaxies at $z\approx$0.24 using a narrowband filter tuned with the redshifted line. The total sky area covered is 706 arcmin$^2$ within a redshift range from 0.234 to 0.252 ($\delta z$=0.018). This corresponds to a volume of 3.9$\times10^3$ Mpc$^3$ when $\Omega_{\rm matter}=0.3$, $\Omega_\Lambda=0.7$, and $H_{\mathrm{0}}$=70 km s$^{-1}$ Mpc$^{-1}$ are adopted. We obtain a sample of 348 H$\alpha$ emitting galaxies whose observed emission-line equivalent widths are greater than 12 \AA. We find an extinction-corrected H$\alpha$ luminosity density of $10^{39.65^{+0.08}_{-0.12}}$ ergs s$^{-1}$ Mpc$^{-3}$. Using the Kennicutt relation between the H$\alpha$ luminosity and star formation rate, the star formation rate density in the covered volume is estimated as $0.036^{+0.006}_{-0.012}$ $M_{\sun}$ yr$^{-1}$ Mpc$^{-3}$. This value is higher by a factor of 3 than the local SFR density.

We study the effects of the outer boundary conditions in neutrino-driven winds on the r-process nucleosynthesis. We perform numerical simulations of hydrodynamics of neutrino-driven winds and nuclear reaction network calculations of the r-process. As an outer boundary condition of hydrodynamic calculations, we set a pressure upon the outermost layer of the wind, which is approaching toward the shock wall. Varying the boundary pressure, we obtain various asymptotic thermal temperature of expanding material in the neutrino-driven winds for resulting nucleosynthesis. We find that the asymptotic temperature slightly lower than those used in the previous studies of the neutrino-driven winds can lead to a successful r-process abundance pattern, which is in a reasonable agreement with the solar system r-process abundance pattern even for the typical proto-neutron star mass Mns ~ 1.4 Msun. A slightly lower asymptotic temperature reduces the charged particle reaction rates and the resulting amount of seed elements and lead to a high neutron-to-seed ratio for successful r-process. This is a new idea which is different from the previous models of neutrino-driven winds from very massive (Mns ~ 2.0 Msun) and compact (Rns ~ 10 km) neutron star to get a short expansion time and a high entropy for a successful r-process abundance pattern. Although such a large mass is sometimes criticized from observational facts on a neutron star mass, we dissolve this criticism by reconsidering the boundary condition of the wind. We also explore the relation between the boundary condition and neutron star mass, which is related to the progenitor mass, for successful r-process.

In this paper we present a numerical study of plasma jets produced by intense laser matter interactions. Through this study we hope to better understand astrophysical jets and their recent experimental simulations in the laboratory. We paid special attention to radiation cooling and the interaction of the jet with ambient gas. Four cases are presented in this paper; two of them deal with the propagation of jets in vacuum, while in the other two the propagation takes place in the ambient gas. Available experimental results are reproduced to good accuracy in the vacuum case. For jets in ambient gas, we find that the existence of the surrounding gas confines the jet into a narrow cylindrical shape so that both the density and temperature of the jet remain high enough for effective radiation cooling. As a result, a collimated plasma jet is formed in these cases. The dimensionless parameters characterizing the laboratory jets and protostellar jets have overlapping domains. We also discuss the cooling lengths for our model and compare them with the corresponding values in the astrophysical jets. A plasma jet in the ambient gas experiment is proposed which is within the reach of present day technology, and can be relevant to astrophysical phenomena.

We reanalyze $r$-process nucleosynthesis in the neutron-rich ejecta from a prompt supernova explosion of a low-mass (11 M$_\odot$) progenitor. A pompt explosion is not yet ruled out as a possibility for low-mass supernova progenitors. Moreover, there is mounting evidence that a new $r$-process site may be required. Hence, we assume that a prompt explosion can occur and make a study of r-process nucleosynthesis in the supernova ejecta. To achieve a prompt explosion we have performed a general relativistic hydrodynamic simulation of adiabatic collapse and bounce using a relativistic nuclear-matter equation of state. The electron fraction $Y_e$ during the collapse was fixed at the initial-model value. The size of the inner collapsing core was then large enough to enable a prompt explosion to occur. Adopting the calculated trajectories of promptly ejected material, we explicitly computed the burst of neutronization due to electron captures on free protons in the photodissociated ejecta after the passage of the shock. The thermal and compositional evolution of the resulting neutron-rich ejecta originating from near the surface of proto-neutron star was obtained. These were used in nuclear reaction network calculations to evaluate the products of $r$-process nucleosynthesis. We find that, unlike earlier studies, the amount of $r$-process material ejected per supernova is quite consistent with observed galactic $r$-process abundances. Furthermore, the computed $r$-process abundances are in good agreement with Solar abundances of $r$-process elements for A$> 100$. This suggests that prompt supernovae are still a viable $r$-process site. Such events may be responsible for the abundances of the heaviest $r$-process nuclei.

We study the equation of state of electron in strong magnetic fields which are greater than the critical value $B_c \simeq 4.4 \times 10^{13}$ Gauss. We find that such a strong magnetic field induces the anisotropic pressure of electron. We apply the result to the neutrino-driven wind in core-collapse supernovae and find that it can produce large entropy per baryon, $S \sim 400 k_B$. This mechanism might successfully account for the production of the heavy nuclei with mass numbers A = 80 -- 250 through the r-process nucleosynthesis.

The vacuum polarization tensor in strong external magnetic fields has been evaluated numerically for various strengths of magnetic fields and momenta of photons under the threshold of the $e^{\pm}$ pair creation. The fitting formula has been obtained which reproduces the calculated results within 10 % of error. The proper time method is employed further to obtain the retarded polarization tensor for finite temperature plasmas.

With a view of applications to the simulations of supernova explosion and proto neutron star cooling, we derive the Boltzmann equations for the neutrino transport with the flavor mixing based on the real time formalism of the nonequilibrium field theory and the gradient expansion of the Green function. The relativistic kinematics is properly taken into account. The advection terms are derived in the mean field approximation for the neutrino self-energy whiles the collision terms are obtained in the Born approximation. The resulting equations take the familiar form of the Boltzmann equation with corrections due to the mixing both in the advection part and in the collision part. These corrections are essentially the same as those derived by Sirera et al. for the advection terms and those by Raffelt et al. for the collision terms, respectively, though the formalism employed here is different from theirs. The derived equations will be easily implemented in numerical codes employed in the simulations of supernova explosions and proto neutron star cooling.

We study the neutrino-driven wind from the proto-neutron star by the general relativistic hydrodynamical simulations. We examine the properties of the neutrino-driven wind to explore the possibility of the r-process nucleosynthesis. The numerical simulations with the neutrino heating and cooling processes are performed with the assumption of the constant neutrino luminosity by using realistic profiles of the proto-neutron star (PNS) as well as simplified models. The dependence on the mass of PNS and the neutrino luminosity is studied systematically. Comparisons with the analytic treatment in the previous studies are also done. In the cases with the realistic PNS, we found that the entropy per baryon and the expansion time scale are neither high nor short enough for the r-process within the current assumptions. On the other hand, we found that the expansion time scale obtained by the hydrodynamical simulations is systematically shorter than that in the analytic solutions due to our proper treatment of the equation of state. This fact might lead to the increase of the neutron-to-seed ratio, which is suitable for the r-process in the neutrino-driven wind. Indeed, in the case of massive and compact proto-neutron stars with high neutrino luminosities, the expansion time scale is found short enough in the hydrodynamical simulations and the r-process elements up to A ~ 200 are produced in the r-process network calculation.

We studied possible modifications of the neutrino - nucleon scattering rate due to the nucleon - nucleon collisions in the hot dense matter which we find in the supernova core. We show that the finite width of the nucleon spectral function induced by the nucleon collisions leads to broadening of the dynamical spin structure function of the nucleon, resulting in the reduction of the rate of neutrino - nucleon scattering via the axial vector current and making the energy exchange between neutrinos and nucleons easier. The reduction rate is relatively large (about 0.6) even at density of about 10^{13}g/cm^{3} and could have a significant impact on the dynamics of the collapse-driven supernova as well as the cooling of the proto neutron star.

In view of the application to supernova simulations, we calculate neutrino reaction rates with nucleons via the neutral and charged currents in the supernova core in the relativistic random phase approximation (RPA) and study their effects on the opacity of the supernova core. The formulation is based on the Lagrangian employed in the calculation of nuclear equation of state (EOS) in the relativistic mean field theory (RMF). The nonlinear meson terms are treated appropriately so that the consistency of the density correlation derived in RPA with the thermodynamic derivative obtained from EOS by RMF is satisfied in the static and long wave length limit. We employ pion and rho meson exchange interactions together with the phenomenological Landau-Migdal parameters for the isospin-dependent nuclear interactions. We find that both the charged and neutral current reaction rates are suppressed from the standard Bruenn's approximate formula considerably in the high density regime. In the low density regime, on the other hand, the vector current contribution to the neutrino-nucleon scattering rate is enhanced in the vicinity of the boundary of the liquid-gas phase transition, while the other contributions are moderately suppressed there also. In the high temperature regime or in the regime where electrons have a large chemical potential, the latter of which is important only for the electron capture process and its inverse process, the recoil of nucleons cannot be neglected and further reduces the reaction rates with respect to the standard approximate formula which discards any energy transfer in the processes. These issues could have a great impact on the neutrino heating mechanism of collapse-driven supernovae.

We have coded a Boltzmann solver based on a finite difference scheme (S_N method) aiming at calculations of neutrino transport in type II supernovae. Close comparison between the Boltzmann solver and a Monte Carlo transport code has been made for realistic atmospheres of post bounce core models under the assumption of a static background. We have also investigated in detail the dependence of the results on the numbers of radial, angular, and energy grid points and the way to discretize the spatial advection term which is used in the Boltzmann solver. A general relativistic calculation has been done for one of the models. We find overall good agreement between the two methods. However, because of a relatively small number of angular grid points (which is inevitable due to limitations of the computation time) the Boltzmann solver tends to underestimate the flux factor and the Eddington factor outside the (mean) ``neutrinosphere'' where the angular distribution of the neutrinos becomes highly anisotropic. This fact suggests that one has to be cautious in applying the Boltzmann solver to a calculation of the neutrino heating in the hot-bubble region because it might tend to overestimate the local energy deposition rate. A comparison shows that this trend is opposite to the results obtained with a multi-group flux-limited diffusion approximation of neutrino transport. The accuracy of the Boltzmann solver can be considerably improved by using a variable angular mesh to increase the angular resolution in the semi-transparent regime.

The large abundance ratio of $^{44}Ti/^{56}Ni$ in Cas A is puzzling. In fact, the ratio seems to be larger than the theoretical constraint derived by Woosley & Hoffman (1991). However, this constraint is obtained on the assumption that the explosion is spherically symmetric, whereas Cas A is famous for the asymmetric form of the remnant. Recently, Nagataki et al. (1997) calculated the explosive nucleosynthesis of axisymmetrically deformed collapse-driven supernova. They reported that the ratio of $^{44}Ti/^{56}Ni$ was enhanced by the stronger alpha-rich freezeout in the polar region. In this paper, we apply these results to Cas A and examine whether this effect can explain the large amount of $^{44}Ti$ and the large ratio of $^{44}Ti/^{56}Ni$. We demonstrate that the conventional spherically symmetric explosion model can not explain the $^{44}$Ti mass produced in Cas A if its lifetime is shorter than $\sim$ 80 years and the intervening space is transparent to the gamma-ray line from the decay of $^{44}$Ti. On the other hand, we show the axisymmetric explosion models can solve the problem. We expect the same effect from a three dimensionally asymmetric explosion, since the stronger alpha-rich freezeout will also occur in that case in the region where the larger energy is deposited.

We have performed 1-dimensional calculations for explosive nucleosynthesis in collapse-driven supernova and investigated its sensitivity to the initial form of the shock wave. We have found the tendency that the peak temperature becomes higher around the mass cut if the input energy is injected more in the form of kinetic energy rather than internal energy. Then, the mass cut becomes larger, and, as a result, neutron-rich matter is less included in the ejecta; this is favorable for producing the observational data compared with a previous model. Our results imply that the standard method to treat various processes for stellar evolution, such as convection and electron capture during the silicon burning stage, are still compatible with the calculation of explosive nucleosynthesis.

Explosive nucleosynthesis under the axisymmetric explosion in Type II supernova has been performed by means of two dimensional hydrodynamical calculations. We have compared the results with the observations of SN 1987A. Our chief findings are as follows: (1) $^{44}Ti$ is synthesized so much as to explain the tail of the bolometric light curve of SN 1987A. We think this is because the alpha-rich freezeout takes place more actively under the axisymmetric explosion. (2) $^{57}Ni$ and $^{58}Ni$ tend to be overproduced compared with the observations. However, this tendency relies strongly on the progenitor's model. We have also compared the abundance of each element in the mass number range $A= 16-73$ with the solar values. We have found three outstanding features. (1) For the nuclei in the range $A=16-40$, their abundances are insensitive to the initial form of the shock wave. This insensitivity is favored since the spherical calculations thus far can explain the solar system abundances in this mass range. (2) There is an enhancement around A=45 in the axisymmetric explosion compared with the spherical explosion fairly well. In particular, $^{44}Ca$, which is underproduced in the present spherical calculations, is enhanced significantly. (3) In addition, there is an enhancement around A=65. This tendency does not rely on the form of the mass cut but of the initial shock wave. This enhancement may be the problem of the overproduction in this mass range, although this effect would be relatively small since Type I supernovae are chiefly responsible for this mass number range.

First results of numerical simulations are presented which compute the dynamical evolution of a neutron star with a mass slightly below the minimum stable mass by means of a new implicit (general relativistic) hydrodynamic code. We show that such a star first undergoes a phase of quasi-static expansion, caused by slow nuclear $\beta$-decays, lasting for about 20 seconds, but then explodes violently. The kinetic energy of the explosion is around $10^{49} erg$, the peak luminosity in electron anti-neutrinos is of order $10^{52} erg/s$, and the thermodynamic conditions of the expanding matter are favorable for r-process nucleosynthesis. These results are obtained for the Harrison-Wheeler equation of state and a simple and, possibly, unrealistic treatment of $\beta$-decay rates and nuclear fission, which were adopted for comparison with previous works. However, we do not expect that the outcome will change qualitatively if more recent nuclear input physics used. Although our study does not rely on a specific scenario ofhow a neutron star starting from a bigger (and stable) masscan reach the dynamical phase, we assume that the final mass-loss event happens on a very short time scale, i.e., on a time scale shorter than a sound-crossing time, by removing a certain amount of mass as an initial perturbation. This assumption implies that the star has no time to adjust its nuclear composition to the new mass mass through a sequence of quasi-equilibria. In the latter case, however, there exists no stable configuration below the minimum mass, because the equation of state of fully catalyzed matter is too soft. Therefore, the dynamics of the explosion will not be too different from what have obtained if different initial perturbations are assumed.

An implicit Lagrangian hydrodynamics code for general relativistic spherical collapse is presented. This scheme is based on an approximate linearized Riemann solver (Roe type scheme). This code is aimed especially at the calculation of the late phase of collapse-driven supernovae and the nascent neutron star, where there is a remarkable contrast between the dynamical time scale of the proto-neutron star and the diffusion time scale of neutrinos, without such severe limitation of the Courant condition at the center of the neutron star. Several standard test calculations have been done. Two other adiabatic simulations have also been done in order to test the performance of the code in the context of the collapse-driven supernovae. It is found that the time step can be extended far beyond the Courant limitation at the center of the neutron star. The details of the scheme and the results of these test calculations are discussed.