Source author record

Ken'ichiro Nakazato

Ken'ichiro Nakazato appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

26works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

26 published item(s)

preprint2022arXiv

Effects of nuclear matter and composition in core-collapse supernovae and long-term proto-neutron star cooling

We study the influence of hot and dense matter in core-collapse supernovae by adopting up-to-date nuclear equation of state (EOS) based on the microscopic nuclear many-body frameworks. We explore effects of EOS based on the Dirac Brueckner Hartree-Fock theory through comparisons with those based on the variational method. We also examine effects of the differences in the composition of nuclei and nucleons by using the same EOS by the variational method but employing two different treatments in computations of nuclear abundances. We perform numerical simulations of core-collapse supernovae adopting the three EOSs. We also perform numerical simulations of the long-term evolution over 70 s of the proto-neutron star cooling. We show that impacts by different modeling of composition are remarkable as in those by different treatments of uniform matter in the gravitational collapse, bounce, and shock propagation. The cooling of proto-neutron star and the resulting neutrino emission are also affected by the compositional difference even if the same treatment in computing uniform matter of EOS.

preprint2022arXiv

Observing Supernova Neutrino Light Curves with Super-Kamiokande. II. Impact of the Nuclear Equation of State

The late-time evolution of the neutrino event rate from supernovae is evaluated for Super-Kamiokande using simulated results of proto-neutron star (PNS) cooling. In the present work we extend the result of Suwa et al. (2019) [arXiv:1904.09996], which studied the dependence on the PNS mass, but focus on the impact of the nuclear equation of state (EOS). We find that the neutrino event rate depends on both the high-density and low-density EOS, where the former determines the radius of the PNS and the latter affects its surface temperature. Based on the present evaluation of the neutrino event rate we propose a new analysis method to extract the time variability of the neutrino average energy taking into account the statistical error in the observation.

preprint2022arXiv

Observing Supernova Neutrino Light Curves with Super-Kamiokande. III. Extraction of Mass and Radius of Neutron Stars from Synthetic Data

Neutrinos are guaranteed to be observable from the next galactic supernova (SN). Optical light and gravitational waves are also observable, but may be difficult to observe if the location of the SN in the galaxy or the details of the explosion are unsuitable. The key to observing the next supernova is to first use neutrinos to understand various physical quantities and then link them to other signals. In this paper, we present Monte Carlo sampling calculations of neutrino events from galactic supernova explosions observed with Super-Kamiokande. The analytical solution of neutrino emission, which represents the long-term evolution of neutrino-light curve from supernovae, is used as a theoretical template. It gives the event rate and event spectrum through inverse beta decay interactions with explicit model parameter dependence. Parameter estimation is performed on these simulated sample data by fitting least squares using the analytical solution. The results show that the mass, radius and total energy of a remnant neutron star produced by a SN can be determined with an accuracy of $\sim 0.1M_\odot$, $\sim 1$ km, and $\sim 10^{51}$ erg, respectively, for a galactic SN at 8 kpc.

preprint2021arXiv

Analytic solutions for neutrino-light curves of core-collapse supernovae

Neutrinos are a guaranteed signal from supernova explosions in the Milky Way, and a most valuable messenger that can provide us with information about the deepest parts of supernovae. In particular, neutrinos will provide us with physical quantities, such as the radius and mass of protoneutron stars (PNS), which are the central engine of supernovae. This requires a theoretical model that connects observables such as neutrino luminosity and average energy with physical quantities. Here, we show analytic solutions for the neutrino-light curve derived from the neutrino radiation transport equation by employing the diffusion approximation and the analytic density solution of the hydrostatic equation for a PNS. The neutrino luminosity and the average energy as functions of time are explicitly presented, with dependence on PNS mass, radius, the total energy of neutrinos, surface density, and opacity. The analytic solutions provide good representations of the numerical models from a few seconds after the explosion and allow a rough estimate of these physical quantities to be made from observational data.

preprint2020arXiv

A New Approach to Mass and Radius of Neutron Stars with Supernova Neutrinos

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.

preprint2020arXiv

