用户名: 密码: 验证码:
A computer program for two-particle intrinsic coefficients of fractional parentage
详细信息查看全文 | 推荐本文 |
摘要
A Fortran 90 program CESOS for the calculation of the two-particle intrinsic coefficients of fractional parentage for several j-shells with isospin and an arbitrary number of oscillator quanta (CESOs) is presented. The implemented procedure for CESOs calculation consistently follows the principles of antisymmetry and translational invariance. The approach is based on a simple enumeration scheme for antisymmetric many-particle states, efficient algorithms for calculation of the coefficients of fractional parentage for j-shells with isospin, and construction of the subspace of the center-of-mass Hamiltonian eigenvectors corresponding to the minimal eigenvalue equal to 3/2 (in 鈩徬?/em>). The program provides fast calculation of CESOs for a given particle number and produces results possessing small numerical uncertainties. The introduced CESOs may be used for calculation of expectation values of two-particle nuclear shell-model operators within the isospin formalism.

Program summary

Program title: CESOS

Catalogue identifier: AELT_v1_0

Program summary URL:

Program obtainable from: CPC Program Library, Queen始s University, Belfast, N. Ireland

Licensing provisions: Standard CPC licence,

No. of lines in distributed program, including test data, etc.: 10鈥?32

No. of bytes in distributed program, including test data, etc.: 61鈥?23

Distribution format: tar.gz

Programming language: Fortran 90

Computer: Any computer with a Fortran 90 compiler

Operating system: Windows XP, Linux

RAM: The memory demand depends on the number of particles A and the excitation energy of the system . Computation of the particle system with the total angular momentum and the total isospin requires around 4 kB of RAM at at , and at .

Classification: 17.18

Nature of problem: The code CESOS generates a list of two-particle intrinsic coefficients of fractional parentage for several j-shells with isospin.

Solution method: The method is based on the observation that CESOs may be obtained by diagonalizing the center-of-mass Hamiltonian in the basis set of antisymmetric A-particle oscillator functions with singled out dependence on Jacobi coordinates of two last particles and choosing the subspace of its eigenvectors corresponding to the minimal eigenvalue equal to 3/2.

Restrictions: One run of the code CESOS generates CESOs for one specified set of () values only. The restrictions on the () values are completely determined by the restrictions on the computation of the single-shell CFPs and two-particle multishell CFPs (GCFPs) [1]. The full sets of single-shell CFPs may be calculated up to the shell (for any particular shell of the configuration); the shell with cannot get full (it is the implementation constraint). The calculation of GCFPs is limited by when (due to the memory constraints); small numbers of particles allow significantly higher excitations. Any allowed values of J and T may be chosen for the specified values of A and . The complete list of allowed values of J and T for the chosen values of A and may be generated by the GCFP program - CPC Program Library, Catalogue Id. AEBI_v1_0. The actual scale of the CESOs computation problem depends strongly on the magnitude of the A and values. Though there are no limitations on A and values (within the limits of single-shell CFPs and multishell CFPs calculation), however the generation of corresponding list of CESOs is the subject of available computing resources. For example, the computing time of CESOs for , at took around 14 hours. The system with , at requires around 15 hours. These computations were performed on Pentium 3 GHz PC with 1 GB RAM [2].

Unusual features: It is possible to test the computed CESOs without saving them to a file. This allows the user to learn their number and approximate computation time and to evaluate the accuracy of calculations.

Additional comments: The program CESOS uses the code from GCFP program for calculation of the two-particle multishell coefficients of fractional parentage.

Running time: It depends on the size of the problem. The particle system with the took around 31 seconds on Pentium 3 GHz PC with 1 GB RAM at and about 2.6 hours at .

References:

[1]

A. Deveikis, A. Juodagalvis, Comp. Phys. Comm. 179 (2008) 607.

[2]

A. Deveikis, A. Fedurtsya, Phys. Part. Nucl. Lett. 8 (2011) 484.

© 2004-2018 中国地质图书馆版权所有 京ICP备05064691号 京公网安备11010802017129号

地址:北京市海淀区学院路29号 邮编:100083

电话:办公室:(+86 10)66554848;文献借阅、咨询服务、科技查新:66554700