Showing 1–50 of 68 results for author: Crawford, L

The highest mass Kepler red giants -- I. Global asteroseismic parameters of 48 stars

Authors: Courtney L. Crawford, Timothy R. Bedding, Yaguang Li, Dennis Stello, Daniel Huber, Jie Yu, K. R. Sreenivas, Tanda Li, Emily F. Kerrison

Abstract: When low- and intermediate-mass stars evolve off the main sequence, they expand and cool into the red giant stages of evolution, which include those associated with shell H burning (the red giant branch), core He burning (the red clump), and shell He burning (the asymptotic giant branch). The majority of red giants have masses $< 2 M_\odot$, and red giants more massive than this are often excluded… ▽ More When low- and intermediate-mass stars evolve off the main sequence, they expand and cool into the red giant stages of evolution, which include those associated with shell H burning (the red giant branch), core He burning (the red clump), and shell He burning (the asymptotic giant branch). The majority of red giants have masses $< 2 M_\odot$, and red giants more massive than this are often excluded from major studies. Here we present a study of the highest-mass stars ($M > 3.0 M_\odot$) in the Kepler sample of 16,000 red giants. We begin by re-estimating their global seismic properties with new light curves, highlighting the differences between using the SAP and PDCSAP light curves provided by Kepler. We use the re-estimated properties to derive new mass estimates for the stars, ending with a final sample of 48 confirmed high-mass stars. We explore their oscillation envelopes, confirming the trends found in recent works such as low mean mode amplitude and wide envelopes. We find, through probabilistic means, that our sample is likely all core He burning stars. We measure their dipole and quadrupole mode visibilities and confirm that the dipole mode visibility tends to decrease with mass. △ Less

Submitted 11 February, 2024; originally announced February 2024.

Comments: 14 pages, 12 figures, 1 table, accepted to MNRAS

Abrasion-fission reactions at intermediate energies

Authors: M. Bowry, O. B. Tarasov, J. S. Berryman, V. Bader, D. Bazin, T. Chupp, H. L. Crawford, A. Gade, E. Lunderberg, A. Ratkiewicz, F. Recchia, B. M. Sherrill, D. Smalley, A. Stolz, S. R. Stroberg, D. Weisshaar, S. Williams, K. Wimmer, J. Yurkon

Abstract: The availability of high-intensity, heavy-ion beams coupled to sensitive, large solid-angleacceptance spectrometers has enabled a detailed examination of the fission fragments produced in induced-fission reactions. The abrasion-fission process involves the formation of projectile-like prefragments in violent nuclear collisions at relative energies in excess of 100 MeV/u. At intermediate energies b… ▽ More The availability of high-intensity, heavy-ion beams coupled to sensitive, large solid-angleacceptance spectrometers has enabled a detailed examination of the fission fragments produced in induced-fission reactions. The abrasion-fission process involves the formation of projectile-like prefragments in violent nuclear collisions at relative energies in excess of 100 MeV/u. At intermediate energies below this threshold, experiments suggest a change in the prefragment kinematic qualities. Information regarding the influence of this transitional phase upon the evolution of nuclei approaching the point of scission is scarce. In this article, data are presented for over 200 nuclei from nickel to palladium produced in abrasion-fission reactions of a 80 MeV/u 238U beam. Cross sections were obtained following yield measurements performed for the principal charge states of the identified fission fragments and a detailed analysis of the ion transmission. A full kinematic analysis of the fission fragments has been performed using the LISE++ software package, where the trajectory of an ion passing through a spectrometer can be reconstructed based upon measurements at the focal plane. The results obtained at the S800 spectrograph are compared with predictions obtained with a three-fission progenitor (3EER) model. Systematic studies of fission-fragment properties continue to provide a valuable experimental benchmark for theoretical efforts directed toward describing this complex decay channel, that is important in the context of planning experiments to explore the neutron-rich region of the nuclear chart at rare-isotope beam facilities. △ Less

Submitted 31 January, 2024; originally announced January 2024.

Comments: 14 pages, 18 figures

Journal ref: Phys. Rev. C 108, 034604 (2023)

A simple method to measure numax for asteroseismology: application to 16,000 oscillating Kepler red giants

Authors: K. R. Sreenivas, Timothy R. Bedding, Yaguang Li, Daniel Huber, Courtney L. Crawford, Dennis Stello, Jie Yu

Abstract: The importance of numax (the frequency of maximum oscillation power) for asteroseismology has been demonstrated widely in the previous decade, especially for red giants. With the large amount of photometric data from CoRoT, Kepler and TESS, several automated algorithms to retrieve numax values have been introduced. Most of these algorithms correct the granulation background in the power spectrum b… ▽ More The importance of numax (the frequency of maximum oscillation power) for asteroseismology has been demonstrated widely in the previous decade, especially for red giants. With the large amount of photometric data from CoRoT, Kepler and TESS, several automated algorithms to retrieve numax values have been introduced. Most of these algorithms correct the granulation background in the power spectrum by fitting a model and subtracting it before measuring numax. We have developed a method that does not require fitting to the granulation background. Instead, we simply divide the power spectrum by a function of the form nu^-2, to remove the slope due to granulation background, and then smooth to measure numax. This method is fast, simple and avoids degeneracies associated with fitting. The method is able to measure oscillations in 99.9% of previously-studied Kepler red giants, with a systematic offset of 1.5 % in numax values that that we are able to calibrate. On comparing the seismic radii from this work with Gaia, we see similar trends to those observed in previous studies. Additionally, our values of width of the power envelope can clearly identify the dipole mode suppressed stars as a distinct population, hence as a way to detect them. We also applied our method to stars with low (0.19--18.35 muHz) and found it works well to correctly identify the oscillations. △ Less

Submitted 15 April, 2024; v1 submitted 30 January, 2024; originally announced January 2024.

Comments: Published in MNRAS

A Vision for the Science of Rare Isotopes

Authors: H. L. Crawford, K. Fossez, S. König, A. Spyrou

Abstract: The field of nuclear science has considerably advanced since its beginning just over a century ago. Today, the science of rare isotopes is on the cusp of a new era with theoretical and computing advances complementing experimental capabilities at new facilities internationally. In this article we present a vision for the science of rare isotope beams (RIBs). We do not attempt to cover the full bre… ▽ More The field of nuclear science has considerably advanced since its beginning just over a century ago. Today, the science of rare isotopes is on the cusp of a new era with theoretical and computing advances complementing experimental capabilities at new facilities internationally. In this article we present a vision for the science of rare isotope beams (RIBs). We do not attempt to cover the full breadth of the field, but rather provide a perspective and address a selection of topics that reflect our own interests and expertise. We focus in particular on systems near the drip lines, where one often finds nuclei that are referred to as "exotic," and where the role of the "nuclear continuum" is only just starting to be explored. An important aspect of this article is the attempt to highlight the crucial connections between nuclear structure and nuclear reactions required to fully interpret and leverage the rich data to be collected in the next years at RIB facilities. Further, we connect the efforts in structure and reactions to key questions of nuclear astrophysics. △ Less

Submitted 14 December, 2023; originally announced December 2023.

HdC and EHe stars through the prism of Gaia DR3: 3D distribution and Gaia's chromatic PSF effects

Authors: P. Tisserand, C. L. Crawford, J. Soon, G. C. Clayton, A. J. Ruiter, I. R. Seitenzahl

Abstract: Upon its release the Gaia DR3 catalogue has led to tremendous progress in multiple fields of astronomy by providing the complete astrometric solution for nearly 1.5 billion sources. We analysed the photometric and astrometric results for Hydrogen-deficient Carbon (HdC), Extreme Helium (EHe), and DYPer type stars to identify any potential biases and thus to select stars suitable for kinematic and s… ▽ More Upon its release the Gaia DR3 catalogue has led to tremendous progress in multiple fields of astronomy by providing the complete astrometric solution for nearly 1.5 billion sources. We analysed the photometric and astrometric results for Hydrogen-deficient Carbon (HdC), Extreme Helium (EHe), and DYPer type stars to identify any potential biases and thus to select stars suitable for kinematic and spatial distribution studies. We investigated the impact of photometric declines in R Coronae Borealis (RCB) stars on Gaia astrometry using information from the Gaia IPD process cross-matched with light curves. We have reached the conclusion that the astrometric fits for numerous RCB stars are not valid due to the Gaia PSF chromaticity effect in both shape and centroid. The astrometric results of all stars with a significant time-dependent colour variation should be similarly affected. RCB stars might thus be promising sources to correct this effect in future Gaia releases. Furthermore, after validating the Gaia astrometric results for 92 stars, we observed that the majority of HdC and EHe stars are distributed across the three old stellar structures, the thick disk, the bulge and the halo. However, we have also uncovered evidence indicating that some of them exhibit orbits characteristic of the thin disk. This is also particularly true for all DYPer type stars studied. Finally, we have produced a list of star memberships for each Galactic substructure, and provided a list of heliocentric radial velocities and associated errors for targets not observed by Gaia DR3. We are beginning to observe a relationship between kinematics, stellar population, and metallicity in RCB and EHe stars. That relation can be explained, within the double degenerate scenario, by the large range in the delay time distribution expected from population synthesis simulations, particularly through the HybCO merger channel. △ Less

Submitted 16 January, 2024; v1 submitted 18 September, 2023; originally announced September 2023.

Comments: 22 pages, 20 figures, accepted for publication in A&A v2: RV section from accompanying paper ArXiv2309.10139 moved into that one +New section added on the comparison between our classification of HdC and EHe stars within a Galactic substructure and the one published lately by Philip Monai et al. (2023), and between both results on dynamical orbital parameters as determined with Galpy simulations

HdC and EHe stars through the prism of Gaia DR3: Evolution of RV amplitude and dust formation rate with effective temperature

Authors: P. Tisserand, C. L. Crawford, J. Soon, G. C. Clayton, A. J. Ruiter, I. R. Seitenzahl

Abstract: The Gaia DR3 release includes heliocentric radial velocity measurements and velocity variability indices for tens of millions of stars observed over 34 months.In this study, we utilise these indices to investigate the intrinsic radial velocity variations of Hydrogen-deficient Carbon (HdC) stars and Extreme Helium (EHe) stars across their large ranges of temperature and brightness. Taking advantage… ▽ More The Gaia DR3 release includes heliocentric radial velocity measurements and velocity variability indices for tens of millions of stars observed over 34 months.In this study, we utilise these indices to investigate the intrinsic radial velocity variations of Hydrogen-deficient Carbon (HdC) stars and Extreme Helium (EHe) stars across their large ranges of temperature and brightness. Taking advantage of the newly defined HdC temperature classes, we examine the evolution of the total velocity amplitude with effective temperature. Additionally, we analyse the variation in the dust production rate of R Coronae Borealis (RCB) stars with temperature using two different proxies for the photometric state of RCB stars: one from Gaia and another from the 2MASS survey. Our observations revealed a trend in the evolution of the maximum radial velocity amplitude across each HdC temperature class. Similarly, we also observed a correlation between stellar temperature and the dust production rate. Interestingly, we possibly observed for the first time some variations of the intrinsic radial velocity amplitude and the dust production rate with HdC temperature class. If confirmed, these variations would indicate that the helium shell-burning giant stage starts with strong atmospheric motions that decrease in strength, up to $\sim$6000 K, before picking up again as the HdC star atmosphere shrinks further in size and reaches warmer temperatures. Moreover, the dust formation rate appears to be much higher in colder RCB stars compared to warmer ones. △ Less

Submitted 16 January, 2024; v1 submitted 18 September, 2023; originally announced September 2023.

