We present XMM-Newton observations of N132D, the X-ray brightest supernova remnant (SNR) in the Large Magellanic Cloud (LMC), using the Reflection Grating Spectrometer (RGS) that enables high-resolution spectroscopy in the soft X-ray band. A dozen emission lines from L-shell transitions of various elements at intermediate charge states are newly detected in the RGS data integrating the ~200-ks on-axis observations. The 0.3-2.0-keV spectra require at least three components of thermal plasmas with different electron temperatures and indicate clear evidence of non-equilibrium ionization (NEI). Our detailed spectral diagnostics further reveal that the forbidden-to-resonance line ratios of O VII and Ne IX are both higher than expected for typical NEI plasmas. This enhancement could be attributed to either resonance scattering or emission induced by charge exchange in addition to a possible contribution from the superposition of multiple temperature components, although the lack of spatial information prevents us from concluding which is most likely.

We propose a concept of multiplexing lobster-eye (MuLE) optics to achieve significant reductions in the number of focal plane imagers in lobster-eye (LE) wide-field X-ray monitors. In the MuLE configuration, an LE mirror is divided into several segments and the X-rays reflected on each of these segments are focused on a single image sensor in a multiplexed configuration. If each LE segment assumes a different rotation angle, the azimuthal rotation angle of a cross-like image reconstructed from a point source by the LE optics identifies the specific segment that focuses the X-rays on the imager. With a focal length of 30 cm and LE segments with areas of 10 x 10 cm^2, ~1 sr of the sky can be covered with 36 LE segments and only four imagers (with total areas of 10 x 10 cm^2). A ray tracing simulation was performed to evaluate the nine-segment MuLE configuration. The simulation showed that the flux (0.5 to 2 keV) associated with the 5-sigma detection limit was ~2 x 10^-10 erg cm^-2 s^-1 (10 mCrab) for a transient with a duration of 100 s. The simulation also showed that the direction of the transient for flux in the range of 14 to 17 mCrab at 0.6 keV was determined correctly with 99.7% confidence limit. We conclude that the MuLE configuration can become an effective on-board device for small satellites for future X-ray wide-field transient monitoring.

In this decade, GeV/TeV gamma-ray observations of several supernova remnants (SNRs) have implied that accelerated particles are escaping from their acceleration sites. However, when and how they escape from the SNR vicinities are yet to be understood. Recent studies have suggested that the particle escape might develop with thermal plasma ages of the SNRs. In this paper, we present a systematic study on time evolution of particle escape using thermal X-ray properties and gamma-ray spectra. We used 38 SNRs which associate with GeV/TeV gamma-ray emissions. We conducted spectral fittings on the gamma-ray spectra using exponential cutoff power law and broken power law models to estimate the exponential cutoff or the break energies, both of which are indicators of particle escape. The plots of the gamma-ray cutoff/break energies over the plasma ages show similar tendencies to those predicted by simple theories of the particle escape under conditions in which a shock is interacting with thin interstellar medium or clouds. The particle escape timescale is estimated as $\sim$100 kyr from decreasing trends of the total energy of the confined protons with the plasma age. The large dispersions of the cutoff/break energies of the data may suggest an intrinsic variety of particle escape environments. This might be the cause of the complicated Galactic cosmic-ray spectral shape measured on Earth.

We report on the results of our detailed analyses on the peculiar recombining plasma of the supernova remnant (SNR) G359.1$-$0.5, and the interacting CO clouds. Combining {\it Chandra} and {\it Suzaku} data, we estimated the ionization state of the plasma with a careful treatment of the background spectrum. The average spectrum showed a remarkably large deviation of the electron temperature ($\sim$0.17 keV) from the initial temperature ($>$ 16 keV), indicating that the plasma is in a highly recombination-dominant state. On the other hand, the recombination timescale $({\it n_{\rm e} t})$ is comparable to those of the other recombining SNRs ($\sim4.2 \times 10^{11}$ cm$^{-3}$ s). We also searched for spatial variation of the plasma parameters, but found no significant differences. Using $^{12}$CO($J$=2--1) data obtained with NANTEN2, we found a new, plausible candidate for the interacting CO cloud, which has a line-of-sight velocity of $\sim -$20 km s$^{-1}$. This indicates that the SNR is located at a distance of $\sim$4 kpc, which is the foreground of the Galactic center, as previously reported. The associated CO cloud does not show clear spatial coincidence with the nearby GeV/TeV emission, indicating that the origins of the GeV/TeV emission are likely unrelated to G359.1$-$0.5.

