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

Basic Information

• Paper ID: 2511.14973
• Title: Generalized one-dimensional nonpolynomial Schrödinger equation for Bose-Einstein condensates with generic transverse confinement
• Authors: Andréia M. Basso, Wesley B. Cardoso (Instituto de Física, Universidade Federal de Goiás, Brazil)
• Classification: cond-mat.quant-gas, nlin.PS
• Submission Date: November 18, 2025
• Paper Link: https://arxiv.org/abs/2511.14973

Abstract

This study performs dimensional reduction for Bose-Einstein condensates (BEC) subject to generalized transverse potential confinement. Starting from the three-dimensional Gross-Pitaevskii equation, the authors employ a variational ansatz method to derive an effective one-dimensional nonpolynomial Schrödinger equation (1D-NPSE) that self-consistently determines transverse width dynamics. The model unifies and generalizes existing formulations for cigar-shaped and funnel-shaped geometries. The method is validated through comprehensive numerical tests, demonstrating excellent agreement with full 3D simulations across various interaction regimes for ground state properties. Finally, real-time simulations of matter-wave scattering at potential barriers validate the model's dynamical robustness, successfully reproducing spatiotemporal evolution and energy-dependent transmission characteristics observed in complete 3D calculations.

Research Background and Motivation

Research Problem

This study addresses the dimensional reduction problem for Bose-Einstein condensates in strongly anisotropic traps, particularly for systems with generic transverse confinement potentials.

Problem Significance

1. Computational Efficiency: Full three-dimensional Gross-Pitaevskii equation (3D-GPE) numerical solutions are computationally expensive, especially for strongly anisotropic systems
2. Physical Understanding: Effective low-dimensional models reveal the essential physical mechanisms of the system, facilitating theoretical analysis
3. Experimental Relevance: Practical BEC experiments typically realize condensates in strongly anisotropic traps, requiring accurate low-dimensional descriptions

Limitations of Existing Methods

1. Geometric Specificity: Existing models primarily target specific geometric configurations (e.g., cigar-shaped for harmonic oscillator potentials 40 or funnel-shaped 65)
2. Lack of Unified Framework: Different confinement potentials require separate derivations, lacking systematic treatment
3. Limited Applicability: Existing models struggle to describe continuously varying confinement geometries

Research Motivation

The authors introduce a parameterized generalized transverse potential $V_\perp(r) = \frac{m\omega_\perp^2 n^2}{8}(\alpha^{4-2n}r^{2n} - 2\beta^{4-n}r^n/r^2)$, using the exponent parameter $n \in \mathbb{R}_+^*$ to unify descriptions of different confinement geometries, establishing a continuously tunable theoretical framework.

Core Contributions

1. Proposed Generalized Transverse Confinement Potential: Introduced a generalized non-harmonic transverse potential characterized by parameter $n$, unifying cigar-shaped ($n=2$) and funnel-shaped ($n=1$) geometries
2. Derived Unified 1D Effective Equation: Through variational methods, derived a one-dimensional nonpolynomial Schrödinger equation (1D-EFG) applicable to continuous values of $n$
3. Obtained Analytical Expression for Transverse Width: Derived analytical solutions for the variational field $\sigma^2$, self-consistently describing transverse width dynamics
4. Comprehensive Numerical Validation: Validated model accuracy across multiple aspects including ground state properties, critical interaction strengths, and dynamical scattering
5. Established Unified Theoretical Framework: Provided systematic theoretical tools for studying BEC under different confinement geometries

Detailed Methodology

Task Definition

Input: Three-dimensional Gross-Pitaevskii equation and generalized transverse confinement potential
Output: Effective one-dimensional nonpolynomial Schrödinger equation and self-consistent solution for transverse width
Constraints: Strong anisotropy condition ($\gamma = \omega_z/\omega_\perp \ll 1$), adiabatic approximation (separation of transverse and longitudinal time scales)

Model Architecture

1. Generalized Transverse Potential Construction

Dimensionless transverse potential: $V_\perp(r) = \frac{n^2}{8}\left(r^{2n} - \frac{2r^n}{r^2}\right)$

Key properties:

• $n=1$: Funnel-shaped, singular at origin
• $n=2$: Harmonic oscillator potential, cigar-shaped
• $n>1$: Potential minimum located at $r_c = \left(\frac{n-2}{n-1}\right)^{1/n}$

2. Variational Ansatz Design

Gaussian-type radial distribution: $\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}}$

Where:

• $f(z,t)$: Normalized axial wave function
• $\sigma(z,t)$: Transverse width variational field
• Normalization factor ensures $2\pi\int_0^\infty |\psi|^2 r dr = |f|^2$

3. Effective Lagrangian Density