Developing an end-to-end simulation framework of supernova neutrino detection

Massive stars can explode as supernovae at the end of their life cycle, releasing neutrinos whose total energy reaches $10^{53}$ erg. Moreover, neutrinos play key roles in supernovae, heating and reviving the shock wave as well as cooling the resulting protoneutron star. Therefore, neutrino detectors are waiting to observe the next galactic supernova and several theoretical simulations of supernova neutrinos are underway. While these simulation concentrate mainly on only the first one second after the supernova bounce, the only observation of a supernova with neutrinos, SN 1987A, revealed that neutrino emission lasts for more than 10 seconds. For this reason, long-time simulation and analysis tools are needed to compare theories with the next observation. Our study is to develop an integrated supernova analysis framework to prepare an analysis pipeline for treating galactic supernovae observations in the near future. This framework deals with the core-collapse, bounce and proto-neutron star cooling processes, as well as with neutrino detection on earth in a consistent manner. We have developed a new long-time supernova simulation in one dimension that explodes successfully and computes the neutrino emission for up to 20 seconds. Using this model we estimate the resulting neutrino signal in the Super-Kamiokande detector to be about 1,800 events for an explosion at 10 kpc and discuss its implications in this paper. We compare this result with the SN 1987A observation to test its reliability.

preprint2020arXiv

Effects of symmetry energy on the radius and tidal deformability of neutron stars in relativistic mean-field model

The radii and tidal deformabilities of neutron stars are investigated in the framework of relativistic mean-field (RMF) model with different density-dependent behaviors of symmetry energy. To study the effects of symmetry energy on the properties of neutron stars, an $ω$ meson and $ρ$ meson coupling term is included in a popular RMF Lagrangian, i.e. the TM1 parameter set, which is used for the widely used supernova equation of state (EoS) table. The coupling constants relevant to the vector-isovector meson, $ρ$, are refitted by a fixed symmetry energy at subsaturation density and its slope at saturation density, while other coupling constants remain the same as the original ones in TM1 so as to update the supernova EoS table. The radius and mass of maximum neutron stars are not so sensitive to the symmetry energy in these family TM1 parameterizations. However, the radii at intermediate mass region are strongly correlated with the slope of symmetry energy. Furthermore, the dimensionless tidal deformabilities of neutron stars are also calculated within the associated Love number. We find that its value at $1.4 M_\odot$ has a linear correlation to the slope of symmetry energy being different from the previous studied. With the latest constraints of tidal deformabilities from GW170817 event, the slope of symmetry energy at nuclear saturation density should be smaller than $60$ MeV in the family TM1 parameterizations. This fact supports the usage of lower symmetry energy slope for the update supernova EoS, which is applicable to simulations of neutron star merger. Furthermore, the analogous analysis are also done within the family IUFSU parameter sets. It is found that the correlations between the symmetry energy slope with the radius and tidal deformability at $1.4 M_\odot$ have very similar linear relations in these RMF models.

preprint2020arXiv

Possibility of rapid neutron star cooling with the realistic equation of state

Whether fast cooling processes occur or not is crucial for the thermal evolution of neutron stars. In particular, the threshold of the direct Urca process, which is one of the fast cooling processes, is determined by the interior proton fraction $Y_p$, or the nuclear symmetry energy. Since recent observations indicate the small radius of neutron stars, a low value is preferred for the symmetry energy. In this study, simulations of neutron star cooling are performed adopting three models for equation of state (EoS): Togashi, Shen, and LS220 EoSs. The Togashi EoS has been recently constructed with realistic nuclear potentials under finite temperature, and found to account for the small radius of neutron stars. As a result, we find that, since the direct Urca process is forbidden, the neutron star cooling is slow with use of the Togashi EoS. This is because symmetry energy of Togashi EoS is lower than those of other EoSs. Hence, in order to account for observed age and surface temperature of isolated neutron stars (INS) with use of the Togashi EoS, other fast cooling processes are needed regardless of the surface composition.

preprint2019arXiv

Influence of density dependence of symmetry energy in hot and dense matter for supernova simulations

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.