Comments: 7 pages, 4 figures, accepted in A&A V2: version accepted, radial velocity section of V1 moved into the accompanying paper ArXiv2309.10148

Statistical Inference on Grayscale Images via the Euler-Radon Transform

Authors: Kun Meng, Mattie Ji, Jinyu Wang, Kexin Ding, Henry Kirveslahti, Ani Eloyan, Lorin Crawford

Abstract: Tools from topological data analysis have been widely used to represent binary images in many scientific applications. Methods that aim to represent grayscale images (i.e., where pixel intensities instead take on continuous values) have been relatively underdeveloped. In this paper, we introduce the Euler-Radon transform, which generalizes the Euler characteristic transform to grayscale images by… ▽ More Tools from topological data analysis have been widely used to represent binary images in many scientific applications. Methods that aim to represent grayscale images (i.e., where pixel intensities instead take on continuous values) have been relatively underdeveloped. In this paper, we introduce the Euler-Radon transform, which generalizes the Euler characteristic transform to grayscale images by using o-minimal structures and Euler integration over definable functions. Coupling the Karhunen-Loeve expansion with our proposed topological representation, we offer hypothesis-testing algorithms based on the chi-squared distribution for detecting significant differences between two groups of grayscale images. We illustrate our framework via extensive numerical experiments and simulations. △ Less

Submitted 27 August, 2023; originally announced August 2023.

Comments: 85 pages, 9 figures

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.

Journal ref: Phys.Lett.B 843 (2023) 138038

Should I Stop or Should I Go: Early Stopping with Heterogeneous Populations

Authors: Hammaad Adam, Fan Yin, Huibin, Hu, Neil Tenenholtz, Lorin Crawford, Lester Mackey, Allison Koenecke

Abstract: Randomized experiments often need to be stopped prematurely due to the treatment having an unintended harmful effect. Existing methods that determine when to stop an experiment early are typically applied to the data in aggregate and do not account for treatment effect heterogeneity. In this paper, we study the early stopping of experiments for harm on heterogeneous populations. We first establish… ▽ More Randomized experiments often need to be stopped prematurely due to the treatment having an unintended harmful effect. Existing methods that determine when to stop an experiment early are typically applied to the data in aggregate and do not account for treatment effect heterogeneity. In this paper, we study the early stopping of experiments for harm on heterogeneous populations. We first establish that current methods often fail to stop experiments when the treatment harms a minority group of participants. We then use causal machine learning to develop CLASH, the first broadly-applicable method for heterogeneous early stopping. We demonstrate CLASH's performance on simulated and real data and show that it yields effective early stopping for both clinical trials and A/B tests. △ Less

Submitted 27 October, 2023; v1 submitted 20 June, 2023; originally announced June 2023.

Comments: NeurIPS 2023 (spotlight)

Microsecond Isomer at the N=20 Island of Shape Inversion Observed at FRIB

Authors: T. J. Gray, J. M. Allmond, Z. Xu, T. T. King, R. S. Lubna, H. L. Crawford, V. Tripathi, B. P. Crider, R. Grzywacz, S. N. Liddick, A. O. Macchiavelli, T. Miyagi, A. Poves, A. Andalib, E. Argo, C. Benetti, S. Bhattacharya, C. M. Campbell, M. P. Carpenter, J. Chan, A. Chester, J. Christie, B. R. Clark, I. Cox, A. A. Doetsch , et al. (41 additional authors not shown)

Abstract: Excited-state spectroscopy from the first Facility for Rare Isotope Beams (FRIB) experiment is reported. A 24(2)-$μみゅー$s isomer was observed with the FRIB Decay Station initiator (FDSi) through a cascade of 224- and 401-keV $γがんま$ rays in coincidence with $^{32}\textrm{Na}$ nuclei. This is the only known microsecond isomer ($1{\text{ }μみゅー\text{s}}\leq T_{1/2} < 1\text{ ms}$) in the region. This nucleus is… ▽ More Excited-state spectroscopy from the first Facility for Rare Isotope Beams (FRIB) experiment is reported. A 24(2)-$μみゅー$s isomer was observed with the FRIB Decay Station initiator (FDSi) through a cascade of 224- and 401-keV $γがんま$ rays in coincidence with $^{32}\textrm{Na}$ nuclei. This is the only known microsecond isomer ($1{\text{ }μみゅー\text{s}}\leq T_{1/2} < 1\text{ ms}$) in the region. This nucleus is at the heart of the $N=20$ island of shape inversion and is at the crossroads of spherical shell-model, deformed shell-model, and ab initio theories. It can be represented as the coupling of a proton hole and neutron particle to $^{32}\textrm{Mg}$, $^{32}\textrm{Mg}+πぱい^{-1} + νにゅー^{+1}$. This odd-odd coupling and isomer formation provides a sensitive measure of the underlying shape degrees of freedom of $^{32}\textrm{Mg}$, where the onset of spherical-to-deformed shape inversion begins with a low-lying deformed $2^+$ state at 885 keV and a low-lying shape-coexisting $0_2^+$ state at 1058 keV. We suggest two possible explanations for the 625-keV isomer in $^{32}$Na: a $6^-$ spherical shape isomer that decays by $E2$ or a $0^+$ deformed spin isomer that decays by $M2$. The present results and calculations are most consistent with the latter, indicating that the low-lying states are dominated by deformation. △ Less

Submitted 26 April, 2023; v1 submitted 22 February, 2023; originally announced February 2023.

Comments: 7 pages, 5 figures, accepted by Physical Review Letters

A Simple Approach for Local and Global Variable Importance in Nonlinear Regression Models

Authors: Emily T. Winn-Nuñez, Maryclare Griffin, Lorin Crawford

Abstract: The ability to interpret machine learning models has become increasingly important as their usage in data science continues to rise. Most current interpretability methods are optimized to work on either (\textit{i}) a global scale, where the goal is to rank features based on their contributions to overall variation in an observed population, or (\textit{ii}) the local level, which aims to detail o… ▽ More The ability to interpret machine learning models has become increasingly important as their usage in data science continues to rise. Most current interpretability methods are optimized to work on either (\textit{i}) a global scale, where the goal is to rank features based on their contributions to overall variation in an observed population, or (\textit{ii}) the local level, which aims to detail on how important a feature is to a particular individual in the data set. In this work, a new operator is proposed called the "GlObal And Local Score" (GOALS): a simple \textit{post hoc} approach to simultaneously assess local and global feature variable importance in nonlinear models. Motivated by problems in biomedicine, the approach is demonstrated using Gaussian process regression where the task of understanding how genetic markers are associated with disease progression both within individuals and across populations is of high interest. Detailed simulations and real data analyses illustrate the flexible and efficient utility of GOALS over state-of-the-art variable importance strategies. △ Less

Submitted 10 August, 2023; v1 submitted 3 February, 2023; originally announced February 2023.

A Spectral Classification System for Hydrogen-deficient Carbon Stars

Authors: Courtney L. Crawford, Patrick Tisserand, Geoffrey C. Clayton, Jamie Soon, Mike Bessell, Peter Wood, D. A. Garcia-Hernandez, Ashley J. Ruiter, Ivo R. Seitenzahl

Abstract: Stellar spectral classification has been highly useful in the study of stars. While there is a currently accepted spectral classification system for carbon stars, the subset of Hydrogen-deficient Carbon (HdC) stars has not been well described by such a system, due predominantly to their rarity and their variability. Here we present the first system for the classification of HdCs based on their spe… ▽ More Stellar spectral classification has been highly useful in the study of stars. While there is a currently accepted spectral classification system for carbon stars, the subset of Hydrogen-deficient Carbon (HdC) stars has not been well described by such a system, due predominantly to their rarity and their variability. Here we present the first system for the classification of HdCs based on their spectra, which is made wholly on their observable appearance. We use a combination of dimensionality reduction and clustering algorithms with human classification to create such a system with eight total classes corresponding to temperature, and an additional second axis corresponding to the carbon molecular band strength. We classify over half of the known sample of HdC stars using this, and roughly calibrate the temperatures of each class using their colors. Additionally, we express trends in the occurrence of certain spectral peculiarities such as the presence of Hydrogen and Lithium lines. We also present three previously unpublished spectra, report the discovery of two new Galactic dustless HdC (dLHdC) stars and additionally discuss one especially unique star that appears to border between the hottest HdCs and the coolest Extreme Helium (EHe) stars. △ Less

Submitted 10 January, 2023; v1 submitted 9 October, 2022; originally announced October 2022.

Comments: 30 pages, 18 figures, submitted to MNRAS January 11, 2023 revision: Changes updated according to referee comments

E2 Rotational Invariants of $0^+_1$ and $2^+_1$ states for $^{106}$Cd: the Emergence of Collective Rotation

Authors: T. J. Gray, J. M. Allmond, R. V. F. Janssens, W. Korten, A. E. Stuchbery, J. L. Wood, A. D. Ayangeakaa, S. Bottoni, B. M. Bucher, C. M. Campbell, M. P. Carpenter, H. L. Crawford, H. David, D. Doherty, P. Fallon, M. T. Febbraro, A. Galindo-Uribarri, C. J. Gross, M. Komorowska, F. G. Kondev, T. Lauritsen, A. O. Macchiavelli, P. Napiorkowsi, E. Padilla-Rodal, S. D. Pain , et al. (7 additional authors not shown)

Abstract: The collective structure of $^{106}$Cd is elucidated by multi-step Coulomb excitation of a 3.849 MeV/$A$ beam of $^{106}$Cd on a 1.1 mg/cm$^2$ $^{208}$Pb target using GRETINA-CHICO2 at ATLAS. Fourteen $E2$ matrix elements were obtained. The nucleus $^{106}$Cd is a prime example of emergent collectivity that possesses a simple structure: it is free of complexity caused by shape coexistence and has… ▽ More The collective structure of $^{106}$Cd is elucidated by multi-step Coulomb excitation of a 3.849 MeV/$A$ beam of $^{106}$Cd on a 1.1 mg/cm$^2$ $^{208}$Pb target using GRETINA-CHICO2 at ATLAS. Fourteen $E2$ matrix elements were obtained. The nucleus $^{106}$Cd is a prime example of emergent collectivity that possesses a simple structure: it is free of complexity caused by shape coexistence and has a small, but collectively active number of valence nucleons. This work follows in a long and currently active quest to answer the fundamental question of the origin of nuclear collectivity and deformation, notably in the cadmium isotopes. The results are discussed in terms of phenomenological models, the shell model, and Kumar-Cline sums of $E2$ matrix elements. The ${\langle 0_2^+ ||E2||2_1^+ \rangle}$ matrix element is determined for the first time, providing a total, converged measure of the electric quadrupole strength, $\langle Q^2 \rangle$, of the first-excited $2_1^+$ level relative to the $0_1^+$ ground state, which does not show an increase as expected of harmonic and anharmonic vibrations. Strong evidence for triaxial shapes in weakly collective nuclei is indicated; collective vibrations are excluded. This is contrary to the only other cadmium result of this kind in $^{114}$Cd by C. Fahlander et al., Nucl. Phys. A485, 327 (1988), which is complicated by low-lying shape coexistence near midshell. △ Less

Submitted 13 September, 2022; originally announced September 2022.

Comments: 7 pages, 4 figures

The Core of $^{25}$F studied by the $^{25}$F(-1p)$^{24}$O reaction

