39%
.1 r... ‘
‘35....

iv!

..
Mu «Ming?

H.‘ - A
.1... EH?
:3.

4...

n 3 Ii
. x . g
.65....
. ‘I ‘5?
Mums“.
I . .I
31%. i.
.A ntu‘gl."
31.65 I .flfi
A 47131

I: €51}

.x 1”“...
4‘)...‘ \...
‘Ptlii‘flcsi‘v
\ $ ‘ . fl... 2!

 

. fifiafiafiég Lgmwé

129v
s

.~ 11.

1 uv .
é «a if... av

mgr:

 

THESIS
Q.

m LIBRARY
Michigan State
University

This is to certify that the
dissertation entitled

DEVELOPMENT OF THE HIGH VELOCITY TRANSIENT

FIELD TECHNIQUE AND MEASUREMENT OF THE 21*
g FACTORS OF 38"‘08

presented by

ANDREW DAVID DAVIES

has been accepted towards fulfillment
of the requirements for the

degree In PHYSICS

 

 

mflwlzfim/LA

Major Professor’s Signature
S’l‘l I 20 o (a
‘ r

Date

 

MSU is an Affirmative Action/Equal Opportunity Institution

-._.—.—--o-..

PLACE IN RETURN BOX to remove this checkout from your record.
To AVOID FINES return on or before date due.
MAY BE RECALLED with earlier due date if requested.

 

DATE DUE

DATE DUE

DATE DUE

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2/05 p:/ClRC/Date0ue.indd-p.1

 

DEVELOPMENT OF THE HIGH VELOCITY
TRANSIENT FIELD TECHNIQUE AND
MEASUREMENT OF THE 2? g FACTORS OF 38"‘08

By

Andrew David Davies

A DISSERTATION

Submitted to
Michigan State University
in partial fulfillment of the requirements
for the degree of

DOCTOR OF PHILOSOPHY
Department of Physics and Astronomy

2006

ABSTRACT

DEVELOPMENT OF THE HIGH VELOCITY TRANSIENT FIELD
TECHNIQUE AND MEASUREMENT OF THE 2f 9 FACTORS OF 38,408

By

Andrew David Davies

The onset Of deformation in the neutron-rich sulfur isotopes has been investigated,
through a new application of the transient field technique for measuring Short-lived
excited-state 9 factors with fast radioactive beams. The 9 factor (the magnetic dipole
moment divided by the Spin Of the nuclear state) can provide insight into the complex
interplay between single particle configurations and collective behavior in a nuclear
state. In the case Of the neutron-rich sulfur isotopes, a transition from spherical tO
deformed shapes has been inferred from first-excited state energy and transition prob—
ability measurements. However, the underlying causes Of the deformation change have
not. yet been deduced. Knowledge of the g(2+) values can greatly constrain nuclear
models and provide insight into the causes Of the deformation change. However, the
measurement of both the Sign and magnitude Of picosecond excited-state g factors
requires the use of the transient field method, which generally requires ion velocities
Of Z U0 (v0 = c/137), or lower. Due to the high velocity Of fragmentation products,
the conventional transient field technique (which has been extensively developed for
stable, Coulomb-barrier energy beams) could not be applied, and challenges posed by
high-velocity radioactive beams needed to be overcome. A new technique, which com-
bines the transient field method with intermediate-energy Coulomb excitation, has
been developed in order to provide access to g(2+) values in fast radioactive secondary
beams. The g(2+) values Of 388 and 408, with N = 22 and 24, were measured in the
first application Of the method at the NSCL. The description Of the technique, results,

and interpretation within a. shell—model context is presented in this dissertation.

for Debbie

iii

ACKNOWLEDGMENTS

This dissertation would certainly not have been possible without the help and
support of a lot of people. I gratefully acknowledge Michigan State University, the
Department Of Physics and Astronomy, the NSCL, the National Science Foundation,
and the APS Division Of Nuclear Physics for financial support for graduate school
and to travel to conferences. My adviser Paul Mantica always ensured I had the
resources necessary to get this project done, and it was a privilege learning how to
be an experimentalist from him. It was the vision and physical insight Of Andrew
Stuchbery that inspired the entire experiment, and his participation in every facet
Of the project despite being located thousands Of miles away is greatly appreciated.
Warren Rogers ultimately was responsible for my decision to go tO graduate school; he
is a true mentor, and has also had a huge impact through the Conference Experience
for Undergraduates. Warren first introduced me to both Paul and Andrew S. — it’s
a small world!

I thank the members Of my thesis committee, Gregers Hansen, Alex Brown,
Thomas Glasmacher, Bhanu Mahanti, Paul Mantica, Andrew Stuchbery, and Michael
Thoennessen, for their guidance and helpful discussions over the years.

The NSCL technical staff are totally outstanding, and literally dozens Of people
helped make the project a reality. I really need to thank Len Morris and Renan Fontus
for their design work, and the N SCL Machine ShOp for producing nearly everything
we needed. John Priller wrote the magnet control script and John Brandon built
the interface module that was inspired by Kelly Davidson. I thank John Yurkon
for the PMT and for detector design advice, and Al Zeller for running the magnet
calculations. Ron Fox helped enormously with the data acquisition, and Wil Mueller
put in a lot of effort to support the use Of SeGA. I needed to call on many others
in order to build and set everything up prior to the run, and I am truly indebted tO

them all.

iv

Theo Mertzimekis deserves a huge amount Of credit for being truly helpful over
the years, especially during the target test run, and together with Sean and Bryan,
made the ANU measurements possible. Paul Davidson and Anna Wilson went above
and beyond to contribute to the success Of the experiment, and I can’t thank them
enough. Many thanks go to the Gamma Group, especially to Heather for her help
with the SeCA stand, and to so many others for help during the setup and with
taking shifts. The other graduate students here at MSU have been great friends and
it has been a pleasure tO spend all these years here in large part because of them.
I especially should thank Chip, Nate, Heather, Russ, Joe, Josh, Ryan, and Pete for
their friendship and for being great study partners for those first years here. Elaine,
Chandana, and Jeremy have been great friends. I have been lucky tO have been part.
Of a research group here that has great chemistry among its members and has made
this time worthwhile.

I want to thank my family for their love and support throughout college and
graduate school. I thank God for them, for making it possible for me to attend
Westmont in the first place, and for bringing Debbie into my life. Debbie, you have
shown me SO much love and patience while going through graduate school at the same

time, and I love you.

‘7

Contents

1 Introduction 1
1.1 Physics Of Neutron-rich Sulfur Isotopes ................. 2
1.1.1 Previous Experimental and Theoretical Work ......... 3

1.1.2 Open Questions .......................... 7

1.2 Electromagnetic Moments Of Nuclei ................... 10
1.2.1 Magnetic Dipole Moments in the Shell Model ......... 13

1.2.2 g(2+) in Even-even Nuclei .................... 17

1.3 Organization of Dissertation ....................... 20
2 'h'ansient Field Technique 22
2.1 Measurement Principles: Stable Beams ................. 25
2.1.1 Angular Distribution from Spin Alignment ........... 26

2.1.2 Transient Field Strength ..................... 32
2.1.3 Measurement Of Spin Precession ................. 37

2.2 g(2+) Of Radionuclides: New Challenges ................ 39
2.3 g(2+) Of Fast Radioactive Beams: New Technique ........... 42
2.3.1 Intermediate—Energy Coulomb Excitation ........... 43

2.3.2 Transient Field Strength at High Velocity ........... 49
2.3.3 Target Design for Optimum Slowing .............. 53

3 Application at NSCL 58
3.1 Isotope Production ............................ 58
3.2 Multilayer Target ............................. 60
3.2.1 Stopping Powers and In—Beam Test ............... 61

3.3 Coulomb Excitation Of 38’408 Isotopes ................. 64
3.3.1 Segmented Germanium Array .................. 64
3.3.2 Maximum Scattering Angle and Alignment ........... 69
3.3.3 Slope of Angular Distribution and Sensitivity ......... 69

3.4 Phoswich Detector ............................ 70
3.4.1 Design and Construction ..................... 71
3.4.2 Characterization Of Phoswich Detector ............. 77
3.4.3 High Count Rate Behavior .................... 81

3.5 Compact Electromagnet ......................... 83
3.5.1 Design and Construction ..................... 83
3.5.2 Hall Probe Measurements .................... 87

3.5.3 Power Supply Control ...................... 89

vi

3.6 Data Acquisition System ......................... 90

3.6.1 Germanium Detector Electronics ................ 90

3.6.2 Phoswich and PIN Detector Electronics ............. 92

3.6.3 Magnet Interface to Datastream ................. 94

3.6.4 Data Acquisition Software .................... 94

3.7 Layout in N3 Vault ............................ 95

4 Results and Conclusions 99
4.1 Particle Identification ........................... 99
4.2 y-ray Sorting ............................... 100
4.2.1 Phoswich, Time Of Flight, and 7 Time Cuts .......... 101

4.2.2 Principles of Doppler Correction ................. 102

4.2.3 Angular Distribution Of '7' Rays ................. 110

4.2.4 Vacuum Deorientation ...................... 112

4.3 Transient Field Effect ........................... 115
4.3.1 Experimental Spin Precession .................. 115

4.3.2 Evaluation Of Transient Field .................. 116

4.3.3 2+ 9 factors of 38’408 ....................... 120

4.4 Interpretation Of Measured g(2+) values ................ 121
4.5 Summary and Outlook .......................... 127
Bibliography 131

vii

List Of Figures

1.1 Region Of the nuclear chart near 38408 ................. 3
1.2 Shell model orbitals ............................ 4
1.3 First-excited 2+ states in sulfur isotopes ................ 5
1.4 Proton and neutron magnetic dipole moments ............. 13
1.5 Schmidt limits ............................... 17
2.1 Experimental principle Of transient field technique ........... 24
2.2 Nuclear orientation examples ....................... 26
2.3 Stopping powers for sulfur in Fe and Cd ................ 36
2.4 Conventional transient field target for inverse kinematics ........ 41
2.5 Projectile and target in intermediate-energy Coulomb excitation . . . 48
2.6 HVTF parametrization in Fe ....................... 50
2.7 Measured TF strengths in Gd hosts ................... 52
2.8 Slowing-down time at different energy ranges .............. 56
3.1 CCF diagram ............................... 59
3.2 PIN spectrum from target test. run ................... 63
3.3 SeGA segment naming convention. ................... 64
3.4 A SeGA detector .............................. 65
3.5 Coordinate system forGe detectors ................... 66
3.6 SeGA configurations for angular distributions .............. 67
3.7 Full SeGA configuration for TF 9 factors ................. 68

viii

3.8

3.9

3.10
3.11
3.12
3.13
3.14
3.15
3.16
3.17
3.18
3.19
3.20
3.21
3.22

4.1
4.2
4.3
4.4
4.5
4.6
4.7
4.8
4.9
4.10
4.11
4.12

Phoswich light components ........................ 72

Phoswich positioning geometry. ..................... 73
Phoswich detector dimensions ...................... 75
Specular reflection in the light guide. .................. 77
Oscilloscope trace Of the phoswich Signals ................ 78
Phoswich electronics diagram ...................... 79
Phoswich response tO radioactive sources ................. 80
Gain shift Of phoswich with rate ..................... 82
POISSON calculation Of magnetic flux lines ............... 86
Plot Of magnetic field vs. current ...................... 88
Plot Of magnetic field along beam axis. ................. 88
SeGA electronics diagram ........................ 91
Trigger logic diagram ........................... 93
Layout in N3 vault ............................. 95
Cross section Of target. chamber ..................... 97
Phoswich Tail vs. Full Energy ...................... 100
Phoswich TOF vs. E ........................... 101
SeGA TDC spectra ............................ 103
Properties Of Compton scattering .................... 105
Components Of SeGA energy resolution ................. 107
Determination of average after-target velocities ............ 108
Comparison Of different Doppler corrections .............. 109
Angular distribution Of y rays ...... 3 ................ 111
Fits to projectile charge states after target ............... 114
Sign convention for detector angles ................... 117
Monte Carlo y-ray Spectrum and spin precession ............ 120
Known sulfur isotope 9 factors ...................... 121

4.13 Shell model and experimental g(2+) values ............... 123

4.14 Proton partitions Of Shell-model wavefunctions ............. 126
4.15 CCF Yields ................................ 128
4.16 Lifetimes Of selected 2+ states ...................... 129

List of Tables

3.1
3.2

4.1
4.2
4.3
4.4
4.5

Production and properties Of the radioactive beams ........... 60
Iron foil magnetization .......................... 85
Precessions and Slopes for each angle pair ................ 117
Nuclear parameters and reaction kinematics .............. 118
Shell model and experimental values compared ............. 122
Theoretical 9 factors compared with experiment ............. 124
Detailed decomposition of g factors .................... 124

xi

Chapter 1

Introduction

Many systems in nature that consist of charged particles possess electromagnetic mO—
ments, including atomic nuclei. The electromagnetic moments of a nucleus are among
its most fundamental properties, and Since the electromagnetic interaction is well
understood, nuclear theories can readily provide predictions of electromagnetic mo-
ments, which are experimental observables. Measurements Of certain electromagnetic
moments stringently test nuclear models because moments are very sensitive tO the
details of the nuclear wavefunction. In part, this is because the static electromagnetic
moments are expectation values of “one-body” Operators acting on single nuclear
states. Since electromagnetic moments probe the properties Of nuclear states, they
are particularly useful in transitional regions of the nuclear chart, where they can be
used to quantify the changes occurring and uncover their underlying causes.

The heavy sulfur isotopes have transitional characteristics, changing from spherical
to deformed as neutrons are added beyond the stable, naturally occurring sulfur
isotopes. The experimental evidence for the changing structure will be summarized in
the following sections, where it will be pointed out that although several observables
have been measured in this region Of the nuclear chart, an understanding Of the
underlying causes Of the shape transition is lacking. The onset Of deformation in

the sulfur isotopes is interesting because it coincides with the crossing of neutron

number N = 20, and continues toward N = 28, which are two “magic” numbers Of
neutrons in the stable nuclei, associated with extra binding energy indicative Of a large
energy gap in the nuclear shell structure. With increasing neutron excess, different
orbits begin to be occupied compared to the Situation at stability, and changes in
the relative spacing of those single-particle orbits occur [1]. Since exotic nuclei are
difficult to produce, experimentally measured properties are Often sparse, providing
no constraint for theoretical predictions, which can vary widely. It is the goal of this
dissertation to describe new measurements on the neutron-rich sulfur isotopes that
give a defining perspective on the evolution Of their shell structure. After giving some
details of the physics Of the neutron-rich sulfur isotopes, an introduction to the nuclear
excited-state magnetic dipole moment will be given, along with the justification Of
how the magnetic moment enables the elucidation Of the underlying contributions to

the structure Of the neutron-rich sulfur isotopes.

1.1 Physics of Neutron-rich Sulfur Isotopes

The isotopes studied in the present work, 38S and 40S, are located in the chart Of the
nuclides in between the N = 20 and N = 28 closed neutron shells, where neutrons
occupy the f7 /2 shell-model orbital. A view of the surrounding nuclides can be seen
in Fig. 1.1, with the stable isotopes Shown in black. Sulfur isotopes have 16 protons,
and the last two protons fill the 231/2 orbit, forming a Z = 16 subshell closure which
is partially evidenced by the fact that sulfur has the largest number of stable isotopes
for any element lighter than calcium. A schematic Of the relevant shell model orbitals
is shown in Fig. 1.2. The heaviest stable isotope Of sulfur is 333820, which has (in
addition to the subshell closure at Z = 16) a full neutron d3/2 orbital which closes
the N = 20 major shell. That situation makes 363 a stable, almost doubly-magic
nucleus, with a high-energy first-excited state Observed at 3.29 MeV. Adding one

neutron pair to 36S makes 38S, a radioactive nuclide with T1/2 = 170.3 m and a

  
  
  

    

Si

26 28

   

22 24

  

20

Figure 1.1: Region Of the nuclear chart near 38*‘lOS, with the stable isotopes in black.
The N = 20,28 closed neutron Shells are shown outlined in black. Sulfur isotopes
have 16 protons.

2.9 MeV fl—decay Q value. Adding one more pair of neutrons makes 40S, the mid-shell
nucleus for the N = 20 — 28 shell, with T1/2 = 8.8 s and a 4.7 MeV fi-decay Q value.
Initial interest in the neutron-rich sulfur isotopes arose when anomalous 46Ca/48Ca
abundance ratios were Observed in meteoritic inclusions [2], and it was suggested that
the neutron-rich sulfur isotopes play an important role in the nucleosynthesis Of heavy
Ca—Ti-Cr isotopes (see Ref. [3] and references therein). Soon thereafter, indications of
transitional properties arose as measurements began to be made on the neutron-rich
sulfur isotopes, and it was recognized that these isotopes also offer a probe of changing
shell structure as the neutron number approaches 28, away from the stability line. A
survey of the experimental study Of these isotopes prior to the current work will be

given in the following section.

1.1.1 Previous Experimental and Theoretical Work

The isotope 388 is accessible with two—neutron transfer onto 368, which is most readily
done using (t,p) reactions. Since the isotopic abundance Of 368 is only 0.02% Of

natural sulfur, the (t, p) measurements were difficult to perform but resulted in the

, 1A,.
2p «77"; 2P In

-‘
4 I ‘~----

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

If 2p
,\‘\\ 28 3/2
20 ‘ 20 lfuz
Id ---- ld3/2
'''''''''''''' 23
2s 16);”:
8 8
1P —:::: ------ 1P l/2
..... 1pm
2 2
ls -- Isl,2
Harmonic HO with Spin-Orbit

Oscillator Interaction

Figure 1.2: Single-particle shell model orbitals relevant to the neutron-rich sulfur
isotopes. The left hand side shows the energy levels given by a harmonic oscillator
(HO) potential, and the right hand Side includes the effect of a negative spin-orbit
force, that Splits each harmonic oscillator level in two, conventionally labeled nfj.
Each level can be occupied by 2 j + 1 particles of each nucleon type. The numbers in
between large gaps in the single-particle levels are the magic numbers Of nucleons.

first studies of 38S in the 1980’s [4, 5, 6]. Other transfer reactions, such as two—proton
pickup from 40Ar in the reaction 40Ar(llB,13N)383 [7], and two-neutron transfer with
heavy ions using 36S(14C,12C)38S and 36S(180,1607)38S reactions [8] also Obtained
low-energy level scheme information. The heavier even-even sulfur isotopes could not
be produced in sufficient quantities to study using nucleon transfer reactions, and
since fusion-evaporation reactions are generally more amenable to neutron-deficient
systems, further progress in studying these isotopes was difficult. It is notable that
Fornal and coworkers were able to glean some '7 rays belonging to 38S from a re-
analysis [9] Of 77 coincidence data from binary reactions Of 36S+160Gd, by gating on
158Gd particles, and this method may hold some promise in reaching the high-energy
levels in some neutron-rich nuclides. However, more direct means of accessing nuclei
in order tO use other experimental probes needed to await prOgrcss in radioactive

beam development.

 

 

 

 

 

 

 

 

