Abstract:We investigate the use of the sinc collocation and harmonic oscillator bases for solving a two-particle system bound by a Gaussian potential described by the radial Schrödinger equation. We analyze the properties of the bound state wave functions by investigating where the basis-state wave functions break down and relate the breakdowns to the infrared and ultraviolet scales for both bases. We propose a correction for the asymptotic infrared region, the long range tails of the wave functions. We compare the calculated bound state eigenvalues and mean square radii obtained within the two bases. From the trends in the numerical results, we identify the advantages and disadvantages of the two bases. We find that the sinc basis performs better in our implementation for accurately computing both the deeply- and weakly-bound states whereas the harmonic oscillator basis is more convenient since the basis-state wave functions are orthogonal and maintain the same mathematical structure in both position and momentum space. These mathematical properties of the harmonic oscillator basis are especially advantageous in problems where one employs both position and momentum space. The main disadvantage of the harmonic oscillator basis as illustrated in this work is the large basis space size required to obtain accurate results simultaneously for deeply- and weakly-bound states. The main disadvantage of the sinc basis could be the numerical challenges for its implementation in a many-body application.

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to compare the effectiveness of using the sinc - configuration basis and the harmonic oscillator basis when solving the bound states of a two - particle system described by the Gaussian potential. Specifically, the researchers analyzed the performance of these two basis functions in the radial Schrödinger equation, explored their advantages and limitations at different energy scales, and proposed a correction method to improve the behavior of the long - range tail. By calculating the eigenvalues and mean - square radii of the bound states, the researchers aimed to identify the applicable ranges of the two basis functions and the challenges in their numerical implementation. ### Background and Objectives of the Paper - **Research Background**: In nuclear structure theory, the three - dimensional harmonic oscillator basis is widely used because it retains rotational symmetry, simplifies the implementation of conservation laws, and has equivalent treatment in position space and momentum space. However, the asymptotic behavior (Gaussian tail) of the harmonic oscillator basis does not match the exponential tail of the bound - state wave function, which limits its application in some problems. - **Research Objectives**: This paper aims to compare the performance of the sinc basis and the harmonic oscillator basis when solving the bound states of a two - particle system described by the Gaussian potential, especially focusing on the characteristics of these basis functions on the infrared (IR) and ultraviolet (UV) scales, and their accuracy in calculating deeply - bound states and weakly - bound states. ### Main Research Contents 1. **Theoretical Framework**: - **Radial Schrödinger Equation**: The researchers started from the radial Schrödinger equation and derived the matrix representations in the sinc basis and the harmonic oscillator basis. - **Selection of Basis Functions**: The mathematical forms of the sinc basis and the harmonic oscillator basis were introduced, including their orthogonality and normalization conditions. - **Characteristic Scales**: The infrared and ultraviolet characteristic scales related to the basis functions were defined, and how these scales affect the accuracy of the wave function was discussed. 2. **Numerical Methods**: - **Harmonic Oscillator Basis**: How to construct the Hamiltonian matrix in the harmonic oscillator basis was introduced in detail, and the eigenvalue problem was solved by numerical diagonalization. - **Sinc Basis**: The configuration method of the sinc basis was introduced, including choosing appropriate mapping functions and step - size parameters to adapt to the asymptotic behavior of the wave function. 3. **Results and Comparisons**: - **Eigenvalues and Mean - Square Radii**: The eigenvalues and mean - square radii of the bound states under different basis functions were calculated, and their accuracies were compared. - **Wave - Function Structure**: By plotting the wave - function diagrams, the performance of the sinc basis and the harmonic oscillator basis at different scales was analyzed, especially their differences in the long - range and short - range regions. - **Error Analysis**: Two indicators, internal error and inter - basis error, were introduced to evaluate the stability of different basis functions in numerical implementation. ### Main Findings - **Advantages of the Sinc Basis**: When calculating deeply - bound states and weakly - bound states, the sinc basis shows higher accuracy, especially in states where long - range physics is important. - **Advantages of the Harmonic Oscillator Basis**: The basis functions of the harmonic oscillator basis are orthogonal and maintain the same mathematical structure in position space and momentum space, which is very advantageous in problems that require simultaneous processing of position and momentum space. - **Limitations**: The harmonic oscillator basis requires a larger basis - space size when simultaneously and accurately describing deeply - bound states and weakly - bound states, while the numerical implementation of the sinc basis in multi - body problems is more complex. ### Conclusions - **Method Improvement**: A method of improving long - range physics by attaching an exponential tail was proposed, thereby increasing the accuracy of the wave function. - **Future Directions**: The possibilities of further optimizing the method were discussed, including choosing more appropriate basis - function parameters and improving numerical algorithms. Through these studies, the paper provides guidance for choosing basis functions suitable for specific problems and offers valuable references for future multi - body system research.

Comparing Sinc and Harmonic Oscillator Basis for Bound States of a Gaussian Interaction

Two rotating particles interacting via two-body Gaussian potential harmonically confined in two spatial dimensions

Basis Set Generation for Quantum Dynamics Simulations Using Simple Trajectory-Based Methods

Non-Stationary Long-Term Dynamics via Selected Incomplete Dual Bases

Application of a Numerical-Basis-state Method to Strong-Field Excitation and Ionization of Hydrogen Atoms

Efficient Representation of Gaussian Fermionic Pure States in Non-Computational Bases

Three-body Forces in Oscillator Bases Expansion

Self‐Consistent Molecular‐Orbital Methods. IX. An Extended Gaussian‐Type Basis for Molecular‐Orbital Studies of Organic Molecules

Spectral Equivalence of Bosons and Fermions in One-Dimensional Harmonic Potentials

Time correlation functions for quantum systems: validating Bayesian approaches for harmonic oscillators and beyond

Two-center harmonic oscillator basis for Skyrme-DFT calculations (I): formalism and Proof of Principle

On a Repulsive Short-Range Potential Influence on the Harmonic Oscillator

The Harmonic Oscillator Potential Perturbed by a Combination of Linear and Non-linear Dirac Delta Interactions with Application to Bose-Einstein Condensation

Assessing the Accuracy of Quantum Dynamics Performed in the Time-Dependent Basis Representation

Hyperspherical explicitly correlated Gaussian approach for few-body systems with finite angular momentum

Barrier penetration in a discrete-basis formalism

All-electron periodic G 0 W 0 implementation with numerical atomic orbital basis functions: Algorithm and benchmarks

Smooth, invariant orthonormal basis for singular potential Schroedinger operators

Finite Space Complete Basis Method: Precision Computation of High-Resolution Spectrum Near Ionization Threshold

Necessity of orthogonal basis vectors for the two-anyon problem in one-dimensional lattice