Search

A new neutron SIMulation program based on the versatile GEANT4 toolkit, neuSIM4, has been developed to describe interactions of neutrons in the NE213 liquid scintillator from 0.1 to 3000 MeV. neuSIM4 is designed to accommodate complicated modern detector geometry setups with multiple scintillator detectors, each of which can be outfitted with more than one photo-multiplier. To address a broad spectrum of neutron energies, two new neutron interaction physics models, KSCIN and NxQMD, have been implemented in GEANT4. For neutrons with energy below 110 MeV, we incorporate a total of eleven neutron induced reaction channels on hydrogen and carbon nuclei, including nine carbon inelastic reaction channels, into KSCIN. Beyond 110 MeV, we implement a neutron induced reaction model, NxQMD, in GEANT4. We use its results as reference to evaluate other neutron-interaction physics models in GEANT4. We find that results from an existing cascade physics model (INCL) in GEANT4 agree very well with the results from NxQMD, and results from both codes agree with new and existing light response data. To connect KSCIN to NxQMD or INCL, we introduce a transition region where the contribution of neuSIM4 linearly decreases with corresponding increased contributions from NxQMD or INCL. To demonstrate the application of the new code, we simulate the light response and performance of a 2 x 2 m2 neutron detector wall array consisting of 25 2m-long scintillation bars. We are able to compare the predicted light response functions to the shape of the experimental response functions and calculate the efficiency of the neutron detector array for neutron energies up to 200 MeV. These simulation results will be pivotal for understanding the performance of modern neutron arrays with intricate geometries, especially in the measurements of neutron energy spectra in heavy-ion reactions.

Properties of the nuclear equation of state (EoS) can be probed by measuring the dynamical properties of nucleus-nucleus collisions. In this study, we present the directed flow ($v_1$), elliptic flow ($v_2$) and stopping (VarXZ) measured in fixed target Sn + Sn collisions at 270 AMeV with the S$\pi$RIT Time Projection Chamber. We perform Bayesian analyses in which EoS parameters are varied simultaneously within the Improved Quantum Molecular Dynamics-Skyrme (ImQMD-Sky) transport code to obtain a multivariate correlated constraint. The varied parameters include symmetry energy, $S_0$, and slope of the symmetry energy, $L$, at saturation density, isoscalar effective mass, $m_{s}^*/m_{N}$, isovector effective mass, $m_{v}^{*}/m_{N}$ and the in-medium cross-section enhancement factor $\eta$. We find that the flow and VarXZ observables are sensitive to the splitting of proton and neutron effective masses and the in-medium cross-section. Comparisons of ImQMD-Sky predictions to the S$\pi$RIT data suggest a narrow range of preferred values for $m_{s}^*/m_{N}$, $m_{v}^{*}/m_{N}$ and $\eta$.

The boundaries of the Chart of Nuclides contain exotic isotopes that possess extreme proton-toneutron asymmetries. Here we report on strong evidence of 9N, one of the most exotic proton-rich isotopes where more than one half of its constitute nucleons are unbound. With seven protons and two neutrons, this extremely proton-rich system would represent the first-known example of a ground-state five-proton emitter. The invariant-mass spectrum of its decay products can be fit with two peaks whose energies are consistent with the theoretical predictions of an open-quantum-system approach, however we cannot rule out the possibility that only a single resonance-like peak is present in the spectrum., Comment: 4 figures, 7pages

Understanding the neutron background is essential for determining the neutron yield from nuclear reactions. In the analysis presented here, the shadow bars are placed in front of neutron detectors to determine the energy dependent neutron background fractions. The measurement of neutron spectra with and without shadow bars is important to determine the neutron background more accurately. The neutron background, along with its sources and systematic uncertainties, are explored with a focus on the impact of background models and their dependence on neutron energy., Comment: 7 pages, 10 figures

Experimental information on fragment emissions is important in understanding the dynamics of nuclear collisions and in the development of transport model simulating heavy-ion collisions. The composition of complex fragments emitted in the heavy-ion collisions can be explained by statistical models, which assume that thermal equilibrium is achieved at collision energies below 100 MeV/u. Our new experimental data together with theoretical analyses for light particles from Sn+Sn collisions at 270 MeV/u, suggest that the hypothesis of thermal equilibrium breaks down for particles emitted with high transfer momentum. To inspect the system's properties in such limit, the scaling features of the yield ratios of particles from two systems, a neutron-rich system of ${}^{132}\mathrm{Sn}+{}^{124}\mathrm{Sn}$ and a nearly symmetric system of ${}^{108}\mathrm{Sn}+{}^{112}\mathrm{Sn}$, are examined in the framework of the statistical multifragmentation model and the antisymmetrized molecular dynamics model. The isoscaling from low energy particles agree with both models. However the observed breakdown of isoscaling for particles with high transverse momentum cannot be explained by the antisymmetrized molecular dynamics model.

In the framework of the Statistical Multifragmentation Model, the nuclear isoscaling analysis is extended to constrain the ratio between the sizes of the decaying sources formed in a collision between two heavy ions. It is found that the ratio between the probabilities of observing n fragments in each event, for each of the sources, follows a scaling law, similar to the traditional nuclear isoscaling. However, the corresponding slope is also sensitive to the sources' sizes. This property is explained analytically using the grand-canonical ensemble. The extent to which our findings are affected by finite size effects and by the deexcitation of the hot primary fragments is also investigated. The scaling turns out to be robust and weakly affected by effects implied by these two aspects. We also find that the Poisson distribution is a fairly good approximation to the above mentioned probabilities, associated with both the primordial fragments, produced at the breakup stage, and the final ones, found at the end of the fragment deexcitation process., Comment: 6 pages, 3 figures

