-
A Bayesian mixture model approach to quantifying the empirical nuclear saturation point
Authors:
C. Drischler,
P. G. Giuliani,
S. Bezoui,
J. Piekarewicz,
F. Viens
Abstract:
The equation of state (EOS) in the limit of infinite symmetric nuclear matter exhibits an equilibrium density, $n_0 \approx 0.16 \, \mathrm{fm}^{-3}$, at which the pressure vanishes and the energy per particle attains its minimum, $E_0 \approx -16 \, \mathrm{MeV}$. Although not directly measurable, the nuclear saturation point $(n_0,E_0)$ can be extrapolated by density functional theory (DFT), pro…
▽ More
The equation of state (EOS) in the limit of infinite symmetric nuclear matter exhibits an equilibrium density, $n_0 \approx 0.16 \, \mathrm{fm}^{-3}$, at which the pressure vanishes and the energy per particle attains its minimum, $E_0 \approx -16 \, \mathrm{MeV}$. Although not directly measurable, the nuclear saturation point $(n_0,E_0)$ can be extrapolated by density functional theory (DFT), providing tight constraints for microscopic interactions derived from chiral effective field theory (EFT). However, when considering several DFT predictions for $(n_0,E_0)$ together, a discrepancy between model classes emerges at high confidence levels that each model prediction's uncertainty cannot explain. How can we leverage these DFT constraints to rigorously benchmark nuclear saturation properties of chiral interactions? To address this question, we present a Bayesian mixture model that combines multiple DFT predictions for $(n_0,E_0)$ using an efficient conjugate prior approach. The inferred posterior distribution for the saturation point's mean and covariance matrix follows a Normal-inverse-Wishart class, resulting in posterior predictives in the form of correlated, bivariate $t$-distributions. The DFT uncertainty reports are then used to mix these posteriors using an ordinary Monte Carlo approach. At the 95% credibility level, we estimate $n_0 \approx 0.157 \pm 0.010 \, \mathrm{fm}^{-3}$ and $E_0 \approx -15.97 \pm 0.40 \, \mathrm{MeV}$ for the marginal (univariate) $t$-distributions. Combined with chiral EFT calculations of the pure neutron matter EOS, we obtain bivariate normal distributions for the nuclear symmetry energy and its slope parameter evaluated at $n_0$: $S_v \approx 32.0 \pm 1.1 \, \mathrm{MeV}$ and $L\approx 52.6\pm 8.1 \, \mathrm{MeV}$ (95%), respectively. Our Bayesian framework is publicly available, so practitioners can readily use and extend our results.
△ Less
Submitted 4 May, 2024;
originally announced May 2024.
-
Toward emulating nuclear reactions using eigenvector continuation
Authors:
C. Drischler,
M. Quinonez,
P. G. Giuliani,
A. E. Lovell,
F. M. Nunes
Abstract:
We construct an efficient emulator for two-body scattering observables using the general (complex) Kohn variational principle and trial wave functions derived from eigenvector continuation. The emulator simultaneously evaluates an array of Kohn variational principles associated with different boundary conditions, which allows for the detection and removal of spurious singularities known as Kohn an…
▽ More
We construct an efficient emulator for two-body scattering observables using the general (complex) Kohn variational principle and trial wave functions derived from eigenvector continuation. The emulator simultaneously evaluates an array of Kohn variational principles associated with different boundary conditions, which allows for the detection and removal of spurious singularities known as Kohn anomalies. When applied to the $K$-matrix only, our emulator resembles the one constructed by Furnstahl et al. [Phys. Lett. B 809, 135719] although with reduced numerical noise. After a few applications to real potentials, we emulate differential cross sections for $^{40}$Ca$(n,n)$ scattering based on a realistic optical potential and quantify the model uncertainties using Bayesian methods. These calculations serve as a proof of principle for future studies aimed at improving optical models.
△ Less
Submitted 19 November, 2021; v1 submitted 18 August, 2021;
originally announced August 2021.
-
A novel experimental setup for rare events selection and its potential application to super heavy elements search
Authors:
Z. Majka,
R. Planeta,
Z. Sosin,
A. Wieloch,
K. Zelga,
M. Adamczyk,
K. Pelczar,
M. Barbui,
S. Wuenschel,
K. Hagel,
X. Cao,
E-J. Kim,
J. Natowitz,
R. Wada,
H. Zheng,
G. Giuliani,
S. Kowalski
Abstract:
The paper presents a novel instrumentation for rare events selection which was tested in our research of short lived super heavy elements production and detection. The instrumentation includes an active catcher multi elements system and dedicated electronics. The active catcher located in the forward hemisphere is composed of 63 scintillator detection modules. Reaction products of damped collision…
▽ More
The paper presents a novel instrumentation for rare events selection which was tested in our research of short lived super heavy elements production and detection. The instrumentation includes an active catcher multi elements system and dedicated electronics. The active catcher located in the forward hemisphere is composed of 63 scintillator detection modules. Reaction products of damped collisions between heavy ion projectiles and heavy target nuclei are implanted in the fast plastic scintillators of the active catcher modules. The acquisition system trigger delivered by logical branch of the electronics allows to record the reaction products which decay via the alpha particle emissions or spontaneous fission which take place between beam bursts. One microsecond wave form signal from FADCs contains information on heavy implanted nucleus as well as its decays.
△ Less
Submitted 17 October, 2018;
originally announced October 2018.
-
Evidence for high excitation energy resonances in the 7 alpha disassembly of $^{28}$Si
Authors:
X. G. Cao,
E. J. Kim,
K. Schmidt,
K. Hagel,
M. Barbui,
J. Gauthier,
S. Wuenschel,
G. Giuliani,
M. R. D. Rodriguez,
S. Kowalski,
H. Zheng,
M. Huang,
A. Bonasera,
R. Wada,
G. Q. Zhang,
C. Y. Wong,
A. Staszczak,
Z. X. Ren,
Y. K. Wang,
S. Q. Zhang,
J. Meng,
J. B. Natowitz
Abstract:
The excitation function for the 7 alpha de-excitation of $^{28}$Si nuclei excited to high excitation energies in the collisions of 35 MeV/nucleon $^{28}$Si with $^{12}$C reveals resonance structures that may indicate the population of high spin toroidal isomers such as those predicted by a number of recent theoretical calculations. This interpretation is supported by extended theoretical analyses.
The excitation function for the 7 alpha de-excitation of $^{28}$Si nuclei excited to high excitation energies in the collisions of 35 MeV/nucleon $^{28}$Si with $^{12}$C reveals resonance structures that may indicate the population of high spin toroidal isomers such as those predicted by a number of recent theoretical calculations. This interpretation is supported by extended theoretical analyses.
△ Less
Submitted 25 April, 2018; v1 submitted 22 January, 2018;
originally announced January 2018.
-
Alpha Conjugate Neck Structures in the Collisions of 35 MeV/nucleon 40Ca with 40Ca
Authors:
K. Schmidt,
X. Cao,
E. J. Kim,
K. Hagel,
M. Barbui,
J. Gauthier,
S. Wuenschel,
G. Giuliani,
M. R. D. Rodrigues,
H. Zheng,
M. Huang,
N. Blando,
A. Bonasera,
R. Wada,
C. Botosso,
G. Liu,
G. Viesti,
S. Moretto,
G. Prete,
S. Pesente,
D. Fabris,
Y. El Masri,
T. Keutgen,
S. Kolwaski,
A. Kumar
, et al. (2 additional authors not shown)
Abstract:
The de-excitation of alpha-conjugate nuclei produced in reactions of 35 MeV/nucleon 40Ca with 40Ca has been investigated. Particular emphasis is placed on examining the dynamics of collisions leading to projectile-like fragment exit channels. A general exploration of the reaction systematics reveals the binary dissipative character of the collisions and a hierarchy effect similar to that seen for…
▽ More
The de-excitation of alpha-conjugate nuclei produced in reactions of 35 MeV/nucleon 40Ca with 40Ca has been investigated. Particular emphasis is placed on examining the dynamics of collisions leading to projectile-like fragment exit channels. A general exploration of the reaction systematics reveals the binary dissipative character of the collisions and a hierarchy effect similar to that seen for heavier systems. Investigation of the subset of events characterized by a total alpha-conjugate mass (alpha particles plus alpha-conjugate fragments) equal to 40 and atomic number equal to 20 reveals a dominance of alpha-conjugate exit channels. The hierarchy effect for these channels leads to the production of alpha-clustered neck structures with potentially exotic geometries and properties.
△ Less
Submitted 22 May, 2017; v1 submitted 1 December, 2016;
originally announced December 2016.
-
The many facets of the (non relativistic) Nuclear Equation of State
Authors:
G. Giuliani,
H. Zheng,
A. Bonasera
Abstract:
A nucleus is a quantum many body system made of strongly interacting Fermions, protons and neutrons (nucleons). This produces a rich Nuclear Equation of State whose knowledge is crucial to our understanding of the composition and evolution of celestial objects. The nuclear equation of state displays many different features; first neutrons and protons might be treated as identical particles or nucl…
▽ More
A nucleus is a quantum many body system made of strongly interacting Fermions, protons and neutrons (nucleons). This produces a rich Nuclear Equation of State whose knowledge is crucial to our understanding of the composition and evolution of celestial objects. The nuclear equation of state displays many different features; first neutrons and protons might be treated as identical particles or nucleons, but when the differences between protons and neutrons are spelled out, we can have completely different scenarios, just by changing slightly their interactions. At zero temperature and for neutron rich matter, a quantum liquid gas phase transition at low densities or a quark-gluon plasma at high densities might occur. Furthermore, the large binding energy of the $α$ particle, a Boson, might also open the possibility of studying a system made of a mixture of Bosons and Fermions, which adds to the open problems of the nuclear equation of state.
△ Less
Submitted 15 January, 2014; v1 submitted 6 November, 2013;
originally announced November 2013.
-
Measurement of the plasma astrophysical S factor for the 3He(D, p)4He reaction in exploding molecular clusters
Authors:
M. Barbui,
W. Bang,
A. Bonasera,
K. Hagel,
K. Schmidt,
J. B. Natowitz,
R. Burch,
G. Giuliani,
M. Barbarino,
H. Zheng,
G. Dyer,
H. J. Quevedo,
E. Gaul,
A. C. Bernstein,
M. Donovan,
S. Kimura,
M. Mazzocco,
F. Consoli,
R. De Angelis,
P. Andreoli,
T. Ditmire
Abstract:
The plasma astrophysical S factor for the 3He(D, p)4He fusion reaction was measured for the first time at temperatures of few keV, using the interaction of intense ultrafast laser pulses with molecular deuterium clusters mixed with 3He atoms. Different proportions of D2 and 3He or CD4 and 3He were mixed in the gas jet target in order to allow the measurement of the cross-section for the 3He(D, p)4…
▽ More
The plasma astrophysical S factor for the 3He(D, p)4He fusion reaction was measured for the first time at temperatures of few keV, using the interaction of intense ultrafast laser pulses with molecular deuterium clusters mixed with 3He atoms. Different proportions of D2 and 3He or CD4 and 3He were mixed in the gas jet target in order to allow the measurement of the cross-section for the 3He(D, p)4He reaction. The yield of 14.7 MeV protons from the 3He(D, p)4He reaction was measured in order to extract the astrophysical S factor at low energies. Our result is in agreement with other S factor parameterizations found in the literature.
△ Less
Submitted 1 July, 2013;
originally announced July 2013.
-
Coulomb corrections to density and temperature of bosons in heavy ion collisions
Authors:
Hua Zheng,
Gianluca Giuliani,
Aldo Bonasera
Abstract:
A recently proposed method, based on quadrupole and multiplicity fluctuations in heavy ion collisions, is modified in order to take into account distortions due to the Coulomb field. This is particularly interesting for bosons produced in heavy ion collisions, such as $d$ and $α$ particles. We derive temperatures and densities seen by the bosons and compare to similar calculations for fermions. Th…
▽ More
A recently proposed method, based on quadrupole and multiplicity fluctuations in heavy ion collisions, is modified in order to take into account distortions due to the Coulomb field. This is particularly interesting for bosons produced in heavy ion collisions, such as $d$ and $α$ particles. We derive temperatures and densities seen by the bosons and compare to similar calculations for fermions. The resulting energy densities agree rather well with each other and with the one derived from neutrons. This suggests that a common phenomenon, such as the sudden opening of many reaction channels and/or a liquid gas phase transition, is responsible for the agreement.
△ Less
Submitted 31 July, 2013; v1 submitted 24 June, 2013;
originally announced June 2013.
-
Coulomb corrections to density and temperature in heavy ion collisions
Authors:
Hua Zheng,
Gianluca Giuliani,
Aldo Bonasera
Abstract:
A recently proposed method, based on quadrupole and multiplicity fluctuations in heavy ion collisions, is modified in order to take into account distortions due to the Coulomb field. The classical and quantum limits for fermions are discussed. In the classical case we find that the temperature determined from ${}^3He$ and ${}^3H$, after the Coulomb correction, are very similar to those obtained fr…
▽ More
A recently proposed method, based on quadrupole and multiplicity fluctuations in heavy ion collisions, is modified in order to take into account distortions due to the Coulomb field. The classical and quantum limits for fermions are discussed. In the classical case we find that the temperature determined from ${}^3He$ and ${}^3H$, after the Coulomb correction, are very similar to those obtained from neutrons within the Constrained Molecular Dynamics (CoMD) approach. In the quantum case, the proton temperature becomes very similar to neutron's, while densities are not sensitive to the Coulomb corrections.
△ Less
Submitted 21 May, 2013;
originally announced May 2013.
-
Radon backgrounds in the DEAP-1 liquid-argon-based Dark Matter detector
Authors:
P. -A. Amaudruz,
M. Batygov,
B. Beltran,
K. Boudjemline,
M. G. Boulay B. Cai T. Caldwell,
M. Chen,
R. Chouinard,
B. T. Cleveland,
D. Contreras,
K. Dering,
F. Duncan,
R. Ford,
R. Gagnon F. Giuliani,
M. Gold V. V. Golovko,
P. Gorel,
K. Graham,
D. R. Grant,
R. Hakobyan,
A. L. Hallin,
P. Harvey,
C. Hearns,
C. J. Jillings,
M. Kuźniak,
I. Lawson,
O. Li
, et al. (27 additional authors not shown)
Abstract:
The DEAP-1 \SI{7}{kg} single phase liquid argon scintillation detector was operated underground at SNOLAB in order to test the techniques and measure the backgrounds inherent to single phase detection, in support of the \mbox{DEAP-3600} Dark Matter detector. Backgrounds in DEAP are controlled through material selection, construction techniques, pulse shape discrimination and event reconstruction.…
▽ More
The DEAP-1 \SI{7}{kg} single phase liquid argon scintillation detector was operated underground at SNOLAB in order to test the techniques and measure the backgrounds inherent to single phase detection, in support of the \mbox{DEAP-3600} Dark Matter detector. Backgrounds in DEAP are controlled through material selection, construction techniques, pulse shape discrimination and event reconstruction. This report details the analysis of background events observed in three iterations of the DEAP-1 detector, and the measures taken to reduce them.
The $^{222}$Rn decay rate in the liquid argon was measured to be between 16 and \SI{26}{\micro\becquerel\per\kilogram}. We found that the background spectrum near the region of interest for Dark Matter detection in the DEAP-1 detector can be described considering events from three sources: radon daughters decaying on the surface of the active volume, the expected rate of electromagnetic events misidentified as nuclear recoils due to inefficiencies in the pulse shape discrimination, and leakage of events from outside the fiducial volume due to imperfect position reconstruction. These backgrounds statistically account for all observed events, and they will be strongly reduced in the DEAP-3600 detector due to its higher light yield and simpler geometry.
△ Less
Submitted 23 April, 2014; v1 submitted 5 November, 2012;
originally announced November 2012.
-
Isospin Dependence of Incomplete Fusion Reactions at 25 Mev/a
Authors:
F. Amorini,
G. Cardella,
G. Giuliani,
M. Papa,
C. Agodi,
R. Alba,
A. Anzalone,
I. Berceanu,
S. Cavallaro,
M. B. Chatterjee,
R. Coniglione,
E. De Filippo,
A. Di Pietro,
E. Geraci,
L. Grassi,
A. Grzeszczuk,
P. Figuera,
E. La Guidara,
G. Lanzalone,
N. Le Neindre,
I. Lombardo,
C. Maiolino,
A. Pagano,
S. Pirrone,
G. Politi
, et al. (7 additional authors not shown)
Abstract:
40Ca+40,48Ca,46Ti reactions at 25 MeV/A have been studied using the 4p CHIMERA detector. An isospin effect on the competition between incomplete fusion and dissipative binary reaction mechanisms has been observed. The probability of producing a compound system is observed to be lower in the case of N=Z colliding systems as compared to the case of reactions induced on the more neutron rich 48Ca t…
▽ More
40Ca+40,48Ca,46Ti reactions at 25 MeV/A have been studied using the 4p CHIMERA detector. An isospin effect on the competition between incomplete fusion and dissipative binary reaction mechanisms has been observed. The probability of producing a compound system is observed to be lower in the case of N=Z colliding systems as compared to the case of reactions induced on the more neutron rich 48Ca target. Predictions based on CoMD-II calculations show that the competition between fusion-like and dissipative reactions, for the selected centrality, can strongly constraint the parameterization of symmetry energy and its density dependence in the nuclear equation of state.
△ Less
Submitted 8 September, 2008;
originally announced September 2008.