Progress in Nuclear Astrophysics of East and Southeast Asia

Nuclear astrophysics is an interdisciplinary research field of nuclear physics and astrophysics, seeking for the answer to a question, how to understand the evolution of the Universe with the nuclear processes which we learn. We review the research activities of nuclear astrophysics in east and southeast Asia which includes astronomy, experimental and theoretical nuclear physics and astrophysics. Several hot topics such as the Li problems, critical nuclear reactions and properties in stars, properties of dense matter, r-process nucleosynthesis and $\nu$-process nucleosynthesis are chosen and discussed in further details. Some future Asian facilities, together with physics perspectives, are introduced.

osynthesis and ν-process nucleosynthesis are chosen and discussed in further details. Some future Asian facilities, together with physics perspectives, are introduced.
Keywords nuclear astrophysics · east and southeast Asia

Introduction
Nuclear astrophysics deals with astronomical phenomena involving atomic nuclei, and therefore it is an interdisciplinary field that consists of astronomy, astrophysics and nuclear physics. We seek for the answer to a question, how to understand the evolution of the Universe with the nuclear processes which we learn. Figure 1 shows a map of nuclides called nuclear chart, where our knowledge about nuclei and various nucleosynthetic processes are indicated.
The important questions in our field include but not limited to : -What is the origin of the elements in the cosmos? -What is the nature of neutron stars and dense nuclear matter? -What are the nuclear reactions that drive the evolution of stars and stellar explosions?
The research involves close collaboration among researchers in various subfields of astronomy, astrophysics, nuclear and particle physics. This includes astronomical observations using telescopes, gravitational wave (GW) detectors, and neutrino detectors; accelerator laboratory experiments using beams of stable or radioactive nuclei, neutrons, and gamma-rays; laboratory analysis of interstellar grains; large scale computer simulations of stellar explosions and nuclei; and theoretical work in nuclear physics and astrophysics. The interdisciplinary research not only answers the fundamental questions in nature but also provides unprecedented opportunities to discover new physics. For example, the decades of collaborative efforts of astronomers, astrophysicists, nuclear and particle physicists successfully established the standard solar model and led to a perfect solution to the long standing solar neutrino problem and a breaking through discovery of the new properties of neutrino beyond the standard model. A new era of astronomy has been opened by the historic multi-messenger observations of the SN1987A supernova (SN) event and the GW170817 neutron star merger (NSM). Precise knowledge of the critical nuclear physics inputs and reliable stellar models are urgently needed to decipher the messages correctly. Nuclear astrophysics community in east and southeast Asia has made tremendous progresses over the past few decades. Since early 1980's, the development of radioactive ion beam facilities (RIB) has opened a new frontier of nuclear astrophysics. Short-lived radioactive nuclei are expected to play the critical roles in explosive nucleosynthesis in stars, galaxies and the universe. Such an RIB facility is the laboratory for astrophysics where one can simulate celestial nuclear reactions of astrophysical interest to seek for the origin of atomic nuclides. One of the recent highlights is the successful approach to the r-process nuclei at RIKEN RI Beam Factory (RIBF). New facilities such as RAON and HIAF are being built to further expand the territory in the unexplored region in the chart of nuclei. The reaction rates of critical reactions, such as the holy grail reaction -12 C(α,γ) 16 O, are the major uncertainties of stellar models, preventing us from understanding the evolution of stars. Jinping Underground Nuclear Astrophysics (JUNA) collaboration developed an underground high-intensity accelerator facility to determine the critical reaction cross sections directly at the Gamow energies, and successfully accomplished the first campaign.
Remarkable progress in supercomputing power with huge memories enables astrophysicists to simulate catastrophically energetic events like supernova explosions or merging binary neutron stars. We predict a variety of nucleosynthetic products in these explosive episodes as well as quasi-static nuclear burning processes in stars. These theoretical predictions motivate spectroscopic astronomical observations of elemental abundances and their distribution in the universe through large scale medium-to-high resolution survey including LAMOST-II, 4MOST, WEAVE, etc. as well as high-dispersion spectrographs equipped on board of ground-based large telescopes like Subaru, Keck, VLT or TMT and the space telescope like Hubble Space Telescope. These astronomy missions are still in middle of the road to reach the oldest metal-free Population III stars that were made from the primeval gas from the big-bang nucleosynthesis (BBN) in the early universe. The ultimate goal is to fully uncover the evolution of the chemical composition in the Universe. Achieving this goal can further bring new insight to the fundamental issues regarding the evolution of dark matter halos.
In this article we first briefly review the nuclear astrophysics programs in five selected countries/regions in east and southeast Asia, and then report the progresses of some hot topics together with some highlights. The future prospect of this field is also provided.

Overview of Asian Activities in Nuclear Astrophysics
Five countries/regions, namely Japan, Korea, mainland of China, Taiwan and Malaysia, are chosen as the examples to showcase the progress achieved in east and southeast Asia.

Japan
Nuclear astrophysics in Japan has made remarkable progress since 1980s as a new interdisciplinary field, largely motivated by the development of the studies in nuclear physics with RI beams.
In nuclear physics, extensive studies of neutron halo, which was first discovered by a Japanese group [5,6] under INS-LBL (Institute for Nuclear Study -Lawrence Berkeley Laboratory) collaboration, were made at the accelerator facility of RIKEN. Search for new super-heavy elements which led later Fig. 1 Nuclear chart and the major nucleosynthetic processes in the universe. Most paths involve unstable nuclei, many of which are either difficult to be studied due to their short lifetimes and low beam intensities, or have not yet been discovered. 226 new isotopes (red boxes) have been synthesized and 84 nuclear masses (blue boxes) have been measured at Asian facilities and adopted in the Atomic Mass Evaluation by 2020 [1,2]. The nuclei shown in purple-red and light blue indicate those with production rates higher than 1 particle per day, given a primary-beam intensity of 1 particle µA at RIBF [3]. Isotopes in the yellow region are predicted based on the KTUY mass formula [4] but not synthesized yet. New facilities such as RAON, HIAF as well as the RIBF upgrade will provide new opportunities to synthesize more isotopes which are critical to the understanding of the nuclear processes in stars, galaxies and the universe.
to the discovery of Nihonium with atomic number Z=113 [7,8], discovery of appearance and disappearance of new and known magic numbers [3,9], and measurement of the half lives of r-process nuclei [10,11,12], respectively, were carried out in this institute. The INS of The University of Tokyo made continuous efforts in developing accelerator technology that would be realized later in HIMAC (Heavy Ion Medical Accelerator in Chiba) and J-PARC (Japan Proton Accelerator Research Complex). In particle physics, celestial neutrinos from SN1987A were detected in KAMIOKANDE, followed later by discovery of the evidence for neutrino-flavor oscillations in atmospheric neutrinos in SUPER-KAMIOKANDE. In astronomy, Subaru Telescope had the first light in 1990s, leading to discovery of several extremely metal-poor (EMP) stars which are enriched with r-process elements. In astrophysics, Japanese X-ray satellites HAKUCHO, TENMA, GINGA, ASCA, and SUZAKU identified new X-ray sources in space.
These discoveries were quickly shared and studied theoretically by qualified research scientists crossing over the neighboring interdisciplinary sub-fields (e.g. [13,14]). Such a coherent integration of activities has been supported by the bi-annual international conference "Origin of Matter and Evolution of Galaxies" (OMEG), which started in collaboration with INS (Institute for Nuclear Study) and RIKEN in 1988 and has been held in Asian countries until today. Together with OMEG, Japan Forum of Nuclear Astrophysics (JaFNA/UKAKUREN), established in 2008, supports world-wide nuclear astrophysics activities as well as China Institute of Nuclear Astrophysics (CINA), and Joint Institute for Nuclear Astrophysics -Center for the Evolution of the Elements (JINA-CEE).
In late 1980s, the developments of physics of unstable nuclei and the RI beam technology enhanced experimental nuclear astrophysics in Japan along with similar efforts in US and Europe. An ISOL-based low-energy RI beam facility was sought in INS, while RIKEN succeeded in providing in-flight type RI beams in early 1990s. Later RIKEN extended the facility to RIBF, the first "third generation" RIB facility, which started regular operation in 2007 and now approaches r-process nuclei. The Center for Nuclear Study (CNS) of University of Tokyo introduced a low-energy in-flight RIB facility CRIB (CNS Radioactive Ion Beam separator) in 2000 [15], which was intended to study directly the astrophysical reactions under explosive burning, and has been used to expand international collaborations especially with researchers from the east Asian countries. Another large-scale facility J-PARC also provides research opportunities for neutron-induced reactions of astrophysical interest and spectroscopy of hyper-nuclei which serves to establish the equation-ofstate (EoS) of neutron stars.
The first target of the experiments in 1980's was for explosive hydrogen burning, which takes place in novae and X-ray bursts [16]. The 15 O(α, γ) 19 Ne(p,γ) 20 Na reactions relevant to the breakout of the rp process were studied by the charged-particle spectroscopy [17]. The onset condition of the hot-CNO cycle was determined by studying the 13 N(p,γ) reaction by emplying the Coulomb dissociation method [18].
Reaction studies by neutron and photon beams were respectively developed by the groups at Tokyo Institute of Technology [19,20] and Konan University [21]. A group at Konan University built an experimental hutch called GACKO in the New SUBARU electron ring to study (γ, n) reactions relevant to p(γ)-process [22] and s-process [23,24]. In connection with the solar neutrino problem and the overproduction problem of BBN Li abundance, precise theoretical calculations were carried out for both reactions 3 H(α, γ) 7 Li and 3 He(α, γ) 7 Be in order to remove nuclear physics uncertainties relevant to these problems [25,26]. Also the 8 B production reaction 7 Be(p,γ) 8 B was investigated by the Coulomb dissociation at RIKEN [27,28,29] and the resonant scattering of 7 Be+p [30] at CRIB. An important stellar reaction 12 C(α,γ) was intensively studied by a group at Kyushu University by developing a sophisticated recoil mass separator [31].
Japanese astronomers are involved in various observations that provide a direct way to identify the site of nucleosynthesis. Recent SN detection is followed up by spectroscopic observations with, for example, the Subaru Telescope [32], which reveal the various types of SN explosions and identify the elements in their ejecta. Merging binary neutron stars, identified for the first time by the GW detection as GW170817 , could be another site of explosive nucleosynthesis. Follow-up multi-wavelength observations give a significant impact on our understanding of the r-process [33]. Recent observations also revealed that novae are the promising sites of 7 Li production [34].
Recent studies with wide-field survey telescopes like Skymapper (e.g., [35]) and large telescopes including the Subaru Telescope (e.g., [36,37]) found EMP stars having Fe abundances more than 10 5 times lower than the solar value. Such observations may constrain the mass distribution of first generations of massive stars, which is a long-standing unsolved problem in astronomy. Detailed chemical abundances of metal-poor stars are also applied to constrain r-process, s-process, and the BBN.
Recently the Gaia mission has precisely measured the position and motion of more than one billion stars. Combined with chemical abundance information for a portion of the samples, the formation process of the Milky Way substructures is being studied (e.g., [38]). Formation of the galactic halo structure should be related to chemical abundances of EMP stars. A group of stars with large excess of r-process elements (e.g., Eu/Fe ratio more than 10 times higher than the solar value) are of great interest. The recent discovery of the ultra-faint dwarf galaxy Reticurum II (Ret. II), whose member stars ubiquitously show enhancement of r-process elements [39,40], has made a hint to interpret such that these r-process-enhanced stars in the halo could have been formed in a small stellar systems affected by an r-process event, and then accreted into the Milky Way. Recent finding of r-process-enhanced stars in a wide range of metallicity in the Milky Way by a China-Japan collaboration team also shed new light on the interpretation of the the r-process.
Large progress in these studies is expected in near future. Detection of GW from the neutron star merger event has a large impact on nuclear astrophysics, in particular on the study of r-process. Further progress is expected by new facilities, KAGRA in Japan that started operation in 2020 and LIGO India that is under construction [41]. Following the large optical surveys of stars in the Milky Way including those with LAMOST, surveys of fainter objects are planned with the Subaru Telescope using a new instrument Prime Focus Spectrometer (PFS) from 2023. For the follow-up high-resolution spectroscopy for such faint objects, next generation extremely large telescopes are required. China, India, and Japan are the members of the Thirty Meter Telescope (TMT) and Korea and Australia are partners of the Giant Magellan Telescope (GMT).