Nuclear Astrophysics is a field at the intersection of nuclear physics and astrophysics, which seeks to understand the nuclear engines of astronomical objects and the origin of the chemical elements. This white paper summarizes progress and status of the field, the new open questions that have emerged, and the tremendous scientific opportunities that have opened up with major advances in capabilities across an ever growing number of disciplines and subfields that need to be integrated. We take a holistic view of the field discussing the unique challenges and opportunities in nuclear astrophysics in regards to science, diversity, education, and the interdisciplinarity and breadth of the field. Clearly nuclear astrophysics is a dynamic field with a bright future that is entering a new era of discovery opportunities., Comment: 96 pages. Submitted to Journal of Physics G

A challenge preventing successful inverse kinematics measurements with heavy nuclei that are not fully stripped is identifying and tagging the beam particles. For this purpose, the HEavy ISotope Tagger (HEIST) has been developed. HEIST utilizes two micro-channel plate timing detectors to measure time of flight, a multi-sampling ion chamber to measure energy loss, and a high purity Ge detector to identify isomer decays and calibrate the isotope identification system. HEIST has successfully identified $^{198}$Pb and other nearby nuclei at energies of about 75 MeV/A. In the experiment discussed, a typical cut containing 89\% of all $^{198}$Pb$^{+80}$ in the beam had a purity of 86\%. We examine the issues of charge state contamination. The observed charge state populations of these ions are presented and are moderately well described by the charge state model GLOBAL., Comment: The following article has been submitted to Review of Scientific Instruments

To efficiently detect energetic light charged particles, it is common to use arrays of energy-loss telescopes involving two or more layers of detection media. As the energy of the particles increases, thicker layers are usually needed. However, carrying out measurements with thick-telescopes may require corrections for the losses due to nuclear reactions induced by the incident particles on nuclei within the detector and for the scattering of incident particles out of the detector, without depositing their full energy in the active material. In this paper, we develop a method for measuring such corrections and determine the reaction and out-scattering losses for data measured with the silicon-CsI(Tl) telescopes of the newly developed HiRA10 array. The extracted efficiencies are in good agreement with model predictions using the GEANT4 reaction loss algorithm for Z=1 and Z=2 isotopes. After correcting for the HiRA10 geometry, a general function that describes the efficiencies from the reaction loss in CsI(Tl) crystals as a function of range is obtained., Comment: 9 pages, 9 figures

Machine Learning (ML) algorithms have been demonstrated to be capable of predicting impact parameter in heavy-ion collisions from transport model simulation events with perfect detector response. We extend the scope of ML application to experimental data by incorporating realistic detector response of the S$\pi$RIT Time Projection Chamber into the heavy-ion simulation events generated from the UrQMD model to resemble experimental data. At 3 fm, the predicted impact parameter is 2.8 fm if simulation events with perfect detector is used for training and testing; 2.4 fm if detector response is included in the training and testing, and 5.8 fm if ML algorithms trained with perfect detector is applied to testing data that has included detector response. The last result is not acceptable illustrating the importance of including the detector response in developing the ML training algorithm. We also test the model dependence by applying the algorithms trained on UrQMD model to simulated events from four different transport models as well as using different input parameters on UrQMD model. Using data from Sn+Sn collisions at E/A=270 MeV, the ML determined impact parameters agree well with the experimentally determined impact parameter using multiplicities, except in the very central and very peripheral regions. ML selects central collision events better and allows impact parameters determination beyond the sharp cutoff limit imposed by experimental methods.

Excited states in 13O were investigated using inelastic scattering of an E/A=69.5-MeV 13O beam off of a 9Be target. The excited states were identified in the invariant-mass spectra of the decay products. Both single proton and sequential two-proton decays of the excited states were examined. For a number of the excited states, the protons were emitted with strong anisotropy where emissions transverse to the beam axis are favored. The measured proton-decay angular distributions were compared to predictions from distorted-wave born-approximation (DWBA) calculations of the spin alignment which was shown to be largely independent of the excitation mechanism. The deduced $^{13}$O level scheme is compared to ab initio no-core shell model with continuum (NCSMC) predictions. The lowest-energy excited states decay isotropically consistent with predictions of strong proton 1s1/2 structure. Above these states in the level scheme, we observed a number of higher-spin states not predicted within the model. Possibly these are associated with rotational bands built on deformed cluster configurations predicted by antisymmetrized molecular dynamics (AMD) calculations. The spin alignment mechanism is shown to be useful for making spin assignments and may have widespread use.

