solve_ivp丢弃复解虚部：狄拉克方程正负电子对产生求解异常

• 问题场景：在求解正负电子对产生的狄拉克方程自由基展开问题时，使用SciPy的solve_ivp积分器数值求解Mocken等人提出的方程组，运行时触发复数转实数丢弃虚部的警告，得到的正负电子对产生概率曲线与文献、Mathematica计算结果均不匹配。

• 解决方案：将初始值数组y0显式声明为复数类型（例如通过numpy.complex128转换），即可避免虚部被自动截断，得到符合预期的计算结果。

• 问题背景：该问题对应SciPy官方的Issue #8453，已存在至少4年，核心原因是solve_ivp默认会根据初始值的类型判断求解问题的域——若初始值为实数，即使方程组会产生复数解，积分器也会强制按实数问题处理，直接丢弃虚部，最终导致数值结果严重偏差。

内容的提问来源于stack exchange，提问作者derdotte