Korea
Since the early 1990s, as the nuclear astrophysics community grew along with new advances of Rare Isotope (RI) accelerators around the world, the number of researchers in nuclear astrophysics has steadily increased also in Korea. During the last two decades, many young scientists carried out their researches in nuclear astrophysics at various accelerator facilities such as Center for Nuclear Study (CNS, Japan), Oak Ridge National Laboratory (ORNL, US), National Superconducting Cyclotron Laboratory (NSCL, US), Argonne National Laboratory (ANL, US), TRIUMF (Canada) and RIKEN (Japan). Many of these researchers currently play an important role in the construction of the Korea RI accelerator, RAON, as well as performing experiments at existing RI accelerator facilities around the world.
Measuring the 7 Be(p,γ) 8 B reaction rate is one of the most important experiments in nuclear astrophysics, because neutrino flux from the β-decay of 8 B, measured by various underground neutrino detectors to test the standard solar model and the neutrino oscillation, showed a correlation with the reaction rate. The first theoretical calculation on the contribution of s and d states of the 7 Be(p,γ) 8 B reaction was performed by Korean researchers [42], and motivated a nuclear astrophysics experiment using RI beams, the Coulomb dissociation of 8 B, investigating the 7 Be(p,γ) 8 B reaction [43].
While there has been no accelerator producing RI beams in Korea, the Korean nuclear astrophysics community explored many indirect studies of important reactions in astrophysics. Some examples are 17 F(p,γ) 18 [52], which are related to astrophysical processes including the primordial BBN, the pp chain, the hot CNO cycle, and the breakout reaction to the rp-process.
A close collaboration between Korean and Japanese nuclear astrophysics communities made experimental studies using RI beams at the CRIB. The elastic scattering excitation functions of 14,15 O(α,α), for example, were measured by bombarding radioactive oxygen beams onto the helium gas target to extract the alpha widths and excitation energies [50]. Both the important 15 O(α,γ) 19 Ne and 14 O(α,p) 17 F reactions proceed through resonance states above the 14,15 O+α thresholds in the explosive process. The 26 Si+p and 25 Al+p elastic scattering experiments [53,54] were also performed to study the 26 Si(p,γ) 27 P and 25 Al(p,γ) 26 Si reactions which are well known as the relevant synthesis mechanism of 26 Al gs .
The RAON facility [55] is under construction and will be the first RIB facility in Korea, providing exotic and unique radioactive isotopes for key experiments in many areas of nuclear physics and nuclear astrophysics. The stable beam commissioning is expected to begin in 2021 and RI beams will be produced by the Korea Broad acceptance Recoil Spectrometer & Apparatus (KoBRA) [56]. Moreover, the Center for Exotic Nuclear Studies, CENS, was recently launched at Institute for Basic Science (IBS). One of the main objec-tives of CENS is to carry out nuclear astrophysics at existing RIB facilities as well as to prepare future experiments at RAON.
RAON is also expected to boost active theory-experiment collaborations in nuclear astrophysics. Theory groups in Korea are performing researches on diverse areas to study the nature of nuclear systems and the properties of nuclei. With a chiral perturbation theory, for example, they could resolve the long-standing 10% discrepancy between theory and experiment in the radiative neutron-proton capture (n+p→d+γ) cross section at threshold [57]. Another examples are pionless and cluster effective field theories for electroweak reactions at low-energy, nuclear lattice and no-core shell model for ab initio description of nuclear matter and light nuclei.
Recently, a new potential dubbed Daejeon16 dedicated to the shell-model and a new energy density functional dubbed KIDS (Korea: IBS-DaeguSKKU) were developed [58]. Neutrino driven heavy-element synthesis as well as the Li problem in the BBN are also under active investigation.

