目录
radial_grid_linear()
radial_grid_log()
radial_grid_exp_transformed()
trapezoid_weights()
radial_grid_mixed()
radial_hamiltonian_matrix()
solve_bound_states_fd()
radial_hamiltonian_matrix_linear_fd5()
solve_bound_states_fd5()
solve_bound_states_fd5_auxlinear()
solve_bound_states_transformed()
build_transformed_hamiltonian()
SCFConfig
SCFResult
run_lsda_x_only()
run_lsda_pz81()
run_lsda_vwn()
HFConfig
HFResult
run_hf_minimal()
HFSCFConfig
HFSCFResult
run_hf_scf_s()
HFSCFGeneralConfig
HFSCFGeneralResult
run_hf_scf()
v_hartree()
OrbitalSpec
default_occupations()
trapz()
normalize_radial_u()
slater_integral_radial()
slater_integral_k0()
SlaterIntegralCache
allowed_k_values()
coupling_factor_ak()
wigner_3j_squared()
exchange_operator_s()
exchange_operator_general()
exchange_operator_general_spin()
vx_dirac()
ex_dirac_density()
lda_c_pz81()
lda_c_vwn()
运行 LSDA X-only 的自洽场计算(球对称、径向)。
cfg (SCFConfig) -- 自洽计算的配置。
verbose (bool)
progress_every (int)
自洽结果(是否收敛、迭代数、能级/波函数、密度与势)。
备注
占据方案以径向通道为单位(指定 \(\ell\) 与通道内序号 \(n_{\text{index}}\))。
对 C 原子,2p 壳层采用 m 平均且自旋极化:上自旋每个 m 的分数占据 \(2/3\),下自旋为 0。
混合策略:密度线性混合(数值稳健,易用)。