- Email: [email protected]
Hauser-feshbach nuclear scattering subroutine liana
C-20 COMPUTER PHYSICS COMMUNICATIONS 1 (1969) 181-190. NORTH-HOLLAND PUBLISHING COMP.. AMSTERDAM
HAUSER-FESHBACH
NUCLEAR
SCATTERING
SUBROUTINE
LIANA
W i l l i a m R. SMITH P h y s i c s Department, Trinity University, San Antonio, Texas, USA Received 22 October 1969
ADAPTATION
SUMMARY
Title o f manuscript: Adaptation of Subroutine SCAT for Use with P r o g r a m LIANA Adaptors' names and a d d r e s s e s : William R. Smith, Physics Department. Trinity University, San Antonio. Texas. USA Title of adaptation (32 c h a r a c t e r s maximum): ADAPT SCAT TO LIANA Adaptation number: 0002 R e f e r e n c e to original program: Catalog~te number: ACQF; Title: SCAT; R e f . in CPC : 1 (1969) 106 Authors of original program: William R. Smith No. of cards required to effect adaptation (including directive cards): 13
PROGRAM
SUMMARY
Title o f manuscript: H a u s e r - F e s h b a c h Nuclear Scattering Subroutine LIANA Authors' names and addresses: William R. Smith, Physics Department, Trinity University, San Antonio, Texas, USA Title o f program (32 c h a r a c t e r s maximum): LIANA Catalogue numbe~. ABKA Computer f o r which the program is designed and others upon which it is operable Computer: CDC 6400, CDC 6600. Installation: Computer Knowledge Corporation, San Antonio, Texas; r e s p . University of Texas, Austin, Texas, USA Operation s y s t e m or monitor under which the program is executed: SCOPE 3.1.6A P r o g r a m m i n g language used: FORTRAN High speed store required: 23,000 words. No. of bits p e r word : 60 Is the program overlaid? No No. of magnetic tapes required: None What other peripherals are used? None No. of cards in combined program and test deck: 571 CPC L i b r a r y subprograms used Catalogue number: ACQF; Title: SCAT; Ref. in CPC: 1 (1969) 106 Adaptation number: 0002; Adaptation title: ADAPT SCAT TO LIANA; Ref. in CPC: i (1969) 181. K e y w o r d s d e s c r i p t i v e o f problem and method o f solution: Nuclear, Reaction, Compound Nucleus, Compound-Elastic, H a u s e r - F e s h b a c h , Differential C r o s s Section, Integrated C r o s s Section, Optical Model, Woods-Saxon, SpinOrbit, Penetrabilities, Width Fluctuation Correction.
C-21 182
W. R. S M I T H
Nature o f physical problem S u b r o u t i n e LIANA c o m p u t e s the s c a t t e r i n g a n g u l a r d i s t r i b u t i o n and i n t e g r a t e d c r o s s s e c t i o n for any twop a r t i c l e n u c l e a r r e a c t i o n which p r o c e e d s t h r o u g h a s t a t i s t i c a l d i s t r i b u t i o n of c o m p o u n d n u c l e u s s t a t e s . T h e s u b r o u t i n e m a y be i n c o r p o r a t e d into o t h e r p r o g r a m s o r e l s e r u n with the d r i v e r p r o g r a m p r o v i d e d . T h e l a t t e r c a l l s s u b r o u t i n e SCAT I1] to obtain the n e c e s s a r y optical m o d e l n u c l e a r p e n e t r a b i l i t i e s . T h e input r e q u i r e d c o n s i s t s of the e n e r g i e s , c h a r g e s , m a s s e s , s p i n s , p a r itics, and optical m o d e l p a r a m e t e r s for all t w o - p a r t icle c h a n n e l s t h r o u g h which the c o m p o u n d n u c l e u s c a n decay.
Method o f solulion T h e H a u s e r - F e s h b a c h t h e o r y [21 g e n e r a l i z e d [31 to include s p i n - o r b i t coupling is e m p l o y e d . E i t h e r the plain H a u s c r - F e s h b a c h m e t h o d i m p l y i n g i n d e p e n d e n c e of c o m p o u n d n u c l e u s f o r m a t i o n and d e c a y or the m o d i fied H a u s c r - F e s h b a c h p r o c e d u r e i n c l u d i n g a c o r r e c tion 141 for an a s s u m e d P o r t e r - T h o m a s d i s t r i b u t i o n of c o m p o u n d n u c l e u s r e s o n a n c e widths [51 m a y be u t i l i z e d .
Resl~qctions on the complexity o f the problem Only p a r t i c l e c h a n n e l s can be t r e a t e d ; g a m m a - r a y c h a n n e l s a r e e x c l u d e d . N o r can c h a n n e l s c o n t a i n i n g t h r e e or m o r e s e p a r a t e p a r t i c l e s be handled. In any c h a n n e l the s p i n - o r b i t p a r t i c l e can h a v e s p i n up to 5
and the o t h e r p a r t i c l e m a y h a v e a n y s p i n . T h e m a x i m u m n u m b e r of p a r t i a l w a v e s n e c e s s a r y in any one c h a n n e l c a n be no g r e a t e r t h a n 31. A total of t e n t y p e s of two p a r t i c l e p a r t i t i o n s , and a total of 600 c h a n n e l s . m a y be s i m u l t a n e o u s l y t r e a t e d . T h e r e s t r i c t i o n s on m a x i m u m s p i n . n u m b e r of p a r t i a l w a v e s , n u m b e r of p a r t i c l e t y p e s , and n u m b e r of c h a n n e l s a r e due s o l e l y to s t o r a g e a l l o t m e n t s and m a y be e a s i l y r e l a x e d .
Typical running time T h e r u n n i n g t i m e is r o u g h l y p r o p o r t i o n a l to t h e n u m b e r of p e n e t r a b i l i t i e s c a l c u l a t e d , s i n c e the r e s t of the c a l c u l a t i o n is s h o r t by c o m p a r i s o n . T h e e a l e u l a t i o n a l t i m e s p e n t by LIANA a l o n e i n c r e a s e s a s t h e v a l u e s of the s p i n . n u m b e r of c h a n n e l s , and n u m b e r of p a r t i a l w a v e s i n c r e a s e . T h e t e s t c a s e r e q u i r e s 13 s e c o n d s to c o m p i l e and 3 s e c o n d s to r u n on a CDC 6400 c o m p u t e r .
References [1] W. R. Smith, C o m p u t e r P h y s . C o m m u n . 1 (1969) 106. [2] W. H a u s e r and H. F e s h b a c h , P h y s . Rev. 87 (1952) 366. [3] G. R. S a t e h l c r . P r o c . P h y s . Soe. (London) A68 (1955) 1041. [4] P. M o l d a u e r , P h y s . Rev. 123 (1961) 968. [5] C. E. P o r t e r and R. G. T h o m a s . P h y s . Rev. 104 {1956) 483.
HAUSER-FESHBACH
NUCLEAR
SCATTERING
SUBROUTINE
LIANA
W i l l i a m R. SMITH P h y s i c s Department, Trinity University, San Antonio, Texas, USA Received 22 October 1969
ADAPTATION
SUMMARY
Title o f manuscript: Adaptation of Subroutine SCAT for Use with P r o g r a m LIANA Adaptors' names and a d d r e s s e s : William R. Smith, Physics Department. Trinity University, San Antonio. Texas. USA Title of adaptation (32 c h a r a c t e r s maximum): ADAPT SCAT TO LIANA Adaptation number: 0002 R e f e r e n c e to original program: Catalog~te number: ACQF; Title: SCAT; R e f . in CPC : 1 (1969) 106 Authors of original program: William R. Smith No. of cards required to effect adaptation (including directive cards): 13
PROGRAM
SUMMARY
Title o f manuscript: H a u s e r - F e s h b a c h Nuclear Scattering Subroutine LIANA Authors' names and addresses: William R. Smith, Physics Department, Trinity University, San Antonio, Texas, USA Title o f program (32 c h a r a c t e r s maximum): LIANA Catalogue numbe~. ABKA Computer f o r which the program is designed and others upon which it is operable Computer: CDC 6400, CDC 6600. Installation: Computer Knowledge Corporation, San Antonio, Texas; r e s p . University of Texas, Austin, Texas, USA Operation s y s t e m or monitor under which the program is executed: SCOPE 3.1.6A P r o g r a m m i n g language used: FORTRAN High speed store required: 23,000 words. No. of bits p e r word : 60 Is the program overlaid? No No. of magnetic tapes required: None What other peripherals are used? None No. of cards in combined program and test deck: 571 CPC L i b r a r y subprograms used Catalogue number: ACQF; Title: SCAT; Ref. in CPC: 1 (1969) 106 Adaptation number: 0002; Adaptation title: ADAPT SCAT TO LIANA; Ref. in CPC: i (1969) 181. K e y w o r d s d e s c r i p t i v e o f problem and method o f solution: Nuclear, Reaction, Compound Nucleus, Compound-Elastic, H a u s e r - F e s h b a c h , Differential C r o s s Section, Integrated C r o s s Section, Optical Model, Woods-Saxon, SpinOrbit, Penetrabilities, Width Fluctuation Correction.
C-21 182
W. R. S M I T H
Nature o f physical problem S u b r o u t i n e LIANA c o m p u t e s the s c a t t e r i n g a n g u l a r d i s t r i b u t i o n and i n t e g r a t e d c r o s s s e c t i o n for any twop a r t i c l e n u c l e a r r e a c t i o n which p r o c e e d s t h r o u g h a s t a t i s t i c a l d i s t r i b u t i o n of c o m p o u n d n u c l e u s s t a t e s . T h e s u b r o u t i n e m a y be i n c o r p o r a t e d into o t h e r p r o g r a m s o r e l s e r u n with the d r i v e r p r o g r a m p r o v i d e d . T h e l a t t e r c a l l s s u b r o u t i n e SCAT I1] to obtain the n e c e s s a r y optical m o d e l n u c l e a r p e n e t r a b i l i t i e s . T h e input r e q u i r e d c o n s i s t s of the e n e r g i e s , c h a r g e s , m a s s e s , s p i n s , p a r itics, and optical m o d e l p a r a m e t e r s for all t w o - p a r t icle c h a n n e l s t h r o u g h which the c o m p o u n d n u c l e u s c a n decay.
Method o f solulion T h e H a u s e r - F e s h b a c h t h e o r y [21 g e n e r a l i z e d [31 to include s p i n - o r b i t coupling is e m p l o y e d . E i t h e r the plain H a u s c r - F e s h b a c h m e t h o d i m p l y i n g i n d e p e n d e n c e of c o m p o u n d n u c l e u s f o r m a t i o n and d e c a y or the m o d i fied H a u s c r - F e s h b a c h p r o c e d u r e i n c l u d i n g a c o r r e c tion 141 for an a s s u m e d P o r t e r - T h o m a s d i s t r i b u t i o n of c o m p o u n d n u c l e u s r e s o n a n c e widths [51 m a y be u t i l i z e d .
Resl~qctions on the complexity o f the problem Only p a r t i c l e c h a n n e l s can be t r e a t e d ; g a m m a - r a y c h a n n e l s a r e e x c l u d e d . N o r can c h a n n e l s c o n t a i n i n g t h r e e or m o r e s e p a r a t e p a r t i c l e s be handled. In any c h a n n e l the s p i n - o r b i t p a r t i c l e can h a v e s p i n up to 5
and the o t h e r p a r t i c l e m a y h a v e a n y s p i n . T h e m a x i m u m n u m b e r of p a r t i a l w a v e s n e c e s s a r y in any one c h a n n e l c a n be no g r e a t e r t h a n 31. A total of t e n t y p e s of two p a r t i c l e p a r t i t i o n s , and a total of 600 c h a n n e l s . m a y be s i m u l t a n e o u s l y t r e a t e d . T h e r e s t r i c t i o n s on m a x i m u m s p i n . n u m b e r of p a r t i a l w a v e s , n u m b e r of p a r t i c l e t y p e s , and n u m b e r of c h a n n e l s a r e due s o l e l y to s t o r a g e a l l o t m e n t s and m a y be e a s i l y r e l a x e d .
Typical running time T h e r u n n i n g t i m e is r o u g h l y p r o p o r t i o n a l to t h e n u m b e r of p e n e t r a b i l i t i e s c a l c u l a t e d , s i n c e the r e s t of the c a l c u l a t i o n is s h o r t by c o m p a r i s o n . T h e e a l e u l a t i o n a l t i m e s p e n t by LIANA a l o n e i n c r e a s e s a s t h e v a l u e s of the s p i n . n u m b e r of c h a n n e l s , and n u m b e r of p a r t i a l w a v e s i n c r e a s e . T h e t e s t c a s e r e q u i r e s 13 s e c o n d s to c o m p i l e and 3 s e c o n d s to r u n on a CDC 6400 c o m p u t e r .
References [1] W. R. Smith, C o m p u t e r P h y s . C o m m u n . 1 (1969) 106. [2] W. H a u s e r and H. F e s h b a c h , P h y s . Rev. 87 (1952) 366. [3] G. R. S a t e h l c r . P r o c . P h y s . Soe. (London) A68 (1955) 1041. [4] P. M o l d a u e r , P h y s . Rev. 123 (1961) 968. [5] C. E. P o r t e r and R. G. T h o m a s . P h y s . Rev. 104 {1956) 483.