Efforts have been made to solve the Dirac equation with axially deformed scalar and vector WoodsSaxon potentials in the coordinate space with the imaginary time step method. The results of the singleparticle energies ...Efforts have been made to solve the Dirac equation with axially deformed scalar and vector WoodsSaxon potentials in the coordinate space with the imaginary time step method. The results of the singleparticle energies thus obtained are consistent with those calculated with the basis expansion method, which demonstrates the feasibility of the imaginary time step method for the relativistic static problems.展开更多
Taking the single neutron levels of 12C in the Fermi sea as examples,the optimization of the imaginary time step(ITS) evolution with the box size and mesh size for the Dirac equation is investigated.For the weakly bou...Taking the single neutron levels of 12C in the Fermi sea as examples,the optimization of the imaginary time step(ITS) evolution with the box size and mesh size for the Dirac equation is investigated.For the weakly bound states,in order to reproduce the exact single-particle energies and wave functions,a relatively large box size is required.As long as the exact results can be reproduced,the ITS evolution with a smaller box size converges faster,while for both the weakly and deeply bound states,the ITS evolutions are less sensitive to the mesh size.Moreover,one can find a parabola relationship between the mesh size and the corresponding critical time step,i.e.,the largest time step to guarantee the convergence,which suggests that the ITS evolution with a larger mesh size allows larger critical time step,and thus can converge faster to the exact result.These conclusions are very helpful for optimizing the evolution procedure in the future self-consistent calculations.展开更多
Single particles moving in a reflection-asymmetric potential are investigated by solving the Schr6dinger equation of the reflectionasymmetric Nilsson Hamiltonian with the imaginary time method in 3D lattice space and ...Single particles moving in a reflection-asymmetric potential are investigated by solving the Schr6dinger equation of the reflectionasymmetric Nilsson Hamiltonian with the imaginary time method in 3D lattice space and the harmonic oscillator basis expansion method. In the 3D lattice calculation, the l2 divergence problem is avoided by introducing a damping function, and the（l2）N term in the non-spherical case is calculated by introducing an equivalent N-independent operator. The efficiency of these numerical techniques is demonstrated by solving the spherical Nilsson Hamiltonian in 3D lattice space. The evolution of the single-particle levels in a reflection-asvmmetric ootential is obtained and discussed bv the above two numerical methods, and their consistencv is shown in the obtained single-particle energies with the differences smaller than 10-4[hω0]展开更多
基金Supported by National Natural Science Foundation of China (10435010, 10775004, 10221003)Major State Basic Research Development Program (2007CB815000)
文摘Efforts have been made to solve the Dirac equation with axially deformed scalar and vector WoodsSaxon potentials in the coordinate space with the imaginary time step method. The results of the singleparticle energies thus obtained are consistent with those calculated with the basis expansion method, which demonstrates the feasibility of the imaginary time step method for the relativistic static problems.
基金supported partially by Guizhou Science and Technology Foundation (Grant No J[2010]2135)the National Basic Research Program of China (Grant No 2007CB815000)the National Natural Science Foundation of China (Grant Nos 10775004, 10947013, and 10975008)
文摘Taking the single neutron levels of 12C in the Fermi sea as examples,the optimization of the imaginary time step(ITS) evolution with the box size and mesh size for the Dirac equation is investigated.For the weakly bound states,in order to reproduce the exact single-particle energies and wave functions,a relatively large box size is required.As long as the exact results can be reproduced,the ITS evolution with a smaller box size converges faster,while for both the weakly and deeply bound states,the ITS evolutions are less sensitive to the mesh size.Moreover,one can find a parabola relationship between the mesh size and the corresponding critical time step,i.e.,the largest time step to guarantee the convergence,which suggests that the ITS evolution with a larger mesh size allows larger critical time step,and thus can converge faster to the exact result.These conclusions are very helpful for optimizing the evolution procedure in the future self-consistent calculations.
基金supported by the National Basic Research Program of China (Grant No. 2013CB834400)the National Natural Science Foundation of China (Grants Nos. 11335002, 11375015, 11461141002, and 11621131001)
文摘Single particles moving in a reflection-asymmetric potential are investigated by solving the Schr6dinger equation of the reflectionasymmetric Nilsson Hamiltonian with the imaginary time method in 3D lattice space and the harmonic oscillator basis expansion method. In the 3D lattice calculation, the l2 divergence problem is avoided by introducing a damping function, and the（l2）N term in the non-spherical case is calculated by introducing an equivalent N-independent operator. The efficiency of these numerical techniques is demonstrated by solving the spherical Nilsson Hamiltonian in 3D lattice space. The evolution of the single-particle levels in a reflection-asvmmetric ootential is obtained and discussed bv the above two numerical methods, and their consistencv is shown in the obtained single-particle energies with the differences smaller than 10-4[hω0]