The large imbalance in the neutron and proton densities in very neutron rich systems increases the nuclear symmetry energy so that it governs many aspects of neutron stars and their mergers. Extracting the density dependence of the symmetry energy therefore constitutes an important scientific objective. Many analyses have been limited to extracting values for the symmetry energy, $S_0$, and its ``derivative'', $L$, at saturation density $\rho_0 \approx 2.6 \times 10^{14}~\mathrm{g/cm^3}$ $\approx 0.16~\mathrm{nucleons/fm^{3}}$, resulting in constraints that appear contradictory. We show that most experimental observables actually probe the symmetry energy at densities far from $\rho_0$, making the extracted values of $S_0$ or $L$ imprecise. By focusing on the densities these observables actually probe, we obtain a detailed picture of the density dependence of the symmetry energy from $0.25\rho_0$ to $1.5\rho_0$. From this experimentally derived density functional, we extract $L_{01}=53.1\pm6.1 MeV$ at $\rho \approx 0.10~\mathrm{fm^{-3}}$, a neutron skin thickness for $^{208}Pb$ of $R_{np} =$ $0.23\pm0.04$ fm, a symmetry pressure at saturation density of $P_0=3.2\pm1.2 MeV/fm^3$ and suggests a radius for a 1.4 solar mass neutron star of $13.1\pm0.6$ km., Comment: 6 pages, 4 figures

Background: $^{132}$Sn+$^{124}$Sn collisions at the beam energy of 270 MeV$/$nucleon have been performed at the Radioactive Isotope Beam Factory (RIBF) in RIKEN to investigate the nuclear equation of state. Reconstructing impact parameter is one of the important tasks in the experiment as it relates to many observables. Purpose: In this work, we employ three commonly used algorithms in machine learning, the artificial neural network (ANN), the convolutional neural network (CNN) and the light gradient boosting machine (LightGBM), to determine impact parameter by analyzing either the charged particles spectra or several features simulated with events from the ultra-relativistic quantum molecular dynamics (UrQMD) model. Method: To closely imitate experimental data and investigate the generalizability of the trained machine learning algorithms, incompressibility of nuclear equation of state and the in-medium nucleon-nucleon cross sections are varied in the UrQMD model to generate the training data. Results: The mean absolute error $\Delta b$ between the true and the predicted impact parameter is smaller than 0.45 fm if training and testing sets are sampled from the UrQMD model with the same parameter set. However, if training and testing sets are sampled with different parameter sets, $\Delta b$ would increase to 0.8 fm. Conclusion: The generalizability of the trained machine learning algorithms suggests that these machine learning algorithms can be used reliably to reconstruct impact parameter in experiment., Comment: 9 pages,9 figures

Many neutron star (NS) properties, such as the proton fraction within a NS, reflect the symmetry energy contributions to the Equation of State that dominate when neutron and proton densities differ strongly. To constrain these contributions at supra-saturation densities, we measure the spectra of charged pions produced by colliding rare isotope tin (Sn) beams with isotopically enriched Sn targets. Using ratios of the charged pion spectra measured at high transverse momenta, we deduce the slope of the symmetry energy to be $42 < L < 117$ MeV. This value is slightly lower but consistent with the $L$ values deduced from a recent measurement of the neutron skin thickness of $^{208}$Pb.

Using the dcQMD transport model, the isoscalar and isovector in-medium potentials of the $\Delta$(1232) baryon are studied and information regarding their effective strength is obtained from a comparison to experimental pion production data in heavy-ion collisions below 800 MeV/nucleon impact energy. The best description is achieved for an isoscalar potential moderately more attractive than the nucleon optical potential and a rather small isoscalar relative effective mass m$^*_\Delta \approx$ 0.45. For the isovector component only a constraint between the potential's strength at saturation and the isovector effective mass difference can be extracted, which depends on quantities such as the slope of the symmetry energy and the neutron-proton effective mass difference. These results are incompatible with the usual assumption, in transport models, that the $\Delta$(1232) and nucleon potentials are equal. The density dependence of symmetry energy can be studied using the high transverse momentum tail of pion multiplicity ratio spectra. Results are however correlated with the value of neutron-proton effective mass difference. This region of spectra is shown to be affected by uncertain model ingredients such as the pion potential or in-medium correction to inelastic scattering cross-sections at levels smaller than 10$\%$. Extraction of precise constraints for the density dependence of symmetry energy above saturation will require experimental data for pion production in heavy-ion collisions below 800 MeV/nucleon impact energy and experimental values for the high transverse momentum tail of pion multiplicity ratio spectra accurate to better than 5$\%$., Comment: 22 pages, 13 figures

In the past two decades, pions created in the high density regions of heavy ion collisions have been predicted to be sensitive at high densities to the symmetry energy term in the nuclear equation of state, a property that is key to our understanding of neutron stars. In a new experiment designed to study the symmetry energy, the multiplicities of negatively and positively charged pions have been measured with high accuracy for central $^{132}$Sn+$^{124}$Sn, $^{112}$Sn+$^{124}$Sn, and $^{108}$Sn+$^{112}$Sn collisions at $E/A=270~\mathrm{MeV}$ with the S$\pi$RIT Time Projection Chamber. While the uncertainties of individual pion multiplicities are measured to 4\%, those of the charged pion multiplicity ratios are measured to 2\%. We compare these data to predictions from seven major transport models. The calculations reproduce qualitatively the dependence of the multiplicities and their ratios on the total neutron to proton number in the colliding systems. However, the predictions of the transport models from different codes differ too much to allow extraction of reliable constraints on the symmetry energy from the data. This finding may explain previous contradictory conclusions on symmetry energy constraints obtained from pion data in Au+Au system. These new results call for better understanding of the differences among transport codes, and new observables that are more sensitive to the density dependence of the symmetry energy., Comment: 8 pages, 4 figures, 1 table (accepted for publication in PLB)

