In the crust of a neutron star, global torsional oscillations could occur in two elastic layers. The outer and inner layers are composed of spherical and cylindrical nuclei and of cylindrical holes (tubes) and spherical holes (bubbles), respectively, while between these two layers, a phase of slab-like (lasagna) nuclei with vanishingly small elasticity is sandwiched. In this work, we update systematic calculations of the eigenfrequencies of the fundamental oscillations in the inner layer by newly allowing for the presence of tubes. We find that the frequencies still depend strongly on the slope parameter of the nuclear symmetry energy, $L$, while being almost independent of the incompressibility of symmetric nuclear matter. We also find that the fundamental frequencies in the inner layer can become smaller than those in the outer layer because the tube phase has a relatively small shear modulus and at the same time dominates the inner layer in thickness. As a result, we can successfully explain not only the quasi-periodic oscillations originally discovered in the observed X-ray afterglow of the giant flare of SGR 1806--20 but also many others recently found by a Bayesian procedure.

We systematically examine torsional shear oscillations of neutron star crusts by newly taking into account the possible presence of the phase of cylindrical nuclei. In this study, we neglect an effect of magnetic fields, under which the shear oscillations can be damped by the magnetic interaction. First, by identifying the low frequency quasi-periodic oscillations (QPOs) observed in the soft-gamma repeaters (SGRs) as the fundamental torsional oscillations, we constrain the slope parameter of the nuclear symmetry energy, $L$, for reasonable values of the star's mass $M$ and radius $R$. Meanwhile, we find that the 1st overtone of torsional oscillations obtained for given $M$ and $R$ can be expressed well as a function of a new parameter $\varsigma\equiv (K_0^4 L^5)^{1/9}$, where $K_0$ is the incompressibility of symmetric nuclear matter. Assuming that the lowest of the QPO frequencies above 500 Hz observed in SGR 1806-20 comes from the 1st overtone, we can constrain the value of $\varsigma$. Then, for each neutron star model, such a value of $L$ as can be obtained from the observed low frequency QPOs translates to the optimal value of $K_0$ via the above constraint on $\varsigma$. Finally, its consistency with allowed values of $K_0$ from empirical giant monopole resonances leads to neutron star models with relatively low mass and large radius, which are qualitatively similar to the prediction in earlier investigations. This result suggests that $L\simeq 58-73$ MeV, even when uncertainties in the neutron superfluid density inside the phase of cylindrical nuclei are allowed for.

With various sets of the parameters that characterize the equation of state (EOS) of nuclear matter, we systematically examine the thickness of a neutron star crust and of the pasta phases contained therein. Then, with respect to the thickness of the phase of spherical nuclei, the thickness of the cylindrical phase, and the crust thickness, we successfully derive fitting formulas that express the ratio of each thickness to the star's radius as a function of the star's compactness, the incompressibility of symmetric nuclear matter, and the density dependence of the symmetry energy. In particular, we find that the thickness of the phase of spherical nuclei has such a strong dependence on the stellar compactness as the crust thickness, but both of them show a much weaker dependence on the EOS parameters. Thus, via determination of the compactness, the thickness of the phase of spherical nuclei as well as the crust thickness can be constrained reasonably, even if the EOS parameters remain to be well-determined.

We consider torsional oscillations that are trapped in a layer of spherical-hole (bubble) nuclear structure, which is expected to occur in the deepest region of the inner crust of a neutron star. Because this layer intervenes between the phase of slab nuclei and the outer core of uniform nuclear matter, torsional oscillations in the bubble phase can be excited separately from usual crustal torsional oscillations. We find from eigenmode analyses for various models of the equation of state of uniform nuclear matter that the fundamental frequencies of such oscillations are almost independent of the incompressibility of symmetric nuclear matter, but strongly depend on the slope parameter of the nuclear symmetry energy $L$. Although the frequencies are also sensitive to the entrainment effect, i.e., what portion of nucleons outside bubbles contribute to the oscillations, by having such a portion fixed, we can successfully fit the calculated fundamental frequencies of torsional oscillations in the bubble phase inside a star of specific mass and radius as a function of $L$. By comparing the resultant fitting formula to the frequencies of quasi-periodic oscillations (QPOs) observed from the soft-gamma repeaters, we find that each of the observed low-frequency QPOs can be identified either as a torsional oscillation in the bubble phase or as a usual crustal oscillation, given generally accepted values of $L$ for all the stellar models considered here.