Mainland of China
The nuclear astrophysics experimental activity in the mainland of China started at China Institute of Atomic Energy (CIAE) in 1990s. The major research facilities for nuclear astrophysics experiments include HI-13 MV tandem accelerator and Beijing Radioactive Isotope Facility (BRIF) -a ISOL facility driven by proton beam at CIAE, and Heavy Ion Research Facility at Lanzhou (HIRFL) , 320kV accelerator and Low Energy Accelerator Facility (LEAF) at Institute of Modern Physics (IMP), Chinese Academy of Sciences. The experimental groups have measured a number of critical reaction cross sections and nuclear properties using stable and radioactive ion beams. For example, the CIAE group has performed a number of indirect measurements of important reactions in stars, such as 7 Be(p,γ) 8 B, 13 C(α,n) 16 O and 12 C(α,γ) 16 O [59,60,61,62]. The IMP group precisely measured the masses of more than 30 short-lived isotopes in the region close to the proton drip line using a cooler storage ring (CSR) to study the explosive nucleosynthesis [63,64,65] and also performed a number of experiments to determine the critical reaction rates using the low energy stable or radioactive beams provided by HIRFL and the CRIB radioactive beam facility at Japan [66,67,68]. The IMP group is also in charge of the atomic mass evaluation (AME), an essential tool for the studies of nuclear physics and nuclear astrophysics [1,2]. The Jinping Underground Nuclear Astrophysics (JUNA) collaboration [69] established an underground accelerator facility for the direct measurements of critical nuclear reaction cross sections at stellar energies. The newly finished back-n white neutron facility at China Spallation Neutron Source (CSNS) [70] and the to-be-finished Shanghai Laser Electron Gamma Source (SLEGS) will offer new opportunities for the studies of neutron or gamma-induced reactions. High Intensity heavy ion Accelerator Facility (HIAF) are being built [71] to measure the critical reactions at stellar energies and study the critical nuclear properties. More details of JUNA and HIAF are available in section 4. Beijing ISOL, a powerful ISOL facility driven by a reactor, is being proposed [72,73,74] for the studies of the r-process nucleosynthesis.
Theoretical nuclear model calculations provide an indispensable approach in the region beyond the current experimental capability. For example, a database of neutrino spectra of the sd-shell nuclei under astrophysical conditions was established with the nuclear shell model [75]. The spherical relativistic continuum Hartree-Bogoliubov (RCHB) theory and mass formula were applied to predict the nuclear masses in the neutron-rich region which is inaccessible at any existing facilities [76,77]. The nuclear beta decay half-lives were systematically studied by fully self-consistent relativistic quasiparticle random phase approximation (QRPA) model [78], and the description of beta-decay half-lives was further improved by either considering tensor force or including many-body correlations beyond QRPA [79,80]. New sets of EoS, for example, Quark Mean Field [81], have been proposed for neutron stars through the confrontation of theoretical calculations with laboratory measurements of nuclear properties and reactions and increasingly accurate astronomy observations.
The Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAM-OST) is the largest Chinese optical telescope constructed in 2008. LAMOST survey [82,83] has finished five-year phase-I survey with the low-resolution mode (R ∼ 1, 800). Currently, it is running phase-II survey in a combination mode with both low-and medium-resolution spectra (R ∼ 7, 500). The LAM-OST survey has obtained over 10 million low-resolution spectra and six million medium-resolution spectra. In the 7th data release (DR7) of LAMOST, seven million stars are released for public use with reliable stellar parameters calculated through an upgraded pipeline [84]. Based on the LAMOST data and follow-up observations with other facilities including the Subaru telescope, Lijiang 2.4m and 1.8m telescopes, etc., great advance has already been made in the field of stellar and Galactic astronomy. The huge database of LAMOST has become an invaluable reservoir of studying rare objects leading to a number of important discoveries. For example, the most Li-rich K giant TYC 429-2097-1 has been discovered to have a Li abundance over 3,000 times higher than the normal amount and provides a great opportunity to investigate the origin and evolution of Li in the Galaxy [85], while Yan et al. [86] presents a uniquely large systematic study of lithium-rich stars and reveals that the distribution of lithium-rich red giants declines steeply with increasing lithium abundance, with an upper limit of around 2.6 dex, whereas the lithium abundances of red clump (RC) stars extend to much higher values; Xing et al. [87] reports the first extremely r-process enhanced (r-II) star exhibiting very low α-elemental abundance, which provides the clearest chemical signature of past accretion events onto the Milky Way; Liu et al. [88] discovers the most massive stellar black hole whose formation poses an extremely challenging problem to current stellar evolution theories. On the other hand, LAMOST also enables systematic studies on large sample of various types of stars in the Galaxy, e.g., the largest metal-poor star searching project, which has discovered over 10,000 candidates with metallicities less than 1% of the sun [89], thus forming the largest catalogue of bright source of metal-poor stars in the world which is very suitable for statistical studies of the Galactic halo and follow-up observation with the existing ground-based telescopes as well. LAMOST is now moving forward to another step to install a high-resolution spectrograph. Though still in the commissioning stage, it has already obtained scientific progresses such as discovering two new Li-rich RC stars [90]. Besides the conventional optical observation, the observation of INSIGHT-hard X-ray modulation telescope helped to confirm the unexpected weak and soft nature of GRB170817A. By collaborating with the MPIE group, astronomers also explore the γ-rays from radioactive isotopes in the Milky Way to probe the supernova progenitors and later evolution of massive stars. A 20kt neutrino observatory (JUNO) will be ready in 2020 to detect all six flavors of neutrinos and anti-neutrinos emitted from the supernova core [91,92,93,94,95].
The stellar modeling program is being developed in the mainland of China. The groups at Yunnan observatories have been working on the progenitors of SNe Ia since 2006. Their detailed 1D stellar evolution calculations suggested that SNe Ia generated from Si-rich WDs and/or C-O-Si WDs are expected to show some particular characteristic features such as the high velocity of Si lines in their spectra. This can be used as the 'smoking gun' for identifying the SN Ia progenitor systems composed of a CO WD and a He star [96]. They are also working on multi-dimensional hydrodynamical simulations of SNe Ia such as the interaction of SN Ia ejecta with a stellar companion star with a hope to predict observables which can be used to compare with the SN Ia observations with different telescopes [97]. This is expected to place strict constraints on SN Ia progenitor models and their explosion mechanism.

Taiwan
Earlier efforts in nuclear astrophysics focused on measurements of short-lived radioactive nuclei in meteorites and understanding their origin [98,99,100,101,102,103,104,105], as well as on experimental side relevant to the nuclear reaction rates in astrophysical environments, e.g., [106,107,108,109]. Since 2018, researches focusing on the theoretical aspects of supernova explosions and the nucleosynthesis of heavy elements have become active in various leading institutes, including the Institute of Physics, Institute of Astronomy and Astrophysics at the Academia Sinica, and the Institute of Astronomy at the National Tsing Hua University. Moreover, as nuclear astrophysics plays a pivotal role in various phenomena that are the prime targets in the newly-opened era of multi-messenger astronomy, collaborative efforts among domestic researchers are under development supported by the National Center for Theoretical Sciences.
As for supernova-related studies, efforts in implementing the efficient neutrino transport scheme -Isotropic Diffusion Source Approximation (IDSA)were made in multidimensional simulations of core-collapse supernova explosions with the open-source code FLASH [110,111,112,113]. With such develop-ment, the dependence of the black-hole forming supernovae on different nuclear EoS and the stellar rotation rates was examined and the expected GW signals were computed [114,115]. Guo, Wu and their collaborators then improved the treatment and consistency of the neutrino interactions in dense matter, including the inverse neutron-decay, muonic charged-current interactions, and nucleon-nucleon Bremsstrahlung [116,117,118,119], which allows to quantify the impact of these interactions on supernova simulations. Meanwhile, in a series of work [120,121], potential explosions of massive stars of ∼ 35 − 50 M triggered by the hadron-quark phase transition were suggested. Such exotic explosions were found to result in unique signatures in both the neutrino signals and the r-process nucleosynthesis yields that can be tested by future neutrino detection and metal-poor star observations. Chen et al. [122,123,124] performed detailed multidimensional long-term simulations of the pair-instability supernovae and magnetar-driven supernovae on even more massive and energetic end. These simulations aim to study the important consequences due to multidimensional effects on the observables.
Regarding the r-process nucleosynthesis in NSMs, various unique signatures of the r-process, including the late-time kilonova lightcurves and the γ-ray emission from merger remnants residing in the Milky Way, were explored [125,126,127]. Future detection of these can shed new lights on our understanding of the nature of r-process in NSMs. Banerjee, Wu and Yuan examined the important question of whether NSMs can be the dominant rprocess sites in the Milky Way [128], which pointed out that the late-time galactic chemical evolution of r-process enrichment can be consistently accounted for by NSMs when considering new effects of neutron star kicks and the inside-out evolution of the Milky Way which was omitted in previous studies and easing the need of additional r-process sources. In parallel, several authors made efforts in studying the onset conditions of neutrino flavor conversions in NSM remnants and their impact on r-process [129,130,131]. Improved numerical modeling of this phenomenon is on-going to further solidify the roles of neutrinos in NSMs.

Malaysia
Nuclear astrophysics and astrophysics as a whole is in its infancy stage of developing into a serious discipline in Malaysia where it started in earnest in the early 1990s on the standard solar model and solar neutrinos. The research focuses on the theoretical aspects of nuclear astrophysics. Currently, University of Malaya and International Islamic University Malaysia are involved in this field.
Amongst the topics of interest in nuclear astrophysics is the study of thermonuclear reaction rates in stellar evolution [132] and on the reliability of the M3Y double folding potentials (DFM) and the α-α double-folding cluster (DFC) potentials in fusion cross sections of light systems of heavy ion reactions. Both potentials have been shown to reproduce the experimental fusion cross section data reasonably well [133]. In fact, the DFC potential reproduces satisfactory fitting to the astrophysical S-factor data and consistent prediction of the astrophysical reaction rates.
In the light of observations of very massive stars (VMS) in the local universe [134], theoretical models were constructed to verify this observation [135]. The focus of this work is mainly on the study of life and death of very massive stars, which are extreme cosmic engines, and enriching their environments with chemically processed material throughout their entire life-time in Ref. [135,136,137]. This work is in collaboration with Keele Astrophysics and Geneva stellar evolution group where both research groups are actively involved in the study of the evolution of massive stars. The current project is a study on the grids of stars with supersolar metallicity which would be useful for studying massive stars in the metal rich environment.
In the area of neutrino astrophysics, the stopping power of matter on neutrinos with oscillation effects is applied to cosmological neutrinos [138] and neutrinos in very massive stars. Recently, neutrino emissivity from thermal processes in VMS in our local universe is investigated. Neutrino luminosity towards the end of the VMS life is higher than in massive stars. This could shed light on the possibility of using the detection of the neutrinos to locate the candidates for the pair instability supernova in our local universe.
The first detection of binary neutron star merger event GW170817 in 2017 by LIGO and Virgo [139] has established limits both on the tidal deformability of the canonical star to be constrained in the range of 70 ≤ Λ 1.4 ≤ 580, and the chirp mass M = 1.188M [139]. The new constraints from binary neutron star merger in terms of tidal polarization open the possibility to revisit the wellestablished model, namely, density-dependent van der Waal model, Skyrme-Hartree-Fock, momentum dependence interaction, and relativistic mean field approaches. Further it is crucial to establish the possible correlations and sensitivity of important nuclear bulk properties, i.e., the symmetry energy, its slope, compressibility, and values of the tidal deformability for the canonical 1.4M [140,141,142,143]. A connection has been established between the finite nuclei and infinite nuclear matter in terms of isospin asymmetry parameter i.e., symmetry energy and its coefficients using relativistic mean-field formalism [144,145]. This result is of considerable importance since due to shell closure over the isotonic chain, will act as an awaiting point in nucleosynthesis of r-process and experimental investigations towards the exotic region of the nuclear chart including the superheavy island.

