HFSZEEMAN95—A program for computing weak and intermediate magneticfield and hyperfineinduced transition rates
(2020) In Computer Physics Communications Abstract
HFSZEEMAN95 is an updated and extended Fortran 95 version of the HFSZEEMAN program (Andersson and Jönsson, 2008). Given relativistic atomic state functions generated by the GRASP2018 package (Fischer et al., 2019), HFSZEEMAN95 together with the accompanying Matlab/GNU Octave program MITHIT allows for: (1) the computation and plotting of Zeeman energy splittings of magnetic fine and hyperfine structure substates as functions of the strength of an external magnetic field, (2) the computation of transition rates between different magnetic fine and hyperfine structure substates in the presence of an external magnetic field and rates of hyperfineinduced transitions in the field free limit, (3) the synthesization of spectral profiles for... (More)
HFSZEEMAN95 is an updated and extended Fortran 95 version of the HFSZEEMAN program (Andersson and Jönsson, 2008). Given relativistic atomic state functions generated by the GRASP2018 package (Fischer et al., 2019), HFSZEEMAN95 together with the accompanying Matlab/GNU Octave program MITHIT allows for: (1) the computation and plotting of Zeeman energy splittings of magnetic fine and hyperfine structure substates as functions of the strength of an external magnetic field, (2) the computation of transition rates between different magnetic fine and hyperfine structure substates in the presence of an external magnetic field and rates of hyperfineinduced transitions in the field free limit, (3) the synthesization of spectral profiles for transitions obtained from (2). With the new features, HFSZEEMAN95 and the accompanying Matlab/GNU Octave program MITHIT are useful for the analysis of observational spectra and to resolve the complex features due to the splitting of the fine and hyperfine levels. Program summary: Program Title: HFSZEEMAN95 Program Files doi: http://dx.doi.org/10.17632/rv2vycs7pg.1 Licensing provisions: GNU General Public License 3 Programming language: Fortran 95, Matlab/GNU Octave Nature of problem: Calculation of transition energies and rates between different magnetic fine and hyperfine structure substates in the presence of an external magnetic field and rates of hyperfine induced transitions in the field free limit. Synthesization of spectral profiles. Solution method: Wave functions for magnetic fine structure substates in the field free case are given by atomic state functions (ASFs). The ASFs are expansions over configuration state functions (CSFs) ΓJM_{J}〉=∑γc_{γ}γJM_{J}〉.The ASFs are computed by the GRASP2018 relativistic atomic structure package (Fischer et al., 2018) and are supposed to be available. Wave functions for magnetic fine structure substates in an external magnetic field are expanded in a basis of ASFs Γ˜M_{J}〉=∑ΓJd_{ΓJ}ΓJM_{J}〉.Wave functions for magnetic hyperfine structure substates in an external magnetic field are expanded in a basis of the combined nuclear and atomic system Γ˜IM_{F}〉=∑ΓJFd_{ΓJF}ΓIJFM_{F}〉,where ΓIJFM_{F}〉 are coupled nuclear and atomic functions ΓIJFM_{F}〉=∑M_{I},M_{J}〈IJM_{I}M_{J}IJFM_{F}〉IM_{I}〉ΓJM_{J}〉.Reduced hyperfine and Zeeman matrix elements, used to construct the total interaction matrix in the given basis, are computed as sums over reduced oneparticle matrix elements of orbitals building the CSFs. By diagonalizing the interaction matrix, Zeeman energy splittings of fine and hyperfine structure substates are obtained together with the expansion coefficients of the basis functions. Transition rates between different magnetic fine and hyperfine structure substates are computed as sums over reduced transition matrix elements between fine structure states weighted by the expansion coefficients of the basis functions and angular factors. Given energies of the magnetic substates along with transition rates, a synthetic spectrum is obtained by convolving the spectral lines with a Gaussian function with a user defined value of the full width half maximum (FWHM). Additional comments including restrictions and unusual features : 1. The complexity of the cases that can be handled is determined by the GRASP2018 package used for the generation of the atomic state functions. 2. The current programs can be used for the calculations of electric dipole (E1), electric quadrupole (E2), magnetic dipole (M1) and magnetic quadrupole (M2) magneticfield and hyperfineinduced transitions, which are caused by the mixing of states due to hyperfine and Zeeman interaction. 3. The present model does not include the necessary nonperturbative treatment of the uncommon case involving neardegeneracies where the radiative width of a fine structure state is of the same order as the hyperfine or magneticfield perturbation; an effect often termed radiation damping (Indelicato et al., 1989 [1]; Robicheaux et al., 1995 [2]; Johnson et al., 1997 [3]).
(Less)
 author
 Li, Wenxian ; Grumer, Jon ^{LU} ; Brage, Tomas ^{LU} and Jönsson, Per
 organization
 publishing date
 2020
 type
 Contribution to journal
 publication status
 published
 subject
 keywords
 Hyperfine structure, Hyperfineinduced transitions, Magnetic field, Magneticfieldinduced transitions, Multiconfiguration Dirac–Hartree–Fock＋Breit, Relativistic atomic wave functions, Unexpected transitions, Zeeman effect
 in
 Computer Physics Communications
 article number
 107211
 publisher
 Elsevier
 external identifiers

 scopus:85079836427
 ISSN
 00104655
 DOI
 10.1016/j.cpc.2020.107211
 language
 English
 LU publication?
 yes
 id
 39b7d42dd17444f3846940e074891730
 date added to LUP
 20200318 13:37:46
 date last changed
 20201229 03:35:28