Authors: H. L. Crawford, M. D. Jones, A. O. Macchiavelli, P. Fallon, D. Bazin, P. C. Bender, B. A. Brown, C. M. Campbell, R. M. Clark, M. Cromaz, B. Elman, A. Gade, J. D. Holt, R. V. F. Janssens, I. Y. Lee, B. Longfellow, S. Paschalis, M. Petri, A. L. Richard, M. Salathe, J. A. Tostevin, D. Weisshaar

Abstract: The $^{25}$F($5/2^+) (-1p) ^{24}$O reaction was studied at the NSCL using the S800 spectrometer. The experimental spectroscopic factor for the ground-state to ground-state transition indicates a substantial depletion of the proton $d_{5/2}$ strength compared to shell-model expectations. Our result supports the findings reported by Tang \textit{et al.}, from their study of the $(p,2p)$ reaction at… ▽ More The $^{25}$F($5/2^+) (-1p) ^{24}$O reaction was studied at the NSCL using the S800 spectrometer. The experimental spectroscopic factor for the ground-state to ground-state transition indicates a substantial depletion of the proton $d_{5/2}$ strength compared to shell-model expectations. Our result supports the findings reported by Tang \textit{et al.}, from their study of the $(p,2p)$ reaction at RIBF. The overlap between the $^{25}$F and $^{24}$O ground-states is considerably less than anticipated if $^{24}$O acted as a robust and rigid doubly-magic core in $^{25}$F. We interpret the results within the framework of the Particle-Vibration Coupling (PVC) of a $d_{5/2}$ proton coupled to a quadrupole phonon of an effective core. This approach provides a good description of the experimental data by requiring an effective $^{24}$O* core with a phonon energy of $\hbarωおめが_2$= 3.2 MeV, and a $B(E2) ~ 2.7$ W.u., softer and more collective than a bare $^{24}$O. Both the Nilsson deformed mean field and the PVC models appear to capture the properties of the effective core of $^{25}$F, suggesting that the additional proton tends to polarize the free, doubly magic $^{24}$O in such a way that it becomes either slightly deformed or a quadrupole vibrator. △ Less

Submitted 7 September, 2022; v1 submitted 5 September, 2022; originally announced September 2022.

Comments: 6 pages, 2 figures

Quenching of Single-Particle Strength in A=15 Nuclei

Authors: B. P. Kay, T. L. Tang, I. A. Tolstukhin, G. B. Roderick, A. J. Mitchell, Y. Ayyad, S. A. Bennett, J. Chen, K. A. Chipps, H. L. Crawford, S. J. Freeman, K. Garrett, M. D. Gott, M. R. Hall, C. R. Hoffman, H. Jayatissa, A. O. Macchiavelli, P. T. MacGregor, D. K. Sharp, G. L. Wilson

Abstract: Absolute cross sections for the addition of $s$- and $d$-wave neutrons to $^{14}$C and $^{14}$N have been determined simultaneously via the ($d$,$p$) reaction at 10 MeV/u. The difference between the neutron and proton separation energies, $ΔでるたS$, is around $-20$ MeV for the $^{14}$C$+$$n$ system and $+8$ MeV for $^{14}$N$+$$n$. The population of the $1s_{1/2}$ and $0d_{5/2}$ orbitals for both system… ▽ More Absolute cross sections for the addition of $s$- and $d$-wave neutrons to $^{14}$C and $^{14}$N have been determined simultaneously via the ($d$,$p$) reaction at 10 MeV/u. The difference between the neutron and proton separation energies, $ΔでるたS$, is around $-20$ MeV for the $^{14}$C$+$$n$ system and $+8$ MeV for $^{14}$N$+$$n$. The population of the $1s_{1/2}$ and $0d_{5/2}$ orbitals for both systems is reduced by a factor of approximately 0.5 compared to the independent single-particle model, or about 0.6 when compared to the shell model. This finding strongly contrasts with results deduced from intermediate-energy knockout reactions between similar nuclei on targets of $^{9}$Be and $^{12}$C. The simultaneous technique used removes many systematic uncertainties. △ Less

Submitted 5 July, 2022; originally announced July 2022.

Comments: 6 pages, 4 figures

Randomness of Shapes and Statistical Inference on Shapes via the Smooth Euler Characteristic Transform

Authors: Kun Meng, Jinyu Wang, Lorin Crawford, Ani Eloyan

Abstract: In this article, we establish the mathematical foundations for modeling the randomness of shapes and conducting statistical inference on shapes using the smooth Euler characteristic transform. Based on these foundations, we propose two chi-squared statistic-based algorithms for testing hypotheses on random shapes. Simulation studies are presented to validate our mathematical derivations and to com… ▽ More In this article, we establish the mathematical foundations for modeling the randomness of shapes and conducting statistical inference on shapes using the smooth Euler characteristic transform. Based on these foundations, we propose two chi-squared statistic-based algorithms for testing hypotheses on random shapes. Simulation studies are presented to validate our mathematical derivations and to compare our algorithms with state-of-the-art methods to demonstrate the utility of our proposed framework. As real applications, we analyze a data set of mandibular molars from four genera of primates and show that our algorithms have the power to detect significant shape differences that recapitulate known morphological variation across suborders. Altogether, our discussions bridge the following fields: algebraic and computational topology, probability theory and stochastic processes, Sobolev spaces and functional analysis, analysis of variance for functional data, and geometric morphometrics. △ Less

Submitted 23 May, 2024; v1 submitted 27 April, 2022; originally announced April 2022.

Comments: 110 pages

Journal ref: Journal of the American Statistical Association, 2024

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.

Comments: 12 pages, 14 figures, accepted for publication in Physical Review C

Weak binding effects on the structure of $^{40}$Mg

Authors: A. O. Macchiavelli, H. L. Crawford, P. Fallon, R. M. Clark, A. Poves

Abstract: While the phenomenon of one- and two-neutron ground-state halo nuclei is well established, the effects of weak binding on nuclear excitation properties remain largely unexplored. Motivated by this question and by recent data in $^{40}$Mg we investigate the coupling of weakly bound (halo) valence neutrons to a core using the known properties of $^{40}$Mg to explore and illustrate possible particle-… ▽ More While the phenomenon of one- and two-neutron ground-state halo nuclei is well established, the effects of weak binding on nuclear excitation properties remain largely unexplored. Motivated by this question and by recent data in $^{40}$Mg we investigate the coupling of weakly bound (halo) valence neutrons to a core using the known properties of $^{40}$Mg to explore and illustrate possible particle-core coupling schemes and their impact on the low-lying excitation spectrum. △ Less

Submitted 10 January, 2022; originally announced January 2022.

Comments: 7 pages, 7 figures

The dawn of a new era for dustless HdC stars with GAIA eDR3

Authors: P. Tisserand, C. L. Crawford, G. C. Clayton, A. J. Ruiter, V. Karambelkar, M. S. Bessell, I. R. Seitenzahl, M. M. Kasliwal, J. Soon, T. Travouillon

Abstract: Decades after their discovery, only four hydrogen-deficient carbon (HdC) stars were known to have no circumstellar dust shell. This is in complete contrast to the $\sim$130 known Galactic HdC stars that are notorious for being heavy dust producers, i.e. the R Coronae Borealis (RCB) stars. Together they form a rare class of supergiant stars that are thought to originate from the merger of CO/He whi… ▽ More Decades after their discovery, only four hydrogen-deficient carbon (HdC) stars were known to have no circumstellar dust shell. This is in complete contrast to the $\sim$130 known Galactic HdC stars that are notorious for being heavy dust producers, i.e. the R Coronae Borealis (RCB) stars. Together they form a rare class of supergiant stars that are thought to originate from the merger of CO/He white dwarf binary systems, otherwise known as the double-degenerate scenario. We searched for new dustless HdC (dLHdC) stars. We used primarily the 2MASS and GAIA eDR3 catalogues to select candidates that were followed-up spectroscopically. We discovered 27 new dLHdC stars, one new RCB star and two new EHe stars. Surprisingly, 20 of the new dLHdC stars share a characteristic of the known dLHdC star HD 148839, having lower atmospheric hydrogen deficiencies. The uncovered population of dLHdC stars exhibit a Bulge-like distribution, like the RCB stars, but show multiple differences from those that indicate they are a different population of HdC stars following its own evolutionary sequence with a fainter luminosity and also a narrow range of effective temperature. We found indication of a current low dust production activity for four of the new dLHdC stars which could be typical RCB stars passing through a transition time. We have evidence for the first time of a large range of absolute magnitudes in the overall population of HdC stars, spanning over 3 mag. In the favoured formation framework, this is explained by a large range in the initial total WD binary mass which leads to a series of evolutionary sequences with distinct maximum brightness and initial temperature. The cold Galactic RCB stars are also noticeably fainter than the Magellanic ones, possibly due to a difference in metallicity resulting in different WD mass ratio. In our Galaxy, there could be as many dLHdC stars as RCB stars. △ Less

Submitted 27 July, 2022; v1 submitted 14 December, 2021; originally announced December 2021.

Comments: 22 pages, 18 figures, accepted in A&A. V3: minor corrections on Table 2 and 3 related to HD stars names and (V-I)0 colours (Fig.10 changed accordingly). V4: changed the colour to light blue for the two new EHe stars indicated in Fig.10, page 13

Journal ref: A&A 667, A83 (2022)

Peculiar Hydrogen-deficient Carbon Stars: Strontium-Rich Stars and the s-Process

Authors: Courtney L. Crawford, Patrick Tisserand, Geoffrey C. Clayton, Bradley Munson

Abstract: Context: R Coronae Borealis (RCB) variables and their non-variable counterparts, the dustless Hydrogen-Deficient Carbon (dLHdC) stars have been known to exhibit enhanced s-processed material on their surfaces, especially Sr, Y, and Ba. No comprehensive work has been done to explore the s-process in these types of stars, however one particular RCB star, U Aqr, has been under scrutiny for its extrao… ▽ More Context: R Coronae Borealis (RCB) variables and their non-variable counterparts, the dustless Hydrogen-Deficient Carbon (dLHdC) stars have been known to exhibit enhanced s-processed material on their surfaces, especially Sr, Y, and Ba. No comprehensive work has been done to explore the s-process in these types of stars, however one particular RCB star, U Aqr, has been under scrutiny for its extraordinary Sr enhancement. Aims: We aim to identify RCB and dLHdC stars that have significantly enhanced Sr abundances, such as U Aqr, and use stellar evolution models to begin to estimate the type of neutron exposure that occurs in a typical HdC star. Methods: We compare the strength of the Sr II 4077 $Å$ spectral line to Ca II H to identify the new subclass of Sr-rich HdCs. We additionally use the structural and abundance information from existing RCB MESA models to calculate the neutron exposure parameter, $τたう$ Results: We identify six stars in the Sr-rich class. Two are RCBs, and four are dLHdCs. We additionally find that the preferred RCB MESA model has a neutron exposure $τたう$ ~ 0.1 mb$^{-1}$, which is lower than the estimated $τたう$ between 0.15 and 0.6 mb$^{-1}$ for the Sr-rich star U Aqr found in the literature. We find trends in the neutron exposure corresponding to He-burning shell temperature, metallicity, and assumed s-processing site. Conclusions: We have found a sub-class of 6 HdCs known as the Sr-rich class, which tend to lie in the halo, outside the typical distribution of RCBs and dLHdCs. We find that dLHdC stars are more likely to be Sr-rich than RCBs, with an occurrence rate of ~13\% for dLHdCs and ~2\% for RCBs. This is one of the first potential spectroscopic differences between RCBs and dLHdCs, along with dLHdCs having stronger surface abundances of $^{18}$O. △ Less

Submitted 14 December, 2021; originally announced December 2021.