SPECIAL TOPICS
Here we discuss several hot topics that exhibit a remarkable progress in nuclear astrophysics on synergy of the three fields of astronomy, nuclear and particle physics, and astrophysics.

Li problems
The Li production in BBN is a longstanding unsolved problem [146]. The standard homogeneous BBN model with cosmological baryon-to-photon ratio η determined by the observations of Cosmic Microwave Background (CMB) anisotropies and precisely measured neutron half-life predicts about a factor of three larger Li abundance than those measured in metal-poor mainsequence turn-off stars, while the other light nuclides 2 H, 3 He and 4 He are in reasonable agreement with observational constraints. The discrepancy of 7 Li is remarkably larger than error sources including non-LTE and 3D effect of stellar atmospheres in calculation of spectral line formation, and is called "Big-bang Li problem" Recent observations extended to the lowest metallicity range ([Fe/H]< −3) suggest systematically lower Li abundances than those at [Fe/H]> −3 [147]. This is a stringent constraint on the mechanism that causes the Li problem.
A number of suggestions have ever been proposed to solve this problem. One is to understand better the effects of diffusive transport and stellar massdependent convection and depletion of Li near the stellar surface of the metalpoor halo stars [148,149]. However these astration effects are hard to be fully understood.
This problem has been extensively investigated in nuclear physics laboratories. Although the production reactions 3 H(α, γ) 7 Li and 3 He(α, γ) 7 Be(e − ,ν e ) 7 Li are reasonably well studied both experimentally [150,151] and theoretically [25,26,152], there are some reactions not well studied yet which might alter the 7 Li abundance. The CIAE group has determined the cross sections of 12 reactions involving lithium [153]. But the BBN calculations showed that the new reaction rates have minimal effect on the final lithium abundances. In fact, 7 Be is much more abundantly produced than 7 Li in the standard BBN with the current accepted η-value. Therefore, the main focus has been placed on the experimental studies of destruction reactions of 7 Be, i.e. 7 Be(n,p) 7 Li and 7 Be(n,α) 4 He. The JAEA experiment [154] of the former reaction reported significant contributions of the decay brunches to the first excited state of 7 Li from the relevant 8 Be resonances, which was not identified in the previous n-ToF experiment at CERN [155]. Hayakawa et al. [156] have successfully measured both transitions, for the first time, at the wide energy range relevant to BBN, applying Trojan Horse Method (THM) by the use of intense 7 Be beam from CRIB at CNS. They also deduced the reaction cross sections for 7 Be(n,α) 4 He as well. This reaction was studied at RCNP [157] by the use of time reverse reaction 4 He(α,n) 7 Be and the p-wave dominance was found. These precise experimental data proved that both destruction reactions of 7 Be are not responsible for the factor of three overproduction in 7 Li abundance.
The 7 Be(d,p) 8 Be and 7 Be(d, 3 He) 6 Li destruction reactions were also studied at RCNP and CIAE, respectively. Their contributions were found to be quite minor [158,159]. An important report from Florida group [160] is that the 7 Be(d,α) reaction might make a remarkable effect on destruction of 7 Be. All other reactions however need to be revisited before making a definite conclusion.
Many other models beyond the standard BBN have also been proposed such as those including non-extensive Tsallis statistics [161], baryon inhomogeneous BBN [162,163,164], supersymmetric particle decay [165,166], sterile neutrinos [167], and so on. Among these, Tsallis statistics is one of the interesting models to solve the big-bang Li problem, and its theoretical establishment has started by taking account of the inhomogeneous primordial magnetic fields [168] still satisfying the observed constraints from CMB [169].
Post-cosmological origin of Li has been studied theoretically in chemical evolution models [170]. Tajitsu et al. [171] have reported the detection of 7 Be, which decays to 7 Li in 53 days, in the ejecta of a nova explosion (Fig. 2). This first discovery was followed by further detection of 7 Be in other novae [34,172] in addition to detection of 7 Li [173]. These observations indicate that nova explosions could be the major production sites of 7 Li in the present Milky Way as predicted theoretically.
Another unsolved question in astronomy is extremely enhanced Li abundances as high as A(Li) = log(n Li /n H ) + 12 > 1.5 in ∼ 1% of low-mass giant stars (e.g. [174]). 7 Be is produced by the 3 He(α, γ) 7 Be reaction in the hydrogen burning shell, conveyed to the stellar surface, and eventually decays to 7 Li, which is called Camelon-Fowler (CF) mechanism [175]. However, since 7 Li is easily destroyed by the 7 Li(p,α) reaction at low temperature, such a high 7 Li abundance level is not explained in the standard theory of low-mass stellar evolution. Recent studies have revealed that most Li-rich objects are found in the core He-burning phase (clump stars) [86,176]. There is a peculiar scenario such that extra mixing is induced by the tidal interaction with companion star and drives Li production. Another scenario is to assume that substellar objects holding high Li abundance were engulfed to form the stellar atmosphere of these stars. It has recently been pointed out that the produced 7 Be is boosted to the stellar surface due to more efficient themohaline mixing if neutrinos had finite magnetic moment to accelerate stellar cooling [177]. These proposed mechanisms of Li enhancement are still under debate.

Critical nuclear reactions and properties during quiescent and explosive burning
Critical reactions in stars are essential parts of the stellar models to understand multi-messages from stars, such as the isotopic/elemental abundances, light/neutrino fluxes or even GW of stars. However, the determination of the very low nuclear reaction rates in stars, and the determination of properties and reactions of very neutron deficient or very neutron rich unstable nuclei are two long standing challenges in nuclear astrophysics [178]. We review some progresses in the precise measurements of the critical nuclear reactions and properties. The r-process nucleosynthesis and the weak interaction process are discussed separately in sec.3.4 and sec.3.5.

Helium and carbon burnings
The 12 C(α, γ) 16 O reaction is the key reaction in the stellar helium burning stage. Competing with 3α process, the 12 C(α, γ) 16 O reaction rate determines the ratio of carbon and oxygen in the universe and affects many other stellar scenarios such as the synthesis of 60 Fe and 26 Al, supernova explosion and the formation of back holes [179]. 12 C(α, γ) 16 O is quoted as the Holy Grail reaction and a reaction rate with accuracy better than 10% is required [180] .
This reaction has been studied directly in the ground level labs using normal and inverse kinematics, respectively, in Japan [31,181]. However, the needed cross sections at around 300 keV is far yet to be studied in the coming years [181]. The CIAE group performed an indirect measurement of the S E2 (300) factor. Comparing to the value constrained by the direct measurement, the uncertainty of S E2 (300) factor is reduced from 56% to 10% , as shown in Fig. 3. The Kyushu group initiated a new program to revisit the nuclear reactions of the He burning using the quality beam from a Tandem accelerator. JUNA collaboration began the direct measurement in the Jinping Underground Laboratory in early 2021(See section 4.3). Another direct measurement of the reaction using an active target TPC is proposed in Korea as well [182]. Important progresses have been made for the triple alpha process in the last decade. Although a non-resonant mechanism was suggested to dominate the process with CDCC method [183], the detailed rate has been derived theoretically [184], which turns out to be inefficient to change the evolution of low-to-intermediate mass stars [185]. The triple alpha process at high temperatures are crucial to explosive burning, such as XRB and supernova explosion [186]. Resonances of 0 + and 2 + were observed at around 10 MeV [187] at RCNP. The 2 + resonance, which is likely a rotational member of the Hoyle state, was confirmed [188]. Recently, this width of the 9.6 MeV 3 − resonance was studied experimentally at RCNP, determining the total decay width for the first time. Their preliminary results shows a reduction of the reaction rate.
The 12 C+ 12 C fusion reaction is the primary reaction in the hydrostatic and explosive C/Ne burnings in massive stars. It is also important for the stellar explosions such as Type Ia supernova and superburst [189,190]. Determination of the 12 C+ 12 C fusion cross section at stellar energies below E c.m. = 3 MeV has been a long standing problem since the first measurement at 1960. While the 12 C( 12 C,n) 23 Mg reaction has been measured directly at Gamow energies [191], the lowest measured energy for the α and p channels has been pushed only down to E c.m. = 2.1 MeV with large uncertainties. Three different kinds of extrapolations have been used to estimate the reaction cross section at lower energies [192,193,194,195,196,197,198]. A recent 24 Mg(α,α ) 24 Mg experiment at RCNP provides important information for the resonances at stellar energies [199]. Extrapolating theory is being developed at Malaysia and Japan. The IMP group is carrying out a direct measurement using the intense carbon beam from the Low Energy Accelerator Facility(LEAF) [200] while the LUNA collaboration is planning an underground measurement [201]. High intensity heavy ion beam, efficient detectors with ultra-low background, theoretical studies of the molecular resonances and fusion reaction theory are needed to achieve a more reliable reaction rate for the astrophysical application.