Quasi-periodic oscillations (QPOs) discovered in soft-gamma repeaters (SGRs) are expected to help us to study the properties of matter in neutron stars. In earlier investigations, we identified the QPOs of frequencies below $\sim100$ Hz observed in giant flares of SGR 1806$-$20 and SGR 1900+14 as the crustal torsional oscillations. For this purpose, we calculated the frequencies of the fundamental torsional oscillations with various angular indices $\ell$, by changing the stellar mass and radius. In this work, we try to explain the additional QPO frequencies recently reported by Huppenkothen et al. (2014a,b) within the same framework as before except that we newly take into account the effect of electron screening, which acts to decrease the frequencies by a small amount. Those QPOs were discovered in two different SGRs, i.e., SGR 1806$-$20 and SGR J1550$-$5418. Then, we find that the newly observed QPO frequency in SGR 1806$-$20 can be still identified as one of the frequencies of the fundamental torsional oscillations, while those in SGR J1550$-$5418 can also be explained in terms of the torsional oscillations although the relevant angular indices are difficult to identify.

Neutrons stars lighter than the Sun are basically composed of nuclear matter of density up to around twice normal nuclear density. In our recent analyses, we showed that possible simultaneous observations of masses and radii of such neutron stars could constrain $\eta\equiv(K_0L^2)^{1/3}$, a combination of the incompressibility of symmetric nuclear matter $K_0$ and the density derivative of the nuclear symmetry energy $L$ that characterizes the theoretical mass-radius relation. In this paper, we focus on the mass-radius constraint of the X-ray burster 4U 1724-307 given by Suleimanov et al. (2011). We therefrom obtain the constraint that $\eta$ should be larger than around 130 MeV, which in turn leads to $L$ larger than around 110, 98, 89, and 78 MeV for $K_0=180$, 230, 280, and 360 MeV. Such a constraint on $L$ is more or less consistent with that obtained from the frequencies of quasi-periodic oscillations in giant flares observed in soft-gamma repeaters.

Neutron stars, produced at the death of massive stars, are often regarded as giant neutron-rich nuclei. This picture is especially relevant for low-mass (below about solar mass) neutron stars, where non-nucleonic components are not expected to occur. Due to the saturation property of nucleonic matter, leading to the celebrated liquid-drop picture of atomic nuclei, empirical nuclear masses and radii can be approximately expressed as function of atomic mass number. It is, however, not straightforward to express masses and radii of neutron stars even in the low-mass range where the structure is determined by a balance between the pressure of neutron-rich nucleonic matter and the gravity. Such expressions would be of great use given possible simultaneous mass and radius measurements. Here we successfully construct theoretical formulas for the masses and radii of low-mass neutron stars from various models that are consistent with empirical masses and radii of stable nuclei. In this process, we discover a new equation-of-state parameter that characterizes the structure of low-mass neutron stars. This parameter, which plays a key role in connecting the mass-radius relation of the laboratory nuclei to that of the celestial objects, could be constrained from future observations of low-mass neutron stars.

Phenomenological approach to inhomogeneous nuclear matter is useful to describe fundamental properties of atomic nuclei and neutron star crusts in terms of the equation of state of uniform nuclear matter. We review a series of researches that we have developed by following this approach. We start with more than 200 equations of state that are consistent with empirical masses and charge radii of stable nuclei and then apply them to describe matter radii and masses of unstable nuclei, proton elastic scattering and total reaction cross sections off unstable nuclei, and nuclei in neutron star crusts including nuclear pasta. We finally discuss the possibility of constraining the density dependence of the symmetry energy from experiments on unstable nuclei and even observations of quasi-periodic oscillations in giant flares of soft gamma-ray 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.

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.

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.