Comments: 8 pages, submitted to A&A

Journal ref: A&A 667, A85 (2022)

Structure of $^{43}$P and $^{42}$Si in a two-level shape-coexistence model

Authors: A. O. Macchiavelli, H. L. Crawford, C. M. Campbell, R. M. Clark, M. Cromaz, P. Fallon, I. Y. Lee, A. Gade, A. Poves, E. Rice

Abstract: Exclusive cross sections for the $^{43}$P$(-1p)^{42}$Si reaction to the lowest $0^+$ and $2^+$ states, measured at NSCL with GRETINA and the S800, are interpreted in terms of a two-level mixing (collective) model of oblate and prolate co-existing shapes. Using the formalism developed for deformed nuclei we calculate the spectroscopic amplitudes and exclusive cross-sections in the strong coupling l… ▽ More Exclusive cross sections for the $^{43}$P$(-1p)^{42}$Si reaction to the lowest $0^+$ and $2^+$ states, measured at NSCL with GRETINA and the S800, are interpreted in terms of a two-level mixing (collective) model of oblate and prolate co-existing shapes. Using the formalism developed for deformed nuclei we calculate the spectroscopic amplitudes and exclusive cross-sections in the strong coupling limit, where for $^{43}$P the schematic wavefunction includes the coupling of the Nilsson [211]$\frac{1}{2}$ proton orbit. Good agreement with the experimental data is obtained when the amplitude of the oblate configuration is $\gtrsim$ 80\%, suggesting that both nuclei are predominantly oblate, in line with theoretical expectations. △ Less

Submitted 16 August, 2021; originally announced August 2021.

Comments: 5 pages, 3 figures

UCGretina GEANT4 Simulation of the GRETINA Gamma-Ray Energy Tracking Array

Authors: L. A. Riley, D. Weisshaar, H. L. Crawford, M. L. Agiorgousis, C. M. Campbell, M. Cromaz, P. Fallon, A. Gade, S. D. Gregory, E. B. Haldeman, L. R. Jarvis, E. D. Lawson-John, B. Roberts, B. V. Sadler, C. G. Stine

Abstract: UCGretina, a GEANT4 simulation of the GRETINA gamma-ray tracking array of highly-segmented high-purity germanium detectors is described. We have developed a model of the array, in particular of the Quad Module and the capsules, that gives good agreement between simulated and measured photopeak efficiencies over a broad range of gamma-ray energies and reproduces the shape of the measured Compton co… ▽ More UCGretina, a GEANT4 simulation of the GRETINA gamma-ray tracking array of highly-segmented high-purity germanium detectors is described. We have developed a model of the array, in particular of the Quad Module and the capsules, that gives good agreement between simulated and measured photopeak efficiencies over a broad range of gamma-ray energies and reproduces the shape of the measured Compton continuum. Both of these features are needed in order to accurately extract gamma-ray yields from spectra collected in in-beam gamma-ray spectroscopy measurements with beams traveling at $v/c \gtrsim 0.3$ at the National Superconducting Cyclotron Laboratory and the Facility for Rare Isotope Beams. In the process of developing the model, we determined that millimeter-scale layers of passive germanium surrounding the active volumes of the simulated crystals must be included in order to reproduce measured photopeak efficiencies. We adopted a simple model of effective passive layers and developed heuristic methods of determining passive-layer thicknesses by comparison of simulations and measurements for a single crystal and for the full array. Prospects for future development of the model are discussed. △ Less

Submitted 17 April, 2021; originally announced April 2021.

Comments: 26 pages, 12 figures

Journal ref: Nucl. Instrum. Methods Phys. Res. A 1003, 165305 (2021)

Hollow-tree Super: a directional and scalable approach for feature importance in boosted tree models

Authors: Stephane Doyen, Hugh Taylor, Peter Nicholas, Lewis Crawford, Isabella Young, Michael Sughrue

Abstract: Current limitations in boosted tree modelling prevent the effective scaling to datasets with a large feature number, particularly when investigating the magnitude and directionality of various features on classification. We present a novel methodology, Hollow-tree Super (HOTS), to resolve and visualize feature importance in boosted tree models involving a large number of features. Further, HOTS al… ▽ More Current limitations in boosted tree modelling prevent the effective scaling to datasets with a large feature number, particularly when investigating the magnitude and directionality of various features on classification. We present a novel methodology, Hollow-tree Super (HOTS), to resolve and visualize feature importance in boosted tree models involving a large number of features. Further, HOTS allows for investigation of the directionality and magnitude various features have on classification. Using the Iris dataset, we first compare HOTS to Gini Importance, Partial Dependence Plots, and Permutation Importance, and demonstrate how HOTS resolves the weaknesses present in these methods. We then show how HOTS can be utilized in high dimensional neuroscientific data, by taking 60 Schizophrenic subjects and applying the method to determine which brain regions were most important for classification of schizophrenia as determined by the PANSS. HOTS effectively replicated and supported the findings of Gini importance, Partial Dependence Plots and Permutation importance within the Iris dataset. When applied to the schizophrenic brain dataset, HOTS was able to resolve the top 10 most important features for classification, as well as their directionality for classification and magnitude compared to other features. Cross-validation supported that these same 10 features were consistently used in the decision-making process across multiple trees, and these features were localised primarily to the occipital and parietal cortices, commonly disturbed brain regions in those with Schizophrenia. It is imperative that a methodology is developed that is able to handle the demands of working with large datasets that contain a large number of features. HOTS represents a unique way to investigate both the directionality and magnitude of feature importance when working at scale with boosted-tree modelling. △ Less

Submitted 7 April, 2021; originally announced April 2021.

Comments: 28 pages, 1 table, 7 figures, PDF format - Submitted to PLOSONE pending review

MSC Class: 62-08

R Coronae Borealis Star Evolution: Simulating 3D Merger Events to 1D Stellar Evolution Including Large Scale Nucleosynthesis

Authors: Bradley Munson, Emmanouil Chatzopoulos, Juhan Frank, Geoffrey C. Clayton, Courtney L. Crawford, Pavel A. Denissenkov, Falk Herwig

Abstract: R Coronae Borealis (RCB) stars are rare hydrogen-deficient carbon-rich variable supergiants thought to be the result of dynamically unstable white dwarf mergers. We attempt to model RCBs through all the relevant timescales by simulating a merger event in Octo-tiger, a 3D adaptive mesh refinement (AMR) hydrodynamics code and mapping the post-merger object into MESA, a 1D stellar evolution code. We… ▽ More R Coronae Borealis (RCB) stars are rare hydrogen-deficient carbon-rich variable supergiants thought to be the result of dynamically unstable white dwarf mergers. We attempt to model RCBs through all the relevant timescales by simulating a merger event in Octo-tiger, a 3D adaptive mesh refinement (AMR) hydrodynamics code and mapping the post-merger object into MESA, a 1D stellar evolution code. We then post-process the nucleosynthesis on a much larger nuclear reaction network to study the enhancement of s-process elements. We present models that match observations or previous studies in most surface abundances, isotopic ratios, early evolution and lifetimes. We also observe similar mixing behavior as previous modeling attempts which result in the partial He-burning products visible on the surface in observations. However, we do note that our sub-solar models lack any enhancement in s-process elements, which we attribute to a lack of hydrogen in the envelope. We also find that the Oxygen-16/Oxygen-18 isotopic ratio is very sensitive to initial hydrogen abundance and increases outside of the acceptable range with a hydrogen mass fraction greater than $10^{-4}$. △ Less

Submitted 2 March, 2021; originally announced March 2021.

Comments: 18 pages, 17 figures, 1 table. To be published in The Astrophysical Journal

Coriolis coupling effects in proton-pickup spectroscopic factors from $^{12}$B

Authors: A. O. Macchiavelli, H. L. Crawford, R. M. Clark, P. Fallon, I. Y. Lee, C. Morse, C. M. Campbell, M. Cromaz, C. Santamaria, J. Chen, C. R. Hoffman, B. P. Kay

Abstract: Spectroscopic factors to low-lying negative-parity states in $^{11}$Be extracted from the $^{12}$B($d$,$^3$He)$^{11}$Be proton-removal reaction are interpreted within the rotational model. Earlier predictions of the $p$-wave proton removal strengths in the strong coupling limit of the Nilsson model underestimated the spectroscopic factors to the $3/2^-_1$ and $5/2^-_1$ states and suggested that de… ▽ More Spectroscopic factors to low-lying negative-parity states in $^{11}$Be extracted from the $^{12}$B($d$,$^3$He)$^{11}$Be proton-removal reaction are interpreted within the rotational model. Earlier predictions of the $p$-wave proton removal strengths in the strong coupling limit of the Nilsson model underestimated the spectroscopic factors to the $3/2^-_1$ and $5/2^-_1$ states and suggested that deviations in the $1^+$ ground state of the odd-odd $^{12}$B due to Coriolis coupling should be further explored. In this work we use the Particle Rotor Model to take into account these effects and obtain a good description of the level scheme in $^{11}$B, with a moderate $K$-mixing of the proton Nilsson levels [110]1/2 and [101]3/2. This mixing, present in the $1^+$ bandhead of $^{12}$B, is key to explaining the proton pickup data. △ Less

Submitted 8 December, 2020; originally announced December 2020.

Comments: 5 Pages, 2 Figures, to be submitted to Phys. Rev. C

Journal ref: Phys. Rev. C 103, 034307 (2021)

Joint analysis of structural connectivity and cortical surface features: correlates with mild traumatic brain injury

Authors: Cailey I. Kerley, Leon Y. Cai, Chang Yu, Logan M. Crawford, Jason M. Elenberger, Eden S. Singh, Kurt G. Schilling, Katherine S. Aboud, Bennett A. Landman, Tonia S. Rex

Abstract: Mild traumatic brain injury (mTBI) is a complex syndrome that affects up to 600 per 100,000 individuals, with a particular concentration among military personnel. About half of all mTBI patients experience a diverse array of chronic symptoms which persist long after the acute injury. Hence, there is an urgent need for better understanding of the white matter and gray matter pathologies associated… ▽ More Mild traumatic brain injury (mTBI) is a complex syndrome that affects up to 600 per 100,000 individuals, with a particular concentration among military personnel. About half of all mTBI patients experience a diverse array of chronic symptoms which persist long after the acute injury. Hence, there is an urgent need for better understanding of the white matter and gray matter pathologies associated with mTBI to map which specific brain systems are impacted and identify courses of intervention. Previous works have linked mTBI to disruptions in white matter pathways and cortical surface abnormalities. Herein, we examine these hypothesized links in an exploratory study of joint structural connectivity and cortical surface changes associated with mTBI and its chronic symptoms. Briefly, we consider a cohort of 12 mTBI and 26 control subjects. A set of 588 cortical surface metrics and 4,753 structural connectivity metrics were extracted from cortical surface regions and diffusion weighted magnetic resonance imaging in each subject. Principal component analysis (PCA) was used to reduce the dimensionality of each metric set. We then applied independent component analysis (ICA) both to each PCA space individually and together in a joint ICA approach. We identified a stable independent component across the connectivity-only and joint ICAs which presented significant group differences in subject loadings (p<0.05, corrected). Additionally, we found that two mTBI symptoms, slowed thinking and forgetfulness, were significantly correlated (p<0.05, corrected) with mTBI subject loadings in a surface-only ICA. These surface-only loadings captured an increase in bilateral cortical thickness. △ Less

Submitted 15 December, 2020; v1 submitted 18 November, 2020; originally announced December 2020.