Explosive Hydrogen burning in Novae, X-ray bursts and Type II Supernovae
Novae and type I X-ray bursts (XRBs) are the most frequently observed explosive stellar events. The observed light curves and isotopic/elemental abundances are important probes to study the compact stellar objects and the related nucleosynthesis. Due to a large systematic uncertainty of the theoretical results from estimated nuclear physics inputs such as nuclear mass and reaction rate [202], experimental knowledge of the constituent reactions that power these violent outbursts are required to improve the uncertainty.
Explosive nucleosynthesis involving the HCNO cycle and the rp-process (see Fig. 1) is thought to be responsible for the production of terrestrial 15 N and the excess 22 Ne seen in many meteorites, as well as the elemental overabundances of O, Ne, Mg, etc. observed in nova ejecta. In order to understand the dynamics of such explosions and the origin of our solar system materials, it is necessary to determine the detailed characteristics of the important nuclear reactions in the explosive process. Important (α,p) reactions on 14 O, 18 Ne, 22 Mg, and 30 S, 15 O(α,γ) 19 Ne and (p,γ) reactions on 13 N, 18 F and 19 Ne were studied with radioactive ion beams at CRIB, RIKEN-RIPS, HIRFL-RIBLL1 and HRIBF [15,50,66,203,204,205]. However, most results are yet limited, and thus further studies are need.
The rp-process is a dominant mechanism of nucleosynthesis in XRB which is ignited in low mass binary systems consisting of a neutron star and a H/He rich companion star (See Fig.1). When the critical temperature and density are reached during XRB, the rp-process is triggered by the 3α-reaction and subsequently powered by a sequence of (p,γ) and (α,p) reactions. This process burns materials in the star into the calcium region. After the process, proton capture and β + -decay occur and produce the nuclei along the proton drip line up to A=100. It has been long considered that 64 Ge, 68 Se, and 72 Kr nuclides are the major waiting-point nuclei along the path of the rp-process. The predictive power of the XRB model is limited by the nuclear uncertainties around these waiting points which can only be removed with precise measurements. Taking the 64 Ge as an example, it was pointed out that the uncertainties of the mass and the proton capture rate of 65 As limited the predictive power of the XRB model. The mass uncertainty was eliminated by a precise mass measurement of 65 As performed using the HIRFL-CSR [65]. An experiment has been performed at RIKEN RIBF using the SAMURAI spectrometer, which measured 65 As and proton decaying in flight from 66 Se produced in one-neutron removal reaction induced by RI beams of 67 Se. The decay information was used to investigate their roles in the proton capture reaction on 65 As which is involved in the breakout of the rp-process starting from one of the waiting point nucleus 64 Ge in type-I XRBs.
The νp-process in type II supernovae might run a little closer to the line of stability because of presence of a small fraction of neutrons which discard the waiting points of the rp-process [206]. The neutrino oscillations, which is discussed in Sec 3.5, the critical nuclear reactions and nuclear properties are essential for the model development. A recent precision mass measurement at HIRFL indicated that the rp-process may go through the Zr-Nb region, suggesting less feasible to have a previously reported Zr-Nb cycle at 84 Mo [63]. 56 Ni(n,p) 56 Co has been identified as the critical reaction which impact the nucleosythesis [206]. This reaction will be studied indirectly by slowing-down beams with a novel technique in the OEDO (Optimized Energy Degrading Optics for RI beam) project, which is installed in RIKEN RIBF by CNS [207].
Recently, astronomers reported observation of a bump near the K edge of nuclides of Z=48-49 in the XRB [208], although the resolution was not good. High-resolution X-ray observations are really awaited in the future.

Radioactive tracers of the galaxy
The long-lived isotopes 26 Al(T 1/2 =0.72 Myr) and 60 Fe(T 1/2 =2.6 Myr) are produced mainly by massive stars. The nucleosynthesis of the long-lived radioactivities is an important constraint on the stellar models. The abundances inferred from gamma-ray astronomy may have important implications for rotationally induced mixing, convection theory, mass loss theory, the initial mass function for massive stars, and the distribution of metals in the galaxy. The present gamma-ray observations determine the 60 Fe/ 26 Al gamma-ray flux ratio to be 0.184±0.042 based on the exponential disk grid maps [209], as shown in Fig.4. But the mean ratio may vary in the range of 0.2 to 0.4 by using different sky maps. Theoretical groups predict the gamma-ray flux ratio at a large variation from 0.1 -1 [210,211,212,213,214]. As the major source of the large uncertainty comes from the reaction rates of the critical reactions, such as the neutron source reactions and those creation and destruction reactions for the unstable isotopes 26 Al and 60 Fe which have not be measured adequately, it is impossible to fully constrain the stellar model with the observed 60 Fe/ 26 Al gamma-ray flux ratio. Fig. 4 60 Fe/ 26 Al flux ratio for the grid of exponential disk models (blue). Including the uncertainties of the fluxes from each spectral fit, the total estimated 60 Fe/ 26 Al flux ratio from exponential disks is 0.184±0.042 [209].
Experimental and theoretical studies have been being carried out to identify the critical reactions and improve their accuracies with various techniques. Nucleosynthesis calculation identified the 59 Fe stellar decay rates as one major uncertainty in the model [215]. This uncertainty has been eliminated recently by an IMP-NSCL collaboration. Their experimental result shows the 60 Fe yield would be reduced by using the new rate, alleviating the tension between some model predictions and the observation [216]. The 59 Fe(n,γ) 60 Fe was restudied by the CIAE and JAEA collaboration using the surrogate method. The 26 Si(p,γ) 27 P and 25 Al(p,γ) 26 Si reactions are well known as the relevant synthesis mechanism of 26 Al gs . The 26 Si+p and 25 Al+p elastic scattering experiments were performed by using CRIB to study the abundance of 26 Al. Six resonance states were observed in the 26 Si+p elastic scattering experiment. The reactions with the isomer state of 26 Al in the stars complicates the nucleosynthesis of 26 Al. Fortunately 26 Al isomer beam becomes available at CRIB and the BRIF facility at CIAE for experiments. 25 Mg(p,γ) 26 Al is the important production reaction for 26 Al. JUNA carried out a direct measurement to determinate the strengths and ground state feeding factors of the crucial resonances.
Another long-lived isotope, 244 Pu(T 1/2 =80.0 Myr), are found in the deep sea sediment together with 60 Fe using AMS [217]. The current observation suggest that it comes from a rare event, such as neutron star merger. The consist understanding of the galactic radioactivity and deep sea sediment would help us understand the origin of the elements heavier than iron discussed latter in Ref. 3.4. It is also interesting to further study the moon samples from the Apollo or Chang-Er missions for finding the clues of the impact of SN events on the solar systems.

Neutron star and EoS of high density nuclear matter
Neutron star(NS) is the smallest and densest stellar object composed of nuclear matter. It is kept in hydrostatic equilibrium only by the pressure produced by the compressed nuclear matter. Nuclear matter in neuron star varies from very low density on the surface to several times the saturation density in the core. The property of nuclear matter is often expressed in terms of the EoS. The symmetry energy term (L) of EoS is one of the key for the high-density nuclear matter. It determines the important NS properties such as pressure, structure, radius, tidal deformability (Λ) and so on (Fig.5). The EoS of NS is also a fundamental component in the models of X-Ray Burst (XRB) and the neutron star merger (NSM), a plausible site for the r-process nucleosynthesis. Thus understanding the NS EoS of dense matter is one of the fundamental quests of both nuclear physics and astrophysics.
Though the nuclear matter cannot be directly examined in laboratories, to indirectly explore the nature of nuclear interactions in neutron-rich situations in finite nuclei is one of the target of experimental nuclear physics. On-going efforts at RIKEN RIBF and RCNP respectively aim at determining the symmetry energy term of EoS from different approaches: studies of multi-particle emission in collisions between neutron-rich nuclei and evaluation of the neutron skin thickness through giant resonance involving proton-neutron exchange. By colliding rare isotope Sn beams with isotopically enriched Sn targets, SπRIT collaboration deduced the slope of the symmetry energy to be 42 MeV< L < 117 MeV [218], which is slightly lower but consistent with the L-values deduced from the measurement of the neutron skin thickness of 208 Pb. Isolated multi-neutron systems and multi-neutron cluster states in neutron-halo nuclei are also investigated at RIBF. With the heavy ion beam covering sub GeV/u range, HIRFL-CEE is being built at IMP to study EoS at about 2 times saturation density via multi observable [219,220]. The inner core of neutron stars is predicted to be composed of hyperon-mixed matter. Its properties can be studied by the nature of hyper nuclei at J-PARC. New opportunity will be available soon at GSI and HIAF by studying the hypernuclear spectroscopy with Heavy-Ion Beams [221].
Lots of work have been done to connect the terrestrial experimental studies consistently with the observational results, like the mass constraint of the most massive pulsars, the tidal deformability constraints from GW170817 by LIGO/Virgo, and the simultaneous estimation of the mass and radius of PSR J0030+0451 by NICER, also incorporating our best knowledge of neutron matter based on ab initio calculations. Many efforts have been made on one of the EOS's main features, i.e., the symmetry energy, which is crucial for interpreting many astrophysical observations related to NSs and NS binaries. The future measurements of the radii of massive pulsars, using X-ray missions like NICER and eXTP, should provide even better constraints on the EOS [222].
In the future, it is important to establish new quantitative results of the EoS that are testable by experiments/observations to understand the nuclear interaction better and to probe the particle degree of freedom in cold, dense matter. It is also helpful to find further opportunities to study EoS from the dynamical processes of NS systems, like NS cooling, pulsar glitches, and short gamma-ray bursts.

