Showing 1–50 of 101 results for author: Kawano, T

Solving one-dimensional penetration problem for fission channel in the statistical Hauser-Feshbach theory

Authors: Toshihiko Kawano, Patrick Talou, Stephane Hilaire

Abstract: We solve the Schrödinger equation for an arbitrary one-dimensional potential energy to calculate the transmission coefficient in the fission channel of compound nucleus reactions. We incorporate the calculated transmission coefficients into the statistical Hauser-Feshbach model calculation for neutron-induced reactions on $^{235,238}$U and $^{239}$Pu. The one-dimensional model reproduces the evalu… ▽ More We solve the Schrödinger equation for an arbitrary one-dimensional potential energy to calculate the transmission coefficient in the fission channel of compound nucleus reactions. We incorporate the calculated transmission coefficients into the statistical Hauser-Feshbach model calculation for neutron-induced reactions on $^{235,238}$U and $^{239}$Pu. The one-dimensional model reproduces the evaluated fission cross section data reasonably well considering the limited number of model parameters involved. A resonance-like structure appears in the transmission coefficient for a double-humped fission barrier shape that includes an intermediate well, which is understood to be a quantum mechanical effect in the fission channel. The calculated fission cross sections for the neutron-induced reactions on $^{235,238}$U and $^{239}$Pu all exhibit a similar structure. △ Less

Submitted 1 July, 2023; originally announced July 2023.

Report number: LA-UR-23-25089

The Los Alamos evaluation of $^{239}$Pu neutron-induced reactions in the fast energy range

Authors: M. R. Mumpower, D. Neudecker, T. Kawano, M. Herman, N. Kleedtke, A. E. Lovell, I. Stetcu, P. Talou

Abstract: A major revision of the evaluation of $^{239}$Pu neutron-induced reaction cross sections is reported in the fast energy range. The evaluation starts at 2.5 keV incident neutron energy and has been extended up to 30 MeV. Several other notable changes are included in this evaluation since the release of ENDF/B-VIII.0 including the adoption of the Standards fission cross section, inclusion of new rad… ▽ More A major revision of the evaluation of $^{239}$Pu neutron-induced reaction cross sections is reported in the fast energy range. The evaluation starts at 2.5 keV incident neutron energy and has been extended up to 30 MeV. Several other notable changes are included in this evaluation since the release of ENDF/B-VIII.0 including the adoption of the Standards fission cross section, inclusion of new radiative capture data of Mosby et al., inclusion of the (n,2n) data of Meot et al., in addition to advances in the treatment of reaction modeling. In contrast to previous evaluation efforts, this evaluation is reproducible with detailed information stored chronologically utilizing a Git repository. The final evaluation results have been compiled into an ENDF-formatted file, which has been processed successfully through NJOY, checked for internal consistency, benchmarked versus older evaluations and validated against a suite of critical assemblies and pulsed-spheres. △ Less

Submitted 6 February, 2023; originally announced February 2023.

Comments: 31 pages, 31 figures, internal Los Alamos report

Report number: LA-UR-22-24121

Consideration of memory of spin and parity in the fissioning compound nucleus by applying the Hauser-Feshbach fission fragment decay model to photonuclear reactions

Authors: Toshihiko Kawano, Amy E. Lovell, Shin Okumura, Hirokazu Sasaki, Ionel Stetcu, Patrick Talou

Abstract: Prompt and $β$-delayed fission observables, such as the average number of prompt and delayed neutrons, the independent and cumulative fission product yields, and the prompt $γ$-ray energy spectra for the photonuclear reactions on $^{235,238}$U and $^{239}$Pu are calculated with the Hauser-Feshbach Fission Fragment Decay (HF$^3$D) model and compared with available experimental data. In the analysis… ▽ More Prompt and $β$-delayed fission observables, such as the average number of prompt and delayed neutrons, the independent and cumulative fission product yields, and the prompt $γ$-ray energy spectra for the photonuclear reactions on $^{235,238}$U and $^{239}$Pu are calculated with the Hauser-Feshbach Fission Fragment Decay (HF$^3$D) model and compared with available experimental data. In the analysis of neutron-induced fission reactions to the case of photo-induced fission, an excellent reproduction of the delayed neutron yields supports a traditional assumption that the photo-fission might be similar to the neutron-induced fission at the same excitation energies regardless of the spin and parity of the fissioning systems. △ Less

Submitted 17 December, 2022; originally announced December 2022.

Report number: LA-UR-22-32277

QRPA calculations for M1 transitions with the noniterative finite amplitude method and the application to neutron radiative capture cross sections

Authors: Hirokazu Sasaki, Toshihiko Kawano, Ionel Stetcu

Abstract: We derive the equations of quasiparticle random-phase approximation (QRPA) based on the finite amplitude method (FAM) with the Hartree-Fock+BCS (HF+BCS) single-particle states, and calculate the magnetic dipole (M1) transition for deformed gadolinium isotopes. Our QRPA calculation shows both large spin-flip transitions in the 5 to 10 MeV excitation energy and the low energy orbital transition that… ▽ More We derive the equations of quasiparticle random-phase approximation (QRPA) based on the finite amplitude method (FAM) with the Hartree-Fock+BCS (HF+BCS) single-particle states, and calculate the magnetic dipole (M1) transition for deformed gadolinium isotopes. Our QRPA calculation shows both large spin-flip transitions in the 5 to 10 MeV excitation energy and the low energy orbital transition that would correspond to the M1 scissors mode observed in nuclear experiments. Then, we calculate neutron capture reactions based on the statistical Hauser-Feshbach theory with the photoabsorption cross sections given by QRPA. We find that the capture cross section is enhanced due to the contribution from the low energy M1 transition although the calculated capture cross section still underestimates the experimental data. This issue in the calculated capture cross section could be improved by uncertainties of low energy E1 transition neglected in our QRPA calculation. △ Less

Submitted 29 November, 2022; originally announced November 2022.

Comments: 9 pages, 6 figures

Report number: LA-UR-22-32329

Nuclear data activities for medium mass and heavy nuclei at Los Alamos

Authors: M. R. Mumpower, T. M Sprouse, T. Kawano, M. W. Herman, A. E. Lovell, G. W. Misch, D. Neudecker, H. Sasaki, I. Stetcu, P. Talou

Abstract: Nuclear data is critical for many modern applications from stockpile stewardship to cutting edge scientific research. Central to these pursuits is a robust pipeline for nuclear modeling as well as data assimilation and dissemination. We summarize a small portion of the ongoing nuclear data efforts at Los Alamos for medium mass to heavy nuclei. We begin with an overview of the NEXUS framework and s… ▽ More Nuclear data is critical for many modern applications from stockpile stewardship to cutting edge scientific research. Central to these pursuits is a robust pipeline for nuclear modeling as well as data assimilation and dissemination. We summarize a small portion of the ongoing nuclear data efforts at Los Alamos for medium mass to heavy nuclei. We begin with an overview of the NEXUS framework and show how one of its modules can be used for model parameter optimization using Bayesian techniques. The mathematical framework affords the combination of different measured data in determining model parameters and their associated correlations. It also has the advantage of being able to quantify outliers in data. We exemplify the power of this procedure by highlighting the recently evaluated 239-Pu cross section. We further showcase the success of our tools and pipeline by covering the insight gained from incorporating the latest nuclear modeling and data in astrophysical simulations as part of the Fission In R-process Elements (FIRE) collaboration. △ Less

Submitted 21 October, 2022; originally announced October 2022.

Comments: 6 pages, 5 figures, Nuclear Data (2022) conference proceedings. Comments welcome!

Report number: LA-UR-22-30993

Collective enhancement in the exciton model

Authors: M. R. Mumpower, D. Nuedecker, H. Sasaki, T. Kawano, A. E. Lovell, M. W. Herman, I. Stetcu, M. Dupuis

Abstract: The pre-equilibrium reaction mechanism is considered in the context of the exciton model. A modification to the one-particle one-hole state density is studied which can be interpreted as a collective enhancement. The magnitude of the collective enhancement is set by simulating the Lawrence Livermore National Laboratory (LLNL) pulsed-spheres neutron-leakage spectra. The impact of the collective enh… ▽ More The pre-equilibrium reaction mechanism is considered in the context of the exciton model. A modification to the one-particle one-hole state density is studied which can be interpreted as a collective enhancement. The magnitude of the collective enhancement is set by simulating the Lawrence Livermore National Laboratory (LLNL) pulsed-spheres neutron-leakage spectra. The impact of the collective enhancement is explored in the context of the highly deformed actinide, 239-Pu. A consequence of this enhancement is the removal of fictitious levels in the Distorted-Wave Born Approximation often used in modern nuclear reaction codes. △ Less

Submitted 21 October, 2022; originally announced October 2022.

Comments: 7 pages, 7 figures. Comments welcome!

Report number: LA-UR-22-30081

Comprehensive Analyses of the Neutrino-Process in the Core-collapsing Supernova

Authors: Heamin Ko, Dukjae Jang, Myung-Ki Cheoun, Motohiko Kusakabe, Hirokazu Sasaki, Xingqun Yao, Toshitaka Kajino, Takehito Hayakawa, Masaomi Ono, Toshihiko Kawano, Grant J. Mathews

Abstract: We investigate the neutrino flavor change effects due to neutrino self-interaction, shock wave propagation as well as matter effect on the neutrino-process of the core-collapsing supernova (CCSN). For the hydrodynamics, we use two models: a simple thermal bomb model and a specified hydrodynamic model for SN1987A. As a pre-supernova model, we take an updated model adjusted to explain the SN1987A em… ▽ More We investigate the neutrino flavor change effects due to neutrino self-interaction, shock wave propagation as well as matter effect on the neutrino-process of the core-collapsing supernova (CCSN). For the hydrodynamics, we use two models: a simple thermal bomb model and a specified hydrodynamic model for SN1987A. As a pre-supernova model, we take an updated model adjusted to explain the SN1987A employing recent development of the $(n,γ)$ reaction rates for nuclei near the stability line $(A \sim 100)$. As for the neutrino luminosity, we adopt two different models: equivalent neutrino luminosity and non-equivalent luminosity models. The latter is taken from the synthetic analyses of the CCSN simulation data which involved quantitatively the results obtained by various neutrino transport models. Relevant neutrino-induced reaction rates are calculated by a shell model for light nuclei and a quasi-particle random phase approximation model for heavy nuclei. For each model, we present abundances of the light nuclei ($^7$Li, $^7$Be, $^{11}$B and $^{11}$C) and heavy nuclei ($^{92}$Nb, $^{98}$Tc, $^{138}$La and $^{180}$Ta) produced by the neutrino-process. The light nuclei abundances turn out to be sensitive to the Mikheyev-Smirnov-Wolfenstein (MSW) region around O-Ne-Mg region while the heavy nuclei are mainly produced prior to the MSW region. Through the detailed analyses, we find that neutrino self-interaction becomes a key ingredient in addition to the MSW effect for understanding the neutrino-process and the relevant nuclear abundances. The normal mass hierarchy is shown to be more compatible with the meteorite data. Main nuclear reactions for each nucleus are also investigated in detail. △ Less