Comments: To be published in Proc SPIE Int Soc Opt Eng. 2021 Feb

Shape Coexistence at Zero Spin in 64Ni Driven by the Monopole Tensor Interaction

Authors: N. Mărginean, D. Little, Y. Tsunoda, S. Leoni, R. V. F. Janssens, B. Fornal, T. Otsuka, C. Michelagnoli, L. Stan, F. C. L. Crespi, C. Costache, R. Lica, M. Sferrazza, A. Turturica, A. D. Ayangeakaa, K. Auranen, M. Barani, P. C. Bender, S. Bottoni, M. Boromiza, A. Bracco, S. Călinescu, C. M. Campbell, M. P. Carpenter, P. Chowdhury , et al. (53 additional authors not shown)

Abstract: The low-spin structure of the semimagic 64Ni nucleus has been considerably expanded: combining four experiments, several 0+ and 2+ excited states were identified below 4.5 MeV, and their properties established. The Monte Carlo shell model accounts for the results and unveils an unexpectedly complex landscape of coexisting shapes: a prolate 0+ excitation is located at a surprisingly high energy (34… ▽ More The low-spin structure of the semimagic 64Ni nucleus has been considerably expanded: combining four experiments, several 0+ and 2+ excited states were identified below 4.5 MeV, and their properties established. The Monte Carlo shell model accounts for the results and unveils an unexpectedly complex landscape of coexisting shapes: a prolate 0+ excitation is located at a surprisingly high energy (3463 keV), with a collective 2+ state 286 keV above it, the first such observation in Ni isotopes. The evolution in excitation energy of the prolate minimum across the neutron N = 40 subshell gap highlights the impact of the monopole interaction and its variation in strength with N. △ Less

Submitted 11 August, 2020; originally announced August 2020.

Journal ref: Physical Review Letters 125, 102502 (2020)

The Core of $^{25}$F in the Rotational Model

Authors: A. O. Macchiavelli, R. M. Clark, H. L. Crawford, P. Fallon, I. Y. Lee, C. Morse, C. M. Campbell, M. Cromaz, C. Santamaria

Abstract: In a recent experiment, carried out at RIBF/RIKEN, the $^{25}$F$(p,2p)$$^{24}$O reaction was studied at 270 MeV/A in inverse kinematics. Derived spectroscopic factors suggest that the effective core of $^{25}$F significantly differs from a free $^{24}$O nucleus. We interpret these results within the Particle-Rotor Model and show that the experimental level scheme of $^{25}$F can be understood in t… ▽ More In a recent experiment, carried out at RIBF/RIKEN, the $^{25}$F$(p,2p)$$^{24}$O reaction was studied at 270 MeV/A in inverse kinematics. Derived spectroscopic factors suggest that the effective core of $^{25}$F significantly differs from a free $^{24}$O nucleus. We interpret these results within the Particle-Rotor Model and show that the experimental level scheme of $^{25}$F can be understood in the rotation-aligned coupling scheme, with its $5/2^+_1$ ground state as the band-head of a decoupled band. The excitation energies of the observed $1/2_1^+$ and $9/2_1^+$ states correlate strongly with the rotational energy of the effective core, seen by the odd proton, and allow us to estimate its $2^+$ energy at $\approx$ 3.2 MeV and a moderate quadrupole deformation, $εいぷしろん_2 \approx 0.15$. The measured fragmentation of the $πぱいd_{5/2}$ single-particle strength is discussed and some further experiments suggested. △ Less

Submitted 3 August, 2020; originally announced August 2020.

Comments: 5 pages, 3 figures. Submitted to Physical Review C

Search for the $1/2^+$ intruder state in $^{35}$P

Authors: M. Salathe, H. L. Crawford, A. O. Macchiavelli, B. P. Kay, C. R. Hoffman, A. D. Ayangeakaa, C. M. Campbell, R. M. Clark, M. Cromaz, P. Fallon, M. D. Jones, S. A. Kuvin, J. Sethi, M. Wiedeking, J. R. Winkelbauer, A. H. Wuosmaa

Abstract: The excitation energy of deformed intruder states (specifically the 2p2h bandhead) as a function of proton number $Z$ along $N=20$ is of interest both in terms of better understanding the evolution of nuclear structure between spherical $^{40}$Ca and the Island of Inversion nuclei, and for benchmarking theoretical descriptions in this region. At the center of the $N=20$ Island of Inversion, the np… ▽ More The excitation energy of deformed intruder states (specifically the 2p2h bandhead) as a function of proton number $Z$ along $N=20$ is of interest both in terms of better understanding the evolution of nuclear structure between spherical $^{40}$Ca and the Island of Inversion nuclei, and for benchmarking theoretical descriptions in this region. At the center of the $N=20$ Island of Inversion, the npnh (where n=2,4,6) neutron excitations across a diminished $N=20$ gap result in deformed and collective ground states, as observed in $^{32}$Mg. In heavier isotones, npnh excitations do not dominate in the ground states, but are present in the relatively low-lying level schemes. With the aim of identifying the expected 2p2h$\otimes\mathrm{s}_{1/2^+}$ state in $^{35}$P, the only $N=20$ isotone for which the neutron 2p2h excitation bandhead has not yet been identified, the $^{36}$S(d,$^3$He)$^{35}$P reaction has been revisited in inverse kinematics with the HELical Orbit Spectrometer (HELIOS) at the Argonne Tandem Linac Accelerator System (ATLAS). While a candidate state has not been located, an upper limit for the transfer reaction cross-section to populate such a configuration within a 2.5 to 3.6\,MeV energy range, provides a stringent constraint on the wavefunction compositions in both $^{36}$S and $^{35}$P. △ Less

Submitted 16 December, 2020; v1 submitted 28 July, 2020; originally announced July 2020.

Journal ref: Phys. Rev. C 102, 064317 (2020)

Generalizing Variational Autoencoders with Hierarchical Empirical Bayes

Authors: Wei Cheng, Gregory Darnell, Sohini Ramachandran, Lorin Crawford

Abstract: Variational Autoencoders (VAEs) have experienced recent success as data-generating models by using simple architectures that do not require significant fine-tuning of hyperparameters. However, VAEs are known to suffer from over-regularization which can lead to failure to escape local maxima. This phenomenon, known as posterior collapse, prevents learning a meaningful latent encoding of the data. R… ▽ More Variational Autoencoders (VAEs) have experienced recent success as data-generating models by using simple architectures that do not require significant fine-tuning of hyperparameters. However, VAEs are known to suffer from over-regularization which can lead to failure to escape local maxima. This phenomenon, known as posterior collapse, prevents learning a meaningful latent encoding of the data. Recent methods have mitigated this issue by deterministically moment-matching an aggregated posterior distribution to an aggregate prior. However, abandoning a probabilistic framework (and thus relying on point estimates) can both lead to a discontinuous latent space and generate unrealistic samples. Here we present Hierarchical Empirical Bayes Autoencoder (HEBAE), a computationally stable framework for probabilistic generative models. Our key contributions are two-fold. First, we make gains by placing a hierarchical prior over the encoding distribution, enabling us to adaptively balance the trade-off between minimizing the reconstruction loss function and avoiding over-regularization. Second, we show that assuming a general dependency structure between variables in the latent space produces better convergence onto the mean-field assumption for improved posterior inference. Overall, HEBAE is more robust to a wide-range of hyperparameter initializations than an analogous VAE. Using data from MNIST and CelebA, we illustrate the ability of HEBAE to generate higher quality samples based on FID score than existing autoencoder-based approaches. △ Less

Submitted 20 July, 2020; originally announced July 2020.

Comments: 13 pages

Prompt-delayed $γがんま$-ray spectroscopy of neutron-rich $^{119,121}$In isotopes

Authors: S. Biswas, A. Lemasson, M. Rejmund, A. Navin, Y. H. Kim, C. Michelagnoli, I. Stefan, R. Banik, P. Bednarczyk, Soumik Bhattacharya, S. Bhattacharyya, E. Clément, H. L. Crawford, G. de France, P. Fallon, G. Frémont, J. Goupil, B. Jacquot, H. J. Li, J. Ljungvall, A. Maj, L. Ménager, V. Morel, R. Palit, R. M. Pérez-Vidal , et al. (1 additional authors not shown)

Abstract: The fusion and transfer induced fission reaction $^{9}$Be($^{238}$U,~f) with 6.2 MeV/u beam energy, using a unique setup consisting of AGATA, VAMOS++ and EXOGAM detectors, was used to populate through the fission process and study the neutron-rich $^{119,121}$In isotopes. This setup enabled the prompt-delayed $γがんま$-ray spectroscopy of isotopes in the time range of $100~\rm{ns} - 200~μみゅー\rm{s}$. In the… ▽ More The fusion and transfer induced fission reaction $^{9}$Be($^{238}$U,~f) with 6.2 MeV/u beam energy, using a unique setup consisting of AGATA, VAMOS++ and EXOGAM detectors, was used to populate through the fission process and study the neutron-rich $^{119,121}$In isotopes. This setup enabled the prompt-delayed $γがんま$-ray spectroscopy of isotopes in the time range of $100~\rm{ns} - 200~μみゅー\rm{s}$. In the odd-$A$ $^{119,121}$In isotopes, indications of a short half-life $19/2^{-}$ isomeric state, in addition to the previously known $25/2^{+}$ isomeric state, were observed from the present data. Further, new prompt transitions above the $25/2^{+}$ isomer in $^{121}$In were identified along with reevaluation of its half-life. The experimental data were compared with the theoretical results obtained in the framework of large-scale shell-model calculations in a restricted model space. The $\langle πぱいg_{9/2} νにゅーh_{11/2};I \arrowvert \hat{\mathcal{H}}\arrowvert πぱいg_{9/2} νにゅーh_{11/2};I\rangle$ two-body matrix elements of residual interaction were modified to explain the excitation energies and the $B(E2)$ transition probabilities in the neutron-rich In isotopes. The (i) decreasing trend of $E(29/2^{+}) - E(25/2^{+})$ in odd-In (with dominant configuration $πぱいg_{9/2}^{-1}νにゅーh_{11/2}^{-2}$ and maximum aligned spin of $29/2^{+}$) and (ii) increasing trend of $E(27/2^{+}) - E(23/2^{+})$ in odd-Sb (with dominant configuration $πぱいg_{7/2}^{+1}νにゅーh_{11/2}^{-2}$ and maximum aligned spin of $27/2^{+}$) with increasing neutron number could be understood as a consequence of hole-hole and particle-hole interactions, respectively. △ Less

Submitted 20 July, 2020; originally announced July 2020.

Comments: Accepted for publication in Phys. Rev. C

Journal ref: Phys. Rev. C 102 014326 (2020)

Modeling R Coronae Borealis Stars: Effects of He-Burning Shell Temperature and Metallicity

Authors: Courtney L. Crawford, Geoffrey C. Clayton, Bradley Munson, Emmanouil Chatzopoulos, Juhan Frank

Abstract: The R Coronae Borealis (RCB) stars are extremely hydrogen-deficient carbon stars which produce large amounts of dust, causing sudden deep declines in brightness. They are believed to be formed primarily through white dwarf mergers. In this paper, we use MESA to investigate how post-merger objects with a range of initial He-burning shell temperatures from 2.1 - 5.4 *10^8 K with solar and subsolar m… ▽ More The R Coronae Borealis (RCB) stars are extremely hydrogen-deficient carbon stars which produce large amounts of dust, causing sudden deep declines in brightness. They are believed to be formed primarily through white dwarf mergers. In this paper, we use MESA to investigate how post-merger objects with a range of initial He-burning shell temperatures from 2.1 - 5.4 *10^8 K with solar and subsolar metallicities evolve into RCB stars. The most successful model of these has subsolar metallicity and an initial temperature near 3 *10^8 K. We find a strong dependence on initial He-burning shell temperature for surface abundances of elements involved in the CNO cycle, as well as differences in effective temperature and radius of RCBs. Elements involved in nucleosynthesis present around 1 dex diminished surface abundances in the 10% solar metallicity models, with the exception of carbon and lithium which are discussed in detail. Models with subsolar metallicities also exhibit longer lifetimes than their solar counterparts. Additionally, we find that convective mixing of the burned material occurs only in the first few years of post-merger evolution, after which the surface abundances are constant during and after the RCB phase, providing evidence for why these stars show a strong enhancement of partial He-burning products. △ Less

Submitted 17 August, 2020; v1 submitted 6 July, 2020; originally announced July 2020.

Comments: 13 pages, 12 figures, accepted to MNRAS on August 15, 2020

Journal ref: MNRAS, Volume 498, Issue 2, October 2020, Pages 2912-2924

Extending the Southern Shore of the Island of Inversion to $^{28}$F

Authors: A. Revel, O. Sorlin, F. M. Marques, Y. Kondo, J. Kahlbow, T. Nakamura, N. A. Orr, F. Nowacki, J. A. Tostevin, C. X. Yuan, 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. (67 additional authors not shown)

Abstract: Detailed spectroscopy of the neutron-unbound nucleus $^{28}$F has been performed for the first time following proton/neutron removal from $^{29}$Ne/$^{29}$F beams at energies around 230 MeV/nucleon. The invariant-mass spectra were reconstructed for both the $^{27}$F$^{(*)}+n$ and $^{26}$F$^{(*)}+2n$ coincidences and revealed a series of well-defined resonances. A near-threshold state was observed… ▽ More Detailed spectroscopy of the neutron-unbound nucleus $^{28}$F has been performed for the first time following proton/neutron removal from $^{29}$Ne/$^{29}$F beams at energies around 230 MeV/nucleon. The invariant-mass spectra were reconstructed for both the $^{27}$F$^{(*)}+n$ and $^{26}$F$^{(*)}+2n$ coincidences and revealed a series of well-defined resonances. A near-threshold state was observed in both reactions and is identified as the $^{28}$F ground state, with $S_n(^{28}$F$)=-199(6)$ keV, while analysis of the $2n$ decay channel allowed a considerably improved $S_n(^{27}$F$)=1620(60)$ keV to be deduced. Comparison with shell-model predictions and eikonal-model reaction calculations have allowed spin-parity assignments to be proposed for some of the lower-lying levels of $^{28}$F. Importantly, in the case of the ground state, the reconstructed $^{27}$F$+n$ momentum distribution following neutron removal from $^{29}$F indicates that it arises mainly from the $1p_{3/2}$ neutron intruder configuration. This demonstrates that the island of inversion around $N=20$ includes $^{28}$F, and most probably $^{29}$F, and suggests that $^{28}$O is not doubly magic. △ Less

Submitted 2 April, 2020; originally announced April 2020.

Comments: 6 pages, 3 figures, Supplemental Material, accepted for publication in Physical Review Letters

ADC Nonlinearity Correction for the MAJORANA DEMONSTRATOR

Authors: N. Abgrall, J. M. Allmond, I. J. Arnquist, F. T. Avignone III, A. S. Barabash, C. J. Barton, F. E. Bertrand, B. Bos, M. Busch, M. Buuck, T. S. Caldwell, C. M. Campbell, Y-D. Chan, C. D. Christofferson, P. -H. Chu, M. L. Clark, H. L. Crawford, C. Cuesta, J. A. Detwiler, A. Drobizhev, D. W. Edwins, Yu. Efremenko, H. Ejiri, S. R. Elliott, T. Gilliss , et al. (42 additional authors not shown)

Abstract: Imperfections in analog-to-digital conversion (ADC) cannot be ignored when signal digitization requirements demand both wide dynamic range and high resolution, as is the case for the Majorana Demonstrator 76Ge neutrinoless double beta decay search. Enabling the experiment's high-resolution spectral analysis and efficient pulse shape discrimination required careful measurement and correction of ADC… ▽ More Imperfections in analog-to-digital conversion (ADC) cannot be ignored when signal digitization requirements demand both wide dynamic range and high resolution, as is the case for the Majorana Demonstrator 76Ge neutrinoless double beta decay search. Enabling the experiment's high-resolution spectral analysis and efficient pulse shape discrimination required careful measurement and correction of ADC nonlinearites. A simple measurement protocol was developed that did not require sophisticated equipment or lengthy data taking campaigns. A slope-dependent hysteresis was observed and characterized. A correction applied to digitized waveforms prior to signal processing reduced the differential and integral nonlinearites by an order of magnitude, eliminating these as dominant contributions to the systematic energy uncertainty at the double-beta decay Q value. △ Less

Submitted 24 March, 2021; v1 submitted 4 March, 2020; originally announced March 2020.

Report number: LA-UR-20-21663

Journal ref: IEEE Transactions on Nuclear Science, vol. 68, no. 3, pp. 359-367, March 2021

Probing the role of proton cross-shell excitations in 70Ni using nucleon knockout reactions

Authors: B. Elman, A. Gade, R. V. F. Janssens, A. D. Ayangeakaa, D. Bazin, J. Belarge, P. C. Bender, B. A. Brown, C. M. Campbell, M. P. Carpenter, H. L. Crawford, B. P. Crider, P. Fallon, A. M. Forney, J. Harker, S. N. Liddick, B. Longfellow, E. Lunderberg, C. J. Prokop, J. Sethi, R. Taniuchi, W. B. Walters, D. Weisshaar, S. Zhu

Abstract: The neutron-rich Ni isotopes have attracted attention in recent years due to the occurrence of shape or configuration coexistence. We report on the difference in population of excited final states in 70Ni following gamma-ray tagged one-proton, one-neutron, and two-proton knockout from 71Cu, 71Ni, and 72Zn rare-isotope beams, respectively. Using variations observed in the relative transition intens… ▽ More The neutron-rich Ni isotopes have attracted attention in recent years due to the occurrence of shape or configuration coexistence. We report on the difference in population of excited final states in 70Ni following gamma-ray tagged one-proton, one-neutron, and two-proton knockout from 71Cu, 71Ni, and 72Zn rare-isotope beams, respectively. Using variations observed in the relative transition intensities, signaling the changed population of specific final states in the different reactions, the role of neutron and proton configurations in excited states of 70Ni is probed schematically, with the goal of identifying those that carry, as leading configuration, proton excitations across the Z = 28 shell closure. Such states are suggested in the literature to form a collective structure associated with prolate deformation. Adding to the body of knowledge for 70Ni, 29 new transitions are reported, of which 15 are placed in its level scheme. △ Less

Submitted 15 September, 2019; originally announced September 2019.

Comments: accepted for publication in Phys. Rev. C

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

Evidence for rigid triaxial deformation in $^{76}$Ge from a model-independent analysis

Authors: A. D. Ayangeakaa, R. V. F. Janssens, S. Zhu, D. Little, J. Henderson, C. Y. Wu, D. J. Hartley, M. Albers, K. Auranen, B. Bucher, M. P. Carpenter, P. Chowdhury, D. Cline, H. L. Crawford, P. Fallon, A. M. Forney, A. Gade, A. B. Hayes, F. G. Kondev, Krishichayan, T. Lauritsen, J. Li, A. O. Macchiavelli, D. Rhodes, D. Seweryniak , et al. (3 additional authors not shown)

Abstract: An extensive, model-independent analysis of the nature of triaxial deformation in $^{76}$Ge, a candidate for neutrinoless double-beta ($0νββ$) decay, was carried out following multi-step Coulomb excitation. Shape parameters deduced on the basis of a rotational-invariant sum-rule analysis provided considerable insight into the underlying collectivity of the ground-state and $γがんま$ bands. Both sequence… ▽ More An extensive, model-independent analysis of the nature of triaxial deformation in $^{76}$Ge, a candidate for neutrinoless double-beta ($0νββ$) decay, was carried out following multi-step Coulomb excitation. Shape parameters deduced on the basis of a rotational-invariant sum-rule analysis provided considerable insight into the underlying collectivity of the ground-state and $γがんま$ bands. Both sequences were determined to be characterized by the same $βべーた$ and $γがんま$ deformation parameter values. In addition, compelling evidence for low-spin, rigid triaxial deformation in $^{76}$Ge was obtained for the first time from the analysis of the statistical fluctuations of the quadrupole asymmetry deduced from the measured $E2$ matrix elements. These newly determined shape parameters are important input and constraints for calculations aimed at providing, with suitable accuracy, the nuclear matrix elements relevant to $0νββ$. △ Less

Submitted 7 September, 2019; originally announced September 2019.

Comments: 6 pages, 5 figures

Journal ref: Phys. Rev. Lett. 123, 102501 (2019)

Experimental Constraint on Stellar Electron-Capture Rates from the ${}^{88}\text{Sr}(t,{}^{3}\text{He}+γがんま){}^{88}\text{Rb}$ reaction at 115 MeV/u

Authors: J. C. Zamora, R. G. T. Zegers, Sam M. Austin, D. Bazin, B. A. Brown, P. C. Bender, H. L. Crawford, J. Engel, A. Falduto, A. Gade, P. Gastis, B. Gao, T. Ginter, C. J. Guess, S. Lipschutz, B. Longfellow, A. O. Macchiavelli, K. Miki, E. Ney, S. Noji, J. Pereira, J. Schmitt, C. Sullivan, R. Titus, D. Weisshaar

Abstract: The Gamow-Teller strength distribution from ${}^{88}$Sr was extracted from a $(t,{}^{3}\text{He}+γがんま)$ experiment at 115 MeV/$u$ to constrain estimates for the electron-capture rates on nuclei around $N=50$, between and including $^{78}$Ni and $^{88}$Sr, which are important for the late evolution of core-collapse supernovae. The observed strength below an excitation energy of 8 MeV was consistent wi… ▽ More The Gamow-Teller strength distribution from ${}^{88}$Sr was extracted from a $(t,{}^{3}\text{He}+γがんま)$ experiment at 115 MeV/$u$ to constrain estimates for the electron-capture rates on nuclei around $N=50$, between and including $^{78}$Ni and $^{88}$Sr, which are important for the late evolution of core-collapse supernovae. The observed strength below an excitation energy of 8 MeV was consistent with zero and below 10 MeV amounted to $0.1\pm0.05$. Except for a very-weak transition that could come from the 2.231-MeV $1^{+}$ state, no $γがんま$ lines that could be associated with the decay of known $1^{+}$ states were identified. The derived electron-capture rate from the measured strength distribution is more than an order of magnitude smaller than rates based on the single-state approximation presently used in astrophysical simulations for most nuclei near $N=50$. Rates based on shell-model and quasiparticle random-phase approximation calculations that account for Pauli blocking and core-polarization effects provide better estimates than the single-state approximation, although a relatively strong transition to the first $1^{+}$ state in $^{88}$Rb is not observed in the data. Pauli unblocking effects due to high stellar temperatures could partially counter the low electron-capture rates. The new data serves as a zero-temperature benchmark for constraining models used to estimate such effects. △ Less

Submitted 13 June, 2019; originally announced June 2019.

Effects of one valence proton on seniority and angular momentum of neutrons in neutron-rich $^{122-131}$Sb$_{51}$ isotopes

Authors: S. Biswas, A. Lemasson, M. Rejmund, A. Navin, Y. H. Kim, C. Michelagnoli, I. Stefan, R. Banik, P. Bednarczyk, S. Bhattacharya, S. Bhattacharyya, E. Clément, H. L. Crawford, G. de France, P. Fallon, G. Frémont, J. Goupil, B. Jacquot, H. J. Li, J. Ljungvall, A. Maj, L. Ménager, V. Morel, R. Palit, R. M. Pérez-Vidal , et al. (36 additional authors not shown)

Abstract: The neutron-rich $^{122-131}$Sb isotopes were produced as fission fragments in the reaction $^{9}$Be($^{238}$U,~f) with 6.2 MeV/u beam energy. An unique setup, consisting of AGATA, VAMOS++ and EXOGAM detectors, was used which enabled the prompt-delayed gamma-ray ($γがんま$) spectroscopy of fission fragments in the time range of 100 ns - 200 $μみゅー$s. New isomers, prompt and delayed transitions were establis… ▽ More The neutron-rich $^{122-131}$Sb isotopes were produced as fission fragments in the reaction $^{9}$Be($^{238}$U,~f) with 6.2 MeV/u beam energy. An unique setup, consisting of AGATA, VAMOS++ and EXOGAM detectors, was used which enabled the prompt-delayed gamma-ray ($γがんま$) spectroscopy of fission fragments in the time range of 100 ns - 200 $μみゅー$s. New isomers, prompt and delayed transitions were established in the even-A $^{122-130}$Sb isotopes. In the odd-A $^{123-131}$Sb isotopes, new prompt and delayed $γがんま$-ray transitions were identified, in addition to the confirmation of the previously known isomers. The half-lives of the isomeric states and the $B(E2)$ transition probabilities of the observed transitions depopulating these isomers were extracted. The experimental data was compared with the theoretical results obtained in the framework of Large-Scale Shell-Model (LSSM) calculations in a restricted model space. Modifications of several components of the shell model interaction were introduced to obtain a consistent agreement with the excitation energies and the $B(E2)$ transition probabilities in neutron-rich Sn and Sb isotopes. The isomeric configurations in Sn and Sb were found to be relatively pure. Further, the calculations revealed that the presence of a single valence proton, mainly in the $g_{7/2}$ orbital in Sb isotopes, leads to significant mixing (due to the $νにゅーπぱい$ interaction) of: (i) the neutron seniorities ($\upsilon_νにゅー$) and (ii) the neutron angular momentum ($I_νにゅー$). The above features have a weak impact on the excitation energies, but have an important impact on the $B(E2)$ transition probabilities. In addition, a constancy of the relative excitation energies irrespective of neutron seniority and neutron number in Sn and Sb was observed. △ Less

Submitted 4 June, 2019; originally announced June 2019.

Comments: 23 pages, 18 figures, 10 tables

Journal ref: Physical Review C 99, 064302 (2018)

Is the structure of 42Si understood?

Authors: A. Gade, B. A. Brown, J. A. Tostevin, D. Bazin, P. C. Bender, C. M. Campbell, H. L. Crawford, B. Elman, K. W. Kemper, B. Longfellow, E. Lunderberg, D. Rhodes, D. Weisshaar

Abstract: A more detailed test of the implementation of nuclear forces that drive shell evolution in the pivotal nucleus \nuc{42}{Si} -- going beyond earlier comparisons of excited-state energies -- is important. The two leading shell-model effective interactions, SDPF-MU and SDPF-U-Si, both of which reproduce the low-lying \nuc{42}{Si}($2^+_1$) energy, but whose predictions for other observables differ sig… ▽ More A more detailed test of the implementation of nuclear forces that drive shell evolution in the pivotal nucleus \nuc{42}{Si} -- going beyond earlier comparisons of excited-state energies -- is important. The two leading shell-model effective interactions, SDPF-MU and SDPF-U-Si, both of which reproduce the low-lying \nuc{42}{Si}($2^+_1$) energy, but whose predictions for other observables differ significantly, are interrogated by the population of states in neutron-rich \nuc{42}{Si} with a one-proton removal reaction from \nuc{43}{P} projectiles at 81~MeV/nucleon. The measured cross sections to the individual \nuc{42}{Si} final states are compared to calculations that combine eikonal reaction dynamics with these shell-model nuclear structure overlaps. The differences in the two shell-model descriptions are examined and linked to predicted low-lying excited $0^+$ states and shape coexistence. Based on the present data, which are in better agreement with the SDPF-MU calculations, the state observed at 2150(13)~keV in \nuc{42}{Si} is proposed to be the ($0^+_2$) level. △ Less

Submitted 10 May, 2019; originally announced May 2019.

Comments: accepted in Physical Review Letters

Journal ref: Phys. Rev. Lett. 122, 222501 (2019)

Interpreting Deep Neural Networks Through Variable Importance

Authors: Jonathan Ish-Horowicz, Dana Udwin, Seth Flaxman, Sarah Filippi, Lorin Crawford

Abstract: While the success of deep neural networks (DNNs) is well-established across a variety of domains, our ability to explain and interpret these methods is limited. Unlike previously proposed local methods which try to explain particular classification decisions, we focus on global interpretability and ask a universally applicable question: given a trained model, which features are the most important?… ▽ More While the success of deep neural networks (DNNs) is well-established across a variety of domains, our ability to explain and interpret these methods is limited. Unlike previously proposed local methods which try to explain particular classification decisions, we focus on global interpretability and ask a universally applicable question: given a trained model, which features are the most important? In the context of neural networks, a feature is rarely important on its own, so our strategy is specifically designed to leverage partial covariance structures and incorporate variable dependence into feature ranking. Our methodological contributions in this paper are two-fold. First, we propose an effect size analogue for DNNs that is appropriate for applications with highly collinear predictors (ubiquitous in computer vision). Second, we extend the recently proposed "RelATive cEntrality" (RATE) measure (Crawford et al., 2019) to the Bayesian deep learning setting. RATE applies an information theoretic criterion to the posterior distribution of effect sizes to assess feature significance. We apply our framework to three broad application areas: computer vision, natural language processing, and social science. △ Less

Submitted 28 April, 2020; v1 submitted 28 January, 2019; originally announced January 2019.

The Age of the Galactic Stellar Halo from Gaia White Dwarfs

Authors: Mukremin Kilic, P. Bergeron, Kyra Dame, N. C. Hambly, N. Rowell, Courtney L. Crawford

Abstract: We use 156 044 white dwarf candidates with $\geq5σしぐま$ significant parallax measurements from the Gaia mission to measure the velocity dispersion of the Galactic disc; $(σしぐま_U,σしぐま_V,σしぐま_W) = (30.8, 23.9, 20.0)$ km s$^{-1}$. We identify 142 objects that are inconsistent with disc membership at the $>5σしぐま$ level. This is the largest sample of field halo white dwarfs identified to date. We perform a detailed mo… ▽ More We use 156 044 white dwarf candidates with $\geq5σしぐま$ significant parallax measurements from the Gaia mission to measure the velocity dispersion of the Galactic disc; $(σしぐま_U,σしぐま_V,σしぐま_W) = (30.8, 23.9, 20.0)$ km s$^{-1}$. We identify 142 objects that are inconsistent with disc membership at the $>5σしぐま$ level. This is the largest sample of field halo white dwarfs identified to date. We perform a detailed model atmosphere analysis using optical and near-infrared photometry and parallaxes to constrain the mass and cooling age of each white dwarf. The white dwarf cooling ages of our targets range from 7 Myr for J1657+2056 to 10.3 Gyr for J1049-7400. The latter provides a firm lower limit of 10.3 Gyr for the age of the inner halo based on the well-understood physics of white dwarfs. Including the pre-white dwarf evolutionary lifetimes, and limiting our sample to the recently formed white dwarfs with cooling ages of $<500$ Myr, we estimate an age of $10.9 \pm 0.4$ Gyr (internal errors only) for the Galactic inner halo. The coolest white dwarfs in our sample also give similar results. For example, J1049-7400 has a total age of 10.9-11.1 Gyr. Our age measurements are consistent with other measurements of the age of the inner halo, including the white dwarf based measurements of the globular clusters M4, NGC 6397, and 47 Tuc. △ Less

Submitted 8 October, 2018; originally announced October 2018.

Comments: MNRAS, in press

Understanding the Low-Energy Enhancement of the $γがんま$-ray Strength Function of $^{56}$Fe

Authors: M. D. Jones, A. O. Macchiavelli, M. Wiedeking, L. A. Bernstein, H. L. Crawford, C. M. Campbell, R. M. Clark, M. Cromaz, P. Fallon, I. Y. Lee, M. Salathe, A. Wiens, A. D. Ayangeakaa, D. L. Bleuel, S. Bottoni, M. P. Carpenter, H. M. Davids, J. Elson, A. Görgen, M. Guttormsen, R. V. F. Janssens, J. E. Kinnison, L. Kirsch, A. C. Larsen, T. Lauritsen , et al. (5 additional authors not shown)

Abstract: A model-independent technique was used to determine the $γがんま$-ray Strength Function ($γがんま$SF) of $^{56}$Fe down to $γがんま$-ray energies less than 1 MeV for the first time with GRETINA using the $(p,p')$ reaction at 16 MeV. No difference was observed in the energy dependence of the $γがんま$SF built on $2^{+}$ and $4^{+}$ final states, supporting the Brink hypothesis. In addition, angular distribution and polari… ▽ More A model-independent technique was used to determine the $γがんま$-ray Strength Function ($γがんま$SF) of $^{56}$Fe down to $γがんま$-ray energies less than 1 MeV for the first time with GRETINA using the $(p,p')$ reaction at 16 MeV. No difference was observed in the energy dependence of the $γがんま$SF built on $2^{+}$ and $4^{+}$ final states, supporting the Brink hypothesis. In addition, angular distribution and polarization measurements were performed. The angular distributions are consistent with dipole radiation. The polarization results show a small bias towards magnetic character in the region of the enhancement. △ Less

Submitted 29 January, 2018; originally announced January 2018.

Journal ref: Phys. Rev. C 97, 024327 (2018)

Variable Prioritization in Nonlinear Black Box Methods: A Genetic Association Case Study

Authors: Lorin Crawford, Seth R. Flaxman, Daniel E. Runcie, Mike West

Abstract: The central aim in this paper is to address variable selection questions in nonlinear and nonparametric regression. Motivated by statistical genetics, where nonlinear interactions are of particular interest, we introduce a novel and interpretable way to summarize the relative importance of predictor variables. Methodologically, we develop the "RelATive cEntrality" (RATE) measure to prioritize cand… ▽ More The central aim in this paper is to address variable selection questions in nonlinear and nonparametric regression. Motivated by statistical genetics, where nonlinear interactions are of particular interest, we introduce a novel and interpretable way to summarize the relative importance of predictor variables. Methodologically, we develop the "RelATive cEntrality" (RATE) measure to prioritize candidate genetic variants that are not just marginally important, but whose associations also stem from significant covarying relationships with other variants in the data. We illustrate RATE through Bayesian Gaussian process regression, but the methodological innovations apply to other "black box" methods. It is known that nonlinear models often exhibit greater predictive accuracy than linear models, particularly for phenotypes generated by complex genetic architectures. With detailed simulations and two real data association mapping studies, we show that applying RATE enables an explanation for this improved performance. △ Less

Submitted 26 August, 2018; v1 submitted 22 January, 2018; originally announced January 2018.

Comments: 28 pages, 5 figures, 1 tables; Supplementary Material

Tropical Sufficient Statistics for Persistent Homology

Authors: Anthea Monod, Sara Kališnik, Juan Ángel Patiño-Galindo, Lorin Crawford

Abstract: We show that an embedding in Euclidean space based on tropical geometry generates stable sufficient statistics for barcodes. In topological data analysis, barcodes are multiscale summaries of algebraic topological characteristics that capture the `shape' of data; however, in practice, they have complex structures that make them difficult to use in statistical settings. The sufficiency result prese… ▽ More We show that an embedding in Euclidean space based on tropical geometry generates stable sufficient statistics for barcodes. In topological data analysis, barcodes are multiscale summaries of algebraic topological characteristics that capture the `shape' of data; however, in practice, they have complex structures that make them difficult to use in statistical settings. The sufficiency result presented in this work allows for classical probability distributions to be assumed on the tropical geometric representation of barcodes. This makes a variety of parametric statistical inference methods amenable to barcodes, all while maintaining their initial interpretations. More specifically, we show that exponential family distributions may be assumed, and that likelihood functions for persistent homology may be constructed. We conceptually demonstrate sufficiency and illustrate its utility in persistent homology dimensions 0 and 1 with concrete parametric applications to human immunodeficiency virus and avian influenza data. △ Less

Submitted 30 June, 2019; v1 submitted 8 September, 2017; originally announced September 2017.

Comments: 31 pages, 5 figures

Report number: MPI MIS 66/2017

Journal ref: SIAM Journal on Applied Algebra and Geometry 3 (2), 337-371 (2019)

Structure of 29F in the Rotation-aligned Coupling Scheme of the Particle-Rotor Model

Authors: A. O. Macchiavelli, P. Fallon, H. L. Crawford, C. M. Campbell, R. M. Clark, M. Cromaz, M. D. Jones, I. Y. Lee, M. Salathe

Abstract: Recent results from RIKEN/RIBF on the low-lying level structure of 29F are interpreted within the Particle-Rotor Model. We show that the experimental data can be understood in the Rotation-aligned Coupling Scheme, with the 5/2+ ground state as the bandhead of a decoupled band. In this picture, the energy of the observed 1/2+1 state correlates strongly with the rotational energy of the core and pro… ▽ More Recent results from RIKEN/RIBF on the low-lying level structure of 29F are interpreted within the Particle-Rotor Model. We show that the experimental data can be understood in the Rotation-aligned Coupling Scheme, with the 5/2+ ground state as the bandhead of a decoupled band. In this picture, the energy of the observed 1/2+1 state correlates strongly with the rotational energy of the core and provides an estimate of the 2+ energy in 28O. Our analysis suggest a moderate deformation, epsilon_2 ~ 0.17, and places the 2+ in 28O at ~ 2.4 MeV. △ Less

Submitted 7 September, 2017; originally announced September 2017.

Comments: 3 pages, 3 figures, To be submitted to Phys. Lett. B

Analysis of Spectroscopic Factors in 11Be and 12Be in the Nilsson Strong Coupling Limit

Authors: A. O. Macchiavelli, H. L. Crawford, C. M. Campbell, R. M. Clark, M. Cromaz, P. Fallon, M. D. Jones, I. Y. Lee, M. Salathe

Abstract: Spectroscopic factors in 10Be, 11Be and 12Be, extracted from (d,p), one neutron knockout, and (p,d) reactions are interpreted within the rotational model. Assuming that the ground state and first excited state of 11Be can be associated with the 1/2[220] and 1/2[101] Nilsson levels, the strong coupling limit gives simple expressions that relate the amplitudes of these wavefunctions (in the spherica… ▽ More Spectroscopic factors in 10Be, 11Be and 12Be, extracted from (d,p), one neutron knockout, and (p,d) reactions are interpreted within the rotational model. Assuming that the ground state and first excited state of 11Be can be associated with the 1/2[220] and 1/2[101] Nilsson levels, the strong coupling limit gives simple expressions that relate the amplitudes of these wavefunctions (in the spherical basis) with the measured cross-sections and derived spectroscopic factors. We obtain good agreement with both the measured magnetic moment of the ground state in 11Be and the reaction data. △ Less

Submitted 18 August, 2017; originally announced August 2017.

Comments: 4 pages, 1 figure

Journal ref: Phys. Rev. C 97, 011302 (2018)

$βべーた$-particle energy-summing correction for $βべーた$-delayed proton emission measurements

Authors: Z. Meisel, M. del Santo, H. L. Crawford, R. H. Cyburt, G. F. Grinyer, C. Langer, F. Montes, H. Schatz, K. Smith

Abstract: A common approach to studying $βべーた$-delayed proton emission is to measure the energy of the emitted proton and corresponding nuclear recoil in a double-sided silicon-strip detector (DSSD) after implanting the $βべーた$-delayed proton emitting ($βべーた$p) nucleus. However, in order to extract the proton-decay energy, the measured energy must be corrected for the additional energy implanted in the DSSD by the… ▽ More A common approach to studying $βべーた$-delayed proton emission is to measure the energy of the emitted proton and corresponding nuclear recoil in a double-sided silicon-strip detector (DSSD) after implanting the $βべーた$-delayed proton emitting ($βべーた$p) nucleus. However, in order to extract the proton-decay energy, the measured energy must be corrected for the additional energy implanted in the DSSD by the $βべーた$-particle emitted from the $βべーた$p nucleus, an effect referred to here as $βべーた$-summing. We present an approach to determine an accurate correction for $βべーた$-summing. Our method relies on the determination of the mean implantation depth of the $βべーた$p nucleus within the DSSD by analyzing the shape of the total (proton + recoil + $βべーた$) decay energy distribution shape. We validate this approach with other mean implantation depth measurement techniques that take advantage of energy deposition within DSSDs upstream and downstream of the implantation DSSD. △ Less

Submitted 18 November, 2016; originally announced November 2016.

Comments: Published in Nuclear Instruments and Methods in Physics Research A

Journal ref: Nuclear Instruments and Methods in Physics Research A 844, 45-52 (2017)

Predicting Clinical Outcomes in Glioblastoma: An Application of Topological and Functional Data Analysis

Authors: Lorin Crawford, Anthea Monod, Andrew X. Chen, Sayan Mukherjee, Raúl Rabadán

Abstract: Glioblastoma multiforme (GBM) is an aggressive form of human brain cancer that is under active study in the field of cancer biology. Its rapid progression and the relative time cost of obtaining molecular data make other readily-available forms of data, such as images, an important resource for actionable measures in patients. Our goal is to utilize information given by medical images taken from G… ▽ More Glioblastoma multiforme (GBM) is an aggressive form of human brain cancer that is under active study in the field of cancer biology. Its rapid progression and the relative time cost of obtaining molecular data make other readily-available forms of data, such as images, an important resource for actionable measures in patients. Our goal is to utilize information given by medical images taken from GBM patients in statistical settings. To do this, we design a novel statistic---the smooth Euler characteristic transform (SECT)---that quantifies magnetic resonance images (MRIs) of tumors. Due to its well-defined inner product structure, the SECT can be used in a wider range of functional and nonparametric modeling approaches than other previously proposed topological summary statistics. When applied to a cohort of GBM patients, we find that the SECT is a better predictor of clinical outcomes than both existing tumor shape quantifications and common molecular assays. Specifically, we demonstrate that SECT features alone explain more of the variance in GBM patient survival than gene expression, volumetric features, and morphometric features. The main takeaways from our findings are thus twofold. First, they suggest that images contain valuable information that can play an important role in clinical prognosis and other medical decisions. Second, they show that the SECT is a viable tool for the broader study of medical imaging informatics. △ Less

Submitted 12 September, 2019; v1 submitted 21 November, 2016; originally announced November 2016.

Comments: 30 pages, 9 figures, 1 table

Journal ref: Journal of the American Statistical Association (2019)

Unexpected distribution of $\nu1f_{7/2}$ strength in the calcium isotopes at $N$=30

Authors: H. L. Crawford, A. O. Macchiavelli, P. Fallon, M. Albers, V. M. Bader, D. Bazin, C. M. Campbell, R. M. Clark, M. Cromaz, J. Dilling, A. Gade, A. T. Gallant, J. D. Holt, R. V. F. Janssens, R. Krücken, C. Langer, T. Lauritsen, I. Y. Lee, J. Menéndez, S. Noji, S. Paschalis, F. Recchia, J. Rissanen, A. Schwenk, M. Scott , et al. (8 additional authors not shown)

Abstract: The calcium isotopes have emerged as an important testing ground for new microscopically derived shell-model interactions, and a great deal of focus has been directed toward this region. We investigate the relative spectroscopic strengths associated with $1f_{7/2}$ neutron hole states in $^{47, 49}$Ca following one-neutron knockout reactions from $^{48,50}$Ca. The observed reduction of strength po… ▽ More The calcium isotopes have emerged as an important testing ground for new microscopically derived shell-model interactions, and a great deal of focus has been directed toward this region. We investigate the relative spectroscopic strengths associated with $1f_{7/2}$ neutron hole states in $^{47, 49}$Ca following one-neutron knockout reactions from $^{48,50}$Ca. The observed reduction of strength populating the lowest 7/2$^{-}_{1}$ state in $^{49}$Ca, as compared to $^{47}$Ca, is consistent with the description given by shell-model calculations based on two- and three-nucleon forces in the neutron $pf$ model space, implying a fragmentation of the $l$=3 strength to higher-lying states. The experimental result is inconsistent with both the GXPF1 interaction routinely used in this region of the nuclear chart and with microscopic calculations in an extended model space including the $\nu1g_{9/2}$ orbital. △ Less

Submitted 5 October, 2016; originally announced October 2016.

Comments: 6 pages, 2 figures

Shape coexistence and the role of axial asymmetry in $^{72}$Ge

Authors: A. D. Ayangeakaa, R. V. F. Janssens, C. Y. Wu, J. M. Allmond, J. L. Wood, S. Zhu, M. Albers, S. Almaras-Calderon, B. Bucher, M. P. Carpenter, C. J. Chiara, D. Cline, H. L. Crawford, H. M. David, J. Harker, A. B. Hayes, C. R. Hoffman, B. P. Kay, K. Kolos, A. Korichi, T. Lauritsen, A. O. Macchiavelli, A. Richard, D. Seweryniak, A. Wiens

Abstract: The quadrupole collectivity of low-lying states and the anomalous behavior of the $0^+_2$ and $2^+_3$ levels in $^{72}$Ge are investigated via projectile multi-step Coulomb excitation with GRETINA and CHICO-2. A total of forty six $E2$ and $M1$ matrix elements connecting fourteen low-lying levels were determined using the least-squares search code, gosia. Evidence for triaxiality and shape coexist… ▽ More The quadrupole collectivity of low-lying states and the anomalous behavior of the $0^+_2$ and $2^+_3$ levels in $^{72}$Ge are investigated via projectile multi-step Coulomb excitation with GRETINA and CHICO-2. A total of forty six $E2$ and $M1$ matrix elements connecting fourteen low-lying levels were determined using the least-squares search code, gosia. Evidence for triaxiality and shape coexistence, based on the model-independent shape invariants deduced from the Kumar-Cline sum rule, is presented. These are interpreted using a simple two-state mixing model as well as multistate mixing calculations carried out within the framework of the triaxial rotor model. The results represent a significant milestone towards the understanding of the unusual structure of this nucleus. △ Less

Submitted 14 March, 2016; originally announced March 2016.

Comments: 7 pages and 7 figures

Journal ref: Phys. Lett. B 754, 254 (2016)