r-process
Origin of radioactive nuclei such as 226 Ra (half-life τ 1/2 =1600y), 232 Th (τ 1/2 =140.5Gy), 238 U (τ 1/2 =44.7Gy) is a long-standing unsolved question in modern science since the discovery of radioactivity by Antoine H. Becquerel in 1896 and the extraction of radioactive substances by Marie Curie in 1898. These atomic nuclides are now thought to be produced in rapid neutron-capture process (rprocess) in extremely neutron-rich environment of explosive phenomena [223]. 80 100 120 Fig. 6 A comparison between the estimated r-process abundance curves with (red) and without (green) taking account of the new experimental data from RIBF [10]. Theoretical abundances are calculated in a canonical SN nucleosynthesis model assuming (n,γ) -(γ,n) equilibrium [224].
Binary neutron star mergers (NSMs) as well as core-collapse supernovae (CC-SNe) and collapsars are the most viable candidates for the origin of r-process elements among many other possible astrophysical sites [224,225]. CCSN includes both neutrino-driven wind (ν-wind) and magneto-hydrodynamic jet (MHDJ) SNe, which leave neutron star as a remnant, and collapsar is a very massive single star collapsing into a black hole [226,227]. Three different aspects are folded in the studies of r-process nucleosynthesis. They are the nuclear properties of extremely neutron-rich short-lived radioactive nuclei, the theoretical modeling of SN explosion, dynamics of NSM and nucleosynthesis, and the spectroscopic observations of r-process elements in kilonovae caused by neutron star mergers as well as in metal-deficient stars to test these theoretical predictions.

Highlights from Experimental Nuclear Physics
The r-process reaction flow path runs on extremely neutron-rich unstable nuclei far from the line of nuclear stability. Recent innovative experimental technique in nuclear physics has reached the production of some of these exotic nuclei. The fission of 238 U is a possible powerful tool for producing very neutron-rich nuclei close to or on the r-process pathway. RIBF has been developing a high-energy and high-intensity 238 U beam, and various secondary beams of unstable nuclei through in-flight fission allows one to study the rprocess nuclei.
The large amounts of lifetime data for β decay have been taken at RIKEN RIBF [10,11,12]. They now cover the r-process pathway below and beyond the N =82 magic number. The impact of the RIKEN data is demonstrated in Fig. 6 in a canonical SN nucleosynthesis model, where a network calculation assuming (n,γ) -(γ,n) equilibrium using RIKEN data available in 2015 is compared in two cases with and without the new experimental lifetimes [10]. One can see the difference clearly. Data at around another waiting point N =50 also are available [228].
In neutron-rich nuclei, the β decay daughter can be in its high lying excited state that emits one or more neutrons, which influences the final abundances. The BRIKEN project has started to measure the neutron emission probability P n , and a first result has already been obtained for a decay of 86 Ge affecting the decay paths to the production of 84 Kr, 85 Rb and 86 Kr.
Nuclei with neutron magic numbers, which are called "waiting points", are expected to play an important role in forming the r-process abundance peaks. Spectroscopic studies at RIBF could access the "doubly-closed" 78 Ni with Z=28 and N =50 for the first time. 78 Ni is a typical seed nucleus for the r-process in some stellar conditions. Its low lying states have been identified by measuring deexcitation γ rays from fast moving 78 Ni nuclei produced by a few nucleon removal reaction in inverse kinematics, where RI beams hit a stable target [229]. Together with lifetime measurements of nuclei in the vicinity of 78 Ni [228], the closed shell nature was found to persist with no clear signal for "shell quenching" [230,231].
Another neutron-rich region of interest is around N =82, responsible for the 2nd peak of the r-process abundance. Nature of palladium isotopes around the r-process path have been studied by measuring γ rays from their excited states populated either by isomer production [232] or nucleon removal reaction [233]. Measured systematic trend of the low-lying state indicates no significant shell quenching for N =82.
Using the low-energy RI beams in the OEDO project, (d,p) reactions as surrogate process of corresponding (n,γ) reactions on r-process nuclei will be studied under the a collaboration involving Asian institutes such as CNS, RIKEN, IMP, and SKKU. The region around 132 Sn with N =82 is a target. The feasibility of the method was already tested successfully for 79 Se(n,γ) at around few MeV of neutron energy.
The KISS project and its future extension KISS-II aims to approach the blank-spot or the "south-east" region of 208 Pb and further "east" beyond N =126, which includes the progenitor parent nuclei for the r-process 3rd peak element. The group of Wako Nuclear Science Center (WNSC) of KEK found the deep-inelastic collisions or multi-nucleon transfers can be more efficient to reach that region compared with the projectile fragmentation or fission usually employed in most of the RI-beam facilities. The project has started by building a system [234].
As illustrated in Fig. 1, fission fragment distribution (FFD) as well as various fission modes, i.e. spontaneous, β-delayed and neutron-capture induced fissions, are known to affect strongly the r-process abundance in NSM and collapsar nucleosynthesis [235,236,237]. It is planned to systematically measure the fission barrier of the neutron-rich actinoids with the SAMURAI setup at RIKEN RIBF. Nuclear data groups at Japan Atomic Energy Agency and Tokyo Institute of Technology successfully describe the FFD and their kinematical properties for long-lived radioactive actinoids and SHEs in a unique method of solving 3D/4D Langevin equation [238,239].
Measurement of nuclear masses is essential to study the nuclear structure date relevant for r-process [240,241]. In addition to the MR-TOF method in the SLOWRI setup for slow RI ions and the attempts to measure TOF in a long flight path, a novel method is being developed. Individually injecting short-lived nuclei are produced and identified in BigRIPS and their masses are determined by realizing the isochronous condition of the Rare RI Ring. The goal is to achieve 10 −6 mass precision for nuclei with a lifetime as short as 1ms, which should greatly improve our knowledge on the r-process network.

Highlights from Astronomical Observations and Theoretical Progress
Follow-up observations of optical-to-infrared light from GW170817, including the contributions of the observations with the Subaru Telescope lead by the J-GEM consortium [242], have provided an evidence that some sort of radioactive nuclei have been produced in binary NSM. The infrared excess over the optical spectra at several days was interpreted as the signature of large opacity due to the lanthanoids. Theoretical modeling of NSM r-process suggests that both the lighter and heavier r-process elements could be produced in the early dynamical ejecta and/or the later wind outflow from the remnant [243,244,225]. However, recent careful re-analysis of the spectra [245] and many other related studies (e.g., [246,247]) have indicated that the amount of lanthanoids is not as large as that expected in theoretical calculation for kilonova associated with GW170817. Meanwhile, the detection of the infrared emission at several tens of days has not yet been able to provide evidence for the presence of 254 Cf or α-decay nuclei. However, such association can be possible with future events [126,248].
Although theoretical model uncertainties still remain, there are large parameter space to make variations of the NSM r-process as a possible astrophysical site of the r-process [225,235,244]. The basic difference between NSM and CCSN is a different emergent event rate and a delay timescale from star formation of their progenitors until the time when r-process occurs. SNe and collapsars are dying single massive stars that culminate their evolution in explosions in a few My. Therefore, they can enrich r-process elements in the interstellar medium (ISM) from the early Galaxy. On the other hand, the majority of NSMs may not contribute to the early Galaxy because they can take cosmologically long timescale due to extremely slow GW radiation to lose kinetic energy of the orbital motion until they merge. The coalescence timescale estimated from observed double pulsars ranges from a few hundred My to the cosmic age or beyond [249] although faster mergers are also possible. Their emergent rate is as low as 0.1-1% of SN rate. Recent theoretical study of Galactic chemical evolution has demonstrated that it may be difficult for NSMs to emerge in metal-deficient region for these reasons [250].
The contribution from NSMs to r-process elements, however, depends on the timescale of metal enrichment and dynamics in chemical evolution models. Recent models include formation of metal-poor stars in small stellar systems like the current dwarf galaxies around the Milky Way, which takes longer timescale for chemical enrichment in general until accreting into the Galactic halo. Such models at least partially reconcile the difficulty of the timescale problem of NSM contribution [251]. Recent discoveries of r-process-enhanced stars in dwarf galaxies like Ret. II are providing useful constraints on such models (see below).
Observations of early generations of stars also made progress in revealing the nature of the r-process and identifying the sites. Chemical compositions of heavy neutron-capture elements show large star-to-star scatter, among which a small fraction of extremely metal-poor (EMP) stars have large excess in rprocess elements like Eu [252]. This indicates that these heavy elements are produced by quite rare events and that the gas clouds were not well mixed in the early universe [253]. Ultra-faint dwarf galaxy named Ret. II which has just thousands of solar masses was discovered, and interestingly most stars seem to be r-process-enhanced [254,255]. This suggests that the progenitor gas cloud was polluted by a few but very efficient r-process events, which could be NSMs, MHDJs or collapsars. It is suggested in a hierarchical structure formation scenario that the Milky Way halo was formed from merging dwarf galaxies which consists of r-process enhanced stars.
China-Japan collaboration team has recently reported the discovery of an r-process-enhanced star with relatively high metallicity compared to those found previously [256]. The metallicity of this object is about 30 times lower than the Sun, but its Eu abundance is as high as the solar value. Whereas such object had not been found previously in the Milky Way, similar objects are known in satellite dwarf galaxies [257]. This discovery suggests that the hierarchical structure formation scenario of the Milky Way is applicable to relatively metal-rich halo stars, as mentioned above, such that the objects were formed in dwarf galaxies first and then accreted later into the Milky Way. Interestingly, objects having similarly large excess of r-process elements have been found among moderately metal-poor stars (−2 <[Fe/H]< −1) in the past two years.
The detailed abundance patterns determined for individual objects provide strong constraint on understanding of r-process nature. A remarkable result found for r-process-enhanced EMP stars is that their abundance patterns of heavy elements are very similar to each other, and also similar to that of the solar-system r-process component, which is called "universality" of the r-process [253]. This empirical fact invited many theoretical models of SN r-process, e.g. [243,258,259], as shown in Fig. 6. Recent calculations have suggested that isotopic abundance patterns as a function of nuclear mass number A are quite different although elemental distributions as a function of atomic number Z are similar to one another among CCSNe, collapsars and NSMs [235,250]. Isotopic separation of r-process elements in EMP halo stars are highly desirable. On the other hand, diversity of abundance patterns of light neutroncapture elements also are found in metal-poor stars [260]. Further theoretical studies of explosive nucleosynthesis and observations of SN and GW events are highly desirable to uncover the origins of neutron-capture elements and the mechanisms to produce diversity in their abundance ratios.