We construct the equation of state (EOS) of dense matter covering a wide range of temperature, proton fraction, and density for the use of core-collapse supernova simulations. The study is based on the relativistic mean-field (RMF) theory, which can provide an excellent description of nuclear matter and finite nuclei. The Thomas--Fermi approximation in combination with assumed nucleon distribution functions and a free energy minimization is adopted to describe the non-uniform matter, which is composed of a lattice of heavy nuclei. We treat the uniform matter and non-uniform matter consistently using the same RMF theory. We present two sets of EOS tables, namely EOS2 and EOS3. EOS2 is an update of our earlier work published in 1998 (EOS1), where only the nucleon degree of freedom is taken into account. EOS3 includes additional contributions from $\Lambda$ hyperons. The effect of $\Lambda$ hyperons on the EOS is negligible in the low-temperature and low-density region, whereas it tends to soften the EOS at high density. In comparison with EOS1, EOS2 and EOS3 have an improved design of ranges and grids, which covers the temperature range $T=0.1$--$10^{2.6}$ MeV with the logarithmic grid spacing $\Delta \log_{10}(T/\rm{[MeV]})=0.04$ (92 points including T=0), the proton fraction range $Y_p=0$--0.65 with the linear grid spacing $\Delta Y_p = 0.01$ (66 points), and the density range $\rho_B=10^{5.1}$--$10^{16}\,\rm{g\,cm^{-3}}$ with the logarithmic grid spacing $\Delta \log_{10}(\rho_B/\rm{[g\,cm^{-3}]}) = 0.1$ (110 points).

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.

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.

An equation of state (EOS) for uniform nuclear matter is constructed at zero and finite temperatures with the variational method starting from the realistic nuclear Hamiltonian composed of the Argonne V18 and UIX potentials. The energy is evaluated in the two-body cluster approximation with the three-body-force contribution treated phenomenologically so as to reproduce the empirical saturation conditions. The obtained energies for symmetric nuclear matter and neutron matter at zero temperature are in fair agreement with those by Akmal, Pandharipande and Ravenhall, and the maximum mass of the neutron star is 2.2 Msolar. At finite temperatures, a variational method by Schmidt and Pandharipande is employed to evaluate the free energy, which is used to derive various thermodynamic quantities of nuclear matter necessary for supernova simulations. The result of this variational method at finite temperatures is found to be self-consistent.

We examine how the properties of inhomogeneous nuclear matter at subnuclear densities depend on the density dependence of the symmetry energy. Using a macroscopic nuclear model we calculate the size and shape of nuclei in neutron star matter at zero temperature in a way dependent on the density dependence of the symmetry energy. We find that for smaller symmetry energy at subnuclear densities, corresponding to larger density symmetry coefficient L, the charge number of nuclei is smaller, and the critical density at which matter with nuclei or bubbles becomes uniform is lower. The decrease in the charge number is associated with the dependence of the surface tension on the nuclear density and the density of a sea of neutrons, while the decrease in the critical density can be generally understood in terms of proton clustering instability in uniform matter.

We propose that spontaneous particle--anti-particle pair creations from the discharged vacuum caused by the strong interactions in dense matter are major sources of $\gamma$-ray bursts. Two neutron star collisions or black hole-neutron star mergers at cosmological distance could produce a compact object with its density exceeding the critical density for pair creations. The emitted anti-particles annihilate with corresponding particles at the ambient medium. This releases a large amount of energy. We discuss the spontaneous $p\bar{p}$ pair creations within two neutron star collision and estimate the exploded energy from $p\bar{p}$ annihilation processes. The total energy could be around $10^{51} - 10^{53}$ erg depending on the impact parameter of colliding neutron stars. This value fits well into the range of the initial energy of the most energetic $\gamma$-ray bursts.

We construct the equation of state (EOS) of nuclear matter at finite temperature and density with various proton fractions within the relativistic mean field (RMF) theory for the use in the supernova simulations. The Thomas-Fermi approximation is adopted to describe the non-uniform matter where we consider nucleus, alpha-particle, proton and neutron in equilibrium. We treat the uniform matter and non-uniform matter consistently using the RMF theory. We tabulate the outcome as the pressure, free energy, entropy etc, with enough mesh points in wide ranges of the temperature, proton fraction, and baryon mass density.

We construct the equation of state (EOS) of nuclear matter using the relativistic mean field (RMF) theory in the wide density, temperature range with various proton fractions for the use of supernova simulation and the neutron star calculations. We first construct the EOS of homogeneous nuclear matter. We use then the Thomas-Fermi approximation to describe inhomogeneous matter, where heavy nuclei are formed together with free nucleon gas. We discuss the results on free energy, pressure and entropy in the wide range of astrophysical interest. As an example, we apply the resulting EOS on the neutron star properties by using the Oppenheimer-Volkoff equation.