（1）短程库仑势的计算

• \rho^\text{s}(\boldsymbol{r}^\text{s})\rightarrow\rho^\text{s}(\boldsymbol{G}^\text{s})
• \phi^\text{s}_\text{S}(\boldsymbol{r}^\text{s})=V_\text{s}(\boldsymbol{G}^\text{s})\rho^\text{s}(\boldsymbol{G}^\text{s})
• \phi_\text{S}^\text{s}(\boldsymbol{G}^\text{s})\rightarrow\phi_\text{S}^\text{s}(\boldsymbol{r}^\text{s})

这部分计算对应的调用为：

call FFT_NEW( DSBX%FFT_STATE, rhoIn, rhoG_DSBX_short )
potG_DSBX_short = rhoG_DSBX_short * htkernel_short
call FFT_NEW( DSBX%FFT_STATE, potG_DSBX_short, potR_DSBX_short )

这里的htkernel_short对应V_\text{s}(\boldsymbol{G}^\text{s})。

（2）长程库仑势的计算（一）

• \rho^\text{s}_\text{m}(\boldsymbol{G}^\text{s})=m(\boldsymbol{G}^\text{s})\rho^\text{s}(\boldsymbol{G}^\text{s})
• \rho^\text{s}_\text{m}(\boldsymbol{G}^\text{s})\rightarrow\rho'^\text{s}_\text{m}(\boldsymbol{G}^\text{s})
• \rho'^\text{s}_\text{m}(\boldsymbol{G}^\text{s})\rightarrow\rho'^\text{s}_\text{m}(\boldsymbol{r}'^\text{s})
• \rho'^\text{s}_\text{m}(\boldsymbol{r}'^\text{s})\rightarrow\rho'_\text{m}(\boldsymbol{r}')

这部分对应的调用为：

rhoG_DSBX_short = rhoG_DSBX_short * maskG_DSBX
call FFT_NEW( DSBX%FFT_STATE, rhoG_DSBX_short, rhoR_DSBX_short )
call DSBX_2_SGBX( rhoR_DSBX_short, rhoR_SGBX )

这里的maskG_DSBX对应m(\boldsymbol{G}^\text{s})。注意到这里的计算流程是\rho^\text{s}_\text{m}(\boldsymbol{G}^\text{s})\rightarrow\rho^\text{s}_\text{m}(\boldsymbol{r}^\text{s})\rightarrow\rho'_\text{m}(\boldsymbol{r}')，与第一篇中不同，但是效果是相同的。

（3）长程库仑势的计算（二）

• \rho'_\text{m}(\boldsymbol{r}')\rightarrow\rho'_\text{m}(\boldsymbol{G})
• \phi'_\text{L}(\boldsymbol{G})=\rho'_\text{m}(\boldsymbol{G})\cdot[V_\text{L}(\boldsymbol{G})/m(\boldsymbol{G})]
• \phi'_\text{L}(\boldsymbol{G})\rightarrow\phi'_\text{L}(\boldsymbol{r}')

这部分对应的调用为：

if (SGBX%active_proc == .TRUE. ) then
    CALL FFT_NEW( SGBX%FFT_STATE, rhoR_SGBX, rhoG_SGBX )
    rhoG_SGBX = rhoG_SGBX * htkernel_long / maskG_SGBX
    CALL FFT_NEW( SGBX%FFT_STATE, rhoG_SGBX, potR_SGBX )
end if

这里的htkernel_long对应的是V_\text{L}(\boldsymbol{G})。注意到由于这部分计算是对SGBX网格进行全局FFT，因此只有第一个网格（ctive_proc标记为.TRUE.）进行相关计算。

（4）长程库仑势的计算（三）

• \phi'_\text{L}(\boldsymbol{r}')\rightarrow\phi'^\text{s}_\text{L}(\boldsymbol{r}'^\text{s})
• W'^\text{s}(\boldsymbol{r}'^\text{s})=M'(\boldsymbol{r}'^\text{s})\phi'^\text{s}_\text{L}(\boldsymbol{r}'^\text{s})
• W'^\text{s}(\boldsymbol{r}'^\text{s})\rightarrow W'^\text{s}(\boldsymbol{G}^\text{s})

这部分对应的调用为：

call SGBX_2_SSBX( potR_SGBX, potR_SSBX )
call GetBufferData( SSBX, potR_SSBX )
potR_SSBX = potR_SSBX * maskX_SSBX
call FFT_NEW( SSBX%FFT_STATE, potR_SSBX, potG_SSBX )

这里的maskX_SSBX对应M'(\boldsymbol{r}'^\text{s})。注意到这里还调用了一个GetBufferData子程序（位于DataExchange_FillBuffer.f90文件中DataExchange_FillBuffer），其作用是将每个小块与附近26个小块交换核心部分数据以填充缓冲区。

（5）长程库仑势的计算（四）

• W'^\text{s}(\boldsymbol{G}^\text{s})\rightarrow W^\text{s}(\boldsymbol{G}^\text{s})
• W^\text{s}(\boldsymbol{G}^\text{s})\rightarrow W^\text{s}(\boldsymbol{r}^\text{s})
• \phi^\text{s}_\text{L}(\boldsymbol{r}^\text{s} )=W^\text{s} (\boldsymbol{r}^\text{s})/M(\boldsymbol{r}^\text{s} )

这部分对应的调用为：

call SSBX_2_DSBX( potG_SSBX, potG_DSBX_long )
call FFT_NEW( DSBX%FFT_STATE, potG_DSBX_long, potR_DSBX_long )
potR_DSBX_long = potR_DSBX_long / maskX_DSBX

这里的maskX_DSBX对应M(\boldsymbol{r}^\text{s})。

（6）将短程与长程库仑势相加

• \phi^\text{s}(\boldsymbol{r}^\text{s})=\phi_\text{S}^\text{s}(\boldsymbol{r}^\text{s})+\phi_\text{L}^\text{s}(\boldsymbol{r}^\text{s})

这部分对应的代码为：

potential = potR_DSBX_short + potR_DSBX_long