preprint2016arXiv

Gravitational-Wave Background from Binary Mergers and Metallicity Evolution of Galaxies

The cosmological evolution of the binary black hole (BH) merger rate and the energy density of the gravitational-wave (GW) background are investigated. To evaluate the redshift dependence of the BH formation rate, BHs are assumed to originate from low-metallicity stars, and the relations between the star formation rate, metallicity and stellar mass of galaxies are combined with the stellar mass function at each redshift. As a result, it is found that when the energy density of the GW background is scaled with the merger rate at the local Universe, the scaling factor does not depend on the critical metallicity for the formation of BHs. Also taking into account the merger of binary neutron stars, a simple formula to express the energy spectrum of the GW background is constructed for the inspiral phase. The relation between the local merger rate and the energy density of the GW background will be examined by future GW observations.

preprint2015arXiv

Spectrum of the Supernova Relic Neutrino Background and Metallicity Evolution of Galaxies

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.

preprint2014arXiv

Self-organized criticality in a spherically closed cellular automaton: Modeling soft gamma repeater bursts driven by magnetic reconnection

A new cellular automaton (CA) model is presented for the self-organized criticality (SOC) in recurrent bursts of soft gamma repeaters (SGRs), which are interpreted as avalanches of reconnection in the magnetosphere of neutron stars. The nodes of a regular dodecahedron and a truncated icosahedron are adopted as spherically closed grids enclosing a neutron star. It is found that the system enters the SOC state if there are sites where the expectation value of the added perturbation is nonzero. The energy distributions of SOC avalanches in CA simulations are described by a power law with a cutoff, which is consistent with the observations of SGR 1806-20 and SGR 1900+14. The power-law index is not universal and depends on the amplitude of the perturbation. This result shows that the SOC of SGRs can be illustrated not only by the crust quake model but also by the magnetic reconnection model.

preprint2013arXiv

A new equation of state for neutron star matter with nuclei in the crust and hyperons in the core

The equation of state for neutron stars in a wide-density range at zero temperature is constructed. The chiral quark-meson coupling model within relativistic Hartree-Fock approximation is adopted for uniform nuclear matter. The coupling constants are determined so as to reproduce the experimental data of atomic nuclei and hypernuclei. In the crust region, nuclei are taken into account within the Thomas-Fermi calculation. All octet baryons are considered in the core region, while only $Ξ^{-}$ appears in neutron stars. The resultant maximum mass of neutron stars is $1.95M_\odot$, which is consistent with the constraint from the recently observed massive pulsar, PSR J1614-2230.

preprint2013arXiv

Imprint of Explosion Mechanism on Supernova Relic Neutrinos

The spectrum and event rate of supernova relic neutrinos are calculated taking into account the dependence on the time it takes for the shock wave in supernova cores to revive. The shock revival time should depend on the still unknown explosion mechanism of collapse-driven supernovae. The contribution of black-hole-forming failed supernovae is also considered. The total event rate is higher for models with a longer shock revival time and/or a failed-supernova contribution. The hardness of the spectrum does not strongly depend on the shock revival time, but the spectrum becomes hard owing to the failed supernovae. Therefore, the shock-revival-time dependence of supernova relic neutrinos has different systematics from the fractions of failed supernovae.

preprint2013arXiv

Possible Constraints on the Density Dependence of the Nuclear Symmetry Energy from Quasiperiodic Oscillations in Soft Gamma Repeaters

We systematically examine the fundamental frequencies of shear torsional oscillations in neutron star crusts in a manner that is dependent on the parameter $L$ characterizing the poorly known density dependence of the symmetry energy. The identification of the lowest quasiperiodic oscillation (QPO) among the observed QPOs from giant flares in soft-gamma repeaters as the $/ell=2$ fundamental torsional oscillations enables us to constrain the parameter $L$ as $L/ge 47.4$ MeV, which is the most conservative restriction on $L$ obtained in the present work that assumes that the mass and radius of the flaring neutron stars range 1.4-1.8 $M_{/odot}$ and 10-14 km. Next, we identify one by one a set of the low-lying frequencies observed in giant flares as the fundamental torsional oscillations. The values of $L$ that can reproduce all the observed frequencies in terms of the torsional oscillations coupled with a part of dripped neutrons via entrainment effects are then constrained as 101.1 MeV $/le L /le$ 131.0 MeV. Alternatively, if only the second lowest frequency observed in SGR 1806$-$20 has a different origin, one obtains relatively low $L$ values ranging 58.0 MeV $/le L /le$ 85.3 MeV, which seem more consistent with other empirical constraints despite large uncertainties.