The first measurement of the $\vec{^3He}(\vec{\gamma},p)d$ process was performed at the High Intensity $\gamma$-ray Source (HI$\gamma$S) facility at Triangle Universities Nuclear Laboratory (TUNL) using a circularly polarized, monoenergetic $\gamma$-ray beam and a longitudinally polarized $^3$He target. The spin-dependent asymmetry and the contribution from the two-body photodisintegration to the $^3$He Gerasimov-Drell-Hearn integrand are extracted and compared with state-of-the-art three-nucleon system calculations at the incident photon energy of 29.0 MeV. The data are in general agreement with the various theoretical predictions based on the Siegert theorem or on explicit inclusion of meson-exchange currents., Comment: 5 pages, 3 figures. arXiv admin note: text overlap with arXiv:1304.5442

We use a Bayesian inference analysis to explore the sensitivity of Taylor expansion parameters of the nuclear equation of state (EOS) to the neutron star dimensionless tidal deformability ($\Lambda$) on 1 to 2 solar masses neutron stars. A global power law dependence between tidal deformability and compactness parameter (M/R) is verified over this mass region. To avoid superfluous correlations between the expansion parameters, we use a correlation-free EOS model based on a recently published meta-modeling approach. We find that assumptions in the prior distribution strongly influence the constraints on $\Lambda$. The $\Lambda$ constraints obtained from the neutron star merger event GW170817 prefer low values of $L_\text{sym}$ and $K_\text{sym}$, for a canonical neutron star with 1.4 solar mass. For neutron star with mass $<1.6$ solar mass, $L_\text{sym}$ and $K_\text{sym}$ are highly correlated with the tidal deformability. For more massive neutron stars, the tidal deformability is more strongly correlated with higher order Taylor expansion parameters.

The SAMURAI Pion Reconstruction and Ion-Tracker Time Projection Chamber (S$\pi$RIT TPC) was designed to enable measurements of heavy ion collisions with the SAMURAI spectrometer at the RIKEN Radioactive Isotope Beam Factory and provide constraints on the Equation of State of neutron-rich nuclear matter. The S$\pi$RIT TPC has a 50.5 cm drift length and an 86.4 cm $\times$ 134.4 cm pad plane with 12,096 pads that are equipped with the Generic Electronics for TPCs readout electronics. The S$\pi$RIT TPC allows excellent reconstruction of particles and provides isotopic resolution for pions and other light charged particles across a wide range of energy losses and momenta. Details of the S$\pi$RIT TPC are presented, along with discussion of the TPC performance based on cosmic ray and experimental data., Comment: 12 pages, 20 figures

Using the waveforms from a digital electronic system, an offline analysis technique on pulse shape discrimination (PSD) has been developed to improve the neutron-gamma separation in a bar-shaped NE-213 scintillator that couples to a photomultiplier tube (PMT) at each end. The new improved method, called the ``valued-assigned PSD'' (VPSD), assigns a normalized fitting residual to every waveform as the PSD value. This procedure then facilitates the incorporation of longitudinal position dependence of the scintillator, which further enhances the PSD capability of the detector system. In this paper, we use radiation emitted from an AmBe neutron source to demonstrate that the resulting neutron-gamma identification has been much improved when compared to the traditional technique that uses the geometric mean of light outputs from both PMTs. The new method has also been modified and applied to a recent experiment at the National Superconducting Cyclotron Laboratory (NSCL) that uses an analog electronic system., Comment: 8 pages, 8 figures

In this paper, we present a software framework, S$\pi$RITROOT, which is capable of track reconstruction and analysis of heavy-ion collision events recorded with the S$\pi$RIT time projection chamber. The track-fitting toolkit GENFIT and the vertex reconstruction toolkit RAVE are applied to a box-type detector system. A pattern recognition algorithm which performs helix track finding and handles overlapping pulses is described. The performance of the software is investigated using experimental data obtained at the Radioactive Isotope Beam Facility (RIBF) at RIKEN. This work focuses on data from $^{132}$Sn + $^{124}$Sn collision events with beam energy of 270 AMeV. Particle identification is established using $\left$ and magnetic rigidity, with pions, hydrogen isotopes, and helium isotopes., Comment: 10 pages, 17 figures

Time projection chambers (TPCs) are widely used in nuclear and particle physics. They are particularly useful when measuring reaction products from heavy ion collisions. Most nuclear experiments at low energy are performed in a fixed target configuration, in which the unreacted beam will pass through the detection volume. As the beam intensity increases, the buildup of positive ions created from the ionization of the detector gas by the beam creates the main source of space charge, distorting the nominal electric field of the TPC. This has a profound effect on the accuracy of the measured momenta of the emitted particles. In this paper we will discuss the magnitude of the effects and construct an observable more appropriate for fixed target experiments to study the effects. We also will present an algorithm for correcting the space charge and some of the implications it has on the momentum determination., Comment: 9 pages, 10 figures

Excited states in $^{14}$O have been investigated both experimentally and theoretically. Experimentally, these states were produced via neutron-knockout reactions with a fast $^{15}$O beam and the invariant-mass technique was employed to isolate the 1$p$ and 2$p$ decay channels and determine their branching ratios. The spectrum of excited states was also calculated with the Shell Model Embedded in the Continuum that treats bound and scattering states in a unified model. By comparing energies, widths and decay branching patterns, spin and parity assignments for all experimentally observed levels below 8 MeV are made. This includes the location of the second 2$^{+}$ state that we find is in near degeneracy with the third 0$^{+}$ state. An interesting case of sequential 2$p$ decay through a pair of degenerate $^{13}$N excited states with opposite parities was found where the interference between the two sequential decay pathways produces an unusual relative-angle distribution between the emitted protons., Comment: 19 pages, 20 figures

