We report multi-wavelength monitoring observations of an M-dwarf flare star AD Leonis with Seimei Telescope (6150--7930 {\AA}), SCAT (Spectroscopic Chuo-university Astronomical Telescope; 3700--7500 {\AA}), NICER (Neutron Star Interior Composition Explorer; 0.2--12.0 keV), and collaborations of OISTER (Optical and Infrared Synergetic Telescopes for Education and Research) program. Twelve flares are detected in total which include ten H$\alpha$, four X-ray, and four optical-continuum flares; one of them is a superflare with the total energy of $\sim$ 2.0$\times$10$^{33}$ erg. We found that (1) during the superflare, the H$\alpha$ emission line full width at 1/8 maximum dramatically increases to 14 {\AA} from 8 {\AA} in the low-resolution spectra (R$\sim$ 2000) accompanied with the large white-light flares, (2) some weak H$\alpha$/X-ray flares are not accompanied with white-light emissions, and (3) the non-flaring emissions show clear rotational modulations in X-ray and H$\alpha$ intensity in the same phase. To understand these observational features, one-dimensional hydrodynamic flare simulations are performed by using the RADYN code. As a result of simulations, we found the simulated H$\alpha$ line profiles with hard and high-energy non-thermal electron beams are consistent with that of the initial phase line profiles of the superflares, while those with more soft- and/or weak-energy beam are consistent with those in decay phases, indicating the changes in the energy fluxes injected to the lower atmosphere. Also, we found that the relation between optical continuum and H$\alpha$ intensity is nonlinear, which can be one cause of the non-white-light flares. The flare energy budget exhibits diversity in the observations and models, and more observations of stellar flares are necessary for constraining the occurrence of various emission line phenomena in stellar flares.

We present optical observations of type Ia supernova (SN) 2019ein, starting at 2 days after the estimated explosion date. The spectra and the light curves show that SN 2019ein belongs to the High-Velocity (HV) and Bload Line groups with relatively rapid decline in the light curves (Delta m15(B) = 1.36 +- 0.02 mag) and the short rise time (15.37 +- 0.55 days). The Si II 6355 velocity, associated with a photospheric component but not with a detached high-velocity feature, reached ~ 20,000 km s-1 at 12 days before the B-band maximum. The line velocity however decreased very rapidly and smoothly toward the maximum light, where it was ~ 13,000 km s-1 as relatively low among HV SNe. This indicates that the speed of the spectral evolution of HV SNe Ia is correlated not only to the velocity at the maximum light, but also to the light curve decline rate like the case for Normal-Velocity (NV) SNe Ia. Spectral synthesis modeling shows that the outermost layer at > 17,000 km s-1 is well described by the O-Ne-C burning layer extending to at least 25,000 km s-1, and there is no unburnt carbon below 30,000 km s-1; these properties are largely consistent with the delayed detonation scenario, and are shared with the prototypical HV SN 2002bo despite the large difference in Delta m15(B). This structure is strikingly different from that derived for the well-studied NV SN 2011fe. We suggest that the relation between the mass of 56Ni (or Delta m15) and the extent of the O-Ne-C burning layer provides an important constraint on the explosion mechanism(s) of HV and NV SNe.

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.

Nuclear matter under the conditions of a supernova explosion unfolds into a rich variety of spatially structured phases, called nuclear pasta. We investigate the role of periodic network-like structures with negatively curved interfaces in nuclear pasta structures, by static and dynamic Hartree-Fock simulations in periodic lattices. As the most prominent result, we identify for the first time the {\it single gyroid} network structure of cubic chiral $I4_123$ symmetry, a well known configuration in nanostructured soft-matter systems, both as a dynamical state and as a cooled static solution. Single gyroid structures form spontaneously in the course of the dynamical simulations. Most of them are isomeric states. The very small energy differences to the ground state indicate its relevance for structures in nuclear pasta.

