搜索

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

(1+1)维非线性薛定谔方程PT对称势函数的数值反演

王扬 徐映红 赵晔丹 张立溥

引用本文:
Citation:

(1+1)维非线性薛定谔方程PT对称势函数的数值反演

王扬, 徐映红, 赵晔丹, 张立溥

Numerical Inversion of PT-Symmetric Potential Functions in (1+1)-Dimensional Nonlinear Schrödinger Equations

WANG Yang, XU Yinghong, ZHAO Yedan, ZHANG Lipu
Article Text (iFLYTEK Translation)
PDF
HTML
导出引用
  • 本文针对一类具有(1+1)维非线性项的薛定谔方程对称势函数的反问题进行了深入研究. 在研究过程中, 我们设定部分边界处的波函数为已知条件, 而将势函数的数值作为未知量进行求解. 通过采用三层差分格式, 我们将原始的非线性方程成功转化为实数域上的非线性方程组. 在此基础上, 我们运用非精确牛顿法对转化后的非线性方程组进行了有效的求解. 研究结果充分证实了该方法在解决此类反问题方面的高效性, 为相关领域的数值求解提供了一种创新的思路和有力的工具.
    This research investigates the inverse problem of reconstructing the PT-symmetric potential in a class of $ (1+1) $-dimensional nonlinear Schrödinger equations. The governing equation is given by:$ \text{i}u_t(x,t) + u_{xx}(x,t) + \alpha\left| u(x,t) \right|^2 u(x,t) + \beta\left| u(x,t) \right|^4 u(x,t) + V_{PT}(x) u(x,t) = 0, $where $ u(x, t) $ denotes the wave function in dimensionless coordinates, and the PT-symmetric potential $ V_{PT}(x) = V(x) + {\bf{i}}W(x) $ consists of a real part $ V(x) $ and an imaginary part $ {\bf{i}}W(x) $, satisfying the symmetry conditions $ V(x) = V(-x) $ and $ W(x) = -W(x) $.In this inverse problem, partial boundary values of the wave function are known, while the potential $ V_{PT}(x) $ is the unknown to be reconstructed. To address this challenge, we construct a three-level finite difference scheme for the corresponding forward problem, discretizing both the wave function and the potential. This approach leads to a nonlinear system of equations that links the known wave data to the unknown potential values. To simplify the computation, we separate the real and imaginary parts of this system and reformulate it as a real-valued nonlinear system of equations.For the numerical solution, we employ an inexact Newton method to iteratively solve the resulting nonlinear system. In each iteration, the Jacobian matrix is approximated numerically. To ensure that the reconstructed potential strictly satisfies the PT-symmetry, a parity correction mechanism is introduced at the end of the iteration process.We conduct numerical experiments under both noise-free (exact data) and noisy (inexact data) conditions. The results demonstrate that in both cases, the proposed method converges within a limited number of iterations and maintains the reconstruction error within the order of $ 10^{-3} $. These findings validate the effectiveness and robustness of the proposed method in solving inverse problems involving PT-symmetric potentials, offering an innovative and practical approach for related numerical applications.
  • 图 1  $ V_{pt}(x) $实部及虚部反演结果 图(a) $ V_{pt}(x) $实部的实际精确值与INTA和mPINNs反演结果的对比; 图(b) $ V_{pt}(x) $虚部的实际精确值与INTA和mPINNs反演结果的对比; 图(c) INTA反演结果实部与虚部的绝对误差

    Fig. 1.  Real and imaginary part inversion results of $ V_{pt}(x) $ (a) Comparison of the real part of $ V_{pt}(x) $ with the inversion results of INTA and mPINNs; (b) Comparison of the imaginary part of $ V_{pt}(x) $ with the inversion results of INTA and mPINNs; (c) Absolute error of the real and imaginary parts of INTA inversion results.

    图 2  误差 (a) INTA每次迭代后的误差error; (b) mPINNs与精确解$ u(x, z) $之间的误差

    Fig. 2.  Errors (a) error after each iteration of INTA; (b) Error between mPINNs and the exact solution $ u(x, z) $.

    图 3  $ V_{pt}(x) $实部及虚部反演结果($ U^* $包含噪声) (a) INTA和mPINNs反演结果与$ V_{pt}(x) $实部的实际精确值的对比; (b) INTA和mPINNs反演结果与$ V_{pt}(x) $虚部的实际精确值的对比; (c) INTA反演结果实部与虚部的绝对误差

    Fig. 3.  Real and imaginary part inversion results of $ V_{pt}(x) $ ($ U^* $ with noise) (a) Comparison of the real part of $ V_{pt}(x) $ with the inversion results of INTA and mPINNs; (b) Comparison of the imaginary part of $ V_{pt}(x) $ with the inversion results of INTA and mPINNs; (c) Absolute error of the real and imaginary parts of INTA inversion results.

    图 4  误差($ U^* $包含噪声) (a) INTA每次迭代后的误差error; (b) mPINNs与精确解$ u(x, z) $之间的误差

    Fig. 4.  Errors($ U^* $ with noise) (a) error after each iteration of INTA; (b) Error between mPINNs and the exact solution $ u(x, z) $.

    表 1  INTA与mPINNs对比

    Table 1.  Comparison between INTA and mPINNs.

    算法 实部误差 虚部误差 运行时间(秒)
    INTA $ 6.1768\times10^{-13} $ $ 7.6891\times10^{-12} $ 61.3
    mPINNs 0.0416 0.01549 113.8
    下载: 导出CSV

    表 2  绝对误差

    Table 2.  Absolute Error.

    x $ 0.1 $ $ 0.2 $ $ 0.3 $ $ 0.4 $ $ 0.5 $ $ 0.6 $ $ 0.7 $ $ 0.8 $ $ 0.9 $
    INTA实部 $ 6.4\text{E-13} $ $ 6.5\text{E-13} $ $ 6.1\text{E-13} $ $ 6.0\text{E-13} $ $ 6.1\text{E-13} $ $ 6.4\text{E-13} $ $ 7.3\text{E-13} $ $ 7.8\text{E-13} $ $ 8.0\text{E-13} $
    mPINNs实部 $ 5.3\text{E-03} $ $ 8.8\text{E-03} $ $ 8.6\text{E-03} $ $ 4.8\text{E-03} $ $ 1.5\text{E-03} $ $ 8.5\text{E-03} $ $ 1.4\text{E-02} $ $ 1.6\text{E-02} $ $ 1.1\text{E-02} $
    INTA虚部 $ 1.4\text{E-12} $ $ 1.6\text{E-12} $ $ 1.3\text{E-12} $ $ 7.5\text{E-13} $ $ 5.2\text{E-13} $ $ 4.7\text{E-13} $ $ 6.8\text{E-13} $ $ 9.5\text{E-13} $ $ 5.5\text{E-13} $
    mPINNs虚部 $ 1.4\text{E-02} $ $ 1.7\text{E-03} $ $ 4.0\text{E-03} $ $ 2.3\text{E-03} $ $ 4.0\text{E-03} $ $ 5.3\text{E-03} $ $ 6.1\text{E-03} $ $ 5.7\text{E-03} $ $ 4.3\text{E-03} $
    下载: 导出CSV

    表 3  INTA与mPINNs对比($ U^* $包含噪声)

    Table 3.  Comparison between INTA and mPINNs ($ U^* $ with noise).

    算法 实部误差 虚部误差 运行时间(秒)
    INTA 0.0043 0.0126 71.5
    mPINNs 0.0098 0.0452 112.3
    下载: 导出CSV

    表 4  绝对误差($ U^* $包含噪声)

    Table 4.  Absolute Error($ U^* $ with noise).

    x $ 0.1 $ $ 0.2 $ $ 0.3 $ $ 0.4 $ $ 0.5 $ $ 0.6 $ $ 0.7 $ $ 0.8 $ $ 0.9 $
    INTA实部 $ 1.9\text{E-04} $ $ 2.3\text{E-03} $ $ 4.2\text{E-03} $ $ 4.4\text{E-03} $ $ 2.8\text{E-03} $ $ 2.8\text{E-03} $ $ 3.1\text{E-03} $ $ 5.2\text{E-03} $ $ 3.2\text{E-03} $
    mPINNs实部 $ 8.9\text{E-03} $ $ 1.1\text{E-02} $ $ 1.2\text{E-02} $ $ 1.1\text{E-02} $ $ 8.5\text{E-03} $ $ 5.9\text{E-03} $ $ 3.5\text{E-03} $ $ 1.9\text{E-03} $ $ 1.3\text{E-03} $
    INTA虚部 $ 3.5\text{E-03} $ $ 2.3\text{E-03} $ $ 2.6\text{E-03} $ $ 8.7\text{E-04} $ $ 2.3\text{E-04} $ $ 1.8\text{E-03} $ $ 2.4\text{E-03} $ $ 3.8\text{E-03} $ $ 1.8\text{E-03} $
    mPINNs虚部 $ 2.4\text{E-03} $ $ 6.2\text{E-03} $ $ 1.0\text{E-02} $ $ 1.6\text{E-02} $ $ 1.6\text{E-02} $ $ 1.6\text{E-02} $ $ 1.5\text{E-02} $ $ 1.2\text{E-02} $ $ 7.1\text{E-03} $
    下载: 导出CSV
  • [1]

    Liu Y P, Gao Y T, Wei G M 2012 Phys. A Stat. Mech. Appl. 391 535Google Scholar

    [2]

    Uthayakumar A, Han Y G, Lee S B 2006 Chaos Solitons Fractals 29 916Google Scholar

    [3]

    Bender C M, Boettcher S 1998 Phys. Rev. Lett. 80 5243Google Scholar

    [4]

    El-Ganainy R, Makris K G, Christodoulides D N, Musslimani Z H 2007 Opt. Lett. 32 2632Google Scholar

    [5]

    Bender C M 2007 Rep. Prog. Phys. 70 947Google Scholar

    [6]

    Baudouin L, Puel J P 2002 Inverse Probl. 18 1537Google Scholar

    [7]

    Avdonin S A, Mikhaylov A S, Mikhaylov V S, Park J C 2021 Inverse Problems 37 035002Google Scholar

    [8]

    Raissi M, Perdikaris P, Karniadakis G E 2019 J. Comput. Phys. 378 686Google Scholar

    [9]

    Zhou Z, Yan Z 2021 Phys. Lett. A 387 127010Google Scholar

    [10]

    Li J, Li B 2021 Commun. Theor. Phys. 73 125001Google Scholar

    [11]

    张坤 2024 理论数学 14 117

    Zhang K 2024 Pure Math. 14 117

    [12]

    Qiu W X, Geng K L, Zhu B W, Liu W, Li J T, Dai C Q 2024 Nonlinear Dyn. 112 10215Google Scholar

    [13]

    Song J, Yan Z 2023 Physica D 448 133729Google Scholar

    [14]

    Wang S, Wang H, Perdikaris P 2021 Comput. Methods Appl. Mech. Eng. 384 113938Google Scholar

    [15]

    Liu Y, Wu R, Jiang Y 2024 J. Comput. Phys. 518 113341Google Scholar

    [16]

    Karniadakis G E, Kevrekidis I G, Lu L, Perdikaris P, Wang S, Yang L 2021 Nat. Rev. Phys. 3 422Google Scholar

    [17]

    Lu L, Meng X, Mao Z, Karniadakis G E 2021 SIAM Rev. 63 208Google Scholar

    [18]

    孙志忠 2022 偏微分方程数值解法 第三版 (北京: 科学出版社) 第 320—349页

    Sun Z Z 2022 Numerical Solution of Partial Differential Equations. 3 rd edn. (Beijing: Science Press), pp 320–349

    [19]

    Dembo R S, Eisenstat S C, Steihaug T 1982 SIAM J. Numer. Anal. 19 400Google Scholar

    [20]

    Göksel İ, Antar N, Bakırtaş İ 2015 Opt. Commun. 354 277Google Scholar

  • [1] 陈韬, 江普庆. 揭示热反射实验中热物性参数的本征关系. 物理学报, doi: 10.7498/aps.73.20241369
    [2] 谢国大, 潘攀, 任信钢, 冯乃星, 方明, 李迎松, 黄志祥. 高阶SF-SFDTD方法在含时薛定谔方程求解中的应用研究. 物理学报, doi: 10.7498/aps.73.20230771
    [3] 王鹏举, 范俊宇, 苏艳, 赵纪军. 基于机器学习构建的环三亚甲基三硝胺晶体势. 物理学报, doi: 10.7498/aps.69.20200690
    [4] 张海峰, 王文旭. 复杂系统重构. 物理学报, doi: 10.7498/aps.69.20200001
    [5] 孙帅令, 冷永刚, 张雨阳, 苏徐昆, 范胜波. 双磁铁多稳态悬臂梁磁力及势能函数分析. 物理学报, doi: 10.7498/aps.69.20191981
    [6] 葛伟宽, 薛纭, 楼智美. 完整力学系统的广义梯度表示. 物理学报, doi: 10.7498/aps.63.110202
    [7] 惠治鑫, 贺鹏飞, 戴瑛, 吴艾辉. 硅功能化石墨烯热导率的分子动力学模拟. 物理学报, doi: 10.7498/aps.63.074401
    [8] 黄启灿, 胡淑娟, 邱春雨, 李宽, 于海鹏, 丑纪范. 基于无导数优化方法的数值模式误差估计. 物理学报, doi: 10.7498/aps.63.149203
    [9] 周耐根, 胡秋发, 许文祥, 李克, 周浪. 硅熔化特性的分子动力学模拟–-不同势函数的对比研究. 物理学报, doi: 10.7498/aps.62.146401
    [10] 周耐根, 洪涛, 周浪. MEAM势与Tersoff势比较研究碳化硅熔化与凝固行为. 物理学报, doi: 10.7498/aps.61.028101
    [11] 沈晶, 沙威, 黄志祥, 陈明生, 吴先良. 含时Schrdinger方程的高阶辛FDTD算法研究. 物理学报, doi: 10.7498/aps.61.190202
    [12] 李会山, 李鹏程, 周效信. 强激光场中模型氢原子的势函数对产生高次谐波强度的影响. 物理学报, doi: 10.7498/aps.58.7633
    [13] 刘 奎, 丁宏林, 张贤高, 余林蔚, 黄信凡, 陈坤基. 量子点浮置栅量子线沟道三栅结构单电子场效应管存储特性的数值模拟. 物理学报, doi: 10.7498/aps.57.7052
    [14] 刘 冬, 王 飞, 黄群星, 严建华, 池 涌, 岑可法. 二维弥散介质温度场的快速重建. 物理学报, doi: 10.7498/aps.57.4812
    [15] 厉江帆, 单树民, 杨建坤, 姜宗福. 失谐量子频率转换系统薛定谔方程的显式解析解. 物理学报, doi: 10.7498/aps.56.5597
    [16] 辛国锋, 陈国鹰, 花吉珍, 赵润, 康志龙, 冯荣珠, 安振峰. 941nm大功率应变单量子阱激光器的波长设计. 物理学报, doi: 10.7498/aps.53.1293
    [17] 张解放, 徐昌智, 何宝钢. 变量分离法与变系数非线性薛定谔方程的求解探索. 物理学报, doi: 10.7498/aps.53.3652
    [18] 李培咸, 郝 跃, 范 隆, 张进城, 张金凤, 张晓菊. 基于量子微扰的AlGaN/GaN异质结波函数半解析求解. 物理学报, doi: 10.7498/aps.52.2985
    [19] 尤云祥, 缪国平. 阻抗障碍物声散射的反问题. 物理学报, doi: 10.7498/aps.51.270
    [20] 刘剑波, 蔡喜平. 一维定态薛定谔方程的宏观模拟解法. 物理学报, doi: 10.7498/aps.50.820
计量
  • 文章访问数:  253
  • PDF下载量:  13
  • 被引次数: 0
出版历程
  • 上网日期:  2025-05-10

/

返回文章
返回