@article{39b7d42dd17444f3846940e074891730, abstract = {<p>HFSZEEMAN95 is an updated and extended Fortran 95 version of the HFSZEEMAN program (Andersson and Jönsson, 2008). Given relativistic atomic state functions generated by the GRASP2018 package (Fischer et al., 2019), HFSZEEMAN95 together with the accompanying Matlab/GNU Octave program MITHIT allows for: (1) the computation and plotting of Zeeman energy splittings of magnetic fine and hyperfine structure substates as functions of the strength of an external magnetic field, (2) the computation of transition rates between different magnetic fine and hyperfine structure substates in the presence of an external magnetic field and rates of hyperfineinduced transitions in the field free limit, (3) the synthesization of spectral profiles for transitions obtained from (2). With the new features, HFSZEEMAN95 and the accompanying Matlab/GNU Octave program MITHIT are useful for the analysis of observational spectra and to resolve the complex features due to the splitting of the fine and hyperfine levels. Program summary: Program Title: HFSZEEMAN95 Program Files doi: http://dx.doi.org/10.17632/rv2vycs7pg.1 Licensing provisions: GNU General Public License 3 Programming language: Fortran 95, Matlab/GNU Octave Nature of problem: Calculation of transition energies and rates between different magnetic fine and hyperfine structure substates in the presence of an external magnetic field and rates of hyperfine induced transitions in the field free limit. Synthesization of spectral profiles. Solution method: Wave functions for magnetic fine structure substates in the field free case are given by atomic state functions (ASFs). The ASFs are expansions over configuration state functions (CSFs) ΓJM<sub>J</sub>〉=∑γc<sub>γ</sub>γJM<sub>J</sub>〉.The ASFs are computed by the GRASP2018 relativistic atomic structure package (Fischer et al., 2018) and are supposed to be available. Wave functions for magnetic fine structure substates in an external magnetic field are expanded in a basis of ASFs Γ˜M<sub>J</sub>〉=∑ΓJd<sub>ΓJ</sub>ΓJM<sub>J</sub>〉.Wave functions for magnetic hyperfine structure substates in an external magnetic field are expanded in a basis of the combined nuclear and atomic system Γ˜IM<sub>F</sub>〉=∑ΓJFd<sub>ΓJF</sub>ΓIJFM<sub>F</sub>〉,where ΓIJFM<sub>F</sub>〉 are coupled nuclear and atomic functions ΓIJFM<sub>F</sub>〉=∑M<sub>I</sub>,M<sub>J</sub>〈IJM<sub>I</sub>M<sub>J</sub>IJFM<sub>F</sub>〉IM<sub>I</sub>〉ΓJM<sub>J</sub>〉.Reduced hyperfine and Zeeman matrix elements, used to construct the total interaction matrix in the given basis, are computed as sums over reduced oneparticle matrix elements of orbitals building the CSFs. By diagonalizing the interaction matrix, Zeeman energy splittings of fine and hyperfine structure substates are obtained together with the expansion coefficients of the basis functions. Transition rates between different magnetic fine and hyperfine structure substates are computed as sums over reduced transition matrix elements between fine structure states weighted by the expansion coefficients of the basis functions and angular factors. Given energies of the magnetic substates along with transition rates, a synthetic spectrum is obtained by convolving the spectral lines with a Gaussian function with a user defined value of the full width half maximum (FWHM). Additional comments including restrictions and unusual features : 1. The complexity of the cases that can be handled is determined by the GRASP2018 package used for the generation of the atomic state functions. 2. The current programs can be used for the calculations of electric dipole (E1), electric quadrupole (E2), magnetic dipole (M1) and magnetic quadrupole (M2) magneticfield and hyperfineinduced transitions, which are caused by the mixing of states due to hyperfine and Zeeman interaction. 3. The present model does not include the necessary nonperturbative treatment of the uncommon case involving neardegeneracies where the radiative width of a fine structure state is of the same order as the hyperfine or magneticfield perturbation; an effect often termed radiation damping (Indelicato et al., 1989 [1]; Robicheaux et al., 1995 [2]; Johnson et al., 1997 [3]).</p>}, author = {Li, Wenxian and Grumer, Jon and Brage, Tomas and Jönsson, Per}, issn = {00104655}, language = {eng}, publisher = {Elsevier}, series = {Computer Physics Communications}, title = {HFSZEEMAN95—A program for computing weak and intermediate magneticfield and hyperfineinduced transition rates}, url = {http://dx.doi.org/10.1016/j.cpc.2020.107211}, doi = {10.1016/j.cpc.2020.107211}, year = {2020}, }