w...
E
.—. 3000» ~23 375 .
E {I
F 2000* g 250 .
N +
v o
”J 1000» I a 125 »
Eu“
0 I I I I I I I m 0 I a RI I I I
14 I6 18 20 22 24 26 28 30 I4 16 I8 20 22 24 26 28 3O
Neutron Number Neutron Number

Figure 1.3: Previously measured 2+ excited state energies (left) and B(E2; 0+ —> 2+)
values (right) for sulfur isotopes, as a function Of neutron number. Data are taken
from Refs. [10, 11] and [12]. Lines are drawn to guide the eye.

Fragmentation reactions can provide high yields Of neutron-rich nuclides, and
this technique made it possible for neutron-rich sulfur isotopes to be studied with
intermediate-energy Coulomb excitation at the N SCL [13], yielding E(2+) and B(E2)
values for the first-excited 2+ states in 38’40’42S [10]. The Coulex measurements put
on solid experimental ground the idea that the neutron-rich sulfur isotopes lie in a
region of rapidly-changing nuclear structure [14, 15, 16]. The N = 28 nucleus 44S was
also studied with intermediate-energy Coulomb excitation, which provided a measure-
ment Of E(2+) and the B(E2) to the first 2+ state [11]. The results are summarized
in Fig. 1.3. From the 2f energies, it appears that the nuclei become more collective
as they approach midshell. This is confirmed by the increasing B(E2) values, which
indicate that the heavy sulfur isotopes are becoming deformed as more neutrons are
added on top Of the spherical 368. Of particular interest is the behavior Of the sulfur
isotopes at the purported closed shell at N = 28, where some collectivity has been
found to remain in 44S [11]. Since the B(E2) value is proportional to the square Of
the proton deformation parameter I32, the Coulomb excitation measurement was only
able to establish that 44S exhibits collectivity, and could not measure the sign Of I32
or tell whether the deformation is static or dynamic in nature.

Along with the Coulomb excitation studies. fl-decay studies Of ~10P were carried

out at the N SCL tO study the level structure Of 4OS [17], and proton inelastic scattering
using a CH2 target Offered another probe Of collectivity and deformation of the 2:5
state in 38S [18].

In-beam 'y-ray Spectroscopy Of neutron-rich sulfur isotopes was carried out fol-
lowing fragmentation reactions [19] at the GANIL facility in France. The workers in
Ref. [19] were able to propose level schemes based on the y—ray spectra, including the
locations Of 4+ states, which in the case Of 403 contradicted the previous study by
fi—decay [17]. Nevertheless, the results Of Shell—model calculations and collective model
calculations that explained the proposed level schemes indicated that the N = 28 shell
closure weakens below Z = 20, and is probably eroded for Z = 16 since the proton
subshell closure that is strong at N = 20 has evidently disappeared at N = 28, ac-
cording tO systematics Of single-particle energies in the K isotopes [20]. Even though
the measurement of higher-lying excited state energies can give an indication Of the
behavior Of a nucleus within the framework of collective models, and can test shell
model calculations, the results Of Ref. [19] can be viewed as somewhat tentative until
further confirmation Of the level schemes is obtained. In general, the experimental
data gathered SO far leave an incomplete picture Of the underlying causes for the
transitioning nuclear properties Of 36443.

The neutron—rich sulfur isotopes have been the subject of several theoretical in-
vestigations that have paralleled the experimental studies. Werner et al. [15] studied
the heavy sulfur isotopes within the Relativistic Mean Field (RMF) framework and
found that the potential energy surface Of the neutron-rich sulfur isotopes is some-
what flat with respect to deformation, in general leading to nearly degenerate prolate
(I32 > 0) and oblate ([32 < 0) Shapes. Other calculations by the same workers us-
ing non-relativistic Skyrme HartreeFock (SkHF) techniques found nearly degenerate
prolate-deformed and oblate-deformed solutions as well. The onset of deformation
in the neutron-rich sulfur isotopes was confirmed by the Coulomb excitation studies

Of Refs. [10, 11] which found [I32] ~ 0.3. Other theoretical studies followed Werner

and collaborators, where RMF studies using improved interactions also confirmed the
theoretical presence Of deformation [21, 22]. Shell model calculations were also under-
taken by Retamosa and coworkers [20], and were later refined by N ummela et al. [23]
and extended by Caurier, Nowacki, and Poves [24] to explain the B(E2) values and
2+ excitation energies Of the neutron-rich sulfur isotopes. The Shell model studies also
found the presence of deformation in the sulfur isotopes, and postulated strong shape
coexistence at the N = 28 nucleus 44S. Attention was beginning to shift to a larger
nuclear physics question: the validity Of the N = 28 magic number for neutron-rich
nuclides with Z < 20. The RMF. SkHF, and Shell model investigations, together with
other studies including angular-momentum projected Generator Coordinate Method
(AMPGCM) calculations of Rodriguez-Guzman et al. [25], all suggested that the
N = 28 shell closure weakens below calcium.

The experimental as well as theoretical effort devoted to this region of the nuclear
chart has been and continues to be significant. The conclusions about the nuclear
structure that have been drawn are varied, largely because the physical interpretation
of level energies and reduced transition probabilities are model-dependent, and many
of the models used in the theoretical studies are insufficiently constrained by the
available data. The lack Of a consistent interpretation Of the available data suggests
that g-factor measurements, which stringently test nuclear models as will be shown
in the next section, are crucial to aid in the understanding Of the interplay between
the migration of single-particle energy levels and the development. Of collectivity in

the neutron-rich sulfur isotopes.

1. 1.2 Open Questions

The changes in nuclear properties that have been Observed for the heavy sulfur iso-
topes have left questions about their causes. In particular, the migration Of the single-
particle energy levels and their occupancies in this region are not well understood.

Since the N = 20 energy gap has been crossed for the sulfur isotopes heavier than

368, it would be expected [5] that the addition of neutron pairs to form each succes-

sive even-even sulfur isotope would dominate the single-particle nature of the nuclear
wavefunction, especially for 38S which lies close to the N = 20 closed shell. On the
other hand, the increase in collectivity that has been Observed would tend to lower
the role Of single-particle effects as protons and neutrons become strongly coupled
together to form excited states. There is a delicate interplay between collectivity and
Single-particle effects in nuclear excitations, and care is needed when drawing con-
clusions from excited-state energies and B(E2) values. In fact, Since the B(E2) is a
reduced matrix element between an initial and a final state, some ambiguities can
arise in the determination of the nuclear properties. The 9 factor is not subject to
such uncertainties Since it is an observable Of the properties Of a single state. For
this reason the onset Of collectivity can be probed using g(2+) values to gain further
insight into the single—particle contributions to the excited-state wavefunctions.

The large energy gap after neutron or proton number 20 is due to the harmonic
oscillator closed shell consisting Of the d5 /2'31 )2-(13 /2 orbits (called the 3d shell, see
Fig. 1.2). The next. level is the f7/2 orbital, which has descended from the orbitals
composing the pf shell due to the spin-orbit interaction, leaving an energy gap above
it that causes the N = 28 magic number. The Shell closure at N = 28 can be called
a “spin-orbit closed shell” because it is primarily due to the spin-orbit force that
there is a gap between the f7/2 orbital and the pf Shell. It had been suggested by
Dobaczewski and collaborators [26] that the spin-orbit force weakens in nuclei with
a large neutron excess — this would result in the raising Of the f7/2 orbit and the
disappearance Of the N = 28 shell closure. Single-particle effects are a result of large
gaps in the single-particle energy levels, and for the Single-particle behavior Of the
neutrons, the N = 20 and N = 28 Shell gaps will clearly be the most important for the
38—448 isotopes. The proton single—particle behavior is most influenced by the proton
31/2-d3/2 gap. As neutrons begin to occupy the f7/2 orbit, therefore changing the

number of valence neutrons. proton-neutron interactions between the f7/2 neutrons

and the sd protons can affect the energies of proton single-particle energy levels, in
particular the proton 31/2—d3/2 gap [27]. It was suggested [28] that the mechanism
of proton-neutron interactions that change strength with increasing neutron valence
number is responsible for the changing structure Of the sulfur isotopes, by eroding the
Z = 16 subshell closure. In contradiction to the assertion in Ref. [28], a weakening Of
the N = 28 shell gap away from stability was recently proposed [19, 24] as the main
cause of deformation in the neutron-rich sulfur isotopes. It will be shown that the 9
factor can separate the proton and neutron contributions to the nuclear wavefunction,
which will enable the effects of the Single-particle energy shifts, and the possible
existence Of shell closures, to be deduced in the sulfur isotopes.

The existence Of shell closures at the traditional magic numbers toward the driplines
is an open question, since many changes in Shell gaps have already been Observed.
For example, the N = 20 shell closure has been found to soften for decreasing proton
number below Z = 16, with clear evidence for non-magic behavior at the Z = 12
nuclide 32Mg [29]. AS was indicated previously, much attention has been focused on
the N = 28 shell gap below calcium, since radioactive beam facilities have now pro-
vided experimental access to this region. It remains an open question as to whether
some nuclides along the N = 28 line can be called “magic” nuclei, most notably
42Si, which has been the subject of conflicting claims from fi-decay studies [30] and
from nucleon knockout measurements [31]. Cottle and Kemper had suggested [32] that
spectroscopy Of Z = 14 42Si would give a definitive answer about the possibility of the
N = 28 magic number below Z = 20, rather than measurements in the neutron-rich
sulfur isotopes; their reasoning is that the effect of the disappearing Z = 16 subshell
closure is strong, but is limited to the sulfur chain. Since an excited-state g-factor
measurement probes the nuclear wavefunction in a way that allows an elucidation Of
neutron and proton behavior (including orbital occupancies and shell gaps), g-factor
measurements Offer a means to search directly for the persistence Of the N = 28 shell

closure in the sulfur isotopes.

9

1.2 Electromagnetic Moments Of Nuclei

Electromagnetic multipole moments are categorized as either electric or magnetic
moments. Electric multipole moments arise from the distribution Of electric charge
within the nucleus, and magnetic moments arise from the distribution Of currents and
magnetization. Calling the charge density within the nucleus p(x’), the potential <I>(x)
outside the region of the nucleus can be expanded in a series of Spherical harmonics.

Following Jackson [33],

1
mi 21W (1 1)
(m_0 __ 54—1” + 1E1," Tl+1 ’ .

where the values Elm are the coefficients Of the expansion, to be calculated from p(x’).
The Elm coefficients are found by equating Eqn. 1.1 with the explicit expression for

the electrostatic potential,

<I>(x)=/ ”(x’) d———-3:r, ’ (1.2)

Ix- x'l

for which the term 1/ [x — x’ I can be expanded using

00 l
[X— x’] 471.2% "112: 214.17. 1+1 lYlmw, ¢,)Ylm(99 (15,) (1.3)

yielding
_ 1 a: I I II I 3,! Y1m(9,¢)) ,
@(X)—47TIZ:§:;—l [fn172(6,¢)7‘ p(X )d .L' T. (1.4)
,m

The quantities in brackets are the Ejm’s, known as the electric multipole moments,

Elm = fr’ln;(6,¢))p(x’)d3r’. (1.5)

Since the nucleus is a quantum system, Eqn. 1.5 must be cast. into a quantum me-

chanical form using p(x) = 2,221 ell/2A2. In a classical picture, magnetic multipole

10

moments can be defined in a manner analogous to Eqn. 1.5, as

in... = / r"n:,,(6,¢)p.u(x’)d3x’, (1.6)

where the magnetic density is pM(r) = —V - M(r), with M(r) = Mc(r) + M3(r‘) is
the sum of contributions from electric currents and intrinsic Spin Of the nucleons [34].
The current density comes from the motion of the protons, leading to contributions
of the form r x p, which is their orbital angular momentum. The spin contributions
arise from contributions of the neutrons as well as the protons, since both particle
types carry intrinsic spin.

For symmetry reasons, the nucleus can not possess Odd-rank static electric mo-
ments or even-rank magnetic moments, and indeed those moments (electric dipole,
magnetic quadrupole, etc.) have been found experimentally to be vanishingly small.
The electric monopole (E0) moment is simply the nuclear electric charge. The only
higher electromagnetic moments experimentally established in nuclei are the mag-
netic dipole (M 1), electric quadrupole (E2), and magnetic octupole (M3), with some
evidence pointing to the existence Of the electric hexadecapole (E4) moment, although
the effects of the higher-order moments rapidly diminish with increasing rank. The
magnetic dipole and electric quadrupole moments are the most important for describ-
ing nuclear structure. The electric quadrupole moment represents the deviation of the
nuclear charge from a Spherical shape. The existence of the magnetic dipole moment.
means that electric currents and magnetization are present within the nucleus, which
indicates that the nucleus possesses internal structure. The classical definition of the
magnetic dipole moment is [1. = f M (xI )d3x' . Quantum mechanically, the contribu-
tions to M(x) from the current and spin terms can be written in terms of the angular
momentum and Spin operators, £7 and s, which will be explained further in the context
of the nuclear shell model.

When a magnetic dipole moment is immersed in a magnetic field, a torque (T =

11

u x b) is exerted on the magnetic moment that causes a rotation with frequency

w = gu7NB, (1.7)

which is the Larmor precession formula for a nucleus in a magnetic field Of strength
B. The usual unit for nuclear magnetic moments is the nuclear magneton, [,IN =
err/2mp 2:: 5051-10"27 J /T (mp is the proton mass). The quantity 9 is the gyromag-
netic mtio, defined as the proportionality constant between the value of the magnetic

moment and the nuclear spin,

[1. = g/INI. (1.8)

It will be Shown that the numerical value of 9 contains information on the nuclear
wavefunction. A measurement Of the Larmor frequency in Eqn. 1.7 is the principle
on which most techniques for measuring magnetic moments are based. The frequency
depends on the product of the field strength and the 9 factor, as can be seen from
Eqn. 1.7. Extraction of the 9 factor therefore relies on the knowledge or measurement
Of the magnetic field strength.

The nucleons (protons and neutrons) are the constituents Of nuclei, and they have
their own complex internal structure. This structure gives the proton and neutron
magnetic dipole moments, which possibly arise from the distribution of quark cur-
rents and magnetization within the nucleons, although at present these are not fully
understood. The nucleon magnetic moments have been measured for the free proton

and neutron, with the resulting nominal values [35]

+2.793 for rotons
#N I) ( 1 . 9)

‘t:
H

-—1.913 IIN for neutrons,

known at present to a precision of A/t/ [1 ~ 10‘7 for the neutron and 10’8 for the

proton. It is the large difference (in both Sign and magnitude. see Fig. 1.4) of the

12

spin [1 4

 

neutron proton

Figure 1.4: Proton and neutron magnetic dipole moment vectors (thick grey arrows)
and the proton and neutron spins (thin black arrows).

neutron and proton magnetic moments that makes the nuclear magnetic moment
such a sensitive probe of the internal arrangement Of neutrons and protons within the
nucleus, since the magnetic moments Of the individual nucleons combine to produce
the total nuclear magnetic moment. A measurement Of the nuclear magnetic moment
therefore gives insight into the particular way neutrons and protons are behaving in a
nuclear state Of interest. It is useful tO interpret the magnetic moment of a nucleus in
terms of a particular model in order to gain insight into the structure of the nucleus.
Two broad classes of models, which generally predict very different 9 factors, are
used to interpret magnetic moments. Microscopic-type models, like the shell model,
focus on the motions Of the individual constituent nucleons. In contrast, macroscopic
collective models focus on the bulk properties of the whole nucleus. Brief introductions

to the different approaches will be given in the following two sections.

1.2.1 Magnetic Dipole Moments in the Shell Model

One Of the most interesting aspects Of the nuclear system is its many-body nature.
In contrast to atomic systems, the nucleons interact very strongly with each other,
and this has interesting consequences. In close analogy to electrons, the nucleons are
spin 1/2 particles which must Obey the Pauli Exclusion Principle. The consequence

of this principle and the nucleon-nucleon interactions within the nucleus are nuclear

l3

properties that Show a periodic Shell structure reminiscent of the atomic structure
seen in the Periodic Table of the Elements. A successful nuclear shell model exists
that can explain many Of these properties for the stable nuclides.

Shell model wavefunctions are Obtained by assuming initially that the nucleons are
moving inside an average potential well created by all of the nucleons, called the “mean
field.” Modern shell model calculations also take into account the residual nucleon-
nucleon interactions. In the extreme single particle Shell model, the Single—particle
wavefunctions can be obtained by solving the SchrOdinger equation with the mean
field potential, and energy levels for particles with quantum numbers (n, f, j) result.
The quantum number n represents the number of nodes in the radial wavefunction, t is
the orbital angular momentum, and j is the total angular momentum that results from
coupling the intrinsic nucleon spin 5 = % to the orbital angular momentum to Obtain
the two possibilities j = E :t i. The nuclear potential generally has a Woods-Saxon
shape, and during the development Of the shell model, a spin-orbit term (—f . s) was
added [36, 37] in order to account for the experimentally-determined “magic numbers”
which correspond to large gaps in the single-particle energy levels. Since nucleons fill
the lowest energy levels first and Obey the Pauli principle, nucleons are added to a
particular (17., 8, j) level until the maximum Of 2 j + 1 protons or neutrons is reached.
Additional nucleons must then begin to occupy higher-lying states. This means that
in some cases, certain nuclides with A nucleons can exhibit an extra degree Of binding
relative to that of the nuclide with A + 1 nucleons. Empirically, these magic nucleon
numbers are found in stable nuclei to be 2, 8, 20. 50, 82, and 126, for either the
proton1 number Z or the neutron number N. Certain nuclei that have magic numbers
of both their protons and neutrons are termed doubly magic (such as 4He and 16O).
which are extremely tightly-bound systems. It is not true that the magic numbers

are universal for all possible combinations Of Z and A, because they depend on the

 

lNo nuclide has been found to exist with as many protons as 126. Such an isotope could lie
beyond the limits Of nuclear existence. or perhaps is simply uI'Ircat-liable with current experimental
techniques.

14

ordering and relative spacing Of the energy levels caused by the nuclear potential well
and nucleon-nucleon interactions, which can change with nucleon number. It is one
of the main goals of nuclear physics to observe and attempt. to explain this change,
in particular for nuclides far removed from stability.

A magnetic dipole moment arises from contributions from the orbital and intrinsic
spins of the protons and neutrons that. make up a particular nucleus, in the same
way that the nuclear spin does (the two quantities are closely connected, as seen in

Eqn. 1.8). The quantum-mechanical magnetic moment operator is given by

u = 141v 2 [art + 988], (1-10)

where the sum runs over all nucleons in the nucleus, for which the orbital 9 factor
g; equals 1 for protons and 0 for neutrons. t is the orbital angular momentum and
s is the intrinsic spin. The spin 9 factors, gs. are identical to the nucleon 9 factors,

Obtained from Eqn. 1.9 by dividing by the nucleon spin. 1/ 2:

+5586 for protons
gs = (1-11)
—3.826 for neutrons.
The magnetic dipole moment is a sensitive probe Of the single particle nature Of
nuclear states, since the Operator involves terms that depend on the coupling Of t and
s, and sums over all the nucleons composing the nucleus. The fact that the proton and
neutron 9 factors are different in both Sign and magnitude means that the nuclear
g factor can change significantly with t’ and j, as well as with particle type, which
means that shell model calculations using different interactions or model spaces can
be stringently tested.
Using the operator in Eqn. 1.10, the magnetic moment. can be Obtained from a

nuclear wavefunction by calculating its expectation value for the wavefunction of the

15

m = I substate,

(Im=l|/I|[m=1). (1.12)

A single-particle state with a particular j gives rise tO a particular magnetic moment
due to the operators t’ and 3 contained in Eqn. 1.10. Since this is the sum of two
angular momenta, one must obey angular momentum algebra in determining the
result. Following Brink and Satchler [38] the addition theorem for the 9 factor of a
state I that is the result Of summing two angular momenta 11 and 12, with 9 factors

91 and 92. IS

“([1 + l) - 12(12 +1)
I(l+1) '

 

1 1
9 = -({11+g2)+ 5(91 -.(12) (1-13)

2

The addition theorem can be applied using gl = 9;, 11 = t”, 92 = gs. 12 = s, and
I = j, to Obtain the g factor Of a. particle in an orbital with angular momentum
j = I :1: % in a compact. form [39]:

gs _ {If

.-: ,i , . 1.14
g} g( 2€+1 ( )

 

Using Eqn. 1.14, simplified predictions can be made for the ground state magnetic mo-
ments of Odd-mass nuclei within the “extreme Single-particle model” which assumes
that all protons and all neutrons form spink0 pairs, and the remaining unpaired nu-
cleon is responsible for the spin and magnetic moment Of the Odd-mass nucleus. The
predicted magnetic moments for Odd-mass nuclei are called the Schmidt limits be-
cause the experimentally measured values generally lie in between the j = [+- 1/2
and j = 5—1/2 values given by Eqn. 1.1-I (see Fig. 1.5). In reality, the magnetic mo-
ments usually deviate from these extreme predictions by a factor Of about 0.6. The
reason for the deviation is that. the model is not realistic enough. Usually, effective g
factors are used to provide a correction for missing interactions in the Hamiltonian,

and 9:“ a: 0.6g?“e is Often used. and there can be corrections for g,r as well. Modern

16

 

 

 

F— 6 l.5 7 o .jzl'l/z"

z . o

a 5 LO 3 g

g 4 0.5 * | 8 0

g 0.0 ~ 2 o g 2

E 3 -0.5 . l I I °

-§ 2 -l.0 » § l l
° D

50 I 'l.5 " : o O

a -20 ~ 3 -

2 0 ' ° j=i+I/2

 

 

 

 

 

 

 

'/2 3’2 5/2 7/2 9/2

Nuclear Spin Nuclear Spin

Figure 1.5: Measured magnetic dipole moments in Odd-proton nuclides (left), Odd-
neutron nuclides (right), and Schmidt predictions (lines). Data are from Ref. [43].

shell model calculations such as those Of Brown [40] or Caurier and coworkers [41]
have prescriptions for including such corrections, and more details can be found in
the literature [42, 39].

The Schmidt prediction in Eqn. 1.14 shows that the 9 factor can change Sign
depending on whether the particle is a proton or a neutron (because Of their respective
spin 9 factors). and depending on whether the orbit is one with j = (+1 / 2 or j = 8—1 / 2
(because of the i Sign in Eqn. 1.14). A measurement of the sign of the g factor is very
useful; this information can Often establish the type Of particle that is dominating the
nuclear structure. The magnitude Of the g’factor can be used to sort out the relative

contributions Of the protons and neutrons to the wavefunction of a nuclear state.

1.2.2 g(2+) in Even-even Nuclei

Even-even nuclei always have zero ground state spin and therefore no ground state
magnetic dipole moment. Excited states in even-even nuclei, however, can have finite
spin and therefore are expected to have magnetic moments. The situation can be
more complicated than in Odd mass nuclei, that can be thought Of as an inert core
with an unpaired nucleon that is responsible for the spin and magnetic moment. In

particular. the first excited state of most even-even nuclei has spin 2. The Spin-2

17

state can be assembled from the single particles in many ways, such as promoting a
pair Of nucleons to a higher-energy orbit, breaking a pair and promoting one nucleon
to another level, or by simply recoupling a pair of nucleons within the same j level
to angular momentum J = 2. In all cases, the addition theorem can be used to
obtain predictions for the g factor Of such configurations, with especially Simple results
for identical particles in the same j orbital, where g(J) = g( j). Therefore, within
the context of the shell model, a 9 factor can be useful to determine the nucleon
configuration for both the protons and the neutrons in an excited state. It is rarely
only a single particle (or a single pair of particles) that make up the excited state
wavefunction; for certain nuclei — most notably those with very large quadrupole
deformation - the participation of a very large number Of nucleons is needed to explain
the large transition probabilities (B(E2) values) that can be 10—100 times larger than
what would be expected from a few Single-particle excitations. In the following, the
example of a hydrodynamical picture is used to give a crude estimation Of a collective

g factor.

Hydrodynamical Model

In contrast to the shell model, which treats the nucleus on a microscopic basis by con-
sidering the interactions among the constituent nucleons, the nucleus can be treated
as nuclear matter, similar to a liquid drop or a rigid body. Such systems, in which all
the nucleons act together, are termed macroscopic or collective. Macroscopic models
such as the liquid drop model [44] and the collective Bohr model [45] can describe
many features of deformed nuclei, which tend to Show characteristic patterns of rota-
tion Or vibration. Rather than describing the details Of the various collective models,
it is more useful in this context to discuss the excited-state 9 factors that result from
these models. In particular, a hydrodynamical treatment of the nucleus as a rotating
charged fluid predicts 9 ~ +Z/A, and this value serves as the canonical “collective”

prediction for excited-state 9 factors although it is not clear that. a real nucleus (even

18

one that is a good rotor) should behave in this fashion. A brief justification for the
Z /A result will be given here.

The classical gyromagnetic ratio of an object is the ratio Of its magnetic moment
to its angular momentum. Using elementary electrodynamics, the expected gyromag-
netic ratio can be calculated for a simplified version of the nucleus: a uniformly-
charged rotating sphere. Let the Sphere have mass Amp and charge Z c distributed
uniformly throughout its volume grR3, where R is the radius Of the sphere. Then,
for a rotational frequency w, the sphere will have angular momentum Iw, with the
moment of inertia I = gMRZ. Therefore, the denominator Of the gyromagnetic ratio
is equal to gAmpR 2w. The magnetic moment of the charged rotating sphere is calcu-
lated by recalling the result for the magnetic moment Of a current loop, #100,, = Ia,
where I is the current and a is the area enclosed by the loop. Using cylindrical coordi-
nates with the z axis along the axis Of rotation, a succession of current loops can form
disks lying in planes Of constant z, with radii corresponding to the cylindrical radius
Of the sphere, r(z) = m. These disks will have magnetic moments dam-“(2),
which can add together to form a sphere with a resulting magnetic moment ”sphere-
The electric charge in a loop Of radius r is 21rrA, where A = pdz dr is the charge per

unit length and p is the volume charge density. SO the current can be written as

_ 27rrA

d] _ 27T/w

 

= rwp dz dr,
using the time for one revolution Of 27r/u), and for a loop Of radius r one Obtains
dmooph) 2 la 2 mopr3dz dr.

The magnetic moment Of the sphere is Obtained by integrating the current loops out

to the cylindrical radius Of the sphere, and then integrating over the z coordinate

19

from 2 = -R to z = +R. This can be expressed as

+12 m
Iisphere = 7rwp / dz / dr 1'3
—R 0
4

= gimp]?5

3 2

(1.15)

By taking the ratio of Eqn 1.15 and the angular momentum one can obtain the

classical gyromagnetic ratio Of a charged rotating sphere:

I‘.sphere_ ézeWRrZ _§_. 8
Jsphere IEATHpWI—t)2 A 2m!)

 

(1.16)

gelass :

Thus the classical result for a charged rotating sphere shows a Z /A dependence,
expressed in units of e/2mp, the classical analog of the nuclear magneton.

The collective models based on classical rotators or surface vibrations tend to pre-
dict 9 factors near +Z/A. For rotational nuclei, this is indicative of how much of the
angular momentum is carried by the protons with respect to the entire nucleus [46].
Vibrational nuclei also have a 9 factor of Z /A that is dominated by the orbital mo-
tion Of the particles [45]. It is very significant that in collective models, the spin
contributions of the nucleons are ignored. In the rare-earth mass region, calculations
performed by Prior and coworkers [47] showed that the contributions from the Spins
of the nucleons are small and generally cancel. It is not clear if the Spin contributions
are small universally across all mass regions and N / Z ratios, since few excited-state

9 factors have been measured Off the line Of stability.

1.3 Organization of Dissertation

The merits of g-factor measurements for elucidating the structure of the nuclear state

wavefunctions have been presented in this Chapter, along with the background Of the

20

sulfur isotopes that are the subject Of this dissertation. Since no existing technique
was clearly suited to measuring g(2+) values of picosecond-lifetime excited states in
nuclei produced as fast fragments (which at present is the most accessible means Of
Obtaining the neutron-rich sulfur nuclei), a new technique needed to be developed.
This thesis concerns the development of a suitable technique and its application to the
transitional isotopes 388 and 408. The method is called the High Velocity Transient
Field (HVTF) technique and is presented in Chapter 2. The experimental realization
of the technique at the National Superconducting Cyclotron Laboratory is described
in Chapter 3. Results of HVTF g-factor measurements of the first-excited states of
neutron-rich 38’403 are reported in Chapter 4, as well as the outlook for further

radioactive fast—beam magnetic moment measurements.

21

Chapter 2

Transient Field Technique

It was recognized by Brady 8; Deutsch [48] early in the development Of the theory
of the angular distribution of radiation that a nucleus may be perturbed by external
fields prior to the emission Of radiation. In the case of a magnetic interaction between
the nuclear magnetic moment and a magnetic field existing at the Site of the nucleus,
the perturbation takes the form Of a Larmor precession Of the nuclear Spin about the
magnetic field axis, with the precession frequency given by Eqn. 1.7. Since nuclear
radiation is emitted in a characteristic angular pattern measured with respect to the
nuclear Spin axis, a precession Of the nuclear spin causes a change in the orientation
Of the angular distribution Of radiation by an amount proportional to the product
Of the field strength and the 9 factor Of the nuclear state, as seen by the following.

Rewriting Eqn. 1.7 using w = Ad/At yields
_ llN

where the negative Sign (due to the interaction H = —;I - B), indicates that the
spin precession proceeds in the direction of the torque r = [l x B. The average
magnetic interaction time is the excited-state lifetime, At = r(2+). Since the g factor

is generally Of order unity, the ratio [t N / n has approximately the value 48 MHz / T. and

22

the 2+ states Of interest have lifetimes in the picosecond range (1 ps = 10‘]2 s), the
magnetic field needs to be Of order kilOTeSla in order to induce a detectable precession
(of order milliradians). Such a field cannot be provided by laboratory magnets, and the
strongest static hyperfine fields present at the Site of nuclei embedded in ferromagnetic
materials are generally ~100 T. It was found, however, that the hyperfine magnetic
field experienced by ions swiftly traversing ferromagnetic materials is generally on
the order Of kilOTesla. Since this phenomenon occurs only while ions are in motion,
it was called the transient field.

The transient field was discovered in the late 1960‘s [49. 50] when it was found that
static-field induced precessions were affected if the ion under study first recoiled into
a ferromagnet and later came to rest, where it experienced a static hyperfine field. It
became apparent that ions precess at a different. rate while Slowing down in the ferro-
magnetic material, compared to the Spin precession expected from the known static
hyperfine field in the material. An intensive study of the properties of the transient
field initiated shortly after the Observation of the transient field phenomenon. The
investigation Of the transient field resulted in the development Of the transient-field
technique to measure short-lived excited-state g factors for states produced in nuclear
reactions, since reaction products Often have a significant amount Of kinetic energy
after the reaction and can recoil into a ferromagnetic medium by a suitable design of
the target. An early review [51] was made by Benczer—Koller, Hass, and Sak in 1980,
which contains an authoritative history Of the progress made in understanding the
behavior Of the transient field, as well as a description Of experimental techniques us-
ing the transient field. For a review encompassing many of the modern development-s,
see Ref. [52]. An outline Of the transient field technique, including a sophisticated
combination of Coulomb excitation in inverse kinematics with the transient field, will
be presented in this Chapter. An overview of a typical experimental arrangement is

shown in Fig. 2.1, which is taken from Ref. [53].

23

Y —detector

 

 

    
 

// ‘ \ / ‘\
I ‘(V |
l \ l
\ \ I
\ \ ’ I
\ ‘ ’ /
\ \ / I
\ \ ’ /
\ \\ / fl
\ /
\\ ’
I
\ \ / I
\ /
\ /
U /
beam recoil nucleus
n I ’ ‘ \\
port detector , / backing \ \
I I \
/ \
/ \
~/ V5 V1 \ \
\ \
/ / \ \
/ / \ \
/ / \ \
I ’ G ‘ \
/ \ \
I / -. \
I / / B \ I
\\ / / “\ II
\
\ _ .- / I ¢ \ \ ‘ v [I

Figure 2.]: Transient field experimental concept. from Ref. [53]. In this arrangement,
the beam enters through the hole in an annular particle detector and impinges on
the target layer Of the multilayer target. The beam can backseatter near 180° after
Coulomb exciting a target nucleus, causing the target nucleus to recoil through the
Fe layer, where it experiences the transient field. The dashed line is the 7-ray angular
distribution, with the direction Of precession shown by the arrows for a magnetic field
directed out of the page, and a positive 9 factor. Placement of a 7-ray detector at
the angle of maximum slope of the angular distribution (as shown) maximizes the
sensitivity of the measurement.

24

2.1 Measurement Principles: Stable Beams

Nearly all Of the techniques for measuring nuclear magnetic moments involve immers-
ing the nucleus in a magnetic field (which may be provided by a laboratory magnet
or via the internal fields in a molecular or solid-state environment), and measuring
the interaction between the magnetic moment of the nuclear state and the external
magnetic field, as expressed in Eqn. 1.7. A natural means Of accomplishing a mea-
surement Of the magnetic interaction (the Larmor precession) that takes place while
the nucleus is in an excited state is to measure a perturbed angular correlationl)
(PAC) Of de—excitation y rays from the nuclear state Of interest; this forms the basis
of the transient field technique. One method, Time-Differential Perturbed Angular
Distribution (TDPAD), measures the intensity variation of the 7 rays emitted into a
particular solid angle as a function of time, which directly results from the nuclear
spin axis “sweeping past” the 'y-ray detector at the Larmor precession frequency. In
practice, TDPAD typically requires a nuclear excited-state lifetime longer than a few
nanoseconds in order to allow at least one nuclear spin revolution to occur. If the
nuclear lifetime is less than a nanosecond, it is not practical to use time-differential
methods, and an integral method is usually used. Time-Integrated Perturbed Angu-
lar Distribution (IPAD), the category to which the transient field technique belongs,
detects not the frequency of precession but the angular degree of precession that a
nuclear Spin turns within a set time interval, Often its lifetime in a particular state.
All PAC / PAD techniques require nuclear spin alignment, in addition to knowledge Of

the field strength in order to extract the g factor.

 

1Another term, used almost synonymously with PAC, is perturbed angular distributions (PAD),
where the former usually refers to measuring the angle between successive '7 rays in a coincidence
measurement, and the latter usually refers to measuring the angular distribution Of "y rays with
respect to the beam axis or the direction Of a particle associated with the population of the state
that emits the r) ray. This dissertation uses the term “angular distribution” to mean the measurement
Of 7 rays with respect to the beam axis, which also corresponds to the associated particle direction.

Prolate Oblate
Alignment Alignment

\.
Tc:

Polarization

  

             

      

     

"Is

In: -2 -1 0 l 2 - -l 0 l 2 - -l 0 l 2
Figure 2.2: Examples of different types of nuclear orientation, shown as the relative

populations Of different magnetic substates for an I = 2 nucleus. Both types Of
alignment are symmetric about the m = 0 substate, while polarization is asymmetric.

I“. i4
‘s

2.1.1 Angular Distribution from Spin Alignment

TO Obtain information on the spin precession that is proportional tO the 9 factor,
excited-state spin alignment is necessary in order to allow a directionally-anisotropic
emission Of 7 rays from the state Of interest. Spin-aligned states (meaning an ensemble
Of nuclei with an unequal population Of magnetic substates, i.e. P(m) = P(—m)
but P(m) 74 P(n), where n is any other magnetic substate; see Fig. 2.2) can be
populated with compound nuclear reactions, direct transfer reactions, or Coulomb
excitation (for example), when there is a preferred axis in the reaction geometry. If
the details of the reaction are known, it is possible to calculate the exact distribution
Of magnetic substates and therefore predict the directional emission of 7 rays from the
excited states formed in the reaction. It is difficult to calculate the m-state population
distribution with much accuracy for most nuclear reactions because the interaction
is not completely understood. However, since the electromagnetic interaction is very
well understood, Coulomb excitation is a case where the m-state populations can be
very accurately calculated. This makes Coulomb excitation (presented in § 2.3.1) an
excellent means of preparing excited states which will decay by 7-ray emission with

a known directional distribution.

26

Fig. 2.1 illustrates the angular distribution for head-on Coulomb excitation Of the
2? state Of an even-even nucleus, where the beam is used to excite a target nucleus.
The experimental arrangement requires the detection of a beam particle scattering
backward at near 180° with respect to the beam axis, causing the excited target
nucleus to recoil forward into the ferromagnetic layer that provides the transient
field. In Such an arrangement, the angular distribution Of 7 rays has a characteristic
four-lobe shape due to the preferred population Of m = 0 substates in a backseattered
geometry. The degree Of spin precession can then be measured by detecting a change in
the 7-ray count rate at a known angle in the angular distribution, for which the slope
Of the angular distribution with respect to 0 is known. Generally, several detectors
are used in order to measure the 7-ray angular distribution and characterize the
slope. The angular distribution measurement. becomes especially important when the
excited state is populated by a compound-nucleus or transfer reaction and the m-
substate population distribution is not known. For Coulomb excitation, the substate
populations can be calculated from the well-known theory and a measurement at
several angles is important mainly to confirm the calculations and to ensure that the

system behaved as expected.

Formalism of Angular Distribution of 7 Rays

A 7 ray is typically emitted by a nucleus as it transitions from one state Of angular
momentum I,- and definite parity, to another state of angular momentum I f, which
also has definite parity. The 7 ray carries away the energy difference between the two
nuclear states, and it is also responsible for the possible changes in angular momentum
and parity. Due to the way that the photon electromagnetic fields couple to angular
momentum, the emission direction of the 7 ray has a relationship to the direction
of the initial nuclear state’s Spin axis [54]. This relationship is expressed analytically
as an angular probability distribution, which only has experimental meaning in the

context Of an ensemble of nuclei rather than only one nucleus. If the nuclei are oriented

27

randomly in space, then all Of their 7—ray angular probability distributions will overlap
to make an isotropic distribution. However, if the ensemble of nuclei is aligned along
a particular axis, then the 7 rays emitted from the ensemble may Show an anisotropy
in space with respect to that axis.

The theory Of the angular distribution Of radiation has been developed through
an examination of the effect the angular momentum change has on the direction of
the emitted radiation. Several references [55, 56, 57, 58] describe this formalism, and
a brief description Of the practical details will be given here. For the case Of axial
symmetry about the beam axis, it is most convenient to express the angular distri-
bution as an expansion in terms Of Legendre polynomials, and within the formalism
Of angular distributions, the coefficient for each term in the series is calculable. The

angular distribution is written as

kmaa:

W(6) = Z akPk(cosd), (2.2)
k=0

and it remains to determine the coefficients ak and the number Of terms in the series.
The emission of 7 radiation is caused by the electromagnetic interaction, which con-
serves parity. The conservation Of parity implies that the emitted radiation Should
not contain information about the direction Of the nuclear spin, but rather should
only have information about the nuclear spin axis. In other words, the angular distri-
bution should be a function Of the angle between the 7 ray and the orientation axis,
but Should not be different if 9 2 0° and t9 = 180° are interchanged. Functions with
that property must be symmetric about 90°. The Legendre polynomials Pk(cos 0) for
Odd k are not symmetric about 90°, but the even-ranked Pk(cos 0) are symmetric.
Therefore, the angular distribution is described only with even values2 Of k. It will be

shown shortly that a 7-ray angular distribution depopulating a 2+ state is completely

 

2The decay of a nucleus by the emission of a 3 particle is governed by the weak interaction, which
violates parity, and Wg(0) is written using the Odd-k Legendre polynomials.

28

described by terms with k = 0, 2, 4, and can be written as

VV(0) = l + GZBQF2Q2 P2(COS 6) + G4B4F4Q4 P4(COS 9), (2.3)

where the coefficients (I), are expressed as factors that separate the various physical
contributions to the final angular distribution. The factors come primarily from de-
tails of the 7 ray itself (F k depends on the multipolarity and electric or magnetic
character of the 7 ray, and the spins of the initial and final states) and the degree Of
initial nuclear alignment (specified by the substate populations P(m), from which the
orientation parameters 3;, are obtained). Also affecting the experimentally-observed
angular distribution is the finite angular resolution of the 7-ray detectors (instead Of
Observing W(9) at an exact angle 9, the finite opening angle of the detector averages
over a range of 6, which can “wash out” some Of the effect), represented by the Qk
factors. An additional effect is due to hyperfine interactions that reduce the orien—
tation of the nucleus if it leaves the target carrying electrons, prior to emitting the
7 ray. The hyperfine interaction introduces “attenuation coefficients” Gk that depend
on the state Of ionization of the projectiles emerging into vacuum [59]. The G k values
will be discussed along with the analysis Of the 7-ray angular distribution Obtained
in this work.

The degree of nuclear alignment is a measure Of the relative populations of the
magnetic substates, and plays a key role in the anisotropy of a 7-ray angular dist-
ribution. A convenient way Of treating an ensemble of oriented states (such as the
excited states resulting from Coulomb excitation) is to use density matrix formalism.
The density matrix is a way to organize all the information necessary to completely
specify the state Of a quantum-mechanical system. In the present application, that
consists Of the normalized populations Of the magnetic substates, using the beam
axis as the quantization axis. The density matrix p can be decomposed into multipole

terms that. affect separate k. terms in the angular distribution series. These multipole

29

terms are called “statistical tensors” pkq, and in the case Of axial symmetry (when
the particles are detected symmetrically about the beam axis) only the q = 0 terms

are needed: 1

pkg = T Z (—)’+m(1—m1m|k0)P(m). (2.4)

m=—I

PPM)

where the symbol :- \/2x—+—. In the angular distribution series, the statistical ten-
sors are contained in the Bk orientation parameters by the relation Bk = \/2k—+Tpk0.
The Clebsch-Gordan coefficients vanish unless k 5 21, Showing that higher-order
terms are not needed to describe the angular distribution completely. The Clebsch-

Gordan coefficients can be written explicitly for I = 2 and k = 2. 4 to obtain

2
1 2
=—-—-—— E 3m -—6Pm, and
p20 3m =¥2l l ( )

2

1 4 2
= 35m - 155122. + 72 P m .
“0 36\/14 mZ=_2[ ] ( )

 

Since P(m) = P(—m) from the definition of alignment, these expressions can be

further simplified to

p20=——(2)l2PP(1)-P0()la and
‘f— (2.6)

p40 = film) _ 4P(1)+ 3P(0)l-

For the case Of 100% prolate alignment, where all of the substate population resides
in the m = :tI states, the theoretical maximum values of the statistical tensors can
be Obtained from Eqn. 2.6 using P(2) = 0.5 and P(l) = P(O ) = 0 to get pl’mf =
0.5345 and pm“: 0.08909. This is the type Of alignment produced in intermediate-
energy Coulomb excitation. The opposite situation occurs when the population mainly
resides in P(O), called oblate alignment; this type Of alignment is produced in “head-
On” Coulomb excitation such as the Situation Shown in Fig. 2.1. For 100% Oblate

alignment, p.36” = —0.2673 and p36“ = +0.2673. It is easy tO see from Eqn. 2.6

30

that for no alignment (equal population Of all substates), )020 and p40 vanish3. In that
case, the angular distribution is isotropic in the frame of the emitting nucleus.

The 7-ray directional distribution coefficients are termed Pk, and depend on k,
the initial spin 1;, the final spin I f, and the multipolarity Of the 7 ray, through a

function called the F-coefl‘icient [58], defined as

I- I ,, I I ,,
Fk(LL’1f1,)=(-1)’f+’i+1kLL'I,-(L L k){’“ L A). (2.7)
1 —1 0 I,- 1,- If

where L and L' are the multipolarities Of the 7 ray, the (. . ) is a Wigner 3— j symbol
and the {. . .} is a 6-j symbol (see Ref. [60] and references therein). For a pure“1
multipole transition, e.g. a pure E2 transition, L = L' = 2. The F -coefficients can be
calculated directly from Eqn. 2.7, and they are also tabulated in the literature (see
e.g. Ref. [61]). For 2+ -—> 0+ transitions, F2 = —0.5976 and F4 = —1.0690.

The corrections to the ak coefficients due to the finite solid angles subtended by
the detectors, the Qk coefficients, have been extensively studied only for cases where
cylindrical 7-ray detectors are oriented with their symmetry axes pointing towards
the target position [62, 63]. AS will be discussed in Chapter 3, the 7-ray detectors
in this work were not oriented in that way, but instead had their symmetry axes
oriented tangentially to spheres centered at the target position. In the tangential
geometry, neither tabulated values of Qk nor analytical formulae known at present.
can be applied. An estimate Of the solid-angle correction factors can be used in lieu of
an analytical calculation, since the germanium detectors were not close-packed around
the target position but were 24.5 cm away. The estimate is Obtained by considering
the detectors to be oriented with their symmetry axes pointing towards the target

position, and performing a calculation in which the diameters of the detectors are

equal to the crystal lengths in the true geometry, and the lengths in the calculation

 

3A useful definition of the % prolate alignment is pgn/pg’éaf.
4For a 7 ray of mixed multipolarity, e.g. an E2+ M 1, see Ref. [58].

31

are set equal to the diameters in the true geometry. The 62,, values that result from
the calculation [64] are Q2 2 0.98 and Q4 = 0.94, which do not differ much from
unity. A more accurate calculation probably requires Monte Carlo techniques, but
is unnecessary since these Qk coefficients are only small corrections to the angular
distribution (on the order Of 2%, since the k = 4 term iS very small), and the statistical

errors of the measurements are generally much larger.

2.1.2 Transient Field Strength

The mechanism Of the transient field interaction is in principle understood, but at
present there are no means to calculate transient field strengths reliably from first
principles. In 1986, Rud and Dybdal made a thorough study [53] Of a mechanism based
on the polarization Of electrons bound tO the moving ions. that was originally proposed
by Borchers et al. [50] as an explanation for their results that discovered the transient
field. The polarized bound electron mechanism provides a coherent picture Of the
general behavior Of the transient field, but Rud and Dybdal admitted that the ability
to Obtain quantitative predictions from the model for the purpose Of making g-factor
measurements was limited, as it remains today. Nevertheless, there are numerous
properties of the transient field that are known empirically. The first property Of the
transient field that was apparent in early studies was that the direction is always
parallel to the applied magnetic field, even if the static hyperfine field in that same
substance is antiparallel. It was also observed that. the transient field depends on
the ion atomic number Z but not the mass number A, and it. also depends on the
species Of ferromagnet (e.g., Fe or Cd), and the velocity at which the ion traverses the
ferromagnet. By using very thin foils, experimenters were able to determine that the
transient field interaction “turns on” within 10 fs Of an ion entering the foil [65]. From
studies using ion probes of different velocities. it was Shown that the transient field is
usually a smoothly-varying function Of velocity, although for light. ions (Z = 8, 9, 10)

at low velocities. matching conditions permit. atomic shell and molecular orbital effects

32

to cause discontinuities in an otherwise smooth velocity dependence [66]. For some
heavy ions at low velocities, matching conditions for 4s electrons cause some field-
strength discontinuity in Fe hosts [67], analogous to the ls matching for the light ions.
Outside of the very low velocity regime, for light ions (Z S 20), the field strength

increases with ion velocity until reaching a maximum at a velocity of

vK=mR am

At much higher velocities, the TF magnitude approaches zero [68]. The quantity
110 = e2/h‘ :4 c/ 137 is the Bohr velocity of atomic electrons, and Z 210 is the velocity Of
the K-shell electrons of an atom of atomic number Z. It is assumed that the transient
field arises from an interaction between the polarized electrons in the ferromagnetic
host material with the ion as it traverses the material, and the velocity Of maximum
transient field strength for light ions indicates that the K -shell electrons Of the ion play
a significant. role. In fact, the magnetic hyperfine interaction between a Is electron and
the nucleus is the largest among all electron-nuclear magnetic hyperfine interactions
due to the high spatial overlap at the site of the nucleus. However, a theory to explain
the interaction fails to describe the transient field strength over a large range Of
ion Z values and velocities. For heavier ions (Z Z 20), the effectiveness of K -shell
polarization has not been fully explored experimentally and the existing data suggest
less efficient polarization Of K -shell electrons.

In general, it is not always possible to perform simultaneous relative measurements
between a reference 9 factor in a neighboring isotope to the nucleus of interest under
identical experimental conditions, in order to calibrate the transient field strength.
Therefore, different parametrizations have been established that are applicable to

different. velocity regimes and ion atomic numbers. Three parametrizations are found

33

in the literature, summarized in [52] as

Linear parametrization [69]: BL]N = aZf—y (2.9)
.,0
I) 0.45
Rutgers parametrization [70]: BRUT = a’Z 1'1 (—) M (2.10)
m
Chalk-River parametrization [71]: BCR = a’IZ—v—eT’av/“O. (2.11)
1,0

The value M in the Rutgers parametrization is the magnetization Of the host material,
and sets of the parameters a, a', a”, and fl exist for different mass regions and velocity
ranges (typically between 1120 and 5m). The Rutgers and Chalk-River parametriza-
tions have been extensively applied to the rare-earth region at typical velocities of
less than 5 v0. The linear parametrization is the most straightforward of the three
parametrizations, but quickly diverges from measured TF strengths as the ion veloc-
ity deviates from the velocity and Z range in which the coefficient a was determined;
an extrapolation to velocities beyond Z 110 would grossly overestimate the magnitude
Of the transient field. The Chalk-River parametrization includes an exponential term
that introduces a stronger velocity dependence, but even though it would Show a
disappearance Of the transient field at a sufficiently-high velocity (dependent on the
values Of a" and )3), it was not intended to be applied outside of the velocity or Z
range for which it was calibrated.

Because BTF is a function of velocity, it is important to know the velocity of the
ions as they travel through the ferromagnetic material. Ions lose energy while travers—
ing matter, primarily by ionizing the material along their trajectories, which causes a
continual decrease in velocity as they move through the material. The rate of energy

loss. called the stopping power, can be expressed by the Bethe-Bloch equation [72],

 

‘ EL 1 2nl-C’C'2I32‘7‘2Trriar _ ’32 _ _ 9 (2.12)
A 732 2 12 2

where z is the projectile atomic number. Z, A. and 6 are the atomic number, mass

34

number, and density of the target material, ,6 and 7 are the projectile velocity and
Lorentz factor, and me is the electron mass. I is the mean excitation energy which
can be approximated as (10 - Z) eV for Z > 16 and Tmax is the maximum energy that
can be given to a free electron in a collision with a particle of mass M and velocity

[3 = We
2mec2I3272
1 + 27me/M + (me/M)?

 

Tmar =

(2.13)

which can be approximated as ’1’me z 2mec2l3272 for 27mg / M << 1. K is a constant
equal to 47rNAr§mec2, where N A is Avogadro’s number, and re is the classical radius
of the electron, e2/me. The Bethe-Bloch equation, and stopping powers in general,
are an active area Of research, with correction factors such as shell corrections of
the target material usually taken into account [73]; for brevity, such details have
been omitted here. In practice, heavy ion stopping calculations are usually done with
computer programs based on the Bethe-Bloch equation or based on empirical stopping
powers. It is important to appreciate the functional form Of the stopping power,
and its dependence on projectile atomic number and target material properties, to

understand the principles of the transient-field targets presented in this Chapter.

'D‘ansient-field Target Principles

After the nuclear excited states are populated in an excitation target layer, the spin-
aligned nuclei pass into a ferromagnetic target layer that provides the transient field.
The ferromagnetic layer is polarized by an external magnetic field, which is usually
provided by a small electromagnet designed to produce a magnetic field Of 500-1500
Gauss (depending on the needs Of the ferromagnetic layer used); the purpose Of the
applied magnetic field is to align the magnetic domains Of the ferromagnet layer and
bring it to its saturation magnetization. Although in principle any ferromagnetic ma-
terial can be used to provide the transient field, two commonly-used ferromagnetic

host materials are iron and gadolinium, which have the most extensively studied tran-

10.0

 

 

 

 

 

E l

a, 7.5 .

E

g 5.0 . Fe

0 ~.

Ga “1......

DD “*4...

.5 2.5 ’ Gd ~~~m""'~~...., ‘

Q1 ~~~~~Ifl~~

9* l

.9

m 0.0 A A . A A . A A
0 500 1000 1500

Sulfur kinetic energy [keV]

Figure 2.3: Ziegler stopping powers for sulfur projectiles in iron and gadolinium, as a
function of energy.

sient field strengths. Iron is a good ferromagnet at room temperature, and has been
well—studied. The use Of Gd to provide the transient field is somewhat advantageous
because ion ranges are longer in Gd than Fe, and thicker targets can be used so more
time can be spent interacting with the transient field. The drawback is that Gd is
not magnetic at room temperature and must be cooled to below its Curie tempera-
ture (TC = 293 K) to enter the ferromagnetic phase, which requires a means to cool
the target (such as liquid nitrogen). However, most recent stable-beam measurements
have been made with Gd due to the advantage of its lower stopping power (compared
to Fe, shown in Fig. 2.3), since the measurements are done with beam energies near
the Coulomb barrier5.

The conventional transient field target usually includes a “stopping layer,” which

consists of a thick, non-magnetic material that provides a perturbation-free environ-

 

5The energy of the electrostatic repulsion of the two positively-charged nuclei, when they
are barely touching (the distance between their centers is the sum of their radii):

82ZIZ2 N 144-le2
R1 + R2 1.25041” + Ag”)

 

 

Ebarrier : MeV.

36

ment in which the nuclei can stop and then 7 decay. The stopping layer ensures that
the Spin precession from the transient field layer is “locked in.” The 7 rays are emit-
ted primarily while the nucleus is at rest, avoiding the complications arising from the
Doppler Shift.

Knowledge Of the transient field strength is necessary tO extract the magnitude of
the 9 factor from a measured spin precession. The transient field is always parallel to
the applied magnetic field direction, so the nuclear spins will precess in a direction
determined by the Sign of the 9 factor, independent of a precise calibration of the
transient field strength. The precession direction is straightforward to measure if the
general shape of the 7-ray angular distribution is known, because a Simple increase or
decrease Of counts in a 7—ray detector at a known angle with respect to the beam axis
will contain this information, as can be seen in Figure 2.1. Reversing the direction
Of the applied magnetic field will reverse the transient field direction, and therefore
the direction that the angular distribution rotates as well. Comparing the counts in
detectors for field up and field down is the main experimental signal that indicates
both the Sign and magnitude of the spin precession in a manner that cancels out most

instrumental asymmetries.

2.1.3 Measurement Of Spin Precession

The number Of counts detected for a particular angle and magnetic field direction,

e.g. up (T) is
N(i9 I) = (in/(6 IIAthbeama

where the efficiency of the detector at the angle i6 is (i, the counting time during
the particular field direction is AtT, and the average beam current during that time

is Ibeam- The sign convention for the detector angle is shown in Fig. 4.10. The double

37

ratio p is used in order to eliminate instrumental asymmetries, as seen here:

 

 

 

 

_ N(+0 T) N(--9 l) _ C+W(+6 DAtTIbeam C—u’r(-91)At11bea1n
N(+91)N(—6 T) C+l'V(+61)At11beam C- ill/(‘6 IIAtTIbea-m

(2.14)

 

 

_ W<+0 T) W(—0 1)
‘ W(+01)W(-61)'

The detector efficiencies, up/ down counting times, and average beam currents cancel,
since the two ratios of the double ratio are each taken with respect to one detector.
The precession effect is then

6 = —————-, (2.15)
which is related to the precession angle by

1 dww)

6
A6 =
(9) d0

S

 

, where S: W (2.16)

is the logarithmic Slope Of the angular distribution at the detector position.

The uncertainties in the quantities p and c are

 

 

4 2
p 6Ni 25p
5 =_ E __ . 6 = ,. 2.17
p Ziizlfm)’ and 6 (1..)2 ( )

Assuming that the counts in the two detectors for each field direction are similar, such
that there are N counts per detector per field direction, and using counting statistics

where 6N z \/ N, the uncertainty on the precession angle can be estimated. Inserting

6N = x/N into 6p gives

(2.18)

 

38

Using 6p 2 p/x/N, the uncertainty in the effect c is

26p _ 1 2p ~ 1
(1+p)3 «WM/03 NV

 

66 = (2.19)

where the last step is permitted since p z 1 for small precessions, which is the same
as the initial assumption that assigns N counts to each detector per field direction.

Since the precession angle is A0 = e / S, the uncertainty in the precession angle is

1
6 z .
A6 ZS r—N'

 

(2.20)

where the uncertainty in the slope factor 8 is generally considerably less than the
statistical uncertainty in c, and is neglected.

The inverse of rim), 5 2N, is Often called the precession sensitivity Of a particular
pair of detectors located at an angle 0. In general, 8(6) is the dominant factor in
determining where to place the detectors for a. transient field experiment, in order to

maximize S 2N.

2.2 g(2+) of Radionuclides: New Challenges

The transient field technique was developed at stable-beam facilities where measure—
ments on most stable isotopes with picosecond 2+ states were carried out. Several
groups began modifying the technique in order to measure excited-state 9 factors in
a limited number of radioactive nuclei. Since it is impractical to fabricate radioactive
targets, especially when the nuclei are short-lived, the radioactive nuclei need to be
produced as projectiles. In one method, the production of the radioactive nuclei was
done with transfer reactions in the first layer Of the transient-field target, and the prin-
ciple challenges are the radioactive decay background from the nuclear reactions. the
uncertainties the production mechanism introduces because of the unknown m-state

populations, and possible feeding from higher-lying states, which lends an uncertainty

39

to the interaction time for the state Of interest. Radioactive Ti isotopes were produced
using a-transfer reactions on stable calcium beams, with 44Ti (T1 /2 = 60 yr) Obtained
using a 40Ca beam [74] and 52Ti (T1 /2 = 1.7 m) Obtained using a 48Ca beam [75]. A
similar method was used to obtain 68Ge (T1 /2 = 271 d) by a-transfer onto a stable
64Zn beam [76]. Inspection Of the nuclear chart shows this method has limited appli-
cability to radioactive isotopes, because it requires a stable target with an isotopic
chain that spans a few units more on the neutron—rich or neutron-deficient side than
a neighboring isotopic chain (this is the case for the Ca isotopes and the Zn isotopes
mentioned here). Some heavier neutron-deficient isotopes can be accessed as the line
Of stability begins to favor larger N /Z ratios, e.g., 10°Cd—illOSn or 124Xe——4128Ba,
and under suitable conditions other candidates may be found.

TO measure 9 factors in radioactive isotopes available directly as beams, a tech-
nique based on Coulomb excitation in inverse kinematics can be used, which is one
of the more recent transient field experimental arrangements that offers many advan-
tages when the nucleus of interest is the projectile (for stable projectile applications,
see e.g. Refs. [77, 78]). The inverse—kinematics arrangement is outlined in Fig. 2.4.
Rather than using the beam to excite a target nucleus under study. a heavy ion Of
interest is accelerated and impinges on a thin target layer consisting of a light nuclear
species such as silicon. In a head-on collision between a heavy projectile and a light
target with the beam energy below the Coulomb barrier, the two nuclei will not touch.
However, due to momentum conservation, the light nucleus will be knocked forward,
followed by the heavy projectile. Since the light target ion is low-Z, it has a long range
in the multilayer target and escapes the stopping layer, where it travels through vac-
uum to be detected in a particle detector subtending a small angular range about the
beam axis. The heavy projectile, which may have been Coulomb-excited in the close
head-on collision, traverses the ferromagnetic layer in an excited state, experiencing
the transient field and precessing. Due to the higher stopping power for the higher—Z

projectile, it stops in the stopping layer and emits a tie-excitation 7 ray which is

40

 

0 Ge Detector
t.

5 x

O

 

 

Incoming
Projectile w
—> i A ' 4.
, 1 Particle
Detector

 

 

 

 

Si Fe Cu

Figure 2.4: View of the conventional transient field target used in conjunction with
Coulomb excitation in inverse kinematics. The incoming heavy projectile is excited
in a head-on collision, sending the light target nucleus forward to be detected in the
particle detector. The heavy projectile stops in the copper stopping layer, where it
can emit a 7 ray in a perturbation-free environment.

detected in coincidence with the light target ion that has recoiled out of the target.
Such an arrangement produces a clean spectrum Of 7 rays from Coulomb excitation
of the heavy projectile on the light target because of the particle coincidence required
(i.e., if the Coulomb excitation occurs in the ferromagnet or stopping layer, no light
particle will be detected in coincidence with the 7 ray). This geometry, which selects
head-on Coulomb excitation events, ensures that the excited nuclear states have a
high degree of spin alignment with the magnetic substate population dominated by
the m = 0 state.

A novel means Of obtaining a radioactive 7°Kr (T1 )2 = 14.8 h) beam at the
Berkeley 88” Cyclotron was developed [79], which used radiochemistry techniques to
separate the long-lived sample produced in a stable cyclotron beam bombardment.
The separated 7°Kr was then introduced into the ion source and accelerated for the ex-
periment (the “re-cyclotron technique”). In this case, the inverse-kinematics Coulomb

excitation technique was used. which overcame the drawbacks of using nuclear reac-

41

tions to populate the excited state. TO solve the problem of the decay background,
the stopping layer was replaced by a moving tape system that periodically removed
the activity. In this way, it was shown that a transient field measurement is feasi-
ble on radioactive ions available with Coulomb-barrier kinetic energies, which are
becoming available as re—accelerated beams at ISOL facilities such as ISAC-II [80],
REX-ISOLDE [81], and HRIBF [82].

Whether the radioactive isotopes are produced by nuclear reactions in the target,
or produced directly as a beam, the intensity (or yield) Of radioactive isotopes is
several orders of magnitude weaker than a stable beam intensity. The low intensity
directly lowers the yield Of 7 rays, and hence increases the statistical uncertainty
of the measurement. More will be said about the factors influencing the uncertainty
later, since this dissertation involves the specific case Of fast radioactive beams, which

come with additional challenges that will be introduced in the following section.

2.3 g(2+) of Fast Radioactive Beams:
New Technique

Most of the details of the first transient field measurements on radioactive isotopes
were specially designed for the nucleus at hand, and something that they all held
in common were the ion velocity ranges used. None Of the radioactive g(2+) mea-
surements were conducted with beam energies significantly different from previous
studies using the same elements, and therefore the well-established parametrizations
of the transient field strength with respect to velocity could be employed. To study
radioactive isotopes produced using in-flight fragmentation, which Offers access to
numerous isotopes away from stability but with higher velocities not encountered in
the previous transient field measurements. the new HVTF technique needed to be
developed .

The HVTF technique consists of three key components: the production Of a beam

42

Of radionuclides with sufficient intensity (intensity requirements are addressed in
§ 2.3.3); intermediate-energy Coulomb excitation to produce spin alignment, using
a thick target designed to degrade the projectile velocity to near Z v0 (target design is
addressed in § 2.3.3); and knowledge of the transient field strength, BTF(Z, II) which
depends on the ion and ferromagnet combination and is a function Of ion velocity that
is usually described by a parametrization (the BTF parametrization is addressed in

§ 2.3.2).

2.3.1 Intermediate—Energy Coulomb Excitation

When nuclei approach each other closely, they are each immersed in the strong electric
field created by the protons Of the other nucleus. Coulomb repulsion causes the nuclei
to deviate from their original paths and follow hyperbolic paths called Rutherford
trajectories. The shapes of the trajectories are hyperbolic because the electric field
is a function Of the distance separating the two nuclei, with a maximum at closest
separation and nearly no effect at large distances (hence the arms of the hyperbola
that asymptotically become straighter as the distance from the vertex increases). If
the nuclei approach each other swiftly, they will experience a strongly-varying electric
field as a function of time. This strong electric field is a time-dependent perturbation
that can cause a nucleus to transition to an excited state in a process known as

Coulomb excitation, Often abbreviated as Couler.

Semiclassical Approach

The strength and time-dependence Of the electromagnetic perturbation is calculable
because the interaction is the most well known Of the fundamental interactions. The

modified° Rutherford trajectory that governs the magnitude and time dependence Of

 

°Since the scattering is inelastic and involves an energy transfer to the excited nucleus, the
classical Rutherford orbit is slightly modified. Since one cannot specify where along the trajectory
the excitation takes place, a symmetrization procedure is performed to account for this modification
by setting up a new trajectory with the beam velocity changed to the average of the velocities before
and after the collision. Fortunately. since most of the excited states studied by Coulomb excitation

43

the electric perturbation is usually treated classically. This is justified because the
de Broglie wavelength Of the heavy ion is small compared to the nuclear scattering
length scale, or half the distance of closest approach in a head-on collision, a0. The

parameter a0 is defined as
Z1 22 e2
Clo = ——+'r—, (2.21
m043 2 )
for a beam velocity [3 = v/c and reduced mass mo of two nuclei with atomic numbers
Z1 and Z2. The ratio aO/A is called the Sommerfeld parameter 7), written explicitly

as
_Z1Z262
_ m3 ,

(2.22)
A classical description of the trajectory is justified when 7) >> 1. The actual tra-
jectory will determine the magnitude Of the electromagnetic interaction, which can
be decomposed into multipole interaction terms between the target’s charge and the
projectile’s electromagnetic moments. These terms are called the “monopOle—2’\-pole”
terms, such as the monopole—monopole term, which causes the Rutherford scattering.
As was pointed out in §1.2, only certain electromagnetic moments exist in nuclei for
symmetry reasons. Therefore, the monopole-dipole term is small since the nucleus has
no static electric dipole moment, and the electric dipole strength that exists due to a
dynamic electric dipole moment is concentrated at high energy in the Giant Dipole
Resonance. The most important term is therefore the interaction of the target’s elec-
tric field and the electric quadrupole moment of the projectile nucleus, which has the
consequence that the majority of the Coulomb excitation strength causes E2 transi-
tions. Magnetic interactions are suppressed by a factor Of (iv/c)2 relative to electric
interactions (which is not entirely negligible in the intermediate-energy regime; how-
ever, the electric quadrupole transitions will still dominate the Coulomb excitation).

For a nucleus subjected to the electromagnetic perturbation, the probability of

 

have only a few MeV of excitation energy, the energy transfer is small compared to the kinetic energy
of the beam. Therefore, the effect of the trajectory modification is not important, especially with
increasing beam energy.

44

Coulomb excitation will depend on the transition matrix element between the initial
and final states, and this probability can be calculated quantum mechanically through
perturbation theory. The total cross section for Coulomb excitation is calculated

semiclassically to be

dacoul = dORuth _ P-
d9 d0 'f’

(2.23)
where Pif is the Coulomb excitation probability for a transition to be induced between
two nuclear states by the electromagnetic perturbation Of a particular Rutherford
trajectory. Therefore, the measurement of the Coulomb excitation cross section yields
information on nuclear matrix elements, which give insight into the nuclear structure
Of the initial and final nuclear states. To extract the matrix element contained in Pif
from the measured excitation cross section, it is essential that the projectile and target
nuclei remain well separated throughout the collision so that the interaction between
them is purely electromagnetic. The usual way to ensure that the nuclei do not interact
with the strong nuclear force — which can happen if any matter overlap occurs — is
to choose the bombarding energy well below the Coulomb barrier. Specifically, a “safe
Coulex” energy can be defined such that the distance of closest approach between the
nuclear centers is the sum Of the nuclear radii plus a few extra fm, where the choice
Of how many additional fm of separation is left to the experimenter’s tolerance of
possible nuclear effects [83].

Coulomb excitation will generally occur only under the condition that the collision
is faster than the time it takes for the nucleus to adjust tO the changing electric field.
This is called the adiabatic condition, and is usually expressed quantitatively as the
ratio between the collision time ta,” z aO/v and the transition time ttmns = h/AE

(with AE = E] - Ei):
aOAE

 

