atomscf.grid.radial_grid_log

atomscf.grid.radial_grid_log(n, rmin, rmax)[源代码]

生成对数(几何)径向网格及其梯形积分权重。

对数网格的定义(对数等差,适用于 Numerov):

\[r_i = r_\min\,\exp\!\left( i\,\Delta x \right),\ \ \Delta x = \frac{\ln(r_\max) - \ln(r_\min)}{N-1}.\]
参数:
  • n (int) -- 网格点数 \(N\),要求 \(N\ge 2\)

  • rmin (float) -- 径向下限 \(r_\min>0\),需严格大于 0 才能取对数。

  • rmax (float) -- 径向上限 \(r_\max\),应满足 \(r_\max > r_\min\)

返回类型:

tuple[ndarray, ndarray]

返回:

  • r (numpy.ndarray) -- 对数径向网格坐标 \(r_i\)

  • w (numpy.ndarray) -- 对应的梯形积分权重 \(w_i\)

备注

  • 对数网格能在小 \(r\) 处加密采样,适合处理库仑势与核附近行为。

  • 该网格满足 ln(r) 等差,适用于 Numerov 方法(参考 codex reply_check_4.md)

  • 若后续采用有限差分离散二阶导数,需使用非均匀网格的差分公式(本包已支持)。