2025-11-29T08:01:19.757167

Generalized one-dimensional nonpolynomial SchrÃ¶dinger equation for Bose-Einstein condensates with generic transverse confinement

Basso, Cardoso

This work presents a dimensional reduction of Bose-Einstein condensates confined by generalized transverse potentials, parametrized by an exponent $n$. Starting from the three-dimensional Gross-Pitaevskii equation, we employ a variational ansatz to derive an effective one-dimensional nonpolynomial SchrÃ¶dinger equation, which self-consistently determines the transverse width dynamics. The model generalizes existing formalisms for cigar- and funnel-shaped geometries. We validate the approach through comprehensive numerical tests, demonstrating excellent agreement with full 3D simulations for ground-state properties across various interaction regimes. Finally, real-time simulations of matter-wave scattering at potential barriers verify the model's dynamical robustness, successfully replicating the spatiotemporal evolution and energy-dependent transmission characteristics observed in full 3D calculations.

academic

Generalized one-dimensional nonpolynomial Schrödinger equation for Bose-Einstein condensates with generic transverse confinement

基本信息

论文ID: 2511.14973
标题: Generalized one-dimensional nonpolynomial Schrödinger equation for Bose-Einstein condensates with generic transverse confinement
作者: Andréia M. Basso, Wesley B. Cardoso (Instituto de Física, Universidade Federal de Goiás, Brazil)
分类: cond-mat.quant-gas, nlin.PS
提交时间: 2025年11月18日
论文链接: https://arxiv.org/abs/2511.14973

摘要

本研究针对受广义横向势约束的玻色-爱因斯坦凝聚体(BEC)进行维度约化。从三维Gross-Pitaevskii方程出发，采用变分ansatz方法推导出有效的一维非多项式薛定谔方程(1D-NPSE)，该方程能自洽地确定横向宽度动力学。该模型统一推广了雪茄型和漏斗型几何的现有形式体系。通过全面的数值测试验证了该方法，在各种相互作用区间的基态性质上与完整3D模拟展现出优异的一致性。最后，势垒处物质波散射的实时模拟验证了模型的动力学鲁棒性，成功复现了完整3D计算中观察到的时空演化和能量依赖的透射特性。

研究背景与动机

研究问题

本研究旨在解决强各向异性阱中玻色-爱因斯坦凝聚体的维度约化问题，特别是针对具有广义横向约束势的系统。

问题重要性

计算效率：完整的三维Gross-Pitaevskii方程(3D-GPE)数值求解计算成本高昂，特别是对于强各向异性系统
物理理解：有效的低维模型能够揭示系统的本质物理机制，便于理论分析
实验相关性：实际实验中的BEC通常在强各向异性阱中实现，需要准确的低维描述

现有方法局限

特定几何限制：现有模型主要针对特定几何构型（如谐振子势的雪茄型40或漏斗型65）
缺乏统一框架：不同约束势需要单独推导，缺少系统性的处理方法
适用范围有限：现有模型难以描述连续变化的约束几何

研究动机

作者引入参数化的广义横向势 $V_\perp(r) = \frac{m\omega_\perp^2 n^2}{8}(\alpha^{4-2n}r^{2n} - 2\beta^{4-n}r^n/r^2)$ ，通过调节指数参数 $n \in \mathbb{R}_+^*$ 来统一描述不同的约束几何，建立一个连续可调的理论框架。

核心贡献

提出广义横向约束势：引入参数 $n$ 表征的广义非谐振横向势，统一了雪茄型( $n=2$ )和漏斗型( $n=1$ )几何
推导统一的1D有效方程：通过变分方法导出适用于连续 $n$ 值的一维非多项式薛定谔方程(1D-EFG)
获得横向宽度解析表达式：导出变分场 $\sigma^2$ 的解析解，自洽描述横向宽度动力学
全面数值验证：在基态性质、临界相互作用强度、动力学散射等多个方面验证模型准确性
建立统一理论框架：为研究不同约束几何下的BEC提供了系统性的理论工具

方法详解

任务定义

输入：三维Gross-Pitaevskii方程及广义横向约束势
输出：有效的一维非多项式薛定谔方程及横向宽度的自洽解
约束：强各向异性条件( $\gamma = \omega_z/\omega_\perp \ll 1$ )，绝热近似(横向和纵向时间尺度分离)

模型架构

1. 广义横向势构造

无量纲化后的横向势： $V_\perp(r) = \frac{n^2}{8}\left(r^{2n} - \frac{2r^n}{r^2}\right)$

关键性质：

$n=1$ ：漏斗型，原点处奇异
$n=2$ ：谐振子势，雪茄型
$n>1$ ：势阱最小值位于 $r_c = \left(\frac{n-2}{n-1}\right)^{1/n}$

2. 变分ansatz设计

采用高斯型径向分布： $\psi(r,t) = \sqrt{\frac{n}{2\pi\Gamma(2/n)}} \exp\left(-\frac{r^n}{2\sigma^2}\right) \frac{f(z,t)}{\sigma^{2/n}}$

其中：

$f(z,t)$ ：归一化轴向波函数
$\sigma(z,t)$ ：横向宽度变分场
归一化因子保证 $2\pi\int_0^\infty |\psi|^2 r dr = |f|^2$

3. 有效拉格朗日密度

通过绝热近似（横向宽度视为常数）得到： $L_{\text{eff}} = \frac{i}{2}\left(f\frac{\partial f^*}{\partial t} - f^*\frac{\partial f}{\partial t}\right) + \frac{1}{2}\left|\frac{\partial f}{\partial z}\right|^2 + [V(z)+1]|f|^2$ $+ \frac{n^2(\sigma^2-1)^2}{8\sigma^{4/n}\Gamma(2/n)}|f|^2 + \frac{gn}{\sigma^{4/n}2^{\frac{n+2}{n}}\Gamma(2/n)}|f|^4$

4. Euler-Lagrange方程

对 $f^*$ 和 $\sigma$ 变分得到：

轴向方程(1D-EFG)： $i\frac{\partial f}{\partial t} = -\frac{1}{2}\frac{\partial^2 f}{\partial z^2} + V(z)f + \frac{n}{\sigma^{4/n}\Gamma(2/n)}\left[\frac{n}{8}(\sigma^2-1)^2 + \frac{g}{2^{2/n}}|f|^2\right]f$