-
Experiences Readying Applications for Exascale
Authors:
Paul T. Bauman,
Reuben D. Budiardja,
Dmytro Bykov,
Noel Chalmers,
Jacqueline Chen,
Nicholas Curtis,
Marc Day,
Markus Eisenbach,
Lucas Esclapez,
Alessandro Fanfarillo,
William Freitag,
Nicholas Frontiere,
Antigoni Georgiadou,
Joseph Glenski,
Kalyana Gottiparthi,
Marc T. Henry de Frahan,
Gustav R. Jansen,
Wayne Joubert,
Justin G. Lietz,
Jakub Kurzak,
Nicholas Malaya,
Bronson Messer,
Damon McDougall,
Paul Mullowney,
Stephen Nichols
, et al. (7 additional authors not shown)
Abstract:
The advent of exascale computing invites an assessment of existing best practices for developing application readiness on the world's largest supercomputers. This work details observations from the last four years in preparing scientific applications to run on the Oak Ridge Leadership Computing Facility's (OLCF) Frontier system. This paper addresses a range of topics in software including programm…
▽ More
The advent of exascale computing invites an assessment of existing best practices for developing application readiness on the world's largest supercomputers. This work details observations from the last four years in preparing scientific applications to run on the Oak Ridge Leadership Computing Facility's (OLCF) Frontier system. This paper addresses a range of topics in software including programmability, tuning, and portability considerations that are key to moving applications from existing systems to future installations. A set of representative workloads provides case studies for general system and software testing. We evaluate the use of early access systems for development across several generations of hardware. Finally, we discuss how best practices were identified and disseminated to the community through a wide range of activities including user-guides and trainings. We conclude with recommendations for ensuring application readiness on future leadership computing systems.
△ Less
Submitted 2 October, 2023;
originally announced October 2023.
-
Augmenting Tactile Simulators with Real-like and Zero-Shot Capabilities
Authors:
Osher Azulay,
Alon Mizrahi,
Nimrod Curtis,
Avishai Sintov
Abstract:
Simulating tactile perception could potentially leverage the learning capabilities of robotic systems in manipulation tasks. However, the reality gap of simulators for high-resolution tactile sensors remains large. Models trained on simulated data often fail in zero-shot inference and require fine-tuning with real data. In addition, work on high-resolution sensors commonly focus on ones with flat…
▽ More
Simulating tactile perception could potentially leverage the learning capabilities of robotic systems in manipulation tasks. However, the reality gap of simulators for high-resolution tactile sensors remains large. Models trained on simulated data often fail in zero-shot inference and require fine-tuning with real data. In addition, work on high-resolution sensors commonly focus on ones with flat surfaces while 3D round sensors are essential for dexterous manipulation. In this paper, we propose a bi-directional Generative Adversarial Network (GAN) termed SightGAN. SightGAN relies on the early CycleGAN while including two additional loss components aimed to accurately reconstruct background and contact patterns including small contact traces. The proposed SightGAN learns real-to-sim and sim-to-real processes over difference images. It is shown to generate real-like synthetic images while maintaining accurate contact positioning. The generated images can be used to train zero-shot models for newly fabricated sensors. Consequently, the resulted sim-to-real generator could be built on top of the tactile simulator to provide a real-world framework. Potentially, the framework can be used to train, for instance, reinforcement learning policies of manipulation tasks. The proposed model is verified in extensive experiments with test data collected from real sensors and also shown to maintain embedded force information within the tactile images.
△ Less
Submitted 19 September, 2023;
originally announced September 2023.
-
AllSight: A Low-Cost and High-Resolution Round Tactile Sensor with Zero-Shot Learning Capability
Authors:
Osher Azulay,
Nimrod Curtis,
Rotem Sokolovsky,
Guy Levitski,
Daniel Slomovik,
Guy Lilling,
Avishai Sintov
Abstract:
Tactile sensing is a necessary capability for a robotic hand to perform fine manipulations and interact with the environment. Optical sensors are a promising solution for high-resolution contact estimation. Nevertheless, they are usually not easy to fabricate and require individual calibration in order to acquire sufficient accuracy. In this letter, we propose AllSight, an optical tactile sensor w…
▽ More
Tactile sensing is a necessary capability for a robotic hand to perform fine manipulations and interact with the environment. Optical sensors are a promising solution for high-resolution contact estimation. Nevertheless, they are usually not easy to fabricate and require individual calibration in order to acquire sufficient accuracy. In this letter, we propose AllSight, an optical tactile sensor with a round 3D structure potentially designed for robotic in-hand manipulation tasks. AllSight is mostly 3D printed making it low-cost, modular, durable and in the size of a human thumb while with a large contact surface. We show the ability of AllSight to learn and estimate a full contact state, i.e., contact position, forces and torsion. With that, an experimental benchmark between various configurations of illumination and contact elastomers are provided. Furthermore, the robust design of AllSight provides it with a unique zero-shot capability such that a practitioner can fabricate the open-source design and have a ready-to-use state estimation model. A set of experiments demonstrates the accurate state estimation performance of AllSight.
△ Less
Submitted 11 November, 2023; v1 submitted 6 July, 2023;
originally announced July 2023.
-
Cross-shell states in $^{15}$C: a test for p-sd interactions
Authors:
J. Lois-Fuentes,
B. Fernández-Domínguez,
X. Pereira-López,
F. Delaunay,
W. N. Catford,
A. Matta,
N. A. Orr,
T. Duguet,
T. Otsuka,
V. Somà,
O. Sorlin,
T. Suzuki,
N. L. Achouri,
M. Assié,
S. Bailey,
B. Bastin,
Y. Blumenfeld,
R. Borcea,
M. Caamaño,
L. Caceres,
E. Clément,
A. Corsi,
N. Curtis,
Q. Deshayes,
F. Farget
, et al. (37 additional authors not shown)
Abstract:
The low-lying structure of $^{15}$C has been investigated via the neutron-removal $^{16}$C$(d,t)$ reaction. Along with bound neutron sd-shell hole states, unbound p-shell hole states have been firmly confirmed. The excitation energies and the deduced spectroscopic factors of the cross-shell states are an important measure of the $[(p)^{-1}(sd)^{2}]$ neutron configurations in $^{15}$C. Our results…
▽ More
The low-lying structure of $^{15}$C has been investigated via the neutron-removal $^{16}$C$(d,t)$ reaction. Along with bound neutron sd-shell hole states, unbound p-shell hole states have been firmly confirmed. The excitation energies and the deduced spectroscopic factors of the cross-shell states are an important measure of the $[(p)^{-1}(sd)^{2}]$ neutron configurations in $^{15}$C. Our results show a very good agreement with shell-model calculations using the SFO-tls interaction for $^{15}$C. However, a modification of the $p$-$sd$ and $sd$-$sd$ monopole terms was applied in order to reproduce the $N=9$ isotone $^{17}$O. In addition, the excitation energies and spectroscopic factors have been compared to the first calculations of $^{15}$C with the $ab~ initio$ self-consistent Green's function method employing the NNLO$_{sat}$ interaction. The results show the sensitivity to the size of the $N=8$ shell gap and highlight the need of going beyond the current truncation scheme in the theory.
△ Less
Submitted 16 February, 2023;
originally announced February 2023.
-
Accelerating reactive-flow simulations using vectorized chemistry integration
Authors:
Nicholas J. Curtis,
Kyle E. Niemeyer,
Chih-Jen Sung
Abstract:
The high cost of chemistry integration is a significant computational bottleneck for realistic reactive-flow simulations using operator splitting. Here we present a methodology to accelerate the solution of the chemical kinetic ordinary differential equations using single-instruction, multiple-data vector processing on CPUs using the OpenCL framework. First, we compared several vectorized integrat…
▽ More
The high cost of chemistry integration is a significant computational bottleneck for realistic reactive-flow simulations using operator splitting. Here we present a methodology to accelerate the solution of the chemical kinetic ordinary differential equations using single-instruction, multiple-data vector processing on CPUs using the OpenCL framework. First, we compared several vectorized integration algorithms using chemical kinetic source terms and analytical Jacobians from the pyJac software against a widely used integration code, CVODEs. Next, we extended the OpenFOAM computational fluid dynamics library to incorporate the vectorized solvers, and we compared the accuracy of a fourth-order linearly implicit integrator -- both in vectorized form and a corresponding method native to OpenFOAM -- with the community standard chemical kinetics library Cantera. We then applied our methodology to a variety of chemical kinetic models, turbulent intensities, and simulation scales to examine a range of engineering and scientific scale problems, including (pseudo) steady-state as well as time-dependent Reynolds-averaged Navier--Stokes simulations of the Sandia flame D and the Volvo Flygmotor bluff-body stabilized, premixed flame. Subsequently, we compared the performance of the vectorized and native OpenFOAM integrators over the studied models and simulations and found that our vectorized approach performs up to 33--35x faster than the native OpenFOAM solver with high accuracy.
△ Less
Submitted 10 May, 2022;
originally announced May 2022.
-
Low-lying single-particle structure of 17C and the N = 14 sub-shell closure
Authors:
X. Pereira-López,
B. Fernández-Domínguez,
F. Delaunay,
N. L. Achouri,
N. A. Orr,
W. N. Catford,
M. Assié,
S. Bailey,
B. Bastin,
Y. Blumenfeld,
R. Borcea,
M. Caamaño,
L. Caceres,
E. Clément,
A. Corsi,
N. Curtis,
Q. Deshayes,
F. Farget,
M. Fisichella,
G. de France,
S. Franchoo,
M. Freer,
J. Gibelin,
A. Gillibert,
G. F. Grinyer
, et al. (36 additional authors not shown)
Abstract:
The first investigation of the single-particle structure of the bound states of 17C, via the d(16C, p) transfer reaction, has been undertaken. The measured angular distributions confirm the spin-parity assignments of 1/2+ and 5/2+ for the excited states located at 217 and 335 keV, respectively. The spectroscopic factors deduced for these states exhibit a marked single-particle character, in agreem…
▽ More
The first investigation of the single-particle structure of the bound states of 17C, via the d(16C, p) transfer reaction, has been undertaken. The measured angular distributions confirm the spin-parity assignments of 1/2+ and 5/2+ for the excited states located at 217 and 335 keV, respectively. The spectroscopic factors deduced for these states exhibit a marked single-particle character, in agreement with shell model and particle-core model calculations, and combined with their near degeneracy in energy provide clear evidence for the absence of the N = 14 sub-shell closure. The very small spectroscopic factor found for the 3/2+ ground state is consistent with theoretical predictions and indicates that the ν1d3/2 strength is carried by unbound states. With a dominant l = 0 valence neutron configuration and a very low separation energy, the 1/2+ excited state is a one-neutron halo candidate.
△ Less
Submitted 11 November, 2020;
originally announced November 2020.
-
Clustering in $^{18}$O -- absolute determination of branching ratios via high-resolution particle spectroscopy
Authors:
S. Pirrie,
C. Wheldon,
Tz. Kokalova,
J. Bishop,
R. Hertenberger,
H. -F. Wirth,
S. Bailey,
N. Curtis,
D. Dell'Aquila,
Th. Faestermann,
D. Mengoni,
R. Smith,
D. Torresi,
A. Turner
Abstract:
The determination of absolute branching ratios for high-energy states in light nuclei is an important and useful tool for probing the underlying nuclear structure of individual resonances: for example, in establishing the tendency of an excited state towards $α$-cluster structure. Difficulty arises in measuring these branching ratios due to similarities in available decay channels, such as (…
▽ More
The determination of absolute branching ratios for high-energy states in light nuclei is an important and useful tool for probing the underlying nuclear structure of individual resonances: for example, in establishing the tendency of an excited state towards $α$-cluster structure. Difficulty arises in measuring these branching ratios due to similarities in available decay channels, such as ($\mathbf{^{18}}$O,$\mathbf{n}$) and ($\mathbf{^{18}}$O,$\mathbf{2n}$), as well as differences in geometric efficiencies due to population of bound excited levels in daughter nuclei. Methods are presented using Monte Carlo techniques to overcome these issues.
△ Less
Submitted 30 November, 2019;
originally announced December 2019.
-
Experimental investigation of $α$-condensation in light nuclei
Authors:
Jack Bishop,
Tzany Kokalova,
Martin Freer,
L Acosta,
M Assie,
S Bailey,
G Cardella,
N Curtis,
E De Filippo,
D Dell'Aquila,
S De Luca,
L Francalanza,
B Gnoffo,
G Lanzalone,
I Lombardo,
N. S. Martorana,
S Norella,
A Pagano,
E. V. Pagano,
M. Papa,
S. Pirrone,
G Politi,
F Rizzo,
P Russotto,
L Quattrocchi
, et al. (7 additional authors not shown)
Abstract:
Method: To examine signatures of this alpha-condensation, a compound nucleus reaction using 160, 280, and 400 MeV 16O beams impinging on a carbon target was used to investigate the 12C(16O,7a) reaction. This permits a search for near-threshold states in the alpha-conjugate nuclei up to 24Mg. Results: Events up to an alpha-particle multiplicity of 7 were measured and the results were compared to bo…
▽ More
Method: To examine signatures of this alpha-condensation, a compound nucleus reaction using 160, 280, and 400 MeV 16O beams impinging on a carbon target was used to investigate the 12C(16O,7a) reaction. This permits a search for near-threshold states in the alpha-conjugate nuclei up to 24Mg. Results: Events up to an alpha-particle multiplicity of 7 were measured and the results were compared to both an Extended Hauser-Feshbach calculation and the Fermi break-up model. The measured multiplicity distribution exceeded that predicted from a sequential decay mechanism and had a better agreement with the multi-particle Fermi break-up model. Examination of how these 7 alpha final states could be reconstructed to form 8Be and 12C(0_2+) showed a quantitative difference in which decay modes were dominant compared to the Fermi break-up model. No new states were observed in 16O, 20Ne, and 24Mg due to the effect of the N-alpha penetrability suppressing the total alpha-particle dissociation decay mode. Conclusion: The reaction mechanism for a high energy compound nucleus reaction can only be described by a hybrid of sequential decay and multi-particle breakup. Highly alpha-clustered states were seen which did not originate from simple binary reaction processes. Direct investigations of near-threshold states in N-alpha systems are inherently impeded by the Coulomb barrier prohibiting the observation of states in the N-alpha decay channel. No evidence of a highly clustered 15.1 MeV state in 16O was observed from (28Si*,12C(0_2+))16O(0_6+) when reconstructing the Hoyle state from 3 alpha-particles. Therefore, no experimental signatures for alpha-condensation were observed.
△ Less
Submitted 19 August, 2019; v1 submitted 11 July, 2019;
originally announced July 2019.
-
Using SIMD and SIMT vectorization to evaluate sparse chemical kinetic Jacobian matrices and thermochemical source terms
Authors:
Nicholas J. Curtis,
Kyle E. Niemeyer,
Chih-Jen Sung
Abstract:
Accurately predicting key combustion phenomena in reactive-flow simulations, e.g., lean blow-out, extinction/ignition limits and pollutant formation, necessitates the use of detailed chemical kinetics. The large size and high levels of numerical stiffness typically present in chemical kinetic models relevant to transportation/power-generation applications make the efficient evaluation/factorizatio…
▽ More
Accurately predicting key combustion phenomena in reactive-flow simulations, e.g., lean blow-out, extinction/ignition limits and pollutant formation, necessitates the use of detailed chemical kinetics. The large size and high levels of numerical stiffness typically present in chemical kinetic models relevant to transportation/power-generation applications make the efficient evaluation/factorization of the chemical kinetic Jacobian and thermochemical source-terms critical to the performance of reactive-flow codes. Here we investigate the performance of vectorized evaluation of constant-pressure/volume thermochemical source-term and sparse/dense chemical kinetic Jacobians using single-instruction, multiple-data (SIMD) and single-instruction, multiple thread (SIMT) paradigms. These are implemented in pyJac, an open-source, reproducible code generation platform. A new formulation of the chemical kinetic governing equations was derived and verified, resulting in Jacobian sparsities of 28.6-92.0% for the tested models. Speedups of 3.40-4.08x were found for shallow-vectorized OpenCL source-rate evaluation compared with a parallel OpenMP code on an avx2 central processing unit (CPU), increasing to 6.63-9.44x and 3.03-4.23x for sparse and dense chemical kinetic Jacobian evaluation, respectively. Furthermore, the effect of data-ordering was investigated and a storage pattern specifically formulated for vectorized evaluation was proposed; as well, the effect of the constant pressure/volume assumptions and varying vector widths were studied on source-term evaluation performance. Speedups reached up to 17.60x and 45.13x for dense and sparse evaluation on the GPU, and up to 55.11x and 245.63x on the CPU over a first-order finite-difference Jacobian approach. Further, dense Jacobian evaluation was up to 19.56x and 2.84x times faster than a previous version of pyJac on a CPU and GPU, respectively.
△ Less
Submitted 4 September, 2018;
originally announced September 2018.
-
An automated target species selection method for dynamic adaptive chemistry simulations
Authors:
Nicholas Curtis,
Kyle Niemeyer,
Chih-Jen Sung
Abstract:
The relative importance index (RII) method for determining appropriate target species for dynamic adaptive chemistry (DAC) simulations using the DRGEP method is developed. The accuracy and effectiveness of this RII method is validated for two fuels: n-heptane and isopentanol. The RII method determines appropriate DRGEP target species solely from the local thermo-chemical state of the simulation, e…
▽ More
The relative importance index (RII) method for determining appropriate target species for dynamic adaptive chemistry (DAC) simulations using the DRGEP method is developed. The accuracy and effectiveness of this RII method is validated for two fuels: n-heptane and isopentanol. The RII method determines appropriate DRGEP target species solely from the local thermo-chemical state of the simulation, ensuring that accuracy will be maintained. Further, the RII method reduces the expertise required of users to select DRGEP target species sets appropriate to the combustion phenomena under consideration. Constant volume autoignition simulations run over a wide range of initial condi- tions using detailed reaction mechanisms for n-heptane and isopentanol show that the RII method is able to maintain accuracy even when traditional static target species sets fail, and are even more accurate than expert-selected target species sets. Additionally, the accuracy and efficiency of the RII method are compared to those of static target species sets in single-cell engine simulations under homogeneous charge compression ignition conditions. For simulations using more stringent DRGEP thresholds, the RII method performs similarly to that of the static target species sets. With a larger DRGEP threshold, the RII method is significantly more accurate than the static target species sets without imposing significant computational overhead. Furthermore, the applicability of the RII method to a DRG-based DAC scheme is discussed.
△ Less
Submitted 5 March, 2018;
originally announced April 2018.
-
Chronology of the three-body dissociation of $^8$He
Authors:
B. Laurent,
F. M. Marqués,
C. Angulo,
N. I. Ashwood,
M. J. G. Borge,
V. Bouchat,
W. N. Catford,
N. M. Clarke,
N. Curtis,
M. Freer,
F. Hanappe,
V. Kinnard,
M. Labiche,
T. Materna,
P. McEwan,
T. Nilsson,
A. Ninane,
G. Normand,
N. A. Orr,
S. D. Pain,
E. Prokhorova,
L. Stuttgé,
C. Timis
Abstract:
The space and time configurations of the dissociation of $^8$He into $^6$He+$n$+$n$, on C and Pb targets, have been explored simultaneously for the first time. The final-state interactions in the $n$-$n$ and $^6$He-$n$ channels are successfully described within a model that considers independent emission of neutrons from a Gaussian volume with a given lifetime. The dissociation on C target exhibit…
▽ More
The space and time configurations of the dissociation of $^8$He into $^6$He+$n$+$n$, on C and Pb targets, have been explored simultaneously for the first time. The final-state interactions in the $n$-$n$ and $^6$He-$n$ channels are successfully described within a model that considers independent emission of neutrons from a Gaussian volume with a given lifetime. The dissociation on C target exhibits a dominant sequential decay through the ground state of $^7$He, consistent with neutrons being emitted from a Gaussian volume of $r_{nn}^{\rm{rms}}=7.3\pm0.6$~fm with a $n$-$n$ delay in the sequential channel of $1400\pm400$~fm/$c$, in agreement with the lifetime of $^7$He. The lower-statistics data on Pb target correspond mainly to direct breakup, and are well described using the $n$-$n$ volume measured, without any $n$-$n$ delay. The validity of the phenomenological model used is discussed.
△ Less
Submitted 12 February, 2019; v1 submitted 4 January, 2018;
originally announced January 2018.
-
An investigation of GPU-based stiff chemical kinetics integration methods
Authors:
Nicholas J. Curtis,
Kyle E. Niemeyer,
Chih-Jen Sung
Abstract:
A fifth-order implicit Runge-Kutta method and two fourth-order exponential integration methods equipped with Krylov subspace approximations were implemented for the GPU and paired with the analytical chemical kinetic Jacobian software pyJac. The performance of each algorithm was evaluated by integrating thermochemical state data sampled from stochastic partially stirred reactor simulations and com…
▽ More
A fifth-order implicit Runge-Kutta method and two fourth-order exponential integration methods equipped with Krylov subspace approximations were implemented for the GPU and paired with the analytical chemical kinetic Jacobian software pyJac. The performance of each algorithm was evaluated by integrating thermochemical state data sampled from stochastic partially stirred reactor simulations and compared with the commonly used CPU-based implicit integrator CVODE. We estimated that the implicit Runge-Kutta method running on a single GPU is equivalent to CVODE running on 12-38 CPU cores for integration of a single global integration time step of 1e-6 s with hydrogen and methane models. In the stiffest case studied---the methane model with a global integration time step of 1e-4 s---thread divergence and higher memory traffic significantly decreased GPU performance to the equivalent of CVODE running on approximately three CPU cores. The exponential integration algorithms performed more slowly than the implicit integrators on both the CPU and GPU. Thread divergence and memory traffic were identified as the main limiters of GPU integrator performance, and techniques to mitigate these issues were discussed. Use of a finite-difference Jacobian on the GPU---in place of the analytical Jacobian provided by pyJac---greatly decreased integrator performance due to thread divergence, resulting in maximum slowdowns of 7.11-240.96 times; in comparison, the corresponding slowdowns on the CPU were just 1.39-2.61 times, underscoring the importance of use of an analytical Jacobian for efficient GPU integration. Finally, future research directions for working towards enabling realistic chemistry in reactive-flow simulations via GPU\slash SIMD accelerated stiff chemical kinetic integration were identified.
△ Less
Submitted 14 February, 2017; v1 submitted 13 July, 2016;
originally announced July 2016.
-
pyJac: analytical Jacobian generator for chemical kinetics
Authors:
Kyle E. Niemeyer,
Nicholas J. Curtis,
Chih-Jen Sung
Abstract:
Accurate simulations of combustion phenomena require the use of detailed chemical kinetics in order to capture limit phenomena such as ignition and extinction as well as predict pollutant formation. However, the chemical kinetic models for hydrocarbon fuels of practical interest typically have large numbers of species and reactions and exhibit high levels of mathematical stiffness in the governing…
▽ More
Accurate simulations of combustion phenomena require the use of detailed chemical kinetics in order to capture limit phenomena such as ignition and extinction as well as predict pollutant formation. However, the chemical kinetic models for hydrocarbon fuels of practical interest typically have large numbers of species and reactions and exhibit high levels of mathematical stiffness in the governing differential equations, particularly for larger fuel molecules. In order to integrate the stiff equations governing chemical kinetics, generally reactive-flow simulations rely on implicit algorithms that require frequent Jacobian matrix evaluations. Some in situ and a posteriori computational diagnostics methods also require accurate Jacobian matrices, including computational singular perturbation and chemical explosive mode analysis. Typically, finite differences numerically approximate these, but for larger chemical kinetic models this poses significant computational demands since the number of chemical source term evaluations scales with the square of species count. Furthermore, existing analytical Jacobian tools do not optimize evaluations or support emerging SIMD processors such as GPUs. Here we introduce pyJac, a Python-based open-source program that generates analytical Jacobian matrices for use in chemical kinetics modeling and analysis. As a demonstration, we first establish the correctness of the Jacobian matrices for kinetic models of hydrogen, methane, ethylene, and isopentanol oxidation, then demonstrate the performance achievable on CPUs and GPUs using pyJac via matrix evaluation timing comparisons.
△ Less
Submitted 19 February, 2017; v1 submitted 10 May, 2016;
originally announced May 2016.
-
Structure of $^{13}$Be probed via secondary beam reactions
Authors:
G. Randisi,
A. Leprince,
H. Al Falou,
N. A. Orr,
F. M. Marqués,
N. L. Achouri,
J. -C. Angélique,
N. Ashwood,
B. Bastin,
T. Bloxham,
B. A. Brown,
W. N. Catford,
N. Curtis,
F. Delaunay,
M. Freer,
E. de Góes Brennand,
P. Haigh,
F. Hanappe,
C. Harlin,
B. Laurent,
J. -L. Lecouey,
A. Ninane,
N. Patterson,
D. Price,
L. Stuttgé
, et al. (1 additional authors not shown)
Abstract:
The low-lying level structure of the unbound neutron-rich nucleus $^{13}$Be has been investigated via breakup on a carbon target of secondary beams of $^{14,15}$B at 35 MeV/nucleon. The coincident detection of the beam velocity $^{12}$Be fragments and neutrons permitted the invariant mass of the $^{12}$Be+$n$ and $^{12}$Be+$n$+$n$ systems to be reconstructed. In the case of the breakup of $^{15}$B…
▽ More
The low-lying level structure of the unbound neutron-rich nucleus $^{13}$Be has been investigated via breakup on a carbon target of secondary beams of $^{14,15}$B at 35 MeV/nucleon. The coincident detection of the beam velocity $^{12}$Be fragments and neutrons permitted the invariant mass of the $^{12}$Be+$n$ and $^{12}$Be+$n$+$n$ systems to be reconstructed. In the case of the breakup of $^{15}$B, a very narrow structure at threshold was observed in the $^{12}$Be+$n$ channel. Contrary to earlier stable beam fragmentation studies which identified this as a strongly interacting $s$-wave virtual state in $^{13}$Be, analysis here of the $^{12}$Be+$n$+$n$ events demonstrated that this was an artifact resulting from the sequential-decay of the $^{14}$Be(2$^+$) state. Single-proton removal from $^{14}$B was found to populate a broad low-lying structure some 0.70 MeV above the neutron-decay threshold in addition to a less prominent feature at around 2.4 MeV. Based on the selectivity of the reaction and a comparison with (0-3)$\hbarω$ shell-model calculations, the low-lying structure is concluded to most probably arise from closely spaced J$^π$=1/2$^+$ and 5/2$^+$ resonances (E$_r$=0.40$\pm$0.03 and 0.85$^{+0.15}_{-0.11}$ MeV), whilst the broad higher-lying feature is a second 5/2$^+$ level (E$_r$=2.35$\pm$0.14 MeV). Taken in conjunction with earlier studies, it would appear that the lowest 1/2$^+$ and 1/2$^-$ levels lie relatively close together below 1 MeV.
△ Less
Submitted 28 February, 2014; v1 submitted 9 December, 2013;
originally announced December 2013.
-
Emergence of the N=16 shell gap in 21O
Authors:
B. Fernandez-Dominguez,
J. S. Thomas,
W. N. Catford,
F. Delaunay,
S. M. Brown,
N. A. Orr,
M. Rejmund,
N. L. Achouri,
H. Al Falou,
N. A. Ashwood,
D. Beaumel,
Y. Blumenfeld,
B. A. Brown,
R. Chapman,
M. Chartier,
N. Curtis,
C. Force,
G. de France,
S. Franchoo,
J. Guillot,
P. Haigh,
F. Hammache,
M. Labiche,
V. Lapoux,
R. C. Lemmon
, et al. (17 additional authors not shown)
Abstract:
The spectroscopy of 21O has been investigated using a radioactive 20O beam and the (d,p) reaction in inverse kinematics. The ground and first excited states have been determined to be Jpi=5/2+ and Jpi=1/2+ respectively. Two neutron unbound states were observed at excitation energies of 4.76 +- 0.10 and 6.16 +- 0.11. The spectroscopic factor deduced for the lower of these interpreted as a 3/2+ leve…
▽ More
The spectroscopy of 21O has been investigated using a radioactive 20O beam and the (d,p) reaction in inverse kinematics. The ground and first excited states have been determined to be Jpi=5/2+ and Jpi=1/2+ respectively. Two neutron unbound states were observed at excitation energies of 4.76 +- 0.10 and 6.16 +- 0.11. The spectroscopic factor deduced for the lower of these interpreted as a 3/2+ level, reveals a rather pure 0d3/2 single-particle configuration. The large energy difference between the 3/2+ and 1/2+ states is indicative of the emergence of the N=16 magic number. For the higher lying resonance, which has a character consistent with a spin-parity assignment of 3/2+ or 7/2-, a 71% branching ratio to the first 2+ state in 20O has been observed. The results are compared with new shell model calculations.
△ Less
Submitted 17 December, 2010;
originally announced December 2010.
-
Migration of nuclear shell gaps studied in the d(24Ne,p gamma)25Ne reaction
Authors:
W. N. Catford,
C. N. Timis,
R. C. Lemmon,
M. Labiche,
N. A. Orr,
B. Fernandez-Dominguez,
R. Chapman,
M. Freer,
M. Chartier,
H. Savajols,
M. Rejmund,
N. L. Achouri,
N. Amzal,
N. I. Ashwood,
T. D. Baldwin,
M. Burns,
L. Caballero,
J. M. Casadjian,
N. Curtis,
G. de France,
W. Gelletly,
X. Liang,
S. D. Pain,
V. P. E. Pucknell,
B. Rubio
, et al. (4 additional authors not shown)
Abstract:
The transfer of neutrons onto 24Ne has been measured using a reaccelerated radioactive beam of 24Ne to study the (d,p) reaction in inverse kinematics. The unusual raising of the first 3/2+ level in 25Ne and its significance in terms of the migration of the neutron magic number from N=20 to N=16 is put on a firm footing by confirmation of this state's identity. The raised 3/2+ level is observed sim…
▽ More
The transfer of neutrons onto 24Ne has been measured using a reaccelerated radioactive beam of 24Ne to study the (d,p) reaction in inverse kinematics. The unusual raising of the first 3/2+ level in 25Ne and its significance in terms of the migration of the neutron magic number from N=20 to N=16 is put on a firm footing by confirmation of this state's identity. The raised 3/2+ level is observed simultaneously with the intruder negative parity 7/2- and 3/2- levels, providing evidence for the reduction in the N=20 gap. The coincident gamma-ray decays allowed the assignment of spins as well as the transferred orbital angular momentum. The excitation energy of the 3/2+ state shows that the established USD shell model breaks down well within the sd model space and requires a revised treatment of the proton-neutron monopole interaction.
△ Less
Submitted 9 April, 2010;
originally announced April 2010.
-
TIARA: a large solid angle silicon array for direct reaction studies with radioactive beams
Authors:
M. Labiche,
W. N. Catford,
R. C. Lemmon,
C. N. Timis,
R. Chapman,
N. A. Orr,
B. Fernandez-Dominguez,
G. Moores,
N. L. Achouri,
N. Amzal,
S. Appleton,
N. I. Ashwood,
T. D. Baldwin,
M. Burns,
L. Caballero,
J. Cacitti,
J. M. Casadjian,
M. Chartier,
N. Curtis,
K. Faiz,
G. de France,
M. Freer,
J. M. Gautier,
W. Gelletly,
G. Iltis
, et al. (17 additional authors not shown)
Abstract:
A compact, quasi-4pi position sensitive silicon array, TIARA, designed to study direct reactions induced by radioactive beams in inverse kinematics is described here. The Transfer and Inelastic All-angle Reaction Array (TIARA) consists of 8 resistive charge division detectors forming an octagonal barrel around the target and a set of double-sided silicon-strip annular detectors positioned at eac…
▽ More
A compact, quasi-4pi position sensitive silicon array, TIARA, designed to study direct reactions induced by radioactive beams in inverse kinematics is described here. The Transfer and Inelastic All-angle Reaction Array (TIARA) consists of 8 resistive charge division detectors forming an octagonal barrel around the target and a set of double-sided silicon-strip annular detectors positioned at each end of the barrel. The detector was coupled to the -ray array EXOGAM and the spectrometer VAMOS at the GANIL Laboratory to demonstrate the potential of such an apparatus with radioactive beams. The 14N(d,p)15N reaction, well known in direct kinematics, has been carried out in inverse kinematics for that purpose. The observation of the 15N ground state and excited states at 7.16 and 7.86 MeV is presented here as well as the comparison of the measured proton angular distributions with DWBA calculations. Transferred l-values are in very good agreement with both theoretical calculations and previous experimental results obtained in direct kinematics.
△ Less
Submitted 13 January, 2010; v1 submitted 20 August, 2009;
originally announced August 2009.
-
B(E1) Strengths from Coulomb Excitation of 11Be
Authors:
N. C. Summers,
S. D. Pain,
N. A. Orr,
W. N. Catford,
J. C. Angelique,
N. I. Ashwood,
V. Bouchat,
N. M. Clarke,
N. Curtis,
M. Freer,
B. R. Fulton,
F. Hanappe,
M. Labiche,
J. L. Lecouey,
R. C. Lemmon,
D. Mahboub,
A. Ninane,
G. Normand,
F. M. Nunes,
N. Soic,
L. Stuttge,
C. N. Timis,
I. J. Thompson,
J. S. Winfield,
V. Ziman
Abstract:
The $B$(E1;$1/2^+\to1/2^-$) strength for $^{11}$Be has been extracted from intermediate energy Coulomb excitation measurements, over a range of beam energies using a new reaction model, the extended continuum discretized coupled channels (XCDCC) method. In addition, a measurement of the excitation cross section for $^{11}$Be+$^{208}$Pb at 38.6 MeV/nucleon is reported. The $B$(E1) strength of 0.1…
▽ More
The $B$(E1;$1/2^+\to1/2^-$) strength for $^{11}$Be has been extracted from intermediate energy Coulomb excitation measurements, over a range of beam energies using a new reaction model, the extended continuum discretized coupled channels (XCDCC) method. In addition, a measurement of the excitation cross section for $^{11}$Be+$^{208}$Pb at 38.6 MeV/nucleon is reported. The $B$(E1) strength of 0.105(12) e$^2$fm$^2$ derived from this measurement is consistent with those made previously at 60 and 64 MeV/nucleon, i n contrast to an anomalously low result obtained at 43 MeV/nucleon. By coupling a multi-configuration description of the projectile structure with realistic reaction theory, the XCDCC model provides for the first time a fully quantum mechanical description of Coulomb excitation. The XCDCC calculations reveal that the excitation process involves significant contributions from nuclear, continuum, and higher-order effects. An analysis of the present and two earlier intermediate energy measurements yields a combined B(E1) strength of 0.105(7) e$^2$fm$^2$. This value is in good agreement with the value deduced independently from the lifetime of the $1/2^-$ state in $^{11}$Be, and has a comparable p recision.
△ Less
Submitted 10 May, 2007; v1 submitted 16 March, 2007;
originally announced March 2007.
-
Structure of 12Be: intruder d-wave strength at N=8
Authors:
S. D. Pain,
W. N. Catford,
N. A. Orr,
J. C. Angelique,
N. I. Ashwood,
V. Bouchat,
N. M. Clarke,
N. Curtis,
M. Freer,
B. R. Fulton,
F. Hanappe,
M. Labiche,
J. L. Lecouey,
R. C. Lemmon,
D. Mahboub,
A. Ninane,
G. Normand,
N. Soic,
L. Stuttge,
C. N. Timis,
J. A. Tostevin,
J. S. Winfield,
V. Ziman,
.
Abstract:
The breaking of the N=8 shell-model magic number in the 12Be ground state has been determined to include significant occupancy of the intruder d-wave orbital. This is in marked contrast with all other N=8 isotones, both more and less exotic than 12Be. The occupancies of the 0 hbar omega neutron p1/2-orbital and the 1 hbar omega, neutron d5/2 intruder orbital were deduced from a measurement of ne…
▽ More
The breaking of the N=8 shell-model magic number in the 12Be ground state has been determined to include significant occupancy of the intruder d-wave orbital. This is in marked contrast with all other N=8 isotones, both more and less exotic than 12Be. The occupancies of the 0 hbar omega neutron p1/2-orbital and the 1 hbar omega, neutron d5/2 intruder orbital were deduced from a measurement of neutron removal from a high-energy 12Be beam leading to bound and unbound states in 11Be.
△ Less
Submitted 16 October, 2005;
originally announced October 2005.
-
Single-neutron transfer from 11Be gs via the (p,d) reaction with a radioactive beam
Authors:
J. S. Winfield,
S. Fortier,
W. N. Catford,
S. Pita,
N. A. Orr,
J. Van de Wiele,
Y. Blumenfeld,
R. Chapman,
S. P. G. Chappell,
N. M. Clarke,
N. Curtis,
M. Freer,
S. Gales,
H. Langevin-Joliot,
H. Laurent,
I. Lhenry,
J. M. Maison,
P. Roussel-Chomaz,
M. Shawcross,
K. Spohr,
T. Suomijarvi,
A. de Vismes
Abstract:
The 11Be(p,d)10Be reaction has been performed in inverse kinematics with a radioactive 11Be beam of E/A = 35.3 MeV. Angular distributions for the 0+ ground state, the 2+, 3.37 MeV state and the multiplet of states around 6 MeV in 10Be were measured at angles up to 16 deg CM by detecting the 10Be in a dispersion-matched spectrometer and the coincident deuterons in a silicon array. Distorted wave…
▽ More
The 11Be(p,d)10Be reaction has been performed in inverse kinematics with a radioactive 11Be beam of E/A = 35.3 MeV. Angular distributions for the 0+ ground state, the 2+, 3.37 MeV state and the multiplet of states around 6 MeV in 10Be were measured at angles up to 16 deg CM by detecting the 10Be in a dispersion-matched spectrometer and the coincident deuterons in a silicon array. Distorted wave and coupled-channels calculations have been performed to investigate the amount of 2+ core excitation in 11Be gs. The use of "realistic" 11Be wave functions is emphasised and bound state form factors have been obtained by solving the particle-vibration coupling equations. This calculation gives a dominant 2s component in the 11Be gs wave function with a 16% [2+ x 1d] core excitation admixture. Cross sections calculated with these form factors are in good agreement with the present data. The Separation Energy prescription for the bound state wave function also gives satisfactory fits to the data, but leads to a significantly larger [2 x 1d] component in 11Be gs.
△ Less
Submitted 2 October, 2000; v1 submitted 20 September, 2000;
originally announced September 2000.
-
One-neutron removal reactions on neutron-rich psd-shell nuclei
Authors:
E. Sauvan,
F. Carstoiu,
N. A. Orr,
J. C. Angelique,
W. N. Catford,
N. M. Clarke,
M. Mac Cormick,
N. Curtis,
M. Freer,
S. Grevy,
C. Le Brun,
M. Lewitowicz,
E. Liegard,
F. M. Marques,
P. Roussel-Chomaz,
M. G. Saint Laurent,
M. Shawcross,
J. S. Winfield
Abstract:
A systematic study of high energy, one-neutron removal reactions on 23 neutron-rich, psd--shell nuclei (Z=5-9, A=12-25) has been carried out. The longitudinal momentum distributions of the core fragments and corresponding single-neutron removal cross sections are reported for reactions on a carbon target. Extended Glauber model calculations, weighted by the spectroscopic factors obtained from sh…
▽ More
A systematic study of high energy, one-neutron removal reactions on 23 neutron-rich, psd--shell nuclei (Z=5-9, A=12-25) has been carried out. The longitudinal momentum distributions of the core fragments and corresponding single-neutron removal cross sections are reported for reactions on a carbon target. Extended Glauber model calculations, weighted by the spectroscopic factors obtained from shell model calculations, are compared to the experimental results. Conclusions are drawn regarding the use of such reactions as a spectroscopic tool and spin-parity assignments are proposed for 15B, 17C, 19-21N, 21,23O, 23-25F. The nature of the weakly bound systems 14B and 15,17C is discussed.
△ Less
Submitted 11 July, 2000;
originally announced July 2000.
-
The Halo of 14Be
Authors:
M. Labiche,
N. A. Orr,
F. M. Marques,
J. C. Angelique,
L. Axelsson,
B. Benoit,
U. C. Bergmann,
M. J. G. Borge,
W. N. Catford,
S. P. G. Chappell,
N. M. Clarke,
G. Costa,
N. Curtis,
A. D'Arrigo,
E. de Goes Brennand,
O. Dorvaux,
G. Fazio,
M. Freer,
B. R. Fulton,
G. Giardina,
S. Grevy,
D. Guillemaud-Mueller,
F. Hanappe,
B. Heusch K. L. Jones,
B. Jonson
, et al. (17 additional authors not shown)
Abstract:
The two-neutron halo nucleus 14Be has been investigated in a kinematically complete measurement of the fragments (12Be and neutrons) produced in dissociation at 35 MeV/nucleon on C and Pb targets. Two-neutron removal cross-sections, neutron angular distributions and invariant mass spectra characteristic of a halo were observed and the electromagnetic (EMD) contributions deduced. Comparison with…
▽ More
The two-neutron halo nucleus 14Be has been investigated in a kinematically complete measurement of the fragments (12Be and neutrons) produced in dissociation at 35 MeV/nucleon on C and Pb targets. Two-neutron removal cross-sections, neutron angular distributions and invariant mass spectra characteristic of a halo were observed and the electromagnetic (EMD) contributions deduced. Comparison with three-body model predictions indicate that the halo wavefunction contains a large 2s1/2^2 admixture. The EMD invariant mass spectrum exhibited a relatively narrow structure near threshold (Edecay=1.8+/-0.1 MeV, Gamma = 0.8+/-0.4 MeV) consistent with a soft-dipole excitation.
△ Less
Submitted 7 June, 2000;
originally announced June 2000.