g = (2.24)

The condition for which Coulomb excitation is probable is g < 1, and if the bombard-

ing energy is tOO low or the excitation energy is too high such that g 2 1, the Coulex

45

probability drops exponentially with g [84]. For typical Coulomb excitation experi-
ments, which run near Coulomb barrier energies, the adiabatic condition limits the
excitation to states with only a few MeV excitation energy. At such energies, a high-Z
nucleus can cause multiple Coulomb excitations of a nucleus under study, especially
if that nucleus is highly deformed, because Of the ease Of inducing E2 transitions in a
nucleus with a large quadrupole moment. In such a situation, it is no longer a weak
perturbative scenario and a system Of coupled-channels equations needs to be solved
to extract the many transition matrix elements involved [85].

When the beam energy is raised into the intermediate-energy range, several addi—
tional effects come into play. The primary effect is the rapid increase in the eccentricity
of the Rutherford orbits, which begin to look like straight lines, Since the Coulomb
repulsion can cause less deviation Of the high-momentum trajectory. The beam en-
ergy is above the Coulomb barrier, so the distance Of closest approach can be smaller
than the sums of the nuclear radii, and the projectile and target can interpenetrate.
Nuclear reactions can then take place, including few-nucleon removal and fragmenta-
tion reactions, which are undesirable because they result in nuclear species changes
rather than inelastic scattering. The experiment can be designed such that only par-
ticles that are deflected within a small angle 0mm; (corresponding to a large enough
minimum impact parameter bmm) are detected. A secondary effect at high beam en-
ergies is the Lorentz transformation of the electric fields Of the projectile and target,
as seen in Fig. 2.5, which increases the concentration Of electric field strength in the
plane perpendicular to the beam axis. This means that the projectile and target dO
not “feel” the influence of each other until they are closer than in the conventional
Coulomb excitation situation. Half the distance of closest approach, a0, is modified
by the Lorentz factor 7 = 1 / W to become a = a0 / 7, which reflects the reduced

effect of the electric fields prior to the scattering. The scattering angle in the center

46

of mass system, 00m, can be described classically by the relation

(1-
b=—°

I

cot. 560,... (2.25)

where b is the impact parameter, shown in Fig. 2.5. The collision time is very short,
Of order b/7v, and during that time there is an increased electric field strength due to
the Lorentz concentration Of the electric fields of the projectile and target. Multiple
Coulomb excitation can no longer occur, since the reaction is much more “sudden.”
The one-step excitation dominates, and can be solved using perturbation theory. The

adiabatic condition from Eqn. 2.24 is now

_a0AE AE l)

"t rhcfi Ear

 

(2.26)

where the distance Of closest approach in a head-on collision is replaced by the appro-
priate separation distance for small-angle scattering, the impact parameter b. From
Eqn. 2.26, one can see that fast collisions with b ~ R1 + R2 (the sum Of the nuclear
radii) can populate states Of order 10 MeV excitation energy. '

As was mentioned previously, intermediateenergy beams can easily cause nu-
clear reactions for impact parameters small enough to allow the nuclei to overlap.
In contrast to a Coulomb excitation cross-section measurement, a HVTF experiment
does not need to strictly ensure “safe Coulex” because the goal is not to extract
electromagnetic transition matrix elements. Therefore, the exact separation to allow
between the projectile and target (which is not a well—established quantity) is not
a major concern. It suflices to choose a maximum scattering angle 9mm such that.
the impact parameter b is larger than the sum of the nuclear radii, where if nuclear
effects set in, they should affect the excitation cross section without impacting the
distribution Of P(m) too severely. Since Eqn. 2.25 is defined for scattering angles in
the center of mass frame, 031?“: must be transformed into the appropriate lab-frame