We use properties of doubly-magic nuclei, ab-initio calculations of low-density neutron matter, and of neutron stars to constrain the parameters of the Skyrme energy-density functional. We find all of these properties can be reproduced within a constrained family of Skyrme parameters. The maximum mass of a neutron star is found to be sensitive to the neutron effective mass. A value of [$ m^{*}_{\rm n}/m](\rho_0) = 0.60-0.65 $ is required to obtain a maximum neutron star mass of 2.1 solar masses. Using the constrained Skyrme functional with the aforementioned effective mass, the predicted radius for a neutron star of 1.4 solar masses is 12.4(1) km and $\Lambda$ = 423(40).

Particle-decaying states of the light nuclei $^{11,12}$N and $^{12}$O were studied using the invariant-mass method. The decay energies and intrinsic widths of a number of states were measured, and the momentum correlations of three-body decaying states were considered. A second 2$p$-decaying 2$^+$ state of $^{12}$O was observed for the first time, and a higher energy $^{12}$O state was observed in the 4$p$+2$\alpha$ decay channel. This 4$p$+2$\alpha$ channel also contains contributions from fission-like decay paths, including $^6$Be$_{g.s.}$+$^{6}$Be$_{g.s.}$. Analogs to these states in $^{12}$O were found in $^{12}$N in the 2$p$+$^{10}$B and 2$p$+$\alpha$+$^6$Li channels. The momentum correlations for the prompt 2$p$ decay of $^{12}$O$_{g.s.}$ were found to be nearly identical to those of $^{16}$Ne$_{g.s.}$, and the correlations for the new 2$^+$ state were found to be consistent with sequential decay through excited states in $^{11}$N. The momentum correlations for the 2$^+_1$ state in $^{12}$O provide a new value for the $^{11}$N ground-state energy. The states in $^{12}$N/$^{12}$O that belong to the $A$=12 isobaric sextet do not deviate from the quadratic isobaric multiplet mass equation (IMME) form., Comment: 16 pages, 16 figures, 7 tables, tables corrected

The LIGO-Virgo collaboration detection of the binary neutron-star merger event, GW170817, has expanded efforts to understand the Equation of State (EoS) of nuclear matter. These measurements provide new constraints on the overall pressure, but do not elucidate its origins, by not distinguishing the contribution to the pressure from symmetry energy which governs much of the internal structure of a neutron star. By combining the neutron star EoS extracted from the GW170817 event and the EoS of symmetric matter from nucleus-nucleus collision experiments, we extract the symmetry pressure, which is the difference in pressure between neutron and nuclear matter over the density region from 1.2$\rho_{0}$ to $4.5\rho_{0}$. While the uncertainties in the symmetry pressure are large, they can be reduced with new experimental and astrophysical results., Comment: arXiv admin note: substantial text overlap with arXiv:1901.07673

The binary neutron-star merger event, GW170817, has cast a new light on nuclear physics research. Using a neutron-star model that includes a crust equation of state (EoS), we calculate the properties of a 1.4 solar-mass neutron star. The model incorporates more than 200 Skyrme energy density functionals, which describe nuclear matter properties, in the outer liquid core region of the neutron star. We find a power-law relation between the neutron-star tidal deformability, $\Lambda$, and the neutron-star radius, R. Without an explicit crust EoS, the model predicts smaller R and the difference becomes significant for stars with large radii. To connect the neutron star properties with nuclear matter properties, we confront the predicted values for $\Lambda$, against the Taylor expansion coefficients of the Skyrme interactions. There is no pronounced correlation between Skyrme parameters in symmetric nuclear matter and neutron star properties. However, we find the strongest correlation between $\Lambda$ and $K_{sym}$, the curvature of the density dependence of the symmetry energy at saturation density. At twice the saturation density, our calculations show a strong correlation between $\Lambda$ and total pressure providing guidance to laboratory nucleus-nucleus collision experiments.

Efficiency corrected single ratios of neutron and proton spectra in central $^{112}$Sn+$^{112}$Sn and $^{124}$Sn+$^{124}$Sn collisions at 120 MeV/u are combined with double ratios to provide constraints on the density and momentum dependencies of the isovector mean-field potential. Bayesian analyses of these data reveal that the isoscalar and isovector nucleon effective masses, $m_s^* - m_v^*$ are strongly correlated. The linear correlation observed in $m_s^* - m_v^*$ yields a nearly independent constraint on the effective mass splitting $\Delta m_{np}^*= (m_n^*-m_p^*)/m_N = -0.05_{-0.09}^{+0.09}\delta$. The correlated constraint on the standard symmetry energy, $S_0$ and the slope, $L$ at saturation density yields the values of symmetry energy $S(\rho_s)=16.8_{-1.2}^{+1.2}$ MeV at a sensitive density of $\rho_s/\rho_0 = 0.43_{-0.05}^{+0.05}$.