Neutron stars are formed in core-collapse supernova explosions, where a large number of neutrinos are emitted. In this paper, supernova neutrino light curves are computed for the cooling phase of protoneutron stars, which lasts a few minutes. In the numerical simulations, 90 models of the phenomenological equation of state with different incompressibilities, symmetry energies, and nucleon effective masses are employed for a comprehensive study. It is found that the cooling timescale is longer for a model with a larger neutron star mass and a smaller neutron star radius. Furthermore, a theoretical expression of the cooling timescale is presented as a function of the mass and radius and it is found to describe the numerical results faithfully. These findings suggest that diagnosing the mass and radius of a newly formed neutron star using its neutrino signal is possible.

We study the influence of density-dependent symmetry energy at high densities in simulations of core-collapse supernovae, black hole formation and proto-neutron star cooling by extending the relativistic mean field (RMF) theory used for the Shen EOS table. We adopt the extended RMF theory to examine the density dependence of the symmetry energy with a small value of the slope parameter $L$, while the original properties of the symmetric nuclear matter are unchanged. In order to assess matter effects at high densities, we perform numerical simulations of gravitational collapse of massive stars adopting the EOS table at high densities beyond $10^{14}$ g/cm$^3$ with the small $L$ value, which is in accord with the experimental and observational constraints, and compare them with the results obtained by using the Shen EOS. Numerical results for 11.2M$_{\odot}$ and 15M$_{\odot}$ stars exhibit minor effects around the core bounce and in the following evolution for 200 ms. Numerical results for 40M$_{\odot}$ and 50M$_{\odot}$ stars reveal a shorter duration toward the black hole formation with a smaller maximum mass for the small $L$ case. Numerical simulations of proto-neutron star cooling over 10 s through neutrino emissions demonstrate increasing effects of the symmetry energy at high densities. Neutrino cooling drastically proceeds in a relatively long timescale with high luminosities and average energies with the small symmetry energy. Evolution toward the cold neutron star is affected because of the different behavior of neutron-rich matter while supernova dynamics around core bounce remains similar in less neutron-rich environments.

We present a concept of the X-ray imaging system with high angular-resolution and moderate sensitivity. In this concept, a two-dimensional detector, i.e., imager, is put at a slightly out-of-focused position of the focusing mirror, rather than just at the mirror focus as in the standard optics, to capture the miniature image of objects. In addition, a set of multi-grid masks (or a modulation collimator) is installed in front of the telescope. We find that the masks work as a coded aperture camera and that they boost the angular resolution of the focusing optics. The major advantage of this concept is that a much better angular resolution an order of 2--3 or more than in the conventional optics is achievable, while a high throughput (large effective area) is maintained, which is crucial in photon-limited high-energy astronomy, because any type of mirrors, including light-weight reflective mirrors, can be employed in our concept. If the signal-to-noise ratio is sufficiently high, we estimate that angular resolutions at the diffraction limit of 4" and 0"4 at about 7 keV can be achieved with a pair of masks at separations of 1~m and 100~m, respectively, at the diffraction limit.

The cooling process of a protoneutron star is investigated with focus on its sensitivity to properties of hot and dense matter. An equation of state, which includes the nucleon effective mass and nuclear symmetry energy at twice the saturation density as control parameters, is constructed for systematic studies. The numerical code utilized in this study follows a quasi-static evolution of a protoneutron star solving the general-relativistic stellar structure with neutrino diffusion. The cooling timescale evaluated from the neutrino light curve is found to be longer for the models with larger effective masses and smaller symmetry energies at high densities. The present results are compared with those for other equations of state and it is found that they are consistent in terms of their dependences on the effective mass and neutron star radius.

