-
Isochronous mass spectrometry at the RIKEN Rare-RI Ring facility
Authors:
D. Nagae,
S. Omika,
Y. Abe,
Y. Yamaguchi,
F. Suzaki,
K. Wakayama,
N. Tadano,
R. Igosawa,
K. Inomata,
H. Arakawa,
K. Nishimuro,
T. Fujii,
T. Mitsui,
T. Yamaguchi,
T. Suzuki,
S. Suzuki,
T. Moriguchi,
M. Amano,
D. Kamioka,
A. Ozawa,
S. Naimi,
Z. Ge,
Y. Yanagisawa,
H. Baba,
S. Michimasa
, et al. (6 additional authors not shown)
Abstract:
A dedicated isochronous storage ring, named the Rare-RI Ring, was constructed at the RI Beam Factory of RIKEN, aiming at precision mass measurements of nuclei located in uncharted territories of the nuclear chart. The Rare-RI Ring employs the isochronous mass spectrometry technique with the goal to achieve a relative mass precision of $10^{-6}$ within a measurement time of less than 1 ms. The perf…
▽ More
A dedicated isochronous storage ring, named the Rare-RI Ring, was constructed at the RI Beam Factory of RIKEN, aiming at precision mass measurements of nuclei located in uncharted territories of the nuclear chart. The Rare-RI Ring employs the isochronous mass spectrometry technique with the goal to achieve a relative mass precision of $10^{-6}$ within a measurement time of less than 1 ms. The performance of the facility was demonstrated through mass measurements of neutron-rich nuclei with well-known masses. Velocity or magnetic rigidity is measured for every particle prior to its injection into the ring, wherein its revolution time is accurately determined. The latter quantity is used to determine the mass of the particle, while the former one is needed for non-isochronicity corrections. Mass precisions on the order of $10^{-5}$ were achieved in the first commissioning, which demonstrates that Rare-RI Ring is a powerful tool for mass spectrometry of short-lived nuclei.
△ Less
Submitted 8 July, 2024;
originally announced July 2024.
-
Spectroscopy of deeply bound orbitals in neutron-rich Ca isotopes
Authors:
P. J. Li,
J. Lee,
P. Doornenbal,
S. Chen,
S. Wang,
A. Obertelli,
Y. Chazono,
J. D. Holt,
B. S. Hu,
K. Ogata,
Y. Utsuno,
K. Yoshida,
N. L. Achouri,
H. Baba,
F. Browne,
D. Calvet,
F. Château,
N. Chiga,
A. Corsi,
M. L. Cortés,
A. Delbart,
J-M. Gheller,
A. Giganon,
A. Gillibert,
C. Hilaire
, et al. (63 additional authors not shown)
Abstract:
The calcium isotopes are an ideal system to investigate the evolution of shell structure and magic numbers. Although the properties of surface nucleons in calcium have been well studied, probing the structure of deeply bound nucleons remains a challenge. Here, we report on the first measurement of unbound states in $^{53}$Ca and $^{55}$Ca, populated from \ts{54,56}Ca($p,pn$) reactions at a beam en…
▽ More
The calcium isotopes are an ideal system to investigate the evolution of shell structure and magic numbers. Although the properties of surface nucleons in calcium have been well studied, probing the structure of deeply bound nucleons remains a challenge. Here, we report on the first measurement of unbound states in $^{53}$Ca and $^{55}$Ca, populated from \ts{54,56}Ca($p,pn$) reactions at a beam energy of around 216 MeV/nucleon at the RIKEN Radioactive Isotopes Beam Factory. The resonance properties, partial cross sections, and momentum distributions of these unbound states were analyzed. Orbital angular momentum $l$ assignments were extracted from momentum distributions based on calculations using the distorted wave impulse approximation (DWIA) reaction model. The resonances at excitation energies of 5516(41)\,keV in $^{53}$Ca and 6000(250)\,keV in $^{55}$Ca indicate a significant $l$\, =\,3 component, providing the first experimental evidence for the $ν0f_{7/2}$ single-particle strength of unbound hole states in the neutron-rich Ca isotopes. The observed excitation energies and cross-sections point towards extremely localized and well separated strength distributions, with some fragmentation for the $��0f_{7/2}$ orbital in $^{55}$Ca. These results are in good agreement with predictions from shell-model calculations using the effective GXPF1Bs interaction and \textit{ab initio} calculations and diverge markedly from the experimental distributions in the nickel isotones at $Z=28$.
△ Less
Submitted 5 July, 2024;
originally announced July 2024.
-
Skeletal Kinetics Reduction for Astrophysical Reaction Networks
Authors:
A. G. Nouri,
Y. Liu,
P. Givi,
H. Babaee,
D. Livescu
Abstract:
A novel methodology is developed to extract accurate skeletal reaction models for nuclear combustion. Local sensitivities of isotope mass fractions with respect to reaction rates are modeled based on the forced optimally time-dependent (f-OTD) scheme. These sensitivities are then analyzed temporally to generate skeletal models. The methodology is demonstrated by conducting skeletal reduction of co…
▽ More
A novel methodology is developed to extract accurate skeletal reaction models for nuclear combustion. Local sensitivities of isotope mass fractions with respect to reaction rates are modeled based on the forced optimally time-dependent (f-OTD) scheme. These sensitivities are then analyzed temporally to generate skeletal models. The methodology is demonstrated by conducting skeletal reduction of constant density and temperature burning of carbon and oxygen relevant to SNe Ia. The 495-isotopes Torch model is chosen as the detailed reaction network. A map of maximum production of $^{56}\text{Ni}$ in SNe Ia is produced for different temperatures, densities, and proton to neutron ratios. The f-OTD simulations and the sensitivity analyses are then performed with initial conditions from this map. A series of skeletal models are derived and their performances are assessed by comparison against currently existing skeletal models. Previous models have been constructed intuitively by assuming the dominance of $α$-chain reactions. The comparison of the newly generated skeletal models against previous models is based on the predicted energy release and $^{44}\text{Ti}$ and $^{56}\text{Ni}$ abundances by each model. The consequences of $\mathtt{y}_e \neq 0.5$ in the initial composition are also explored where $\mathtt{y}_e$ is the electron fraction. The simulated results show that $^{56}\text{Ni}$ production decreases by decreasing $\mathtt{y}_e$ as expected, and that the $^{43}\text{Sc}$ is a key isotope in proton and neutron channels toward $^{56}\text{Ni}$ production. It is shown that an f-OTD skeletal model with 150 isotopes can accurately predict the $^{56}\text{Ni}$ abundance in SNe Ia for $\mathtt{y}_e \lesssim 0.5$ initial conditions.
△ Less
Submitted 26 April, 2024;
originally announced April 2024.
-
Cross Interpolation for Solving High-Dimensional Dynamical Systems on Low-Rank Tucker and Tensor Train Manifolds
Authors:
Behzad Ghahremani,
Hessam Babaee
Abstract:
We present a novel tensor interpolation algorithm for the time integration of nonlinear tensor differential equations (TDEs) on the tensor train and Tucker tensor low-rank manifolds, which are the building blocks of many tensor network decompositions. This paper builds upon our previous work (Donello et al., Proceedings of the Royal Society A, Vol. 479, 2023) on solving nonlinear matrix differenti…
▽ More
We present a novel tensor interpolation algorithm for the time integration of nonlinear tensor differential equations (TDEs) on the tensor train and Tucker tensor low-rank manifolds, which are the building blocks of many tensor network decompositions. This paper builds upon our previous work (Donello et al., Proceedings of the Royal Society A, Vol. 479, 2023) on solving nonlinear matrix differential equations on low-rank matrix manifolds using CUR decompositions. The methodology we present offers multiple advantages: (i) It delivers near-optimal computational savings both in terms of memory and floating-point operations by leveraging cross algorithms based on the discrete empirical interpolation method to strategically sample sparse entries of the time-discrete TDEs to advance the solution in low-rank form. (ii) Numerical demonstrations show that the time integration is robust in the presence of small singular values. (iii) High-order explicit Runge-Kutta time integration schemes are developed. (iv) The algorithm is easy to implement, as it requires the evaluation of the full-order model at strategically selected entries and does not use tangent space projections, whose efficient implementation is intrusive. We demonstrate the efficiency of the presented algorithm for several test cases, including a nonlinear 100-dimensional TDE for the evolution of a tensor of size $70^{100} \approx 3.2 \times 10^{184}$ and a stochastic advection-diffusion-reaction equation with a tensor of size $4.7 \times 10^9$.
△ Less
Submitted 10 June, 2024; v1 submitted 19 March, 2024;
originally announced March 2024.
-
Reconstructing Blood Flow in Data-Poor Regimes: A Vasculature Network Kernel for Gaussian Process Regression
Authors:
Shaghayegh Z. Ashtiani,
Mohammad Sarabian,
Kaveh Laksari,
Hessam Babaee
Abstract:
Blood flow reconstruction in the vasculature is important for many clinical applications. However, in clinical settings, the available data are often quite limited. For instance, Transcranial Doppler ultrasound (TCD) is a noninvasive clinical tool that is commonly used in the clinical settings to measure blood velocity waveform at several locations on brain's vasculature. This amount of data is gr…
▽ More
Blood flow reconstruction in the vasculature is important for many clinical applications. However, in clinical settings, the available data are often quite limited. For instance, Transcranial Doppler ultrasound (TCD) is a noninvasive clinical tool that is commonly used in the clinical settings to measure blood velocity waveform at several locations on brain's vasculature. This amount of data is grossly insufficient for training machine learning surrogate models, such as deep neural networks or Gaussian process regression. In this work, we propose a Gaussian process regression approach based on physics-informed kernels, enabling near-real-time reconstruction of blood flow in data-poor regimes. We introduce a novel methodology to reconstruct the kernel within the vascular network, which is a non-Euclidean space. The proposed kernel encodes both spatiotemporal and vessel-to-vessel correlations, thus enabling blood flow reconstruction in vessels that lack direct measurements. We demonstrate that any prediction made with the proposed kernel satisfies the conservation of mass principle. The kernel is constructed by running stochastic one-dimensional blood flow simulations, where the stochasticity captures the epistemic uncertainties, such as lack of knowledge about boundary conditions and uncertainties in vasculature geometries. We demonstrate the performance of the model on three test cases, namely, a simple Y-shaped bifurcation, abdominal aorta, and the Circle of Willis in the brain.
△ Less
Submitted 14 March, 2024;
originally announced March 2024.
-
Mass, spectroscopy and two-neutron decay of $^{16}$Be
Authors:
B. Monteagudo,
F. M. Marqués,
J. Gibelin,
N. A. Orr,
A. Corsi,
Y. Kubota,
J. Casal,
J. Gómez-Camacho,
G. Authelet,
H. Baba,
C. Caesar,
D. Calvet,
A. Delbart,
M. Dozono,
J. Feng,
F. Flavigny,
J. -M. Gheller,
A. Giganon,
A. Gillibert,
K. Hasegawa,
T. Isobe,
Y. Kanaya,
S. Kawakami,
D. Kim,
Y. Kiyokawa
, et al. (43 additional authors not shown)
Abstract:
The structure and decay of the most neutron-rich beryllium isotope, $^{16}$Be, has been investigated following proton knockout from a high-energy $^{17}$B beam. Two relatively narrow resonances were observed for the first time, with energies of $0.84(3)$ and $2.15(5)$ MeV above the two-neutron decay threshold and widths of $0.32(8)$ and $0.95(15)$ MeV respectively. These were assigned to be the gr…
▽ More
The structure and decay of the most neutron-rich beryllium isotope, $^{16}$Be, has been investigated following proton knockout from a high-energy $^{17}$B beam. Two relatively narrow resonances were observed for the first time, with energies of $0.84(3)$ and $2.15(5)$ MeV above the two-neutron decay threshold and widths of $0.32(8)$ and $0.95(15)$ MeV respectively. These were assigned to be the ground ($J^π=0^+$) and first excited ($2^+$) state, with $E_x=1.31(6)$ MeV. The mass excess of $^{16}$Be was thus deduced to be $56.93(13)$ MeV, some $0.5$ MeV more bound than the only previous measurement. Both states were observed to decay by direct two-neutron emission. Calculations incorporating the evolution of the wavefunction during the decay as a genuine three-body process reproduced the principal characteristics of the neutron-neutron energy spectra for both levels, indicating that the ground state exhibits a strong spatially compact dineutron component, while the 2$^+$ level presents a far more diffuse neutron-neutron distribution.
△ Less
Submitted 30 January, 2024;
originally announced January 2024.
-
First Exploration of Monopole-Driven Shell Evolution above the N = 126 shell closure: new Millisecond Isomers in 213Tl and 215Tl
Authors:
T. T. Yeung,
A. I. Morales,
J. Wu,
M. Liu,
C. Yuan,
S. Nishimura,
V. H. Phong,
N. Fukuda,
J. L. Tain,
T. Davinson,
K. P. Rykaczewski,
R. Yokoyama,
T. Isobe,
M. Niikura,
Zs. Podolyak,
G. Alcala,
A. Algora,
J. Agramunt,
C. Appleton,
H. Baba,
R. Caballero-Folch,
P. Calvino,
M. P. Carpenter,
I. Dillmann,
A. Estrade
, et al. (30 additional authors not shown)
Abstract:
Isomer spectroscopy of heavy neutron-rich nuclei beyond the N=126 closed shell has been performed for the first time at the Radioactive Isotope Beam Factory of the RIKEN Nishina Center. New millisecond isomers have been identified at low excitation energies, 985.3(19) keV in 213Tl and 874(5) keV in 215Tl. The measured half-lives of 1.34(5) ms in 213Tl and 3.0(3) ms in 215Tl suggest spins and parit…
▽ More
Isomer spectroscopy of heavy neutron-rich nuclei beyond the N=126 closed shell has been performed for the first time at the Radioactive Isotope Beam Factory of the RIKEN Nishina Center. New millisecond isomers have been identified at low excitation energies, 985.3(19) keV in 213Tl and 874(5) keV in 215Tl. The measured half-lives of 1.34(5) ms in 213Tl and 3.0(3) ms in 215Tl suggest spins and parities 11/2- with the single proton-hole configuration h11/2 as leading component. They are populated via E1 transitions by the decay of higher-lying isomeric states with proposed spin and parity 17/2+, interpreted as arising from a single s1/2 proton hole coupled to the 8+ seniority isomer in the (A+1)Pb cores. The lowering of the 11/2- states is ascribed to an increase of the h11/2 proton effective single-particle energy as the second g9/2 orbital is filled by neutrons, owing to a significant reduction of the proton-neutron monopole interaction between the h11/2 and g9/2 orbitals. The new ms-isomers provide the first experimental observation of shell evolution in the almost unexplored N>126 nuclear region below doubly-magic 208Pb.
△ Less
Submitted 25 April, 2024; v1 submitted 12 January, 2024;
originally announced January 2024.
-
A DEIM Tucker Tensor Cross Algorithm and its Application to Dynamical Low-Rank Approximation
Authors:
Behzad Ghahremani,
Hessam Babaee
Abstract:
We introduce a Tucker tensor cross approximation method that constructs a low-rank representation of a $d$-dimensional tensor by sparsely sampling its fibers. These fibers are selected using the discrete empirical interpolation method (DEIM). Our proposed algorithm is referred to as DEIM fiber sampling (DEIM-FS). For a rank-$r$ approximation of an $\mathcal{O}(N^d)$ tensor, DEIM-FS requires access…
▽ More
We introduce a Tucker tensor cross approximation method that constructs a low-rank representation of a $d$-dimensional tensor by sparsely sampling its fibers. These fibers are selected using the discrete empirical interpolation method (DEIM). Our proposed algorithm is referred to as DEIM fiber sampling (DEIM-FS). For a rank-$r$ approximation of an $\mathcal{O}(N^d)$ tensor, DEIM-FS requires access to only $dNr^{d-1}$ tensor entries, a requirement that scales linearly with the tensor size along each mode. We demonstrate that DEIM-FS achieves an approximation accuracy close to the Tucker-tensor approximation obtained via higher-order singular value decomposition at a significantly reduced cost. We also present DEIM-FS (iterative) that does not require access to singular vectors of the target tensor unfolding and can be viewed as a black-box Tucker tensor algorithm. We employ DEIM-FS to reduce the computational cost associated with solving nonlinear tensor differential equations (TDEs) using dynamical low-rank approximation (DLRA). The computational cost of solving DLRA equations can become prohibitive when the exact rank of the right-hand side tensor is large. This issue arises in many TDEs, especially in cases involving non-polynomial nonlinearities, where the right-hand side tensor has full rank. This necessitates the storage and computation of tensors of size $\mathcal{O}(N^d)$. We show that DEIM-FS results in significant computational savings for DLRA by constructing a low-rank Tucker approximation of the right-hand side tensor on the fly. Another advantage of using DEIM-FS is to significantly simplify the implementation of DLRA equations, irrespective of the type of TDEs. We demonstrate the efficiency of the algorithm through several examples including solving high-dimensional partial differential equations.
△ Less
Submitted 8 January, 2024;
originally announced January 2024.
-
Future-proofing geotechnics workflows: accelerating problem-solving with large language models
Authors:
Stephen Wu,
Yu Otake,
Daijiro Mizutani,
Chang Liu,
Kotaro Asano,
Nana Sato,
Hidetoshi Baba,
Yusuke Fukunaga,
Yosuke Higo,
Akiyoshi Kamura,
Shinnosuke Kodama,
Masataka Metoki,
Tomoka Nakamura,
Yuto Nakazato,
Taiga Saito,
Akihiro Shioi,
Masahiro Takenobu,
Keigo Tsukioka,
Ryo Yoshikawa
Abstract:
The integration of Large Language Models (LLMs) like ChatGPT into the workflows of geotechnical engineering has a high potential to transform how the discipline approaches problem-solving and decision-making. This paper delves into the innovative application of LLMs in geotechnical engineering, as explored in a hands-on workshop held in Tokyo, Japan. The event brought together a diverse group of 2…
▽ More
The integration of Large Language Models (LLMs) like ChatGPT into the workflows of geotechnical engineering has a high potential to transform how the discipline approaches problem-solving and decision-making. This paper delves into the innovative application of LLMs in geotechnical engineering, as explored in a hands-on workshop held in Tokyo, Japan. The event brought together a diverse group of 20 participants, including students, researchers, and professionals from academia, industry, and government sectors, to investigate practical uses of LLMs in addressing specific geotechnical challenges. The workshop facilitated the creation of solutions for four different practical geotechnical problems as illustrative examples, culminating in the development of an academic paper. The paper discusses the potential of LLMs to transform geotechnical engineering practices, highlighting their proficiency in handling a range of tasks from basic data analysis to complex, multimodal problem-solving. It also addresses the challenges in implementing LLMs, particularly in achieving high precision and accuracy in specialized tasks, and underscores the need for expert oversight. The findings demonstrate LLMs' effectiveness in enhancing efficiency, data processing, and decision-making in geotechnical engineering, suggesting a paradigm shift towards more integrated, data-driven approaches in this field. This study not only showcases the potential of LLMs in a specific engineering domain, but also sets a precedent for their broader application in interdisciplinary research and practice, where the synergy of human expertise and artificial intelligence redefines the boundaries of problem-solving.
△ Less
Submitted 14 December, 2023;
originally announced December 2023.
-
Constraining nucleon effective masses with flow and stopping observables from the S$π$RIT experiment
Authors:
C. Y. Tsang,
M. Kurata-Nishimura,
M. B. Tsang,
W. G. Lynch,
Y. X. Zhang,
J. Barney,
J. Estee,
G. Jhang,
R. Wang,
M. Kaneko,
J. W. Lee,
T. Isobe,
T. Murakami,
D. S. Ahn,
L. Atar,
T. Aumann,
H. Baba,
K. Boretzky,
J. Brzychczyk,
G. Cerizza,
N. Chiga,
N. Fukuda,
I. Gasparic,
B. Hong,
A. Horvat
, et al. (30 additional authors not shown)
Abstract:
Properties of the nuclear equation of state (EoS) can be probed by measuring the dynamical properties of nucleus-nucleus collisions. In this study, we present the directed flow ($v_1$), elliptic flow ($v_2$) and stopping (VarXZ) measured in fixed target Sn + Sn collisions at 270 AMeV with the S$π$RIT Time Projection Chamber. We perform Bayesian analyses in which EoS parameters are varied simultane…
▽ More
Properties of the nuclear equation of state (EoS) can be probed by measuring the dynamical properties of nucleus-nucleus collisions. In this study, we present the directed flow ($v_1$), elliptic flow ($v_2$) and stopping (VarXZ) measured in fixed target Sn + Sn collisions at 270 AMeV with the S$π$RIT Time Projection Chamber. We perform Bayesian analyses in which EoS parameters are varied simultaneously within the Improved Quantum Molecular Dynamics-Skyrme (ImQMD-Sky) transport code to obtain a multivariate correlated constraint. The varied parameters include symmetry energy, $S_0$, and slope of the symmetry energy, $L$, at saturation density, isoscalar effective mass, $m_{s}^*/m_{N}$, isovector effective mass, $m_{v}^{*}/m_{N}$ and the in-medium cross-section enhancement factor $η$. We find that the flow and VarXZ observables are sensitive to the splitting of proton and neutron effective masses and the in-medium cross-section. Comparisons of ImQMD-Sky predictions to the S$π$RIT data suggest a narrow range of preferred values for $m_{s}^*/m_{N}$, $m_{v}^{*}/m_{N}$ and $η$.
△ Less
Submitted 8 December, 2023;
originally announced December 2023.
-
Low-Rank Solution Operator for Forced Linearized Dynamics with Unsteady Base Flows
Authors:
Alireza Amiri-Margavi,
Hessam Babaee
Abstract:
Understanding the linear growth of disturbances due to external forcing is crucial for flow stability analysis, flow control, and uncertainty quantification. These applications typically require a large number of forward simulations of the forced linearized dynamics, often in a brute-force fashion. When dealing with simple steady-state or periodic base flows, there exist powerful and cost-effectiv…
▽ More
Understanding the linear growth of disturbances due to external forcing is crucial for flow stability analysis, flow control, and uncertainty quantification. These applications typically require a large number of forward simulations of the forced linearized dynamics, often in a brute-force fashion. When dealing with simple steady-state or periodic base flows, there exist powerful and cost-effective solution operator techniques. Once constructed, these operators can be used to determine the response to various forcings with negligible computational cost. However, these methods are not applicable to problems with arbitrarily time-dependent base flows. This paper develops and investigates reduced-order modeling with time-dependent bases (TDBs) to build low-rank solution operators for forced linearized dynamics with arbitrarily time-dependent base flows. In particular, we use forced optimally time-dependent decomposition (f-OTD), which extracts the time-dependent correlated structures of the flow response to various excitations.
We also demonstrate that in the case of a steady-state mean flow subject to harmonic forcing, the f-OTD subspace converges to the dominant resolvent analysis modes. The demonstration includes four cases: a toy model, the Burgers equation, the 2D temporally evolving jet, and two-dimensional decaying isotropic turbulence. In these cases, we demonstrate the utility of the low-rank operator for (i) identifying the excitation that leads to maximum amplification, and (ii) reconstructing the full-state flow without incurring additional cost.
△ Less
Submitted 10 December, 2023; v1 submitted 10 November, 2023;
originally announced December 2023.
-
Validation of the $^{10}\mathrm{Be}$ Ground-State Molecular Structure Using $^{10}\mathrm{Be}(p,pα)^{6}\mathrm{He}$ Triple Differential Reaction Cross-Section Measurements
Authors:
P. J. Li,
D. Beaumel,
J. Lee,
M. Assié,
S. Chen,
S. Franchoo,
J. Gibelin,
F. Hammache,
T. Harada,
Y. Kanada-En'yo,
Y. Kubota,
S. Leblond,
P. F. Liang,
T. Lokotko,
M. Lyu,
F. M. Marqués,
Y. Matsuda,
K. Ogata,
H. Otsu,
E. Rindel,
L. Stuhl,
D. Suzuki,
Y. Togano,
T. Tomai,
X. X. Xu
, et al. (36 additional authors not shown)
Abstract:
The cluster structure of the neutron-rich isotope $^{10}$Be has been probed via the $(p,pα)$ reaction at 150 MeV/nucleon in inverse kinematics and in quasifree conditions. The populated states of $^{6}$He residues were investigated through missing mass spectroscopy. The triple differential cross-section for the ground-state transition was extracted for quasifree angle pairs ($θ_{p}$, $θ_α$) and co…
▽ More
The cluster structure of the neutron-rich isotope $^{10}$Be has been probed via the $(p,pα)$ reaction at 150 MeV/nucleon in inverse kinematics and in quasifree conditions. The populated states of $^{6}$He residues were investigated through missing mass spectroscopy. The triple differential cross-section for the ground-state transition was extracted for quasifree angle pairs ($θ_{p}$, $θ_α$) and compared to distorted-wave impulse approximation reaction calculations performed in a microscopic framework using successively the Tohsaki-Horiuchi-Schuck-Röpke product wave-function and the wave-function deduced from Antisymmetrized Molecular Dynamics calculations. The remarkable agreement between calculated and measured cross-sections in both shape and magnitude validates the molecular structure description of the $^{10}$Be ground-state, configured as an $α$-$α$ core with two valence neutrons occupying $π$-type molecular orbitals.
△ Less
Submitted 21 November, 2023;
originally announced November 2023.
-
Level Structures of $^{56,58}$Ca Cast Doubt on a doubly magic $^{60}$Ca
Authors:
S. Chen,
F. Browne,
P. Doornenbal,
J. Lee,
A. Obertelli,
Y. Tsunoda,
T. Otsuka,
Y. Chazono,
G. Hagen,
J. D. Holt,
G. R. Jansen,
K. Ogata,
N. Shimizu,
Y. Utsuno,
K. Yoshida,
N. L. Achouri,
H. Baba,
D. Calvet,
F. Château,
N. Chiga,
A. Corsi,
M. L. Cortés,
A. Delbart,
J. -M. Gheller,
A. Giganon
, et al. (58 additional authors not shown)
Abstract:
Gamma decays were observed in $^{56}$Ca and $^{58}$Ca following quasi-free one-proton knockout reactions from $^{57,59}$Sc beams at $\approx 200$ MeV/nucleon. For $^{56}$Ca, a $γ$ ray transition was measured to be 1456(12) keV, while for $^{58}$Ca an indication for a transition was observed at 1115(34) keV. Both transitions were tentatively assigned as the $2^+_1 \rightarrow 0^+_{gs}$ decays, and…
▽ More
Gamma decays were observed in $^{56}$Ca and $^{58}$Ca following quasi-free one-proton knockout reactions from $^{57,59}$Sc beams at $\approx 200$ MeV/nucleon. For $^{56}$Ca, a $γ$ ray transition was measured to be 1456(12) keV, while for $^{58}$Ca an indication for a transition was observed at 1115(34) keV. Both transitions were tentatively assigned as the $2^+_1 \rightarrow 0^+_{gs}$ decays, and were compared to results from ab initio and conventional shell-model approaches. A shell-model calculation in a wide model space with a marginally modified effective nucleon-nucleon interaction depicts excellent agreement with experiment for $2^+_1$ level energies, two-neutron separation energies, and reaction cross sections, corroborating the formation of a new nuclear shell above the $N$ = 34 shell. Its constituents, the $0f_{5/2}$ and $0g_{9/2}$ orbitals, are almost degenerate. This degeneracy precludes the possibility for a doubly magic $^{60}$Ca and potentially drives the dripline of Ca isotopes to $^{70}$Ca or even beyond.
△ Less
Submitted 13 July, 2023;
originally announced July 2023.
-
Silicon tracker array for RIB experiments at SAMURAI
Authors:
A. I. Stefanescu,
V. Panin,
L. Trache,
T. Motobayashi,
H. Otsu,
A. Saastamoinen,
T. Uesaka,
L. Stuhl,
J. Tanaka,
D. Tudor,
I. C. Stefanescu,
A. E. Spiridon,
K. Yoneda,
H. Baba,
M. Kurokawa,
Y. Togano,
Z. Halasz,
M. Sasano,
S. Ota,
Y. Kubota,
D. S. Ahn,
T. Kobayashi,
Z. Elekes,
N. Fukuda,
H. Takeda
, et al. (27 additional authors not shown)
Abstract:
This work describes a silicon tracker system developed for experiments with proton-rich radioactive ion beams at the SAMURAI superconducting spectrometer of RIBF at RIKEN. The system is designed for accurate angular reconstruction and atomic number identification of relativistic heavy ions and protons which are simultaneously produced in reactions motivated by studies of proton capture reactions o…
▽ More
This work describes a silicon tracker system developed for experiments with proton-rich radioactive ion beams at the SAMURAI superconducting spectrometer of RIBF at RIKEN. The system is designed for accurate angular reconstruction and atomic number identification of relativistic heavy ions and protons which are simultaneously produced in reactions motivated by studies of proton capture reactions of interest for nuclear astrophysics. The technical characteristics of the tracking array are described in detail as are its performance in two pilot experiments. The physics justification for such a system is also presented.
△ Less
Submitted 13 July, 2023;
originally announced July 2023.
-
Searching for universality of dineutron correlation at the surface of Borromean nuclei
Authors:
A. Corsi,
Y. Kubota,
J. Casal,
M. Gomez-Ramos,
A. M. Moro,
G. Authelet,
H. Baba,
C. Caesar,
D. Calvet,
A. Delbart,
M. Dozono,
J. Feng,
F. Flavigny,
J. -M. Gheller,
J. Gibelin,
A. Giganon,
A. Gillibert,
K. Hasegawa,
T. Isobe,
Y. Kanaya,
S. Kawakami,
D. Kim,
Y. Kiyokawa,
M. Kobayashi,
N. Kobayashi
, et al. (43 additional authors not shown)
Abstract:
The dineutron correlation is systematically studied in three different Borromean nuclei near the neutron dripline, 11Li, 14Be and 17B, via the (p, pn) knockout reaction measured at the RIBF facility in RIKEN. For the three nuclei, the correlation angle between the valence neutrons is found to be largest in the same range of intrinsic momenta, which can be associated to the nuclear surface. This re…
▽ More
The dineutron correlation is systematically studied in three different Borromean nuclei near the neutron dripline, 11Li, 14Be and 17B, via the (p, pn) knockout reaction measured at the RIBF facility in RIKEN. For the three nuclei, the correlation angle between the valence neutrons is found to be largest in the same range of intrinsic momenta, which can be associated to the nuclear surface. This result reinforces the prediction that the formation of the dineutron is universal in environments with low neutron density, such as the surface of neutron-rich Borromean nuclei.
△ Less
Submitted 12 July, 2023;
originally announced July 2023.
-
Decay Pattern of Pygmy States Observed in Neutron-Rich 26 Ne
Authors:
J. Gibelin,
D. Beaumel,
T. Motobayashi,
Y. Blumenfeld,
N. Aoi,
H. Baba,
Z. Elekes,
S. Fortier,
N. Frascaria,
N. Fukuda,
T. Gomi,
K. Ishikawa,
Y. Kondo,
T. Kubo,
V. Lima,
T. Nakamura,
A. Saito,
Y. Satou,
J. -A. Scarpaci,
E. Takeshita,
S. Takeuchi,
T. Teranishi,
Y. Togano,
A. M. Vinodkumar,
Y. Yanagisawa
, et al. (1 additional authors not shown)
Abstract:
Coulomb excitation of the exotic neutron-rich nucleus Ne26 on a Pb208 target was measured at 58 MeV/u in order to search for low-lying E1 strength above the neutron emission threshold. This radioactive beam experiment was carried out at the RIKEN Accelerator Research Facility. Using the invariant mass method in the Ne25+n channel, we observe a sizable amount of E1 strength between 6 and 10 MeV exc…
▽ More
Coulomb excitation of the exotic neutron-rich nucleus Ne26 on a Pb208 target was measured at 58 MeV/u in order to search for low-lying E1 strength above the neutron emission threshold. This radioactive beam experiment was carried out at the RIKEN Accelerator Research Facility. Using the invariant mass method in the Ne25+n channel, we observe a sizable amount of E1 strength between 6 and 10 MeV excitation energy. By performing a multipole decomposition of the differential cross section, a reduced dipole transition probability of B(E1)=0.49+-0.16e2fm2 is deduced, corresponding to 4.9+-1.6% of the Thomas-Reiche-Kuhn sum rule. For the first time, the decay pattern of low-lying strength in a neutron-rich nucleus is measured. The extracted decay pattern is not consistent with several mean-field theory descriptions of the pygmy states.
△ Less
Submitted 11 July, 2023;
originally announced July 2023.
-
Intruder configurations in $^{29}$Ne at the transition into the island of inversion: Detailed structure study of $^{28}$Ne
Authors:
H. Wang,
M. Yasuda,
Y. Kondo,
T. Nakamura,
J. A. Tostevin,
K. Ogata,
T. Otsuka,
A. Poves,
N. Shimizu,
K. Yoshida,
N. L. Achouri,
H. Al Falou,
L. Atar,
T. Aumann,
H. Baba,
K. Boretzky,
C. Caesar,
D. Calvet,
H. Chae,
N. Chiga,
A. Corsi,
H. L. Crawford,
F. Delaunay,
A. Delbart,
Q. Deshayes
, et al. (71 additional authors not shown)
Abstract:
Detailed $γ$-ray spectroscopy of the exotic neon isotope $^{28}$Ne has been performed for the first time using the one-neutron removal reaction from $^{29}$Ne on a liquid hydrogen target at 240~MeV/nucleon. Based on an analysis of parallel momentum distributions, a level scheme with spin-parity assignments has been constructed for $^{28}$Ne and the negative-parity states are identified for the fir…
▽ More
Detailed $γ$-ray spectroscopy of the exotic neon isotope $^{28}$Ne has been performed for the first time using the one-neutron removal reaction from $^{29}$Ne on a liquid hydrogen target at 240~MeV/nucleon. Based on an analysis of parallel momentum distributions, a level scheme with spin-parity assignments has been constructed for $^{28}$Ne and the negative-parity states are identified for the first time. The measured partial cross sections and momentum distributions reveal a significant intruder $p$-wave strength providing evidence of the breakdown of the $N=20$ and $N=28$ shell gaps. Only a weak, possible $f$-wave strength was observed to bound final states. Large-scale shell-model calculations with different effective interactions do not reproduce the large $p$-wave and small $f$-wave strength observed experimentally, indicating an ongoing challenge for a complete theoretical description of the transition into the island of inversion along the Ne isotopic chain.
△ Less
Submitted 28 June, 2023;
originally announced June 2023.
-
Remark on the local well-posedness of compressible non-Newtonian fluids with initial vacuum
Authors:
Hind Al Baba,
Bilal Al Taki,
Amru Hussein
Abstract:
We discuss in this short note the local-in-time strong well-posedness of the compressible Navier-Stokes system for non-Newtonian fluids on the three dimensional torus. We show that the result established recently by Kalousek, Mácha, and Nečasova in \doi{10.1007/s00208-021-02301-8} can be extended to the case where vanishing density is allowed initially. Our proof builds on the framework developed…
▽ More
We discuss in this short note the local-in-time strong well-posedness of the compressible Navier-Stokes system for non-Newtonian fluids on the three dimensional torus. We show that the result established recently by Kalousek, Mácha, and Nečasova in \doi{10.1007/s00208-021-02301-8} can be extended to the case where vanishing density is allowed initially. Our proof builds on the framework developed by Cho, Choe, and Kim in \doi{10.1016/j.matpur.2003.11.004} for compressible Navier-Stokes equations in the case of Newtonian fluids. To adapt their method, special attention is given to the elliptic regularity of a challenging nonlinear elliptic system. We show particular results in this direction, however, the main result of this paper is proven in the general case when elliptic $W^{2,p}$-regularity is imposed as an assumption. Also, we give a finite time blow-up criterion.
△ Less
Submitted 8 April, 2024; v1 submitted 29 May, 2023;
originally announced May 2023.
-
Oblique projection for scalable rank-adaptive reduced-order modeling of nonlinear stochastic PDEs with time-dependent bases
Authors:
M. Donello,
G. Palkar,
M. H. Naderi,
D. C. Del Rey Fernández,
H. Babaee
Abstract:
Time-dependent basis reduced order models (TDB ROMs) have successfully been used for approximating the solution to nonlinear stochastic partial differential equations (PDEs). For many practical problems of interest, discretizing these PDEs results in massive matrix differential equations (MDEs) that are too expensive to solve using conventional methods. While TDB ROMs have the potential to signifi…
▽ More
Time-dependent basis reduced order models (TDB ROMs) have successfully been used for approximating the solution to nonlinear stochastic partial differential equations (PDEs). For many practical problems of interest, discretizing these PDEs results in massive matrix differential equations (MDEs) that are too expensive to solve using conventional methods. While TDB ROMs have the potential to significantly reduce this computational burden, they still suffer from the following challenges: (i) inefficient for general nonlinearities, (ii) intrusive implementation, (iii) ill-conditioned in the presence of small singular values, and (iv) error accumulation due to fixed rank. To this end, we present a scalable method based on oblique projections for solving TDB ROMs that is computationally efficient, minimally intrusive, robust in the presence of small singular values, rank-adaptive, and highly parallelizable. These favorable properties are achieved via low-rank approximation of the time discrete MDE. Using the discrete empirical interpolation method (DEIM), a low-rank decomposition is computed at each iteration of the time stepping scheme, enabling a near-optimal approximation at a fraction of the cost. We coin the new approach TDB-CUR since it is equivalent to a CUR decomposition based on sparse row and column samples of the MDE. We also propose a rank-adaptive procedure to control the error on-the-fly. Numerical results demonstrate the accuracy, efficiency, and robustness of the new method for a diverse set of problems.
△ Less
Submitted 16 August, 2023; v1 submitted 7 May, 2023;
originally announced May 2023.
-
Learning stiff chemical kinetics using extended deep neural operators
Authors:
Somdatta Goswami,
Ameya D. Jagtap,
Hessam Babaee,
Bryan T. Susi,
George Em Karniadakis
Abstract:
We utilize neural operators to learn the solution propagator for the challenging chemical kinetics equation. Specifically, we apply the deep operator network (DeepONet) along with its extensions, such as the autoencoder-based DeepONet and the newly proposed Partition-of-Unity (PoU-) DeepONet to study a range of examples, including the ROBERS problem with three species, the POLLU problem with 25 sp…
▽ More
We utilize neural operators to learn the solution propagator for the challenging chemical kinetics equation. Specifically, we apply the deep operator network (DeepONet) along with its extensions, such as the autoencoder-based DeepONet and the newly proposed Partition-of-Unity (PoU-) DeepONet to study a range of examples, including the ROBERS problem with three species, the POLLU problem with 25 species, pure kinetics of the syngas skeletal model for $CO/H_2$ burning, which contains 11 species and 21 reactions and finally, a temporally developing planar $CO/H_2$ jet flame (turbulent flame) using the same syngas mechanism. We have demonstrated the advantages of the proposed approach through these numerical examples. Specifically, to train the DeepONet for the syngas model, we solve the skeletal kinetic model for different initial conditions. In the first case, we parametrize the initial conditions based on equivalence ratios and initial temperature values. In the second case, we perform a direct numerical simulation of a two-dimensional temporally developing $CO/H_2$ jet flame. Then, we initialize the kinetic model by the thermochemical states visited by a subset of grid points at different time snapshots. Stiff problems are computationally expensive to solve with traditional stiff solvers. Thus, this work aims to develop a neural operator-based surrogate model to solve stiff chemical kinetics. The operator, once trained offline, can accurately integrate the thermochemical state for arbitrarily large time advancements, leading to significant computational gains compared to stiff integration schemes.
△ Less
Submitted 23 February, 2023;
originally announced February 2023.
-
STRASSE: A Silicon Tracker for Quasi-free Scattering Measurements at the RIBF
Authors:
H. N. Liu,
F. Flavigny,
H. Baba,
M. Boehmer,
U. Bonnes,
V. Borshchov,
P. Doornenbal,
N. Ebina,
M. Enciu,
A. Frotscher,
R. Gernhäuser,
V. Girard-Alcindor,
D. Goupillière,
J. Heuser,
R. Kapell,
Y. Kondo,
H. Lee,
J. Lehnert,
T. Matsui,
A. Matta,
T. Nakamura,
A. Obertelli,
T. Pohl,
M. Protsenko,
M. Sasano
, et al. (13 additional authors not shown)
Abstract:
STRASSE (Silicon Tracker for RAdioactive nuclei Studies at SAMURAI Experiments) is a new detection system under construction for quasi-free scattering (QFS) measurements at 200-250 MeV/nucleon at the RIBF facility of the RIKEN Nishina Center. It consists of a charged-particle silicon tracker coupled with a dedicated thick liquid hydrogen target (up to 150-mm long) in a compact geometry to fit insi…
▽ More
STRASSE (Silicon Tracker for RAdioactive nuclei Studies at SAMURAI Experiments) is a new detection system under construction for quasi-free scattering (QFS) measurements at 200-250 MeV/nucleon at the RIBF facility of the RIKEN Nishina Center. It consists of a charged-particle silicon tracker coupled with a dedicated thick liquid hydrogen target (up to 150-mm long) in a compact geometry to fit inside large scintillator or germanium arrays. Its design was optimized for two types of studies using QFS: missing-mass measurements and in-flight prompt $γ$-ray spectroscopy. This article describes (i) the resolution requirements needed to go beyond the sensitivity of existing systems for these two types of measurements, (ii) the conceptual design of the system using detailed simulations of the setup and (iii) its complete technical implementation and challenges. The final tracker aims at a sub-mm reaction vertex resolution and is expected to reach a missing-mass resolution below 2 MeV in $σ$ for $(p,2p)$ reactions when combined with the CsI(Na) CATANA array.
△ Less
Submitted 23 January, 2023;
originally announced January 2023.
-
Multiple Mechanisms in Proton-Induced Nucleon Removal at $\sim$100 MeV/Nucleon
Authors:
T. Pohl,
Y. L. Sun,
A. Obertelli,
J. Lee,
M. Gomez-Ramos,
K. Ogata,
K. Yoshida,
B. S. Cai,
C. X. Yuan,
B. A. Brown,
H. Baba,
D. Beaumel,
A. Corsi,
J. Gao,
J. Gibelin,
A. Gillibert,
K. I. Hahn,
T. Isobe,
D. Kim,
Y. Kondo,
T. Kobayashi,
Y. Kubota,
P. Li,
P. Liang,
H. N. Liu
, et al. (26 additional authors not shown)
Abstract:
We report on the first proton-induced single proton- and neutron-removal reactions from the neutron-deficient $^{14}$O nucleus with large Fermi-surface asymmetry $S_n-S_p$ = 18.6 MeV at $\sim$100 MeV/nucleon, a widely used energy regime for rare-isotope studies. The measured inclusive cross sections and parallel momentum distributions of the $^{13}$N and $^{13}$O residues are compared to the state…
▽ More
We report on the first proton-induced single proton- and neutron-removal reactions from the neutron-deficient $^{14}$O nucleus with large Fermi-surface asymmetry $S_n-S_p$ = 18.6 MeV at $\sim$100 MeV/nucleon, a widely used energy regime for rare-isotope studies. The measured inclusive cross sections and parallel momentum distributions of the $^{13}$N and $^{13}$O residues are compared to the state-of-the-art reaction models, with nuclear structure inputs from many-body shell-model calculations. Our results provide the first quantitative contributions of multiple reaction mechanisms including the quasifree knockout, inelastic scattering and nucleon transfer processes. It is shown that the inelastic scattering and nucleon transfer, usually neglected at such energy regime, contribute about 50% and 30% to the loosely bound proton and deeply bound neutron removal, respectively. These multiple reaction mechanisms should be considered in analyses of inclusive one-nucleon removal cross sections measured at intermediate energies for quantitative investigation of single-particle strengths and correlations in atomic nuclei.
△ Less
Submitted 27 April, 2023; v1 submitted 10 January, 2023;
originally announced January 2023.
-
Unbound states in 17C and p-sd shell-model interactions
Authors:
Sunji Kim,
Jongwon Hwang,
Yoshiteru Satou,
Nigel A. Orr,
Takashi Nakamura,
Yosuke Kondo,
Julien Gibelin,
Nadia Lynda Achouri,
Thomas Aumann,
Hidetada Baba,
Franck Delaunay,
Pieter Doornenbal,
Naoki Fukuda,
Naohito Inabe,
Tadaaki Isobe,
Daisuke Kameda,
Daiki Kanno,
Nobuyuki Kobayashi,
Toshio Kobayashi,
Toshiyuki Kubo,
Sylvain Leblond,
Jenny Lee,
Francisco Miguel Marqués,
Ryogo Minakata,
Tohru Motobayashi
, et al. (20 additional authors not shown)
Abstract:
Unbound states in 17C were investigated via one-neutron removal from a 18C beam at an energy of 245 MeV/nucleon on a carbon target. The energy spectrum of 17C, above the single-neutron decay threshold, was reconstructed using invariant mass spectroscopy from the measured momenta of the 16C fragment and neutron, and was found to exhibit resonances at Er=0.52(2), 0.77(2), 1.36(1), 1.91(1), 2.22(3) a…
▽ More
Unbound states in 17C were investigated via one-neutron removal from a 18C beam at an energy of 245 MeV/nucleon on a carbon target. The energy spectrum of 17C, above the single-neutron decay threshold, was reconstructed using invariant mass spectroscopy from the measured momenta of the 16C fragment and neutron, and was found to exhibit resonances at Er=0.52(2), 0.77(2), 1.36(1), 1.91(1), 2.22(3) and 3.20(1) MeV. The resonance at Er=0.77(2) MeV [Ex=1.51(3) MeV] was provisionally assigned as the second 5/2+ state. The two resonances at Er=1.91(1) and 3.20(1) MeV [Ex=2.65(2) and 3.94(2) MeV] were identified, through comparison of the energies, cross sections and momentum distributions with shell-model and eikonal reaction calculations, as p-shell hole states with spin-parities 1/2- and 3/2-, respectively. A detailed comparison was made with the results obtained using a range of shell-model interactions. The YSOX shell-model Hamiltonian, the cross-shell part of which is based on the monopole-based universal interaction, was found to provide a very good description of the present results and those for the neighbouring odd-A carbon isotopes - in particular for the negative parity cross-shell states.
△ Less
Submitted 25 December, 2022;
originally announced December 2022.
-
Isoscaling in central Sn+Sn collisions at 270 MeV/u
Authors:
J. W. Lee,
M. B. Tsang,
C. Y. Tsang,
R. Wang,
J. Barney,
J. Estee,
T. Isobe,
M. Kaneko,
M. Kurata-Nishimura,
W. G. Lynch,
T. Murakami,
A. Ono,
S. R. Souza,
D. S. Ahn,
L. Atar,
T. Aumann,
H. Baba,
K. Boretzky,
J. Brzychczyk,
G. Cerizza,
N. Chiga,
N. Fukuda,
I. Gasparic,
B. Hong,
A. Horvat
, et al. (39 additional authors not shown)
Abstract:
Experimental information on fragment emissions is important in understanding the dynamics of nuclear collisions and in the development of transport model simulating heavy-ion collisions. The composition of complex fragments emitted in the heavy-ion collisions can be explained by statistical models, which assume that thermal equilibrium is achieved at collision energies below 100 MeV/u. Our new exp…
▽ More
Experimental information on fragment emissions is important in understanding the dynamics of nuclear collisions and in the development of transport model simulating heavy-ion collisions. The composition of complex fragments emitted in the heavy-ion collisions can be explained by statistical models, which assume that thermal equilibrium is achieved at collision energies below 100 MeV/u. Our new experimental data together with theoretical analyses for light particles from Sn+Sn collisions at 270 MeV/u, suggest that the hypothesis of thermal equilibrium breaks down for particles emitted with high transfer momentum. To inspect the system's properties in such limit, the scaling features of the yield ratios of particles from two systems, a neutron-rich system of ${}^{132}\mathrm{Sn}+{}^{124}\mathrm{Sn}$ and a nearly symmetric system of ${}^{108}\mathrm{Sn}+{}^{112}\mathrm{Sn}$, are examined in the framework of the statistical multifragmentation model and the antisymmetrized molecular dynamics model. The isoscaling from low energy particles agree with both models. However the observed breakdown of isoscaling for particles with high transverse momentum cannot be explained by the antisymmetrized molecular dynamics model.
△ Less
Submitted 5 November, 2022;
originally announced November 2022.
-
Adaptive sparse interpolation for accelerating nonlinear stochastic reduced-order modeling with time-dependent bases
Authors:
Mohammad Hossein Naderi,
Hessam Babaee
Abstract:
Stochastic reduced-order modeling based on time-dependent bases (TDBs) has proven successful for extracting and exploiting low-dimensional manifold from stochastic partial differential equations (SPDEs). The nominal computational cost of solving a rank-$r$ reduced-order model (ROM) based on time-dependent basis, a.k.a. TDB-ROM, is roughly equal to that of solving the full-order model for $r$ rando…
▽ More
Stochastic reduced-order modeling based on time-dependent bases (TDBs) has proven successful for extracting and exploiting low-dimensional manifold from stochastic partial differential equations (SPDEs). The nominal computational cost of solving a rank-$r$ reduced-order model (ROM) based on time-dependent basis, a.k.a. TDB-ROM, is roughly equal to that of solving the full-order model for $r$ random samples. As of now, this nominal performance can only be achieved for linear or quadratic SPDEs -- at the expense of a highly intrusive process. On the other hand, for problems with non-polynomial nonlinearity, the computational cost of solving the TDB evolution equations is the same as solving the full-order model. In this work, we present an adaptive sparse interpolation algorithm that enables stochastic TDB-ROMs to achieve nominal computational cost for generic nonlinear SPDEs. Our algorithm constructs a low-rank approximation for the right hand side of the SPDE using the discrete empirical interpolation method (DEIM). The presented algorithm does not require any offline computation and as a result the low-rank approximation can adapt to any transient changes of the dynamics on the fly. We also propose a rank-adaptive strategy to control the error of the sparse interpolation. Our algorithm achieves computational speedup by adaptive sampling of the state and random spaces. We illustrate the efficiency of our approach for two test cases: (1) one-dimensional stochastic Burgers' equation, and (2) two-dimensional compressible Navier-Stokes equations subject to one-hundred-dimensional random perturbations. In all cases, the presented algorithm results in orders of magnitude reduction in the computational cost.
△ Less
Submitted 24 November, 2022; v1 submitted 14 July, 2022;
originally announced July 2022.
-
Border of the Island of Inversion: Unbound states in $^{29}$Ne
Authors:
M. Holl,
S. Lindberg,
A. Heinz,
Y. Kondo,
T. Nakamura,
J. A. Tostevin,
H. Wang,
T. Nilsson,
N. L. Achouri,
H. Al Falou,
L. Atar,
T. Aumann,
H. Baba,
K. Boretzky,
C. Caesar,
D. Calvet,
H. Chae,
N. Chiga,
A. Corsi,
H. L. Crawford,
F. Delaunay,
A. Delbart,
Q. Deshayes,
P. Díaz Fernández,
Z. Dombrádi
, et al. (67 additional authors not shown)
Abstract:
The nucleus $^{29}$Ne is situated at the border of the island of inversion. Despite significant efforts, no bound low-lying intruder $f_{7/2}$-state, which would place $^{29}$Ne firmly inside the island of inversion, has yet been observed. Here, the first investigation of unbound states of $^{29}$Ne is reported. The states were populated in $^{30}\mathrm{Ne}(p,pn)$ and $^{30}\mathrm{Na}(p,2p)$ rea…
▽ More
The nucleus $^{29}$Ne is situated at the border of the island of inversion. Despite significant efforts, no bound low-lying intruder $f_{7/2}$-state, which would place $^{29}$Ne firmly inside the island of inversion, has yet been observed. Here, the first investigation of unbound states of $^{29}$Ne is reported. The states were populated in $^{30}\mathrm{Ne}(p,pn)$ and $^{30}\mathrm{Na}(p,2p)$ reactions at a beam energy of around $230$ MeV/nucleon, and analyzed in terms of their resonance properties, partial cross sections and momentum distributions. The momentum distributions are compared to calculations using the eikonal, direct reaction model, allowing $\ell$-assignments for the observed states. The lowest-lying resonance at an excitation energy of 1.48(4) MeV shows clear signs of a significant $\ell$=3-component, giving first evidence for $f_{7/2}$ single particle strength in $^{29}$Ne. The excitation energies and strengths of the observed states are compared to shell-model calculations using the sdpf-u-mix interaction
△ Less
Submitted 11 February, 2022;
originally announced February 2022.
-
A First Glimpse at the Shell Structure beyond $^{54}$Ca: Spectroscopy of $^{55}$K, $^{55}$Ca, and $^{57}$Ca
Authors:
T. Koiwai,
K. Wimmer,
P. Doornenbal,
A. Obertelli,
C. Barbieri,
T. Duguet,
J. D. Holt,
T. Miyagi,
P. Navrátil,
K. Ogata,
N. Shimizu,
V. Somà,
Y. Utsuno,
K. Yoshida,
N. L. Achouri,
H. Baba,
F. Browne,
D. Calvet f,
F. Château,
S. Chen,
N. Chiga,
A. Corsi,
M. L. Cortés,
A. Delbart,
J. -M. Gheller
, et al. (58 additional authors not shown)
Abstract:
States in the $N=35$ and 37 isotopes $^{55,57}$Ca have been populated by direct proton-induced nucleon removal reactions from $^{56,58}$Sc and $^{56}$Ca beams at the RIBF. In addition, the $(p,2p)$ quasi-free single-proton removal reaction from $^{56}$Ca was studied. Excited states in $^{55}$K, $^{55}$Ca, and $^{57}$Ca were established for the first time via in-beam $γ$-ray spectroscopy. Results f…
▽ More
States in the $N=35$ and 37 isotopes $^{55,57}$Ca have been populated by direct proton-induced nucleon removal reactions from $^{56,58}$Sc and $^{56}$Ca beams at the RIBF. In addition, the $(p,2p)$ quasi-free single-proton removal reaction from $^{56}$Ca was studied. Excited states in $^{55}$K, $^{55}$Ca, and $^{57}$Ca were established for the first time via in-beam $γ$-ray spectroscopy. Results for the proton and neutron removal reactions from $^{56}$Ca to states in $^{55}$K and $^{55}$Ca for the level energies, excited state lifetimes, and exclusive cross sections agree well with state-of-the-art theoretical calculations using different approaches. The observation of a short-lived state in $^{57}$Ca suggests a transition in the calcium isotopic chain from single-particle dominated states at $N=35$ to collective excitations at $N=37$.
△ Less
Submitted 7 February, 2022;
originally announced February 2022.
-
Skeletal Reaction Models for Methane Combustion
Authors:
Yinmin Liu,
Hessam Babaee,
Peyman Givi,
Harsha Chelliah,
Daniel Livescu,
Arash Nouri
Abstract:
A local-sensitivity-analysis technique is employed to generate new skeletal reaction models for methane combustion from the foundational fuel chemistry model (FFCM-1). The sensitivities of the thermo-chemical variables with respect to the reaction rates are computed via the forced-optimally time dependent (f-OTD) methodology. In this methodology, the large sensitivity matrix containing all local s…
▽ More
A local-sensitivity-analysis technique is employed to generate new skeletal reaction models for methane combustion from the foundational fuel chemistry model (FFCM-1). The sensitivities of the thermo-chemical variables with respect to the reaction rates are computed via the forced-optimally time dependent (f-OTD) methodology. In this methodology, the large sensitivity matrix containing all local sensitivities is modeled as a product of two low-rank time-dependent matrices. The evolution equations of these matrices are derived from the governing equations of the system. The modeled sensitivities are computed for the auto-ignition of methane at atmospheric and high pressures with different sets of initial temperatures, and equivalence ratios. These sensitivities are then analyzed to rank the most important (sensitive) species. A series of skeletal models with different number of species and levels of accuracy in reproducing the FFCM-1 results are suggested. The performances of the generated models are compared against FFCM-1 in predicting the ignition delay, the laminar flame speed, and the flame extinction. The results of this comparative assessment suggest the skeletal models with 24 and more species generate the FFCM-1 results with an excellent accuracy.
△ Less
Submitted 28 August, 2023; v1 submitted 26 January, 2022;
originally announced January 2022.
-
Scalable In Situ Compression of Transient Simulation Data Using Time-Dependent Bases
Authors:
Shaghayegh Zamani Ashtiani,
Mujeeb R. Malik,
Hessam Babaee
Abstract:
Large-scale simulations of time-dependent problems generate a massive amount of data and with the explosive increase in computational resources the size of the data generated by these simulations has increased significantly. This has imposed severe limitations on the amount of data that can be stored and has elevated the issue of input/output (I/O) into one of the major bottlenecks of high-perform…
▽ More
Large-scale simulations of time-dependent problems generate a massive amount of data and with the explosive increase in computational resources the size of the data generated by these simulations has increased significantly. This has imposed severe limitations on the amount of data that can be stored and has elevated the issue of input/output (I/O) into one of the major bottlenecks of high-performance computing. In this work, we present an in situ compression technique to reduce the size of the data storage by orders of magnitude. This methodology is based on time-dependent subspaces and it extracts low-rank structures from multidimensional streaming data by decomposing the data into a set of time-dependent bases and a core tensor. We derive closed-form evolution equations for the core tensor as well as the time-dependent bases. The presented methodology does not require the data history and the computational cost of its extractions scales linearly with the size of data -- making it suitable for large-scale streaming datasets. To control the compression error, we present an adaptive strategy to add/remove modes to maintain the reconstruction error below a given threshold. We present four demonstration cases: (i) analytical example, (ii) incompressible unsteady reactive flow, (iii) stochastic turbulent reactive flow, and (iv) three-dimensional turbulent channel flow.
△ Less
Submitted 10 January, 2022;
originally announced January 2022.
-
Reduced Order Modeling of Turbulence-Chemistry Interactions using Dynamically Bi-Orthonormal Decomposition
Authors:
Aidyn Aitzhan,
Arash G. Nouri,
Peyman Givi,
Hessam Babaee
Abstract:
The performance of the dynamically bi-orthogonal (DBO) decomposition for the reduced order modeling of turbulence-chemistry interactions is assessed. DBO is an on-the-fly low-rank approximation technique, in which the instantaneous composition matrix of the reactive flow field is decomposed into a set of orthonormal spatial modes, a set of orthonormal vectors in the composition space, and a factor…
▽ More
The performance of the dynamically bi-orthogonal (DBO) decomposition for the reduced order modeling of turbulence-chemistry interactions is assessed. DBO is an on-the-fly low-rank approximation technique, in which the instantaneous composition matrix of the reactive flow field is decomposed into a set of orthonormal spatial modes, a set of orthonormal vectors in the composition space, and a factorization of the low-rank correlation matrix. Two factors which distinguish between DBO and the reduced order models (ROMs) based on the principal component analysis (PCA) are: (i) DBO does not require any offline data generation; and (ii) in DBO the low-rank composition subspace is time-dependent as opposed to static subspaces in PCA. Because of these features, DBO can adapt on-the-fly to intrinsic and externally excited transient changes in state of the transport variables. For demonstration, simulations are conducted of a non-premixed CO/H2 flame in a temporally evolving jet. The GRI-Mech 3.0 model with 53 species is used for chemical kinetics modeling. The results are appraised via a posteriori comparisons against data generated via full-rank direct numerical simulation (DNS) of the same flame, and the PCA reduction of the DNS data. The DBO also yields excellent predictions of various statistics of the thermo-chemical variables.
△ Less
Submitted 6 January, 2022;
originally announced January 2022.
-
Reduced order modeling with time-dependent bases for PDEs with stochastic boundary conditions
Authors:
Prerna Patil,
Hessam Babaee
Abstract:
Low-rank approximation using time-dependent bases (TDBs) has proven effective for reduced-order modeling of stochastic partial differential equations (SPDEs). In these techniques, the random field is decomposed to a set of deterministic TDBs and time-dependent stochastic coefficients. When applied to SPDEs with non-homogeneous stochastic boundary conditions (BCs), appropriate BC must be specified…
▽ More
Low-rank approximation using time-dependent bases (TDBs) has proven effective for reduced-order modeling of stochastic partial differential equations (SPDEs). In these techniques, the random field is decomposed to a set of deterministic TDBs and time-dependent stochastic coefficients. When applied to SPDEs with non-homogeneous stochastic boundary conditions (BCs), appropriate BC must be specified for each of the TDBs. However, determining BCs for TDB is not trivial because: (i) the dimension of the random BCs is different than the rank of the TDB subspace; (ii) TDB in most formulations must preserve orthonormality or orthogonality constraints and specifying BCs for TDB should not violate these constraints in the space-discretized form. In this work, we present a methodology for determining the boundary conditions for TDBs at no additional computational cost beyond that of solving the same SPDE with homogeneous BCs. Our methodology is informed by the fact the TDB evolution equations are the optimality conditions of a variational principle. We leverage the same variational principle to derive an evolution equation for the value of TDB at the boundaries. The presented methodology preserves the orthonormality or orthogonality constraints of TDBs. We present the formulation for both the dynamically bi-orthonormal (DBO) decomposition as well as the dynamically orthogonal (DO) decomposition. We show that the presented methodology can be applied to stochastic Dirichlet, Neumann, and Robin boundary conditions. We assess the performance of the presented method for linear advection-diffusion equation, Burgers' equation, and two-dimensional advection-diffusion equation with constant and temperature-dependent conduction coefficient.
△ Less
Submitted 7 November, 2022; v1 submitted 28 December, 2021;
originally announced December 2021.
-
First application of mass measurement with the Rare-RI Ring reveals the solar r-process abundance trend at A=122 and A=123
Authors:
H. F. Li,
S. Naimi,
T. M. Sprouse,
M. R. Mumpower,
Y. Abe,
Y. Yamaguchi,
D. Nagae,
F. Suzaki,
M. Wakasugi,
H. Arakawa,
W. B. Dou,
D. Hamakawa,
S. Hosoi,
Y. Inada,
D. Kajiki,
T. Kobayashi,
M. Sakaue,
Y. Yokoda,
T. Yamaguchi,
R. Kagesawa,
D. Kamioka,
T. Moriguchi,
M. Mukai,
A. Ozawa,
S. Ota
, et al. (19 additional authors not shown)
Abstract:
The Rare-RI Ring (R3) is a recently commissioned cyclotron-like storage ring mass spectrometer dedicated to mass measurements of exotic nuclei far from stability at Radioactive Isotope Beam Factory (RIBF) in RIKEN. The first application of mass measurement using the R3 mass spectrometer at RIBF is reported. Rare isotopes produced at RIBF, $^{127}$Sn, $^{126}$In, $^{125}$Cd, $^{124}$Ag, $^{123}$Pd,…
▽ More
The Rare-RI Ring (R3) is a recently commissioned cyclotron-like storage ring mass spectrometer dedicated to mass measurements of exotic nuclei far from stability at Radioactive Isotope Beam Factory (RIBF) in RIKEN. The first application of mass measurement using the R3 mass spectrometer at RIBF is reported. Rare isotopes produced at RIBF, $^{127}$Sn, $^{126}$In, $^{125}$Cd, $^{124}$Ag, $^{123}$Pd, were injected in R3. Masses of $^{126}$In, $^{125}$Cd, and $^{123}$Pd were measured whereby the mass uncertainty of $^{123}$Pd was improved. This is the first reported measurement with a new storage ring mass spectrometery technique realized at a heavy-ion cyclotron and employing individual injection of the pre-identified rare nuclei. The latter is essential for the future mass measurements of the rarest isotopes produced at RIBF. The impact of the new $^{123}$Pd result on the solar $r$-process abundances in a neutron star merger event is investigated by performing reaction network calculations of 20 trajectories with varying electron fraction $Y_e$. It is found that the neutron capture cross section on $^{123}$Pd increases by a factor of 2.2 and $β$-delayed neutron emission probability, $P_\mathrm{1n}$, of $^{123}$Rh increases by 14\%. The neutron capture cross section on $^{122}$Pd decreases by a factor of 2.6 leading to pileup of material at $A=122$, thus reproducing the trend of the solar $r$-process abundances. The trend of the two-neutron separation energies (S$_\mathrm{2n}$) was investigated for the Pd isotopic chain. The new mass measurement with improved uncertainty excludes large changes of the S$_\mathrm{2n}$ value at $N=77$. Such large increase of the S$_\mathrm{2n}$ values before $N=82$ was proposed as an alternative to the quenching of the $N=82$ shell gap to reproduce $r$-process abundances in the mass region of $A=112-124$.
△ Less
Submitted 16 February, 2022; v1 submitted 9 December, 2021;
originally announced December 2021.
-
Investigation of the ground-state spin inversion in the neutron-rich 47,49Cl isotopes
Authors:
B. D. Linh,
A. Corsi,
A. Gillibert,
A. Obertelli,
P. Doornenbal,
C. Barbieri,
S. Chen,
L. X. Chung,
T. Duguet,
M. Gómez-Ramos,
J. D. Holt,
A. Moro,
P. Navrátil,
K. Ogata,
N. T. T. Phuc,
N. Shimizu,
V. Somà,
Y. Utsuno,
N. L. Achouri,
H. Baba,
F. Browne,
D. Calvet,
F. Château,
N. Chiga,
M. L. Cortés
, et al. (61 additional authors not shown)
Abstract:
A first gamma-ray study of 47,49Cl spectroscopy was performed at the Radioactive Isotope Beam Factory with 50Ar projectiles at 217 MeV/nucleon, impinging on the liquid hydrogen target of the MINOS device. Prompt de-excitation gamma-rays were measured with the NaI(Tl) array DALI2+. Through the one-proton knockout reaction 50Ar(p,2p), a spin assignment could be determined for the low-lying states of…
▽ More
A first gamma-ray study of 47,49Cl spectroscopy was performed at the Radioactive Isotope Beam Factory with 50Ar projectiles at 217 MeV/nucleon, impinging on the liquid hydrogen target of the MINOS device. Prompt de-excitation gamma-rays were measured with the NaI(Tl) array DALI2+. Through the one-proton knockout reaction 50Ar(p,2p), a spin assignment could be determined for the low-lying states of 49Cl from the momentum distribution obtained with the SAMURAI spectrometer. A spin-parity J = 3/2+ is deduced for the ground state of 49Cl, similar to the recently studied N = 32 isotope 51K.
△ Less
Submitted 7 October, 2021;
originally announced October 2021.
-
Physics-informed neural networks for improving cerebral hemodynamics predictions
Authors:
Mohammad Sarabian,
Hessam Babaee,
Kaveh Laksari
Abstract:
Determining brain hemodynamics plays a critical role in the diagnosis and treatment of various cerebrovascular diseases. In this work, we put forth a physics-informed deep learning framework that augments sparse clinical measurements with fast computational fluid dynamics (CFD) simulations to generate physically consistent and high spatiotemporal resolution of brain hemodynamic parameters. Transcr…
▽ More
Determining brain hemodynamics plays a critical role in the diagnosis and treatment of various cerebrovascular diseases. In this work, we put forth a physics-informed deep learning framework that augments sparse clinical measurements with fast computational fluid dynamics (CFD) simulations to generate physically consistent and high spatiotemporal resolution of brain hemodynamic parameters. Transcranial Doppler (TCD) ultrasound is one of the most common techniques in the current clinical workflow that enables noninvasive and instantaneous evaluation of blood flow velocity within the cerebral arteries. However, it is spatially limited to only a handful of locations across the cerebrovasculature due to the constrained accessibility through the skull's acoustic windows. Our deep learning framework employs in-vivo real-time TCD velocity measurements at several locations in the brain and the baseline vessel cross-sectional areas acquired from 3D angiography images, and provides high-resolution maps of velocity, area, and pressure in the entire vasculature. We validated the predictions of our model against in-vivo velocity measurements obtained via 4D flow MRI scans. We then showcased the clinical significance of this technique in diagnosing the cerebral vasospasm (CVS) by successfully predicting the changes in vasospastic local vessel diameters based on corresponding sparse velocities measurements. The key finding here is that the combined effects of uncertainties in outlet boundary condition subscription and modeling physics deficiencies render the conventional purely physics-based computational models unsuccessful in recovering accurate brain hemodynamics. Nonetheless, fusing these models with clinical measurements through a data-driven approach ameliorates predictions of brain hemodynamic variables.
△ Less
Submitted 25 August, 2021;
originally announced August 2021.
-
Isomeric states in neutron-rich nuclei around $N = 40$
Authors:
K. Wimmer,
F. Recchia,
S. M. Lenzi,
S. Riccetto,
T. Davinson,
A. Estrade,
C. J. Griffin,
S. Nishimura,
V. Phong,
P. -A. Söderström,
O. Aktas,
M. Al-Aqeel,
T. Ando,
H. Baba,
S. Bae,
S. Choi,
P. Doornenbal,
J. Ha,
L. Harkness-Brennan,
T. Isobe,
P. R. John,
D. Kahl,
G. Kiss,
I. Kojouharov,
N. Kurz
, et al. (15 additional authors not shown)
Abstract:
Neutron-rich nuclei in the vicinity of the $N=40$ island of inversion are characterized by shell evolution and exhibit deformed ground states. In several nuclei isomeric states have been observed and attributed to excitations to the intruder neutron $1g_{9/2}$ orbital. In the present study we searched for isomeric states in nuclei around $N=40$, $Z=22$ produced by projectile fragmentation at RIBF.…
▽ More
Neutron-rich nuclei in the vicinity of the $N=40$ island of inversion are characterized by shell evolution and exhibit deformed ground states. In several nuclei isomeric states have been observed and attributed to excitations to the intruder neutron $1g_{9/2}$ orbital. In the present study we searched for isomeric states in nuclei around $N=40$, $Z=22$ produced by projectile fragmentation at RIBF. Delayed $γ$ rays were detected by the EURICA germanium detector array. High statistics data allowed for an updated decay scheme of $^{60}$V. The lifetime of an isomeric state in $^{64}$V was measured for the first time in the present experiment. A previously unobserved isomeric state was discovered in $^{58}$Sc. The measured lifetime suggests a parity changing transition, originating from an odd number of neutrons in the $1g_{9/2}$ orbital. The nature of the isomeric state in $^{58}$Sc is thus different from isomers in the less exotic V and Sc nuclei.
△ Less
Submitted 25 June, 2021;
originally announced June 2021.
-
A Quantum Inspired Approach to Exploit Turbulence Structures
Authors:
Nikita Gourianov,
Michael Lubasch,
Sergey Dolgov,
Quincy Y. van den Berg,
Hessam Babaee,
Peyman Givi,
Martin Kiffner,
Dieter Jaksch
Abstract:
Understanding turbulence is the key to our comprehension of many natural and technological flow processes. At the heart of this phenomenon lies its intricate multi-scale nature, describing the coupling between different-sized eddies in space and time. Here we introduce a new paradigm for analyzing the structure of turbulent flows by quantifying correlations between different length scales using me…
▽ More
Understanding turbulence is the key to our comprehension of many natural and technological flow processes. At the heart of this phenomenon lies its intricate multi-scale nature, describing the coupling between different-sized eddies in space and time. Here we introduce a new paradigm for analyzing the structure of turbulent flows by quantifying correlations between different length scales using methods inspired from quantum many-body physics. We present results for interscale correlations of two paradigmatic flow examples, and use these insights along with tensor network theory to design a structure-resolving algorithm for simulating turbulent flows. With this algorithm, we find that the incompressible Navier-Stokes equations can be accurately solved within a computational space reduced by over an order of magnitude compared to direct numerical simulation. Our quantum-inspired approach provides a pathway towards conducting computational fluid dynamics on quantum computers.
△ Less
Submitted 4 July, 2022; v1 submitted 10 June, 2021;
originally announced June 2021.
-
Skeletal Model Reduction with Forced Optimally Time Dependent Modes
Authors:
A. G. Nouri,
H. Babaee,
P. Givi,
H. K. Chelliah,
D. Livescu
Abstract:
Sensitivity analysis with forced optimally time dependent (f-OTD) modes is introduced and its application for skeletal model reduction is demonstrated. f-OTD expands the sensitivity coefficient matrix into a lowdimensional, time dependent, orthonormal basis which captures directions of the phase space associated with most dominant sensitivities. These directions highlight the instantaneous active…
▽ More
Sensitivity analysis with forced optimally time dependent (f-OTD) modes is introduced and its application for skeletal model reduction is demonstrated. f-OTD expands the sensitivity coefficient matrix into a lowdimensional, time dependent, orthonormal basis which captures directions of the phase space associated with most dominant sensitivities. These directions highlight the instantaneous active species and reaction paths. Evolution equations for the orthonormal basis and the projections of sensitivity matrix onto the basis are derived, and the application of f-OTD for skeletal reduction is described. In this framework, the sensitivity matrix is modeled, stored in a factorized manner, and never reconstructed at any time during the calculations. For demonstration purposes, sensitivity analysis is conducted of constant pressure ethylene-air burning in a zero-dimensional reactor and new skeletal models are generated. The flame speed, the ignition delay, and the extinction curve of resulted models are compared against some of the existing skeletal models. The results demonstrate the ability of f-OTD approach to eliminate unimportant reactions and species in a systematic, efficient and accurate manner.
△ Less
Submitted 3 March, 2021;
originally announced March 2021.
-
Probing the Symmetry Energy with the Spectral Pion Ratio
Authors:
J. Estee,
W. G. Lynch,
C. Y. Tsang,
J. Barney,
G. Jhang,
M. B. Tsang,
R. Wang,
M. Kaneko,
J. W. Lee,
T. Isobe,
M. Kurata-Nishimura,
T. Murakami,
D. S. Ahn,
L. Atar,
T. Aumann,
H. Baba,
K. Boretzky,
J. Brzychczyk,
G. Cerizza,
N. Chiga,
N. Fukuda,
I. Gasparic,
B. Hong,
A. Horvat,
K. Ieki
, et al. (38 additional authors not shown)
Abstract:
Many neutron star (NS) properties, such as the proton fraction within a NS, reflect the symmetry energy contributions to the Equation of State that dominate when neutron and proton densities differ strongly. To constrain these contributions at supra-saturation densities, we measure the spectra of charged pions produced by colliding rare isotope tin (Sn) beams with isotopically enriched Sn targets.…
▽ More
Many neutron star (NS) properties, such as the proton fraction within a NS, reflect the symmetry energy contributions to the Equation of State that dominate when neutron and proton densities differ strongly. To constrain these contributions at supra-saturation densities, we measure the spectra of charged pions produced by colliding rare isotope tin (Sn) beams with isotopically enriched Sn targets. Using ratios of the charged pion spectra measured at high transverse momenta, we deduce the slope of the symmetry energy to be $42 < L < 117$ MeV. This value is slightly lower but consistent with the $L$ values deduced from a recent measurement of the neutron skin thickness of $^{208}$Pb.
△ Less
Submitted 11 March, 2021;
originally announced March 2021.
-
Quasi-free Neutron Knockout Reaction Reveals a Small $s$-orbital Component in the Borromean Nucleus $^{17}$B
Authors:
Z. H. Yang,
Y. Kubota,
A. Corsi,
K. Yoshida,
X. -X. Sun,
J. G. Li,
M. Kimura,
N. Michel,
K. Ogata,
C. X. Yuan,
Q. Yuan,
G. Authelet,
H. Baba,
C. Caesar,
D. Calvet,
A. Delbart,
M. Dozono,
J. Feng,
F. Flavigny,
J. -M. Gheller,
J. Gibelin,
A. Giganon,
A. Gillibert,
K. Hasegawa,
T. Isobe
, et al. (51 additional authors not shown)
Abstract:
A kinematically complete quasi-free $(p,pn)$ experiment in inverse kinematics was performed to study the structure of the Borromean nucleus $^{17}$B, which had long been considered to have neutron halo. By analyzing the momentum distributions and exclusive cross sections, we obtained the spectroscopic factors for $1s_{1/2}$ and $0d_{5/2}$ orbitals, and a surprisingly small percentage of 9(2)$\%$ w…
▽ More
A kinematically complete quasi-free $(p,pn)$ experiment in inverse kinematics was performed to study the structure of the Borromean nucleus $^{17}$B, which had long been considered to have neutron halo. By analyzing the momentum distributions and exclusive cross sections, we obtained the spectroscopic factors for $1s_{1/2}$ and $0d_{5/2}$ orbitals, and a surprisingly small percentage of 9(2)$\%$ was determined for $1s_{1/2}$. Our finding of such a small $1s_{1/2}$ component and the halo features reported in prior experiments can be explained by the deformed relativistic Hartree-Bogoliubov theory in continuum, revealing a definite but not dominant neutron halo in $^{17}$B. The present work gives the smallest $s$- or $p$-orbital component among known nuclei exhibiting halo features, and implies that the dominant occupation of $s$ or $p$ orbitals is not a prerequisite for the occurrence of neutron halo.
△ Less
Submitted 6 February, 2021;
originally announced February 2021.
-
Spectroscopy of $^{33}$Mg with knockout reactions
Authors:
D. Bazin,
N. Aoi,
H. Baba,
J. Chen,
H. Crawford,
P. Doornenbal,
P. Fallon,
K. Li,
J. Lee,
M. Matsushita,
T. Motobayashi,
H. Sakurai,
H. Scheit,
D. Steppenbeck,
R. Stroberg,
S. Takeuchi,
H. Wang,
K. Yoneda,
C. X. Yuan
Abstract:
The structure of $^{33}$Mg was investigated by means of two knockout reactions, one-neutron removal from $^{34}$Mg and one-proton removal from $^{34}$Al. Using comparative analysis of the population of observed excited states in the residual $^{33}$Mg, the nature of these states can be deciphered. In addition, the long-standing controversy about the parity of the $^{33}$Mg ground state is resolved…
▽ More
The structure of $^{33}$Mg was investigated by means of two knockout reactions, one-neutron removal from $^{34}$Mg and one-proton removal from $^{34}$Al. Using comparative analysis of the population of observed excited states in the residual $^{33}$Mg, the nature of these states can be deciphered. In addition, the long-standing controversy about the parity of the $^{33}$Mg ground state is resolved using momentum distribution analysis, showing a clear signature for negative parity. Partial cross section measurements are compared with the results of eikonal reaction theory combined with large-scale shell model calculations of this complex nucleus located in the island of inversion, where configuration mixing plays a major role.
△ Less
Submitted 28 January, 2021;
originally announced January 2021.
-
Persistence of the ${Z=28}$ shell gap in ${A=75}$ isobars: Identification of a possible ${(1/2^-)}$ $μ$s isomer in ${^{75}}$Co and $β$ decay to ${^{75}}$Ni
Authors:
S. Escrig,
A. I. Morales,
S. Nishimura,
M. Niikura,
A. Poves,
Z. Y. Xu,
G. Lorusso,
F. Browne,
P. Doornenbal,
G. Gey,
H. -S. Jung,
Z. Li,
P. -A. Söderström,
T. Sumikama,
J. Taprogge,
Zs. Vajta,
H. Watanabe,
J. Wu,
A. Yagi,
K. Yoshinaga,
H. Baba,
S. Franchoo,
T. Isobe,
P. R. John,
I. Kojouharov
, et al. (18 additional authors not shown)
Abstract:
Background: The evolution of shell structure around doubly magic exotic nuclei is of great interest in nuclear physics and astrophysics. In the `southwest' region of $^{78}$Ni, the development of deformation might trigger a major shift in our understanding of explosive nucleosynthesis. To this end, new spectroscopic information on key close-lying nuclei is very valuable.
Purpose: We intend to me…
▽ More
Background: The evolution of shell structure around doubly magic exotic nuclei is of great interest in nuclear physics and astrophysics. In the `southwest' region of $^{78}$Ni, the development of deformation might trigger a major shift in our understanding of explosive nucleosynthesis. To this end, new spectroscopic information on key close-lying nuclei is very valuable.
Purpose: We intend to measure the isomeric and $β$ decay of $^{75}$Co, with one-proton and two-neutron holes relative to $^{78}$Ni, to access new nuclear structure information in $^{75}$Co and its $β$-decay daughters $^{75}$Ni and $^{74}$Ni.
Methods: The nucleus $^{75}$Co is produced in relativistic in-flight fission reactions of $^{238}$U at the Radioactive Ion Beam Factory in the RIKEN Nishina Center. Its isomeric and $β$ decay are studied exploiting the BigRIPS and EURICA setups.
Results: We obtain partial $β$-decay spectra for $^{75}$Ni and $^{74}$Ni, and report a new isomeric transition in $^{75}$Co. The energy [$E_γ=1914(2)$ keV] and half-life [$t_{1/2}=13(6)$ $μ$s] of the delayed $γ$ ray lend support for the existence of a $J^π=(1/2^-)$ isomeric state at 1914(2) keV. A comparison with PFSDG-U shell-model calculations provides a good account for the observed states in $^{75}$Ni, but the first calculated $1/2^-$ level in $^{75}$Co, a prolate $K=1/2$ state, is predicted about 1 MeV below the observed $(1/2^-)$ level.
Conclusions: The spherical-like structure of the lowest-lying excited states in $^{75}$Ni is proved. In the case of $^{75}$Co, the results suggest that the dominance of the spherical configurations over the deformed ones might be stronger than expected below $^{78}$Ni. Further experimental efforts to discern the nature of the $J^π=(1/2^-)$ isomer are necessary.
△ Less
Submitted 23 June, 2021; v1 submitted 15 January, 2021;
originally announced January 2021.
-
On-the-fly Reduced Order Modeling of Passive and Reactive Species via Time-Dependent Manifolds
Authors:
Donya Ramezanian,
Arash G. Nouri,
Hessam Babaee
Abstract:
One of the principal barriers in developing accurate and tractable predictive models in turbulent flows with a large number of species is to track every species by solving a separate transport equation, which can be computationally impracticable. In this paper, we present an on-the-fly reduced order modeling of reactive as well as passive transport equations to reduce the computational cost. The p…
▽ More
One of the principal barriers in developing accurate and tractable predictive models in turbulent flows with a large number of species is to track every species by solving a separate transport equation, which can be computationally impracticable. In this paper, we present an on-the-fly reduced order modeling of reactive as well as passive transport equations to reduce the computational cost. The presented approach seeks a low-rank decomposition of the species to three time-dependent components: (i) a set of orthonormal spatial modes, (ii) a low-rank factorization of the instantaneous species correlation matrix, and (iii) a set of orthonormal species modes, which represent a low-dimensional time-dependent manifold. Our approach bypasses the need to solve the full-dimensional species to generate high-fidelity data - as it is commonly performed in data-driven dimension reduction techniques such as the principle component analysis. Instead, the low-rank components are directly extracted from the species transport equation. The evolution equations for the three components are obtained from optimality conditions of a variational principle. The time-dependence of the three components enables an on-the-fly adaptation of the low-rank decomposition to transient changes in the species. Several demonstration cases of reduced order modeling of passive and reactive transport equations are presented.
△ Less
Submitted 11 January, 2021;
originally announced January 2021.
-
Shape changes in the mirror nuclei $^{70}$Kr and $^{70}$Se
Authors:
K. Wimmer,
W. Korten,
P. Doornenbal,
T. Arici,
P. Aguilera,
A. Algora,
T. Ando,
H. Baba,
B. Blank,
A. Boso,
S. Chen,
A. Corsi,
P. Davies,
G. de Angelis,
G. de France,
J. -P. Delaroche,
D. T. Doherty,
J. Gerl,
R. Gernhäuser,
M. Girod,
D. Jenkins,
S. Koyama,
T. Motobayashi,
S. Nagamine,
M. Niikura
, et al. (13 additional authors not shown)
Abstract:
We studied the proton-rich $T_z=-1$ nucleus $^{70}$Kr through inelastic scattering at intermediate energies in order to extract the reduced transition probability, $B(E2;\;0^+ \rightarrow 2^+)$. Comparison with the other members of the $A=70$ isospin triplet, $^{70}$Br and $^{70}$Se, studied in the same experiment, shows a $3σ$ deviation from the expected linearity of the electromagnetic matrix el…
▽ More
We studied the proton-rich $T_z=-1$ nucleus $^{70}$Kr through inelastic scattering at intermediate energies in order to extract the reduced transition probability, $B(E2;\;0^+ \rightarrow 2^+)$. Comparison with the other members of the $A=70$ isospin triplet, $^{70}$Br and $^{70}$Se, studied in the same experiment, shows a $3σ$ deviation from the expected linearity of the electromagnetic matrix elements as a function of $T_z$.
At present, no established nuclear structure theory can describe this observed deviation quantitatively. This is the first violation of isospin symmetry at this level observed in the transition matrix elements. A heuristic approach may explain the anomaly by a shape change between the mirror nuclei $^{70}$Kr and $^{70}$Se contrary to the model predictions.
△ Less
Submitted 6 January, 2021;
originally announced January 2021.
-
Computing Sensitivities in Evolutionary Systems: A Real-Time Reduced Order Modeling Strategy
Authors:
Michael Donello,
Mark Carpenter,
Hessam Babaee
Abstract:
We present a new methodology for computing sensitivities in evolutionary systems using a model-driven low-rank approximation. To this end, we formulate a variational principle that seeks to minimize the distance between the time derivative of the reduced approximation and sensitivity dynamics. The first-order optimality condition of the variational principle leads to a system of closed-form evolut…
▽ More
We present a new methodology for computing sensitivities in evolutionary systems using a model-driven low-rank approximation. To this end, we formulate a variational principle that seeks to minimize the distance between the time derivative of the reduced approximation and sensitivity dynamics. The first-order optimality condition of the variational principle leads to a system of closed-form evolution equations for an orthonormal basis and corresponding sensitivity coefficients. This approach allows for the computation of sensitivities with respect to a large number of parameters in an accurate and tractable manner by extracting correlations between different sensitivities on the fly. The presented method requires solving forward evolution equations, sidestepping the restrictions imposed by forward/backward workflow of adjoint sensitivities. For example, the presented method, unlike the adjoint equation, does not impose any I/O load and can be used in applications in which real time sensitivities are of interest. We demonstrate the utility of the method for three test cases: (1) computing sensitivity with respect to model parameters in the Rossler system (2) computing sensitivity with respect to an infinite-dimensional forcing parameter in the chaotic Kuramoto-Sivashinsky equation and (3) computing sensitivity with respect to reaction parameters for species transport in a turbulent reacting flow.
△ Less
Submitted 27 December, 2020;
originally announced December 2020.
-
Multi-fidelity Bayesian Neural Networks: Algorithms and Applications
Authors:
Xuhui Meng,
Hessam Babaee,
George Em Karniadakis
Abstract:
We propose a new class of Bayesian neural networks (BNNs) that can be trained using noisy data of variable fidelity, and we apply them to learn function approximations as well as to solve inverse problems based on partial differential equations (PDEs). These multi-fidelity BNNs consist of three neural networks: The first is a fully connected neural network, which is trained following the maximum a…
▽ More
We propose a new class of Bayesian neural networks (BNNs) that can be trained using noisy data of variable fidelity, and we apply them to learn function approximations as well as to solve inverse problems based on partial differential equations (PDEs). These multi-fidelity BNNs consist of three neural networks: The first is a fully connected neural network, which is trained following the maximum a posteriori probability (MAP) method to fit the low-fidelity data; the second is a Bayesian neural network employed to capture the cross-correlation with uncertainty quantification between the low- and high-fidelity data; and the last one is the physics-informed neural network, which encodes the physical laws described by PDEs. For the training of the last two neural networks, we use the Hamiltonian Monte Carlo method to estimate accurately the posterior distributions for the corresponding hyperparameters. We demonstrate the accuracy of the present method using synthetic data as well as real measurements. Specifically, we first approximate a one- and four-dimensional function, and then infer the reaction rates in one- and two-dimensional diffusion-reaction systems. Moreover, we infer the sea surface temperature (SST) in the Massachusetts and Cape Cod Bays using satellite images and in-situ measurements. Taken together, our results demonstrate that the present method can capture both linear and nonlinear correlation between the low- and high-fideilty data adaptively, identify unknown parameters in PDEs, and quantify uncertainties in predictions, given a few scattered noisy high-fidelity data. Finally, we demonstrate that we can effectively and efficiently reduce the uncertainties and hence enhance the prediction accuracy with an active learning approach, using as examples a specific one-dimensional function approximation and an inverse PDE problem.
△ Less
Submitted 18 December, 2020;
originally announced December 2020.
-
Symmetry energy investigation with pion production from Sn+Sn systems
Authors:
G. Jhang,
J. Estee,
J. Barney,
G. Cerizza,
M. Kaneko,
J. W. Lee,
W. G. Lynch,
T. Isobe,
M. Kurata-Nishimura,
T. Murakami,
C. Y . Tsang,
M. B. Tsang,
R. Wang,
D. S. Ahn,
L. Atar,
T. Aumann,
H. Baba,
K. Boretzky,
J. Brzychczyk,
N. Chiga,
N. Fukuda,
I. Gasparic,
B. Hong,
A. Horvat,
K. Ieki
, et al. (55 additional authors not shown)
Abstract:
In the past two decades, pions created in the high density regions of heavy ion collisions have been predicted to be sensitive at high densities to the symmetry energy term in the nuclear equation of state, a property that is key to our understanding of neutron stars. In a new experiment designed to study the symmetry energy, the multiplicities of negatively and positively charged pions have been…
▽ More
In the past two decades, pions created in the high density regions of heavy ion collisions have been predicted to be sensitive at high densities to the symmetry energy term in the nuclear equation of state, a property that is key to our understanding of neutron stars. In a new experiment designed to study the symmetry energy, the multiplicities of negatively and positively charged pions have been measured with high accuracy for central $^{132}$Sn+$^{124}$Sn, $^{112}$Sn+$^{124}$Sn, and $^{108}$Sn+$^{112}$Sn collisions at $E/A=270~\mathrm{MeV}$ with the S$π$RIT Time Projection Chamber. While the uncertainties of individual pion multiplicities are measured to 4\%, those of the charged pion multiplicity ratios are measured to 2\%. We compare these data to predictions from seven major transport models. The calculations reproduce qualitatively the dependence of the multiplicities and their ratios on the total neutron to proton number in the colliding systems. However, the predictions of the transport models from different codes differ too much to allow extraction of reliable constraints on the symmetry energy from the data. This finding may explain previous contradictory conclusions on symmetry energy constraints obtained from pion data in Au+Au system. These new results call for better understanding of the differences among transport codes, and new observables that are more sensitive to the density dependence of the symmetry energy.
△ Less
Submitted 13 December, 2020;
originally announced December 2020.
-
$\boldsymbol{N=32}$ shell closure below calcium: Low-lying structure of $^{50}$Ar
Authors:
M. L. Cortés,
W. Rodriguez,
P. Doornenbal,
A. Obertelli,
J. D. Holt,
J. Menéndez,
K. Ogata,
A. Schwenk,
N. Shimizu,
J. Simonis,
Y. Utsuno,
K. Yoshida,
L. Achouri,
H. Baba,
F. Browne,
D. Calvet,
F. Château,
S. Chen,
N. Chiga,
A. Corsi,
A. Delbart,
J-M. Gheller,
A. Giganon,
A. Gillibert,
C. Hilaire
, et al. (56 additional authors not shown)
Abstract:
Low-lying excited states in the $N=32$ isotope $^{50}$Ar were investigated by in-beam $γ$-ray spectroscopy following proton- and neutron-knockout, multi-nucleon removal, and proton inelastic scattering at the RIKEN Radioactive Isotope Beam Factory. The energies of the two previously reported transitions have been confirmed, and five additional states are presented for the first time, including a c…
▽ More
Low-lying excited states in the $N=32$ isotope $^{50}$Ar were investigated by in-beam $γ$-ray spectroscopy following proton- and neutron-knockout, multi-nucleon removal, and proton inelastic scattering at the RIKEN Radioactive Isotope Beam Factory. The energies of the two previously reported transitions have been confirmed, and five additional states are presented for the first time, including a candidate for a 3$^-$ state. The level scheme built using $γγ$ coincidences was compared to shell-model calculations in the $sd-pf$ model space, and to ab initio predictions based on chiral two- and three-nucleon interactions. Theoretical proton- and neutron-knockout cross sections suggest that two of the new transitions correspond to $2^+$ states, while the previously proposed $4^+$ state could also correspond to a $2^+$ state.
△ Less
Submitted 21 November, 2020;
originally announced November 2020.
-
Development and operation of an electrostatic time-of-flight detector for the Rare RI storage Ring
Authors:
D. Nagae,
Y. Abe,
S. Okada,
S. Omika,
K. Wakayama,
S. Hosoi,
S. Suzuki,
T. Moriguchi,
M. Amano,
D. Kamioka,
Z. Ge,
S. Naimi,
F. Suzaki,
N. Tadano,
R. Igosawa,
K. Inomata,
H. Arakawa,
K. Nishimuro,
T. Fujii,
T. Mitsui,
Y. Yanagisawa,
H. Baba,
S. Michimasa,
S. Ota,
G. Lorusso
, et al. (6 additional authors not shown)
Abstract:
An electrostatic time-of-flight detector named E-MCP has been developed for quick diagnostics of circulating beam and timing measurement in mass spectrometry at the Rare-RI Ring in RIKEN. The E-MCP detector consists of a conversion foil, potential grids, and a microchannel plate. Secondary electrons are released from the surface of the foil when a heavy ion hits it. The electrons are accelerated a…
▽ More
An electrostatic time-of-flight detector named E-MCP has been developed for quick diagnostics of circulating beam and timing measurement in mass spectrometry at the Rare-RI Ring in RIKEN. The E-MCP detector consists of a conversion foil, potential grids, and a microchannel plate. Secondary electrons are released from the surface of the foil when a heavy ion hits it. The electrons are accelerated and deflected by 90$^\circ$ toward the microchannel plate by electrostatic potentials. A thin carbon foil and a thin aluminum-coated mylar foil were used as conversion foils. We obtained time resolutions of 69(1) ps and 43(1) ps (standard deviation) for a $^{84}$Kr beam at an energy of 170 MeV/u when using the carbon and the aluminum-coated mylar foils, respectively. A detection efficiency of approximately 90% was obtained for both foils. The E-MCP detector equipped with the carbon foil was installed inside the Rare-RI Ring to confirm particle circulation within a demonstration experiment on mass measurements of nuclei around $^{78}$Ge produced by in-flight fission of uranium beam at the RI Beam Factory in RIKEN. Periodic time signals from circulating ions were clearly observed. Revolution times for $^{78}$Ge, $^{77}$Ga, and $^{76}$Zn were obtained. The results confirmed successful circulation of the short-lived nuclei inside the Rare-RI Ring.
△ Less
Submitted 3 November, 2020;
originally announced November 2020.
-
Surface localization of the dineutron in $^{11}$Li
Authors:
Y. Kubota,
A. Corsi,
G. Authelet,
H. Baba,
C. Caesar,
D. Calvet,
A. Delbart,
M. Dozono,
J. Feng,
F. Flavigny,
J. -M. Gheller,
J. Gibelin,
A. Giganon,
A. Gillibert,
K. Hasegawa,
T. Isobe,
Y. Kanaya,
S. Kawakami,
D. Kim,
Y. Kikuchi,
Y. Kiyokawa,
M. Kobayashi,
N. Kobayashi,
T. Kobayashi,
Y. Kondo
, et al. (42 additional authors not shown)
Abstract:
The formation of a dineutron in the nucleus $^{11}$Li is found to be localized to the surface region. The experiment measured the intrinsic momentum of the struck neutron in $^{11}$Li via the $(p,pn)$ knockout reaction at 246 MeV/nucleon. The correlation angle between the two neutrons is, for the first time, measured as a function of the intrinsic neutron momentum. A comparison with reaction calcu…
▽ More
The formation of a dineutron in the nucleus $^{11}$Li is found to be localized to the surface region. The experiment measured the intrinsic momentum of the struck neutron in $^{11}$Li via the $(p,pn)$ knockout reaction at 246 MeV/nucleon. The correlation angle between the two neutrons is, for the first time, measured as a function of the intrinsic neutron momentum. A comparison with reaction calculations reveals the localization of the dineutron at $r\sim3.6$ fm. The results also support the density dependence of dineutron formation as deduced from Hartree-Fock-Bogoliubov calculations for nuclear matter.
△ Less
Submitted 9 October, 2020;
originally announced October 2020.
-
The S$π$RIT Time Projection Chamber
Authors:
J. Barney,
J. Estee,
W. G. Lynch,
T. Isobe,
G. Jhang,
M. Kurata-Nishimura,
A. B. McIntosh,
T. Murakami,
R. Shane,
S. Tangwancharoen,
M. B. Tsang,
G. Cerizza,
M. Kaneko,
J. W. Lee,
C. Y. Tsang,
R. Wang,
C. Anderson,
H. Baba,
Z. Chajecki,
M. Famiano,
R. Hodges-Showalter,
B. Hong,
T. Kobayashi,
P. Lasko,
J. Łukasik
, et al. (15 additional authors not shown)
Abstract:
The SAMURAI Pion Reconstruction and Ion-Tracker Time Projection Chamber (S$π$RIT TPC) was designed to enable measurements of heavy ion collisions with the SAMURAI spectrometer at the RIKEN Radioactive Isotope Beam Factory and provide constraints on the Equation of State of neutron-rich nuclear matter. The S$π$RIT TPC has a 50.5 cm drift length and an 86.4 cm $\times$ 134.4 cm pad plane with 12,096…
▽ More
The SAMURAI Pion Reconstruction and Ion-Tracker Time Projection Chamber (S$π$RIT TPC) was designed to enable measurements of heavy ion collisions with the SAMURAI spectrometer at the RIKEN Radioactive Isotope Beam Factory and provide constraints on the Equation of State of neutron-rich nuclear matter. The S$π$RIT TPC has a 50.5 cm drift length and an 86.4 cm $\times$ 134.4 cm pad plane with 12,096 pads that are equipped with the Generic Electronics for TPCs readout electronics. The S$π$RIT TPC allows excellent reconstruction of particles and provides isotopic resolution for pions and other light charged particles across a wide range of energy losses and momenta. Details of the S$π$RIT TPC are presented, along with discussion of the TPC performance based on cosmic ray and experimental data.
△ Less
Submitted 21 May, 2020;
originally announced May 2020.