Nuclear matter under astrophysical conditions is explored with time-dependent and static Hartree-Fock calculations. The focus is in a regime of densities where matter segregates into liquid and gaseous phases unfolding a rich scenario of geometries, often called nuclear pasta shapes (e.g. spaghetti, lasagna). Particularly the appearance of the different phases depending on the proton fraction and the transition to uniform matter are investigated. In this context the neutron background density is of special interest, because it plays a crucial role for the type of pasta shape which is built. The study is performed in two dynamical ranges, one for hot matter and one at temperature zero to investigate the effect of cooling.

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 present simulations of neutron-rich matter at subnuclear densities, like supernova matter, with the time-dependent Hartree-Fock approximation at temperatures of several MeV. The initial state consists of $\alpha$ particles randomly distributed in space that have a Maxwell-Boltzmann distribution in momentum space. Adding a neutron background initialized with Fermi distributed plane waves the calculations reflect a reasonable approximation of astrophysical matter. This matter evolves into spherical, rod-like, and slab-like shapes and mixtures thereof. The simulations employ a full Skyrme interaction in a periodic three-dimensional grid. By an improved morphological analysis based on Minkowski functionals, all eight pasta shapes can be uniquely identified by the sign of only two valuations, namely the Euler characteristic and the integral mean curvature. In addition, we propose the variance in the cell density distribution as a measure to distinguish pasta matter from uniform matter.

We illustrate the formation of a thin-long structure of heavy nuclei by three-nucleus simultaneous collisions within time-dependent density functional theory. The impact parameter dependence for such formation is systematically demonstrated through clarifications of the difference between binary and ternary collision events. A new method for producing thin-long heavy nuclei in the laboratory is suggested, as well as the possible formation of the thin-long structure in hot dense matter such as that encountered in core collapse supernovae.

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.

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.

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.

Screening effects of electrons on inhomogeneous nuclear matter, which includes spherical, slablike, and rodlike nuclei as well as spherical and rodlike nuclear bubbles, are investigated in view of possible application to cold neutron star matter and supernova matter at subnuclear densities. Using a compressible liquid-drop model incorporating uncertainties in the surface tension, we find that the energy change due to the screening effects broadens the density region in which bubbles and nonspherical nuclei appear in the phase diagram delineating the energetically favorable shape of inhomogeneous nuclear matter. This conclusion is considered to be general since it stems from a model-independent feature that the electron screening acts to decrease the density at which spherical nuclei become unstable against fission and to increase the density at which uniform matter becomes unstable against proton clustering.

Recent deep Chandra LETG+HRC-S observations suggest the possibility that RX J1856.5--3754 is a compact star whose radiation radius is 3.8--8.2 km. In this paper, we systematically calculate the mass-radius relations of quark stars within the bag model. Assuming that RX J1856.5--3754 is a pure quark star, we derive an upper limit on its mass for various sets of the bag-model parameters.

Formation and decay of rod-like and slab-like nuclei, as may be encountered in the outer parts of neutron stars, are examined within a zero-temperature stability analysis with respect to perturbations inducing fission and proton clustering. The nonspherical nuclei are found to be stable against these perturbations, implying the possibility that the size of the stellar region containing such nuclei exceeds the equilibrium prediction.

Geometrical structure of matter at subnuclear densities is investigated in the presence of a degenerate gas of neutrinos as encountered in stellar collapse. The crystalline phases with spherical, cylindrical and planar nuclei as well as with spherical and cylindrical nuclear bubbles are considered by using a compressible liquid-drop model. This model allows for uncertainties in the lepton fraction $Y_{\rm L}$ in addition to those in the nuclear surface tension $E_{\rm surf}$ and in the proton chemical potential in bulk neutron matter $\mu_{\rm p}^{(0)}$. The phase diagrams obtained at zero temperature show that only the phases with rod-like and slab-like nuclei occur at typical values of $Y_{\rm L}$, $E_{\rm surf}$ and $\mu_{\rm p}^{(0)}$, whereas the bubble phases, especially with spherical bubbles, are at best expected at hypothetically low values of $Y_{\rm L}$ and/or $E_{\rm surf}$. For the rod-like and slab-like nuclei, thermally induced displacements are calculated from their respective elastic constants. It is found that at temperatures appropriate to supernova cores, thermal fluctuations would destroy the layered lattice of slab-like nuclei almost independently of the nuclear models and of the degree of the neutrino degeneracy.