preprint2013arXiv

Stellar Core Collapse with Hadron-Quark Phase Transition

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.

preprint2013arXiv

Supernova Neutrino Light Curves and Spectra for Various Progenitor Stars: From Core Collapse to Proto-neutron Star Cooling

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.

preprint2012arXiv

Effect of superfluidity on neutron star crustal oscillations

We consider how superfluidity of dripped neutrons in the crust of a neutron star affects the frequencies of the crust's fundamental torsional oscillations. A nonnegligible superfluid part of dripped neutrons, which do not comove with nuclei, act to reduce the enthalpy density and thus enhance the oscillation frequencies. By assuming that the quasi-periodic oscillations observed in giant flares of soft gamma repeaters arise from the fundamental torsional oscillations and that the mass and radius of the neutron star is in the range of 1.4 < M/M_SUN < 1.8 and 10 km < R < 14 km, we constrain the density derivative of the symmetry energy as 100 MeV < L < 130 MeV, which is far severer than the previous one, L > 50 MeV, derived by ignoring the superfluidity.

preprint2012arXiv

Probing the Equation of State of Nuclear Matter via Neutron Star Asteroseismology

We general relativistically calculate the frequency of fundamental torsional oscillations of neutron star crusts, where we focus on the crystalline properties obtained from macroscopic nuclear models in a way depending on the equation of state of nuclear matter. We find that the calculated frequency is sensitive to the density dependence of the symmetry energy, but almost independent of the incompressibility of symmetric nuclear matter. By identifying the lowest-frequency quasi-periodic oscillation in giant flares observed from soft gamma-ray repeaters as the fundamental torsional mode and allowing for the dependence of the calculated frequency on stellar models, we provide a lower limit of the density derivative of the symmetry energy as $L\simeq 50$ MeV.

preprint2011arXiv

Curvature effect on nuclear pasta: Is it helpful for gyroid appearance?

In supernova cores and neutron star crusts, nuclei are thought to deform to rodlike and slablike shapes, which are often called nuclear pasta. We study the equilibrium properties of the nuclear pasta by using a liquid drop model with curvature corrections. It is confirmed that the curvature effect acts to lower the transition densities between different shapes. We also examine the gyroid structure, which was recently suggested as a different type of nuclear pasta by analogy with the polymer systems. The gyroid structure investigated in this paper is approximately formulated as an extension of the periodic minimal surface whose mean curvature vanishes. In contrast to our expectations, we find from the present approximate formulation that the curvature corrections act to slightly disfavor the appearance of the gyroid structure. By comparing the energy corrections in the gyroid phase and the hypothetical phases composed of d-dimensional spheres, where d is a general dimensionality, we show that the gyroid is unlikely to belong to a family of the generalized dimensional spheres.

preprint2011arXiv

Hyperon Matter and Black Hole Formation in Failed Supernovae

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

preprint2010arXiv

Exploring Hadron Physics in Black Hole Formations: a New Promising Target of Neutrino Astronomy

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.

preprint2010arXiv

Impact of Quarks and Pions on Dynamics and Neutrino Signal of Black Hole Formation in Non-rotating Stellar Core Collapse

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.

preprint2008arXiv

Oscillation and Future Detection of Failed Supernova Neutrinos from 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.

preprint2007arXiv

Numerical Study on Stellar Core Collapse and Neutrino Emission: Probe into the Spherically Symmetric Black Hole Progenitors with 3 - 30Msun Iron Cores

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.

preprint2005arXiv

Gravitational Collapse and Neutrino Emission of Population III Massive Stars

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