Submitted 16 August, 2022; v1 submitted 24 March, 2022; originally announced March 2022.

Journal ref: Astrophys. J. 937, 116 (2022)

Non-iterative finite amplitude methods for E1 and M1 giant resonances

Authors: Hirokazu Sasaki, Toshihiko Kawano, Ionel Stetcu

Abstract: The finite amplitude method (FAM) is a very efficient approach for solving the fully self-consistent random-phase approximation (RPA) equations. We use FAM to rederive the RPA matrices for general Skyrme-like functionals, calculate the electric dipole (E1) and the magnetic dipole (M1) giant resonances, and compare the results with available experimental and evaluated data. For the E1 transitions i… ▽ More The finite amplitude method (FAM) is a very efficient approach for solving the fully self-consistent random-phase approximation (RPA) equations. We use FAM to rederive the RPA matrices for general Skyrme-like functionals, calculate the electric dipole (E1) and the magnetic dipole (M1) giant resonances, and compare the results with available experimental and evaluated data. For the E1 transitions in heavy nuclei, the calculations reproduce well the resonance energy of the photoabsorption cross sections. In the case of M1 transitions, we show that the residual interaction does not affect the transition strength of double-magic nuclei, which suggests that the spin terms in the Skyrme force currently neglected in the present computation could improve the agreement between FAM and experimental data. △ Less

Submitted 26 February, 2022; originally announced February 2022.

Comments: 15 pages, 6 figures

Report number: LA-UR-22-21601

Discrete Anomaly Matching for the Pouliot Type Dualities

Authors: Teruhiko Kawano

Abstract: We compute the 't Hooft anomalies of discrete symmetries in the Pouliot type dual theories and check their anomaly matching conditions. The Pouliot type dual theories we will consider in this paper are two dual pairs; the dual pair of ${\cal N}=1$ supersymmetric theories of a $Spin(7)$ gauge theory with spinors and a $SU(N_f-4)$ gauge theory with a symmetric tensor, fundamentals and singlets, and… ▽ More We compute the 't Hooft anomalies of discrete symmetries in the Pouliot type dual theories and check their anomaly matching conditions. The Pouliot type dual theories we will consider in this paper are two dual pairs; the dual pair of ${\cal N}=1$ supersymmetric theories of a $Spin(7)$ gauge theory with spinors and a $SU(N_f-4)$ gauge theory with a symmetric tensor, fundamentals and singlets, and the other dual pair of ${\cal N}=1$ supersymmetric theories of a $Spin(10)$ gauge theory with a spinor and vectors and a $SU(N_f-5)$ gauge theory with a symmetric tensor, fundamentals and singlets. We will show that the both pairs satisfy the discrete anomaly matching conditions. △ Less

Submitted 12 April, 2022; v1 submitted 8 February, 2022; originally announced February 2022.

Comments: 32 pages, typos corrected

Beta-delayed fission in the coupled Quasi-particle Random Phase Approximation plus Hauser-Feshbach approach

Authors: M. R. Mumpower, T. Kawano, T. M. Sprouse

Abstract: Beta-delayed neutron emission and $β$-delayed fission ($β$df) probabilities were calculated for heavy, neutron-rich nuclei using the Los Alamos coupled Quasi-Particle Random Phase Approximation plus Hauser-Feshbach (QRPA+HF) approach. In this model, the compound nucleus is initially populated by $β$-decay and is followed through subsequent statistical decays taking into account competition between… ▽ More Beta-delayed neutron emission and $β$-delayed fission ($β$df) probabilities were calculated for heavy, neutron-rich nuclei using the Los Alamos coupled Quasi-Particle Random Phase Approximation plus Hauser-Feshbach (QRPA+HF) approach. In this model, the compound nucleus is initially populated by $β$-decay and is followed through subsequent statistical decays taking into account competition between neutrons, $γ$-rays and fission. The primary output of these calculations includes branching ratios along with neutron and $γ$-ray spectra. We find a relatively large region of heavy nuclides where the probability of $β$df is near 100%. For a subset of nuclei near the neutron dripline, delayed neutron emission and the probability to fission are both large which leads to the possibility of multi-chance $β$df (mc-$β$df). We comment on prospective neutron-rich nuclei that could be probed by future experimental campaigns and provide a full table of branching ratios in ASCII format in the supplemental material for use in various applications. △ Less

Submitted 8 January, 2022; originally announced January 2022.

Comments: 12 pages, 7 figures

Report number: LA-UR-21-31648

Emergence of Lie group symmetric classical spacetimes in canonical tensor model

Authors: Taigen Kawano, Naoki Sasakura

Abstract: We analyze a wave function of a tensor model in the canonical formalism, when the argument of the wave function takes Lie group invariant or nearby values. Numerical computations show that there are two phases, which we call the quantum and the classical phases, respectively. In the classical phase, fluctuations are suppressed, and there emerge configurations which are discretizations of the class… ▽ More We analyze a wave function of a tensor model in the canonical formalism, when the argument of the wave function takes Lie group invariant or nearby values. Numerical computations show that there are two phases, which we call the quantum and the classical phases, respectively. In the classical phase, fluctuations are suppressed, and there emerge configurations which are discretizations of the classical geometric spaces invariant under the Lie group symmetries. This is explicitly demonstrated for the emergence of $S^n\ (n=1,2,3)$ for $SO(n+1)$ symmetries by checking the topological and the geometric (Laplacian) properties of the emerging configurations. The transition between the two phases has the form of splitting/merging of distributions of variables, resembling a matrix model counterpart, namely, the transition between one-cut and two-cut solutions. However this resemblance is obscured by a difference of the mechanism of the distribution in our setup from that in the matrix model. We also discuss this transition as a replica symmetry breaking. We perform various preliminary studies of the properties of the phases and the transition for such values of the argument. △ Less

Submitted 20 September, 2021; originally announced September 2021.

Comments: 32 pages, 17 figures

Report number: YITP-21-99

Journal ref: Prog Theor Exp Phys (2022)

Angular momentum removal by neutron and $γ$-ray emissions during fission fragment decays

Authors: I. Stetcu, A. E. Lovell, P. Talou, T. Kawano, S. Marin, S. A. Pozzi, A. Bulgac

Abstract: We investigate the angular momentum removal from fission fragments (FFs) through neutron and $γ$-ray emission, where we find that about half the neutrons are emitted with angular momenta $\ge 1.5\hbar$ and that the change in angular momentum after the emission of neutrons and statistical $γ$ rays is significant, contradicting usual assumptions. Per fission event, in our simulations, the neutron an… ▽ More We investigate the angular momentum removal from fission fragments (FFs) through neutron and $γ$-ray emission, where we find that about half the neutrons are emitted with angular momenta $\ge 1.5\hbar$ and that the change in angular momentum after the emission of neutrons and statistical $γ$ rays is significant, contradicting usual assumptions. Per fission event, in our simulations, the neutron and statistical $γ$-ray emissions change the spin of the fragment by 3.5 -- 5~$\hbar$, with a large standard deviation comparable to the average value. Such wide angular momentum removal distributions can hide any underlying correlations in the fission fragment initial spin values. Within our model, we reproduce data on spin measurements from discrete transitions after neutron emissions, especially in the case of light FFs. The agreement further improves for the heavy fragments if one removes from the analysis the events that would produce isomeric states. Finally, we show that while in our model the initial FF spins do not follow a saw-tooth like behavior observed in recent measurements, the average FF spin computed after neutron and statistical $γ$ emissions exhibits a shape that resembles a saw tooth. This suggests that the average FF spin measured after statistical emissions is not necessarily connected with the scission mechanism as previously implied. △ Less

Submitted 29 November, 2021; v1 submitted 9 August, 2021; originally announced August 2021.

Comments: 9 pages, 6 figures; version accepted for publication

Report number: LA-UR-21-23498, NT@UW-21-07

Journal ref: Phys. Rev. Lett. 127, 222502 (2021)

Influence of non-statistical properties in nuclear structure on emission of prompt fission neutrons

Authors: Toshihiko Kawano, Shin Okumura, Amy E. Lovell, Ionel Stetcu, Patrick Talou

Abstract: The Hauser-Feshbach Fission Fragment Decay (HF$^3$D) model is extended to calculate the prompt fission neutron spectrum (PFNS) for the thermal neutron induced fission on $^{235}$U, where the evaporated neutrons from all possible fission fragment pairs are aggregated. By studying model parameter sensitivities on the calculated PFNS, as well as non-statistical behavior of low-lying discrete level sp… ▽ More The Hauser-Feshbach Fission Fragment Decay (HF$^3$D) model is extended to calculate the prompt fission neutron spectrum (PFNS) for the thermal neutron induced fission on $^{235}$U, where the evaporated neutrons from all possible fission fragment pairs are aggregated. By studying model parameter sensitivities on the calculated PFNS, as well as non-statistical behavior of low-lying discrete level spin distribution, we conclude that discrepancies between the aggregation calculation and the experimental PFNS seen at higher neutron emission energies can be attributed to both the primary fission fragment yield distribution and the possible high spin states that are not predicted by the statistical theory of nuclear structure. △ Less

Submitted 2 April, 2021; originally announced April 2021.