0. The relativistic formulae for this transformation are given here. and can also be

47

:1) \]]])/—/e

 

4i

. //l\\ [M
WI)

Figure 2.5: The electric field lines from the projectile and target nuclei are contracted
along the beam axis by Lorentz effects for intermediate-energy beams. The projectile
has Z1 protons and mass A1, and the target has 22 protons and mass A2. The
impact parameter b is the distance separating the lines of initial trajectories, and the
projectile scattering angle is 6.

 

 

found in standard texts [86]. For a. projectile with mass number A and kinetic energy

T, the Lorentz factor 7 is
T

.__, . 2.27
+ AmN c2 ( )

7:1

with m N02 = 931.494 MeV. Using the ratio of the projectile mass to the target mass,
p = Mp/Mt, the Lorentz factors for the projectile and target. can be defined for the

center of mass frame,

 

 

”Y + P
”I = . 2.28
l (1+2m+p2)‘/2 ( )
and
1+ '
’72 = ”>10 (2.29)

(1+ 2m + p2)1/2‘
Eqns. 2.28 and 2.29 can be used in the expression for the lab frame scattering angle

0Lab that corresponds to the center of mass scattering angle 66",:

 

1 sin 6
tan OLab = — ° cm

. 2.30
)2 ((3089017) + Alp'Tl/Alt’l‘Z) ( )

48

2.3.2 Transient Field Strength at High Velocity

Few studies have been undertaken to study the action Of the transient field at ve-
locities v > Z v0 for ions with Z S 16, although a great number Of experiments
have been conducted at ion velocities between the L and K Shell electron veloci-
ties (éZvo < v < Z '00). Measurements Of Spin precessions using stable beams with
v at Z120 have been achieved for example using sulfur projectiles [87] at the Max—
Planck-Institut fiir Kernphysik in Heidelberg, and a higher-Z Cr beam (Z = 24) at
752 MeV, corresponding to v = 21120 was produced with the UNILAC accelerator at
Darmstadt [88], lending hope to the idea that high velocity measurements are pos-
sible. However, radioactive beams produced by fragmentation have velocities on the
order Of 0.3-0.5 c, and the K shell electron velocity is Z co z 3276" which is much
less than 0.3 c (especially for ions with Z = 16 in this work). There is no a priori
reason why any Of the conventional TF parametrizations could be extended into this
velocity regime; indeed, it would seem that no TF interaction is possible at such high
velocities. Nevertheless, there has been great interest in the use Of fast radioactive
beams for transient field g factor measurements, and several groups have reported
their approaches to the problem [89, 90]. The primary concern is that there would
be a sufficient precession effect to make a HVTF technique feasible for radioactive
beams, which Often are only available with low intensity.

A parametrization Of the transient field strength is needed in order tO extract
9 factors from possible spin precessions. The parametrization can then be used tO—
gether with the other experimental details such as the intermediate—energy Coulomb
excitation cross section and the degree Of spin alignment, the angular distribution
Slope parameters at the detector locations, and the 7-ray detector efficiencies to esti-
mate the experimental sensitivity tO the spin precession expected in the high velocity

regime.

~19

 

 

 

 

2 I I r
._. IO Fe
r3 tNe
:93 'Mg
1 _ oSi ‘
M
a *8
“.4
JJ
CD
0 I I 1
0 1 2
v/Zvo

Figure 2.6: TF strengths for light ions in Fe (data points) and high velocity TF
parametrization that results from a fit using the functional form described in the text
(solid line), from Ref. [91]. Also Shown is a linear parametrization used in this mass
and velocity range.

Parametrization for High Velocity Light Ions

A parametrization was proposed by Stuchbery [91] after an examination Of the avail-
able data for transient field strengths for ions with 6 5 Z S 16, traversing Fe or
Gd foils, at velocities greater than $2140. A brief description of this parametrization,
which was used to plan the experiment on 38’“OS and analyze the resulting data, will
be given here. In principle, Since high-velocity light ions will exist primarily in high
charge states (e.g., q = Z — 1) while traversing matter, the physical situation is sim-
plified relative to heavier ions, and it is plausible that a model-based parametrization
can be found to describe the transient field behavior. The basic physical principle for
the model-based parametrization, here called the ANU parametrization, is that the
transient field strength for a particular ion as a function of velocity should follow the
hydrogen-like charge-state distribution for that ion as a function of velocity, except
that the peak of the field strength should occur at v = Z v0. This is a reasonable
assumption because it is well-established that. the K shell electrons are responsible

for the majority of the transient field effect. for light ions [53]. The magnetic field at

50

the site Of the nucleus, due to a Is electron in particular [92], is
31, = 16.7Z3R(Z) T, (2.31)

where the relativistic correction factor, B(Z) 2 [1 + (Z /84)2'5], is not appreciably
different from. 1 for light ions. This expression is Often used as a starting point in

describing transient field strengths, which are generally described as [51, 52, 53]
3,, o< {1.0, Z)Fl.<v, 2.. 298142). (2.32)

In the above expression, {13 is the polarization Of the Is electron carried by the
ion and F113 is the fraction Of ions that are carrying one 13 electron (this is also
called the “single K —vacancy fraction”). For the development Of the parametrization,
Stuchbery used the Observation that 513 is virtually velocity-independent and wrote
the polarization factor {13 as {18(Z) = AgZpé, where Ag and Pg are parameters that
can be different for each ferromagnetic host material. The K -vacancy fraction F113
was used to describe the velocity-dependence Of the transient field, using a function
that reproduced the behavior Of the hydrogen-like charge state fractions for light
ions traversing Fe and Cd foils, calculated with the program LISE [93, 94]. The

adopted function has the form f (x) = axe—bx

2
, where a: = E /A is the kinetic energy
per nucleon. Using this functional form together with the Z—dependent form for the

polarization parameter and Bis, the AN U parametrization is
1 I i 4
Bt,.(v,Z) = AZP(u/Zvo)Qe—?(L/ZIO) , (2.33)

with the parameters A and P determined for Fe and Cd hosts by fitting the available
data on transient field strengths. For Fe hosts. summarized in Fig. 2.6, the value of

. 2 182(5) T with P = 3 was Obtained [91].

 

 

 

 

 

 

 

I T I
'3'
3 30-— —
NN
\ 20— J
a 10- _
I3
0' . ' ‘
o 1 2

v/Zvo

Figure 2.7: Compilation Of previously-known TF strengths for light ions in Gd (solid
symbols) and high velocity TF parametrization (solid line), and newer measurements
that test the parametrization (Open symbols), from Ref. [95]. A linear parametrization
(dashed line) is also shown.

Stable Beam Transient Field Strength Measurement

When the light-ion HVTF parametrization was proposed, it was apparent that very
few data on transient field strengths beyond 12 = Z 220 exist, and that more experimen-
tal work was necessary to validate the parametrization. An experiment was performed
at the 14UD Pelletron facility at The Australian National University as a collaboration
between researchers at the NSCL and at ANU, to measure transient field strengths
for Mg ions in Gd, at velocities above Z 00 [95]. Gadolinium was chosen as the ferro-
magnet because there are fewer “high-velocity” field strength measurements reported
in the literature (where high velocity means I) > fiZvO) for gadolinium compared to
iron. In addition, a v = 1.25Zv0 measurement (the highest velocity with which a TF
measurement has been made for ions with 6 S Z S 16) has been reported for neon in
Gd, with large error bars. In Figs. 2.6 and 2.7, the previously—known TF strengths are

shown as filled points, and the v > Z140 Ne point leaves open to questioning whether

the transient-field strength is increasing or decreasing at v = 1.25Z-00. A confirmation
of the proposed parametrization in gadolinium would provide strong evidence for the
validity Of the proposed iron parametrization, which was based on a fit to more data
points than gadolinium. In order to ensure the validity of the result, a low-velocity
point for Mg in Gd was measured in the region where the transient field has been well
characterized, 1) ~ éZvo, and the transient field strength extracted from that mea-
surement agreed with systematics of known transient field strengths. The new AN U
measurements are Shown in Fig. 2.7 as open points. The new measurements consist
of two high-velocity and one low-velocity Mg TF strength measurements, along with
an additional point taken with a Ne probe near éZvO [95]. The previously-measured
transient field strengths (filled points) and the ANU parametrization [91] (solid line)
are shown for comparison. Good agreement is seen between the new high velocity data
and the parametrization for Gd hosts, and since the transient field behavior is Fe is
similar, the Fe parametrization should also hold. Further experiments are planned to
measure high velocity transient field strengths in Fe at AN U for additional confirma-

tion.

2.3.3 Target Design for Optimum Slowing

The concentration Of transient field strength near Z en provides a stringent constraint
on the experiment design. Slowing the fast radioactive projectiles to near Coulomb-
barrier velocities is the most critical aspect to be addressed in the choice of the
experimental parameters, i. e. primary beam energy, primary target thickness, wedge
degrader thickness, and HVTF target layer thicknesses, because without a proper
Slowing of the projectiles, there is no transient-field induced spin precession. Accom-
plishing this slowing in an Optimum fashion involves several compromises to ensure a
successful experiment. Some aspects Of the compromises that affect the experimental
sensitivity to spin precession can be quantified by figure Of merit expressions involv-

ing the number Of 7-ray counts and the slope Of the 7-ray angular distribution at

53

the 7-ray detector positions, as discussed in § 2.1.3. To first order, however, these
quantities only impact the amount Of beam time required to achieve a statistically
significant result. Of greater importance is the necessity to Slow the beam down into
the velocity range where BT F can act at all.

The requirement that the projectiles be Slowed to v x Z v0 (which corresponds to
~5 MeV/ nucleon for the heavy sulfur isotopes) is nearly the equivalent to requiring
the ions to Slow as much as possible without stopping in the HVTF target. Slowing
ions in matter necessarily increases the energy spread of the projectiles. A compromise
must therefore be made between how much slowing can be done in the fragment
separator (and the fragment transmission that lowers as the momentum spread Of the
beam increases) and the amount Of slowing that gets left for the HVTF target itself.
The practical limit (Emu) — 0E > 0 Should be obeyed, where a E is the full width at
half maximum Of the energy distribution for the emerging ions. An estimate of the
increase in energy straggling is given by Gottschalk, where for an incoming energy

spread oi, the outgoing energy spread is approximately [96]

00 = —— - 0,, (2.34)

where S (E) = —dE/d:r is the energy—dependent stopping power for the ion/ material
combination of interest. Using Eqn. 2.34 and an incoming energy spread, the com-
promise between slowing in the separator and slowing in the HVTF target can be
quantified. In general, a simulation Of the fragment separator and ion optics is neces-
sary to calculate beam properties after production, separation and transport.

An additional constraint on the experiment design is the slowing down time. In the
HVTF target, the slowing must be done quickly because after the Coulomb excitation
event, the states begin decaying. Decays that occur in the non-ferromagnetic layer

accomplish nothing because there is no transient field present. However, since the

2

target? it is advantageous to use a thick,

cross section for Coulomb excitation goes as Z

54

high-Z Coulomb excitation layer to ensure enough excitation yield. The thickness of
the Coulex layer increases the transit time of the ions prior to their arrival in the
ferromagnet, so a compromise needs to be made between excitation yield and the
need to ensure a large interaction time with BTF prior to the decay governed by
the excited-state lifetime 7'. It turns out that iron has a higher stopping power than
many common Coulex targets such as gold, so if iron is used, it is beneficial to place
the transition to the ferromagnetic layer at a point where the ion velocity is still
beyond 22120; there is no transient field initially, but the ions are quickly slowed into
the velocity range where BTF is large. To determine whether a particular nuclear
excited-state 9 factor can be measured using the HVTF method, the slowing time
needs to be compared to the excited-state lifetime.

The slowing time can be estimated using tabulated ranges of ions in matter,
available, for instance, with the program SRIM [97]. As a benchmark, the time it
takes for ions with 6 3 Z s 20 to slow in a single step from 100 MeV/nucleon
down to 5 MeV/ nucleon in aluminum (a typical degrader material), is shown in the
left panel of Fig. 2.8. These values were calculated for each proton number Z (for
simplicity using each isotope’s stable mass value) by looking up the ranges in A] for
100 MeV / nucleon, R 2(100), and subtracting the range of the ion with 5 MeV / nucleon,
RZ(5), to get the total slowing distance, R g(AE) = Rz(100) — RZ(5). To obtain
the slowing down time, the energy loss was approximated linearly through the target
and the velocity was approximated as 21(1) = W. The inverse velocity as a
function of energy was integrated through the distance R g(AE). As a comparison
to this “one step” slowing, the times were calculated for an iron transient-field layer,
assuming that the ions have first been slowed by some means to 20 MeV/nucleon
prior to the iron layer. The right side of Fig. 2.8 shows the resulting times for ions
to slow from 20 to 5 MeV/ nucleon in iron, using range tables for ions from carbon to
calcium and the linear energy loss approximation. It can be seen that the times range

from ~2-5 ps, meaning that excited states with lifetimes on the order of a few ps

55

 

 

 

 

 

 

 

 

150 ......... 6 r a f T - r
.—. Time to slow from _ 5 _ Time to slow from
a 100 MeV/u to 5 MeV/u a 20 MeV/u to 5 MeV/u
“3 '00 0 . in aluminum ‘5‘ 4 _ in iron
g 0 . é o
g...
bl) . O E; 3 . 0 o
.5 . 0 . .E . . o . .
g 50 0 o . . g 2 . . .
"r13 :75 l .
0 1 nnnnn 1 - l A 1 0 L 1 - 1 1 1 L4 1 1
6 8 IO 12 l4 16 18 20 6 8 10 12 I4 16 I8 20
Projectile Proton Number Projectile Proton Number

Figure 2.8: Approximate time for ions with 6 3 Z S 20 to slow down to
5 MeV/nucleon, in one step from 100 MeV/nucleon in Al (left) and as a second
step from 20 MeV/nucleon in Fe (right); see text for calculation details. Note the
change in vertical scale between the two scenarios.

are accessible, if the transit time in the Coulomb excitation layer is also kept to a
few ps. Since the slowing from 100 MeV/nucleon —— the typical production energy
—— has been shown to be much greater than tens of picoseconds, it is clear that
the Coulomb excitation layer cannot be responsible for slowing the ions from 100 to
20 MeV/ nucleon, and the bulk of the slowing must be done first in the separator.
Therefore, the general technique of transient field measurements on fast beams
is a “three-step slowing” process, where the ions are slowed using degraders in the
fragment separator, and then slowed in a two layer HVTF target, where a high-Z
Coulex layer slows the ions while populating the excited states, and then the transient
field layer slows the ions to near or below Z‘vo. It has been shown that an iron
transient-field layer can accomplish slowing from 20 MeV/ nucleon to 5 MeV / nucleon
in 2—5 ps and is therefore a better choice than gadolinium, which has a lower stopping
power and therefore a longer stopping time. It is possible, however, to use a single-
layer HVTF target where gadolinium (Z = 64) serves as both a Coulex target and the
ferromagnetic host. This choice may have a smaller effective time between the average
excitation and the onset of BTF under certain experimental conditions, but it is a

compromise between transient field interaction time and the Coulomb excitation cross

56

 

section, in addition to the added complexity of requiring a cooling mechanism for the
gadolinium.

The beam intensity requirements for the fast-beam transient field method depend
on the Coulex yield, integrated BTF, state lifetime T, 'y-ray detection efficiency and
angular distribution, and detector positions. Intermediate-energy Coulomb excitation
cross sections are usually of order tens of millibarns (depending on the B(E2) value,
target Z, beam energy, and maximum scattering angle 0mm), and thick targets of
order hundreds of mg/cm2 can be used. The 7 rays are emitted with a probability
that decays exponentially in time after excitation, meaning that a certain fraction
will be emitted prior to reaching the ferromagnet (or prior to the projectile velocity
being near Z110 where BTF has any effect). The usual figure of merit to maximize in
conventional TF experiment planning is S 2N, but a different figure of merit should
be used for fast beams due to some different conditions. Due to the long stopping
time for fast beams, the only '7 rays that contribute to the number N are those that
are Doppler corrected properly into a peak area that can be accurately obtained. For
practical reasons that will be explained in Chapter 4, those 7 rays are emitted after
the projectiles leave the target and are flying with constant velocities. Therefore, the
number N needs to be reduced by a factor that depends on the excited state lifetime,
7', and the transit time through the target, T: N ——> N 6‘77 T. The figure of merit that

can be used to determine the required beam intensity is therefore
82¢2l'lx'ae'T/T. (2.35)

with q‘) = AQ/g representing the g-factor dependent term (a larger 9 factor requires
less time to measure), S is the slope and W is the angular distribution at a particular

detector angle, and a is the excitation cross section.

Chapter '3

Application at NSCL

The radioactive 38408 isotopes under study were produced using in-fiight fragmen-
tation at the NSCL’s Coupled Cyclotron Facility. The high kinetic energy of the
resulting fragments created challenges to the g(2+) measurement that needed to be
overcome as described in Chapter 2. The experimental technique and equipment used
to accomplish the g(2+) measurement. for 38‘408 at the NSCL are described in this

Chapter.

3. 1 Isotope Production

Two different stable isotopes were accelerated as primary beams to produce the de-
sired radioactive sulfur isotopes. To produce 388, enriched 40Ar gas was introduced
into a superconducting Electron Cyclotron Resonance (ECR) ion source, leading to
40Ar7+ ions that were guided into the K500 superconducting cyclotron and acceler-
ated to 12 MeV/ nucleon. The accelerated 40Ar7+ ions were sent through the coupling
line and radially injected near the center of the K1200 cyclotron, where a thin carbon
foil stripped them of their remaining electrons. The K1200 was then used to accelerate
the 40Ar18+ ions to 140 MeV/nucleon, at which energy they were extracted and di-

rected onto a thin target of beryllium located close to the exit of the K1200 cyclotron.

 

 

 

 

 

 

 

 

 

Figure 3.1: Schematic of the coupled cyclotrons and A1900 fragment separator.

Fragmentation reactions occurred inside the target, creating a wide range of nuclides
lighter than 40Ar, by removing several nucleons and leaving the products flying for-
ward with a similar velocity to that of the incoming primary beam. For instance, 388
was created through the removal of two protons from the 40Ar projectiles.

The fragment separator, called the A1900 [98], was necessary to separate the de-
sired fragment from among the variety of products. The A1900 consists of a symmet-
ric arrangement of four superconducting dipole magnets, and several superconducting
quadrupole triplet magnets, and has a large acceptance and high resolving power. The
target box after the K1200 cyclotron is located at the object position of the A1900,
which works by separating the reaction products according to their momentum to
charge ratios. A diagram of the A1900 is shown in Fig. 3.1, taken from Ref. [98]. At
the second intermediate image, there is a large dispersion of the fragments in the hor-
izontal direction. At this position a wedge-shaped degrader and a slit system can be
used to selectively enhance the transmission of the species of interest. The second half
of the A1900 then focuses the fragments back for transmission to the experimental
areas.

The 408 ions could not be created by removing nucleons from ‘mAr, so in a separate
run, a primary beam of 48Ca was created by heating a sample of enriched 48Ca oxide

material to ~875° C in an oven connected to the ECR ion source. The 48Ca vapor

59

Table 3.1: Production and properties of the radioactive beams.

 

 

 

Primary beam 9Be Target A1900 Settings Secondary beam

Ion Intensity Thickness Bplg Bp3,4 Ion E Intensity Purity
(pnA) (mg/CH9) (Tm) (Tm) (MeV) (PPS) (‘76)

40.41 25 987 3.5303 2.2050 388 1547.5 2.105 >99

480a 15 1034 3.6580 2.2872 408 1582.5 2104 >95

 

 

 

 

 

 

was ionized, and the ions in the 8+ charge state were accelerated through the coupled
K500-K1200 system to 140 MeV/ nucleon, and the resulting fully-stripped 48Ca ions
impinged on a Be target. The production of 403 occurred through the removal of four
protons and four neutrons from 48Ca.

The primary target thicknesses and A1900 magnetic rigidities are summarized
in Table 3.1. An acrylic plexiglass wedge degrader, 971 mg/cm2 thick, as well as a
0.5% momentum slit, were employed in the second intermediate image position. The
resulting kinetic energies of the 38’408 fragments were approximately 40 MeV / nucleon,
which was found to be the minimum energy to which the secondary beams could be
slowed before the transmission through the beamlines would suffer due to energy
straggling in the slowing down process. The HVTF target was therefore specifically

designed for 38"10S ions of energy ~40 MeV/ nucleon.

3.2 Multilayer Target

Using stopping powers of Ziegler and the fragment separator simulator LISE [93, 94],
a multilayer target was designed to perform the HVTF measurement on the 38"10S
isotopes. The target materials and thicknesses needed to be chosen carefully to ensure
a. successful measurement, as described in § 2.3.3. The program LISE was used to
estimate the secondary beam energy for 388 fragments produced by the fragmentation

of 140 MeV/nuclear 4OAr on a 9Be target. A pure beam was desirable for the first

HVTF measurement to minimize unwanted background. but in general it is not a

60

stringent requirement. In the case of the 9Be(40Ar,38S)X fragmentation reaction it
was possible to set the A1900 fragment separator settings to obtain a > 99% pure
388 secondary beam with an energy of ~40 MeV/ nucleon. Using this beam energy
as a starting point, it remained to design an HVTF target, assemble the target by
rolling the target layers to the desired thicknesses, and then to confirm the expected
energy loss by testing it in-beam prior to the g-factor measurement. The goal was
to choose a target that would slow the projectiles down to as low of a velocity as
possible without stopping, or around 100-200 MeV.

The Coulomb excitation target layer was chosen to be gold due to its high atomic
number, Z = 79. Since A11 is monoisotopic in stable form, there will be less target
7 rays from different levels in the target nuclei Coulomb excited by the beam. The
malleable nature of Au is also beneficial, and allows it to be easily rolled to the
proper thickness. For the transient field layer, iron was chosen as the ferromagnet
since it does not require cooling with liquid nitrogen or other means, and as was
shown in Chapter 2, it has a high stopping power, which is important to ensure that
the projectiles slow as quickly as possible. Natural iron has 4 isotopes, 54’56’57’58Fe,
and the most abundant isotope, 56Fe, has a ’y-ray transition at 846 keV, which is near

the 903 keV first-excited state energy in 40S.

3.2.1 Stopping Powers and In—Beam Test

The Au and Fe thicknesses were chosen based on calculations with the program TRIM,
which uses the stopping powers of Ziegler and coworkers [99]. The Au thickness was
chosen to be 355 mg/cm2, which was calculated to slow the sulfur ions to kinetic en-
ergies of ~780 MeV. For a comparison, calculations using the ATIMA [100] stopping
powers implemented in the program LISE predicted a kinetic energy of ~810 MeV
when Lindhard-Sorensen (LS) corrections [101] were included, and ~780 MeV with-
out LS corrections. The optimum Fe thickness was determined to be in the range

100—110 mg/cm2 as a result. for which predictions of final energies of 130, 216, and

61

131 MeV were obtained using stopping powers of Ziegler, ATIMA+LS, and ATIMA,
respectively, for 110 mg/cm2 Fe.

The two foils were purchased from Goodfellow Cambridge Limited, a 25 mm x
25 mm Au foil with a nominal thickness of 0.250 mm and purity 99.99570, and a
50 mm x 50 mm Fe foil of nominal thickness 0.125 mm and purity 99.99970. The Au
foil was rolled to 355 mg/cm2 using the metal roller in the NSCL detector lab, by
sandwiching the Au foil between clean strips of stainless steel (by folding one strip
in half). The sandwich of steel and Au was passed through the roller several times,
reversing the orientation of the sandwich each time in order to reduce any systematic
thickness variation across the Au. The Fe foil was cut into 4 squares of ~25 mm2,
and three of the pieces were rolled to thicknesses of 90, 100, and 110 mg/cmz.

It was essential to test the two-layer target in-beam to choose the optimum Au / Fe
target combination, for two reasons. The stopping powers are not known to high
accuracy for the slowing down of 388 from ~1500 MeV to ~150 MeV, as shown
above. In addition, the exact initial beam energy would be unknown prior to the
actual beam extraction from the cyclotrons. For the target test, a primary beam of
40Ar was accelerated to 140 MeV/nucleon and impinged upon a 945 mg/cm2 Be
production target. An acrylic wedge degrader of thickness 971 mg/cm2 was used at
the dispersive image of the A1900 to separate the desired ions from other species
produced in the fragmentation reaction. The secondary beam of 38S was transmitted
through the second half of the A1900, with the magnets set to a magnetic rigidity of
Bpg = 2.2053 Tm, corresponding to a kinetic energy of 1548 MeV. The Au target was
held on a rotating target ladder in order to change the effective thickness according
to t’ = t/ cos 0, where the normal target thickness is t and the thickness along the
beam axis is t’.

Two particle detectors were used, a silicon PIN diode of 980 pm thickness (for
measuring the after-target energy of the 388 ions) and a plastic scintillator phoswich

detector (for identification of possible fragmentation products in the thick target),

62

 

250 . . r T
200 r
150 ~

100]

Counts / MeV

50-

    

 

 

 

0 50 100 150 200 250 300
Particle Energy [MeV]

Figure 3.2: Representative PIN spectrum taken during the in-beam target test run,

using a 355 mg/cm2 Au target layered with a 110 mg/cm2 Fe target, positioned
normal to the beam axis. The PIN energy is calibrated (see text), and the 383 centroid
energy is 140 MeV.

consisting of a 0.75 mm layer of fast BC-400 plastic scintillator and a 2” layer of slow
BC-444 plastic scintillator. The phoswich was viewed by a single photomultiplier
tube (PMT), which was coupled to the scintillators by an acrylic light guide. The
light guide and PMT had a thin layer of optical grease to couple them. The slow
scintillator was attached to the light guide with optical epoxy, and the fast plastic
was simply placed in contact with the slow scintillator and held in place by a thin
wrapping of aluminized mylar for light tightness. The scintillators were disks with 4”
diameters, and the PMT was a Hamamatsu PMT / base assembly with a 2” diameter.
More details about the operation of phoswich detectors and the description of the
actual phoswich detector built for the HVTF g-factor measurements will be given in
§ 3.4.

The PIN detector was calibrated using a 228Th (1 source as well as AE sig-
nals from attenuated primary beam having a known Bp, which were sent through
the PIN without the AuFe target in the beamline, from which the calibration E =

(0.2253/ channel 5 X + 3.683) MeV was deduced to convert from channel number X to

63

energy. The PIN detector measurements showed that the 388 energy was ~800 MeV
after the Au layer, and that the ions exited the Au and 110 mg/cm2 Fe foil combi-
nation in a distribution between 80 and 200 MeV, with a peak value at ~140 MeV,
seen in Fig. 3.2. The phoswich detector showed that the majority of the projectiles
emerged from the target without undergoing nuclear reactions. The test established
that the Au and 110 mg/cm2 HVTF target was appropriate for slowing the beam

into the peak of the transient field strength.

3.3 Coulomb Excitation of 38ADS Isotopes

3.3.1 Segmented Germanium Array

The Segmented Germanium Array (SeGA) is an array of n-type high-purity germa-
nium solid-state 'y-ray detectors [102]. Each germanium crystal is cylindrical, with a.
length of 8 cm and a diameter of 7 cm. Each crystal has a closed-end coaxial geom-
etry with a lithium-diffused central contact to detect the full energy signal from the
crystal. The outer surface of the crystal is an ion-implanted p—type contact, which is
etched into four azimuthal segments and eight longitudinal segments, giving a total
of 32 segments per crystal. Fig. 3.3 shows the segmentation geometry, and Fig. 3.4
shows a view of a complete SeGA detector, which has the crystal cryostat mounted
at a 45° angle to the liquid nitrogen dewar in order to allow detector positioning

with the fine segmentation axis facing the beam axis for Doppler correction purposes,

ABCDE FGH

 

A2 A1

I—n—I—d-J --lb— In—

 

 

A3 A4

 

 

 

 

 

 

 

 

 

 

L

Figure 3.3: SeGA detector crystal segmentation naming convention.

 

64

 

 

ct; Derr‘

 

\W/

 

 

 

 

EJ

Figure 3.4: The SeGA detector consists of a liquid nitrogen dewar and a thin-walled (<
1 mm) aluminum cryostat containing the single Ge crystal. The Field Effect Transistor
(F ET) compartment is situated near the crystal compartment but outside of the
cryostat so that the FETs are at room temperature. The preamplifiers are located in
the angled section on the dewar.

which are discussed in §4.2.2. The signals from the segments can give information 011
the 'y-ray interaction points within the crystal volume, and a measure of the energy
deposited by the 7 ray in each segment that was hit. The segment energy signals
were calibrated relative to the central contact energy using a procedure outlined in
Ref. [103]. Conventional electronics were used and the ”y-ray interaction points were

assigned to the centers of any hit segments.

Angular Distribution Configuration

A stand was constructed at the NSCL for the purposes of SeGA 'y-ray angular distri-
bution measurements using intermediate—energy beams [104]. The stand was designed
so that SeGA would be sensitive to the angular distribution as a function of the polar
angle 0. Since the scattered projectiles are detected symmetrically about the beam
axis, the 7—ray angular distribution is insensitive to the azimuthal angle 05. However,
the precession induced by the transient field occurs about the applied magnetic field
axis, which is oriented along the vertical axis. Therefore, the sensitivity of the Ge
detectors to the precession is maximized in the plane perpendicular to the magnetic
field, and decreases to zero sensitivity at 45 = 0° or 1800 (Ge detectors located directly

above or below the beam axis). The Cartesian coordinate axes are shown in Fig. 3.5.

65

Magnetic field axis
+x, (1) = 0°
A

Detector

I
I .
a s
I N‘
I “
. s
t “_
a ‘.
0 ‘~
‘._
.
A ~-
~‘
. .‘
‘~ 9 9
I ‘Q
s
.
‘Q
.
y. .
u. ,
n. I
‘s I
.
. I
‘~ 0

 

 

Beam axis
+2, 0 = 0°

Figure 3.5: Relation between beam axis. magnetic field axis, and the Cartesian and
spherical polar coordinate systems used in this dissertation.

 

For the TF g-factor measurements, sensitivity also increases with the slope of the ”y-
ray angular distribution. Intermediate-energy Coulomb excitation populates excited
states predominantly in the m = 21:] magnetic substates, and in the specific case of a
2+ excited state transitioning to a 0+ ground state, the transition is necessarily dom-
inated by Am = :l:2 (since the 0+ state has only the substate m. = 0). A Am = i2,
2+ ->0+ angular distribution has its greatest slope near 0mm = 25° and 6mm = 155°,
where 0mm is the angle in the projectile frame measured from the downstream beam
axis. The shape of the angular distribution for Coulomb-excited intermediate-energy
beams always has a zero point in the slope at 0 = 90° in the projectile frame, so it is
advantageous to avoid 90° and cluster the detectors at angles near the beam axis.
The SeGA stand design for intermediate-energy beam 'y-ray angular distributions
takes into account the large degree of Lorentz boosting that focuses the y rays forward,
so that even an isotropic 'y ray in the projectile frame will have fewer counts in the
backward hemisphere and more counts at forward angles in the lab frame. In the
standard SeGA angular distribution configuration, the detectors in the horizontal