We study a hierarchical triple system with the Kozai-Lidov mechanism, and analyse the cumulative shift of periastron time of a binary pulsar by the emission of gravitational waves. Time evolution of the osculating orbital elements of the triple system is calculated by directly integrating the first-order post-Newtonian equations of motion. The Kozai-Lidov mechanism will bend the evolution curve of the cumulative shift when the eccentricity becomes large. We also investigate the parameter range of mass and semi-major axis of the third companion with which the bending of the cumulative-shift curve could occur within 100 years.

162173 Ryugu, the target of Hayabusa2, has a round shape with an equatorial ridge, which is known as a spinning top-shape. A strong centrifugal force is a likely contributor to Ryugu's top-shaped features. Observations by Optical Navigation Camera onboard Hayabusa2 show a unique longitudinal variation in geomorphology; the western side of this asteroid, later called the western bulge, has a smooth surface and a sharp equatorial ridge, compared to the other side. Here, we propose a structural deformation process that generated the western bulge. Applying the mission-derived shape model, we employ a finite element model technique to analyze the locations that experience structural failure within the present shape. Assuming that materials are uniformly distributed, our model shows the longitudinal variation in structurally failed regions when the spin period is shorter than ~3.75 h. Ryugu is structurally intact in the subsurface region of the western bulge while other regions are sensitive to structural failure. We infer that this variation is indicative of the deformation process that occurred in the past, and the western bulge is more relaxed structurally than the other region. Our analysis also shows that this deformation process might occur at a spin period between ~3.5 h and ~3.0 h, providing the cohesive strength ranging between ~4 Pa and ~10 Pa.

The Optical Navigation Camera (ONC-T, ONC-W1, ONC-W2) onboard Hayabusa2 are also being used for scientific observations of the mission target, C-complex asteroid 162173 Ryugu. Science observations and analyses require rigorous instrument calibration. In order to meet this requirement, we have conducted extensive inflight observations during the 3.5 years of cruise after the launch of Hayabusa2 on 3 December 2014. In addition to the first inflight calibrations by Suzuki et al. (2018), we conducted an additional series of calibrations, including read-out smear, electronic-interference noise, bias, dark current, hot pixels, sensitivity, linearity, flat-field, and stray light measurements for the ONC. Moreover, the calibrations, especially flat-fields and sensitivities, of ONC-W1 and -W2 are updated for the analysis of the low-altitude (i.e., high-resolution) observations, such as the gravity measurement, touchdowns, and the descents for MASCOT and MINERVA-II payload releases. The radiometric calibration for ONC-T is also updated in this study based on star and Moon observations. Our updated inflight sensitivity measurements suggest the accuracy of the absolute radiometric calibration contains less than 1.8% error for the ul-, b-, v-, Na-, w-, and x-bands based on star calibration observations and ~5% for the p-band based on lunar calibration observations. The radiance spectra of the Moon, Jupiter, and Saturn from the ONC-T show good agreement with the spacecraft-based observations of the Moon from SP/SELENE and WAC/LROC and with ground-based telescopic observations for Jupiter and Saturn.

The future detection projects of gravitational waves are expected to have the sensitivity of detecting the primordial gravitational waves, which may be useful to get new insight into the very early universe. It is essential to analyze the evolution equation of the gravitational waves to estimate the present field strength of the primordial gravitational waves. In this paper, we obtain analytic solutions of the gravitational wave equation in the presence of non-relativistic matters as well as the cosmological constant. Although it is difficult to obtain the solution directly, we find that the equation for the square of the amplitudes has a simple polynomial solution. This quantity, which is directly related to the energy density of the gravitational waves, turns out to be useful to construct analytic solutions for the amplitudes by using Weierstrass's elliptic functions.

