We investigate the cooling of neutron stars with relativistic and non-relativistic models of dense nuclear matter. We focus on the effects of uncertainties originated from the nuclear models, the composition of elements in the envelope region, and the formation of superfluidity in the core and the crust of neutron stars. Discovery of $2 M_\odot$ neutron stars PSR J1614-2230 and PSR J0343+0432 has triggered the revival of stiff nuclear equation of state at high densities. In the mean time, observation of a neutron star in Cassiopeia A for more than 10 years has provided us with very accurate data for the thermal evolution of neutron stars. Both mass and temperature of neutron stars depend critically on the equation of state of nuclear matter, so we first search for nuclear models that satisfy the constraints from mass and temperature simultaneously within a reasonable range. With selected models, we explore the effects of element composition in the evenlope region, and the existence of superfluidity in the core and the crust of neutron stars. Due to uncertainty in the composition of particles in the envelope region we obtain a range of cooling curves that can cover substantial region of observation data.

We study the thermal evolution of neutron stars in the presence of hyperons or kaons in the core. Our results indicate that the nucleon and hyperon direct Urca processes play crucial roles for the cooling of neutron stars. The presence of hyperons drives fast cooling mechanisms in two ways: 1) it allows the hyperon direct Urca prior to the nucleon direct Urca, 2) and it makes the nucleon direct Urca more feasible by reducing the neutron Fermi momentum. We found that the neutron star equation of state (EOS) with hyperons can be consistent with both mass and temperature observations. We also found that the neutron star EOS with kaon condensation can be consistent with observations, even though the cooling behavior is seldom useful to identify or isolate the effect of kaon condensation.

We consider the so called hyperon puzzle of neutron star (NS). We employ Skyrme force models for the description of in-medium nucleon-nucleon, nucleon-Lambda hyperon ($N\Lambda$), and Lambda-Lambda ($\Lambda\Lambda$) interactions. A phenomenological finite-range force for the $\Lambda\Lambda$ interaction is considered as well. Equation of state (EoS) of NS matter is obtained in the framework of density functional theory, and Tolman-Oppenheimer-Volkoff equations are solved to obtain the mass-radius relations of NSs. It has been generally known that the existence of hyperons in the NS matter is not well supported by the recent discovery of large-mass NSs ($M \simeq 2 M_\odot$) since hyperons make the EoS softer than the one without them. For the selected interaction models, $N\Lambda$ interactions reduce the maximum mass of NS by about 30~\%, while $\Lambda\Lambda$ interactions can give about 10~\% enhancement. Consequently, we find that some Skyrme force models predict the maximum mass of NS consistent with the observation of $2 M_\odot$ NSs, and at the same time satisfy observationally constrained mass-radius relations.

We investigate nuclear matter equations of state in neutron star with kaon condensation. It is generally known that the existence of kaons in neutron star makes the equation of state soft so that the maximum mass of neutron star is not likely to be greater than 2.0 $M_{\odot}$, the maximum mass constrained by current observations. With existing Skyrme force model parameters, we calculate nuclear equations of state and check the possibility of kaon condensation in the core of neutron stars. The results show that even with the kaon condensation, the nuclear equation of state satisfies both the maximum mass and the allowed ranges of mass and radius.

In this work, masses and radii of neutron stars are considered to investigate the effect of nuclear symmetry energy to the astrophysical observables. A relativistic mean field model with density-dependent meson-baryon coupling constants is employed in describing the equation of state of dense nuclear matter, and the density dependencies of the symmetry energies are quoted from the recent phenomenological formulae obtained from the heavy ion data at subnuclear saturation densities. Since hyperons can take part in the $\beta$-equilibrium of the dense matter inside neutron stars, we include hyperons in our estimation and their roles are discussed in combination with that of the nuclear symmetry energy.

The proton-proton fusion reaction, $pp\to de^+\nu$, is studied in pionless effective field theory (EFT) with di-baryon fields up to next-to leading order. With the aid of the di-baryon fields, the effective range corrections are naturally resummed up to the infinite order and thus the calculation is greatly simplified. Furthermore, the low-energy constant which appears in the axial-current-di-baryon-di-baryon contact vertex is fixed through the ratio of two- and one-body matrix elements which reproduces the tritium lifetime very precisely. As a result we can perform a parameter free calculation for the process. We compare our numerical result with those from the accurate potential model and previous pionless EFT calculations, and find a good agreement within the accuracy better than 1%.

The total cross section for radiative neutron capture on a proton, $np \to d \gamma$, is evaluated at big bang nucleosynthesis (BBN) energies. The electromagnetic transition amplitudes are calculated up to next-to leading order within the framework of pionless effective field theory with dibaryon fields. We also calculate the $d\gamma\to np$ cross section and the photon analyzing power for the $d\vec{\gamma}\to np$ process from the amplitudes. The values of low energy constants that appear in the amplitudes are estimated by a Markov Chain Monte Carlo analysis using the relevant low energy experimental data. Our result agrees well with those of other theoretical calculations except for the $np\to d\gamma$ cross section at some energies estimated by an R-matrix analysis. We also study the uncertainties in our estimation of the $np\to d\gamma$ cross section at relevant BBN energies and find that the estimated cross section is reliable to within $\sim$1% error.