plane were placed as seen on the left. side of Fig. 3.6. The back—angle detectors were

66

 

Figure 3.6: The detectors in the horizontal plane for the intermediate-energy ’y-ray
angular distribution measurements (left) and the TF g—factor measurements (right).
The beam axis, indicated by the arrow, determines 0 = 0°.

located at 0 = 127° and there are also detectors at 78°. The TF g-factor measurements
are done with slower beams, and the Lorentz effect is not as pronounced as it is for the
angular distribution measurements with full-energy beams. Therefore, an alternate
configuration was designed in which the same stand can have the 127° detectors
“flipped” so that the crystals are located at 6 = 151°, close to the optimum slope.
Part of the TF configuration is seen on the right-hand side of Fig. 3.6, which shows the
detectors located in the horizontal plane. Because of the geometry of the SeGA liquid
nitrogen dewars, it was not possible to include the 78° detectors in this configuration.
Since the 78° detectors are close to 90°, they would ofier very little sensitivity to the
spin precession and their removal would not impact a g-factor measurement.

The full SeGA configuration for the TF g—factor measurement is shown in Fig. 3.7.
There are two detectors for each position in 0. located at azimuthal angles 4) and
(15+ 180°. For the actual experiment. one of the detectors at 0 = 90° was not installed
because it would be mostly blocked by the magnet return yoke (discussed in § 3.5).
On the side without the magnet yoke, the 90° detector was installed. Another SeGA
detector normally included in this configuration, the d) = 180° detector, could not

be installed due to space constraints caused by the magnet support. frame. As was

67

 

Figure 3.7: Drawing of the full SeGA configuration designed for the TF g-factor
measurements.

previously discussed. the <15 = 180° detector would have zero precession sensitivity,
and was left out. In total. 14 detectors were used for the g-factor measurements:
six pairs at 9 = 29°,40°.60°, 139°, 147°. 151°. and single detectors at 0 = 24° and
0 = 90°.

Prior to the experiment. the detector positions were measured with a theodolite
system to provide accurate position information for Doppler correction. The measure-
ments determined the relative coordinates between the center of the stand (target
position) and several reference points on the outside of the SeGA crystal cryostats.
Detector positions were measured — instead of simply using coordinates from a CAD
drawing — to account accurately for the flexing of the SeGA stand when the detec-
tors were mounted. The locations of each germanium crystal segment relative to the
outside of the crystal cryostats were previously catalogued for each SeGA detector by
Miller [105], and this information was combined with the theodolite detector position
measurements to find the segment positions relative to the SeGA stand target posi-
tion. For more information on the theodolite measurement procedure and the SeGA

stand construction, see Ref. [104].

68

3.3.2 Maximum Scattering Angle and Alignment

A maximum scattering angle was chosen for the Coulex process, based on the re-
lation between the impact parameter and center-of-mass scattering angle given in
Eqn. 2.25, and transformed into the lab frame using Eqns. 2.28-2.30. The resulting
lab angle was 9mm = 55°, corresponding to an impact parameter of touching spheres
plus ~1 fm at the entrance of the gold layer. The plan was to enhance the Coulomb
excitation yield above that of a smaller angle with a greater separation between pro-
jectile and target (of course, soon after entering the gold layer, the projectiles lose
energy and the projectile-target separation distance increases for a fixed scattering
angle). A similar approach has been used in other transient field measurements uti-
lizing Coulomb excitation, although at much lower energies (see, e.g., Ref. [106]).
With the theory of Coulomb excitation, implemented in the program GKINT [107],
the statistical tensors were calculated for 381408 scattering on gold to the maximum
angle of 55°, taking into account the slowing-down process in the gold layer which
changes the adiabacity parameter and therefore the statistical tensors. Cross-section
weighted statistical tensors for decays occurring after the transient field layer were
obtained for the two isotopes, with 388 having p20 = 0.3016 and p40 = 0.0326, and
4OS having p20 = 0.3647 and p40 = 0.0435, indicating an effective alignment of 56%
for 388 and 68% for 408, for the decays occurring in vacuum beyond the target. More
details of the GKINT calculations as well as the measured angular distributions will

be presented along with the data analysis in § 4.2.3.

3.3.3 Slope of Angular Distribution and Sensitivity