Report number: {LA-UR-21-22687

Journal ref: Phys. Rev. C 104, 014611 (2021)

Energy Dependent Calculations of Fission Product, Prompt, and Delayed Neutron Yields for Neutron Induced Fission on $^{235}$U, $^{238}$U, and $^{239}$Pu

Authors: S. Okumura, T. Kawano, A. E. Lovell, T. Yoshida

Abstract: We perform energy dependent calculations of independent and cumulative fission product yields for $^{235}$U, $^{238}$U, and $^{239}$Pu in the first chance fission region. Starting with the primary fission fragment distributions taken from available experimental data and analytical functions based on assumptions for the excitation energy and spin-parity distributions, the Hauser-Feshbach statistica… ▽ More We perform energy dependent calculations of independent and cumulative fission product yields for $^{235}$U, $^{238}$U, and $^{239}$Pu in the first chance fission region. Starting with the primary fission fragment distributions taken from available experimental data and analytical functions based on assumptions for the excitation energy and spin-parity distributions, the Hauser-Feshbach statistical decay treatment for fission fragment de-excitation is applied to more than 1,000 fission fragments for the incident neutron energies up to 5 MeV. The calculated independent fission product yields are then used as an input of $β$-decay to produce the cumulative yield, and summation calculations are performed. Model parameters in these procedures are adjusted by applying the Bayesian technique at the thermal energy for $^{235}$U and $^{239}$Pu and in the fast energy range for $^{238}$U. The calculated fission observable quantities, such as the energy-dependent fission yields, and prompt and delayed neutron yields, are compared with available experimental data. We also study a possible impact of the second chance fission opening on the energy dependence of the delayed neutron yield by extrapolating the calculation. △ Less

Submitted 1 February, 2021; originally announced February 2021.

Report number: LA-UR-21-20820

Fission Fragment Decay Simulations with the CGMF Code

Authors: P. Talou, I. Stetcu, P. Jaffke, M. E. Rising, A. E. Lovell, T. Kawano

Abstract: The CGMF code implements the Hauser-Feshbach statistical nuclear reaction model to follow the de-excitation of fission fragments by successive emissions of prompt neutrons and $γ$ rays. The Monte Carlo technique is used to facilitate the analysis of complex distributions and correlations among the prompt fission observables. Starting from initial configurations for the fission fragments in mass, c… ▽ More The CGMF code implements the Hauser-Feshbach statistical nuclear reaction model to follow the de-excitation of fission fragments by successive emissions of prompt neutrons and $γ$ rays. The Monte Carlo technique is used to facilitate the analysis of complex distributions and correlations among the prompt fission observables. Starting from initial configurations for the fission fragments in mass, charge, kinetic energy, excitation energy, spin, and parity, $Y(A,Z,KE,U,J,π)$, CGMF samples neutron and $γ$-ray probability distributions at each stage of the decay process, conserving energy, spin and parity. Nuclear structure and reaction input data from the RIPL library are used to describe fission fragment properties and decay probabilities. Characteristics of prompt fission neutrons, prompt fission gamma rays, and independent fission yields can be studied consistently. Correlations in energy, angle and multiplicity among the emitted neutrons and $γ$ rays can be easily analyzed as a function of the emitting fragments. △ Less

Submitted 20 November, 2020; originally announced November 2020.

Report number: LA-UR-20-21264

Journal ref: Comput. Phys. Commun. 269, 108087 (2021)

Extension of the Hauser-Feshbach Fission Fragment Decay Model to Multi-Chance Fission

Authors: A. E. Lovell, T. Kawano, S. Okumura, I. Stetcu, M. R. Mumpower, P. Talou

Abstract: The Hauser-Feshbach fission fragment decay model, $\mathtt{HF^3D}$, which calculates the statistical decay of fission fragments, has been expanded to include multi-chance fission, up to neutron incident energies of 20 MeV. The deterministic decay takes as input pre-scission quantities - fission probabilities and the average energy causing fission - and post-scission quantities - yields in mass, ch… ▽ More The Hauser-Feshbach fission fragment decay model, $\mathtt{HF^3D}$, which calculates the statistical decay of fission fragments, has been expanded to include multi-chance fission, up to neutron incident energies of 20 MeV. The deterministic decay takes as input pre-scission quantities - fission probabilities and the average energy causing fission - and post-scission quantities - yields in mass, charge, total kinetic energy, spin, and parity. From these fission fragment initial conditions, the full decay is followed through both prompt and delayed particle emissions, allowing for the calculation of prompt neutron and $γ$ properties, such as multiplicity and energy distributions, both independent and cumulative fission yields, and delayed neutron observables. In this work, we describe the implementation of multi-chance fission into the $\mathtt{HF^3D}$ model, and show an example of prompt and delayed quantities beyond first-chance fission, using the example of neutron-induced fission on $^{235}$U. This expansion represents significant progress in consistently modeling the emission of prompt and delayed particles from fissile systems. △ Less

Submitted 26 October, 2020; originally announced October 2020.

Comments: 15 pages, 11 figures, submitted to PRC

Report number: LA-UR-20-27423

Unified description of the coupled-channels and statistical Hauser-Feshbach nuclear reaction theories for low energy neutron incident reactions

Authors: Toshihiko Kawano

Abstract: We incorporate the coupled-channels optical model into the statistical Hauser-Feshbach nuclear reaction theory, where the scattering matrix is diagonalized by performing the Engelbrecht-Weidenmüller transformation. This technique has been implemented in the coupled-channels optical model code ECIS by J. Raynal, and we extend this method so that all the open channels in a nucleon-induced reaction o… ▽ More We incorporate the coupled-channels optical model into the statistical Hauser-Feshbach nuclear reaction theory, where the scattering matrix is diagonalized by performing the Engelbrecht-Weidenmüller transformation. This technique has been implemented in the coupled-channels optical model code ECIS by J. Raynal, and we extend this method so that all the open channels in a nucleon-induced reaction on a deformed nucleus can be calculated consistently. △ Less

Submitted 26 September, 2020; originally announced September 2020.

Report number: LA-UR 20-27465

The Baghdad Atlas: A relational database of inelastic neutron-scattering $(n,n'γ)$ data

Authors: A. M. Hurst, L. A. Bernstein, T. Kawano, A. M. Lewis, K. Song

Abstract: A relational database has been developed based on the original ($n,n'γ$) work carried out by A. M. Demidov $et$ $al$., at the Nuclear Research Institute in Baghdad, Iraq [$"Atlas$ $of$ $Gamma$-$Ray$ $Spectra$ $from$ $the$ $Inelastic$ $Scattering$ $of$ $Reactor$ $Fast$ $Neutrons"$, Nuclear Research Institute, Baghdad, Iraq (Moscow, Atomizdat 1978)] for 105 independent measurements comprising 76 ele… ▽ More A relational database has been developed based on the original ($n,n'γ$) work carried out by A. M. Demidov $et$ $al$., at the Nuclear Research Institute in Baghdad, Iraq [$"Atlas$ $of$ $Gamma$-$Ray$ $Spectra$ $from$ $the$ $Inelastic$ $Scattering$ $of$ $Reactor$ $Fast$ $Neutrons"$, Nuclear Research Institute, Baghdad, Iraq (Moscow, Atomizdat 1978)] for 105 independent measurements comprising 76 elemental samples of natural composition and 29 isotopically-enriched samples. The information from this Atlas includes: $γ$-ray energies and relative intensities; nuclide and level data corresponding to the residual nucleus and meta data associated with the target sample that allows for the extraction of the flux-weighted ($n,n'γ$) cross sections for a given transition relative to a defined value. The optimized angular-distribution-corrected fast-neutron flux-weighted partial $γ$-ray cross section for the production of the 846.8-keV $2^{+}_{1} \rightarrow 0^{+}_{\rm gs}$ $γ$-ray transition in $^{56}$Fe, determined to be $\langle σ_γ \rangle = 143(29)$ mb, is used for this purpose. However, different values for the adopted cross section can be readily implemented to accommodate user preference based on revised determinations of this quantity. The Atlas ($n,n'γ$) data has been compiled into a series of CSV-style ASCII data sets and a suite of Python scripts have been developed to build and install the database locally. The database can then be accessed directly through the SQLite engine, or using alternative methods such as the Jupyter Notebook Python-browser interface. Several examples exploiting different interaction methodologies are distributed with the complete software package. △ Less

Submitted 15 January, 2021; v1 submitted 29 January, 2020; originally announced January 2020.

Comments: 24 pages, 14 figures

Characterizing $r$-Process Sites through Actinide Production

Authors: Erika M. Holmbeck, Rebecca Surman, Anna Frebel, G. C. McLaughlin, Matthew R. Mumpower, Trevor M. Sprouse, Toshihiko Kawano, Nicole Vassh, Timothy C. Beers

Abstract: Of the variations in the elemental abundance patterns of stars enhanced with $r$-process elements, the variation in the relative actinide-to-lanthanide ratio is among the most significant. We investigate the source of these actinide differences in order to determine whether these variations are due to natural differences in astrophysical sites, or due to the uncertain nuclear properties that are a… ▽ More Of the variations in the elemental abundance patterns of stars enhanced with $r$-process elements, the variation in the relative actinide-to-lanthanide ratio is among the most significant. We investigate the source of these actinide differences in order to determine whether these variations are due to natural differences in astrophysical sites, or due to the uncertain nuclear properties that are accessed in $r$-process sites. We find that variations between relative stellar actinide abundances is most likely astrophysical in nature, owing to how neutron-rich the ejecta from an $r$-process event may be. Furthermore, if an $r$-process site is capable of generating variations in the neutron-richness of its ejected material, then only one type of $r$-process site is needed to explain all levels of observed relative actinide enhancements. △ Less

Submitted 23 January, 2020; originally announced January 2020.

Comments: 12 pages, 5 figures, NPA-IX proceedings

Journal ref: Journal of Physics: Conference Series 1668, 012020 (2020)

Reference Database for Photon Strength Functions

Authors: S. Goriely, P. Dimitriou, M. Wiedeking, T. Belgya, R. Firestone, J. Kopecky, M. Krticka, V. Plujko, R. Schwengner, S. Siem, H. Utsunomiya, S. Hilaire, S. Peru, Y. S. Cho, D. M. Filipescu, N. Iwamoto, T. Kawano, V. Varlamov, R. Xu

Abstract: Photon strength functions describing the average response of the nucleus to an electromagnetic probe are key input information in the theoretical modelling of nuclear reactions. Consequently they are important for a wide range of fields such as nuclear structure, nuclear astrophysics, medical isotope production, fission and fusion reactor technologies. They are also sources of information for wide… ▽ More Photon strength functions describing the average response of the nucleus to an electromagnetic probe are key input information in the theoretical modelling of nuclear reactions. Consequently they are important for a wide range of fields such as nuclear structure, nuclear astrophysics, medical isotope production, fission and fusion reactor technologies. They are also sources of information for widely used reaction libraries such as the IAEA Reference Input Parameter Library and evaluated data files such as EGAF. In the past two decades, the amount of reaction gamma-ray data measured to determine photon strength functions has grown rapidly. Different experimental techniques have led to discrepant results and users are faced with the dilemma which (if any) of the divergent data to adopt. We report on a coordinated effort to compile and assess the existing experimental data on photon strength functions from the giant dipole resonance region to energies below the neutron separation energy. The assessment of the discrepant data at energies around or below the neutron separation energy has been possible only in a few cases where adequate information on the model-dependent analysis and estimation of uncertainties was available. In the giant dipole resonance region, we adopt the recommendations of the new IAEA photonuclear data library. We also present global empirical and semi-microscopic models that describe the photon strength functions in the entire energy region and reproduce reasonably well most of the experimental data. The compiled experimental photon strengths and recommended model calculations are available from the PSF database hosted at the IAEA (URL:www-nds.iaea.org/PSFdatabase). △ Less

Submitted 15 October, 2019; originally announced October 2019.

Comments: IAEA Coordinated Research Project on Generating Reference Database for Photon Strength Functions (F41032)

Journal ref: Eur. Phys. J. A (2019) 55:172

IAEA Photonuclear Data Library 2019

Authors: T. Kawano, Y. S. Cho, P. Dimitriou, D. Filipescu, N. Iwamoto, V. Plujko, X. Tao, H. Utsunomiya, V. Varlamov, R. Xu, R. Capote, I. Gheorghe, O. Gorbachenko, Y. L. Jin, T. Renstrøm, K. Stopani, Y. Tian, G. M. Tveten, J. M. Wang, T. Belgya, R. Firestone, S. Goriely, J. Kopecky, M. Krtička, R. Schwengner , et al. (2 additional authors not shown)

Abstract: Photo-induced reaction cross section data are of importance for a variety of current or emerging applications, such as radiation shielding design and radiation transport analyses, calculations of absorbed dose in the human body during radiotherapy, physics and technology of fission reactors (influence of photo-reactions on neutron balance) and fusion reactors (plasma diagnostics and shielding), ac… ▽ More Photo-induced reaction cross section data are of importance for a variety of current or emerging applications, such as radiation shielding design and radiation transport analyses, calculations of absorbed dose in the human body during radiotherapy, physics and technology of fission reactors (influence of photo-reactions on neutron balance) and fusion reactors (plasma diagnostics and shielding), activation analyses, safeguards and inspection technologies, nuclear waste transmutation, medical isotope production and astrophysical applications. Since the release of the IAEA Photonuclear Data Library in 1999 however, new experimental data as well as new methods to assess the reliability of experimental cross sections have become available. Theoretical models and input parameters used to evaluate photo-induced reactions have improved significantly over the years. In addition, new measurements of partial photoneutron cross sections using mono-energetic photon beams and advanced neutron detection systems have been performed allowing for the validation of the evaluations and assessments of the experimental data. Furthermore, technological advances have led to the construction of new and more powerful gamma-beam facilities, therefore new data needs are emerging. We report our coordinated efforts to address these data needs and present the results of the new evaluations of more than 200 nuclides included in the new updated IAEA Photonuclear Data Library, where the photon energy goes up to 200 MeV. We discuss the new assessment method and make recommendations to the user community in cases where the experimental data are discrepant and the assessments disagree. In addition, in the absence of experimental data, we present model predictions for photo-induced reaction cross section on nuclides of potential interest to medical radioisotope production. △ Less

Submitted 1 August, 2019; originally announced August 2019.

Report number: LA-UR-19-26964

Neutrino self-interaction and MSW effects on the supernova neutrino-process

Authors: Heamin Ko, Myung-Ki Cheoun, Eunja Ha, Motohiko Kusakabe, Takehito Hayakawa, Hirokazu Sasaki, Toshitaka Kajino, Masa-aki Hashimoto, Masaomi Ono, Mark D. Usang, Satoshi Chiba, Ko Nakamura, Alexey Tolstov, Ken'ichi Nomoto, Toshihiko Kawano, Grant J. Mathews

Abstract: We calculate the abundances of $^{7}$Li, $^{11}$B, $^{92}$Nb, $^{98}$Tc, $^{138}$La, and $^{180}$Ta produced by neutrino $(ν)$ induced reactions in a core-collapse supernova explosion. We consider the modification by $ν$ self-interaction ($ν$-SI) near the neutrinosphere and the Mikheyev-Smirnov-Wolfenstein effect in outer layers for time-dependent neutrino energy spectra. Abundances of $^{7}$Li an… ▽ More We calculate the abundances of $^{7}$Li, $^{11}$B, $^{92}$Nb, $^{98}$Tc, $^{138}$La, and $^{180}$Ta produced by neutrino $(ν)$ induced reactions in a core-collapse supernova explosion. We consider the modification by $ν$ self-interaction ($ν$-SI) near the neutrinosphere and the Mikheyev-Smirnov-Wolfenstein effect in outer layers for time-dependent neutrino energy spectra. Abundances of $^{7}$Li and heavy isotopes $^{92}$Nb, $^{98}$Tc and $^{138}$La are reduced by a factor of 1.5-2.0 by the $ν$-SI. In contrast, $^{11}$B is relatively insensitive to the $ν$-SI. We find that the abundance ratio of heavy to light nucleus, $^{138}$La/$^{11}$B, is sensitive to the neutrino mass hierarchy, and the normal mass hierarchy is more likely to be consistent with the solar abundances. △ Less

Submitted 23 November, 2019; v1 submitted 28 February, 2019; originally announced March 2019.

Journal ref: Astrophys. J. Lett. 891, L24 (2023)

Unified Coupled-Channels and Hauser-Feshbach Model Calculation for Nuclear Data Evaluation

Authors: Toshihiko Kawano

Abstract: We present an overview of the coupled-channels optical model and the Hauser-Feshbach theory code CoH$_3$, which focuses on the nuclear reaction calculations in the keV to tens of MeV region with special attention to the nuclear deformation. The code consists of three major sections that undertake the one-body potential mean-field theory, the coupled-channels optical model, and the Hauser-Feshbach… ▽ More We present an overview of the coupled-channels optical model and the Hauser-Feshbach theory code CoH$_3$, which focuses on the nuclear reaction calculations in the keV to tens of MeV region with special attention to the nuclear deformation. The code consists of three major sections that undertake the one-body potential mean-field theory, the coupled-channels optical model, and the Hauser-Feshbach statistical decay. There are other complementary segments to perform the whole nuclear reaction calculations, such as the direct/semidirect radiative capture process, pre-equilibrium process, and prompt fission neutron emission. △ Less

Submitted 17 January, 2019; originally announced January 2019.

Report number: LA-UR 19-20315

Transmission coefficients in compound-nucleus reaction theory

Authors: Y. Alhassid, G. F. Bertsch, P. Fanto, T. Kawano

Abstract: A recent article (D.A. Brown, et al., Phys. Rev. C98 024616 (2018)) rproposed a modification of the cross section formula used in practical calculations of compound nucleus reactions. We discuss the main concepts and approximations of statistical reaction theory and conclude that the standard practical implementations of the cross section formula remain the preferred choice. A recent article (D.A. Brown, et al., Phys. Rev. C98 024616 (2018)) rproposed a modification of the cross section formula used in practical calculations of compound nucleus reactions. We discuss the main concepts and approximations of statistical reaction theory and conclude that the standard practical implementations of the cross section formula remain the preferred choice. △ Less

Submitted 14 January, 2019; originally announced January 2019.

Comments: To be submitted to Phys. Rev. C

Journal ref: Phys. Rev. C 99, 024621 (2019)

Accretion geometry of the black-hole binary Cygnus X-1 from X-ray polarimetry

Authors: M. Chauvin, H. -G. Florén, M. Friis, M. Jackson, T. Kamae, J. Kataoka, T. Kawano, M. Kiss, V. Mikhalev, T. Mizuno, N. Ohashi, T. Stana, H. Tajima, H. Takahashi, N. Uchida, M. Pearce

Abstract: Black-hole binary (BHB) systems comprise a stellar-mass black hole and a closely orbiting companion star. Matter is transferred from the companion to the black hole, forming an accretion disk, corona and jet structures. The resulting release of gravitational energy leads to emission of X-rays. The radiation is affected by special/general relativistic effects, and can serve as a probe of the proper… ▽ More Black-hole binary (BHB) systems comprise a stellar-mass black hole and a closely orbiting companion star. Matter is transferred from the companion to the black hole, forming an accretion disk, corona and jet structures. The resulting release of gravitational energy leads to emission of X-rays. The radiation is affected by special/general relativistic effects, and can serve as a probe of the properties of the black hole and surrounding environment, if the accretion geometry is properly identified. Two competing models describe the disk-corona geometry for the hard spectral state of BHBs, based on spectral and timing measurements. Measuring the polarization of hard X-rays reflected from the disk allows the geometry to be determined. The extent of the corona differs between the two models, affecting the strength of relativistic effects (e.g., enhancement of polarization fraction and rotation of polarization angle). Here, we report observational results on linear polarization of hard X-ray (19-181 keV) emission from a BHB, Cygnus X-1, in the hard state. The low polarization fraction, <8.6% (upper limit at 90% confidence level), and the alignment of the polarization angle with the jet axis show that the dominant emission is not influenced by strong gravity. When considered together with existing spectral and timing data, our result reveals that the accretion corona is either an extended structure, or is located far from the black hole in the hard state of Cygnus X-1. △ Less

Submitted 24 December, 2018; originally announced December 2018.

Comments: 44 pages, 15 figures, 1 table

Journal ref: Nature Astronomy volume 2, pages 652-655 (2018)

High resolution measurement of tagged two-neutron energy and angle correlations in Cf-252(sf)

Authors: P. F. Schuster, M. J. Marcath, S. Marin, S. D. Clarke, M. Devlin, R. C. Haight, R. Vogt, P. Talou, I. Stetcu, T. Kawano, J. Randrup, S. A. Pozzi

Abstract: Background: Spontaneous fission events emit prompt neutrons correlated with one another in emission angle and energy. Purpose: We explore the relationship in energy and angle between correlated prompt neutrons emitted from 252Cf spontaneous fission. Methods: Measurements with the Chi-Nu array provide experimental data for coincident neutrons tagged with a fission chamber signal with 10 degree angu… ▽ More Background: Spontaneous fission events emit prompt neutrons correlated with one another in emission angle and energy. Purpose: We explore the relationship in energy and angle between correlated prompt neutrons emitted from 252Cf spontaneous fission. Methods: Measurements with the Chi-Nu array provide experimental data for coincident neutrons tagged with a fission chamber signal with 10 degree angular resolution and 1 ns timing resolution for time-of-flight energy calculations. The experimental results are compared to simulations produced by the fission event generators CGMF, FREYA, and MCNPX-POLIMI IPOL(1)=1. Results: We find that the measurements and the simulations all exhibit anisotropic neutron emission, though differences exist between fission event generators. Conclusions: This work shows that the dependence of detected neutron energy on the energy of a neutron detected in coincidence, although weak, is non-negligible, indicating that there may be correlations in energy between two neutrons emitted in the same fission event. △ Less

Submitted 16 November, 2018; originally announced November 2018.

Journal ref: Phys. Rev. C 100, 014605 (2019)

Number of Particles in Fission Fragments

Authors: Marc Verriere, Nicolas Schunck, Toshihiko Kawano

Abstract: In current simulations of fission, the number of protons and neutrons in a given fission fragment is almost always obtained by integrating the total density of particles in the sector of space that contains the fragment. Because of the antisymmetry of the many-body wave function of the whole nucleus, this procedure systematically gives noninteger numbers of particles in the fragments. We introduce… ▽ More In current simulations of fission, the number of protons and neutrons in a given fission fragment is almost always obtained by integrating the total density of particles in the sector of space that contains the fragment. Because of the antisymmetry of the many-body wave function of the whole nucleus, this procedure systematically gives noninteger numbers of particles in the fragments. We introduce a novel sampling method to estimate rigorously the probability of finding $Z$ protons and $N$ neutrons in a fission fragment without resorting to projectors, which can sometimes give unwieldy results. When applied on standard Hartree-Fock-Bogoliubov many-body states, we show that our approach reproduces indeed the results of full particle-number projection. We then estimate the charge and mass number dispersion of several scission configurations in $^{240}$Pu with and without pairing correlations included. We show that odd-even effects in the charge probability naturally occur within our approach, which could explain the well-known odd-even staggering of charge distributions. Our method is applicable either in static calculations of scission configurations such as, e.g., in the macroscopic-microscopic approach or energy density functional theory, but also in explicitly time-dependent density functional theory simulations of fission. △ Less

Submitted 17 August, 2020; v1 submitted 13 November, 2018; originally announced November 2018.

Report number: LLNL-JRNL-760851

Journal ref: Phys. Rev. C 100, 024612 (2019)

Angular momentum of fission fragments

Authors: G. F. Bertsch, T. Kawano, LM. Robledo

Abstract: The strong deformation present immediately after scission has consequences for the angular momentum population of the fragments as well as the angular distribution of their decay radiation. We find that the usual spin-cutoff parameterization describes very well the angular momentum distribution associated with the deformation of the fragments at the scission point. Depending on the deformation, it… ▽ More The strong deformation present immediately after scission has consequences for the angular momentum population of the fragments as well as the angular distribution of their decay radiation. We find that the usual spin-cutoff parameterization describes very well the angular momentum distribution associated with the deformation of the fragments at the scission point. Depending on the deformation, its contribution can comparable to the thermal contribution to the angular momentum of the newly formed fragments. The $M$-distribution of the angular momentum is highly polarized and gives rise to large anisotropies in the subsequent gamma cascade. We treat in detail a typical gamma cascade in a daughter nucleus, following usual model assumptions except for the anisotropy of the initial state. In principle, the observed anisotropy can provide information on the relative amounts of deformation and thermal energy present at the scission point. △ Less

Submitted 31 October, 2018; originally announced October 2018.

Comments: 5 pages with 4 figures

Journal ref: Phys. Rev. C 99, 034603 (2019)

Actinide Production in Neutron-Rich Ejecta of a Neutron Star Merger

Authors: Erika M. Holmbeck, Rebecca Surman, Trevor M. Sprouse, Matthew R. Mumpower, Nicole Vassh, Timothy C. Beers, Toshihiko Kawano

Abstract: The rapid-neutron-capture ("r") process is responsible for synthesizing many of the heavy elements observed in both the solar system and Galactic metal-poor halo stars. Simulations of r-process nucleosynthesis can reproduce abundances derived from observations with varying success, but so far fail to account for the observed over-enhancement of actinides, present in about 30% of r-process-enhanced… ▽ More The rapid-neutron-capture ("r") process is responsible for synthesizing many of the heavy elements observed in both the solar system and Galactic metal-poor halo stars. Simulations of r-process nucleosynthesis can reproduce abundances derived from observations with varying success, but so far fail to account for the observed over-enhancement of actinides, present in about 30% of r-process-enhanced stars. In this work, we investigate actinide production in the dynamical ejecta of a neutron star merger and explore if varying levels of neutron richness can reproduce the actinide boost. We also investigate the sensitivity of actinide production on nuclear physics properties: fission distribution, beta-decay, and mass model. For most cases, the actinides are over-produced in our models if the initial conditions are sufficiently neutron-rich for fission cycling. We find that actinide production can be so robust in the dynamical ejecta that an additional lanthanide-rich, actinide-poor component is necessary in order to match observations of actinide-boost stars. We present a simple actinide-dilution model that folds in estimated contributions from two nucleosynthetic sites within a merger event. Our study suggests that while the dynamical ejecta of a neutron star merger is a likely production site for the formation of actinides, a significant contribution from another site or sites (e.g., the neutron star merger accretion disk wind) is required to explain abundances of r-process-enhanced, metal-poor stars. △ Less

Submitted 18 October, 2018; v1 submitted 17 July, 2018; originally announced July 2018.

Comments: 15 pages, 19 figures

Californium-254 and kilonova light curves

Authors: Y. Zhu, R. T. Wollaeger, N. Vassh, R. Surman, T. M. Sprouse, M. R. Mumpower, P. Moller, G. C. McLaughlin, O. Korobkin, T. Kawano, P. J. Jaffke, E. M. Holmbeck, C. L. Fryer, W. P. Even, A. J. Couture, J. Barnes

Abstract: Neutron star mergers offer unique conditions for the creation of the heavy elements and additionally provide a testbed for our understanding of this synthesis known as the $r$-process. We have performed dynamical nucleosynthesis calculations and identified a single isotope, $^{254}$Cf, which has a particularly high impact on the brightness of electromagnetic transients associated with mergers on t… ▽ More Neutron star mergers offer unique conditions for the creation of the heavy elements and additionally provide a testbed for our understanding of this synthesis known as the $r$-process. We have performed dynamical nucleosynthesis calculations and identified a single isotope, $^{254}$Cf, which has a particularly high impact on the brightness of electromagnetic transients associated with mergers on the order of 15 to 250 days. This is due to the anomalously long half-life of this isotope and the efficiency of fission thermalization compared to other nuclear channels. We estimate the fission fragment yield of this nucleus and outline the astrophysical conditions under which $^{254}$Cf has the greatest impact to the light curve. Future observations in the middle-IR which are bright during this regime could indicate the production of actinide nucleosynthesis. △ Less

Submitted 25 June, 2018; originally announced June 2018.

Comments: 9 pages, 5 figures

Report number: LA-UR-18-25559

Canonical tensor model through data analysis -- Dimensions, topologies, and geometries --

Authors: Taigen Kawano, Dennis Obster, Naoki Sasakura

Abstract: The canonical tensor model (CTM) is a tensor model in Hamilton formalism and is studied as a model for gravity in both classical and quantum frameworks. Its dynamical variables are a canonical conjugate pair of real symmetric three-index tensors, and a question in this model was how to extract spacetime pictures from the tensors. We give such an extraction procedure by using two techniques widely… ▽ More The canonical tensor model (CTM) is a tensor model in Hamilton formalism and is studied as a model for gravity in both classical and quantum frameworks. Its dynamical variables are a canonical conjugate pair of real symmetric three-index tensors, and a question in this model was how to extract spacetime pictures from the tensors. We give such an extraction procedure by using two techniques widely known in data analysis. One is the tensor-rank (or CP, etc.) decomposition, which is a certain generalization of the singular value decomposition of a matrix and decomposes a tensor into a number of vectors. By regarding the vectors as points forming a space, topological properties can be extracted by using the other data analysis technique called persistent homology, and geometries by virtual diffusion processes over points. Thus, time evolutions of the tensors in the CTM can be interpreted as topological and geometric evolutions of spaces. We have performed some initial investigations of the classical equation of motion of the CTM in terms of these techniques for a homogeneous fuzzy circle and homogeneous two- and three-dimensional fuzzy spheres as spaces, and have obtained agreement with the general relativistic system obtained previously in a formal continuum limit of the CTM. It is also demonstrated by some concrete examples that the procedure is general for any dimensions and topologies, showing the generality of the CTM. △ Less

Submitted 20 June, 2018; v1 submitted 12 May, 2018; originally announced May 2018.

Comments: 44 pages, 16 figures, minor corrections

Report number: YITP-18-38

Journal ref: Phys. Rev. D 97, 124061 (2018)

Excitation functions for (p,x) reactions of niobium in the energy range of E$_{\text{p}}$ = 40-90 MeV

Authors: Andrew S. Voyles, Lee A. Bernstein, Eva R. Birnbaum, Jonathan W. Engle, Stephen A. Graves, Toshihiko Kawano, Amanda M. Lewis, Francois M. Nortier

Abstract: A stack of thin Nb foils was irradiated with the 100 MeV proton beam at Los Alamos National Laboratory's Isotope Production Facility, to investigate the $^{93}$Nb(p,4n)$^{90}$Mo nuclear reaction as a monitor for intermediate energy proton experiments and to benchmark state-of-the-art reaction model codes. A set of 38 measured cross sections for $^{\text{nat}}$Nb(p,x) and $^{\text{nat}}$Cu(p,x) rea… ▽ More A stack of thin Nb foils was irradiated with the 100 MeV proton beam at Los Alamos National Laboratory's Isotope Production Facility, to investigate the $^{93}$Nb(p,4n)$^{90}$Mo nuclear reaction as a monitor for intermediate energy proton experiments and to benchmark state-of-the-art reaction model codes. A set of 38 measured cross sections for $^{\text{nat}}$Nb(p,x) and $^{\text{nat}}$Cu(p,x) reactions between 40-90 MeV, as well as 5 independent measurements of isomer branching ratios, are reported. These are useful in medical and basic science radionuclide productions at intermediate energies. The $^{\text{nat}}$Cu(p,x)$^{56}$Co, $^{\text{nat}}$Cu(p,x)$^{62}$Zn, and $^{\text{nat}}$Cu(p,x)$^{65}$Zn reactions were used to determine proton fluence, and all activities were quantified using HPGe spectrometry. Variance minimization techniques were employed to reduce systematic uncertainties in proton energy and fluence, improving the reliability of these measurements. The measured cross sections are shown to be in excellent agreement with literature values, and have been measured with improved precision compared with previous measurements. This work also reports the first measurement of the $^{\text{nat}}$Nb(p,x)$^{82\text{m}}$Rb reaction, and of the independent cross sections for $^{\text{nat}}$Cu(p,x)$^{52\text{g}}$Mn and $^{\text{nat}}$Nb(p,x)$^{85\text{g}}$Y in the 40-90 MeV region. The effects of $^{\text{nat}}$Si(p,x)$^{22,24}$Na contamination, arising from silicone adhesive in the Kapton tape used to encapsulate the aluminum monitor foils, is also discussed as a cautionary note to future stacked-target cross section measurements. \emph{A priori} predictions of the reaction modeling codes CoH, EMPIRE, and TALYS are compared with experimentally measured values and used to explore the differences between codes for the $^{\text{nat}}$Nb(p,x) and $^{\text{nat}}$Cu(p,x) reactions. △ Less

Submitted 21 June, 2018; v1 submitted 18 April, 2018; originally announced April 2018.

Comments: 34 pages, submitted to NIM-B

Report number: LA-UR-18-22980

Journal ref: Nuclear Instruments and Methods in Physics Research B, 429 (2018) 53-74

Orbital solution leading to an acceptable interpretation for the enigmatic gamma-ray binary HESS J0632+057

Authors: Yuki Moritani, Takafumi Kawano, Sho Chimasu, Akiko Kawachi, Hiromitsu Takahashi, Jumpei Takata, Alex C. Carciofi

Abstract: High-dispersion spectroscopic monitoring of HESS J0632+057 has been carried out over four orbital cycles in order to search for orbital modulation, covering the entire orbital phase. We have measured radial velocity of H$α$ emission line with the method introduced by Shafter et al. (1986), which has been successfully applied to some Be stars. The velocity is seen to increase much earlier than expe… ▽ More High-dispersion spectroscopic monitoring of HESS J0632+057 has been carried out over four orbital cycles in order to search for orbital modulation, covering the entire orbital phase. We have measured radial velocity of H$α$ emission line with the method introduced by Shafter et al. (1986), which has been successfully applied to some Be stars. The velocity is seen to increase much earlier than expected for the orbital period of 315 days, and much more steeply than expected at around "apastron". The period of the H$α$ modulation is found to be as $308^{+26}_{-23}$ days. We have also analyzed Swift/XRT data from 2009 to 2015 to study the orbital modulation, selecting the data with good statistics ($\geq$ 30 counts). With additional two-year data to the previous works, the orbital period has been updated to $313^{+11}_{-8}$ days, which is consistent with the previous X-ray periods and the spectroscopic one. The past XMM-Newton and Chandra observations prefer the period of 313 days. With the new period, assuming that H$α$ velocities accurately trace the motion of the Be star, we have derived a new set of the orbital parameters. In the new orbit, which is less eccentric ($e \sim 0.6$), two outbursts occur after apastron, and just after periastron. Besides, the column density in bright phase ($4.7^{+0.9}_{-08}\times10^{21}\;\mathrm{cm^{-2}}$) is higher than in faint phase ($2.2\pm0.5\times10^{21}\;\mathrm{cm^{-2}}$). These facts suggest that outbursts occur when the compact object passes nearby/through the Be disk. The mass function implies that mass of the compact object is less than 2.5 $\mathrm{M_{sun}}$ assuming that the mass of the Be star is 13.2--18.2 $\mathrm{M_{sun}}$ (Aragona at al. 2010) unless the inclination is extremely small. The photon index indicates that the spectra becomes softer when the system is bright. These suggest that the compact object is a pulsar. △ Less

Submitted 11 April, 2018; originally announced April 2018.

Comments: 15 pages, 7 figures, accepted to PASJ

In-flight performance of the Canadian Astro-H Metrology System

Authors: Luigi C. Gallo, Alexander Koujelev, Stephane Gagnon, Timothy Elgin, Martin Guibert, Ryo Iizuka, Manabu Ishida, Kosei Ishimura, Naoko Iwata, Taro Kawano, Casey Lambert, Franco Moroso, Shiro Ueno, Takayuki Yuasa

Abstract: The Canadian Astro-H Metrology System (CAMS) on the Hitomi X-ray satellite is a laser alignment system that measures the lateral displacement (X/Y) of the extensible optical bench (EOB) along the optical axis of the hard X-ray telescopes (HXTs). The CAMS consists of two identical units that together can be used to discern translation and rotation of the deployable element along the axis. This pape… ▽ More The Canadian Astro-H Metrology System (CAMS) on the Hitomi X-ray satellite is a laser alignment system that measures the lateral displacement (X/Y) of the extensible optical bench (EOB) along the optical axis of the hard X-ray telescopes (HXTs). The CAMS consists of two identical units that together can be used to discern translation and rotation of the deployable element along the axis. This paper presents the results of in-flight usage of the CAMS during deployment of the EOB and during two observations (Crab and G21.5-0.9) with the HXTs. The CAMS was extremely important during the deployment operation by providing real-time positioning information of the EOB with micrometer scale resolution. In this work, we show how the CAMS improves data quality coming from the hard X-ray imagers. Moreover, we demonstrate that a metrology system is even more important as the angular resolution of the telescope increases. Such a metrology system will be an indispensable tool for future high resolution X-ray imaging missions. △ Less

Submitted 28 March, 2018; originally announced March 2018.

Comments: 46 pages, 24 figures. Accepted for publication in JATIS

Journal ref: J. Astron. Telesc. Instrum. Syst. 4(2), 021405 (2018)

The PoGO+ view on Crab off-pulse hard X-ray polarisation

Authors: M. Chauvin, H. -G. Florén, M. Friis, M. Jackson, T. Kamae, J. Kataoka, T. Kawano, M. Kiss, V. Mikhalev, T. Mizuno, H. Tajima, H. Takahashi, N. Uchida, M. Pearce

Abstract: The linear polarisation fraction and angle of the hard X-ray emission from the Crab provide unique insight into high energy radiation mechanisms, complementing the usual imaging, timing and spectroscopic approaches. Results have recently been presented by two missions operating in partially overlapping energy bands, PoGO+ (18-160 keV) and AstroSat CZTI (100-380 keV). We previously reported PoGO+ r… ▽ More The linear polarisation fraction and angle of the hard X-ray emission from the Crab provide unique insight into high energy radiation mechanisms, complementing the usual imaging, timing and spectroscopic approaches. Results have recently been presented by two missions operating in partially overlapping energy bands, PoGO+ (18-160 keV) and AstroSat CZTI (100-380 keV). We previously reported PoGO+ results on the polarisation parameters integrated across the light-curve and for the entire nebula-dominated off-pulse region. We now introduce finer phase binning, in light of the AstroSat CZTI claim that the polarisation fraction varies across the off-pulse region. Since both missions are operating in a regime where errors on the reconstructed polarisation parameters are non-Gaussian, we adopt a Bayesian approach to compare results from each mission. We find no statistically significant variation in off-pulse polarisation parameters, neither when considering the mission data separately nor when they are combined. This supports expectations from standard high-energy emission models. △ Less

Submitted 21 February, 2018; originally announced February 2018.

Comments: 6 pages, 3 figures, 5 tables. Accepted for publication in MNRAS Letters

$β$-delayed fission in $r$-process nucleosynthesis

Authors: M. R. Mumpower, T. Kawano, T. M. Sprouse, N. Vassh, E. M. Holmbeck, R. Surman, P. Moller

Abstract: We present $β$-delayed neutron emission and $β$-delayed fission calculations for heavy, neutron-rich nuclei using the coupled Quasi-Particle Random Phase Approximation plus Hauser-Feshbach (QRPA+HF) approach. From the initial population of a compound nucleus after $β$-decay, we follow the statistical decay taking into account competition between neutrons, $γ$-rays, and fission. We find a region of… ▽ More We present $β$-delayed neutron emission and $β$-delayed fission calculations for heavy, neutron-rich nuclei using the coupled Quasi-Particle Random Phase Approximation plus Hauser-Feshbach (QRPA+HF) approach. From the initial population of a compound nucleus after $β$-decay, we follow the statistical decay taking into account competition between neutrons, $γ$-rays, and fission. We find a region of the chart of nuclides where the probability of $β$-delayed fission is $\sim100$%, that likely prevents the production of superheavy elements in nature. For a subset of nuclei near the neutron dripline, neutron multiplicity and the probability of fission are both large, leading to the intriguing possibility of multi-chance $β$-delayed fission, a new decay mode for extremely neutron-rich heavy nuclei. In this new decay mode, $β$-decay can be followed by multiple neutron emission leading to subsequent daughter generations which each have a probability to fission. We explore the impact of $β$-delayed fission in rapid neutron capture process ($r$-process) nucleosynthesis in the tidal ejecta of a neutron star--neutron star merger and show that it is a key fission channel that shapes the final abundances near the second $r$-process peak. △ Less

Submitted 12 February, 2018; originally announced February 2018.

Comments: 9 pages, 5 figures, submitted

Report number: LA-UR-18-21034

$^{235}$U(n, f) Independent Fission Product Yield and Isomeric Ratio Calculated with the Statistical Hauser-Feshbach Theory

Authors: Shin Okumura, Toshihiko Kawano, Patrick Talou, Patrick Jaffke, Satoshi Chiba

Abstract: We have developed a Hauser-Feshbach fission fragment decay model, HF$^3$D, which can be applied to the statistical decay of more than 500 primary fission fragment pairs (1,000 nuclides) produced by the neutron induced fission of $^{235}$U. The fission fragment yield $Y(A)$ and the total kinetic energy TKE are model inputs, and we estimate them from available experimental data for the $^{235}$U(n… ▽ More We have developed a Hauser-Feshbach fission fragment decay model, HF$^3$D, which can be applied to the statistical decay of more than 500 primary fission fragment pairs (1,000 nuclides) produced by the neutron induced fission of $^{235}$U. The fission fragment yield $Y(A)$ and the total kinetic energy TKE are model inputs, and we estimate them from available experimental data for the $^{235}$U(n$\rm_{th}$,f) system. The model parameters in the statistical decay calculation are adjusted to reproduce some fission observables, such as the neutron emission multiplicity $\overlineν$, its distribution $P(ν)$, and the mass dependence $\overlineν(A)$. The calculated fission product yield and isomeric ratio are compared with experimental data. We show that the calculated independent fission product yield $Y_I(A)$ at the thermal energy reproduces the experimental data well, while the calculated isomeric ratios tend to be lower than the Madland-England model prediction. The model is extended to higher incident neutron energies up to the second chance fission threshold. We demonstrate for the first time that most of the isomeric ratios stay constant, although the production of isomeric state itself changes as the incident energy increases. △ Less

Submitted 4 February, 2018; originally announced February 2018.

Correlated Prompt Fission Data in Transport Simulations

Authors: P. Talou, R. Vogt, J. Randrup, M. E. Rising, S. A. Pozzi, J. Verbeke, M. T. Andrews, S. D. Clarke, P. Jaffke, M. Jandel, T. Kawano, M. J. Marcath, K. Meierbachtol, L. Nakae, G. Rusev, A. Sood, I. Stetcu, C. Walker

Abstract: Detailed information on the fission process can be inferred from the observation, modeling and theoretical understanding of prompt fission neutron and $γ$-ray~observables. Beyond simple average quantities, the study of distributions and correlations in prompt data, e.g., multiplicity-dependent neutron and \gray~spectra, angular distributions of the emitted particles, $n$-$n$, $n$-$γ$, and $γ$-$γ$~… ▽ More Detailed information on the fission process can be inferred from the observation, modeling and theoretical understanding of prompt fission neutron and $γ$-ray~observables. Beyond simple average quantities, the study of distributions and correlations in prompt data, e.g., multiplicity-dependent neutron and \gray~spectra, angular distributions of the emitted particles, $n$-$n$, $n$-$γ$, and $γ$-$γ$~correlations, can place stringent constraints on fission models and parameters that would otherwise be free to be tuned separately to represent individual fission observables. The FREYA~and CGMF~codes have been developed to follow the sequential emissions of prompt neutrons and $γ$-rays~from the initial excited fission fragments produced right after scission. Both codes implement Monte Carlo techniques to sample initial fission fragment configurations in mass, charge and kinetic energy and sample probabilities of neutron and $γ$~emission at each stage of the decay. This approach naturally leads to using simple but powerful statistical techniques to infer distributions and correlations among many observables and model parameters. The comparison of model calculations with experimental data provides a rich arena for testing various nuclear physics models such as those related to the nuclear structure and level densities of neutron-rich nuclei, the $γ$-ray~strength functions of dipole and quadrupole transitions, the mechanism for dividing the excitation energy between the two nascent fragments near scission, and the mechanisms behind the production of angular momentum in the fragments, etc. Beyond the obvious interest from a fundamental physics point of view, such studies are also important for addressing data needs in various nuclear applications. (See text for full abstract.) △ Less

Submitted 1 February, 2018; v1 submitted 29 September, 2017; originally announced October 2017.

Comments: 39 pages, 57 figure files, published in Eur. Phys. J. A, reference added this version

Report number: LA-UR-17-28181, LLNL-JRNL-738697

Journal ref: Eur. Phys. J. A (2018) 54: 9

Shedding new light on the Crab with polarized X-rays

Authors: M. Chauvin, H. -G. Florén, M. Friis, M. Jackson, T. Kamae, J. Kataoka, T. Kawano, M. Kiss, V. Mikhalev, T. Mizuno, N. Ohashi, T. Stana, H. Tajima, H. Takahashi, N. Uchida, M. Pearce

Abstract: Strong magnetic fields, synchrotron emission, and Compton scattering are omnipresent in compact celestial X-ray sources. Emissions in the X-ray energy band are consequently expected to be linearly polarized. X-ray polarimetry provides a unique diagnostic to study the location and fundamental mechanisms behind emission processes. The polarization of emissions from a bright celestial X-ray source, t… ▽ More Strong magnetic fields, synchrotron emission, and Compton scattering are omnipresent in compact celestial X-ray sources. Emissions in the X-ray energy band are consequently expected to be linearly polarized. X-ray polarimetry provides a unique diagnostic to study the location and fundamental mechanisms behind emission processes. The polarization of emissions from a bright celestial X-ray source, the Crab, is reported here for the first time in the hard X-ray band (~20-160 keV). The Crab is a complex system consisting of a central pulsar, a diffuse pulsar wind nebula, as well as structures in the inner nebula including a jet and torus. Measurements are made by a purpose-built and calibrated polarimeter, PoGO+. The polarization vector is found to be aligned with the spin axis of the pulsar for a polarization fraction, PF = (20.9 $\pm$ 5.0)%. This is higher than that of the optical diffuse nebula, implying a more compact emission site, though not as compact as, e.g., the synchrotron knot. Contrary to measurements at higher energies, no significant temporal evolution of phase-integrated polarisation parameters is observed. The polarization parameters for the pulsar itself are measured for the first time in the X-ray energy band and are consistent with observations at optical wavelengths. △ Less

Submitted 28 June, 2017; originally announced June 2017.

Comments: Accepted for publication in Scientific Reports, www.nature.com/srep

Estimation of M1 scissors mode strength for deformed nuclei in the medium to heavy mass region by statistical Hauser-Feshbach model calculations

Authors: M. R. Mumpower, T. Kawano, J. L. Ullmann, M. Krtička, T. M. Sprouse

Abstract: Radiative neutron capture is an important nuclear reaction whose accurate description is needed for many applications ranging from nuclear technology to nuclear astrophysics. The description of such a process relies on the Hauser-Feshbach theory which requires the nuclear optical potential, level density and $γ$-strength function as model inputs. It has recently been suggested that the M1 scissors… ▽ More Radiative neutron capture is an important nuclear reaction whose accurate description is needed for many applications ranging from nuclear technology to nuclear astrophysics. The description of such a process relies on the Hauser-Feshbach theory which requires the nuclear optical potential, level density and $γ$-strength function as model inputs. It has recently been suggested that the M1 scissors mode may explain discrepancies between theoretical calculations and evaluated data. We explore statistical model calculations with the strength of the M1 scissors mode estimated to be dependent on the nuclear deformation of the compound system. We show that the form of the M1 scissors mode improves the theoretical description of evaluated data and the match to experiment in both the fission product and actinide regions. Since the scissors mode occurs in the range of a few keV $\sim$ a few MeV, it may also impact the neutron capture cross sections of neutron-rich nuclei that participate in the rapid neutron capture process of nucleosynthesis. We comment on the possible impact to nucleosynthesis by evaluating neutron capture rates for neutron-rich nuclei with the M1 scissors mode active. △ Less

Submitted 22 June, 2017; originally announced June 2017.

Comments: 11 pages, 9 figures

Report number: LA-UR-17-23922

White Paper on Nuclear Data Needs and Capabilities for Basic Science

Authors: F. G. Kondev, M. Thoennessen, J. Batchelder, T. Kawano, J. Kelley, E. McCutchan, M. Smith, A. Sonzogni, I. Thompson

Abstract: The Workshop on Nuclear Data Needs and Capabilities for Basic Science was held at the University of Notre Dame on 10-11 August 2016. The purpose of this targeted workshop was to assemble and prioritize the needs of the nuclear physics research community for data sets, services and capabilities in areas including nuclear structure, nuclear reactions, nuclear astrophysics, fundamental interactions,… ▽ More The Workshop on Nuclear Data Needs and Capabilities for Basic Science was held at the University of Notre Dame on 10-11 August 2016. The purpose of this targeted workshop was to assemble and prioritize the needs of the nuclear physics research community for data sets, services and capabilities in areas including nuclear structure, nuclear reactions, nuclear astrophysics, fundamental interactions, neutrino physics and nuclear theory. An overview of nuclear data needs and capabilities identified at this meeting are summarized in the present document. △ Less

Submitted 12 May, 2017; originally announced May 2017.

Comments: This white paper and copies of the presentations can be found on the meeting website at http://meetings.nscl.msu.edu/2016ND_workshop/html/home.html

Calibration and performance studies of the balloon-borne hard X-ray polarimeter PoGO+

Authors: M. Chauvin, M. Friis, M. Jackson, T. Kawano, M. Kiss, V. Mikhalev, N. Ohashi, T. Stana, H. Takahashi, M. Pearce

Abstract: Polarimetric observations of celestial sources in the hard X-ray band stand to provide new information on emission mechanisms and source geometries. PoGO+ is a Compton scattering polarimeter (20-150 keV) optimised for the observation of the Crab (pulsar and wind nebula) and Cygnus X-1 (black hole binary), from a stratospheric balloon-borne platform launched from the Esrange Space Centre in summer… ▽ More Polarimetric observations of celestial sources in the hard X-ray band stand to provide new information on emission mechanisms and source geometries. PoGO+ is a Compton scattering polarimeter (20-150 keV) optimised for the observation of the Crab (pulsar and wind nebula) and Cygnus X-1 (black hole binary), from a stratospheric balloon-borne platform launched from the Esrange Space Centre in summer 2016. Prior to flight, the response of the polarimeter has been studied with polarised and unpolarised X-rays allowing a Geant4-based simulation model to be validated. The expected modulation factor for Crab observations is found to be $M_{\mathrm{Crab}}=(41.75\pm0.85)\%$, resulting in an expected Minimum Detectable Polarisation (MDP) of $7.3\%$ for a 7 day flight. This will allow a measurement of the Crab polarisation parameters with at least $5σ$ statistical significance assuming a polarisation fraction $\sim20\%$ $-$ a significant improvement over the PoGOLite Pathfinder mission which flew in 2013 and from which the PoGO+ design is developed. △ Less

Submitted 22 March, 2017; originally announced March 2017.

Comments: Accepted for publication in Nuclear Inst. and Methods in Physics Research, A

Black hole spin of Cyg X-1 determined from the softest state ever observed

Authors: Takafumi Kawano, Chris Done, Shin'ya Yamada, Hiromitsu Takahashi, Magnus Axelsson, Yasushi Fukazawa

Abstract: We show the softest ever spectrum from Cyg~X-1, detected in 2013 with Suzaku. This has the weakest high energy Compton tail ever seen from this object, so should give the cleanest view of the underlying disk spectrum, and hence the best determination of black hole spin from disk continuum fitting. Using the standard model of a disk with simple non-thermal Comptonisation to produce the weak high en… ▽ More We show the softest ever spectrum from Cyg~X-1, detected in 2013 with Suzaku. This has the weakest high energy Compton tail ever seen from this object, so should give the cleanest view of the underlying disk spectrum, and hence the best determination of black hole spin from disk continuum fitting. Using the standard model of a disk with simple non-thermal Comptonisation to produce the weak high energy tail gives a high spin black hole. However, we get a significantly better fit by including an additional, low temperature thermal Comptonisation component, which allows a much lower black hole spin. Corroboration of the existence of an additional Compton component comes from the frequency dependent hard lags seen in the rapid variability in archival high/soft state data. These can not be explained if the continuum is a single non-thermal Comptonisation component, but are instead consistent with a radially stratified, multi zone Comptonisation spectrum, where the spectrum is softer further from the black hole. A complex multi-zone Comptonisation continuum is required to explain both spectra and timing together, and this has an impact on the derived black hole spin. △ Less

Submitted 20 January, 2017; originally announced January 2017.

Comments: 15 pages, 5 figures

Exit-Channel Suppression in Statistical Reaction Theory

Authors: G. F. Bertsch, T. Kawano

Abstract: Statistical reaction theories such as Hauser-Feshbach assume that branching ratios follow Bohr's compound nucleus hypothesis by factorizing into independent probabilities for different channels. Corrections to the factorization hypothesis are known in both nuclear theory and quantum transport theory, particularly an enhanced memory of the entrance channel. We apply the Gaussian orthogonal ensemble… ▽ More Statistical reaction theories such as Hauser-Feshbach assume that branching ratios follow Bohr's compound nucleus hypothesis by factorizing into independent probabilities for different channels. Corrections to the factorization hypothesis are known in both nuclear theory and quantum transport theory, particularly an enhanced memory of the entrance channel. We apply the Gaussian orthogonal ensemble to study a complementary suppression of exit channel branching ratios. The effect can be large enough to provide information on the number of exit channels. The effect is demonstrated for the branching ratios of neutron-induced reactions on a $^{235}$U target. △ Less

Submitted 2 November, 2017; v1 submitted 1 January, 2017; originally announced January 2017.

Comments: 8 pages, 1 Figure, to be published in Physical Review Letters

Journal ref: Phys. Rev. Lett. 119, 222504 (2017)

Neutron-gamma competition for $β$-delayed neutron emission

Authors: Matthew Mumpower, Toshihiko Kawano, Peter Moller

Abstract: We present a coupled Quasi-particle Random Phase Approximation and Hauser-Feshbach (QRPA+HF) model for calculating delayed particle emission. This approach uses microscopic nuclear structure information which starts with Gamow-Teller strength distributions in the daughter nucleus, and then follows the statistical decay until the initial available excitation energy is exhausted. Explicitly included… ▽ More We present a coupled Quasi-particle Random Phase Approximation and Hauser-Feshbach (QRPA+HF) model for calculating delayed particle emission. This approach uses microscopic nuclear structure information which starts with Gamow-Teller strength distributions in the daughter nucleus, and then follows the statistical decay until the initial available excitation energy is exhausted. Explicitly included at each particle emission stage is $γ$-ray competition. We explore this model in the context of neutron emission of neutron-rich nuclei and find that neutron-gamma competition can lead to both increases and decreases in neutron emission probabilities, depending on the system considered. A second consequence of this formalism is a prediction of more neutrons on average being emitted after $β$-decay for nuclei near the neutron dripline compared to models that do not consider the statistical decay. △ Less

Submitted 5 August, 2016; originally announced August 2016.

Comments: 9 pages, 6 figures, submitted

Report number: LA-UR-16-25845

Late Time Emission of Prompt Fission Gamma Rays

Authors: P. Talou, T. Kawano, I. Stetcu, J. P. Lestone, E. McKigney, M. B. Chadwick

Abstract: The emission of prompt fission $γ$ rays within a few nanoseconds to a few microseconds following the scission point is studied in the Hauser-Feshbach formalism applied to the deexcitation of primary excited fission fragments. Neutron and $γ$-ray evaporations from fully accelerated fission fragments are calculated in competition at each stage of the decay, and the role of isomers in the fission pro… ▽ More The emission of prompt fission $γ$ rays within a few nanoseconds to a few microseconds following the scission point is studied in the Hauser-Feshbach formalism applied to the deexcitation of primary excited fission fragments. Neutron and $γ$-ray evaporations from fully accelerated fission fragments are calculated in competition at each stage of the decay, and the role of isomers in the fission products, before $β$-decay, is analyzed. The time evolution of the average total $γ$-ray energy, average total $γ$-ray multiplicity, and fragment-specific $γ$-ray spectra, is presented in the case of neutron-induced fission reactions of $^{235}$U and $^{239}$Pu, as well as spontaneous fission of $^{252}$Cf. The production of specific isomeric states is calculated and compared to available experimental data. About 7% of all prompt fission $γ$ rays are predicted to be emitted between 10 nsec and 5 $μ$sec following fission, in the case of $^{235}$U and $^{239}$Pu $(n_{\rm th},f)$ reactions, and up to 3% in the case of $^{252}$Cf spontaneous fission. The cumulative average total $γ$-ray energy increases by 2 to 5% in the same time interval. Finally, those results are shown to be robust against significant changes in the model input parameters. △ Less

Submitted 1 July, 2016; originally announced July 2016.

Comments: 11 pages, 14 figures

Report number: LA-UR-16-24045

Journal ref: Phys. Rev. C 94, 064613 (2016)

Optimising a balloon-borne polarimeter in the hard X-ray domain: from the PoGOLite Pathfinder to PoGO+

Authors: Maxime Chauvin, Miranda Jackson, Takafumi Kawano, Mózsi Kiss, Merlin Kole, Victor Mikhalev, Elena Moretti, Hiromitsu Takahashi, Mark Pearce

Abstract: PoGOLite is a balloon-borne hard X-ray polarimeter dedicated to the study of point sources. Compton scattered events are registered using an array of plastic scintillator units to determine the polarisation of incident X-rays in the energy range 20 - 240 keV. In 2013, a near circumpolar balloon flight of 14 days duration was completed after launch from Esrange, Sweden, resulting in a measurement o… ▽ More PoGOLite is a balloon-borne hard X-ray polarimeter dedicated to the study of point sources. Compton scattered events are registered using an array of plastic scintillator units to determine the polarisation of incident X-rays in the energy range 20 - 240 keV. In 2013, a near circumpolar balloon flight of 14 days duration was completed after launch from Esrange, Sweden, resulting in a measurement of the linear polarisation of the Crab emission. Building on the experience gained from this Pathfinder flight, the polarimeter is being modified to improve performance for a second flight in 2016. Such optimisations, based on Geant4 Monte Carlo simulations, take into account the source characteristics, the instrument response and the background environment which is dominated by atmospheric neutrons. This paper describes the optimisation of the polarimeter and details the associated increase in performance. The resulting design, PoGO+, is expected to improve the Minimum Detectable Polarisation (MDP) for the Crab from 19.8% to 11.1% for a 5 day flight. Assuming the same Crab polarisation fraction as measured during the 2013 flight, this improvement in MDP will allow a 5σ constrained result. It will also allow the study of the nebula emission only (Crab off-pulse) and Cygnus X-1 if in the hard state. △ Less

Submitted 14 June, 2016; originally announced June 2016.

Statistical Hauser-Feshbach theory with width fluctuation correction including direct reaction channels for neutron induced reaction at low energies

Authors: T. Kawano, R. Capote, S. Hilaire

Abstract: A model to calculate particle-induced reaction cross sections with statistical Hauser-Feshbach theory including direct reactions is given. The energy average of scattering matrix from the coupled-channels optical model is diagonalized by the transformation proposed by Engelbrecht and Weidenmüller. The ensemble average of $S$-matrix elements in the diagonalized channel space is approximated by a mo… ▽ More A model to calculate particle-induced reaction cross sections with statistical Hauser-Feshbach theory including direct reactions is given. The energy average of scattering matrix from the coupled-channels optical model is diagonalized by the transformation proposed by Engelbrecht and Weidenmüller. The ensemble average of $S$-matrix elements in the diagonalized channel space is approximated by a model of Moldauer [Phys.Rev.C {\bf 12}, 744 (1975)] using newly parametrized channel degree-of-freedom $ν_a$ to better describe the Gaussian Orthogonal Ensemble (GOE) reference calculations. Moldauer approximation is confirmed by a Monte Carlo study using randomly generated $S$-matrix, as well as the GOE three-fold integration formula. The method proposed is applied to the $^{238}$U(n,n') cross section calculation in the fast energy range, showing an enhancement in the inelastic scattering cross sections. △ Less

Submitted 28 April, 2016; originally announced April 2016.

Observation of polarised hard X-ray emission from the Crab by the PoGOLite Pathfinder

Authors: M. Chauvin, H. -G. Florén, M. Jackson, T. Kamae, T. Kawano, M. Kiss, M. Kole, V. Mikhalev, E. Moretti, G. Olofsson, S. Rydström, H. Takahashi, A. Iyudin, M. Arimoto, Y. Fukazawa, J. Kataoka, N. Kawai, T. Mizuno, F. Ryde, H. Tajima, T. Takahashi, M. Pearce

Abstract: We have measured the linear polarisation of hard X-ray emission from the Crab in a previously unexplored energy interval, 20-120 keV. The introduction of two new observational parameters, the polarisation fraction and angle stands to disentangle geometrical and physical effects, thereby providing information on the pulsar wind geometry and magnetic field environment. Measurements are conducted usi… ▽ More We have measured the linear polarisation of hard X-ray emission from the Crab in a previously unexplored energy interval, 20-120 keV. The introduction of two new observational parameters, the polarisation fraction and angle stands to disentangle geometrical and physical effects, thereby providing information on the pulsar wind geometry and magnetic field environment. Measurements are conducted using the PoGOLite Pathfinder - a balloon-borne polarimeter. Polarisation is determined by measuring the azimuthal Compton scattering angle of incident X-rays in an array of plastic scintillators housed in an anticoincidence well. The polarimetric response has been characterised prior to flight using both polarised and unpolarised calibration sources. We address possible systematic effects through observations of a background field. The measured polarisation fraction for the integrated Crab light-curve is ($18.4^{+9.8}_{-10.6}$)%, corresponding to an upper limit (99% credibility) of 42.4%, for a polarisation angle of ($149.2\pm16.0)^\circ$. △ Less

Submitted 30 November, 2015; v1 submitted 9 November, 2015; originally announced November 2015.

Comments: 6 pages, 5 figures. Accepted for publication in MNRAS Letters (In this revision the polarisation angle is correctly stated as 149.2$\pm$16.0 degrees and not 142.2$\pm$16.0 degrees (typo corrected in abstract and on page 4)

The first demonstration of the concept of "narrow-FOV Si/CdTe semiconductor Compton camera"

Authors: Yuto Ichinohe, Yuusuke Uchida, Shin Watanabe, Ikumi Edahiro, Katsuhiro Hayashi, Takafumi Kawano, Masanori Ohno, Masayuki Ohta, Shin'ichiro Takeda, Yasushi Fukazawa, Miho Katsuragawa, Kazuhiro Nakazawa, Hirokazu Odaka, Hiroyasu Tajima, Hiromitsu Takahashi, Tadayuki Takahashi, Takayuki Yuasa

Abstract: The Soft Gamma-ray Detector (SGD), to be deployed onboard the {\it ASTRO-H} satellite, has been developed to provide the highest sensitivity observations of celestial sources in the energy band of 60-600~keV by employing a detector concept which uses a Compton camera whose field-of-view is restricted by a BGO shield to a few degree (narrow-FOV Compton camera). In this concept, the background from… ▽ More The Soft Gamma-ray Detector (SGD), to be deployed onboard the {\it ASTRO-H} satellite, has been developed to provide the highest sensitivity observations of celestial sources in the energy band of 60-600~keV by employing a detector concept which uses a Compton camera whose field-of-view is restricted by a BGO shield to a few degree (narrow-FOV Compton camera). In this concept, the background from outside the FOV can be heavily suppressed by constraining the incident direction of the gamma ray reconstructed by the Compton camera to be consistent with the narrow FOV. We, for the first time, demonstrate the validity of the concept using background data taken during the thermal vacuum test and the low-temperature environment test of the flight model of SGD on ground. We show that the measured background level is suppressed to less than 10\% by combining the event rejection using the anti-coincidence trigger of the active BGO shield and by using Compton event reconstruction techniques. More than 75\% of the signals from the field-of-view are retained against the background rejection, which clearly demonstrates the improvement of signal-to-noise ratio. The estimated effective area of 22.8~cm$^2$ meets the mission requirement even though not all of the operational parameters of the instrument have been fully optimized yet. △ Less

Submitted 25 September, 2015; originally announced September 2015.

Comments: 40 pages, 12 figures, accepted for publication in Nuclear Inst. and Methods in Physics Research, A