Nuclear electroweak response and ν-astrophysics
Nuclear electroweak processes affect various masses of stellar evolution from the main sequence stage until the white dwarfs or gravitational collapse of massive stars [261,262,263]. Neutrinos play an essential role in successful explosion of CCSNe by heating the gas behind the shock [264] and affect various explosive nucleosyntheses. Since solar neutrinos are the direct observable to exhibit nuclear fusion inside invisible solar interior, their observed flux serves to construct the standard solar model [265]. Nuclear electroweak response has been studied in laboratory experiments by measuring nuclear β-decays and hadronic charge-exchange (CEX) and photo-nuclear reaction cross sections because neutrinos are thus invaluable messengers of weak processes in celestial events. The CEX and photo-induced reactions are analogous to charged and neutral current neutrino-induced reactions, respectively.
Systematic theoretical calculations of the neutrino-nucleus reaction cross sections have been performed based on sophisticated theoretical models such as nuclear shell models using new generations of Hamiltonian based on the progress in physics of exotic nuclei [266,267] and quasi-particle random phase approximation (QRPA) [268,269]. The results very well explain the experimental data on 12 C [270] and 56 Fe [271,272], and these models have been extensively applied to provide unmeasured ν-nucleus reaction cross sections of astrophysical interest.
Neutrino experiments in Asian countries/regions have contributed remarkably to determine the mixing angles, i.e., θ 12 (Super-KAMIOKANDE, Kam-LAND), θ 23 (Super-KAMIOKANDE, K2K, T2K) and θ 13 (T2K, Daya Bay, RENO) although mass hierarchy and CP-violation phase are still unknown. Upon these experimental efforts, it is the recent focus to determine the unknown oscillation parameters in terms of SN expolosions that provide threeflavor neutrinos and their anti-particles [273].
The neutrinos emitted from proto-neutron stars propagate through the stellar interior where several nuclei are produced by ν-nucleus interactions. The nuclides heavier than iron are produced in the r-and νp-processes in neutrinodriven winds [274,275,276,277] and the ν-process in outer layers [278,279,280,281]. Neutrinos also play an important role in the NSM outflows where a wide mass rage of r-process nuclei are produced [243,282,283] The 7 Li and 11 B are particularly important [281]. A recent improved and more consistent calculation of stellar evolution from He star to SN explosion has found a dependence of their yields on stellar metallicity [284].
It was found for the first time that remarkable amount of radioactive isotope 98 Tc is produced via charged-current reactions byν e although the other nuclei are produced mainly by the ν e -nucleus interactions. Thus, 98 Tc is characterized as a special probe of SNν e spectrum [285]. Extensive calculations showed that another short-lived radioactive isotope 92 Nb and stable ones 138 La and 180 Ta also are produced in the ν-process. These nuclei are useful to constrain still unknown mass hierarchy [286,287,288]. Fig. 7 shows an example of the 7 Be and 92 Nb abundances as a function of the Lagrangian mass M r . The high-density MSW effect [289,290] occurs near the bottom of the He-layer. On the other hand, the collective flavor oscillation induced by coherent ν-ν scattering [291,292] occurs inside the iron core. Therefore, 7 Be which is mainly produced in the He-layer is affected by both MSW and collective oscillation effects, while 92 Nb which is predominantly produced in the O/Ne/Mg-layer is affected by only the collective oscillation effect. Neutrino mass hierarchy is thus imprinted in each isotopic abundance in a different specific manner. Almost all these nuclei are ejected in SN explosions, from which the siliconcarbide pre-solar grains called SN grains condensate. Meteoritic measurement of SN grains is highly desirable to find a piece of evidence for the neutrino mass hierarchy [288,293].
Neutrino also holds the keys in the early universe because the cosmic expansion rate depends on the generations of light neutrino family [294]. Several short-baseline neutrino oscillation experiments reported an anomaly in the measured flux which could be explained by fission-induced reactions or assuming the non-active fourth-generation neutrino which is called "sterile neutrino" [295]. Reactor neutrino experiments at Daya Bay (China), RENO (South Korea) and JSNS 2 (Japan) are going to look for the signal of sterile neutrinos in sub keV region.
The dark energy and dark matter are the other key ingredients of the universe. Dark matter candidates include exotic weakly interacting massive particles (WIMPs) such as SUSY particles, right-handed massive neutrinos and axions. These hypothetical particles could be directly detected via the elastic scattering off the nuclei or atomic electrons through the weak interaction. Several experiments for WIMP hunting are underway at CDEX (China), PandaX (China), COSINE (South Korea), XMASS (Japan), NEWAGE (Japan), and PICO-LON (Japan).

Future Facilities
New facilities such as large telescope, neutrino observatory, next generation of Radioactive Ion Beam facilities and underground nuclear astrophysics experiments are being built or planned. More exciting discoveries are awaiting for us in coming decades. Here we briefly report some progresses of the facilities and related scientific problems.  [284]. Colors distinguish γ-process region at high temperature, O/Ne/Mg, O/C, C/He, and He-rich layers for ν-process region as labeled. 7 Be is mainly produced in the outer He-layer, while 92 Nb is produced also in the central hot region via the γ-process as well as ν-process.

Multi-messenger Astronomy
The first successful detection of the neutrinos from SN1987A with KAMIOKANDE-III, IMB and Baksan experiments opened the door to neutrino astronomy, which was followed by successive discovery of neutrino oscillation. This event of the century enriched our understanding of explosive nucleosynthesis and the roles of neutrino interactions in SNe. Supernovae provide heat and matter into space and drive the chemical and dynamical evolution of galaxies in the cycle of star birth -stellar evolution -explosion -next generation of star birth [296]. One of the current targets in neutrino astronomy at the underground laboratories is to detect relic SN neutrinos [297], which hints to determine the SN event rate in the Milky Way. In Japan, Super-KAMIOKANDE detector is in and Xγ-rays are respectively the messengers of gravitational, weak, electromagnetic, and electromagnetic plus strong nuclear forces that make a variety of physical and chemical processes in the universe, galaxies and stars.
operation with improved sensitivity by adding gadolinium into the Cherenkov medium made of light water. Next project Hyper-KAMIOKANDE with remarkably upgraded size of water Cherenkov detector has started. In China, Jiangmen Underground Neutrino Observatory (JUNO) is under construction. GW astronomy will link to neutrino astronomy. GW170817 made a big impact on the studies of r-process nucleosynthesis and EoS of neutron stars. Simultaneous detection of both GW and neutrinos in the future would provide more details how the elementary particle and nuclear processes go on in catastrophic SN explosions as well as merging phenomena of relativistic compact objects. KAGRA in Japan has started operation 1 , LIGO-India currently under construction [41], and the space-based GW detection projects TianQin [298] and Taiji [299] are planned in China.
Optical-to-infrared spectroscopic observations of SN1987A and GW170817 also gave fascinating insights into nucleosynthesis of heavy nuclei and dynamics ejecta. SN1987A appeared very close to us in Large Magellanic Cloud at a distance of 140kpc, and both multi-wave length photons and neutrinos were successfully detected. Although GW170817 appeared at 40Mpc which is relatively close and bright enough for GW detection and spectroscopic observations, it was too far to observe neutrinos. Future follow-up spectroscopic observations require deeper observations. High-resolution spectroscopic observation of faint stars in the Milky Way is planned with Subaru Telescope linking to large op-tical survey with LAMOST. Construction of next generation extremely large ground-based telescopes like Thirty Meter Telescope (TMT), Giant Magellan Telescope (GMT) and the space telescope like JWST are planned as international collaboration projects including Japan, Korea, China, India and Australia.
The last piece of multi-messengers is the messenger of strong nuclear force and electromagnetism in the atomic and nuclear processes in the cosmos. Space X-Ray satellite mission XRISM (X-Ray Imaging and Spectroscopy Mission) is planned at ISAS in Japan as post HITOMI project after the successful X-ray satellite missions of HAKUCHO, TENMA, GINGA, ASKA, and SUZAKU in the past. The enhanced X-ray Timing and Polarimetry mission (eXTP) is planned to be launched in 2027 to study the state of matter under extreme conditions of density, gravity and magnetism by the collaboration between Chinese and European institutions. Technology in analyzing chemical compositions of meteorites or returned sample from asteroids like Ryugu in Hayabusa-II project 2 would bring new insight how our solar-system formed and what the primordial solar material was made of. These might bring a piece of evidence for solving the mystery why optical chirality of amino acids on Earth is broken to be only left-handed [300], which hints the creation of life. However, as for the origin of elements in the cosmos, these are still limited to nearby solar system. Complementary spectroscopic observations of various ages of stars await to seek for the cosmic and galactic evolution of atomic nuclides from the big-bang era until the solar system formation.
It relies highly on the theoretical progress, too, in modeling the big-bang universe and high-energy phenomena like SNe and NSMs including numerical simulations of explosive nucleosynthesis there, as illustrated in Fig. 8. Newgeneration supercomputers such as Japanese Fugaku and Chinese Tianhe and Sunway Taihu-Light have enormously big memory and extremely high-speed computing power so that they can simulate cosmic, galactic and stellar evolution and associated nucleosynthetic phenomena by integrating all observed information and knowledge to be brought into multi-messenger astronomy and astrophysics.