The light output produced by light ions (Z<=4) in CsI(Tl) crystals is studied over a wide range of detected energies (E<=300 MeV). Energy-light calibration data sets are obtained with the 10 cm crystals in the recently upgraded High-Resolution Array (HiRA10). We use proton recoil data from 40,48Ca + CH2 at 28 MeV/u, 56.6 MeV/u, 39 MeV/u and 139.8 MeV/u and data from a dedicated experiment with direct low-energy beams. We also use the punch through points of p, d, and t particles from 40,48Ca + 58,64Ni, 112,124Sn collisions reactions at 139.8 MeV/u. Non-linearities, arising in particular from Tl doping and light collection efficiency in the CsI crystals, are found to significantly affect the light output and therefore the calibration of the detector response for light charged particles, especially the hydrogen isotopes. A new empirical parametrization of the hydrogen light output, L(E,Z=1,A), is proposed to account for the observed effects. Results are found to be consistent for all 48 CsI(Tl) crystals in a cluster of 12 HiRA10 telescopes.

The LIGO-Virgo collaboration ground-breaking detection of the binary neutron-star merger event, GW170817, has expanded efforts to understand the Equation of State (EoS) of nuclear matter. These measurements provide new constraints on the overall pressure, but do not, by itself, elucidate its microscopic origins, including the pressure arising from the symmetry energy, that governs much of the internal structure of a neutron star. To correlate microscopic constraints from nuclear measurements to the GW170817 constraints, we calculate neutron star properties with more than 200 Skyrme energy density functionals that describe properties of nuclei. Calculated neutron-star radii (R) and the tidal deformabilities which show a strong correlation with pressure at twice saturation density. By combining the neutron star EoS extracted from the GW170817 event and the EoS of symmetric matter from nucleus-nucleus collision experiments, we extract the density dependence of the symmetry pressure from 1.2 to 4.5 times saturation density. While the uncertainties in the symmetry pressure are large, they can be reduced with new experimental and astrophysical results., Comment: 10 pages, 3 figures. arXiv admin note: text overlap with arXiv:1807.06571

The structure of the extremely proton-rich nucleus $^{11}_{~8}$O$_3$, the mirror of the two-neutron halo nucleus $^{11}_{~3}$Li$_8$, has been studied experimentally for the first time. Following two-neutron knockout reactions with a $^{13}$O beam, the $^{11}$O decay products were detected after two-proton emission and used to construct an invariant-mass spectrum. A broad peak of width $\sim$3\,MeV was observed. Within the Gamow coupled-channel approach, it was concluded that this peak is a multiplet with contributions from the four-lowest $^{11}$O resonant states: $J^{\pi}$=3/2$^-_1$, 3/2$^-_2$, 5/2$^+_1$, and 5/2$^+_2$. The widths and configurations of these states show strong, non-monotonic dependencies on the depth of the $p$-$^9$C potential. This unusual behavior is due to the presence of a broad threshold resonant state in $^{10}$N, which is an analog of the virtual state in $^{10}$Li in the presence of the Coulomb potential. After optimizing the model to the data, only a moderate isospin asymmetry between ground states of $^{11}$O and $^{11}$Li was found., Comment: 6 pages, three figures

Neutron transfer reactions with fast secondary beams of $^{17}$Ne, $^{15}$O, and $^9$C have been studied with the HiRA and CAESAR arrays. Excited states of $^{18}$Ne, $^{16}$O, and $^{10}$C in the continuum have been identified using invariant-mass spectroscopy. The best experimental resolution of these states is achieved by selecting events where the decay fragments are emitted transverse to the beam direction. We have confirmed a number of spin assignments made in previous works for the negative-parity states of $^{18}$Ne. In addition we have found new higher-lying excited states in $^{16}$O and $^{18}$Ne, some of which fission into two ground-state $^8$Be fragments. Finally for $^{10}$C, a new excited state was observed. These transfer reactions were found to leave the remnant of the $^9$Be target nuclei at very high excitation energies and maybe associated with the pickup of a deeply-bound $^9$Be neutron., Comment: 21 pages, 18 figures

To relate constraints from nuclear physics to the tidal deformabilities of neutron stars, we construct a neutron star model that accepts input from a large collection of Skyrme density functions to calculate properties of 1.4 solar-mass neutron stars. We find that restricting this set of Skyrme to density functions that describe nuclear masses, isobaric analog states, and low energy nuclear reactions does not sufficiently restrict the predicted neutron-star radii and the tidal deformabilities. However, pressure constraints on the EoS around twice saturation density ($2\times2.74\times10^{14}g/cm^3$), obtained from high energy nucleus-nucleus collisions, does constrain predicted tidal deformabilities with uncertainties smaller than those obtained from the analysis of GW170817. We also found that the density-pressure constraint on the EoS obtained from a recent analysis of the neutron-star merger event agree very well with the density pressure constraints obtained from nuclear physics experiments published in 2002., Comment: arXiv admin note: substantial text overlap with arXiv:1807.06571