We present an X-ray study of the GeV gamma-ray supernova remnant (SNR) HB 21 with Suzaku. HB 21 is interacting with molecular clouds and the faintest in the GeV band among known GeV SNRs. We discovered strong radiative recombination continua of Si and S from the center of the remnant, which provide the direct evidence of a recombining plasma (RP). The total emission can be explained with the RP and ionizing plasma components. The electron temperature and recombination timescale of the RP component were estimated as 0.17 (0.15-0.18) keV and 3.2 (2.0-4.8) $\times$ 10$^{11}$ s cm$^{-3}$, respectively. The estimated age of the RP (RP age; $\sim$ 170 kyr) is the longest among known recombining GeV SNRs, because of very low density of electrons ($\sim$ 0.05 cm$^{-3}$). We have examined dependencies of GeV spectral indices on each of RP ages and SNR diameters for nine recombining GeV SNRs. Both showed possible positive correlations, indicating that both the parameters can be good indicators of properties of accelerated protons, for instance, degree of escape from the SNR shocks. A possible scenario for a process of proton escape is introduced; interaction with molecular clouds makes weaker magnetic turbulence and cosmic-ray protons escape, simultaneously cooling down the thermal electrons and generate an RP.

A proto-neutron star (PNS) is a newly formed compact object in a core collapse supernova. In this Letter, the neutrino emission from the cooling process of a PNS is investigated using two types of nuclear equation of state (EOS). It is found that the neutrino signal is mainly determined by the high-density EOS. The neutrino luminosity and mean energy are higher and the cooling time scale is longer for the softer EOS. Meanwhile, the neutrino mean energy and the cooling time scale are also affected by the low-density EOS because of the difference in the population of heavy nuclei. Heavy nuclei have a large scattering cross section with neutrinos owing to the coherent effects and act as thermal insulation near the surface of a PNS. The neutrino mean energy is higher and the cooling time scale is longer for an EOS with a large symmetry energy at low densities, namely a small density derivative coefficient of the symmetry energy, $L$.

We present a novel Lyman alpha (Ly$\alpha$) radiative transfer code, SEURAT, where line scatterings are solved adaptively with the resolution of the smoothed particle hydrodynamics (SPH). The radiative transfer method implemented in SEURAT is based on a Monte Carlo algorithm in which the scattering and absorption by dust are also incorporated. We perform standard test calculations to verify the validity of the code; (i) emergent spectra from a static uniform sphere, (ii) emergent spectra from an expanding uniform sphere, and (iii) escape fraction from a dusty slab. Thereby we demonstrate that our code solves the Ly$\alpha$ radiative transfer with sufficient accuracy. We emphasise that SEURAT can treat the transfer of Ly$\alpha$ photons even in highly complex systems that have significantly inhomogeneous density fields. The high adaptivity of SEURAT is desirable to solve the propagation of Ly$\alpha$ photons in the interstellar medium of young star-forming galaxies like Ly$\alpha$ emitters (LAEs). Thus, SEURAT provides a powerful tool to model the emergent spectra of Ly$\alpha$ emission, which can be compared to the observations of LAEs.

We present results obtained from a Suzaku observation of the accretion powered X-ray pulsar GX 1+4. Broad-band continuum spectrum of the pulsar was found to be better described by a simple model consisting of a blackbody component and an exponential cutoff power-law than the previously used compTT continuum model. Though the pulse profile had a sharp dip in soft X-rays ($<$10 keV), phase-resolved spectroscopy confirmed that the dimming was not due to increase in photoelectric absorption. Phase-sliced spectral analysis showed the presence of a significant spectral modulation beyond 10 keV except for the dip phase. A search for the presence of cyclotron resonance scattering feature in the Suzaku spectra yielded a negative result. Iron K-shell (K$_\alpha$ and K$_\beta$) emission lines from nearly neutral iron ions ($<$Fe III) were clearly detected in the source spectrum. A significant K$_\alpha$ emission line from almost neutral Ni atoms was detected for the first time in this source. We estimated the iron abundance of $\sim$80 % of the solar value and Ni/Fe abundance ratio of about two times of the solar value. We searched for a iron Ly$_\alpha$ emission line and found a significant improvement in the spectral fitting by inclusion of this line.

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.