The sensitivity of the SeGA detectors to the spin precession of the 7-ray angular dist-
ribution depends on the logarithmic slope S (19) = [IV-dag: of the angular distribution
at the SeGA angles, and the number of '7 rays counted at those angles. Differentiation

of the angular distribution with respect to 0 yields, for the detectors located in the

69

plane perpendicular to the magnetic field,

8(0) = — (a2P21(cos 0) + (141341 (cos 9))/W(0), (3.1)

where P21 and P41 are the associated Legendre polynomials. Since the precession
occurs about the magnetic field axis, the spherical polar coordinate unit vector does
not correspond to the direction of rotation for detectors located out of the plane
perpendicular to the magnetic field. Instead of 3/60, the directional derivative of
the angular distribution W(6l) should be taken about the magnetic field axis, in
the direction of the angular coordinate in a cylindrical coordinate system (where
the cylinder axis coincides with the magnetic field axis). The logarithmic slope for
a detector located at coordinates (9,05) is S(l9.¢) = 5(0) sin (I), with 8(0) given by
Eqn. 3.1. In the design of the SeGA stand, care was taken to place detectors as close
as physically possible to the beam axis (corresponding to the smallest possible 0),
while keeping the crystal axes tangent to spheres centered on the target (to take
advantage of the angular sensitivity of the segmentation). The detectors placed at
0 = 29° and 9 = 151° (the two pairs shown in Fig. 3.6) represent the best physically-
achievable detector placement consistent with a 4” beampipe. The four detectors at
0 = :l:29° and j: 151° dominate the total sensitivity of the array, since they are located
near the angles of maximum S (9) and sin 05 = 1 in the plane perpendicular to the

magnetic field.

3.4 Phoswich Detector

A beam particle detector was designed for use in detecting the 38308 ions exiting
the HVTF target. The selected phoswich detector served as the trigger to look for
particle-7 coincidences between scattered particles near the beam axis and '7 rays

de—exciting the nuclear level under investigation.

70

3.4.1 Design and Construction

There were several parameters used to decide on the design of the phoswich detector.
Of prime importance was the ability of the detector to handle high event rates. Since
the detector is situated over an angular range from 0° to some angle 0m”, it has to
detect a very high fraction of the incident beam, on the order of 105 particles per
second in the case of 388. Such a high rate prohibits the use of solid—state detectors
which would succumb to radiation damage. Plastic scintillators have very fast fluo—
rescent times of order nanoseconds, and can survive particle rates of 105/s. Another
important characteristic for the particle detector is the ability to identify particles.
The most common technique for particle identification is the AE-E method, where
two detectors are used to measure an energy loss (AE) and a residual energy (E).
According to the Bethe—Bloch formula, charged particles generally deposit an amount
of energy that is proportional to their atomic number Z 2. The residual energy left
after passing through the energy loss detector is collected by a second detector, which
is thick enough to stop the particles. The total energy is proportional to the particle
mass. To employ a AE-E method with plastic scintillators, it is convenient to use a
“phoswich” design. The term phoswich comes from the phrase “phosphor sandwich,”
meaning multiple layers of scintillating materials. Using two different scintillators that
have very different light emission characteristics, the light emitted as charged parti-
cles pass through one layer and stop in the next layer can be detected with a single
photomultiplier tube and separated using pulse-shape analysis.

The basic phoswich design chosen was a fast / slow plastic scintillator design sim-
ilar to that previously employed for intermediate-energy Coulomb excitation at the
NSCL [10, 108]. A 0.75 mm layer of fast BC—400 plastic scintillator and a 2” thick
layer of slow BC-444 plastic scintillator were supplied by Bicron Corporation. The or-
ganic scintillators are polyvinyltoluene—based materials with densities of 1.032 g/cm3,
primarily consisting of hydrogen and carbon. The two scintillators have very different

time characteristics. BC-400 has a fast decay time of 2.4 ns and BC-444 has a decay

71

 

 

 

 

 

 

,, .
5 1. I
E -, . 80444
a I. . o
e I 304 i
in Be_am_§
A r": g
Total Light
Time

Figure 3.8: Components of phoswich light from the different types of plastic scintil-
lator. The light from the fast BC-400 is shown as the dashed line, the slower light
from the BC—444 layer is shown as the dotted line, and the solid line is their total.
Pulse-shape analysis can be used to discriminate between particles that penetrate the
fast plastic layer or stop in it.

time of 180 ns [109]. When a charged particle traverses both layers of the phoswich,
light is emitted by both layers with two characteristic time profiles, shown in Fig. 3.8.
Analyzing different time regions of the pulse reveals the amounts of “fast” and “slow”
light, providing a measurement of AE-E. Since the particles being detected by the
phoswich are heavy ions of Z ~ 16 and kinetic energy E ~ 150 MeV, they are weakly-
penetrating particles. The consequence is that the ion range in plastic is only around
40 pm. In order to perform a AE-E measurement, a fast plastic layer of thickness
< 40pm would be required. It was decided instead to use a “thick” fast plastic layer
of 0.75 mm, which is a macroscopic thickness that is rugged enough to be handled
and that will stop all the heavy ions before they reach the slow plastic. If nuclear
reactions occur inside the target and result in lower-Z ions accompanied by protons
or alpha particles, then the slow scintillator will emit light due to the penetration of
such lower-Z particles through the fast plastic layer. The beam purity for the 38,408

was projected to be high, and a test run proved that, for 38S, the purity was > 99%.

72

 

9/

 

 

V

XCITl

 

 

 

 

 

Figure 3.9: Geometry affecting difference in arrival times for 7 rays emitted from
the target and 7 rays emitted from the phoswich. The phoswich position affects the
particle flight time from the target as well as the the 7 ray flight time from the
phoswich to the Ge detector.

This eased the requirement for the phoswich particle identification from unambiguous
(Z, A) identification to a discrimination between heavy and light ions.

An appropriate diameter for the plastic scintillators, which were disk-shaped,
needed to be determined. The role of the detector was to limit the maximum particle
scattering angle to 0mm, accomplished by placing the detector at a distance 6’ from

the target position such that.

Tphos
i 3

 

tan Omar : (3.2)

where Tphos is the phoswich radius. The dimension E was selected to place the detector
at a sufficient distance from the target to allow the possibility to distinguish between
7 rays originating from Coulomb excitation in the target and 7 rays originating from
the beam hitting the phoswich detector. The time response of the SeGA detectors is
about 25 ns [110], and the distance 6 needed to satisfy the condition that the arrival
times of 7 rays coming from the target and the phoswich had a separation of 25 ns,
as presented in Fig. 3.9. The SeGA detector at 29° was selected for the purpose of

the calculation since it is one of the detectors located closest to the phoswich. Beam

73

particles were assumed to leave the target with 230 MeV kinetic energy, corresponding
to a velocity of 3.404 cm/ns for a nucleus of mass number 38. The velocity is essentially
an upper limit for particles exiting the HVTF target. The propagation time of a 7 ray
between the target and the 29° detector was calculated. The time for the particle to
travel from the target position to the phoswich, and then for a 7 ray to travel from the
phoswich to the 29° detector was also calculated and the relationship was determined

for the difference in 7 ray arrival times:

 

x + 21.43) cm \/(11.88 cm)2 + (x cm)2 24.5 cm
3.404 cm/ns 30 cm/ns 30 cm/ns'

 

 

At, = ( (3.3)
The distance 6 = (:r. + 21.43) cm was calculated from this relationship to be ~80 cm
for At7 > 25 ns. Using Eqn. 3.2 with 6' = 80 cm and 0mm- = 55°. a radius of
7.7 cm results. This corresponds to a detector diameter of 6.07”. and for ease of
manufacturing a 6” diameter was chosen. The 80 cm distance between the particle
detector and the germanium array lowers the 7-ray flux that the germanium detectors
see from the stopped radioactive beam in the phoswich, compared to conventional
transient field arrangements that often employ small silicon particle detectors located
a few cm from the target.

After choosing the scintillator radial dimensions, the next step was to design
the light collection. system and choose a photomultiplier tube (PMT). A fast PMT
with high rate capability and good energy and time resolution is required for the
application of high-rate heavy ion particle detection, and the Electron Tubes, Ltd.
9821B phototube was selected. This PMT has 12 linear-focused BeCu dynodes and
a 3” (67 mm active) diameter. The linear-focused design generally has very good
timing characteristics and pulse linearity. From an efficiency standpoint, the larger
the PMT diameter the better, because from the conservation of photon flux per unit
solid area, the maximum fraction of light that can be collected is the ratio of the PMT

active area to the area of the scintillator cross section. The light collection efficiency,

74

K

 

r—-1
T 98213 PMT PMT Base
t
(u

. ' 80-808

6.98

 

 

 

 

 

 

 

 

Figure 3.10: Side view of the phoswich detector, including the scintillators, light guide,
photomultiplier tube, and voltage-divider base. Dimensions are in centimeters.

therefore, has a theoretical upper limit of the ratio of the squares of the radii of the
two components, or (r p AlT/rphos)2v which is ~0.19.

Light collection uniformity, i.e., the dependence of the light collection on particle
position, was a higher priority for the light collection system than total efficiency. A
non-uniform light collection could mean that particles of the same energy but a dif-
ferent implantation position will appear to have different energies. The large diameter
of the scintillators means that light emitted near the edge of the detector is proba-
bly less efficiently collected than light emitted directly on the axis of the cylindrical
scintillator. In order to maximize the light collection uniformity, a light guide was
designed that tapered from the 6” scintillator diameter to the PMT diameter over a
distance on the order of the scintillator diameter. The shallow angle of the light guide
enhances the collection of light emitted near the edge of the scintillator. Placing the
PMT 6” back from the scintillator also ensures that the PMT views the entire scin-
tillator volume in a more uniform way than it would if a short light guide was used.
The phoswich assembly is depicted in Fig. 3.10. The light guide was made of polished
acrylic and joined by the manufacturer to the BC-444 slow plastic scintillator with

BC-600 optical epoxy. The PMT was attached to the light guide with Tra—Bond F113

75

optical epoxy. The thin fast plastic layer was secured in place by aluminized mylar
foil. No optical coupling was used in between the fast and slow plastic layers (such
as optical grease) because any material would introduce a dead layer of inherently
non-uniform thickness.

The type of reflective wrapping to use for the phoswich was considered carefully,
since wrapping was needed to hold the fast. and slow pieces of scintillator together. It
is important to have a very thin and uniform layer of material on the front face of the
fast scintillator, to avoid attenuating the weakly-penetrating heavy ions. Aluminized
mylar was chosen since it is highly reflective and serves to reflect light from the
scintillators back towards the PMT. Three options were considered for the sides of
the scintillators and light guide: no wrapping, a specular reflecting material. or a
diffuse reflecting material. To use no wrapping is to rely on internal reflection alone
for the light collection. Internal reflection reaches 100% when the condition of total
internal reflection (TIR) is met. TIR occurs when light is incident on a boundary
between regions of different indices of refraction, and the angle of incidence is less
than a critical angle defined by

0,”, = sin—1 :71) (3.4)
where no is the refractive index of the light guide and n1 is the refractive index
of the surrounding medium, in this case vacuum. Light will “leak” out of the light
guide as the angle of incidence increases beyond 00,”, in accordance with the Fresnel
equations [111]. There will be position—dependent light collection because of the cor—
relation between the angle of incidence and the amount of reflected light that reaches
the PMT. A better solution is to wrap the phoswich with a reflective material to
try to catch light that leaks out of the light guide. A specular reflective material is
one that obeys the Reflection Law (where the angle of reflection equals the angle of

incidence), and is usually shiny in appearance. There is also a large degree of corre-

76

 

 

 

Figure 3.11: A light ray undergoing specular reflection in the light guide. Points along
the ray trajectory are labeled numerically in order of occurrence. This trajectory
results in the light ultimately being reflected away from the PMT because of the
angled light guide.

lation between the light that reaches the PMT and its initial emission position when
specular reflectors are employed. In fact, since the light guide gets narrower from the
scintillator to the PMT, some light can be reflected away from the PMT, as can be
seen in Fig. 3.11. To minimize the correlation between initial position and the fraction
of collected light, a diffuse reflector (white Teflon tape) was selected. Diffuse reflectors

obey Lambert’s Law, where the distribution of reflected light has the form

W(¢) 0( cos ()5, (3.5)

with ¢ denoting the reflection angle. Eqn. 3.5 demonstrates that the reflected light
is independent of the angle of incidence, and is peaked around 0° (normal to the
surface). In this case, a large amount of mixing occurs for the light as it is reflected
in the light guide and there is less correlation between the initial light position and

the amount that ultimately reaches the PMT.

3.4.2 Characterization of Phoswich Detector

The light output of the phoswich is the sum of two components, a fast and slow light

pulse. The relative magnitude of the fast and slow light is determined by how much

77

 

 

Signal

 

 

 

JA+ J

 

Time (100 ns/div)

Figure 3.12: Oscilloscope trace showing the two phoswich signals with a relative delay'
between them, and the common QDC gate. The gate intersects the signals in different
regions and the comparison of the two regions accomplishes the pulse-shape analysis.

ionization is caused by the charged particle, electron, 7 ray, or neutron as it trav-
els through the scintillator. In order to analyze the relative magnitudes of the fast
and slow light, pulse-shape discrimination was accomplished by charge-integrating
different parts of the phoswich pulse using a Charge-to—Digital Converter (QDC).
The QDC used was a CAEN Model V792 Charge-to Digital Converter, which has
32 channels that integrate the incoming charge falling within a selected time win-
dow. Pulseshape discrimination is possible using the V792 module by splitting the
phoswich signal into two and delaying one relative to the other. The full pulse shape is
integrated in one channel of the QDC and only the tail of the same pulse is integrated
in another channel. An oscilloscope trace of the two signals is shown in Fig. 3.12. The
actual scintillation light pattern is preserved during the amplification process in the
photomultiplier tube, but the electronics beyond the PMT anode has a particular
RC time constant as it responds to the current produced by the PMT, stretching the
signal in time.

The phoswich detector PMT was biased to -1700V using a Canberra High Voltage
Power Supply. The PMT base, which was designed and built at the NSCL, had two

outputs. One output is positive and originates from the first dynode of the PMT and

78

 

VME BUS

I f FULL

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

CAEN V262 NSCL GATE CAEN V792 510““
I/O MODULE GENERATOR QDC
‘ TAIL A CABLE
SIGNAL DELAY
CABLE P
DELAY ASSIVE
SPLITTER
A
TENNELEC NSCL FAST 3:3:
TC455 CFD ‘ INVERTING AMP
: ANODE
] SIGNAL
v I
TRIGGER DYN ODE "-‘ll
SUBSYSTEM SIGNAL PHOSWICH DETECTOR

 

 

 

Figure 3.13: Electronics used in the phoswich subsystem. Two signals originating from
the PMT base were used for timing and energy analysis. The positive signal from the
first dynode of the PMT was inverted and discriminated with the TC-455 CFD for
timing and triggering purposes. The signal after all 12 dynode stages was passively
split into two and sent to the charge-integrating ADC (QDC) with a relative delay

between them. Both energy signals were integrated with a common gate started by
the CAEN V262 I/ O Module.

79

 

   

Tail Signal

01 particles

 

 

 

Full Signal

Figure 3.14: Spectrum taken using a 241Am a source and a 6000 7-ray source. The
7 rays are primarily along the slow (tail signal) axis and the (1 particles are along the
fast (full signal) axis.

was used for timing. The other output is the negative anode signal which comes after
amplification by all 12 dynode stages. The positive dynode signal was sent through
an NSCL fast inverting shaping amplifier and was used as input for a Tennelec TC—
455 Constant Fraction Discriminator, which produced the timing signal. A rotary
attenuator and the fast amplifier gain could be adjusted to control the amplitude
of the input to the CFD. The negative anode signal was passively split into two
components, and a cable delay of 100 ns was used for one of the signals to allow pulse-
shape discrimination. A master live gate signal that was created by the logical AND of
the CFD signal and a computer NOT-BUSY signal opened a QDC integration gate for
300 us. The management of the master live gate was done using the NSCLDAQ [112]
readout software that communicated via the LEMO outputs of the VME V262 I/O
module.

Radioactive sources that emit different types of radiation were used to characterize
the phoswich detector. The goal was to ensure that the different types of radiation
could be distinguished by pulse-shape discrimination, and it was found that there
was a discernible difference between the radiation types. The a particles emitted by
typical 0 sources such as 228Th or 241Am generally have ~5 MeV of energy, with a.

maximum of ~8 MeV for one a emitted by 228Th. Such particles can not penetrate

80

the 0.75 mm fast plastic and therefore only give rise to fast scintillation. The 228Th
chain, however, also includes fi-decaying daughters as well as de—excitation 7 rays from
the resulting nuclei, and those radiations have greater penetrating power, leading to
slow scintillation (in the case of 7 rays, almost exclusively slow scintillation). The
sources used were 137Cs, which emits a single 7 ray of energy 662 keV; °°Co, which
emits two 7 rays of energies 1172 keV and 1333 keV; 90Sr, which emits 73 particles
with an endpoint of 546 keV, and subsequently the 90Y daughter emits [3 particles
with an endpoint energy of 2.3 MeV; 241Am, which emits a 5 MeV a accompanied by
a 60 keV 7 ray; and 228Th, which has a complicated decay chain including (7, fl, and
7 radiation. A representative spectrum from the 60Co and 241Am sources is shown in
Fig. 3.14.

Signals consisting of purely fast pulses should in principle only contribute to the
“full signal” and not to the “tail” signal, whereas signals consisting of only slow pulses
should be identical in both the “full” and “tail” integrated regions. In practice, there is
always some of each type present in each integrated region, which causes a departure
from a 45° angle between the “full” and “tail” signals when they are plotted against.

each other.

3.4.3 High Count Rate Behavior

A test was undertaken to determine the effect of count rate on the phoswich perfor-
mance. The test was performed at two PMT bias voltages. -1700V and ~2000V. A
241Am alpha source with an activity of 100 uCi was used for the test. The 0 particles
from the decay of 241Am have energies of z5.6 MeV and either directly populate
the ground state of 237Np or an excited state which emits a 60 keV 7 ray. Neither
a particle can penetrate the 0.75 mm BC-400 layer and therefore only cause fast
scintillation. The resolution of the plastic scintillator is not sufficient to separate the
two a particles that are separated by 60 keV.

The test was conducted in vacuum to eliminate attenuation of the 0 particles on

81

 

 

 

 

 

 

 

 

250 . 900 .
8 PMT Bias -1700v PMT Bias 200W
5 225 800 b
.56: II.- .. ]" .
g I I - 700 L .
5 2m . l l . I
5% " 500 -
Q.

175 ‘ 1 500 . .

0.0 1.5 3.0 4.5 0.0 1.5 3.0 4.5

Master Gate Rate [105 counts /s]

Master Gate Rate [105 counts /s]

Figure 3.15: Centroid of 241Am (1 particles as a function of a particle rate resulting
from changing the source-detector distance, for a PMT bias of -1700 V (left) and
-2000 V (right). Note the change in vertical scales.

the path from the source to the detector. The 100 ”Ci 2“Am source was fixed to
the blanked-off end of a 10” beampipe and the phoswich detector was moved in 17
steps from 2 to 74 cm from the source, resulting in count rates from 420.000 to 5,000
master gates per second.

At each source distance and PMT bias voltage setting. charge-integrated spectra
were recorded and a two-dimensional figure was plotted showing the full pulse spec-
trum versus the tail spectrum. The characteristics of the full pulse spectrum, which is
most sensitive to the fast plastic scintillator. was the main observable of the test. The
peak centroids were determined and the behavior with -1700V PMT bias is shown
in Fig. 3.15. The peak centroids shift as a function of count rate, with a 5% shift
observed at 150,000 counts per second and a 10% shift observed at 400.000 counts
per second. The PMT behavior at -2000V was similar to that at -1700V. with a 5%
gain shift at 150,000 counts per second and a 10% shift at 400,000 counts per second.

An appreciable deviation of gain was observed with rates in excess of 105 a par-
ticles per second, and since the sulfur projectiles in the g-factor measurement. stop in
the fast plastic layer, the same qualitative rate-dependent behavior would be expected

when sulfur ions are implanted into the phoswich. At the highest rate of 420.000 par-

82

ticles per second, the a peak quality was visibly diminished as the peak shifted down
towards the lower level threshold. To guarantee undiminished performance of the
phoswich detector at the 105/s rate expected for the 388 fragments, a Pb mask was
designed to block the angular range 0°-2.5°, where the Rutherford scattering cross
section is very high but the Coulomb excitation probability is low due to the large
impact parameter. No mask was used in the 408 runs since the particle rate was an

order of magnitude lower and pileup and gain shift effects were not appreciable.

3.5 Compact Electromagnet

An external magnetic field that polarizes the ferromagnetic target layer was required
to utilize the transient field. It is reasonable to expect that the transient field effect
scales with the magnetization of the ferromagnetic host material (the magnetization
dependence is explicitly written in the Rutgers parametrization, Eqn. 2.10). Since
ferromagnetic materials exhibit a saturation magnetization for a particular applied
magnetic field strength Bsat: nothing is gained by applying a stronger magnetic field
than Bsat- For polycrystalline iron, the magnetic field strength required to fully align
all magnetic domains is on the order of 10’1 T. or ~1000 G [113], and can be lower de-
pending on the iron sample and foil geometry. A magnetic field strength of ~1000 G is
readily available with permanent magnets, bttt because the execution of a TF g-factor
measurement relies on reversing the magnetic field direction periodically in order to
increase the measured effect. and cancel most systematic errors, an electromagnet is
usually used. This section will describe the parameters of the electromagnet solution

designed for TF g-factor measurements at the NSCL.

3.5.1 Design and Construction

Several physical constraints were included in the electromagnet design. The magnet

and target chamber needed to fit inside the SeGA stand. which holds the SeGA

83

detectors with their crystal centers 24.5 cm away from the target position. Taking
into account the crystal cryostat diameters, the magnet and target chamber were
constrained into a cylindrical space with a diameter of ~40 cm. The target foils
had dimensions 26 x 26 mm2, chosen to accommodate the projected beamspot size
realized for experiments run in the N3 vault. The width of the magnet poles was
designed to be no smaller than the target width to ensure full target foil saturation.
The target foils were held in a target frame designed to secure them with a minimal
amount of material on all sides. The magnet pole tips were machined with grooves to
accommodate the target frame, and therefore the magnet “pole gap” was fixed at 32
mm.

Besides the physical dimensions imposed by the target frame and the SeGA stand,
it was important to design the magnet such that it produced a magnetic field sufficient
to magnetize the target foil with a minimum stray fringing field. A strongly—varying
magnetic field can damage Ge detectors, which can experience a Penning discharge,
or a spark at the surface of the germanium crystal which is held at several kV during
normal operation. It is unlikely that. a Penning discharge will occur in a magnetic field
of only 1000 G. A strong magnetic field may also have an effect on the beam itself
as it passes through the magnet fringe field. The steering of the beam, called beam
bending, is usually a concern for low-energy beams and not for intermediate—energy
beams. For the purpose of the magnet design, the fringing field was minimized by
keeping the pole dimension short along the beam axis direction, and by choosing a
“C”-magnet design that uses a return yoke to channel the magnetic flux in a well-
defined region. The pole tips were also tapered from 1” thick along the beam axis to
0.25” thick where they held the target frame, to increase the magnetic flux density
near the target foil. The saturation magnetic field could then be attained without
introducing a large region of magnetic field along the beam direction for the beam to
pass through, which would be the case if the pole tips were symmetric in dimension

along the beam axis and perpendicular to it.

84

Table 3.2: Results of iron foil magnetization measurements [114].

 

Applied Field [Gauss] Magnetization [%]

 

367 91
487 97
621 99

 

 

 

 

To establish the magnetic field required to saturate the iron foil, the magnetiza-
tion of a piece of the same iron used for the HVTF target was measured using an AC
magnetometer at Rutgers [115], with the results shown in Table 3.2. Those measure-
ments determined the magnetization of the foil as a function of applied field, and it
was concluded that the iron foil would be fully magnetized with an applied field of
greater than 700 G.

Simulations were run using the code POISSON [116] to model the magnet design
and to estimate the number of ampere—turns needed to produce a minimum magnetic
field of 700 G at the center of the gap. The calculation included the small gaps where
the magnet yoke inside and outside the vacuum was separated by the thinly-machined
regions of the aluminum target chamber flanges (see § 3.7). A plot of the magnetic flux
is shown in Fig. 3.16. The target value was determined to be 1200 ampere-turns. The
coils were to be operated by an in-house power supply, which could deliver 20 Amps
at 20 Volts. The choice was made to use a large-diameter copper wire, AWG #14
which has a resistance of 2.524 I? per 1000 ft, and run 4.3 A through 280 turns of
the 14 gauge wire to get the 1200 ampere-turns. The power dissipated per coil was
calculated to be ~20 watts. The cross sectional area of 280 turns of this wire was
calculated to be a little over 1” square, which would fit in the available space on the
top and bottom flanges of the target chamber.

The iron pole tips and return yoke pieces were fabricated by the N SCL machine
shop from high-quality low-carbon steel. The return yoke was designed to fit on the

outside of the 4” cross target chamber, and was kept close to the side of the chamber

85

 
 
     

[I [I l \\
- Mu:
//:I ‘. ’::,;t] m:[::[,:, \\\\\ \(

Figure 3.16: Graphical results of a POISSON calculation showing the behavior of a
C—magnet including a small air gap near the coil.

in order not to protrude into the space near the SeGA detectors. The return yoke
and pole tips were made in seven pieces. Two small vertical pieces attached to the
top and bottom flanges with screws, and they extended 1.5” in the vertical direction
to accommodate the magnet coils and coil support plates. Two horizontal sections of
the return yoke connected the coil sections to the vertical return yoke section, which
were all 1 in2. The target frame was held by the upper pole tip with a set screw and
rested in a groove in the lower pole piece. The target was installed by attaching it
to the upper pole piece before assembling the magnet onto the target chamber. Since
the return yoke was detachable, the target could be accessed by removing the top
flange of the target chamber after disconnecting the yoke at the outer vertical piece.

The magnet coils were made out of AWG #14 enameled copper conductor using
a coil winding station available at the NSCL. To keep the turns secure and to prevent
electrical shorts, epoxy (miller-stephenson Versamid 140 and EPON® Resin 815C)
was applied to the coils after each layer accumulated. In order for the coil to be
removable from the form after being epoxied, all surfaces of the form were covered by
3M Teflon tape, to which epoxy does not easily adhere. Four coils were wound, with
290, 280, 260, and 259 turns.

The coils were tested using a bench top Power Designs, Inc. Model TW347D

0-16 V, 5 A power supply. Current was driven through the coils to test their heat

86

dissipation and to measure their current-voltage relationships. For the 290—turn coil,
a voltage of 3.6 V was required to sustain a 5.55 A current after the power supply
and coil were run for one hour to allow the operational temperature to be reached,
for which the coil resistance would plateau. The temperature of the coil was deter-
mined to be 165°F using a RayTek MiniTemp MT gun-type infrared thermometer. To
ensure safe operation, temperature switches were attached to the coils with thermally-
conductive epoxy. The devices were model STO-150 switches that trip at 145-155°F,
and they were wired into the control system to switch off the current should the tem-
perature exceed that range, which was established by safety regulations. The vault
pressurized air supply was used to blow air across the coils, using a regulator to obtain
a pressure of 20 PSI through 1 / 4” diameter nylon tubing. During normal operation,

the coil temperatures were ~110°F with a current of 4.5 A (1.1 kG).

3.5.2 Hall Probe Measurements

An F.W. Bell Model 5080 Gaussmeter (a Hall probe) was used to measure the mag-
netic field strength as a function of magnet current. The Au and Fe target foils were
secured to the target frame and the magnet chamber was fully assembled. The Hall
probe was held with the sensitive area inserted in a notch near the pole tip, ensuring
that it was held in a reproducible location. The magnet current was cycled at discrete
values :l:I and the magnetic field was recorded. A plot of the measured .field strength
is shown in Fig. 3.17. A magnet current of 4.5 A, corresponding to a 1.1 kG field at
the pole tips, was chosen to ensure full target polarization during the run.

The magnetic field was also measured along the beam axis to characterize the
fringing field. For the measurement, a current of 5.5 A was run through the coils and
no target foil was in place. The Hall probe was insensitive to the magnetic field at a
distance of 30 cm from the target along the beam axis. The results of the measure-
ments are shown in Fig. 3.18. In contrast to the pole-gap field measurement shown

in Fig. 3.17, the measurement shown in Fig. 3.18 has a large systematic uncertainty

87

1.5 ......................

 

 

 

 

 

 

l DDU'
6 1.0 - EDD
54 [:1
V . Cl
2 0.5
.2 '
LL.
.2 0.0
3 .
go —0.5 - DD
2 DD
-l.0 - C]
. mg .
1:1
-15 ---..-.-.1- ...........
-6 —4 -2 0 2 4 6
Magnet Current (amps)

Figure 3.17: Plot of magnetic field (in kiloGauss) measured with a Hall probe for
various values of the current (in amps) running through the coils.

750

 

500 . .. ~

250] L 4
l
o ; 4.4 w «Iv '” 1 .

0 5 10 15 20 25 30 35 40
Distance along beamline (cm)

 

 

Magnetic Field (G)

 

 

 

Figure 3.18: Plot of magnetic field as a function of distance along the beamline,
measured from the target chamber entrance flange, without the target in place.

88

because the position and angle of the Hall probe could not be held rigidly in place, in
a reproducible manner. This measurement gives at most a qualitative understanding
of the magnetic field along the beam axis, since the presence of the iron target would
be expected to narrow this magnetic field distribution near the target position. How-
ever, the Hall probe test firmly established that the magnetic field is negligible at the

Ge detectors, at a ~25 cm radius from the target.

3.5.3 Power Supply Control

The external magnetic field direction needed to be reversed at preset intervals (on the
order of minutes) during the transient field g-factor measurements. in order to reduce
systematic errors. Since the power supply used to provide current to the magnet coils
was an NSCL—built four—quadrant power supply, it could output. a current that is
proportional to the control system input voltage, including the polarity. It was possible
to convey information from the N SCL EPICS control system into the datastream by
interfacing the magnet control system with the data acquisition (DAQ) system, and in
that way the data could be sorted by field direction, which was reversed periodically.

The power supply used was the NSCL 20/20 power supply (maximum output
20 Volts and 20 amps). The power supply is designed to provide a constant current
that corresponds to the “set value” that is proportional to the voltage applied from
the control system; the power supply will automatically adjust the voltage necessary
to run the desired current through the attached load. Large-diameter current leads
were used to connect to the magnet coil leads in order to avoid transmission losses
since the power supply was located ~5 m from the electromagnet. A VME module
containing an on-board computer running VxWorks was used to run a control script.
written for the purpose of controlling the electromagnet power supply as well as
supplying the status signals to be included in the datastream. The script. also had
several user-selectable options to control its functionality. The user interface to the

control script was provided through a Panellylate application that could run on a PC

89

computer. The magnet could be run in a “swapping” or a “non-swapping” mode,
whereby either the magnet current would swap after a predefined number of seconds
from I —> —I (where the current I is set by the user), or the magnet current could

simply be held constant at the value I.

3.6 Data Acquisition System

The DAQ system consisted of separate subsystems for the germanium array and the
phoswich detector, which were developed and tested separately and then integrated
together. The model for a DAQ system at the NSCL is to run nuclear electronics
modules, standardized according to the N 1M, CAMAC, or VME standards, and access
the data by mapping the CAMAC addresses into the VME address space managed by
a DAQ computer connected to the VME crate. In this instance, the SeGA electronics
were CAMAC modules and the phoswich and PIN electronics were VME modules.
Trigger logic was performed using NIM modules, and the CAEN VME I/O module
shown as part of the phoswich subsystem in Fig. 3.13 was the master trigger manager
for the DAQ system that interfaced between the VME bus and the DAQ computer

that hosted the DAQ software. An overview of the DAQ electronics will be given here.

3.6.1 Germanium Detector Electronics

The SeGA subsystem consisted of electronics channels to process signals from the
central contacts as well as the detector segments, meaning that each detector has 33
channels of energy information. Each SeGA detector has a dual output of the central
contact energy, one of which is used for timing. The timing signals from each Ge
detector were processed by TC-455 NIM Constant Fraction Discriminator modules,
the outputs of which are combined using a logical OR gate in a NIM multiplicity logic
unit. The GE-OR signal was used for the trigger logic which will be described in the

next section. The Ge CFD signals were also fanned out to the SeGA Phillips 7186

90

 

CAMAC

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I I r if
PHILLIPS LECROY 2551 ORTEC AD413 PHILLIPS
718611 TDC SCALER 8K ADC 716411 ADC

CAEN C469 ORTEC 572
DELAY AMPLIFIER
TENNELEC ORTEC 863 CENTRAL
TC455 CFD _‘ TFA CONTACT
.. SeGA :‘——'_i_j;‘_: NSCL
TRIGGER DETECTOR i 1": SHAPERS
SUBSYSTEM SE GMENTS

 

 

 

Figure 3.19: Electronics used in the SeGA subsystem. One SeGA crystal produces a
full-energy central contact signal as well as 32 segment energy signals. The central
contact signal is split internally, and used for timing and energy analysis. The segment
signals are shaped and digitized independently of the central contact signal. The
trigger subsystem starts the TDC upon registering a valid trigger condition, and
individual SeGA detectors can then stop the TDC. The DAQ system determines
which SeGA detectors have above-threshold events by first reading out the TDC
module, and then reading out the central contact and segment ADC modules for the
SeGA detectors that have sent stop signals to the TDC.

Time-tO—Digital Converter (TDC) module that was started by the master trigger. The
TDC spectra were 800 ns wide and were used to discriminate the germanium spectra
in time, to correct for random coincidences. The central contact energy signals for each
SeGA detector that registered a 7-ray hit were digitized using Ortec 413 CAMAC 8k
Analog-to—Digital Converters (ADCs). The segment energy signals were shaped using
NSCL-built shaping amplifiers and digitized with Philips 7164 CAMAC ADCs.

91

3.6.2 Phoswich and PIN Detector Electronics

The phoswich detector electronics were presented in § 3.4.2, and the same electronics
arrangement shown in Fig. 3.13 was used during the TF g-factor measurements. The
trigger condition was a coincidence between the phoswich detector and the Ge array,
in order for the SeGA electronics only to be read out when both a 7 ray and an
associated beam particle were present. The phoswich subsystem provided the particle
trigger signal that was one half Of the coincidence logic shown in Fig. 3.20. The dynode
signal Of the phoswich was inverted by a fast amp and then discriminated using a TC-
455 CF D. The CF D output Signal opened a 300 ns gate called PARTICLE with which
an OR Of the SeGA CFD signals called GE—OR could satisfy a logical AND condition.
The result of the AND gate, GE—OR- and-PARTICLE, was a 340 as wide gate that could
overlap in time with a 40 ns—wide DELAYED—PARTICLE gate that occurred at the end
of the 300 ns PARTICLE gate. The logical AND of the DELAYED—PARTICLE gate and
the GE—OR- and-PARTICLE was the master trigger for the DAQ system. The timing Of
the trigger was with respect to the particle, regardless of the timing behavior Of the
germanium detectors. which was still preserved in the TDC spectrum for each SeGA
detector.

During the beam tuning, a 488 ,um silicon PIN detector was used to diagnose the
secondary fragments using a AE—TOF method. The PIN was located in the N3 vault
near the target position, as will be described in § 3.7, and was retracted from the beam
line during the TF measurements. During the beam tuning time, the DAQ trigger
was changed to PIN singles OR downscaled phoswich singles, and the PIN signal was
shaped and amplified by an NSCL Quad Shaper and digitized using a CAEN V785
VME Peak-Sensing ADC.

92

phoswich CFD

starts PARTICLE gate U

I
PARTlC LE

 

I ,
_|<———— 300ns——>|

 

I
' I
I
DELAYED-PARTICLE ; i r—

GE-OR

 

 

 

 

 

 

 

 

 

GE-OR and PARTICLE

 

 

I-
I
t-

 

 

| 300+40 ns f T"

1

 

 

 

 

 

7J—

(GE-OR and PARTICLE)
and DELAYED-PARTICLE
is master gate

Figure 3.20: Trigger logic used in the experiment. The phoswich CFD signal started
a 300 ns gate (PARTICLE) and a 40 ns gate that started at the end of the particle
gate called DELAYED—PARTICLE. The logical OR of the SeGA CFD signals (GE—OR)
started a 340 ns gate (CE—OR-and-PARTICLE) if it overlapped with the PARTICLE
gate, and the logical AND Of this gate with the DELAYED-PARTICLE gate created a
master trigger that was correlated with the phoswich timing, avoiding any jitter in
the germanium timing.

93

3.6.3 Magnet Interface to Datastream

In order to embed the status of the magnetic field direction in the datastream for
analysis, the magnet control script running in the control system VME crate controlled
coils in the Modicon/PLC to send high or low logic signals (0 V and 24 V) to the data
acquisition system. A system Of three signals were used, two for the electromagnet
polarity and one for the swapping status. Prior to the polarity Swapping, the “swap”
signal was set to true, and it remained true for a small delay after swapping to ensure
the magnetic field of the electromagnet had completely reversed before the swap was
considered finished. The three signals needed to be interfaced into the DAQ system to
be recorded as bits in the datastream (1=true, 0=false), and a hardware interface was
built to electrically isolate the three signals to avoid noise from the control system and
to adapt the 24 V control signals to the CAEN V262 VME I/ O Module used in the
DAQ system. The interface module consisted Of mechanical relays that would Open

or close current sources with the proper impedances to produce NIM logic Signals.

3.6.4 Data Acquisition Software

The software used for data acquisition was the NSCLDAQ Production Readout soft-
ware, which provides a framework for the integration of distinct data acquisition
subsystems into one coherent DAQ system. The SeGA electronics and data acqui-
sition code had been developed for previous experiments with SeGA coupled to the
S800 spectrograph, in the framework Of Classic Readout, an Older version Of the DAQ
software. The readout system for the phoswich detector was written as an event seg-
ment within Production Readout, and the Older SeGA Readout software was decou-
pled from the S800 and formed into another event segment for Production Readout,
using a functionality that encapsulates Older software to make it compatible with
Production Readout. The two separate subsystems (phoswich and SeGA) were then

combined, along with the magnet status flags, into one datastream structure.

94

LN2 dewar —. . - — filling system control

phoswich \ filling manifold
' /— 92” chamber endcap

 

 

 

 

 

 

 

 

 

 

 

10” pot
> turbo pump
"1.9. / [— quadrupole
\ ,,,,,,,,, way,»
A exhaust fan . ‘ 4‘
., SeGA -
h.
I I s... .....
‘ leads
. PLC
NAQ electronics magnet power cabinet
DAQ computer supply L

 

 

 

Figure 3.21: Arrangement of equipment in the N3 vault. The existing control system
interface via the PLC cabinet was used to control the magnet power supply and
interface with the DAQ electronics. The building liquid nitrogen supply was used
to fill the Ge detectors after first accumulating in a large dewar. Data acquisition
was run using electronics inside the vault, and data were recorded onto disks via an
ethernet connection.

3.7 Layout in N3 Vault

The layout in N3 is shown in Fig. 3.21, which shows the components used in the
experiment. The beamline was aligned using an Optical reticule mounted on the en-
trance flange inside the 92” chamber, as well as a target permanently fixed to the wall
inside the N3 vault beamdump to define the beam axis. An in-place 92” chamber me-
chanical roughing pump was used tO back a 900 liter/s Alcatel 5900 turbomolecular
pump that was installed on the bottom of the 92” chamber entrance pot to provide
pumping for the entire experimental beamline beyond a gate valve located after the
quadrupole doublet inside the N3 vault. The system could provide vacuum better
than 10‘5 torr, necessary to ensure the slowed projectiles did not scatter on their
80 cm flight path between the target and phoswich.

A 10” pot was installed upstream Of the target position. The pot contained a

5”-thick copper beamstop on which a phosphorescent paint was applied. The beam-

95

stop could be remotely raised to stop the beam, and a CCD camera could view the
beamspot for tuning. The Si PIN diode was installed in the pot and raised into the
beam for tuning, and to measure the transmission from the A1900 to the N3 vault.
In addition, the PIN could be used to estimate the transmission through the HVTF
target, by comparing counts in the PIN to counts in the phoswich with the PIN low-
ered out Of the beam. At the entrance to the pot, a collimator was installed in order
to prevent a beam halo from hitting the target frame or electromagnet pole tips. The
collimator consisted of a thick plate Of tungsten with a 1” diameter hole drilled in the
center.

Downstream Of the 10” pot, a section Of beampipe containing a periscopic mirror
was used to view the target position. A phosphor-coated aluminum plate on a circular
target frame was inserted in the horizontal section of the target chamber upstream of
the target position, to serve as a target viewer during beam tuning. Since the vacuum
needed to be broken to insert the viewer and also to remove the viewer after tuning,
it consumed some time during the experiment, and was only used once for the 38S
run and once for the 40S run.

The SeGA stand was clamped to a translation stage that was bolted to the vault
floor just downstream Of the 92” chamber rails. Two Optical targets were held in place
on the SeGA stand using aluminum alignment bars that were built for that purpose.
The same alignment bars were used when the stand was first completed to hold the
targets for measurements that located the SeGA detectors with respect to the targets
by means Of a theodolite system (this formed the database Of detector coordinates
used for the Doppler correction, which is discussed in §4.2.2). The translation stage
could move the entire stand (with respect to the beam axis) up, down, left, right,
forward, and backward, but it could not rotate about the vertical axis, the beam
axis, or the axis perpendicular to the two. In order to rotate the stand about the
vertical, the clamps holding it to the translation stage were loosened and the stand

was manually rotated. Rotations about other axes could be accomplished by inserting

96

Shims beneath the stand in order to level it, but by measuring with a mechanical level
this was seen to be unnecessary because the stand was level.

A 4” beampipe cross design was chosen for the target chamber. The 4” cross was
purchased and the necessary modifications were made by the NSCL machine shop.
The target chamber is shown in Fig. 3.22. The target chamber downstream pipe was
stepped to a 6” diameter as soon as possible after exiting the array, in order to provide
an Opening angle of ~9.5° for scattered beam particles. Below this angle, essentially
all of the Rutherford cross section is contained, ensuring that a negligible fraction
(~ 10’”) of the beam particles would stop near the array. The 6” diameter section is
welded to an adapter flange for connecting to an 8” beampipe that holds the phoswich
detector downstream of the array. The vertical section Of the 4” cross was used to
hold the electromagnet. The chamber was designed to allow the electromagnet coils
and return yoke to be located outside of the vacuum. In order for the magnet pole
tips to enter the vacuum Of the target chamber, a “magnetic feedthrough” rather

than a “physical feedthrough” was used. The top and bottom target chamber flanges

L.

 

 

 

 

 

 

 

 

 

.1

 

Figure 3.22: Cross section Of the HVTF target chamber, with the beam entering from
the left. The target is held in a target frame between the pole tips. The magnet yoke
is separated into a piece that attaches inside the vacuum chamber and a piece that
attaches to the outside of the chamber. The magnet coils are located outside Of the
vacuum.

97

were machined very thin (on both the outside and the vacuum side) in a 1 in2 region
to match the yoke area. On the outside Of the flange, the vertical magnet yoke piece
was attached with screws, and the magnet coil fit snugly around it, held in place by
the horizontal yoke piece. On the vacuum side of the chamber flanges, the magnet
pole tips were attached with screws. In this manner, the magnetic field carried by the
iron yoke could “jump” across the thinly machined section Of the flange without the
need for an actual vacuum feedthrough. The cross section Of the target chamber, in
Fig. 3.22, shows the internal arrangement of the target chamber.

An 8” beampipe, attached to the exit Of the target chamber, was the final section
of beampipe and housed the phoswich detector. The phoswich signals were brought.

out of the vacuum using a feedthrough plate that. was insulated from the beamline.

98

Chapter 4

Results and Conclusions

A separate run was performed for each sulfur isotope. For each secondary beam,
phoswich-7 coincidences as well as downscaled phoswich events were recorded, during
which the electromagnet ran at 4.5 amps (1.1 kG near the pole tips), switching polarity
every 600 3. Data were collected for 81 hours for the 38S fragments and 68 hours for
the 40S fragments. This Chapter presents the data analysis and results Of the two

measurements.

4.1 Particle Identification

Secondary beam particle identification was performed using the A1900 focal plane
detectors for energy loss and time-Of-flight information, verifying that the dominant
component of the beam was the desired sulfur isotope. When the ions were sent to the
experimental endstation, 7 rays from the radioactive decay Of the projectiles (which
stop in the phoswich detector) were used to confirm the particle ID. Known 7 rays
in 38C] and 38Ar were Observed during the 38S run and 7 rays in 40C1 and 40Ar were
Observed during the 40S run. After implanting beam into the phoswich, the beamstop
was inserted and the phoswich CFD threshold was raised in order to cut out triggers

coming from radioactive decay of the beam particles. Those signals were much lower

99

 

Tail Signal (a.u.)
“c

 

 

 

 

Full Signal (a.u.)

Figure 4.1: Phoswich tail signal vs. the full energy signal, with counts in a log scale
indicated by grayscale shade. A lower-level threshold is used in this plot for better
contrast. Most counts are for particles that only cause fast scintillation (the darker
band with low values of tail light). Particles that punch through the fast. plastic layer
will show up away from the dark band with larger values Of tail light.

in energy than the beam particles, but it was desirable to remove them as a source

Of random particle-7 coincidences.

4.2 7-ray Sorting

Particle-7 coincidences were recorded, and both lab—frame and Doppler-corrected
spectra were Obtained. Particle energy, time of flight, and SeGA TDC gates were
used to produce true coincidence spectra that were corrected for random particle-7
coincidences, and Doppler correction was Optimized for the 7 rays emitted after the
projectiles traversed the entire ferromagnet layer. The angular distribution of 7 rays
was analyzed and compared to the calculated angular distribution using the theory
Of Coulomb excitation, and the count rate changes upon magnetic field reversal were
obtained for each of the six pairs of detectors at symmetric angles 21:9 with respect
to the beam direction. Calculations were performed to determine the transient field
strength BTF: and g(2+) values were Obtained for 38408. The details of each step

will be presented in this section.

100

 

 

 

Time of F light (a.u)
<

 

 

 

Phoswich Energy (a.u.)

Figure 4.2: Time Of flight between the cyclotron RF and the phoswich, versus phoswich
energy, with counts shown in a log scale indicated by grayscale shade. A lower-
level threshold is used in this plot for better contrast. The three main features are:
(a) events triggering the slow plastic (light particles), (b) events having particle-7
coincidences primarily in the fast plastic only, and (c) downscaled particle Singles (no
particle-7 coincidence).

4.2.1 Phoswich, Time Of Flight, and 7 Time Cuts

A fast/Slow 2D spectrum. shown in Fig. 4.1. was Obtained by plotting the tail-
integrated QDC channel versus the total-pulse integrated QDC channel. The dark
band located along the “full signal” axis, corresponding to small “tail signal” values,
are the sulfur projectiles. Any deviation from the heavy-ion band is caused by an
increased proportion Of slow scintillation light, indicating punch-through events that
are caused by low-Z particles. In addition to the 2D fast/Slow spectrum, the phoswich .
energy signal (the full pulse shape integration) was plotted versus the time-Of-flight
between the cyclotron RF and the phoswich pulse, to form a 2D E—TOF spectrum
for the particles. Three distinct features are evident in the spectrum (a representa-
tive spectrum is shown in Fig. 4.2). By gating on different regions of the phoswich
2D fast/slow spectrum, it was established that the upper feature (labeled “a” in
Fig. 4.2) corresponds tO those events causing a large amount of slow scintillation in
the phoswich, meaning that. they are not heavy ions (those would stop in the thin
fast plastic layer) but are light particles coming from nuclear reactions in the target

or other material such as the Image 2 wedge. The middle region (labeled “b”) cor-

101

responds to phoswich events triggered by sulfur projectiles that are in coincidence
with 7 rays, and the lower region (labeled “c”) is made up of downscaled phoswich
events (the downscaler module introduces an additional delay that shifts this group
Of events away from the phoswich-7 coincidences).

The germanium spectra were produced by gating on the sulfur region of the
phoswich 2D fast / slow spectrum and the particle-7 coincidence region Of the phoswich—
TOF 2D spectrum. The 7-ray spectra were corrected for random coincidences, be-
ginning by projecting Out the TDC spectra corresponding to particles in the 2D
fast / slow and TOP spectra that came in coincidence with 7 rays in the energy range
1000—1500 keV for 38S runs and 700-1200 keV for 40S runs. The 38S TDC spectra
are shown in Fig. 4.3, which includes the raw (ungated) TDC spectra superimposed
with the phoswich- and energy-gated TDC spectra. For each SeGA detector, gates
were placed around the prompt 7-ray peak in the TDC Spectrum (shown as the filled
channels in Fig. 4.3), and gates were placed to each side Of the prompt peak, in order
tO produce germanium spectra with good coincidences and with random coincidences.
Normalized random coincidence Spectra were then subtracted from the coincidence
spectra for the prompt region, yielding SeGA spectra that were analyzed to extract
the count differences at each angle for magnetic field up and magnetic field down.
Random—subtracted spectra were created for both the lab-frame and the Doppler-
corrected spectra. Separate spectra were created for magnetic field up and magnetic
field down data; these spectra could be added together to Obtain an unperturbed

angular distribution since the precession angle is small.

4.2.2 Principles of Doppler Correction

One of the principal challenges posed by in-beam 7-1'ay studies is the Doppler shift
Of the 7 rays. For intermediate-energy beams, the Doppler shift is very pronounced.

For a projectile with velocity 73 = 1.1/c. a 7 ray emitted at an angle 9 with respect to

102

 

7.5x104

 

 

 

 

 

 

 

 

 

 

 

  

 

 

 

  

 

 

 

7.5x 104

 

 

 

 

 

 

 

 

 

 

 

 

  

 

 

 

 

1.0x104 1.0on4

1.0x103 1.011103

1.0x102 1.011102 ..

I.0x10' 1.0x10' 'i

7.5on4 7.5x104

1.0on4 1.011104 .

1.071103 1.0x103

1.0on2 I.0x102

1.0x10' 1.0xto'

7.51104 9:90,, 7.5x104

1.0x104 1.0x104

I.0on3 1.0xI03

1.01102 1.011102

1.0x10' ' ' ‘ ‘ * 1.0x10' ”
0 1000 2000 3000 4000

7.51104 931;”, 7.51104 e=+151°

1.0x10“ 1.07110“ ’

1.0on3 1.0x103

1.0x102 1.0x102

1.0x10' ' ‘ ~ ‘ 1.0x10' ‘ *
0I000200030004000 01000200030004000

Figure 4.3: SeGA time spectra for 38S showing the raw spectra (upper curves) and
the resulting spectra when gated on the 2D phoswich heavy ion region, the phoswich—
TOF good particle—7 coincidence region, and SeGA 7-ray energies between 1000 and
1600 keV (lower curves). The shaded region shows the channels used to make the
prompt 7-ray TDC gate.

103

the projectile direction will have an energy

E,“C\/1_——II_2

Elab =
7 1 — ficos9

 

(4.1)

where Efymc is the energr Of the 7 ray in the rest frame Of the projectile. Doppler
correction is performed event-by-event on the detected 7-ray energy Eff” in order to
recover the “true” 7-ray energy. From Eqn. 4.1 it can be seen that the beam velocity 9
and the 7-ray emission angle 9 are required. Since the projectiles only undergo small-
angle scattering in intermediate-energy Coulomb excitation, the 7-ray emission angle
can safely be measured with respect to the beam axis. In other words, the angle 9
used in Eqn. 4.1 is the 7-ray detection angle measured with SeGA.

SeGA detectors give sub—crystal position resolution for the 7-ray interaction point
in order to improve the Doppler correction over that of a conventional germanium
detector. For 7 rays with an energy Of ~ 1 MeV, the primary type Of interaction in
germanium is Compton scattering with a mean free path Of 3 cm as calculated from
the germanium mass attenuation coefficient [117]. Since the crystal size is larger than
the mean free path, it is possible that more than one interaction will occur within
the crystal, and hits might be recorded in more than one segment. A method to
assign the “first hit” segment and therefore choose the corresponding angle to use in
the Doppler correction formula must be implemented. The method used to choose
the first-hit segment for the SeGA detectors was the following empirically—justified
procedure: the segment with the highest energy signal is selected to be the first-hit if
the total 7-ray energy is above 500 keV, and for 7 rays below 500 keV total energy,
the lowest-energy segment is chosen.

In Compton scattering, a photon imparts some of its energy to an electron and
then continues traveling through the material (perhaps leaving the material if the
interaction occurs near a boundary edge). The electron and lower-enery photon

leave the scattering site with an angle between them that governs how much energy

104

 

 

 

 

 

 

 

 

120 E 1.5
R 2
3'3. 90 P ‘ >4
2 E” 1 *
O0
5 60 » . g
31 E l
as E 0.5 *
b 30 ~ I c:
> 00
< "s
O ‘ . . LE 0 . . .
0 0.5 1 I5 2 0 0.5 l 1.5 2
Initial gamma-ray energy [MeV] Initial gamma energy [MeV]

Figure 4.4: Klein—Nishina differential cross section for the Compton scattering Of
1 MeV 7 rays (left), and the final 7-ray energy as a function Of incoming 7—ray energy
for Compton scattering at the energy-dependent average scattering angle (right).

is transferred from the photon to the electron. The relationship between the Compton

scattering angle 90 and the outgoing photon energy can be written as

E — E”
f — 1+ (ET/7726262)“ —00890)’

 

(4.2)

where E23! are the initial and final photon energies and me is the electron mass.
The angle- and energy-dependent cross section for Compton scattering is called the

Klein-Nishina equation [118]:

2
‘10 ’"(21Ef Ei Ef ~2
_:__ _ __ 4.
d9 2 E? (Ef-l- Ei Sln 9C , ( 3)

where m = e2/mec2 is the classical radius of the electron, and 9C- is the usual angle
between the initial and final photon directions. Inserting 4.2 into 4.3, one Obtains
the scattering cross section as a function of 9, shown on the left side of Fig. 4.4, for
1 MeV 7 rays. For Compton scattering at the average 90 for a given incoming photon
energy E,, which can be calculated using Eqn. 4.3, the final photon energy is plotted
in Fig. 4.4. On average, an incoming 7 ray will deposit most Of its energy in the first

interaction, at least for 7 rays Of more than Z500 keV. Another interaction becomes

105

likely for a photon with an energy less than 500 keV as it travels through germanium,
and that is absorption through the photoelectric effect. Since a photoelectric absorp-
tion is the most likely event for the lower-energy 7 rays, resulting in the disappearance
of the 7 ray and all Of its energy deposited in the segment containing the interaction,
it is not as likely that the highest-energy segment results from Compton scattering.
Following the assignment of the first-hit segment, the angle 9 in Eqn. 4.1 can be
calculated based on the known coordinates Of the segment with respect to the target.
position. The detection angle as well as the beam velocity 6 are used in Eqn. 4.1 to
perform the Doppler correction. The angular width Of a segment, which is 1 cm wide
and located 21 cm from the target, is ~ 27°. The effect of this angular uncertainty A9
and any uncertainty in the beam velocity Afi on the width of the Doppler-corrected
7-ray peak can be seen by partial differentiation of Eqn. 4.1 with respect to 9 and I3,

yielding

AE, 2_ Esme 2 2
(7.) - f—I-Iacosol‘m

cos9—fi 2 ’2
+ ((1—132)(1—13c039))(AB) (4.4)

(A570,)?
+ — ,
E7

 

where a term accounting for the intrinsic energy resolution of the germanium has been
added. As was pointed out by Glasmacher [119], the relative width of the Doppler-
corrected peak scales directly with the uncertainties A9 and A6 and not with the
7—ray energy. For the case of the 38’40S measurements presented here, the particle
velocity ranges from 6 ~ 0.28 at the front of the target to 13 ~ 0.08 at the back Of the
target. The A6 Of 0.2 means that the A6 term in Eqn. 4.4 will be the dominating term
if a constant )6 value is used for the Doppler correction. In the usual intermediate-
energy Coulomb excitation studies, there is not a large change in velocity through the

target and the use of a constant )6 such as the mid—target velocity is acceptable. In the

106

 

 

 

 

 

 

 

 

 

 

0.05 0.020 . f
5 0
'g '04 ’ 0.015 I
'8‘ 0.03 ~
g 0.010
>. 0.02 .
E.“
E 0.01 > 0.005 f

0.00 ‘ 0.000 ‘ ‘ ‘

0 45 90 135 180 0 45 90 I35 I80

SeGA Angle [deg] SeGA Angle [deg]

Figure 4.5: Square of the fractional energy resolution using Eqn. 4.4 with A9 = 10°
for two different A6 values: A6 = 0.2 (left) and A6 = 0.1 (right). where the solid
line represents the term in Eqn. 4.4 involving A6, and the other two curves are the
A9 term and the total.

HVTF case, the impact on (AE / E )2 due to two values Of A6 is displayed in Fig. 4.5.
The panel on the left represents the actual case if a single beta is used. and the panel
on the right represents a factor of two decrease in A6 to get A6 = 0.1, where it is
apparent that the A9 and A6 terms are on equal footing. In the present context, the
long Doppler tail that results from the emission of 7 rays during the slowing—down
process is not the feature Of interest, as will be discussed further. Instead, the portion
of the spectrum resulting from 7-ray emission after the projectiles have left the target
is more important, and the Doppler correction was performed using a 6 determined
event-by-event, in order to Optimize the reconstruction of the in-vacuum flight peak.

Event-by-event Doppler correction was performed by using the SeGA first-hit
segment positions and the after-target particle velocities determined from a transfor-
mation of the phoswich energy signals into velocities. The event-by-event 6 values
resulted in an improved Doppler correction that was Optimized for the 7 rays emitted
after the projectile experienced the full transient field. The average after-target ve-
locities were extracted using the lab-frame centroids Of the Doppler-shifted projectile

2+ 7-ray peaks. The Doppler shift formula for 7 rays emitted in flight,
E0 = T'Elab (1 - I3 00859 (4.5)

107

 

3”S 1292 keV inverse centroids

4O - -
in Lab frame vs. c050 S 903 keV Inverse centrOIds

8,5e-04 1.11 Lab frame vs. COSO

 

 

    

     

 

 

 

 

 

 

' ' 1.25e—03 . .
310W *6-43 ' 6'05 Slope: -9.79le-05
.7 Mercer“: 77726-04 Intercept: I.I mm
3 8.0e-04 ' 1.15c-03 :
as.
s
g 7.5e-04 . 1.05e-03 »
7.0e-04 e e * 9.50c-04 . . .
-l -05 0 0.5 I -l -05 0 0.5 I

cos 0 cost)

Figure 4.6: Inverse peak centroids versus cos9 for 38S (left) and 40S (right) 7 rays
in the lab frame, with the slopes and intercepts Obtained by a linear fit tO the data
shown in the plots. The average after-target velocities are the ratios Of the slopes and
intercepts (see text).

was used to determine the average 6 Of the moving source. since the 7-ray energy
was measured at different angles. Solving Eqn. 4.5 for E10,, and inverting both sides

Of the equation gives

1 _ 7 — 76 cos9 7 _ 76cos9 (4 6)
E100 E0 En E0 ' '

It can be seen that a plot Of l/Elab versus cos9 will be a line of Slope —67/E0 with
an intercept at 7/E0; the quotient. Of those two quantities is the beam velocity 6.
Fig. 4.6 shows the results of a least-squares fit to the inverse centroids in the lab
frame versus cos9 for 38S and 40S.

The values of 6 extracted from the fits are 638 = 0.083(1) and 6‘") = 0.088(1),
where the quoted errors are 30% larger than those given by the errors Of the fit
parameters. The 638 and 640 values were used to fix the average velocity after the
transformation from phoswich energy channel to 6 value. The velocity distribution
was assumed to have a x/E functional form, and it remained to choose scaling factors
that would result in a realistic velocity distribution (i.e.. ensure that low phoswich

channels did not give 6 ——+ 0). From the phoswich energy-TOP plot. the kinematical

1 08

 

450

 

        
 

 

 

 

 

 

 

Lab Frame 56Fe 1238 keV 200 ‘ Lab Frame
300 - / 385 56Fe 846 keV
/ 100 v 40S
150 I /
O 0
v=0.08c Frame v=0.09c Frame
3 300 _ 561:6 388 100 _ 561:8 408
,3 / / /
C
a 150
C)
0 0

Projectile Frame

200

   
 
 
 

403 903 keV
/

56Fe
100

 

 

 

 

 

 

O A A A A l A A _A_ A O A A A A A A
1100 1350 1600 650 900 1150

Energy [keV] Energy [keV]

 

Figure 4.7: Laboratory frame (top), single-6 Doppler corrected (middle), and event-
by-event Doppler corrected spectra (bottom) for 383 (left) and 40S (right) at 9 = 40°.
The bottom spectrum is Optimized for decays occurring in vacuum by detecting the
kinetic energy Of the particle in the phoswich detector and using that velocity as the
value Of 6 in Eqn. 4.5.

109

curve of the particle distribution was analyzed to extract the scaling factors for use

in the following transformation.

E phos + B

(312/103 = A C I

(4.7)

where Ephos is the channel number Of the phoswich full-energy signal and A, B,
and C are constants chosen to produce resulting velocity distributions that had the
proper average velocities and minima and maxima of 6 ~0.05 and 0.12. In order
to check that the velocity distribution used for the Doppler correction was correct,
additional data were taken at the end Of the experiment with the phoswich position
shifted by :tl5 cm. The change in time Of flight was measured, confirming that the
adopted projectile velocity distribution was correct. Examples Of Doppler-corrected
7—ray spectra are shown in Fig. 4.7, and the difference between using a single value

of 6 and using a 6 determined event-by-event is clearly seen.

4.2.3 Angular Distribution of 7 Rays

Gamma-ray peak areas were obtained at each angle by integrating the counts in the
Ge spectra between limits set close to the high and low side Of the 2+ peaks, and
subtracting counts summed in background regions on the low and high energy sides of
the peak, normalized to the widths of the summing regions. For integrating the counts
in the magnet up/down Spectra, the summation limits were set using the summed
up+down spectra in order to ensure that no bias was introduced for the up and down
counting results, since the precession was deduced by comparing the counts in the up
versus the down spectra.

The angular distributions of the 1292 keV 38S 7 rays and the 903 keV 40S 7 rays are
shown in Fig. 4.8. The angular distribution has been transformed into the projectile
frame by correcting the number Of counts for the Lorentz boost Of the solid angle

and for the angle transformation between the laboratory and projectile frames. The

110

1.5
1.25 -

1 .
0.75 -
0.5

0.25 . . - - .
-180 -90 0 90 130 -l80 -90 0 90 180

6nuc [deg] Onuc [deg]

 

 

W(eflUC)

 

 

 

 

"v" .H

 

 

 

Figure 4.8: Angular correlation Of 38S (left) and 40S (right) 7 rays. The solid (dotted)
line shows the calculated angular correlation in the frame Of the projectile nucleus
including (not including) the effect of vacuum deorientation. The SeGA angles have
been transformed into the projectile frame.

former effect is more important, and can be expressed as [56]

1—52

(19 =
lab (1 + I3 cos 0mm)

 

2 dQn-uc, (4.8)

with the resulting angular distribution transformed as

4910.5
dQnuc .

 

I’Vnuc(0nucl : fl]lrlab(91ab) (4'9)

The transformation needed to Obtain the SeGA angles in the projectile frame is [120]

,8 — COS glob

cos 9 = .
mic [3 COS 9101, - 1

 

aim

The solid lines in Fig. 4.8 are calculations Of the angular distributions using the
program GKINT [107]. GKINT divides the target layers into thin slices and uses the
stopping powers Of Ziegler and coworkers [99] to compute the kinetic energy as a
function Of depth using a Runge-Kutta routine. Using the theory Of Coulomb excita-
tion [121], excitation cross sections in each slice are Obtained by integrating over the

phoswich scattering angle acceptance, and the kinetic energy of the particles emerging

111

from the target is calculated as a function Of excitation depth. The flight time from
excitation to target exit is used to determine the number Of 7 rays emitted during
target transit and in vacuum, using the 2+ lifetimes. The in-vacuum emission rates
are used tO create weighted-average values for the angular distribution coefficients
(a2 and a4) that are detected in the vacuum-flight peaks, including the effects Of
vacuum deorientation, which will be described in the next section. For 38S, the values
a2 = —0.3596 and a4 = —0.0597 were Obtained, and for 40S, the values a2 = —0.3850
and a4 = -0.0653 were obtained. Good agreement can be seen between the calculated
angular distribution and the measured values, and it is important to note that the

solid line is not a fit to the data except for the overall scaling.

4.2.4 Vacuum Deorientation

The 38S and 40S 7-ray peaks at 1292 keV and 903 keV in the Doppler-corrected
spectra appear sharp because those 7 rays were emitted by nuclei in flight through
vacuum after leaving the target, rather than during the slowing down process. Once
an ion exits the target, its velocity does not change until it is detected in the phoswich.
Therefore, the velocity detected in the phoswich detector will correspond to the ion
velocity at which the 7 ray was emitted, and those 7 rays will be Doppler corrected to
the true transition energy. It is advantageous from the point Of view Of HVTF g-factor
measurements to have a long-lived excited state (i.e., several picoseconds) for which
a large fraction Of decays occur after the ion leaves the target. However, from LISE
calculations as well as explicit charge-state measurements undertaken at ANU [122],
it was established that the sulfur ions were carrying 1-3 electrons upon emerging from
iron at velocities near Z v0. Therefore, recoil-in-vacuum (RIV) effects needed to be
examined to determine the extent Of their impact on the g-factor measurement.
Recoil-in-vacuum effects are interactions between the nucleus Of an ion and the
atomic electrons. The dominant interaction between orbital electrons and the atomic

nucleus is the magnetic dipole interaction between the nuclear magnetic moment and

112

the magnetic field created at the site of the nucleus by the ns electrons, with the Is
electron field (from Eqn. 2.31) being most significant since the fields for higher orbitals
fall Off with 1/n3. Additionally, it is important to note that only unpaired electrons
have a magnetic field contribution at the site of the nucleus, because the fields from
paired electrons cancel each other. Therefore, RIV effects are most important for
hydrogen- and lithium-like ions, which is most likely the situation for projectiles in
the HVTF regime or in any other scenario where projectiles are slowed to v x Z 710.
The nuclear angular momentum I and the electronic angular momentum J will couple
to a total angular momentum F = I + J, which is the conserved quantity for the free
ions recoiling through vacuum. The nuclear and electronic spins will precess about F,
and this serves to deorient the nuclear Spins from their initial orientation in Coulomb
excitation. While traversing a solid material, the ion-electron collision times are so
fast that electrons are constantly hopping on and Off the nucleus. This means that
nO spin precession can occur about the vector F because each new electron also
has a randomly-oriented spin direction. Upon emerging from material into vacuum,
however, an ion left with a few electrons attached quickly reaches its electronic ground
state. The static electronic configuration has a unique orientation for J , and hence
RIV effects can happen.

The deorientation is experimentally observed through an attenuation Of the an-
gular distribution, which is manifested by a reduction of the angular distribution
parameters a2 and a4 in a specific way. The reduction parameters, called attenua-
tion coefficients, are denoted 02 and G4. Goldring [59] discussed these attenuation
coefficients for a variety of scenarios, depending on the state Of ionization Of the
ions, by calculating the interaction between the nuclear 9 factor and the magnetic
fields created by the atomic electrons, which can be a complicated configuration. The

‘cold ionization’ scenario, where the ions are carrying only a few electrons, applies to

113

 

 

 

 

  
  

   
 
  

 

  
 

    

 

 

 

 

 

 

 

o
A 10 ' . r 0 . .
§- 4-937371% . ....................... '0 ............ ./04IIHII< ...............
g I ...................
--- ' 1 -l
g '0 E31319: "r 10 0.137 H-like -‘
c . 6 5' 5 e W r"/’

H-l'k _ ......_..._
g '04 r I e q— 16 g 104 . q: 16 __........_‘
3 q = :3 —- . q =15 _,
OD q = "”“"""""‘ q = 14 _________________
0.005 B

E: 3 0.0038are q=l3 —-~ / are q=l3 ——
U 10' m ‘ 10‘3 / L L

0.65 0.8 0.95 1.1 0.65 0.8 0.95 1.1

Reduced velocity (v/Zvo) Reduced velocity (v/Zvo)

Figure 4.9: Charge-state fractions for sulfur projectiles exiting iron foils near the sulfur
K -shell electron velocity, measured at the ANU 14UD Pelletron [122]. The velocity
0.083 c is marked on the left, corresponding to the 38S ions, and the velocity 0.088 c
is marked on the right panel, corresponding to the 40S ions.

high-velocity ions recoiling into vacuum. The attenuation coefficients in this case are

k(k +1) (.57)2
(21,-+1)21+(w*r)2’

 

Gk(w7')=1— (4.11)

and they depend on the lifetime of the excited state (T), and the Larmor frequency
(w) for the free-ion hyperfine field. Since the hydrogen- and lithium-like hyperfine
fields are so large, 027' approaches infinity and the attenuation reaches a “hard core”

value Of

E 1 — bk. (4.12)

The actual attenuation for the angular distribution will depend on the fraction Of ions
in the hydrogen— and lithium-like states, which depends on the velocity with which
the ions leave the target material, the ion atomic number Z, and the target material
itself. A measurement Of the charge state distribution of sulfur ions emerging out
of iron foils at velocities near the sulfur K —shell electron velocity was undertaken at
the ANU 14UD Pelletron [122]. Fits to the resulting charge state distributions are

plotted in Fig. 4.9. The attenuation coefficients for an ensemble of nuclei in different

114

few-electron charge states can be written as

Gk = 1 — (QH + QLilbk~ (4-13)

with the I = 2 values of ()2 = 0.24 and b4 2 0.8. The measured charge-state fractions
of Q H + QLz- = 0.435 yield attenuation coefficients Cg = 0.90 and G4 = 0.65. From
the angular distribution plots in Fig. 4.8, it can be seen that the overall effect of
these coefficients on the angular distributions are small; in fact, at the SeGA angles
used in the TF measurements, there is no sensitivity to the RIV effects except for
the slope changes they cause. The RIV effects calculated with the program GKINT
were included in the a2 and a4 coefficients used to calculate the slope factors S (9),

and uncertainties of 10% were assigned to the S factors.

4.3 Transient Field Effect

The transient field effect is observed through the Larmor precession of the angular
distribution. The total precession angle can be small due to the short interaction time
with the transient field. Therefore, an established analysis method using ratios of
counts in the germanium detectors is employed [51], rather than fits to the perturbed
angular distribution. In the following subsections the details Of the transient field

effect are given, including the extracted g(2+) values of 38’408.

4.3.1 Experimental Spin Precession

The precession angle A9 was extracted from the measured change in the angular
distribution for the two magnetic field directions, using the procedure described in
§ 2.1.3. The logarithmic slopes Of the 7 ray angular distribution W(9) were evaluated
at the SeGA angles 9 in the frame of the projectile due to Lorentz effects. The Lorentz

boost of the 7 ray angular distribution does not affect. the experimental precession

115

effect c, since the factors dfljab/dflnuc that affect N (9) will cancel in the double ratio
p. However, the spin precession occurs in the projectile frame, and it should be noted
that the SeGA Observation angles in the laboratory frame are not located at the
same angles in the projectile frame. Therefore, the slope factors were evaluated at the
projectile-frame SeGA angles, Obtained by transforming the SeGA angles according
to Eqn. 4.10, when calculating c / S. The S factors are shown in Table 4.1, along with

the extracted spin precessions for each angle pair.

4.3.2 Evaluation of Transient Field

The transient field strength was evaluated with the program GKIN T using the ANU
HVTF parametrization described in Chapter 2, and the experimental incoming beam
energies. The kinematical calculations, performed to Obtain 7-ray angular distribu-
tions weighted by excitation depth in the target and exponential decay in time, also
resulted in the knowledge of particle velocity profiles through the target. The particle
velocity profiles were integrated through the target (from the velocity at excitation,
to the velocity upon exiting the target) and the transient field interaction was cal-
culated using the TF parametrization. The particles statistically emit de—excitation
7 rays, which carry the information of the transient field interaction up until the time
of 7—ray emission. Therefore, the Doppler shifts of the 7 rays, which are essentially
measures Of the particle depths within the target at the times Of emission, reveal a
lineshape-dependent spin precession, in principle. The lineshape dependence is due to
the different respective transient-field interaction times for particles that emit 7 rays
at different depths in the target. The particles that emit 7 rays after leaving the
target have the longest transient-field interaction time and therefore the largest Spin
precession, A9, which is deduced by analyzing the 7 rays in the vacuum flight peak.
Therefore, the spin precession per unit 9 factor, A9/g, was evaluated in GKINT for
decays occurring in vacuum.

Results of the GKINT calculations are shown in Table 4.2. Spin precessions per

116

Table 4.1: Precession results obtained for each angle pair from the double ratios p

and 6, and the slope factors S.
388 408

 

detectors
9° (0° 6 S A9 [mrad] e S A9 [mrad]

21:29 90 -50(20) +0.79(8) —64(27) -17(35) +0.86(9) -20(41)
4151 90 +26(23) -O.80(8) -33(29) -4(35) -0.88(9) +5(40)
1139 46 -10(23) 054(5) +19(43) +7(28) 058(6) -12(48)
i147 143 -7(27) -0.48(5) +15(57) -24(31) -0.52(5) +46(60)
3:40 131 4109) +0.50(5) -83(39) +26(32) +0.53(5) +49(60)
i60 61 -18(15) +0.26(3) -69(58) +1(28) +0.27(3) +4(103)

 

 

 

 

 

 

 

average: —43( 15) +5(21)
Beam
¢=90°< Q -0° ¢=270°
+y if _
+x

 

+2 9=0°

Figure 4.10: Sign convention used for pair Of SeGA detectors at symmetric polar
angles from the beam axis. When seen from above, positive 9 is defined as increasing
in a counterclockwise fashion from the downstream beam axis.

117

Table 4.2: Nuclear parameters and reaction kinematics. (Etc) and ("0236) are the ion
kinetic energies and velocities at the entrance and exit of the iron layer. The effective
transient-field interaction time tefl‘ and the spin precession per unit 9 factor, A9/g,
are evaluated for ions that decay after leaving the target. The ion entrance and exit
velocities, teg, and A9/g were calculated with the program GKINT [107].

 

 

Isotope 1(23‘) (E1) (Ea) (vi/2m) (ve/Zvo) teff A49/9
[ps] [MeV] [MeV] [ps] [mrad]

388 4.9 762 123 1.75 0.71 2.98 —330
408 21 782 145 1.73 0.75 2.99 —339

 

 

 

unit 9 factor, (15 2 A9/ g, Of -330 mrad for 38S and -339 mrad for 40S were Obtained
from the GKIN T analysis for ions that decay after leaving the target. The effects Of the
parametrization and the 7-ray lineshape on the magnitude Of A9/g are investigated

in the following section.

Effects of BTF Uncertainty and Doppler Tail Effects

The uncertainty due to the parametrization was investigated by applying a quenching
factor to the v > Z 00 region of the BT 17(2)) function, since it is the high-velocity side
that does not have as much experimental verification as the low-velocity Side, which
was fit to a large number Of experimental measurements over the 6 S Z S 16 range.
Quenching the transient field strength on the high-velocity side did not appreciably
change the (15 value, because the field strength is small as v moves away from Z v0,
and the projectiles spend the least time at the higher velocities (as the ions slow,
more time is spent at the lower velocities where the transient field strength is higher).
It was seen that the parametrization does not have a significant impact on the spin
precession because the majority Of the precession comes from the 21 ~ Zvo region
where the field strength is well-characterized. There have been several measurements
using sulfur ions in the 1%2220 S v S Z110 velocity range, with a measurement at
v % éZvo in Fe and two measurements at v = Z v0 and v = éZvo in Gd, which

make up part Of the data set that was fit. to Obtain the ANU BTF parametrization.

118

Rather than attempting to deduce BTF(Z = 16; Fe) directly from only the sulfur
transient field measurements, it is more accurate to find a global BT F behavior, since
the transient field is a general phenomenon for light ions at high velocities as found in
Ref. [91]. Treating BT F as a global function is justifiable at high velocities (near Z vo)
because the projectile energy is well above the region where atomic shell effects and
possible molecular orbital effects come into play due to energy matching. Therefore
the uncertainty due to the global parametrization does not play a significant role in
the present measurements.

Since the 7-ray lineshape contains information on the velocity profile of the par-
ticles at the time of 7—ray emission, and the transient field is strongly velocity—
dependent, a study was undertaken to determine the effects of the Doppler tails on
the A9/ g value and the S values, particularly as a function of detector angle. Using
the GKINT routines, a Monte Carlo code was developed to analyze the effect Of the
lineshape on the precession effect [123]. The code generated incoming particle ener-
gies from a distribution that resulted in an outgoing particle energy spread resembling
the measured particle energy distribution shown in Fig. 3.2. Coulomb excitation cross
sections were calculated for randomly-picked interaction depths (hence energies) and
scattering angles within the phoswich acceptance, and were accepted or rejected in
a Monte Carlo fashion based on the Coulomb excitation cross section at that energy
and scattering angle. After excitation, the particle was traced through the target, and
a simulated 7 ray could be emitted along its path at any time after excitation, ac-
cording tO the lifetime of the state. Using the particle energy after leaving the target,
Doppler correction was performed in the same manner as the experimental analysis,
for each SeGA detector, and a Doppler-corrected 7-ray lineshape was generated. The
02 and a4 parameters were also calculated as a function of lineshape to see the differ-
ence between the angular distribution of the Doppler tails and the vacuum peaks. It
was concluded that the spin precession per unit 9 factor is dominated by the vacuum

flight peak (accounting for 92% Of the precession Of the total lineshape) and that the

119

 

 

    

 

1000 E 1

'8 . Mrs- Spin Precession (g=1)
"‘ I ‘1.»

> E. '

o r: 100 5

£2 1

:2 g :

t:

3 8 ’

U i 10 -
.E ' .
Q. '1
m I,

 

 

l I 1 1 1 n 1 ,7
1350 1400 1450 1500 1550 1600 1650
Lab Ge Energy [keV]

 

Figure 4.11: Doppler-shifted 7-ray spectrum modeled using Monte Carlo tech-
niques [123], overlayed with the lineshape-dependent spin precession per unit 9 factor.
The detected spin precession is dominated by contributions from the in-vacuum flight
peak (see text).

Doppler tails do not impact the g factor measurement. The Monte Carlo modeling of
the 7-ray lineshapes yielded A9 / g values in agreement with the GKINT values shown
in Table 4.2.

4.3.3 2+ 9 factors of 38ADS

The 9 factors for the 2? levels of 38S and 40S are deduced to be +0.13(5) and —0.01(6),
respectively, using the Spin precessions from a weighted average Of the six SeGA angle
pairs, shown in Table 4.1, and the calculated transient field strength in the spin pre-
cession per unit g factor shown in Table 4.2, which has been shown to be independent
of the TF parametrization. The uncertainty in BT F1 Obtained from the change upon
quenching the high velocity TF strength discussed in the previous section, did not
change the final uncertainty when added in quadrature with the statistical errors. All
the available g(2+) values for sulfur isotopes are shown in Fig. 4.12. The errors in
the 9 factors are dominated by the 7 ray counting statistics. Interpretations of the

g-factor results for 38,408 are given in the following section, within the framework Of

120

 

0.8

 

 

 

 

0.4 " + W ............................... ZA
9 I
U I
«‘3 O 1
Of)
-0.4
Schmidt va2 ..............................................
-0.8

14 16 18 20 22 24 26 28 30
Neutron Number

Figure 4.12: Comparison Of previously-measured 9 factors Of stable sulfur isotopes
(Open points) and the present radioactive isotope results (filled points). Also shown
is the Schmidt value for f7/2 neutrons, the orbital that begins to fill for N 2 20, and

the usual collective Z /A value, representing the two naive predictions for the single
particle and the collective effects, respectively.

the shell model.

4.4 Interpretation of Measured g(2+) values

Shell model calculations were undertaken in order to interpret the excited-state g fac-
tors. The shell model calculations were performed with the program OXBASH [124],
using an sd-p f model space where valence protons were restricted to the sd Shell
and valence neutrons were restricted to the pf shell. Calculations were performed for
1132820, 333822, and [118824, and their isotones fli’gArgO, ‘ngrgg and §§Ar24. The Nowacki
interaction, which was developed in Refs. [20, 23] for application to neutron-rich nu-
clides near N = 28, was used for the shell-model calculations. The effective interaction
consists of USD matrix elements [125] for the proton-proton interactions, the KB’ ma-
trix elements [126] for the neutron-neutron interactions, and the Kahana, Lee, and

Scott G matrix [127] for the proton-neutron interactions.

121

Table 4.3: Shell model and experimental values of excitation energies and B (E2; 0? —+

2?) values for 36'38’40S and 38"10:42Ar nuclides, as well as predicted quadrupole mo-
ments, and deformation parameters extracted from the predicted quadrupole mo—
ments, predicted B(E2) values, and measured B(E2) values.

 

 

 

 

Nuclide E31, 1331‘." Q(2f) B(E2)Ih B(E2)€XP 135‘ mg] leexpl
(MeV) (MeV) e fm2 e2 fm4
363 3.426 3.291 -ll.6 141 104(28) +0.21 0.20 017(2)
383 1.531 1.292 —9.7 268 235(30) +0.17 0.26 025(2)
40S 0.980 0.904 —19.3 473 334(36) +0.33 0.34 028(2)
38Ar 2.018 2.167 +4.3 178 130(10) —0.07 0.19 016(1)
40Ar 1.371 1.461 +8.5 263 330(40) —0.13 0.22 025(2)
42Ar 1.243 1.208 +5.2 360 430(10) —-0.08 0.25 028(3)

 

 

The calculations gave first-excited 2+ state energies, reduced transition proba-
bilities to the first 2+ states, and quadrupole moments Of the 2+ states of 36’38’4OS
and 38"10:42Ar. Standard effective charges of ep ~ 1.5 and en ~ 0.5 were used. The
E(2+) and B(E2) values are in good agreement with experimental values, and are
compared in Table 4.3. In general, the excited-state energies are reproduced within
about 200 keV, and the trends in the B(E2) values are correctly reproduced.

The g(2+) values were obtained from the calculated shell-model wavefunctions, us-
ing bare nucleon 9 factors. The calculated 9 factors are compared to the experimental
values in Fig. 4.13. The 38"1°Ar values were recently measured using the conventional
transient field method [128, 129]. The 36S g(2+) value has not been measured due to
its high excitation energy and Short lifetime, although 36S is stable. The 42Ar g(2+)
value has not been measured, since 42Ar is radioactive. The calculations reproduce
the g-factor behavior of the sulfur and argon isotopes fairly well, and therefore some
insight can be gained by examining the shell model results in more detail.

The individual contributions of the neutrons and protons are separated in Ta-

ble 4.4, and one can see that a near cancelation Of the neutron and proton 9 factors

122

 

 

 

 

 

 

 

 

 

1.5 _ F m I 4_ I I I .1
C sulfur isotopes (Z=16) j; argon isotopes (Z=18):
. t .. .
._ \ .1... .1
H 1'0 : \ :: .
g C \\ :: 2
“:3 0.5: I. —Z/A: ¥ *Z/A:
or) I- \ 4+ f «
: \5 :: \
0.0 _ = r‘ ,_ \ r ,.
. 4.. \r”
-O 5 P l 1 l .. 1 1 1
' 20 22 24 20 22 24
Neutron number Neutron number

Figure 4.13: Comparison Of measured 9 factors and shell model predictions (see text).
The nearest stable nuclides to the sulfur isotopes for which g-factor results are avail-
able, the argon isotopes (Refs. [128] and [129]), are shown on the right hand side. and
the results of the current work are shown on the left.

is causing the small values for the N = 22 and N = 24 isotones. At N = 20, the 2+
wavefunction of 36S consists exclusively Of proton excitations within this model space,
leading to a large positive 9 factor that comes primarily from coupling an 31/2 proton
and a d3/2 proton to spin 1 = 2. As neutrons begin to occupy the f7/2 single-particle
orbital, there are still positive contributions Of proton excitations in the 21+ states of
38S and 40S. The proton excitations effectively cancel the neutron contributions that
would naively dominate the excited-state wavefunctions due to the close proximity of
the N = 20 shell closure.

TO investigate the neutron shell structure further, the shell model calculations were
repeated with the valence neutrons restricted to the f7/2 orbital (i.e.. not allowed to
be excited to the other orbitals in the f p shell that lie above the N = 28 gap). The
results Of the two shell model calculations are displayed in Table 4.5. There were two
substantial differences in the g(2+) value for 38S. The neutron g factor resembles more
closely the f7/2 Schmidt value, as would be expected when the neutrons are restricted
to the f7/2 orbital. However, an additional effect is that the proton 9 factor becomes

less positive, due to fewer proton excitations contributing to the 2+ wavefunction.

123

Table 4.4: Theoretical 9 factors compared with experiment. The measured 38Ar and
40Ar 9 factors are from Refs. [128] and [129], respectively.

 

 

N “Glide ggrl‘oton gitgutron gth 9 exp
ffgsgo +1.15 0 +1.15

£33322 +0298 —0301 -—0.0026 +0.13(5)
[(3324 +0276 —0241 +0035 —001(6)
18Ar20 +0308 0 +0308 +024(12)
1143mm, +0164 —0364 —0200 —002(4)
18mg, +0220 —0280 -0.060

 

 

 

Table 4.5: Detailed decomposition of theoretical 9 factors compared with experiment.
The 9 factors that result from restricting valence neutrons to the f7 2 orbital are in
the rows labeled with the “SDF” model space, and the 9 factors from the full f p shell
calculation are in the rows labeled with the “SDPF” model space.

 

 

 

NUClide IVIOdel g $1601.01] gillflflltl‘oll g”! gexp
orbital spin total
‘13:;820 SDPF 0.9667 0.1868 +1.15 0 +1.15

16822 SDPF 0.2248 0.0733 +0298 —0301 —0003 +0.13(5)
SDF 0.087 0.049 +0136 —0494 —0.358

16824 SDPF 0.2249 0.0513 +0276 —0241 +0035 -0.01(6)
SDF 0.2489 0.0695 +03184 —04037 —00853

11%sz SDPF 1.1509 —0.8424 +0308 0 +0308 +024(12)

313mm SDPF 0.3109 —0.1470 +0164 —0364 —0200 —002(4)
SDF 0.2278 —00922 01356 —04310 —02954

43mg, SDPF 0.2629 —0.0431 +0220 —0.280 —0060
SDF 0.2421 -—0.0254 0.2167 -—04167 —02000

 

 

 

124

The proton excitations seem to be very closely linked to neutron excitations into the
higher orbitals in the fp shell, and the resulting theoretical g factor changes from
-0.358 (with neutrons restricted to f7/2), to —0.003 (when neutrons are allowed to
spread into the full f p shell).

Details of the contributing components to the 2+ wavefunctions of 38S and 40S
are Shown graphically in Fig. 4.14. The contributions of the main proton partitions
are shown for both the restricted and full- f p shell model calculations. There are eight
valence protons available to partition among the three sd orbitals (d5/2, 31/2, and
113/2) beyond the Z = 8 closed shell at oxygen. The different combinations shown in
Fig. 4.14 give rise to different g factors. depending upon how the different particles are
coupled together. and for the simple cases. g factors are shown next to the relevant
partition. Not Shown are Off-diagonal components, such as 1/ f7 /2 — uf5/2 or ads/2 —
7rd3/2 terms which have A8 = 0 and can therefore contribute to the magnetic moment.
The difference between the restricted and full calculation is more pronounced in 38S
than in 4OS. The N = 22 nucleus 38S has one pair of valence neutrons, and the two
neutrons were found by the full f p shell calculations to be partitioned as 63.2% f72/2,
24.4% f71/2p31/2, and 4.2% f732p32/2, taking into account partitions that contribute
1% or more to the total wavefunction. In the 2f“ wavefunction of 38S, it is more
favorable for a proton to be excited from s1/2 to d3 /2 when the neutrons are also
allowed to be excited. Such a strong influence Of the neutron occupancies on the
proton configuration indicates that they are strongly coupled together, which is a
prerequisite for the onset of deformation. When two additional neutrons are added to
make 40S, the shell model calculations show that they are partitioned as 47% f74/2,
25% f73/2p31/2, and 8.8% f72/2p32/2, again only accounting for partitions that contribute
1% or more to the total wavefunction. There is not as large of a change in the 40S
proton configuration, which was already composed of several proton excitations to

begin with than in the case Of 38S. Closer inspection Of Table 4.5 shows that the

change in the 40S 9 factor when neutrons are restricted tO f7/2 is due almost entirely

125

ds/z S1/2 d3/2 proton partitions, 38S

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1
6 2 0 , '{0:0$1631.33}???{0:80:03z0:9:on3103;033:9319}:01910319101910: -
=0
6 1 l g" .
..........
6 0 2 ' ‘
4 2 2 ‘
5 2 1
f ‘~:
5 1 2 V 7’2 :
full vfp I
4 A A 1 A A A l A A A l
0 20 40 60

Fraction of 2+ Wavefunction [%]

d5/2 S1/2 d3/2 proton partitions, 408

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6 O 2 gp=0.08 ‘
6 l 1 gp=l 1 ‘
6 2 O ‘
5 2 l ‘
4 2 2 ‘
va2 5:23:35: ,
full vfp I
40 6O

Fraction Of 2+ Wavefunction [%]

Figure 4.14: Calculated components of the 2+ wavefunctions of 38S (left) and 40S
(right). The main proton partitions are labeled and the horizontal bars represent the
fraction Of the 2Jr wavefunctions. Two different shell model results are shown. The
hatched bars are results Of a restricted shell model calculation that forced the valence
neutrons to occupy the f7/ 2 orbital, and the solid bars Show the results Of the full f p
shell calculation that allows neutrons to occupy higher-lying orbitals. The Schmidt
values for particular proton partitions are Shown for selected cases.

126

to the change in the neutron g factor contribution.

The 403 g(2+) result is interesting because it is nearly zero, although 40S is as-
sumed to be deformed according to the E(2+) and B(E2) systematics. The shell
model calculations for 40S predict a quadrupole moment that leads to a deformation
parameter of 62 z +0.3, which agrees with the value extracted from the B(E2) mea-
surement [10]. Therefore the onset Of deformation is confirmed by the shell model
calculations that explain the g-factor results. However, the magnetic behavior of the
deformed 40S nucleus is manifested in an unprecedented way, since in general, de-
formed nuclei have been found to exhibit 9 factors that are the result Of the orbital
motion of the protons alone (since the neutrons have 93 = 0), with intrinsic spin con-
tributions Of the neutrons and protons usually assumed to be very small in collective
models. Traditionally, it has been assumed that deformed nuclei should always have
9 factors near Z /A, since without exception the deformed, stable, even-even nuclei
with known 9 factors (whether they are the A ~ 150 rare earth elements, or deformed
light N = Z nuclei) have g(2+) z Z /A. Large neutron spin contributions to the g(2+)
value in a deformed nucleus, as is seen in 40S, have not been reported previously in
the literature. The Z /A collective prediction for g(2+) is not generally regarded to
be well-justified theoretically, but it is generally held as an empirical rule. There is
now evidence that Shell effects, which lead to large neutron spin contributions to the

9 factor, can overcome the tendency for deformed nuclei to have gyromagnetic ratios

Of +Z/A.

4.5 Summary and Outlook

The successful measurement of the 2‘] g factors of 38S and 40S demonstrates that
the transient. field technique can be applied successfully to fast fragment beams. Shell
model results indicate that the transitional neutron-rich sulfur isotopes Show shell

effects indicative Of active neutron shell gaps around the f7/2 orbit. The filling Of the

127

 

 
 
 

,Brocess

”19E“ path
—— 50
Yield #/s

I> 10,:
l> 10
l> 10
> .015
l> 10

 
 
  

 

Proton Number

 

 

Neutron Number

 

Figure 4.15: Calculated yields of the Coupled Cyclotron Facility. The 408 measure-
ment reported in this dissertation indicates that fast-fragment g factors of short lived
states can be measured with rates of 104 particles per second.

f7/2 neutron orbit, however, affects the protons through an interaction that causes the
31/2-(13/2 gap to change, as has been the subject of much recent experimental [130] as
well as theoretical [27] interest. In addition, it was observed that the neutron occupa-
tion of higher-lying states across the N = 28 gap has a profound effect on the proton
occupancies, indicating that the neutrons and protons are strongly coupled. Since
the g factors contain information on orbital components of the wavefunction, they
offer a sensitive and unique means to assess how shifts in the single-particle energy
spectrum and orbital occupancies occur. The transition from single-particle behavior
to collectivity that occurs in the neutron-rich sulfur isotopes has been investigated
by measurements Of their g(2+) values, allowing a microscopic study of the dramatic
deformation change, since the 9 factors provide stringent constraints on the nuclear
wavefunctions of 388 and 40S.

Earlier measurements on the stable 32’348 isotopes [131] can be compared to the
present work, as shown in Fig. 4.12. The similar magnitudes of the error bars, given

the disparity in the beam intensities for the stable versus the radioactive sulfur mea-

128

 

1 m : t T F T :
l
‘9 I
10 I Cf .59 3°". ”'0,

E. ~ 3' , a. .5
7 .7: or!!! W.» I,» "Ge ii 1
D: f. , 5 ii i‘. :2 J ""4 “'0 _ . i
2? ‘Mat‘ .4" 9!”?2 «Zn»; «

-l , _I '-’ ’.' '._. i 1'

.a l i- 110. ' i! f, 5: '0 I I I

3 E fig-«33!. ’ .‘f'Ti
0. 1 [r i ca . ,7"; l
’ '5 v’Ni :
0.0 l

 

 

 

10 20 30 40 50 60 70 80 90
Mass Number

Figure 4.16: Estimates of first 2+ lifetimes (using the Grodzins relation, for levels
with energies less that 3 MeV) in radioactive even-even nuclides. Isotopic chains are
labeled and connected by lines.

surements, shows that the HVT F technique can take advantage of the high projectile
energies and utilize thick ferromagnet layers, which increases the interaction time
with BTF- In fact, most of the interaction time is Spent near Z v0 as the fast beam
slows through the peak Of the transient field strength, whereas most stable-beam
measurements are moving away from Zvo as they slow down in the target. For the
specific cases of the stable sulfur isotopes shown in Fig. 4.12, the interaction time with
the transient field is low due to the short excited state lifetimes Of less than 0.5 ps,
but the beam currents available were on the order of 1011 particles per second. The
4OS measurement demonstrates that the HVTF method is applicable to beams with
intensities as low as 104 particles per second. The yield Of secondary fragments at the
Coupled Cyclotron Facility is shown in Fig. 4.15, and a large region Of the nuclear
chart is available through fragmentation, at the 104 pps rate.

The transient field measurement Of g(2+) with fast. projectiles depends on slowing
the projectiles down in a timely manner, governed by the 2+ lifetime T. The T(2+)

value Of a particular isotope provides the severest constraint on the feasibility of an

129

HVTF measurement. In general, the time to slow ions in the target material will
depend on the ion atomic number Z and the incoming projectile velocity. The final
velocity should be near Z m, which means higher-Z ions need less slowing than lower
Z ions, although the behavior of the transient field at high velocities for higher Z ions
is not well studied, and the maximum transient field strength may no longer be at
the velocity Z v0. Since stopping powers are increasing functions of the projectile Z ,
it appears that moving higher in the periodic table is advantageous from a projectile
velocity point of view, but much work is still needed to characterize the transient
field for Z > 16 ions at and above the K -shell electron velocity. This is a difficult
task to undertake at a radioactive beam facility, but may be well-suited for stable-
beam facilities since in general stable isotopes have known g(2+) values and are
available with high beam intensities. An estimation of 2? lifetimes (using the Grodzins
relation [12]) in radioactive even-even nuclei with Z S 32 is shown in Fig. 4.16 for
states with an excitation energy of 3.3 MeV or less. It can be seen that a large number
of states have favorable lifetimes for the HVTF method. One such state, the 2? state
in the deformed N = 20 nucleus 32Mg, has a similar energ, B(E2) value, and lifetime
to 408. Shell effects, such as those observed in the neutron-rich sulfur isotopes, are
not expected to be strong (since the N = 20 shell closure has disappeared at Z = 12),
and a g(2+) measurement is highly desirable to investigate the magnetic properties

of a well-deformed nucleus located far from stability.

130

Bibliography

[1] A. de—Shalit and M. Goldhaber, Phys. Rev. 9221211 (1953).

[2] T. Lee, D. A. Papanastassiou, and G. J. Wasserburg, Astrophys. J. 2202L21
(1978).

 

[3] O. Sorlin, D. Guillemaud-lVIueller, A. C. Mueller, et al., Phys. Rev. C 4722941
(1993).

[4] N. J. Davis, J. A. Kuehner, A. A. Pilt, A. J. Trude], M. C. Vetterli, C. Bamber,
E. K. Warburton, J. W. Olness, and S. Raman, Phys. Rev. C 322713 (1985).

[5] J. W. Olness, E. K. Warburton, J. A. Becker, D. J. Decrnan, E. A. Henry, L. G.
Mann, and L. Ussery, Phys. Rev. C 3422049 (1986).

[6] E. K. Warburton, D. E. Alburger, and G. Wang, Phys. Rev. C 362429 (1987).

[7] L. K. Fifield et al., Nucl. Phys. A 4172534 (1984).

[8] W. A. Mayer et al., Z. Phys. A 319:287 (1984).

[9] B. Fornal, R. H. Mayer, I. G. Bearden, et al., Phys. Rev. C 4922413 (1994).
[10] H. Scheit, T. Glasmacher, B. A. Brown, et al., Phys. Rev. Lett. 7723967 (1996).

[11] T. Glasmacher, B. A. Brown, M. J. Chromik, et al., Phys. Lett. B 3952163
(1997)

[12] S. Raman, C. W. Nestor, and P. Tikkanen, At. Data Nucl. Data Tables 78:1
(2001)

[13] H. Scheit, Ph.D. thesis, Michigan State University (1998).
[14] P. Moller and J. R. Nix, At. Data Nucl. Data Tables 262165 (1981).

[15] T. R. Werner, J. A. Sheikh, W. Nazarewicz, M. R. Strayer, A. S. Umar, and
M. Misu, Phys. Lett. B 3352259 (1994).

[16] T. R. Werner, J. A. Sheikh, M. Misu, W. Nazarewicz, J. Rikovska, K. Heeger,
A. S. Umar, and M. R. Strayer, Nucl. Phys. A 5972327 (1996).

[17] J. A. Winger, P. F. Mantica, R. M. Ronningen, and M. A. Caprio. Phys. Rev.
0' 642064318 (2001).

131

 

[18] J. H. Kelley, T. Suomijarvi, S. E. Hirzebruch, et al., Phys. Rev. C 5621206
(1997)

[19] D. Sohler, Z. Dombrédi, J. Timar, et al., Phys. Rev. C 662054302 (2002).

[20] J. Retamosa, E. Caurier, F. Nowacki, and A. Poves, Phys. Rev. C 5521266
(1997)

[21] D. Hirata, K. Sumiyoshi, B. V. Carlson, H. Toki, and I. Tanihata, Nucl. Phys.
A 6092131 (1996).

[22] G. A. Lalazissis et al., Phys. Rev. C 602014310 (1999).

[23] S. Nummela, P. Baumann, E. Caurier, P. Dessagne, A. Jokinen, et al., Phys.
Rev. C 63:044316 (2001).

[24] E. Caurier, F. Nowacki, and A. Poves, Nucl. Phys. A 742214 (2004).

[25] R. Rodriguez-Guzman, J. L. Egido, and L. M. Robledo, Phys. Rev. C 652024304
(2002)

[26] J. Dobaczewski, I. Hammamoto, W. Nazarewicz, and J. A. Sheikh, Phys. Rev.
Lett. 722981 (1994).

[27] T. Otsuka, T. Suzuki, R. Fujimoto, H. Grawe, and Y. Akaishi, Phys. Rev. Lett.
952232502 (2005).

[28] P. D. Cottle and K. W. Kemper, Phys. Rev. C 5823761 (1998).

[29] T. Motobayashi, Y. Ikeda, Y. Ando, et al., Phys. Lett. B 34629 (1995).

[30] S. Grévy, J. C. Angelique, P. Baumann, et al., Phys. Lett. B 5942252 (2004).
[31] J. Fridmann, I. Wiedenhover, A. Gade, et al., Nature 4352922 (2005).

[32] P. D. Cottle and K. W. Kemper, Phys. Rev. C 662061301 (2002).

[33] J. D. Jackson, Classical Electrodynamics (John Wiley and Sons, New York.
1975)

[34] J. M. Blatt and V. F. Weisskopf, Theoretical Nuclear Physics (Wiley. New York.
1952)

[35] S. Eidelman, K. G. Hayes, K. A. Olive, et al., Phys. Lett. B 59221 (2004).
Particle Data Group.

[36] M. G. Mayer, Phys. Rev. 7521969 (1949).
[37] O. Haxel, J. H. Jensen, and H. E. Suess, Phys. Rev. 7521766 (1949).

[38] D. M. Brink and G. R. Satchler, Angular Momentum. (Clarendon Press. Oxford,
1968)

132

[39] B. Castel and I. S. Towner, Modern Theories of Nuclear Moments (Clarendon
Press, Oxford, 1990).

[40] B. A. Brown, Prog. Part. Nucl. Phys. 472517 (2001).

[41] E. Caurier, G. Martinez-Pinedo, F. Nowacki, A. Poves. and A. P. Zuker, Rev.
Mod. Phys. 772427 (2005).

[42] R. J. Blin-Stoyle, Rev. Mod. Phys. 28:75 (1956).
[43] N. J. Stone, At. Data Nucl. Data Tables 90275 (2005).
[44] N. Bohr, Nature 1372344 (1936).

[45] A. Bohr and B. Mottelson, Nuclear Structure, vol. 2 (Benjamin, New York,
1975)

[46] A. E. Stuchbery, Nucl. Phys. A 5892222 (1995).
[47] O. Prior, F. Boehm, and S. G. Nilsson, Nucl. Phys. A 1102257 (1968).
[48] E. L. Brady and M. Deutsch. Phys. Rev. 782558 (1950).

[49] R. R. Borchers, J. D. Bronson, D. E. Murnick, and L. Grodzins, Phys. Rev.
Lett. 1721099 (1966).

[50] R. R. Borchers, B. Herskind. J. D. Bronson, L. Grodzins. R. Kalish. and D. E.
Murnick, Phys. Rev. Lett. 202424 (1968).

[51] N. Benczer-Koller, M. Hass. and J. Sak, Annu. Rev. Nucl. Part. Sci. 30:53
(1980)

[52] K.-H. Speidel, O. Kenn, and F. Nowacki, Prog. Part. Nucl. Phys. 49:91 (2002).
[53] N. Rud and K. Dybdal, Phys. Scr. 342561 (1986).

[54] R. D. Gill, Gamma-Ray Angular Correlations (Academic Press, New York,
1975)

[55] T. Yamazaki, Nucl. Data A 321 (1967).
[56] A. E. Stuchbery, Nucl. Phys. A 723269 (2003).

[57] H. Morinaga and T. Yamazaki, In-Beam Gamma-Ray Spectroscopy (North-
Holland, Amsterdam, 1976).

[58] K. S. Krane, in Low Temperature Nuclear Orientation. edited by N. J. Stone
and H. Postma (Elsevier Science. 1986). p. 31.

[59] G. Goldring. in Heavy Ion Collisions. edited by R. Bock (North-Holland. Am-
sterdam, 1982), vol. 3, p. 484.

133

 

 

[60] R. E. Tuzan, P. Burkhardt, and D. Secrest, Comp. Phys. Comm. 1122112
(1998).

[61] H. Frauenfelder and R. M. Steffen, in Alpha-, Beta-, and Gamma-Ray Spec-
troscopy, edited by K. Siegbahn (North-Holland, Amsterdam, 1965).

[62] M. J. L. Yates, in Alpha-, Beta-, and Gamma—Ray Spectroscopy, edited by
K. Siegbahn (North—Holland, Amsterdam, 1965).

[63] D. C. Camp and A. L. Van Lehn, Nucl. Inst. 62 Meth. 762192 (1969).
[64] A. E. Stuchbery, private communication.

[65] U. Reuter, F. Hagelberg, S. Kremeyer, H.-J. Simonis, K.-H. Speidel, M. Knopp,
W. Karle, J. Cub, P. N. Tandon, and J. Gerber, Phys. Lett. B 230216 (1989).

[66] P. C. Zalm, A. Holthuizen. J. A. G. de Raedt. and G. van Middelkoop, Phys.
Lett. B 692157 (1977).

[67] A. E. Stuchbery, C. G. Ryan, H. Ohnuma.. G. B. Beard, and H. H. Bolotin,
Phys. Rev. C 272434 (1983).

[68] M. B. Goldberg, W. Knauer, G. J. Kumbartzki, K. H. Speidel. J. C. Adloff, and
J. Gerber, Hyperfine Interact. 4:262 (1978).

[69] J. L. Eberhardt, R. E. Horstmann, P. C. Zalm, H. A. Doubt, and G. van Mid-
delkoop, Hyperfine Interact. 3:195 (1977).

[70] N. K. B. Shu, D. Melnik, J. M. Brennan, W. Semmler, and . Benczer-Koller,
Phys. Rev. C 2121828 (1980).

[71] H. R. Andrews, O. Hausser, D. Ward, P. Taras, R. Nicole, J. Keinonen,
P. Skensved, and B. Haas, Nucl. Phys. A 3832509 (1982).

[72] C. Caso et al., Eur. Phys. Jour. C321 (1998).
[73] J. F. Ziegler, J. Appl. Phys. 8521249 (1999).

[74] S. Schielke, K.-H. Speidel, O. Kenn, J. Leske, N. Gemein, M. Offer, Y. Y.
Sharon, L. Zamick, J. Gerber, and P. Maier-Komor, Phys. Lett. B 5672153
(2003).

[75] K.-H. Speidel, J. Leske, S. Schielke, S. C. Bedi, O. Zell, P. Maier-Komor, S. J. Q.
Robinson, Y. Y. Sharon, and L. Zamick, Phys. Lett. B 6332219 (2006).

[76] O. Kenn, K.-H. Speidel, R. Ernst, S. Schielke, S. Wagner. J. Gerber. P. Maier-
Komor, and F. Nowacki, Phys. Rev. C 652034308 (2002).

[77] K.-H. Speidel, N. Benczer-Koller. G. Kumbartzki, et al.. Phys. Rev. C 5722181
(1998).

[78] J. Holden, N. Benczer-Koller, G. Jakob, et al.. Phys. Rev. C 632024315 (2001).

134

 

[79] G. Kumbartzki, J. R. Cooper, N. Benczer~Koller, et al., Phys. Lett. B 5912213
(2004).

[80] R. E. Laxdal, Nucl. Inst. Meth. Phys. Res. B 2042400 (2003).

[81] O. Kester, T. Sieber, S. Emhofer, et al., Nucl. Inst. Meth. Phys. Res. B 204220
(2003).

[82] C. J. Gross, J. Phys. C 31281639 (2005).
[83] D. H. Feng, A. R. Barnett, and L. J. B. Goldfarb, Phys. Rev. C 1321151 (1976).
[84] A. N. F. Aleixo and C. A. Bertulani, Nucl. Phys. A 5052448 (1989).

[85] T. Czosnyka, D. Cline, L. Hasselgren, C. Y. Wu, R. M. Diamond, H. Kluge,
’ C. Roulet, E. K. Hulet, R. W. Lougheed, and C. Baktash, Nucl. Phys. A 4582123
(1986).

[86] H. Goldstein, Classical Mechanics (Addison-Wesley, 1980).
[87] J. Cub, M. Bussas, K. H. Speidel, et al., Z. Phys. A 34521 (1993).
[88] U. Grabowy, K. H. Speidel, et al.. Z. Phys. A 3592377 (1997).

[89] O. Kenn, K.—H. Speidel, R. Ernst, J. Gerber, P. Maier-Komor, N. Benczer-
Koller, G. Kumbartzki, and F. Becker, Nucl. Inst. Meth. Phys. Res. B 1712589
(2000).

[90] H. Ueno, W. Sato, H. Ogawa, et al., Hyperfine Interact. 1362211 (2001).
[91] A. E. Stuchbery, Phys. Rev. C 692064311 (2004).
[92] H. Kopfermann. Nuclear Moments (Academic Press, New York, 1958).

[93] O. Tarasov, D. Bazin, M. Lewitowicz, and O. Sorlin, Nucl. Phys. A 7012661
(2002)

[94] D. Bazin, O. Tarasov, M. Lewitowicz, and O. Sorlin, Nucl. Inst. Meth. Phys.
Res. A 4822307 (2002).

[95] A. E. Stuchbery. A. N. Wilson, P. M. Davidson, A. D. Davies, T. J. Mertzimekis,
S. N. Liddick, B. E. Tomlin, and P. F. Mantica, Phys. Lett. B 611281 (2005).

[96] J. F. Ziegler, SRIM Special — Degrader Blocks (2000), Documentation of com-
puter program SRIM.

[97] J. F. Ziegler. Computer program SRIM, Tech. rep., IBM-Research, New York
(2000).

[98] D. J. Morrissey, B. M. Sherrill, M. Steiner, A. Stolz. and I. Wiedenhoever, Nucl.
Inst. Meth. Phys. Res. B 204290 (2003).

[99] J. F. Ziegler, J. P. Biersack, and U. Littmark, in The stopping and ranges of
ions in matter, edited by J. F. Ziegler (Pergamon, New York, 1985), vol. 1.

[100] H. Geissel and C. Scheidenberger, Nucl. Inst. Meth. Phys. Res. B 1362114
(1998).

[101] J. Lindhard and A. H. Sorensen, Phys. Rev. A 5322443 (1996).

[102] W. F. Mueller, J. A. Church, T. Glasmacher, D. Gutknecht, G. Hackman, P. G.
Hansen, Z. Hu, K. L. Miller, and P. Quirin, Nucl. Inst. Meth. Phys. Res. A
4662492 (2001).

[103] Z. Hu, T. Glasmacher, W. F. Mueller, and I. Wiedenhéver, Nucl. Inst. Meth.
Phys. Res. A 4822715 (2002).

[104] H. Zwahlen, Ph.D. thesis, Michigan State University (2005).

[105] K. L. Miller, T. Glasmacher, C. Campbell, L. Morris, W. F. Mueller, and
E. Strahler, Nucl. Inst. Meth. Phys. Res. A 4902140 (2002).

[106] G. Jakob. N. Benczer-Koller, G. Kumbartzki, et al., Phys. Rev. C 652024316
(2002)

[107] A. E. Stuchbery, Some notes on the program CKIN T: Transient-field g—factor
kinematics at intermediate energies, Department of Nuclear Phyiscs, The Aus-
tralian National University, report no. ANU-P/ 1678 (2005).

[108] J. A. Church, C. M. Campbell, D.-C. Dinca, et al., Phys. Rev. C 722054320
(2005)

[109] BC—400 andBC-M4 Plastic Scintillator Data Sheets (2005), Saint Gobain Ce-
ramics & Plastics, Inc., Newbury, OH.

[110] W. F. Mueller, private communication.
[111] E. Hecht, Optics (Addison-Wesley, 1998).

[112] R. Fox, E. Kasten, K. Orji, C. Bolen, C. Maurice, and J. Venema, IEEE Trans.
Nucl. Sci. 512571 (2004).

[113] C. Kittel, Introduction to Solid State Physics (Wiley, New York, 1996).
[114] N. Benczer—Koller, private communication.

[115] A. Piqué, J. M. Brennan, R. Darling, R. Tanczyn, D. Ballon, and N. Benczer-
Koller. Nucl. Inst. Meth. Phys. Res. A 2792579 (1989).

[116] K. Halbach and R. F. Holsinger, Particle Accelerators 7:213 (1976).

[117] J. H. Hubbell and S. M. Seltzer, Tables of X-Ray Mass Attenuation Coeffi-
cients and Mass Energy-Absorption Coefi‘icients, Tech. rep., National Institute
of Standards and Technology, Gaithersburg, MD (2004).

136

[118] O. Klein and Y. Nishina, Z. Physik 522853 (1929).
[119] T. Glasmacher, Annu. Rev. Nucl. Part. Sci. 4821 (1998).

[120] D. Pelte and D. Schwalm, in Heavy Ion Collisions, edited by R. Bock (North-
Holland, Amsterdam, 1982), vol. 3, p. 1.

[121] C. A. Bertulani, A. E. Stuchbery, T. J. Mertzimekis, and A. D. Davies, Phys.
Rev. 0 682044609 (2003).

[122] A. E. Stuchbery, P. M. Davidson, and A. N. Wilson, Nucl. Inst. Meth. Phys.
Res. B (2005), in press.

[123] A. E. Stuchbery and A. D. Davies, to be published.

[124] B. A. Brown, A. Etchegoyen, N. S. Godwin, W. D. M. Rae, W. A. Richter, W. E.
Ormand, E. K. Warburton, J. S. Winfield, L. Zhao, and C. H. Zimmerman,
Ozbash for Windows PC, Michigan State University, report no. MSU-NSCL
1289 (2004).

[125] B. H. Wildenthal, Prog. Part. Nucl. Phys. 1125 (1984).
[126] A. Poves and A. Zuker, Phys. Rep. 702235 (1981).
[127] S. Kahana, H. C. Lee, and C. K. Scott, Phys. Rev. 1802956 (1969).

[128] K.-H. Speidel, S. Schielke, J. Leske, J. Gerber, P. Maier-Kormor, S. J. Q. Robin-
son, Y. Y. Sharon, and L. Zamick, Phys. Lett. B 6322207 (2006).

[129] E. Stefanova, N. Benczer-Koller, G. Kumbartzki, et al., Phys. Rev. C 722014309
(2005).

[130] A. Gade, B. A. Brown, D. Bazin, et al., Evolution of the E(1/2)-E{3/2) en-
ergy spacing in odd-mass K, Cl, and P Isotopes for N=20~28, NSCL Preprint
MSUCL—1309 (2006).

[131] P. C. Zalm, A. Holthuizen, J. A. G. de Raedt, and G. van Middelkoop, Nucl.
Phys. A 3152133 (1979).

137

 

 

   

“11111111111ill]