The LIGO-Virgo collaboration ground-breaking detection of the binary neutron-star merger event, GW170817, has intensified efforts towards the understanding of the equation of state (EoS) of nuclear matter. In this letter, we compare directly the density-pressure constraint on the EoS obtained from a recent analysis of the neutron-star merger event to density-pressure constraints obtained from nuclear physics experiments. To relate constraints from nuclear physics to the radii and the tidal deformabilities of neutron stars, we use a large collection of Skyrme density functionals that describe properties of nuclei to calculate properties of 1.4 solar maass neutron stars. We find that restricting this set of Skyrme equations of state to density functionals that describe nuclear masses, isobaric analog states, and low energy nuclear reactions does not sufficiently restrict the predicted neutron-star radii and the tidal deformabilities. Including pressure constraints on the EoS around twice saturation density, obtained from high energy nucleus-nucleus collisions, does constrain predicted radii and tidal deformabilities to be consistent with the results obtained from the analysis of GW170817. We discuss how new measurements of nucleus-nucleus collisions can improve these constraints on the EoS to be more restrictive than the current constraints from the GW170817 merger event., Comment: 11 pages, 3 figures

The density dependence of the nuclear symmetry energy governs important aspects of very neutron rich systems such as heavy nuclei and their collisions, neutron stars and their mergers. Many analyses of experimental data have generated constraints on the symmetry energy and its first derivative at saturation density, rho_0 ~ 2.7xE14 g/cm3. We show that each analysis does not accurately constrain the symmetry energy at rho_0, but rather at a lower density, rho_s, that is most sensitively probed by that analysis. Using published constraints on the symmetry energy and its first derivative at rho_0, we constrain the symmetry energy within a density range of 0.25rho_0 to 0.75rho_0 that is relevant to the inner crusts of neutron stars and to the neutrino-sphere of core-collapse supernovae. With appropriate data, the symmetry energy can be similarly constrained over a wider density range., Comment: 12 pages, 4 figures

In this work, two particular properties of the position-sensitive, thick silicon detectors (known as the "E" detectors) in the High Resolution Array (HiRA) are investigated: the thickness of the dead layer on the front of the detector, and the overall thickness of the detector itself. The dead layer thickness for each E detector in HiRA is extracted using a measurement of alpha particles emitted from a $^{212}$Pb pin source placed close to the detector surface. This procedure also allows for energy calibrations of the E detectors, which are otherwise inaccessible for alpha source calibration as each one is sandwiched between two other detectors. The E detector thickness is obtained from a combination of elastically scattered protons and an energy-loss calculation method. Results from these analyses agree with values provided by the manufacturer., Comment: Accepted for publication in Nuclear Instruments and Methods in Physics Research

The properties of the kinetic energy spectra of light isotopes produced in the breakup of a nuclear source and during the deexcitation of its products are examined. The initial stage, at which the hot fragments are created, is modeled by the Statistical Multifragmentation Model, whereas the Weisskopf-Ewing evaporation treatment is adopted to describe the subsequent fragment deexcita- tion, as they follow their classical trajectories dictated by the Coulomb repulsion among them. The energy spectra obtained are compared to available experimental data. The influence of the fusion cross-section entering into the evaporation treatment is investigated and its influence on the qual- itative aspects of the energy spectra turns out to be small. Although these aspects can be fairly well described by the model, the underlying physics associated with the quantitative discrepancies remains to be understood., Comment: 7 pages, 4 figures

We investigate correlations in the fragment momentum distribution due to the propagation of fragments under the influence of their mutual Coulomb field, after the breakup of an excited nuclear source.The magnitude of the effects on the nuclear temperatures obtained from such distributions is estimated with the help of a simple approach in which a charged fragment interacts with a homogeneous charged sphere. The resuslts are used to correct the temperatures obtained from the asymptotic momentum distributions of fragments produced by a Monte-Carlo simulation in which the system's configuration at breakup is provided by the canonical version of the Statistical Multifragmentation Model. In a separate calculation, the dynamics of this many-particle charged system is followed in a molecular dynamics calculation until the fragments are far away from the breakup volume. The results suggest that, although the magnitude of the corrections is similar in both models, many-particle correlations present in the second approach are non-negligible and should be taken into account in order to minimize ambiguities in such studies., Comment: 7 pages, 4 figures

Pion energy spectra are presented for central collisions of neutron-rich 132Sn+124Sn and neutron-deficient 108Sn+112Sn systems using simulations with Boltzmann-Uehling-Uhlenbeck transport model. These calculations, which incorporate isospin-dependent mean field potentials for relevant baryons and mesons, display a sensitivity to the pion spectra that could allow significant constraints on the density dependence of the symmetry energy and its mean field potential at supra-saturation densities. The predicted sensitivity increases with the isospin asymmetry of the total system and decreases with incident energy.

A simple but novel driver system has been developed to operate the wire gating grid of a Time Projection Chamber (TPC). This system connects the wires of the gating grid to its driver via low impedance transmission lines. When the gating grid is open, all wires have the same voltage allowing drift electrons, produced by the ionization of the detector gas molecules, to pass through to the anode wires. When the grid is closed, the wires have alternating higher and lower voltages causing the drift electrons to terminate at the more positive wires. Rapid opening of the gating grid with low pickup noise is achieved by quickly shorting the positive and negative wires to attain the average bias potential with N-type and P-type MOSFET switches. The circuit analysis and simulation software SPICE shows that the driver restores the gating grid voltage to 90% of the opening voltage in less than 0.20 $\mu$s. When tested in the experimental environment of a time projection chamber larger termination resistors were chosen so that the driver opens the gating grid in 0.35 $\mu$s. In each case, opening time is basically characterized by the RC constant given by the resistance of the switches and terminating resistors and the capacitance of the gating grid and its transmission line. By adding a second pair of N-type and P-type MOSFET switches, the gating grid is closed by restoring 99% of the original charges to the wires within 3 $\mu$s., Comment: 18 pages, 10 figures, 2 tables. Submitted to Nuclear Instruments & Methods in Physics Resarch A