A new table of the nuclear equation of state (EOS) based on realistic nuclear potentials is constructed for core-collapse supernova numerical simulations. Adopting the EOS of uniform nuclear matter constructed by two of the present authors with the cluster variational method starting from the Argonne v18 and Urbana IX nuclear potentials, the Thomas-Fermi calculation is performed to obtain the minimized free energy of a Wigner-Seitz cell in non-uniform nuclear matter. As a preparation for the Thomas-Fermi calculation, the EOS of uniform nuclear matter is modified so as to remove the effects of deuteron cluster formation in uniform matter at low densities. Mixing of alpha particles is also taken into account following the procedure used by Shen et al. (1998, 2011). The critical densities with respect to the phase transition from non-uniform to uniform phase with the present EOS are slightly higher than those with the Shen EOS at small proton fractions. The critical temperature with respect to the liquid-gas phase transition decreases with the proton fraction in a more gradual manner than in the Shen EOS. Furthermore, the mass and proton numbers of nuclides appearing in non-uniform nuclear matter with small proton fractions are larger than those of the Shen EOS. These results are consequences of the fact that the density derivative coefficient of the symmetry energy of our EOS is smaller than that of the Shen EOS.

The spectrum of the supernova relic neutrino (SRN) background from past stellar collapses including black hole formation (failed supernovae) is calculated. The redshift dependence of the black hole formation rate is considered on the basis of the metallicity evolution of galaxies. Assuming the mass and metallicity ranges of failed supernova progenitors, their contribution to SRNs is quantitatively estimated for the first time. Using this model, the dependences of SRNs on the cosmic star formation rate density, shock revival time and equation of state are investigated. The shock revival time is introduced as a parameter that should depend on the still unknown explosion mechanism of core collapse supernovae. The dependence on equation of state is considered for failed supernovae, whose collapse dynamics and neutrino emission are certainly affected. It is found that the low-energy spectrum of SRNs is mainly determined by the cosmic star formation rate density. These low-energy events will be observed in the Super-Kamiokande experiment with gadolinium-loaded water.

Nuclear equation of state (EoS) plays an important role in understanding the formation of compact objects such as neutron stars and black holes. The true nature of the EoS has been a matter of debate at any density range not only in the nuclear physics but also in the astronomy and astrophysics. We have constructed a database of EoSs by compiling data from the literature. Our database contains the basic properties of the nuclear EoS of symmetric nuclear matter and of pure neutron matter. It also includes the detailed information about the theoretical models, for example the adopted methods and assumptions in individual models. The novelty of the database is to consider new experimental probes such as the symmetry energy, its slope relative to the baryon density, and the incompressibility, which enables the users to check their model dependences. We demonstrate the performance of the EOSDB through the examinations of the model dependence among different nuclear EoSs. It is reveled that some theoretical EoSs, which is commonly used in astrophysics, do not satisfactorily agree with the experimental constraints.

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 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.

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 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 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.

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.

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 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 consider brane gas models based on type II string theories and analyze the mass, the Ramond-Ramond charge and the charge on moduli fluctuations of branes wrapping over cycles of a compactified space in the four-dimensional Einstein frame. A six-dimensional torus and Calabi-Yau threefolds are considered for the Kaluza-Klein reduction. A large volume of the compactified space and a weak string coupling gives rise to point particles of the wrapped branes which have a light mass and a small charge of the Ramond-Ramond flux and of the moduli fluctuations, while the particles become very heavy in the string frame. We find that the masses and the charges satisfy the sea-saw like dual relations which become time-independent in the four-dimensional Einstein frame.

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.

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.

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 consider a compactification with a six-dimensional torus in the type II brane gas models. We show that the dilaton and the scale of each cycle of the internal space are fixed in the presence of NS5-brane and Kaluza-Klein monopoles as well as D-branes with the gauge fields. We can construct various models that lead to fixed moduli by using T-duality transformations.

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 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 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 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 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 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.

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 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 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.

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.

We study the influence of the symmetry energy of the equation of state on the thermal evolution of protoneutron stars and the properties of supernova neutrinos by the numerical simulations after the protoneutron stars are formed. As for the equation of state (EOS) of nuclear matter, we take two EOS's with different symmetry energies obtained by the relativistic mean field theory. We find the symmetry energy plays the essential role on the evolution of lepton profiles and the neutrino fluxes.