# Jobname
jobname  sample_us

# Atom number & Electron configuration
#   (calc_type : nonrel, srel, rel)
#   (spin_type : restricted, polarized)
atom_number      58   58.0
electron_config  19   srel  restricted
    1s   2.0   1
    2s   2.0   1
    2p   6.0   1
    3s   2.0   1
    3p   6.0   1
    3d  10.0   1
    4s   2.0   1
    4p   6.0   1
    4d  10.0   1
    4f   1.0   1
    5s   2.0   1
    5p   6.0   1
    5d   1.0   1
    5f   0.0   0
    6s   2.0   1
    6p   0.0   1
    6d   0.0   0
    7s   0.0   0
    7p   0.0   0

# Exchange-correlation potential 
#   (xc_type : ldapz81, ldapw91, ggapbe, xlda, none)
xc_potential   ldapw91

# Exchange-correlation potential 
#   (xc_type : ldapz81, ldapw91, ggapbe, xlda, none)
xc_potential   ggapbe

# Pseudopotential (pp_type : nc, general, none)
#                 (local : orbital, special)
pseudo_potential   general
   local     orbital    d
   orbitals   3
    5s       us      1
             2.20    0.0     1.1
    5p       us      1
             2.20    0.0     1.1
    5d       nc      1
             2.50    0.0     1.1
 phi_np       6     
 phi_gmesh   500     5.5    40.0
 qps_np       5     
 qps_gmesh   500    18.0    40.0

# Pseudopotential (pp_type : nc, general, none)
#                 (local : orbital, special)
pseudo_potential   general
   local     special    polynomial
             rcut       2.50
             vcut       0.00
   orbitals   6
    5s       us      1
             2.20    0.0     1.1
    5p       us      1
             2.20    0.0     1.1
    5d       nc      1
             2.50    0.0     1.1
    4f       us      1
             2.30    0.0     1.1
    6s       us      1
             2.50    0.0     1.1
    6p       us      1
             2.50    0.0     1.1
 phi_np       6     
 phi_gmesh   500     5.5    40.0
 qps_np       5     
 qps_gmesh   500    18.0    40.0

# Pseudopotential (pp_type : nc, general, none)
#                 (local : orbital, special)
pseudo_potential   general
   local     special    bhs
             rcut       1.50
             vcut       0.00
             lambda     6.50
   orbitals   6
    5s       us      1
             2.20    0.0     1.1
    5p       us      1
             2.20    0.0     1.1
    5d       nc      1
             2.50    0.0     1.1
    4f       us      1
             2.30    0.0     1.1
    6s       us      1
             2.50    0.0     1.1
    6p       us      1
             2.50    0.0     1.1
 phi_np       6     
 phi_gmesh   500     5.5    40.0
 qps_np       5     
 qps_gmesh   500    18.0    40.0

# Deficit charge
deficit_charge  angular
lmax_qps   6
   0    1.00
   1    1.02
   2    1.04
   3    1.06
   4    1.08
   5    1.10
   6    1.12

# Deficit charge
deficit_charge  noopt 

# Core correction (pcc_type : polynomial, sphbes, none)
core_correction  pcc  polynomial
   rcut        1.85
   np           4
   gmesh       500   9.0  40.0

# Core correction (pcc_type : polynomial, sphbes, none)
core_correction  pcc  polynomial
   ratio       2.00    
   np           4
   gmesh       500   9.0  40.0    

# Solve PP spin
solve_pp_spin
   spin  polarized  2.0
   mix    0.30

# Solve PP spin
solve_pp_spin
   spin  polarized  manual
     5s   1.0   1.0      
     5p   3.0   3.0
     5d   1.0   0.0      
     4f   1.0   0.0
     6s   1.0   1.0      
     6p   0.0   0.0
   mix    0.30

# Solve PP spin
solve_pp_spin
   spin  restricted  manual
     5s   2.0      
     5p   6.0
     5d   1.0      
     4f   1.0
     6s   2.0      
     6p   0.0
   mix    0.3

# Mesh (mesh_type : manual, standard)
rmesh   standard

# Mesh (mesh_type : manual, standard)
rmesh   manual
   num       1000
   range     1.d-6   60.d0

# Logarithmic derivative (logderi_type : default, manual, none)
logderi   none

# Logarithmic derivative (logderi_type : default, manual, none)
logderi   manual
   nonlocal  separable
   rcut      4.0
   num       701
   range    -5.0  +2.0

# Fourier transform (fourier_type : default, manual, none)
fourier   manual
   num       300
   range     0.05  15.0

# Convergence criteria
conv
   loop_conv    250
   dee_conv     5.d-9

# Order of calculation
order_pc         5
order_diff       4
order_nec        2
order_coeff      5

# Epsilon
eps_check      1.d-10
eps_de         1.d-10

# Switches
sw_calc_hubbard  0
sw_write_rho     1
sw_write_us      1
sw_write_qps     1
sw_write_pcc     1
sw_write_pp      1
sw_write_gncpp   1
sw_write_gncpp2  1
sw_write_ciaopp  1
sw_write_sol     1
sw_write_logderi 1
sw_write_fourier 1
sw_with_core     1
sw_with_ae       1
sw_with_dipole   1
sw_debug         0

# End of input data
end