Through adiabatic approximation (treating transverse width as constant): $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 Equations

Variational equations for $f^*$ and $\sigma$ yield:

Axial Equation (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$

Transverse Width Solution: $\sigma_\pm^2 = \frac{n-2 \pm \sqrt{\frac{4^{2n-1}}{n}g|f|^2(n-1) + n^2}}{2(n-1)}$

Physical solution corresponds to positive branch $\sigma_+^2$ (negative branch leads to complex chemical potential).

Technical Innovations

1. Parametric Unification: Continuous interpolation of different confinement geometries through single parameter $n$, avoiding separate derivations for each potential
2. Self-Consistent Variational Field: Analytical solution for $\sigma^2$ depends on local density $|f|^2$, achieving self-consistent description of transverse-longitudinal coupling
3. Analytical Tractability: Despite being nonpolynomial, the equation retains analytical structure facilitating theoretical analysis
4. Correct Asymptotic Behavior: As $n\to\infty$, $\sigma^2\to 1$, automatically reducing to cubic nonlinearity model

Experimental Setup

Numerical Methods

• Algorithm: Second-order Split-Step Fourier method 69
• Imaginary Time Evolution: Ground state solution (gradient flow minimization)
• Real Time Evolution: Dynamical process simulation

System Parameters

• Axial Potential: $V(z) = \gamma^2 z^2/2$, anisotropy parameter $\gamma^2 = 0.1$
• Parameter $n$: Integer values $n = 1, 2, 3, 4, 5, 6$ examined
• Interaction Strength: $g = 1, 10, 100$ (repulsive) and $g = -0.5$ (attractive)

Comparison Methods

1. Full 3D-GPE (Equation 5): Reference benchmark
2. 1D Cubic Model (1D-CN) (Equation 17): Simplified model with fixed transverse width
3. 1D-EFG Model (Equation 10): Proposed generalized effective model

Evaluation Metrics

Static Properties:

• Axial density distribution $\rho(z) = |f(z)|^2$
• Chemical potential $\mu$ (Equation 12)
• Expectation value of longitudinal length squared $\langle z^2\rangle$
• Critical interaction strength $g_c$ (collapse threshold)

Dynamical Properties:

• Transmission coefficient $T$ (Equations 25-26)
• Kinetic/potential energy ratio $U/E$

Experimental Results

Main Results

1. Ground State Axial Density (Figure 3)

• Trend: Increasing $n$ and $g$ leads to significant axial stretching
• Mechanism: Stronger transverse confinement and repulsive interactions force atoms to distribute along the axis

2. Density Peak Comparison (Figure 4)

Conclusion: 1D-EFG maintains accuracy in strong interaction regime, while 1D-CN is only effective for weak interactions.

3. Chemical Potential and $\langle z^2\rangle$ (Figure 5)

• Nonlinear Effects: Most pronounced at low $n$
• Convergence Behavior: Approaches linear dependence for large $n$
• Accuracy Degradation: Discrepancy between 1D-EFG and 3D increases under strong repulsive interactions (enhanced transverse-longitudinal coupling violates weak coupling assumption)

4. Critical Interaction Strength (Table I)

Findings:

• 1D-EFG correctly predicts collapse phenomenon (1D-CN completely fails)
• Relative error increases with $n$
• Physical essence correct; quantitative accuracy depends on $n$

Dynamical Scattering Experiments

Experimental Setup

• Initial State: Ground state from imaginary time evolution in displaced harmonic potential
• Dynamics: Potential removed at $t=0$, momentum $p=3$ imparted
• Scattering Potential: Gaussian barrier $V_b(z) = A\exp(-z^2/B^2)$, $A=5$, $B=\sqrt{2}$
• Interaction: $g=-0.5$ (attractive)

Axial Density Evolution (Figure 6)

• Qualitative Consistency: Similar dynamics exhibited for all $n$ values
• Interference Effects: More pronounced interference with barrier at low $n$

Transmission Coefficient (Figure 7)

• Monotonicity: $T$ increases with $n$
• Accuracy: 1D-EFG highly consistent with 3D-GPE
  • $n=2$: Relative deviation 0.094% (1D-EFG) vs 0.153% (1D-CN)
  • $n=5$: 0.108% vs 0.168%
  • $n=10$: 0.119% vs 0.122%

Energy Ratio $U/E$ (Figure 8)

• Trend: Increases with $n$
• Accuracy: 1D-EFG consistent with 3D across all $n$ values, 1D-CN shows notable deviation at low $n$

Ablation Studies

Role of Transverse Width Variational Field (Figure 2)

• $\sigma_+^2$ Branch: Physical solution, increases with $g|f|^2$
• $\sigma_-^2$ Branch: Non-physical (negative values lead to complex chemical potential)
• Asymptotic Behavior: As $n\to\infty$, $\sigma^2\to 1$, convergence rate inversely proportional to $g|f|^2$

Cubic Model Comparison

Fixed $\sigma=1$ cubic model fails in:

1. Collapse prediction ($g_c$ nonexistent)
2. Density peak severely deviates under strong interactions
3. Dynamical transmission coefficient shows larger discrepancies

Experimental Findings

1. Confinement Geometry Effects: Enhanced transverse confinement (increasing $n$) results in:
   • Axial stretching
   • Increased transmission probability
   • Weakened nonlinear effects
2. Interaction Strength Dependence:
   • Weak Interactions: All models perform well
   • Strong Interactions: Only 1D-EFG maintains accuracy
3. Time Scale Separation Validity: Adiabatic approximation remains valid during dynamical processes
4. Dimensional Reduction Precision: Quantitative accuracy depends on $n$ and $g$, but physical essence always correct

Related Work

Classification of Dimensional Reduction Methods

1. Variational Methods

• Salasnich et al. (2002) 40: 1D-NPSE for cigar-shaped condensates
• dos Santos et al. (2019) 65: Funnel-shaped geometry
• This Work: Unified generalization to arbitrary $n$

2. Adiabatic Approximation Methods

• Jackson et al. (1998) 39: Transverse-longitudinal mode decoupling
• Mateo & Delgado (2008-2009) 50,51: Disk-shaped and tube-shaped condensates
• This work employs similar assumptions while retaining transverse compressibility

3. Special Geometry Extensions

• Tubular 64, Periodic Transverse Potential 46, Spatially Modulated Nonlinearity 60
• Anisotropic Harmonic Traps 52, Vortex Dynamics 41,51
• Binary Mixtures 57, Spin-Orbit Coupling 59,61

Advantages of This Work

1. Parametric Continuity: Single parameter $n$ interpolates different geometries, facilitating systematic study
2. Analytical Tractability: Retains analytical solutions for transverse width
3. Comprehensive Validation: Full testing of static and dynamical properties
4. Practical Utility: Applicable to non-ideal confinement potentials in experiments

Conclusions and Discussion

Main Conclusions

1. Theoretical Contribution: Successfully derived parameterized generalized 1D effective equation, unifying transverse confinement for $n\in\mathbb{R}_+^*$
2. Numerical Validation: Ground state properties, critical interaction strengths, and dynamical scattering highly consistent with 3D-GPE
3. Model Superiority: 1D-EFG significantly outperforms cubic model with fixed transverse width
4. Physical Insight: Reveals systematic influence of confinement geometry parameter $n$ on condensate properties

Limitations

1. Strong Interaction Constraint: Transverse-longitudinal coupling strengthens for large $g$, adiabatic approximation fails, accuracy decreases
2. High-$n$ Quantitative Deviation: Critical interaction strength $g_c$ relative error increases with $n$
3. Integer $n$ Verification: Theory applies to continuous $n$, but numerical validation limited to integer values
4. Single-Component System: Does not address multi-component, spin-orbit coupling, and other complex scenarios

Future Directions

1. Extension to Other Systems:
   • Multi-component BEC
   • Spin-orbit coupling systems
   • Fermi-Bose mixtures
2. Non-Adiabatic Effects: Study impact of transverse excitations in fast dynamical processes
3. Experimental Validation: Comparison with actual BEC experimental data
4. Method Optimization: Improve variational ansatz for enhanced accuracy in strong interaction regime

In-Depth Evaluation

Strengths

1. Methodological Innovation (★★★★★)

• Unified Framework: First continuous description of different confinement geometries through single parameter
• Analytical Tractability: Retains analytical expressions for transverse width, balancing accuracy and operability
• Theoretical Rigor: Clear variational derivation, explicit physical assumptions

2. Experimental Sufficiency (★★★★☆)

• Multi-Dimensional Validation: Static properties (density, chemical potential, $\langle z^2\rangle$) and dynamical properties (scattering)
• Broad Parameter Space: Covers $n=1$ to $6$, interaction strength spanning 3 orders of magnitude
• Quantitative Comparison: Systematic comparison with full 3D-GPE and simplified cubic model
• Limitation: Lacks verification for non-integer $n$, relatively limited dynamical scenarios

3. Result Convincingness (★★★★☆)

• High Precision: Dynamical transmission coefficient relative deviation <0.2%
• Physical Consistency: Correctly predicts collapse phenomena, asymptotic behavior
• Honest Limitations: Clearly indicates accuracy degradation in strong interaction regime

4. Writing Clarity (★★★★★)

• Logical Structure: Theory derivation → numerical validation → dynamical applications, clear progression
• Rich Figures: 8 comprehensive figures displaying results
• Mathematical Expression: Detailed formula derivation, clearly defined notation

Weaknesses

1. Method Limitations

• Adiabatic Approximation Dependence: May fail under strong interactions or rapid dynamics
• Fixed Ansatz Form: Gaussian radial distribution may not suit all scenarios
• Quantitative Accuracy: Larger critical interaction strength deviation at high $n$

2. Experimental Setup

• Limited Dynamical Scenarios: Only barrier scattering examined, missing solitons, vortices, etc.
• Fixed Parameter Choice: $\gamma^2=0.1$ fixed, lacking systematic study of anisotropy parameter effects
• Insufficient Error Analysis: Limited numerical error and convergence analysis

3. Theoretical Analysis

• Incomplete Asymptotic Analysis: Missing detailed discussion of $n\to 0$ and $n\to\infty$ limits
• Stability Analysis: Only linear stability judged by imaginary time convergence, lacking Bogoliubov spectrum analysis
• Energy Functional: Insufficient discussion of effective energy functional properties

Impact

1. Contribution to Field (★★★★☆)

• Theoretical Tool: Provides flexible dimensional reduction framework for BEC community
• Methodological Inspiration: Parametric unification approach generalizable to other physical systems
• Experimental Guidance: Helps understand non-ideal confinement potential effects on BEC

2. Practical Value (★★★★☆)

• Computational Efficiency: 1D equation orders of magnitude faster than 3D-GPE
• Parameter Scanning: Facilitates systematic study of confinement geometry effects
• Experimental Design: Guides optimization of trap parameters

3. Reproducibility (★★★★☆)

• Strengths: Clear equations, standard numerical methods (Split-Step)
• Weaknesses: No public code, some parameters (grid resolution) not specified

Applicable Scenarios

Most Suitable

1. Strongly Anisotropic BEC: $\gamma \ll 1$ cigar-shaped, funnel-shaped geometries
2. Weak to Moderate Interactions: $|g| \lesssim 10$ regime
3. Adiabatic Dynamics: Processes with time scales far exceeding transverse oscillation periods
4. Parameter Optimization Studies: Scenarios requiring systematic confinement geometry parameter scanning

Unsuitable

1. Strong Interaction Regime: Accuracy decreases for $|g| \gg 10$
2. Weakly Anisotropic Systems: $\gamma \sim 1$ requires 2D or 3D treatment
3. Fast Non-Adiabatic Processes: Significant transverse excitation dynamics
4. Complex Many-Body Effects: Beyond mean-field description

Potential Improvements

1. Higher-Order Corrections: Perturbative corrections including transverse excitations
2. Adaptive Ansatz: Dynamically adjust radial distribution based on local density
3. Machine Learning Assistance: Neural networks fitting optimal variational fields
4. Experimental Data Fitting: Reverse-engineer effective parameters from measured density distributions

Selected References

Foundational Works

• 40 Salasnich et al., Phys. Rev. A 65, 043614 (2002): Original 1D-NPSE for cigar-shaped BEC
• 39 Jackson et al., Phys. Rev. A 58, 2417 (1998): Adiabatic approximation methods
• 65 dos Santos et al., J. Phys. B 52, 245301 (2019): Funnel-shaped geometry

Methodology Related

• 51 Mateo & Delgado, Ann. Phys. 324, 709 (2009): Dimensional reduction for vortex dynamics
• 68 Salasnich et al., Phys. Rev. A 66, 043603 (2002): Adiabatic approximation theoretical foundation
• 69 Yang, Nonlinear Waves (2010): Split-Step numerical method

Experimental Background

• 5 Anderson et al., Science 269, 198 (1995): First BEC realization (87Rb)
• 70 Nguyen et al., Nat. Phys. 10, 918 (2014): Barrier scattering experiments

Overall Assessment

This is a high-quality theoretical physics paper making substantial contributions to BEC dimensional reduction. By introducing parameterized generalized transverse potential, the authors successfully establish a unified theoretical framework describing different confinement geometries, solving the problem of requiring separate derivations for each potential in existing methods. The variational derivation is rigorous, numerical validation comprehensive, and results convincing.

Main Highlights: (1) High theoretical innovation, (2) Tight integration of analytical and numerical methods, (3) Clear physical picture. Main Limitations: (1) Limited accuracy in strong interaction regime, (2) Relatively limited dynamical scenarios, (3) Lack of non-integer $n$ verification.

This work provides practical theoretical tools for the BEC community, with expected impact in cold atom physics and nonlinear optics. Subsequent research should extend to multi-component systems and non-adiabatic dynamics.

Recommendation Index: ★★★★☆ (4.5/5)