Computer Physics Communications 11(1976)249—256 © North-Holland Publishing Company
QUASI-HOUND STATE WAVEFUNCTIONS R.M. DeVRIES Nuclear Structure Research Laboratory ~‘, University of Rochester, Rochester, New York 14627, USA Received 15 March 1976
PROGRAM SUMMARY Title ofprogram: QUASI-BOUND STATE WAVEFUNC— TIONS Catalogue number: ABMQ Program obtainable from: CPC Program Library, Queen’s University of Belfast, N. Ireland (see application form in this issue)
Computer: IBM 360/65; Installation: University of Rochester
Nature of the physical problems A program has been written which calculates the solution to the time independent Schrddinger equation for particles trapped by a Coulomb barrier but whose energy is greater than zero. Method For fixedofsolution Woods—Saxon radius and diffuseness, the depth of
Program language used: FORTRAN
the nuclear potential is iteratively found starting from a radius where the total potential energy exceeds the particle energy. Using this nuclear potential the quasi-bound wavefunction is generated from the boundary conditions at r = of
High speed store required: 175K (WATFIV compiler)
G(r)/r where G(r) is the irregular Coulomb wavefunction
No. of bits in a word: 64 (double precision) Overlay structure: None
Restriction on the complexity of the problem The program is dimensioned for a maximum of 400 mesh points and a maximum angular momentum of 10. Both these
No. of magnetic tapes required: None
limits are easily increased if desired.
Other peripherals used: Card reader, line printer Number of cards in cvmbined program and test deck: 628
Typical running time 212Po 208Pb + a), 7 sec using the WATFIV For the test case ( compiler — this includes compilation time.
Operating system: OS MVT Release 21.8
Card punching mode: IBM(029) EBCDIC Keywords: Nuclear, wave function, quasi-bound, unbound, Schrodingers equation, resonant state, reaction, Woods— Saxon. *
Supported by the National Science Foundation,
Unusual features of the program Because the solution is exactly on resonance no searching on the phase shift is required. Therefore this technique works even if the resonance widths are extremely narrow (1010 eV for example).
R.M. De Vries/Quasi-bound state wavefunctions
250
LONG WRITE-UP 1. Introduction
The well-known solution for large r (VN ~ 0) is
In nuclear reactions, states are often excited, which may be considered to consist of a core nucleus plus a single nucleon or cluster (an alpha particle, for example). Fig. 1 exemplifies the radial potentials involved; the orbital angular momentum is ignored in fig. 1, but is easily included. Discrete states at various energies En may exist as shown. The states are classifled as ‘bound’ if En <0 (in fig. 1), ‘quasi-bound’ if 0Eb. The solution of the radial Schrodinger equation for each of the binding conditions requires somewhat different techniques. The method of solution for bound states is well known [11. Unbound states and related reaction mechanism calculations have also been studied [2]. The quasi-bound state solutions are usually found by searching on the phase shift (discussed in section 2) and requiring that it go through iij23. A different technique is employed here which is particularly useful if the state involved is very narrow.
~I•~Sifl (kr + 6), (2) where 6 is the ‘phase shift’ which for a sharp state with given L, and Ea~is 6 = fir. The term ‘resonance’ comes from scattering studies. If Coulomb forces are included, the solution to the Schrodinger equation (again for larger) isa linear combination of the regular and irregular Coulomb functions which look like, respectively, / zZe2 F sin ~kr in (kr) +
2. Theory
— —--
and G
cos IIkr
+
VN(r)) i~i(r)= En ~(r).
.
-~
0
(1)
Eb
tions 0
/
/
I
I
R~
R0
r
-~
/
/ i’Z__—~i VNucIear + VCouIOmb
Lii
(3)
v
corresponds to a resonance situation is required. For relative angular momentum zero and in the absence of Coulomb forces, the Schrödinger equation becomes m V2
zZe —h---in (kr) + ~2
Clearly our quasi-bound resonant state wavefunction must behave like G(r) at large radii [boundary condition (a)] The other boundary conditions are: (b) ri~fi-÷ 0 as r 0, and (c) i~tishall have a specified number of internal nodes. It is not possible to satisfy all these conditions with any arbitrary VN(r). If VN(r) is taken to have a Woods—Saxon form
The solution to the Schrodinger equation which
/
—
______
VN(r)
R’/ 1 (4) exp 1~r , then for fixed radius R and diffuseness a, it will be =
~
“
‘
—
tions are satisfied. Of course, could vary any necessary to vary V0 theone three boundary condithe other constants inuntil the problem rather than V of 0. The problem then becomes how to start from some initial (guess) value of V0 and iterate to the final (unique) value for which the three boundary condiare satisfied. One common technique is to start the numerical solution to the Schrodinger equation at the origin, vary the depth V0 until the proper number of (internal) nodes are obtained inside of a selected surface radius, and then integrate on out to large r. V0 is then varied in smaller steps until 6 = 1/2 or equivalently the ratio of the wavefunction in the intenor to exterior (large r) is maximized [31 Vincent [2] starts at large r withG(r) and integrates inward. Likewise an outward mtegration starting with condition (b) above is compared with the inward integra.
Fig. 1. Energetics of a (1 = 0) quasi-bound resonant state systern as a function of radius between the core and particle (cluster).
.
.
.
R.M. Dc Vries/Quasi-bound state wavefunctions
tion
at some matching radius. The difference between the two solutions can be related to a change in V which starts the next iteration. This method may also be used for unbound states. A related but different method is used here. Notice in fig. 1 that between R1 and R0, En < Vtotal, ~ that by starting an inward integration in this region, the usual methods of solving the Schrödinger equa-
Ic’
0 2
6
~
C
~
~
82
251
~ ~
tionV0 foruntil bound states may employed, varying conditions (b)be and (c) abovethereby are satisfied.
~
The problem with this technique is what to use for
I0~~ 116
I
starting values since the usual bound state starting (r~a,~) values are [1] i,D(rmax) exp(—krm~), (5)
10-s Normal (~‘-e~) Starting Ca,ditions
with k = \/~)~z2 and p the reduced mass. Our positive En energy rules out these starting values; however, it is demonstrated in fig. 2 that almost any starting values may in fact be used and only
I0~
Arbitrary starti~gI Conditions IO~
2
4
6
8
0
2
4
6
8
20
the region of only ~tiwilltobe in error. are outermost using this technique obtain V Since we 0 at the ongin, we can use purely arbitrary starting conditions. Obviously one should start close to R0 (rather than
R(fm) 51Cr = ~ v + ~i (E~= — 5.0 Fig. 2. Bound state system of MeV) comparing integration results using normal (asymptot ically correct) and arbitrary (~i(rniax)= 0.97 X 10—30, ~(rmax — 0.1 fm) = 1.0 X 10—30) starting conditions. A few fermis inside of rmax the two methods yield identical results. Each was normalized to unity,
Ri). Having established V0 (and, therefore, VN) we now perform a new integration inwards from Rmax (well beyond R0 in fig. I) using initial values of G(r). The result is our quasi-bound state wavefunction. It is found that the first wavefunction (satisfying condi-
0~
I0~ —
~
I0~ 212
o~
Po (g. s.) = 208 Pb + a
0
(‘.1 -~-
jot
_____
—
100
poshive negative
-—
Ic-I
2 icI0~
“ 4
8
12
6
20
24
28
—
32
‘
%i”T\I~”YfliT\r~~Y. \[ 36
40
44
48
52
56
R208pb+~(fm) Fig. 3. Quasi-bound resonant
state wavefunction normalized to G(r)/r for large r.
60
252
R.M. De Vries/Quasi-bound state wavefunctions
tions (b) and (c) but with arbitrary starting values at rmax) is in perfect registry with the final result except, of course, for the arbitrary normalization of the for-
Table 1 Input Card No. Parameter Description
men wavefunction. A typical final result is shown in fig. 3. The overall normalization of i,li has two different forms in this program: (a) iP(r) = G(r)/r for large r, and
1
2
R~
(b)
f
i,1i2(r) r2 dr = 1.
3
0
The R~cutoff radius is taken as the first exterior node (see fig. 3, for an example), but any radius in this region yields the same normalization if En is sufficiently below Eb. The program outputs ui normalized in both of these ways.
ICON(12) control integers — not used — see section 3
1211
ALPHA DRF RMAX ALVAL ECM
run label step size maximum radius I (angular momentum) En (center-of-mass energy)
8X,13A4 9F8.4
FM
particle (or cluster) mass (in AMU) particle (or cluster) atomic number
Z FMA ZA RY
3. Program 4
Besides a main program there are four subroutines, The main program handles input of control integers (unused, but available for individual disk writing erations), and other (but not all) input data. It also
FZ
~
VR
processes the normalization (b) above. 3.1. Subroutine FORMF
RY ARI VSOR
This subroutine is a modified version of a subroutine in the program DWUCK [4] Two cards of input parameters are read in first. The Coulomb potential and VN(r) are generated and then the I (angular momentum),j (total angular momentum) and s (spin of the particle or cluster) values are read. These potentials and the calculated rm~value are used to find the required V 0 value by calling subroutine BIND. The step size in this search is fixed at 0.1 fm. FORMF is called (by MAIN) again to generate VN(r) with the correct V0 and using Rmax and the chosen step size. .
3.2. Subroutine BIND
This routine is taken from the program DWUCK [4] As discussed in section 2, it determines the correct V0 by integration until boundary conditions (b) and (c) are met. The method is that of ref. [1].
Format
__________________________________________________
5
FN FL
FJ2
9F8.4
core nucleus mass core nucleus atomic number Coulomb radius Rc=RYFMA’~’3if RY<0 Rc=RY.(FMA~’3+FM”3) if RY >0 type of potential (only volume Woods—Saxon wells included, but FZ eases user modification) V 0 (guess value) 3,ifRY = VR*VTRIAL <0 R = RY*(FMAI /~ + FML’3) R= RY*FMAi~’ if RY > 0 a (diffuseness) Thomas—Fermi spin orbit strength number of radial nodes
9F8.4
(not includingr = 0, °“) angular momentum, 1, of the state — same as card ALVAL, card 2 twice the total angular momentum of the state
FSS
twice the spin of the particle
VTRIAL
first guess for V0 — if 0., VTRIAL = 60. is used
Note: The input of this program follows that used in the Distorted Wave Born Approximation programs DWUCK (P.D. Kunz, unpublished) and LOLA (R.M. DeVries, unpublished).
R.M. De Vries/Quasi-bound state wavefunctions
253
Table 2 Error messages Error message
Subroutine
Cause
RMAX CHANGED TO
MAIN
Rmax < R
0 linputas> 10
Rmax STOP
FORMF BIND
E0 ~ after 16 tries search for V0 (correct) has not converged
STOP STOP
MAXIMUM L VALUE 15 10 YOUR STATE IS UNBOUND OR NEARLY SO-SORRY BUT WE CANNOT HELP YOU BOUND STATE SEARCH FAILS TO CONVERGE IN 16 INTERATIONS
Note: if Rmax/DRF > 400, Rinax
=
Result =
> 400*DRF
> 400DRF.
3.3. Subroutine COU
to one table 2.
This routine was written by Wills [5] and calculates the regular and irregular Coulomb wavefunctions at Rm~and R~~j~DRF when DRF is the step size.
‘.
The program error messages are listed in
3.4. Subroutine OUTER The final resulting quasi-bound wavefunction is generated from starting values of G(r)/r which are obtained from subroutine COU. The method of integration is Fox—Goodwin.
6. Description of the test run = 2o8P~~ + a is calculated. The final result is shown 212Poin The wavefunction appropriate to the state fig. 3 and was used to determine absolute alpha reduced widths from alpha decay and alpha transfer reactions [6] The program has also been tested for proton states.
4. Program notes
Acknowledgements
The program must be run with at least forty bit words for sufficient accuracy. If a resonance profile is required (i.e. the resonance width is non-negligible) the R max starting conditions are easily changed to F(r) sin 4, + G(r) cos 4, since subroutine COU calculates both G(r) and F(r). The additional phase 4, can also be used to make certain that the ratio of ui(intenor) to 4,(exterior) is maximum, as required by a resonance.
C.M. Vincent, I. Towner and H.T. Fortune all provided valuable help to this study.
5. Input/output The input cards are listed in table 1. Normal output is illustrated at the end of this paper. The list of ~ values are the interations in subroutine BIND attempting to find V 0. The inverse of the integration of 4, per normalization (b) above is labeled ‘normalization .~
.
References Eli
S.M. Perez, Report No. 38, Nuclear Physics Theoretical Group, Oxford University; see also ref. [41. [21 C.M. Vincent,.private comrnunication; see also H.T. Fortune and C.M. Vincent, Phys. Rev. 185 (1969) 1401. [3) D.H. Youngblood, R.L. Kozub, R.A. Kenefick and J.C. Hiebert, Phys. Rev. C2 (1970) 477. [41 P.D. Kunz, program DWUCK, unpublished. [5] J .G. Wills, Thesis, University of Indiana, unpublished; see also J.G. Wills, J. Comput. Phys. 8 (1971) 162. [61 R.M. DeVries, Shapi~a,W.G. avies,Lett. G.C.35 Ball, J.S. 835; Forster and W. D. McLatcl~ie, Phys.DRev. (1975) W.G. Davies, R.M. DeVties, G.C. Ball, J.S. Forster, R.E. Warner, D. Shapira and J. Take (to be published).
R.M. DeVries/Quasi-bound state Wavefunct ions
254
N 10
16 * 16 16 16 16 * 16 16
N
16
0 a 0 0~ 0
II
16 16 16 16 * 16 16 16 16 16
N II N 1.1 .4 0
0
U)
100 00
trIo
16 16
0
o o
*
0 0 .‘.
0
II
•
16
O
~
*on~on *
*05 5-lID 0 I-IN
5-4
5-
*
N _ .4
16 * * 000 16 ~ 00 00 * ~ 00 *
ZN =0 04 0 • 0
O
* 0
~0
I
1.-I us .4 0’
•
U)
54 * N 1!) on on
U)
1
•
co 0~O
5-1.4 .4~ on
5.4~4
.4 0.0 0
p4..4.
.
N
II tt 04 US
*0.404 *N~b1 16 16 16
5- 5-
In ~
~0 a
It H
nuscnvv
000000000
10 N ~D 5- 0 ~ 0 I’S ID N Ill ID I’, N NO IDmO,onotnonoO
•
.
a
.-0 I
~~N.16.OOOO III
It It
II II II It II It II It It
.4
on 0
1414F~14 NNNNN 000000000
05* 10 5- ~ ~ N
00 00 00 00
00 00 00 00
. . U, 5-
.~‘.
~ UE’t
onon
.
It II
.
It II U) .45)
r U
5l~
0
10~Z on on 5-4 .4 on .4 04 N .4*
16
0.
.4
.4
0 0 5.1
.~
0I’5Q .4 54 NO
.4.1-410*0 405r.1*
0 ~
55 05-0
010.4*54
.
00 100
5.4 .4
U
•
00
1~CN 5-N . .
on
U 0. 16 0 .-IN *0. rn54’*Z
•
5- 5-
0 N
16 16 16
*5-4
•
0 0’
*0 0
ID
•
Z**Z**Z**
5” *1110 H.~ 0
l.a.
•
,.~.-
16 .16~’1
*
16 0’*.4 H 54 N 55 N• 16 16 54. .4
•
m.-,-v-,-’-.-.
II II II II II II II II II
I 5-
*~us
on N
5-)t00
16
* 16
H40
monO on onr—o ID
‘-
16 a 16 on 16 16 16
000000000 000000000 000000000 000000000
NNNND3NNN
5-
U)
54
**t*m*~~*
II II
NN
U H
U
•
.4 on 0
* ~N * II It II * 5-4 *01!) *0.45.
•
II II II II II II II II II
,
*5-U)N *I0~
1.
00 00 00 100 . S 00
0
• U)
54
I—
onus
Z 0 H 5’I .4
H
z
II II
p4 N N l~ • *
•
S
.4
*N00 16 •
U,
0
•
ono
*000
ID
I—
~0 0
16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 * * 16 * 16 16 16 000 * 000
000000000 000000000 000.000000 ~ ID ID N N F- F- P- F-
.4 5-4 .4
5.40
~.
on
5-4 on on
54~ on 0.54
0
*
NO 4’) 10 1!) N O~ - ~ U’) ~ 4’) N 0’ U) 10 IfS N I’S N (‘~N
N 10 4’) N N
.
I!) 10 0~ N N ‘.0 10 5N N N N e —‘ i-r-
‘.0 N N ,-
It II II II II II It tI II
R.M. DeVries/Quasi-bound state wavefunctions
255 5-s.
—
N
5-.5-5-
—
o00000000oeeooo000 000000000000000000 1010 sO 10l0 so 5010.0 sO sOsO 10111101010 on 000 5- N N 00000000000 N 5-5-5-0000000—5-5- N 5-5IItIIIIIIIIIII 000100010010101000000101001010
N
oeNeonep0OeN.O.VSNm.-*000105-5-*105-o.4onoMNoVS5-a.O.oO.onU)rnO.
on *0 VS IfS N N0 050*5-I’) II) 005-5-001* N 1115-00110
N 5-.- N 1~45-5-- N 010 N 00 V~IN on N 05.4.40 N 514 ‘5010 I’. I’) Oi5- 5-11)100*010* N II’) 05 I’I* 01O NO SO 005-105-5~-N 5- on IOU) N N 0105-5- N N 5-10 0,,-1005-5-5-.N5-Inp0usmop0U)NN5-5-o5-0VS105-01P-N10N.4NO1001NOSON5-onNVSN5-10010111 •VS.s..nalN.5-.-eI.45-No5-0O.U)NNN01O10rnVSVS10NrnON10USIAN5-ONIC5-O10oNI115-ON r.N4*5-I’~5-.4410 at,son55-ons.p05-#5-10on,-0suIonN.-,-5- 010 .NIn—on* * 10 one. 5’ essononf-i • • • • S • • • 00000000000000
I
I
I
•
I
•
• • S • • • •S*•~••~S e. 00000000000000000000000000000000000
•S~IS~.S•.
•
• •
I
S •
•••
•
••I•S•
1111111
N 10 1010111 1010101010 SOlO lb VS 5111010 on on N N N 5-5--
00000005-5-.00000000000000000000000000000000000000000000000000
N
5-.-
5-5-5-5-5-N
5-5-5-—
P1
111111 liii 11111111 0101000000000001010010001010101000100100001 01001010100011001010100000101 105-NO051001N0000N10NION*N10NO5-N5-N01101010N01100m1005-snOr.NN*OION 005-5-5-010 IA NO.4.41001 N NO.5-- 1010 N on 5-015-ON In 055- N IA 1000 NO .-e.~s4 US N — N N05- 00 eOISON10N.05-onInu’)5-0505r45-NeoNer-01ON.-on*s0005-01NOOUSmN105-O.5-o5-ONNO 10 N 001 N 015- US 5-se on In 5-10 N 1000 IA IA 5-I N 05000 P.* 5- .- P.. P1105-01 IA N N 05-10100010110 on 00 ‘0 N VS 05055-N IA In N 5-rn N 5-10 N C4 5-0*01105-5-1010*105-5-10 N N 5- N 100.405-5-10 N 0101001011001 1001Np0N10ms.10N5-05-05-10N10N01O.mNO.01N1005-o1010N01’55-0~sOP4e.-N10105-N.0 t4e. 5-5-N 15-5-on on on .on 5-15-5-105-50 N 505-01105-N ‘5-5-050 10 N05-N* ~ 105-5-ON on 5- l’410 • S • S • S • •tS••S•• • 00000000000000000000000000000000000000006000000000 I I I I I
S
I
•......i
S
1114111
.o us us 10.4 on on ~i N.- 5-5-5-0000,005-5-5-5-5N 5~ 5-~ P45-5-~5- N ~8~~eeooeeeooeooeoooooeeoeoeeoeooeooeoeoeo .0100000000000000100001000001010101001000 ~l)t&I!) ~ ON 415-N SO N 5-000 N 5-105-01 N 01.01 5-’~ 01100015-N .‘5-05O5-055-0Ifl1010 N 5-05-N 505-5-5-010505-01 N 5-01005-N N 5-VS N N 010’5-10 010100110*0100 N 105-5-010 IA 5-0505-- l~00*5-1010 VS 5-5105-0005-10 OUSNPI NNOu, 11505- 10 00 100101010 ‘0.11 ON 10100. N11110 511*5-10 l*mONNOOIOonOS NNN*rn.VS5-01*015-51505U50.”5-VS1010ON1010N5-N5-5-N5-10N10N5-N10VSl~VSNNso05-10*10 01010005-ONOO5-Q5-NNSOON105-011001050NNU’)0.-0N105-0100001a5-5-105-105-5010510105-5-0*01 N N 00 N 10*rn’10001 N5- 105-100105-05100105-NO 00001001 I~5-N10 05-10.105OVS5-5-N5-NNN10.-105-*5-105-N5-5-Non5-5-105-N5-5-.010U55-5-01N5-Ø’*I’5.-111Nonon14N S S S • •S,...S.SS.S . 00000000000000000000000000000000000000006000000000 I P I I I Ililill N N N 101010101010101010101010*105-5-N N (45-55-5-000 0~0 05-5-5-5-5N.- 5-5-5-5-.’. N 5-5-5N
A~
00000000000000000000000000000000000000000000000000 00000001001
000100100000
00
I
I
I I
I
I I
00
050
I II I I I I I I I I I ~ 0001 00010001 5~5105-1001101005-0 N 010 N 1005-010101 ~ N 50N .4*105-10
01 01 0 0 01~0 01 01 00 0 0001 0 N 5-5-105005-10 N N N 1000 101051 N 5-055-N *101001001 NON105N 1010100 N 5-1001010 N N N 5-N N N ~0 use us N l~0100011001001 N NO.4 N 5-015-5- *0*015-ON 050’ N 01105-5-10000050105- N10 05 N InN *010005,01005-N 010105-10
N 1010015-N 105-5-005015-0105-*10015-10U510m10*
011005-N s 100 N 10105-5-010 *054fl10105-5-100N0m5-0101O.10*N105-010105-01O10,105-5-O.*mm10010105-O.01ONUSN5-5 O.N5-05-5-100N5-rnIANNIANUS105-1010100O.05-N5-N01105-00N5-I’1,010 OIS~5105-105-5-5-In10In InN 01010 N01*sN 10 5115-105-10 P405-.N 5-N 5-10* N 5-5-5-05 N 105-5-10 Nm *105-515-N 5-5-NO. 5S.SS.
S
00000000000000000000000000000000000000000000000000 I I I I I
I
N N N N 101010101010101010101010105-5-N N N 55-0000005-5-’ o 0000000000000000000000000000000000000000000000000
5- N
I
I
I
I
I
I
5-5-5-’5-5-5-5-5-5-5-
I I I I I I I I II I I I I I I I I 01001001010 0~0100101001000 NONVSO.00015-.5-100*10*NNNO01O.NN1000N5-O*5-N.-Qs015-N10N10*5-.1A1010N10.5o N 10010 NO.0 5-1010005-105-N 1050105-NO 0105-N 10105-0105-1010 ~5 N05- ~ 1001051010 P4005-
I
NO.010*N10105-5-In5-5-a5-015-10O5-5-NmN05O*O.O10010001NNOoao~105-01NNoIn5-N ,105-N5-5-0IOsONOUtN1000—05-5-5-N ONNNN0510NOeON10ONON10SO,.,.O.5-ON.OmN N5-0N100N005-101010.-OrnNUSN10N..us 10NN10rn5-.45-NO.N0O0’5-NN~Io105-oN1010010 — N 505-5- P41001101010 N 0.-N 100.4*100010 N N N 5-100’ 0105- P45-101010 N N N I1~1010011001 (‘410 N5ON 5-5-010 *P4 —5-101044105-105-10 NO.105-0 .4 P15-N 10 IN (‘45-5-0’ N U’) 5-5-10 (‘45- ~ *5N55-5P1500000000000000000000000000000000000000000000000000
I I I I I I 1111111 0NNN50101010’010501010101010105-5-5-NN5-5-5-5-00eO0O,-5-5-.-,~N.-5-5-5-.-,.N.-5-
51 ~
~
15 0 54 0
50 50
5-,~ 000000000 000000000 000 00 0000000000 • 00 I I 00p0 I I I I 00000 I I’ I I I 000000 t I I I I I oooooooooooooooooooo~ooooooo~c~ 0000000C)C~ 001010101000101 1000’505-5O5-055-NIANm5-5—(405-5-Nl’)N *m5-05-05ON*05(’-N10mO.5-*$NNION.4OSNO. C’ N 0100 N 10 N* NO.05-1005-100*5-S N 01054’) N 11105-SN N 050* N 100100 N 1010 N 10 NON 000’ N *00’ 511 (‘45-1054’) 5-011010105-N N U’) 10*10105-5-105-105-010 1-rn 5- * 5-10.4001010 0 5-NO 10 N 105-050 1001NOONNOUS5-10N*N01NN1000N1010N5-0*101050*505-Q5-10N10O10S.1OO10N01NN5-0 5-N10015-N100N5-NN50mN0O.N5-0N0NNO5-*O.50N50ew505-N.4*101-~qmN50oIn*o*.0O105-N100110NOO05-5-05-05-P~O1005-50*5-5050O5-10N10~N0ONe10*InN5-NN,-InO.* 5-P4’01NIn*5-NN*N1-IN50NNN5-10N01*N5~N*5-N’5-5-01N105-5-5-N5-**5-P15-.-5-5-N_
øo0oeOoeo~ocSoo~ooooooooeoooaoooooeoooooo0eooeooeee I I I I I I 0NNN5050505050501010s01$)in*4,5-5-5-NNN5-5-5-000eo O0O0OO0O00O0O000000000O0000o0eo0000o00e00oeoooeo~
III -
05-.5-5-’5-
I
Ills
‘5-5-5-5-5-N5-.5-
— 10
PIIJ_IIIIIIIIIIII’III
0
00.0InO.5-10*N1005NON5-0N5-5-r010O.1-10NNN5-0m0OrnNIn-00 01010000010000000000101000101001001010101O00101010101~.~010100001010101~rI~
10 50 0
0InO’0NO.005-d1010In0100O0IS.P4ooNooNmNN050m101-5-oNNP410(-,~m50NmmrnNo. ‘N 5-**U)01ONON105-105-N’05-005N5-010O10N*S~O0S’)N10P110.4100510NN10NQNIn’0O 0N10NNP.N5-10N 5-Inm5-5-N005-NNO.P1N05-*’N5-O5-10Npi,.N.4rn05*1041N5-..~NN05~5-_~~ 1000 P100110*0 N (115-5- III 0005-104’) 5-100’10050110* N ON 10105-1001 N 5- N (S~ ON 051510 P1101055-N 5NNIn50510 5-10100505-5-10*105-01-5-105-5-5- ONONO 010 *5-05N N10P1 ONO.0010 010I5~1005-5-InO5flO50
0
5-N5-5-ON’0IfS*I)5-*N0NNNO5-5-U’)N5-IflP1’.0*rnN5-5-05N1010NN...rn~.mP4*...m5-m
5.4 sI
50
o
10 N 5-
0
00000000000000000000000000000000000000000000000000
I I I I I I IIIIIII 0 1-I ON50%05050’010s0,010s01010101010m5-mNNN.—5-00000e 05-5-5-5N5-5-5-5-N5-5-.-.54 ~ tI,IIIIII.IIIIIIIII 50 0a0010000000O0010010001000001000000010010000000(~00001000100 0’ O*N5-’0*5-M0O.05-NNN*10101005-N50O*rn10.’._on05-10rn1010eP4*IA0IA10_10’0N’0P1_ 5.. In10NO.QO05-NN*105-OO10N5-ON1010e00NN_10105-NN50*5-.~050’0p~Nrn010NmP4~ 50 OMNN5-N*010*O.0s’.oN105-1010ON105-5-5-.-O.NOONON105-W5-*N’0OIAO_*O.*_O 5-ON01100N010N10105-In*5-.4010NN5-5-’0NOO5-5-5-P1 54 *5-ON5-*NUS*5-O10O050mO’0**N100O.5010m*0550eNN10N1010105-1051~5-O.NeNrnmO.o .4 N5-N001N*N’010N5-s~NP1Na10N10O.*NO5-(’4...No5-,0O_mNe*In,...4’)Qe_...eNNO.’0 50 N1050**O.10 50P1.,.~’0O’0O.5-OUS,.,.0 5-5-05 N 501010 US* *5-105-05 C’! OP4 015-5-10 P15-In (115-0105- P15-5-5-050 ~ (‘1 ,5 5*10 N105-14 0. • S O ~I 0’ I 1.1 I P 1111111 0 • 54 ‘0 mm 0’
01 10 5-
N 50 0 0
54 ISS
5-s-
OOOOOgO~e~e ~OgOgoO~00eOCeoeog$~gggo 000000000 00000 000000000 0000000000000000000000000000000000000j000010’c)000 0 000000000 O.NIn5-5-O.N105-5-01NIA5-s01NIA5-501NIn5-.01N’ln5-’.o.NIn5-.-05r..us5-,.~O.NInrn5-e.NInrn ‘a
S
.•.•S•••.
005-N5-10*1010N0005Os.NP45-*IA1010N001Oe..P1on1010VS10N0o01e..~4N5-*VSSO’0No05 ~5-5-5-5-5-~5-5-N N (‘I NN N N 545-1 N N N (‘45-5-I 545-5-5-5-5-5-5-5-5-
50
0
I—I 10
50 5.4 .4 50 S a
R.M, De Vries/Quasj-bound state wai’efunctions
256
5-5-5-5-5-5-5-5-5-5-5-5-N P15-5-10*10101010101050 N N N N N N N 0000100000000000005 500000000000000000000000000000000000000000000000000 III lip III I 111111111 I~I 1111111 III 111111 III I 111111 II 0001010000000000000000 ~ 01000000000000001000—010000000-00 5-50051- ‘N p0010050010015-5N N 5005-05* 10501001M0’OO SON P405-055-105- .01 N 5-100 OIlS N 05105-5-P1e50100505-5-50.u’)010.40505.4IAeIA5-.-N0VS105-s-01100.-1-N*.5-5-5-NlfllA,P15-IN 011005-001-N 1-00 05 Ne45-*00 05N5115-N 10N* P400-10 0100 5-5- 055- *.NN P11- I’~1- 101-50 P155-. 1-1010* P110010510100*0001 N 5-105-10 IAN 5-05055-5-05015- N 05105-001 VS ON 5-505- N 105-I N 5050 P15-N5-~5-*5-.N5-NIN5-10IAP1OP105105-5-(A50N05-05005-5-OVS105-(A0tOS’..100105-101001N ON10OS*055-1010N055-10’N015-010505-IA10005-0*O5-505-001-NNNONr5-IA5-5-10105-USO 5-5-5-5-101010101010*105-105-5-5-105-N15N. P5-N 5-N 1015-P4.-s..-. 01015-N 5-5-ID ID 105-5-5-5-1010 N N 05
SetS
S.
000000000000000000000000001000000000000000000000000 I I I I I 05-.-
1111111
N P45-5-10101010 US 10501050 N N N 5- N 1-1-000000000000000001
5-P4.-5-5-5-5-5-.-
00000000000000000000000000000000000000000000000000 I
III
I
PIll
II
1111111
0J11~11!5 ~ 0N5-01050P4L05-In055-N01N010&ANU’)1005-P1*.,.N10050f410
I
III
PIll
1111
I
III
I
01000i-105005-00r450*0
5-1-N ON 0 IAN 005 5- 05-5-5-10001 00 105- sON 105-10 100101-5-01-011055-lAOs. 05-10015-005-0 P101-IA 0 IA 5-1050101010101-5-1-5-5-010 N 05-5-5-N INN 001 NO NON N10’5005-Ifl5-0110010105-55-5-N NN5- 015-5-055-1005010015-5-ON 50010150015-10 ~ 010511.11010015-NO 115-ON N 5-01010 ~ 1110 5010I115-5-0N5005-5-.5-IA5-5-*10N*5-1000IA5-IA50r.0N5-*NN05N10105-SO05NO%0*050 *0 01.410 P15-N 0101101-01010010 IA 10105- N N 1- P45- N P4105- P11-5-05-01 IA 010505’ N 5- N 1005-5- N 00 — N5-~ *5- P45-10105-105-105.10—10(1105-5-010 N N 105-N 5-5-0105-5-5-5-111111105-5-5-5-1010-N IA . 00000000000000000000000000000000000000000000000000 I I I I I I I~I I I I I 05505-5-5-105-5-5-5-5’. 5- C’! P15-5-10 101010101010101010000000000000000000000000000000000000000900 I I I I I I I I I I I I I I I I I I P I
N N N 1-5-1- 00000
0~0 0005
00101 I P I I I p I II , I 010101-10105-100 P41-01010 (‘41)50105-105001010105-0 P4101-SO 1-rn10505-500SO5-05-N
I
1041105 P15-10 N N1000055-5-Vl5-MVS101-5-*NO0IA510105-IA*N....5-5-1005100101-~450...01N*ON05-5-SO5-05 105-105-0510100105-N.’N1001101-105-5-0105105-Cfl5-0150 m05055-*05105a01-5-100 01’5- 5010 IA 010010 NUSO1-05501-0%1000150NN5-P4010105-005-P45-0N1005-P1105-5-10005-(AO5-0505-N0005P15-100500’ 1005 N (‘1 IAN INN 505-1000 P41001-5-05105-05-101-005-01010101-5-1-105-5-01-0100550 00(11 01- N10 051.545 NO.5010 P41-5-10 1- 5-001 105-5-01 015- 5-1-5-5- 10P405-5-5-O1*0 0~ 015-5-001 105-5-105-5-5-5P45-N 5115-5115-105-10 N 01105-010 P45- 1-105-N5-5-50105-5-5-5-IA IA 111105-N P410105-0 •0
0010000000100000000000000000000000000000000000000000 000010 1111111 000000100000 I I P I I 00000000000000000000000000000000000000000000000000 00005-” IIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIII,III1III4III 5- 5-NIN5-5-N P45-5-1010*1010501050101-N1-1-1-1-1-0 0000010000000010001000000000010000000000 01000’0000010000 000.5-5101- 5- 00010 P110 1-105- 00 5-05 1010 5-1-OP1*5- 101010 05-0N51005- ,..0IflSO 105-105-10 105-05*5-10 5-5-0NN50*505005N*10P105001011-1-NN105-0105-010105-10505-*P15-~~~ 1-N015-10055-5-N0N1050011010N’0**50050~105-10P1U’)OIA5-IA105-10105-10... 10115-5-N 005- 105-OUIP1 10 011- 0510 105-10 05- 5-10 P15-N P1545001-0 015-10 5- *1- 00 01-5-1111010 10005005- 10 100100 N P105-100501-5-505p4 010005-1-5-005-105001105P15-00011-505-105-1050 P1105-0 P4005 101-105-1-N05- 051- P150 P11005005-1-0500001101010 P110100105- 01105-505-5-01-10 N “5-1-C’4* 5-5N (‘15-5-1-105-5-105-0’ (AN 505-50 (‘41- P15-10 P40510 P45-0105-N .-r5-050105-05-*101010 P410 P410105-5-
00000000000000000000000000000000000000000000000000 I
•
I
1111111
5-00005-5-55-P4P15-,”P1P15-5-1010*1014t105050101-NNN1-1-N0000010 500000000000 000000000000000000000000000000000000000O 00005-00000 01010 11111111 III N0N0P1010100110105-0015-00N010010N100501-10N10100050501010.1-NP110NN50100110505-5IlItIPlI 5-N05-O.1O.5-Q5-N10105-P100N5-1-(A0015-0IA05-501-00010050100P1N105-*01000U’)’05-O.0550 0501-10©5-50500100110N05-5-**P41111050P110*1005-N011005-5-5-P15-1-50105-1001 10001 010 1010ONN5-1000ItS*05-051010100P1P1NP1*5-0,”010~N5-1001101001N1010IA1-00101001010O.50O P105-010051010 ‘0OIA1-050*10010105-5-0N5-05100550505001501001N05050510N1-NON101-0 010c501-I.05-5-051050P405P15-0P15-P105-Nl11011010100105-5-01001P40501050101050100*10* 5-5-N5 5-1-105-5-100 InNS.. P150 P105-5-In P15-In P15-0105-P15-55-0550105-10 P1.10101010 NO. P110 * 5- ~ 0~
00000000000000000000000000000000000000000000000000 I I I I I I
50 50 0 0
5-i 01
N 10
1111111
5-000055-5- 5- 5-5- 5-5-P1N 5-5- 10 *10 10101010 50101- 1-1- 1-1-1-N 0000 050000 000% 000000000000000000000000000000000000000000000 1_I1IIIIIIIIIIIIIIIIIIIIIII-IIIIIIIIIIIII__II.III 01010000100101010000000001001010 0010100001000 01000100001000,,.00000 00501010 5010 5-,P1 01100 P11001 P45-50015005-1-1014)10105-501-1-N5-P405-N5-5-5-10010N1-5-5-1005P1111*05-5-50*5105-~5-0105-5010P40500IA10O50 InP1m*C’ 5-01A0N50010100*50*0 505-U550rn5-5-101005-1010101-OO.NNNO’0105-005-01005-0” 1000U’)5-5-1-5-ON0-5-5151-01010 1-100’O.10P1051-*10105-50In5’~5-N05P1In05-N05-1001051010105-0050O5-e.O0OC’N~”10O10O. 0 05 (‘11-1010 P11-055-05 P1 N P15001501010 P1 N N010 (‘4100105-5-10N0lnIn5-050510,0105105-50Ifl5-0050U500500~I15P1 105-05-,-10*1010501-P1N05*50N 5-O*Q%5-5-1000105-5010N505-1-10N1010545*10015-10005O.5-
5-5-
•
o0oo00e0o00000000000000000000000000000000000000000 5-rnN5-5-C’N10.4NP110P10N0P1O.m5-10P15-105-5-C’In5-N5-5-5-011-10P4P1N10101010P4NP15-* I I I I I I 11111111 _Q0O05-5-5-,-5-5-5-,~.5-P15-5-5-**10101050~10501-N1-NNN00O0000O0OO05OOOOO 00000000000 0000 000000 0005,t0000000 0 00005-0000000 I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I
I
I
I
I
00000 I I I P
I
m 5-N_10*01-05-NN0’05-104r0050N015-N05N5-0110105-5-N1-N5-5050O.05-1005-P101105-N
5.4
50 .4 O .4 o 50 N 0 I-I 14
U N 0’ N 54 5-. 4 01 N 50 0’ 0 .01 I I-I 5) 0
50P40N01001100P11-01-P10505-0*N5-5-5-5-051-5-5-O.ON1-05-011001-00-105-NP15-5-’)1005*10 5-5-01105-*(”! N’,.P1rn 5050 05-OO.5-IAO’010*0101*mO0100100101-10*N100mN0105-010010101-001-1010505-01015-*10 *5-04SS 0105-N NI’SP1N0N01-010N 1-105- sO5-P5-Nl111-001001N1000 5-O*U’)_Om’0m_5-N5-10055-0**010000*010105-10105-105001005N0N015010 5-5-5-5-5-011-101010 m,om5m055-.’5-5-sON5505-5-5-50* P1P45-50’1-IAN05N10USUs 105-US 5-5-105-IN 000P110O.0m01105-100*505-5-N05N10NI’4*1005-105-0100P1000*’0010N00 0N050I15P10105-00N00
00000000000000000000000000000000000000000000000000 I
I
I
I
I
I
1111
III
5-005-5-5-5-5-55-5-5-5-NP1mm1010 101010501010101-N1-5-1-NN000050000000100000 000000000000000000000000000000000000000000000000 11111111 III I tlIIIIIIIIlIIIIl 11111111 III I 1111111
O 5-~ P105O.50mN5-*OO505-5-O_O0101010P101P15-N10505-05-00*550*1-10N10P15-~5-1-~4P1OOP1 N10P1_*P41-5-*1-’0e5-P1P1100015-1010*1-5-P41010N101-01-105-0050550015-05N5-10100 10051000015-*.5-0N5-mN0510P1101010015-N01’405-5-P11-P150101001051005-*105-0510500150N O5P15-0O0101050105-05-5-00501O5-0500110010m5-5-100010001-**5010105005-O*5-10~IA* 5-5-00* 01001010N015-10P101*05000101-0P150 1010*05— P15-0501-1-505010*105-5-5-5-5-5-105-1-5-5-105-5-5-5010 P1 P15-5-5-1-10 ~ 5-5- ~ —510 N 05 P15-05 * 05000 P15-1010 P11050100500 P4505-015-5-0’ P41000,” N 0*0*0 P1* P1101010105P~1-5-1-5-1-1-100 0 S Ir50P15-5-100105-5-01100P4111010500000000000000000000000000000000000000000000000000 I
I-
I
I
I
I
02000000002000000
ooeeoeeeog~oooooo~gS~8 000000000
1111111
~ 000000 oee 00o000000000000000 5-011-105-5-0’ 1- U’) 5-5-011-105-5-051-105-5-011-105-5-051-105-5-055-105-5-011IA 5-5-011- US 5- 5-~051-105— 5-5-5-
5-l’5-5-5-5-
5-P1NNP4P1NP1NNNP1f4P4rn
Il5-5-5-5-
5-5-5-5-5-I5