Equilibrium phase diagrams for neutron star matter at subnuclear densities are obtained at zero temperature. Spherical, rod-like and slab-like nuclei as well as spherical and rod-like nuclear bubbles are taken into account by using a compressible liquid-drop model. This model is designed to incorporate uncertainties in the nuclear surface tension and in the proton chemical potential in a gas of dripped neutrons. The resultant phase diagrams show that for typical values of these quantities, the phases with rod-like nuclei and with slab-like nuclei occur in the form of Coulomb lattice at densities below a point where the system becomes uniform. Thermal fluctuations leading to displacements of such nuclei from their equilibrium positions are considered through explicit evaluations of their elastic constants; these fluctuations can be effective at destroying the layered lattice of slab-like nuclei in the temperature region typical of matter in the neutron star crust.

The influence of the presence of hyperons in dense hadronic matter on the quantum nucleation of quark matter is examined at low temperatures relevant to neutron star cores. We calculate the equation of state and the composition of matter before and after deconfinement by using a relativistic mean-field theory and an MIT bag model, respectively; the case in which hyperons are present in the hadronic system is considered, together with the case of the system without hyperons. We find that strangeness contained in hyperons acts to reduce a density jump at deconfinement as well as a lepton fraction in the hadronic phase. As a result of these reductions, a quark matter droplet being in a virtual or real state has its effective mass lightened and its electric charge diminished into nearly zero. The Coulomb screening of leptons on the droplet charge, which has significance to the droplet growth after nucleation in the absence of hyperons, is thus shown to be of little consequence. If the effective droplet mass is small enough to become comparable to the height of the potential barrier, the effect of relativity brings about an exponential increase in the rate of droplet formation via quantum tunneling, whereas the role played by energy dissipation in decelerating the droplet formation, dominant for matter without hyperons, becomes of less importance. For matter with and without hyperons, we estimate the overpressure needed to form the first droplet in the star during the compression due to stellar spin-down or mass accretion from a companion star. The temperature at which a crossover from the quantum nucleation to the Arrhenius-type thermal nucleation takes place is shown to be large compared with the temperature of matter in the core.

Rates for nucleation of two-flavor quark matter in a neutron star core, originally composed of nuclear matter in $\beta$ equilibrium, are calculated at zero temperature by a quantum tunneling analysis incorporating the electrostatic energy. We find that a nucleated droplet would develop into bulk matter due to electron screening effects.

Transitions of nuclear compositions in the crust of a neutron star induced by stellar spin-down are evaluated at zero temperature. We construct a compressible liquid-drop model for the energy of nuclei immersed in a neutron gas, including pairing and shell correction terms, in reference to the known properties of the ground state of matter above neutron drip density, $4.3 \times 10^{11} g cm^{-3}$. Recent experimental values and extrapolations of nuclear masses are used for a description of matter at densities below neutron drip. Changes in the pressure of matter in the crust due to the stellar spin-down are calculated by taking into account the structure of the crust of a slowly and uniformly rotating relativistic neutron star. If the initial rotation period is $\sim$ ms, these changes cause nuclei, initially being in the ground-state matter above a mass density of about $3 \times 10^{13} g cm^{-3}$, to absorb neutrons in the equatorial region where the matter undergoes compression, and to emit them in the vicinity of the rotation axis where the matter undergoes decompression. Heat generation by these processes is found to have significant effects on the thermal evolution of old neutron stars with low magnetic fields; the surface emission predicted from this heating is compared with the $ROSAT$ observations of X-ray emission from millisecond pulsars and is shown to be insufficient to explain the observed X-ray luminosities.