DFMDEF18: A C-code for the double folding interaction potential of a spherical nucleus with deformed nucleus

I. I. Gontchar, M. V. Chushnyakova

    Research output: Contribution to journalArticle

    Abstract

    This is a new version of the DFMDEF code published earlier. The new version is designed to obtain the nucleus–nucleus potential between two nuclei (one of which may be deformed) by using the double folding model (DFM). In particular the code enables to find the Coulomb barrier. The new version allows the user to employ his (her) own charge, proton, and neutron density distributions. The main functionalities of the original code (the nucleus–nucleus potential as a function of the incident angle and the distance between the centers of mass of colliding nuclei; the Coulomb barrier characteristics) have not been modified. New version program summary Program title: DFMDEF18 Program Files doi: http://dx.doi.org/10.17632/y9pwd2p24z.1 Licensing provisions: CC0 1.0 Programming language: C Journal reference of previous version: Comp. Phys. Comm. 184 (2013) 172 Does the new version supersede the previous version? Yes Nature of problem: The code calculates in a semimicroscopic way the bare interaction potential between two colliding nuclei one of which can be deformed and axially symmetric. The potential is evaluated as a function of the center of mass distance and the angle between the axis of symmetry and the beam direction. The heights and the positions of the Coulomb barriers are found. Dependence of the barrier parameters upon the characteristics of the effective NN forces (like, e.g. the range of the exchange part of the nuclear term) as well as upon the parameters of the density distributions can be investigated. Method of solution: The nucleus–nucleus potential is calculated using the double folding model with the Coulomb and the effective M3Y NN interactions. For the direct parts of the Coulomb and the nuclear terms, the Fourier transform method is used. In order to calculate the exchange parts, the density matrix expansion method is applied. Reason for new version: Many users asked us how to implement their own density distributions in the code. Now this option has been added. Summary of revisions: 1. Additional features of DFMDEF18: Projectile and target densities as input files In the DFMDEF [1, 2], only the Woods–Saxon profile for the charge, proton, and neutron density distributions was used. In the new version, DFMDEF18, the user can use the same profile, but there is an option to provide six input files <inp_rhoP_Z.c>, <inp_rhoP_N.c>, <inp_rhoP_q.c>, <inp_rhoT_Z.c>, <inp_rhoT_N.c>, and <inp_rhoT_q.c>similar to what was done in [3]. In these files the proton (_Z), neutron (_N), and charge (_q) density distributions are defined for the projectile (P) and target (T) nuclei as functions of the distance from the nucleus center (r) and the zenith angle (the_deg). We will refer to this set of six files as to “rho-input files”. The technical explanation of the rho-input files might be found in file <Program_changes.txt>(subsection 2.11). The values of the densities required for producing the interaction potential are found by means of the following interpolating polynomial. In order to avoid the overlap in notations for the zenith angle, instead of the traditional θ, we use here Θ. ρr,Θ=ρ[Formula presented]+ρ[Formula presented]+ρ[Formula presented]+ρ[Formula presented]+ρ[Formula presented]−1.Since the fractional difference does not exceed 0.5% we believe the approximation (2) can beaccepted. The original version of the code [1] allowed performing the calculations of the strong (nuclear) term of the nucleus–nucleus interaction potential in two different ways. Namely, the phenomenological Woods–Saxon parametrization [4, 5] and semimicroscopic double folding calculations could be used. In order to use the Woods–Saxon parametrization it is necessary to know the deformation parameters of the target nucleus. If the code uses the rho-input files for the densities provided by the user, such parameters can be absent. Therefore the Woods–Saxon parametrization for the potential is removed from the present version of the code. The Woods–Saxon approximation of the calculated DFM nuclear term of the potential is also removed for the same reason. Thus in the present code only the options concerning the double folding calculations remain. 2. The program The code consists now of 8 files and one header file. It reads the data from 8 input files and prints the results into two output files. For specific details regarding the changes in each source file see the file <Program_changes.txt>. The main input file has been split into two files: <inp_dfpdef.c> and <inp_dens.c>. Their description as well as the description of the output file might be also found in the file <Program_changes.txt>. References [1] I. I. Gontchar, M. V. Chushnyakova, Comput. Phys. Commun. 184 (2013) 172.[2] I. I. Gontchar, D. J. Hinde, M. Dasgupta, C. R. Morton, J. O. Newton, Phys. Rev. C 73 (2006) 034610.[3] I. I. Gontchar, M. V. Chushnyakova, Comput. Phys. Commun. 206 (2016) 97.[4] I.I. Gontchar, M. Dasgupta, D. J. Hinde, R. D. Butt, A. Mukherjee, Phys. Rev. C 65 (2002) 034610.[5] I.I. Gontchar, M. Dasgupta, D. J. Hinde, J. O. Newton, Phys. At. Nucl. 69 (2006) 1428 Appendix TEST RUN OUTPUT Input file <inp_dfpdef.c> [Figure presented] Output file <out_dfmsph.c> [Figure presented]

    Original languageEnglish
    Pages (from-to)414-417
    Number of pages4
    JournalComputer Physics Communications
    Volume222
    DOIs
    Publication statusPublished - 1 Jan 2018

    Fingerprint

    files
    folding
    Neutrons
    Protons
    Projectiles
    nuclei
    interactions
    Charge density
    Computer programming languages
    Fourier transforms
    Polynomials
    density distribution
    zenith
    neutrons
    newton
    center of mass
    projectiles
    output
    C (programming language)
    licensing

    Keywords

    • Coulomb barrier
    • Density dependent NN forces
    • Double folding model
    • M3Y-interaction
    • Nucleus–nucleus collision

    ASJC Scopus subject areas

    • Hardware and Architecture
    • Physics and Astronomy(all)

    Cite this

    DFMDEF18 : A C-code for the double folding interaction potential of a spherical nucleus with deformed nucleus. / Gontchar, I. I.; Chushnyakova, M. V.

    In: Computer Physics Communications, Vol. 222, 01.01.2018, p. 414-417.

    Research output: Contribution to journalArticle

    @article{2f63cfa7f9fa434bbdc0e95374758efd,
    title = "DFMDEF18: A C-code for the double folding interaction potential of a spherical nucleus with deformed nucleus",
    abstract = "This is a new version of the DFMDEF code published earlier. The new version is designed to obtain the nucleus–nucleus potential between two nuclei (one of which may be deformed) by using the double folding model (DFM). In particular the code enables to find the Coulomb barrier. The new version allows the user to employ his (her) own charge, proton, and neutron density distributions. The main functionalities of the original code (the nucleus–nucleus potential as a function of the incident angle and the distance between the centers of mass of colliding nuclei; the Coulomb barrier characteristics) have not been modified. New version program summary Program title: DFMDEF18 Program Files doi: http://dx.doi.org/10.17632/y9pwd2p24z.1 Licensing provisions: CC0 1.0 Programming language: C Journal reference of previous version: Comp. Phys. Comm. 184 (2013) 172 Does the new version supersede the previous version? Yes Nature of problem: The code calculates in a semimicroscopic way the bare interaction potential between two colliding nuclei one of which can be deformed and axially symmetric. The potential is evaluated as a function of the center of mass distance and the angle between the axis of symmetry and the beam direction. The heights and the positions of the Coulomb barriers are found. Dependence of the barrier parameters upon the characteristics of the effective NN forces (like, e.g. the range of the exchange part of the nuclear term) as well as upon the parameters of the density distributions can be investigated. Method of solution: The nucleus–nucleus potential is calculated using the double folding model with the Coulomb and the effective M3Y NN interactions. For the direct parts of the Coulomb and the nuclear terms, the Fourier transform method is used. In order to calculate the exchange parts, the density matrix expansion method is applied. Reason for new version: Many users asked us how to implement their own density distributions in the code. Now this option has been added. Summary of revisions: 1. Additional features of DFMDEF18: Projectile and target densities as input files In the DFMDEF [1, 2], only the Woods–Saxon profile for the charge, proton, and neutron density distributions was used. In the new version, DFMDEF18, the user can use the same profile, but there is an option to provide six input files , , , , , and similar to what was done in [3]. In these files the proton (_Z), neutron (_N), and charge (_q) density distributions are defined for the projectile (P) and target (T) nuclei as functions of the distance from the nucleus center (r) and the zenith angle (the_deg). We will refer to this set of six files as to “rho-input files”. The technical explanation of the rho-input files might be found in file <Program_changes.txt>(subsection 2.11). The values of the densities required for producing the interaction potential are found by means of the following interpolating polynomial. In order to avoid the overlap in notations for the zenith angle, instead of the traditional θ, we use here Θ. ρr,Θ=ρ[Formula presented]+ρ[Formula presented]+ρ[Formula presented]+ρ[Formula presented]+ρ[Formula presented]−1.Since the fractional difference does not exceed 0.5{\%} we believe the approximation (2) can beaccepted. The original version of the code [1] allowed performing the calculations of the strong (nuclear) term of the nucleus–nucleus interaction potential in two different ways. Namely, the phenomenological Woods–Saxon parametrization [4, 5] and semimicroscopic double folding calculations could be used. In order to use the Woods–Saxon parametrization it is necessary to know the deformation parameters of the target nucleus. If the code uses the rho-input files for the densities provided by the user, such parameters can be absent. Therefore the Woods–Saxon parametrization for the potential is removed from the present version of the code. The Woods–Saxon approximation of the calculated DFM nuclear term of the potential is also removed for the same reason. Thus in the present code only the options concerning the double folding calculations remain. 2. The program The code consists now of 8 files and one header file. It reads the data from 8 input files and prints the results into two output files. For specific details regarding the changes in each source file see the file <Program_changes.txt>. The main input file has been split into two files: and . Their description as well as the description of the output file might be also found in the file <Program_changes.txt>. References [1] I. I. Gontchar, M. V. Chushnyakova, Comput. Phys. Commun. 184 (2013) 172.[2] I. I. Gontchar, D. J. Hinde, M. Dasgupta, C. R. Morton, J. O. Newton, Phys. Rev. C 73 (2006) 034610.[3] I. I. Gontchar, M. V. Chushnyakova, Comput. Phys. Commun. 206 (2016) 97.[4] I.I. Gontchar, M. Dasgupta, D. J. Hinde, R. D. Butt, A. Mukherjee, Phys. Rev. C 65 (2002) 034610.[5] I.I. Gontchar, M. Dasgupta, D. J. Hinde, J. O. Newton, Phys. At. Nucl. 69 (2006) 1428 Appendix TEST RUN OUTPUT Input file [Figure presented] Output file [Figure presented]",
    keywords = "Coulomb barrier, Density dependent NN forces, Double folding model, M3Y-interaction, Nucleus–nucleus collision",
    author = "Gontchar, {I. I.} and Chushnyakova, {M. V.}",
    year = "2018",
    month = "1",
    day = "1",
    doi = "10.1016/j.cpc.2017.09.008",
    language = "English",
    volume = "222",
    pages = "414--417",
    journal = "Computer Physics Communications",
    issn = "0010-4655",
    publisher = "Elsevier",

    }

    TY - JOUR

    T1 - DFMDEF18

    T2 - A C-code for the double folding interaction potential of a spherical nucleus with deformed nucleus

    AU - Gontchar, I. I.

    AU - Chushnyakova, M. V.

    PY - 2018/1/1

    Y1 - 2018/1/1

    N2 - This is a new version of the DFMDEF code published earlier. The new version is designed to obtain the nucleus–nucleus potential between two nuclei (one of which may be deformed) by using the double folding model (DFM). In particular the code enables to find the Coulomb barrier. The new version allows the user to employ his (her) own charge, proton, and neutron density distributions. The main functionalities of the original code (the nucleus–nucleus potential as a function of the incident angle and the distance between the centers of mass of colliding nuclei; the Coulomb barrier characteristics) have not been modified. New version program summary Program title: DFMDEF18 Program Files doi: http://dx.doi.org/10.17632/y9pwd2p24z.1 Licensing provisions: CC0 1.0 Programming language: C Journal reference of previous version: Comp. Phys. Comm. 184 (2013) 172 Does the new version supersede the previous version? Yes Nature of problem: The code calculates in a semimicroscopic way the bare interaction potential between two colliding nuclei one of which can be deformed and axially symmetric. The potential is evaluated as a function of the center of mass distance and the angle between the axis of symmetry and the beam direction. The heights and the positions of the Coulomb barriers are found. Dependence of the barrier parameters upon the characteristics of the effective NN forces (like, e.g. the range of the exchange part of the nuclear term) as well as upon the parameters of the density distributions can be investigated. Method of solution: The nucleus–nucleus potential is calculated using the double folding model with the Coulomb and the effective M3Y NN interactions. For the direct parts of the Coulomb and the nuclear terms, the Fourier transform method is used. In order to calculate the exchange parts, the density matrix expansion method is applied. Reason for new version: Many users asked us how to implement their own density distributions in the code. Now this option has been added. Summary of revisions: 1. Additional features of DFMDEF18: Projectile and target densities as input files In the DFMDEF [1, 2], only the Woods–Saxon profile for the charge, proton, and neutron density distributions was used. In the new version, DFMDEF18, the user can use the same profile, but there is an option to provide six input files , , , , , and similar to what was done in [3]. In these files the proton (_Z), neutron (_N), and charge (_q) density distributions are defined for the projectile (P) and target (T) nuclei as functions of the distance from the nucleus center (r) and the zenith angle (the_deg). We will refer to this set of six files as to “rho-input files”. The technical explanation of the rho-input files might be found in file <Program_changes.txt>(subsection 2.11). The values of the densities required for producing the interaction potential are found by means of the following interpolating polynomial. In order to avoid the overlap in notations for the zenith angle, instead of the traditional θ, we use here Θ. ρr,Θ=ρ[Formula presented]+ρ[Formula presented]+ρ[Formula presented]+ρ[Formula presented]+ρ[Formula presented]−1.Since the fractional difference does not exceed 0.5% we believe the approximation (2) can beaccepted. The original version of the code [1] allowed performing the calculations of the strong (nuclear) term of the nucleus–nucleus interaction potential in two different ways. Namely, the phenomenological Woods–Saxon parametrization [4, 5] and semimicroscopic double folding calculations could be used. In order to use the Woods–Saxon parametrization it is necessary to know the deformation parameters of the target nucleus. If the code uses the rho-input files for the densities provided by the user, such parameters can be absent. Therefore the Woods–Saxon parametrization for the potential is removed from the present version of the code. The Woods–Saxon approximation of the calculated DFM nuclear term of the potential is also removed for the same reason. Thus in the present code only the options concerning the double folding calculations remain. 2. The program The code consists now of 8 files and one header file. It reads the data from 8 input files and prints the results into two output files. For specific details regarding the changes in each source file see the file <Program_changes.txt>. The main input file has been split into two files: and . Their description as well as the description of the output file might be also found in the file <Program_changes.txt>. References [1] I. I. Gontchar, M. V. Chushnyakova, Comput. Phys. Commun. 184 (2013) 172.[2] I. I. Gontchar, D. J. Hinde, M. Dasgupta, C. R. Morton, J. O. Newton, Phys. Rev. C 73 (2006) 034610.[3] I. I. Gontchar, M. V. Chushnyakova, Comput. Phys. Commun. 206 (2016) 97.[4] I.I. Gontchar, M. Dasgupta, D. J. Hinde, R. D. Butt, A. Mukherjee, Phys. Rev. C 65 (2002) 034610.[5] I.I. Gontchar, M. Dasgupta, D. J. Hinde, J. O. Newton, Phys. At. Nucl. 69 (2006) 1428 Appendix TEST RUN OUTPUT Input file [Figure presented] Output file [Figure presented]

    AB - This is a new version of the DFMDEF code published earlier. The new version is designed to obtain the nucleus–nucleus potential between two nuclei (one of which may be deformed) by using the double folding model (DFM). In particular the code enables to find the Coulomb barrier. The new version allows the user to employ his (her) own charge, proton, and neutron density distributions. The main functionalities of the original code (the nucleus–nucleus potential as a function of the incident angle and the distance between the centers of mass of colliding nuclei; the Coulomb barrier characteristics) have not been modified. New version program summary Program title: DFMDEF18 Program Files doi: http://dx.doi.org/10.17632/y9pwd2p24z.1 Licensing provisions: CC0 1.0 Programming language: C Journal reference of previous version: Comp. Phys. Comm. 184 (2013) 172 Does the new version supersede the previous version? Yes Nature of problem: The code calculates in a semimicroscopic way the bare interaction potential between two colliding nuclei one of which can be deformed and axially symmetric. The potential is evaluated as a function of the center of mass distance and the angle between the axis of symmetry and the beam direction. The heights and the positions of the Coulomb barriers are found. Dependence of the barrier parameters upon the characteristics of the effective NN forces (like, e.g. the range of the exchange part of the nuclear term) as well as upon the parameters of the density distributions can be investigated. Method of solution: The nucleus–nucleus potential is calculated using the double folding model with the Coulomb and the effective M3Y NN interactions. For the direct parts of the Coulomb and the nuclear terms, the Fourier transform method is used. In order to calculate the exchange parts, the density matrix expansion method is applied. Reason for new version: Many users asked us how to implement their own density distributions in the code. Now this option has been added. Summary of revisions: 1. Additional features of DFMDEF18: Projectile and target densities as input files In the DFMDEF [1, 2], only the Woods–Saxon profile for the charge, proton, and neutron density distributions was used. In the new version, DFMDEF18, the user can use the same profile, but there is an option to provide six input files , , , , , and similar to what was done in [3]. In these files the proton (_Z), neutron (_N), and charge (_q) density distributions are defined for the projectile (P) and target (T) nuclei as functions of the distance from the nucleus center (r) and the zenith angle (the_deg). We will refer to this set of six files as to “rho-input files”. The technical explanation of the rho-input files might be found in file <Program_changes.txt>(subsection 2.11). The values of the densities required for producing the interaction potential are found by means of the following interpolating polynomial. In order to avoid the overlap in notations for the zenith angle, instead of the traditional θ, we use here Θ. ρr,Θ=ρ[Formula presented]+ρ[Formula presented]+ρ[Formula presented]+ρ[Formula presented]+ρ[Formula presented]−1.Since the fractional difference does not exceed 0.5% we believe the approximation (2) can beaccepted. The original version of the code [1] allowed performing the calculations of the strong (nuclear) term of the nucleus–nucleus interaction potential in two different ways. Namely, the phenomenological Woods–Saxon parametrization [4, 5] and semimicroscopic double folding calculations could be used. In order to use the Woods–Saxon parametrization it is necessary to know the deformation parameters of the target nucleus. If the code uses the rho-input files for the densities provided by the user, such parameters can be absent. Therefore the Woods–Saxon parametrization for the potential is removed from the present version of the code. The Woods–Saxon approximation of the calculated DFM nuclear term of the potential is also removed for the same reason. Thus in the present code only the options concerning the double folding calculations remain. 2. The program The code consists now of 8 files and one header file. It reads the data from 8 input files and prints the results into two output files. For specific details regarding the changes in each source file see the file <Program_changes.txt>. The main input file has been split into two files: and . Their description as well as the description of the output file might be also found in the file <Program_changes.txt>. References [1] I. I. Gontchar, M. V. Chushnyakova, Comput. Phys. Commun. 184 (2013) 172.[2] I. I. Gontchar, D. J. Hinde, M. Dasgupta, C. R. Morton, J. O. Newton, Phys. Rev. C 73 (2006) 034610.[3] I. I. Gontchar, M. V. Chushnyakova, Comput. Phys. Commun. 206 (2016) 97.[4] I.I. Gontchar, M. Dasgupta, D. J. Hinde, R. D. Butt, A. Mukherjee, Phys. Rev. C 65 (2002) 034610.[5] I.I. Gontchar, M. Dasgupta, D. J. Hinde, J. O. Newton, Phys. At. Nucl. 69 (2006) 1428 Appendix TEST RUN OUTPUT Input file [Figure presented] Output file [Figure presented]

    KW - Coulomb barrier

    KW - Density dependent NN forces

    KW - Double folding model

    KW - M3Y-interaction

    KW - Nucleus–nucleus collision

    UR - http://www.scopus.com/inward/record.url?scp=85031320971&partnerID=8YFLogxK

    UR - http://www.scopus.com/inward/citedby.url?scp=85031320971&partnerID=8YFLogxK

    U2 - 10.1016/j.cpc.2017.09.008

    DO - 10.1016/j.cpc.2017.09.008

    M3 - Article

    AN - SCOPUS:85031320971

    VL - 222

    SP - 414

    EP - 417

    JO - Computer Physics Communications

    JF - Computer Physics Communications

    SN - 0010-4655

    ER -