Future Facilities in Experimental Nuclear Physics
We have several world-class facilities for radioactive ion beams existing in east and southeast Asia to study important nuclear astrophysics problems. Moreover, there are several upgrades and constructions toward more powerful facilities. Among them, upgrades of RIKEN RI Beam Factory (RIBF) and J-PARC in Japan, Rare Isotope Accelerator complex for ON-line experiments (RAON) in Korea, and High-Intensity heavy ion Accelerator Facility (HIAF) in China are introduced in this section.
RIBF facility has been in operation for more than ten years with its world's best capability of producing rare isotopes far from the stability. A variety of experimental studies related to nuclear astrophysics have been performed with RI beams of them (or these isotopes), as briefly described in previous sections. The upgrade of RIBF is being planned to continuously provide the best beam quality to our nuclear physics community in the world, together with upcoming facilities, such as FRIB in USA, FAIR in Germany, RAON and HIAF. A new charge-stripping scheme with specially-designed ring devices will transmit beam ions with multiple charge states, and recover the decrease of the beam in the two existing strippers. As a result, one order of magnitude higher beam intensity is expected [301]. Installation of a new superconducting LINAC, other improvements of the accelerators and the fragment separator BigRIPS will increase the accessibility of more neutron/proton-rich nuclei including isotopes related to the r-process.
J-PARC is a proton accelerator facility, and the intensive high-energy proton beams are used to produce high-flux secondary beams of neutrons, neutrinos and mesons such as kaons and pions. Astrophysics oriented investigations such as (n,γ) type reactions involved in the s-process are studied as well as properties of hyper-nuclei, hyperon-nucleon scattering. Neutrinos produced in the pion beam are also used to study neutrino oscillation with the KAMIOKA neutrino detector. The current extension plan for the hadron beamline will lead to systematic studies of strangeness nuclear physics and hence shed light on the hyperon-mixed matter.
One of the main scientific objectives of the RAON accelerator facility is to perform nuclear astrophysics experiments using RI beams. While RAON is designed to produce RI beams using known techniques including spectrometer, Isotope Separation On-Line (ISOL) and In-flight Fragmentation (IF), the most remarkable configuration is to utilize the combination of ISOL and IF system to produce more exotic species of isotopes. The method is very unique to provide neutron-rich beams with a greater intensity than that achieved by any existing facilities. Expected RI beams and their intensities are shown in Fig.9. Details of the beam production methods are described in Ref. [55].
The KoBRA at RAON is a multi-purpose mass separator for studying nuclear structures and reactions of exotic nuclei in the energy range of less than 30 MeV/u, as well as nuclear astrophysics experiments. The KoBRA project is divided into two stages. The first stage of KoBRA, will be utilized to produce RI beams using stable beams for the early phase of RAON. KoBRA will be extended by having more beam line components and a spectrometer, which is the second stage, for the study of astrophysically interesting reactions including 15 O(α,γ) 19 Ne and 14 O(α,p) 17 F reactions.
HIAF is a new accelerator facility in China, aiming to provide high intensity heavy ion beams from a LINAC with a maximum current up to 1 emA. With a booster, the beam energies can reach 9.3 GeV/u for proton and 0.8 GeV/u for 238 U. The lower energy terminal after the LINAC will be mainly used to produce new isotopes/elements using multiple nucleon transfer reactions or fusion reactions for the studies of nuclear reaction and spectroscopy. The research program with the higher energy beam from the booster will focus on the production of neutron-rich isotopes, precise mass measurements of exotic nuclei, charge exchange reactions, properties of hypernucleus, EoS of nuclear matter and the QCD phase diagram. The phase-I construction will be completed by 2025. China Initiated Accelerator Driven Subcritical System (CIADS) is also being constructed adjacent to HIAF to transmute nuclear waste using a powerful proton beam from the LINAC. In phase-II, an ISOL terminal will be built at CIADS to feed the neutron rich isotopes into HIAF. More critical r-process nuclei will be produced by the projectile fragmentation reaction or multiple nucleon transfer reaction induced by the neutron-rich beams. Furthermore, successful collaborations among different facilities are crucial since the teamwork will lead the project development in a very efficient way and improve the performance of the experimental instruments. New experimental techniques such as silicon detector arrays, γ-ray detector arrays and active target time projection chambers are critical for the success of the facilities. Collaborative efforts among the facilities will not only be very helpful to solve many challenging problems during the developments of novel technologies in each institute, but also offer new possibilities to operate expensive arrays by sharing the cost.

JUNA: Jinping Underground Nuclear Astrophysics experiment
The determination of the very low nuclear reaction rates in stars is a long standing challenge in nuclear astrophysics. In order to measure the extremely small nuclear cross section in the Gamow window, one of the most effective methods is to take the advantage of the rock with a thickness of several kilometers to shield the cosmic rays and the corresponding backgrounds. In 2014, the second phase of the China JinPing underground Laboratory (CJPL-II) was initialized. The CJPL-II is 2400 meters deep and listed as the deepest underground laboratory around the world. Jinping Underground lab for Nuclear Astrophysics (JUNA) collaboration is taking the advantage of the ultra-low background of CJPL-II and the high intensity accelerator to measure the critical reactions , 12 C(α,γ) 16 O, 13 C(α,n) 16 O, 19 F(p,α) 16 O and 25 Mg(p,γ) 26 Al at their stellar energies [69]. The layout of JUNA and other CJPL-II projects can be found in Fig. 10. JUNA is managed by CIAE, jointly supported by NSFC, CNNC and CAS. Development of high intensity accelerator have been completed. It is compatible to deliver proton and 4 He + beams with an intensity up to 10 emA and 4 He 2+ beam with an intensity up to 2 emA. A number of test experiments have been carried out at the ground level. The first beam was delivered in the Jinping Underground lab in Dec. 2020 and the experimental campaign began in Jan. 2021 to provide more accurate reaction rates for the studies such as s-process, helium burning in stars and the galactic radioactivities. JUNA collaboration is also planning to add a heavy ion accelerator with higher energies to study the critical reactions, such as 12 C+ 12 C, in the advanced burning phases of stars.

Conclusion
By detecting photons with the electromagnetic interaction, neutrinos with the weak interaction, cosmic rays with the strong interaction, GW with the gravitational interaction, multi-messenger astronomy has greatly enriched our ways of perceiving the universe. The successful multi-messenger observation of GW170817 has advanced greatly our knowledge of r-process and the EoS of nuclear matter. The collaborative efforts of LAMOST and Subaru discovered an r-process-enhanced star with relatively high metallicity, suggesting a hierarchical structure formation scenario of the Milky Way. The progresses made at nuclear facilities, such as discoveries of new isotopes and precise measurements of nuclear mass and lifetime, precise S(E2) factor of 12 C(α,γ) 16 O and measurements of the difficult reactions involving 7 Be, have eliminated some critical nuclear uncertainties, and provide a more solid foundation for the development of stellar models. The JUNA collaboration is likely to achieve some breakthroughs in the direct measurements of the critical cross sections at stellar energies with the intense beams in Jinping deep underground lab. New method of producing even more neutron-rich isotopes is being pursued and likely to open a new door towards the r-process studies. Future observatories, such as TMT and JUNO, will bring the human eyesight even further in the Universe and deeper inside of stars. As the new nuclear facilities such as RAON and HIAF and better technologies become available, we will continue to challenge the two long standing problems, the determination of the very low nuclear reaction rates in stars and the determination of properties and reactions of very neutron deficient or very neutron rich unstable nuclei. Complimentary to the other reaction library projects in US and Europe, theoretical collaborations such as A3LIB are developing global databases for the astrophysical applications with predictable errors. With the accurate nuclear reaction rates, more powerful computing resources, and advances in theory aspects, stellar models will become more realistic and reliable and eventually reveal the truths of stars carried by the multi-messengers.
Interdisciplinary feature of nuclear astrophysics demands the close collaborations among astronomers, astrophysicists and nuclear physicists and among the facilities. As we demonstrated in the paper, NO single facility or model will answer all the quests in our field. How to be successful in nuclear astrophysics? Here are the advises from Willy Fowler: seek for truth, work hard and help people.

Authors' contributions
All authors contributed equally to all aspects of the manuscript. All authors read and approved the final manuscript.