Detector efficiency determination is essential to correct the measured yields and extract reliable cross sections of particles emitted in nuclear reactions. We investigate the efficiencies for measuring the full energies of light charged particle in arrays of CsI crystals employed in particle detection arrays such as HiRA, LASSA and MUST2. We perform these simulations with a GEANT4 Monte Carlo transport code implemented in the NPTool framework. Both Coulomb multiple scattering and nuclear reactions within the crystal can significantly reduce the efficiency of detecting the full energy of high energy particles. The calculated efficiencies decrease exponentially as a function of the range of the particle and are quite similar for both the hydrogen ($p, d, t$) and helium ($^3$He, $\alpha$) isotopes. The use of a close-packed array introduces significant position dependent efficiency losses at the interior boundaries between crystals that need to be considered in the design of an array and in the efficiency corrections of measured energy spectra.

KATANA - the Krakow Array for Triggering with Amplitude discrimiNAtion - has been built and used as a trigger and veto detector for the S$\pi$RIT TPC at RIKEN. Its construction allows operating in magnetic field and providing fast response for ionizing particles, giving the approximate forward multiplicity and charge information. Depending on this information, trigger and veto signals are generated. The article presents performance of the detector and details of its construction. A simple phenomenological parametrization of the number of emitted scintillation photons in plastic scintillator is proposed. The effect of the light output deterioration in the plastic scintillator due to the in-beam irradiation is discussed., Comment: 14 pages, 11 figures

Directed and elliptic flows of neutrons and light charged particles were measured for the reaction 197Au+197Au at 400 MeV/nucleon incident energy within the ASY-EOS experimental campaign at the GSI laboratory. The detection system consisted of the Large Area Neutron Detector LAND, combined with parts of the CHIMERA multidetector, of the ALADIN Time-of-flight Wall, and of the Washington-University Microball detector. The latter three arrays were used for the event characterization and reaction-plane reconstruction. In addition, an array of triple telescopes, KRATTA, was used for complementary measurements of the isotopic composition and flows of light charged particles. From the comparison of the elliptic flow ratio of neutrons with respect to charged particles with UrQMD predictions, a value \gamma = 0.72 \pm 0.19 is obtained for the power-law coefficient describing the density dependence of the potential part in the parametrization of the symmetry energy. It represents a new and more stringent constraint for the regime of supra-saturation density and confirms, with a considerably smaller uncertainty, the moderately soft to linear density dependence deduced from the earlier FOPI-LAND data. The densities probed are shown to reach beyond twice saturation., Comment: 25 pages, 22 figures, accepted for publication in Physical Review C

Using covariance analysis, we quantify the correlations between the interaction parameters in a transport model and the observables commonly used to extract information of the Equation of State of Asymmetric Nuclear Matter in experiments. By simulating $^{124}$Sn+$^{124}$Sn, $^{124}$Sn+$^{112}$Sn and $^{112}$Sn+$^{112}$Sn reactions at beam energies of 50 and 120 MeV per nucleon, we have identified that the nucleon effective mass splitting are most strongly correlated to the neutrons and protons yield ratios with high kinetic energy from central collisions especially at high incident energy. The best observable to determine the slope of the symmetry energy, L, at saturation density is the isospin diffusion observable even though the correlation is not very strong ($\sim$0.7). Similar magnitude of correlation but opposite in sign exists for isospin diffusion and nucleon isoscalar effective mass. At 120 MeV/u, the effective mass splitting and the isoscalar effective mass also have opposite correlation for the double n/p and isoscaling p/p yield ratios. By combining data and simulations at different beam energies, it should be possible to place constraints on the slope of symmetry energy (L) and effective mass splitting with reasonable uncertainties., Comment: 13 pages, 1 figure

Nuclear energy released by splitting Uranium and Thorium isotopes into two, three, four, up to eight fragments with nearly equal size are studied. We found that the energy released come from equally splitting the $^{235,238}$U and $^{230,232}$Th nuclei into to three fragments is largest. The statistical multifragmentation model is employed to calculate the probability of different breakup channels for the excited nuclei. Weighing the the probability distributions of fragments multiplicity at different excitation energies for the $^{238}$U nucleus, we found that an excitation energy between 1.2 and 2 MeV/u is optimal for the $^{235}$U, $^{238}$U, $^{230}$Th and $^{232}$Th nuclei to release nuclear energy of about 0.7-0.75 MeV/u.

The 9Be(28Mg,27Na) one-proton removal reaction with a large proton separation energy of Sp(28Mg)=16.79 MeV is studied at intermediate beam energy. Coincidences of the bound 27Na residues with protons and other light charged particles are measured. These data are analyzed to determine the percentage contributions to the proton removal cross section from the elastic and inelastic nucleon removal mechanisms. These deduced contributions are compared with the eikonal reaction model predictions and with the previously measured data for reactions involving the re- moval of more weakly-bound protons from lighter nuclei. The role of transitions of the proton between different bound single-particle configurations upon the elastic breakup cross section is also quantified in this well-bound case. The measured and calculated elastic breakup fractions are found to be in good agreement., Comment: Phys. Rev. C 2014 (accepted)