Next Article in Journal
Modifications for the Differential Evolution Algorithm
Previous Article in Journal
RGB Image Encryption through Cellular Automata, S-Box and the Lorenz System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Algebraic DVR Approaches Applied to Piecewise Potentials: Symmetry and Degeneracy

by
Eduardo Suárez
1,
Rubén Dario Santiago-Acosta
2 and
Renato Lemus
1,*
1
Instituto de Ciencias Nucleares, UNAM A.P. 70-543, Circuito Exterior, C. U., Mexico City 04510, Mexico
2
Tecnológico de Monterrey, Escuela de Ingeniería y Ciencias, Carr. al Lago de Guadalupe Km. 3.5, Monterrey 64849, Mexico
*
Author to whom correspondence should be addressed.
Symmetry 2022, 14(3), 445; https://doi.org/10.3390/sym14030445
Submission received: 19 December 2021 / Revised: 21 January 2022 / Accepted: 16 February 2022 / Published: 23 February 2022
(This article belongs to the Section Mathematics)

Abstract

:
Algebraic discrete variable representation (DVR) methods that have been recently proposed are applied to describe 1D and 2D piecewise potentials. First, it is shown that it is possible to use a DVR approach to describe 1D square well potentials testing the wave functions with exact results. Thereafter, Morse and Pöschl-Teller (PT) potentials are described with multistep piecewise potentials in order to explore the sensibility of the potential to reproduce the transition from a pure square well energy pattern to an anharmonic energy spectrum. Once the properties of the different algebraic DVR approaches are identified, the 2D square potential as a function of the potential depth is studied. We show that this system displays natural degeneracy, accidental degeneracy and systematic accidental degeneracy. The latter appears only for a confined potential, where the symmetry group is identified and irreducible representations are constructed. One particle confined in a rectangular well potential with commensurate sides is also analyzed. It is proved that the systematic accidental degeneracy appearing in this system is removed for finite potential depth.

1. Introduction

The quantum box, also known as the 1D infinite square well, is the simplest quantum mechanical system used in textbooks to exemplify the solution of the Schrödinger equation. This system has only bound states and its generalization to 2D and 3D Cartesian infinite square well potentials is straightforward [1,2,3,4]. Although the energy spectra are indeed easy to obtain, the 2D and 3D cases present a degeneracy pattern not explained by the geometrical point group associated with square and cubic geometrical symmetries; the degeneracy degree is expected to correspond to the dimension of one of the irreducible representations (irrep) of the symmetry group. In situations where the latter is not fulfilled, degeneracy is identified as accidental. However the appearance of systematic accidental degeneracy signals the presence of hidden symmetries, a fact that always deserves further physical insight [5]. In such cases, a new symmetry group is expected to exist where the accidental degeneracy is rendered normal. The non-relativistic hydrogen atom is representative of this situation where the S O ( 3 ) geometrical symmetry should be extended to a four-dimensional S O ( 4 ) group, in which its generators correspond to angular momentum and Runge–Lenz components [6,7,8]. Surprisingly, the considerably more simple system of one particle in an impenetrable 2D square box presents a similar situation where the geometrical symmetry group C 4 v had to be extended to the semidirect group G = T C 4 v , where T is a one-dimensional compact continuous group [9]. The fact that this group contains discrete as well as continuum elements makes the construction of its irreps a non trivial task, which from a formal point of view has a close analogy with the construction of irreps for space groups [10]. A similar situation is present for the case of a cubic box, where the extended group turns out to be G = T O h , where T is a two dimensional compact continuous group for which generators correspond to linear combinations of one-dimensional Hamiltonians [11].
An interesting situation emerges when the potential depth becomes finite either for the square or cubic well potentials. In this study, we shall prove that, in the high region of the spectrum, degenerate levels split to render normal degeneracy with respect to the geometrical symmetry group, albeit with the possibility of remaining an accidental degeneracy as energy decreases. Hence, the symmetry group becomes the expected geometrical group C 4 v or O h for the square and cubic potentials, respectively, although with possible presence of accidental degeneracy, which is completely removed by enough decreases in the potential depth. This behavior is also present for a particle in an impenetrable rectangular well potential with commensurate sides [12].
Although analyses of the square, rectangular and cubic impenetrable well potentials have been presented in a more or less degree of formality [9,11,12], their study for finite wells has not been worked out. One of the reasons for this may be that the solutions cannot be obtained in a straightforward manner but also because their study requires the use of the machinery of group representation theory, a background beyond the scope of standard textbooks. In contrast, 1D piecewise constant potentials do not present degeneracy for the bound states, an ideal situation to simplify the analytic manipulations for illuminating features of both the postulates and properties of eigenfunctions; the study of piecewise potentials provides the background to estimate, in a qualitative form, the behavior of wave functions for arbitrary potentials. In particular, a half-square well potential, symmetric square well and the delta function potentials provide common systems in textbooks to illustrate the concepts of quantum mechanic. Even though the general solutions in each region of the potential are easily obtained, the eigenvalues are determined only by solving transcendental equations after imposing the matching conditions in wave functions and derivatives. When the piecewise potential presents two or three regions, this approach is feasible, but as regions increases, the procedure become somewhat cumbersome in 1D systems, but much worse in 2D and 3D potentials. Hence, alternative approaches that avoid the explicit matching of the wave functions and their derivatives are welcome.
A powerful method that allows the solutions of piecewise potentials to be obtained is the discrete variable representation approach (DVR). The feature of this approach is that it provides a discrete basis where the coordinate is diagonal and, consequently, the representation of the potential too. The DVR method has a long story, with original ideas proposed in the 1960s [13,14] and becoming widely used in the 1980s with the works of J.C. Light and others [15,16,17,18,19,20,21,22,23]. However, this method has been independently developed with different names: quadrature discretization method [24,25,26,27,28,29], Lagrange mesh method [30,31,32,33] and configuration localized states [34,35,36]. The main ingredient of these methods is the use of orthonormal polynomials with the Gaussian quadrature method to establish a grid associated with the zeros of the polynomials. The common features and differences between the DVR methods are provided by Baye [33]. The fundamentals of the DVR methods can be identified in 1D systems, but its extension to higher dimensions was immediate. However it is convenient to highlight two routes to extend the method to higher dimensions. One possibility consists in considering systems that intrinsically are defined in 3D, such as the solution of the non-relativistic hydrogen atom. Another possibility for proceeding to higher dimensions consists in taking direct products of 1D systems. This situation is present in the description of the vibrational degrees of freedom where for each internal coordinate an oscillator is associated [37], although the direct product has also been employed to determine ro-vibrational spectra [38]. The direct product basis has the advantage of simplicity, but it has the disadvantage of the large basis dimension to obtain converged energy levels. It is possible to overcome this problem at least in part by the use of the Lanczos algorithm, which does not require storage of the Hamiltonian matrix [39,40], a method that combined with a procedure to exclude functions that do not significantly contribute to the physical description allows the computational cost to be reduced [41].
Recently, an alternative DVR method based on dynamical groups U ( n + 1 ) for n D systems was proposed. These groups emerge when the number of bosons n ^ associated with the harmonic oscillator is constrained by introducing an extra scalar boson s ( s ) , imposing a constant total of number of bosons N ^ = n ^ + s s [42]. In this manner, number N characterizes the Hilbert space where the physical description is carried out. This scheme to construct a dynamical group has two advantages; on one hand, the dynamical group is identified with a compact unitary group and, on the other hand, three physical group chains are identified. The key to derive a DVR method relies on the connection of the U ( n + 1 ) algebraic scheme with the configuration space in the N large limit. This connection leads us to associate the U ( n ) chain with energy representation, the S O ( n + 1 ) chain with coordinate representation and the S O ¯ ( n + 1 ) chain with momentum representation. The diagonalization of the coordinate and momentum in the U ( n ) basis provides bases where the coordinate and momentum are diagonal. In this scheme, it is through the transformation brackets connecting the different bases that makes it possible to obtain the representation of the Hamiltonian matrix in a simple form in terms of diagonal matrices, a feature characterizing a DVR approach. This approach has been called the unitary group approach (UGA), and it has been presented in detail for n D potentials [43,44,45,46,47,48]. A natural alternative to this approach consists in projecting the harmonic oscillator basis to a subspace of finite dimension. The discrete variable representation —called the HO-DVR approach due to the fact that it is purely defined in the space of Harmonic oscillator wave functions— is obtained through the diagonalization of the algebraic representation of operators r 2 and p 2 . The resulting discrete basis exactly corresponds to localized states related with the zeros of the polynomials associated with the solutions of the harmonic oscillator. Both the U ( n + 1 ) -UGA and HO-DVR methods are defined in terms of the harmonic oscillator basis, but they provide different discrete representation bases and consequently one method may be more suitable for a specific problem [49]. In fact, in the description of the Stark effect in the non-relativistic Hydrogen atom, U ( 4 ) -UGA turns out to be more suitable for describing the breaking of spherical symmetry due to the electric field [50], while in all tested potentials where the angular momentum is conserved, HO-DVR offers better convergence [49].
Although the Harmonic oscillator basis has proved to be useful in applications from atoms to quarks [51], it is not expected to be suitable for any potential. A basis to be used in a given problem depends on the boundary conditions and, consequently, the alternative algebraic DVR methods associated with different bases are needed in order to improve convergence. In this venue, two more algebraic DVR approaches for 1D systems have been proposed, one associated with the Morse potential (M-DVR) and the other one to the Pöschl-Teller potential (PT-DVR), both based on the realization of the coordinate and momentum in terms of the generators of the U ( 1 , 1 ) dynamical group [49]. As expected, the M-DVR approach is more appropriate for asymmetric potentials while the PT-DVR method represents an improvement to symmetric potentials where bound and continuous energies are present. These methods, although constrained to 1D systems, provide useful bases to deal with multidimensional problems such as n D square well potentials with finite wells.
In this contribution, we present a study of one particle under 1D and 2D piecewise potentials. This work has two goals. First, we are interested in studying convergence in these kind of systems. This is a nontrivial question given the fact that potentials that are not analytic do not provide exponential convergence [52]. On the other hand, the algebraic DVR approach provide us with a simple and powerful method to study symmetry breaking T C 4 v C 4 v appearing when the walls of a confined particle in a 2D square well potential become finite. We intend to show that this simple system exemplifies the concepts of natural, accidental and systematic accidental degeneracy, which constitute fundamental concepts in establishing a connection between symmetry and degeneracy. In addition, we present the study of a particle confined in a rectangular well potential with commensurate sides, a system where a systematic accidental degeneracy also appears, and it is removed when the depth of the potential becomes finite.
This article is organized as follows. Section 2 is devoted to presenting the salient features of 1D-DVR approaches to be used in the description of piecewise potentials, presenting as an example the HO-DVR approach based on the harmonic oscillator basis. Alternative methods are provided in the Appendix A. In Section 3, the analysis of convergence for a symmetric and asymmetric well potentials using the different 1D algebraic DVR approaches is presented. In the same section, the transition from the spectrum characterizing the square well potentials to an anharmonic spectrum provided by PT and Morse potentials is also analyzed. In Section 4, the analysis of the energy spectrum for square well potential as a function of the potential depth in terms of the direct product of 1D systems is presented. Section 5 is devoted to presenting a similar analysis for a particle in a rectangular well potential with commensurate sides. Finally, in Section 6, the summary and our conclusions are presented.

2. Algebraic DVR Methods in 1D Systems

The aim of this section is to provide the basic ingredients characterizing algebraic DVR methods that will be applied to obtain the solutions of systems defined in terms of piecewise potentials. We shall take as an example the HO-DVR method. Additional methods known as PT-DVR, M-DVR and S U ( 2 ) -UGA are presented in the Appendix A. The basic idea of these methods consists in diagonalizing the matrix representation of the coordinate and momentum to generate discrete variable representation bases, which allows the matrix representation of the Hamiltonian to be expressed in terms of diagonal matrices for both kinetic and potential terms. The relevant aspect of this approach is that the matrix elements are obtained with purely algebraic means, avoiding the intrinsic numerical problems of integration.
The 1D harmonic oscillator with reduced mass μ and frequency ω is defined through the following Hamiltonian:
H ^ h . o . = 1 2 μ p 2 + ω 2 μ 2 q 2 ,
with eigenfunctions given by the following [1]:
ψ n ( q ) = N n e q 2 2 λ 0 2 H n ( q / λ 0 ) ,
where H n ( x ) stands for Hermite polynomials, with λ 0 = / ( ω μ ) and normalization constant N n = ( n ! 2 n λ 0 π ) 1 / 2 . Introducing the bosonic operators { a , a } with the usual realization:
a = 1 2 1 λ 0 q i λ 0 p ; a = 1 2 1 λ 0 q + i λ 0 p ,
the Hamiltonian takes the following simple form:
H ^ F o c k h . o . = ω ( n ^ + 1 / 2 ) ; n ^ = a a ,
with eigenkets [1]
| n = 1 n ! ( a ) n | 0 .
Using this algebraic representation, the matrix elements for the coordinate and momenta are easily obtained
n | q | n = λ 0 2 n + 1 δ n , n + 1 + n δ n , n 1 ,
n | p | n = i 2 λ 0 n + 1 δ n , n + 1 n δ n , n 1 .
If we now consider the projection to a finite space of the following:
L N h . o . = { | n , n = 0 , 1 , , N 1 } ,
and carry out their diagonalization, we obtain the discrete representation of the coordinate and momenta given by the following eigenvectors:
| q i = n = 0 N 1 n | q i | n ,
| p i = n = 0 N 1 n | p i | n ,
where matrices T = | | n | q i | | and W = | | n | p i | | correspond to the coefficients that define the basis transformations (9) and (10) where coordinates and momentum are diagonal.
q | q i = q i | q i , p | p i = p i | p i .
Equation (11) establishes the discrete variable representation characterized by the following:
q j | V ( q ) | q i = V ( q i ) δ i , j , p j | G ( p ) | p i = G ( p i ) δ i , j
for any function of either coordinates V ( q ) or momentum G ( p ) . The matrix representation of a Hamiltonian associated with a general potential V ( q ) in the harmonic oscillator basis thus takes the following form:
H = W Λ ( p ) W + T Λ ( q ) T ,
where the diagonal matrices are given by | | Λ ( p ) | | = ( p i 2 / 2 μ ) δ i j and | | Λ ( q ) | | = V ( q i ) δ i j in momentum and coordinate representations, respectively. In Appendix A, alternative methods following the same idea are presented.

3. 1D Square Well Potentials

Here, we start considering a particle with mass μ inside a symmetric square well potential defined as follows:
V ( x ) = { 0 ; x [ a / 2 , a / 2 ] V 0 ; x > a / 2 , x < a / 2 ,
and depicted in Figure 1. Because of the invariance of the potential under inversion, the solutions are expected to have good parity. Following the usual procedure of imposing the boundary matching relations of both continuity and derivatives for the functions at the limits of the regions, we have the following for the even functions:
Ψ I e ( x ) = A e α x , Ψ I I e ( x ) = A e α a / 2 sec ( k a / 2 ) cos ( k x ) , Ψ I I I e ( x ) = A e α x ,
with the definitions α 2 = 2 μ E / 2 and k 2 = 2 μ V o / 2 α 2 . The energies are obtained by the following transcendental equation.
tan ( k a / 2 ) = α / k .
On the other hand, for the odd functions, we have the following:
Ψ I o ( x ) = A e α x , Ψ I I o ( x ) = A e α a / 2 csc ( k a / 2 ) sin ( k x ) , Ψ I I I o ( x ) = A e α x ,
with the corresponding transcendental equation for energies
tan ( k a / 2 ) = k / α .
In both cases, amplitude A is determined by normalization. A particular case of academic interest is the confined particle corresponding to infinite potential barriers. In this case, the eigenvalues take the following form:
E n = 2 π 2 2 μ a 2 n 2 ,
with eigenfunctions carrying good parity:
Ψ n ( e ) ( x ) = 2 a cos n π a x ; n = 1 , 3 , 5 ,
Ψ n ( o ) ( x ) = 2 a sin n π a x ; n = 2 , 4 , 6 ,
for even and odd parity, respectively.
The procedure to obtain the solutions for the potential (14), although relatively easy to obtain, has two inconveniences: On one hand, energies are obtained by transcendental equations and, on the other hand, this approach becomes very cumbersome when additional piecewise steps in the potential are involved. We now proceed to obtain the bound states of potential (14) using algebraic DVR approaches presented in Appendix A. Because of the symmetry of the potential, the appropriate methods satisfying the boundary conditions of potential (14) are the HO-DVR, PT-DVR and S U ( 2 ) -UGA approaches. The Hamiltonian matrix representation in HO-DVR and S U ( 2 ) -UGA-DVR approaches takes the general form:
H = Λ ( E ) + T Λ ( q ) T ,
where Λ ( E ) is the diagonal contribution of the harmonic oscillator contribution with an explicit expression depending on the method in the following form:
HO−DVR : | | Λ ( E ) | | = ω n + 1 2 δ n n ,
U ( 2 ) UGA : | | Λ ( E ) | | = ω n + 1 2 ϵ n 2 N δ n n ,
where ϵ = 1 for n < N / 2 and ϵ = 0 for n N / 2 in accordance with (A36). The diagonal matrix Λ ( q ) corresponds to the modified potential in the coordinate representation given by the following.
HO−DVR : | | Λ ( q ) | | = ω V ( x i ) 1 2 μ ω 2 q i 2 δ i j ,
U ( 2 ) UGA : | | Λ ( q ) | | = ω V ζ 2 N 1 N μ ω 2 ζ 2 δ ζ ζ .
Accordingly the transformation coefficients involved in (20) are defined as follows.
HO−DVR : T = | | n | q i | | ; E q u a t i o n ( 9 )
U ( 2 ) UGA : T = | | [ N ] n | [ N ] ζ | | ; E q u a t i o n ( A 34 ) .
Regarding the PT-DVR method, the Hamiltonian matrix representation takes the form (13), where we have for the diagonal matrix Λ ( p ) in the momentum representation, | | Λ ( p ) | | = p i 2 2 μ δ i j , while for the diagonal matrix for the potential in the position representation, we have | | Λ ( q ) | | = V ( q i ) δ i j . In this case, both transformation coefficients defined in (A15) are needed. In every case, the potential in discretized representation V ( x i ) corresponds to the square well potential defined in (14). Diagonalization is carried out in a truncated basis specified by N, the dimension of the finite projected space. Each dimension space provides a particular set of discretized points that are distributed along the piecewise potential in an inhomogeneous form. This feature does not allows a smooth convergence, as we shall see.
In order to test the methods, we start the application of DVR methods to the square well potential with infinite walls, which in practice were taken to have a height of V 0 = 10 8 cm 1 with the length of the box at a = 4 Å. In all cases of this contribution, the mass was taken to be μ = 1.6726219 × 10 27 kg, an approximation to the proton mass. The comparisons between the spectrum generated by the three approaches (HO-DVR, PT-DVR and S U ( 2 ) -UGA) and the exact energies are displayed in Table 1. The criterion to obtain convergence in energies for HO-DVR and PT-DVR methods was to reach less than 1% of errors in all states. Both methods are similar in both quality and basis dimension. However, S U ( 2 ) -UGA does not provide reasonable results even when considering a much larger basis of N = 2483 . Increasing the basis dimension does not result in a stable calculation, albeit the chosen basis dimension corresponds to a minimum in the rms. Since S U ( 2 ) -UGA provides a discretized basis but is not localized, these results mean that localized discrete bases are indeed the best option to describe piecewise potentials.
We now consider the case of a 1D square well potential with finite potential V 0 = 1000 cm 1 and length a = 4 Å. This system presents 10 bound states that can be obtained in exact form by solving transcendental Equations (15) and (16). The comparisons between the exact energies and the ones obtained using DVR methods are provided in Table 2. As expected, for S U ( 2 ) -UGA, a much larger basis is needed to obtain comparable errors. In contrast, the HO-DVR and PT-DVR approaches provide a pretty good description of the system. The criterion to obtain convergence in the energies for both HO-DVR and PT-DVR methods was again to reach less than 1% of error in all states. In order to obtain a more detailed comparison between the HO-DVR and PT-DVR methods, in Figure 2, a plot of the errors is depicted. This shows that the PT-DVR approach provides a better uniform description in most of the states. Only in the last state was the error observed to increase over the error provided by the HO-DVR approach.
Convergence in energy does not imply convergence in wave functions. Wave functions are tested with fidelity. Fidelity F α ( N ) for a given state | Ψ α parameterized with the basis dimension N is defined as the overlap between consecutive eigenstates separated by Δ N :
F α ( N ) = Ψ α N + Δ N | Ψ α N .
For the basis dimension given in Table 2, we have calculated the fidelity for each state taking Δ N = 20 . The results are displayed in Figure 3. The better quality of the wave functions provided by the PT-DVR methods is manifested in the closer values to unity. In order to show the impact of fidelity on the quality of the wave functions, we display in Figure 4 the calculation of wave functions q | Ψ 1 , q | Ψ 5 and q | Ψ 10 provided by both methods, indicating fidelity. In order to identify the best approach, at the right of the wave functions, the absolute value of the difference between the exact and calculated wave function is also shown. From these results, it is clear that PT-DVR provides a better description in both energies and wave functions. This is a somewhat expected result, given the fact that the PT basis carries information about the continuum. From these results, it is clear that PT-DVR provides a better description in both energies and wave functions.
One feature of this system is the peculiar N-dependence of convergence. In Figure 5, we present as an example the error dependence of basis dimension N for state | Ψ 5 . The oscillations become evident as well as the fact that each method displays a characteristic pattern. We thus have to be careful in selecting the appropriate dimension N: it may be the case that increasing the basis does not imply a better convergence since, in these systems, an exponential convergence is not present [52], and a DVR method does not imply an variational approach [53]. Hence, in order to take full advantage of the DVR methods, a study of the errors should be conducted. It is only through this analysis that N can be appropriately chosen to obtain the best description.

Multistep Piecewise Potentials

The potential displayed in Figure 1 may be compared with a realistic PT-potential. In a similar manner, the half-space square well potential may be compared with the Morse potential, as shown in Figure 6. This comparison is based on the same boundary conditions as well as the presence of continuum states. However, the corresponding spectra are completely different: While PT and Morse potentials have an anharmonic spectrum characterized by a negative slope in a Birge–Sponer plot involving difference E v + 1 E v = ω e 2 x e ω e ( v + 1 ) vs. number of state, the square well potentials present a positive slope α in E n + 1 E n = α ( 2 n + 1 ) . This fact results in the conclusion that anharmonicity is not only a consequence of the presence of the continuum but is fundamental result of the form of the potential curve in particular regions. In order to elucidate the region of the potential responsable to the anharmonic behavior, we propose to represent both the PT and Morse potentials in terms of multipiecewise potentials, as displayed in Figure 7, where in the limit of infinite steps, true potentials are expected to be recovered. Our goal is to use the DVR approach to study the number of steps (degree of smoothness) to reproduce the anharmonicity pattern characterizing the molecular spectra.
First, we consider the PT potential approximated by a piecewise potential following the potential curve. We have divided the range of the potential in equals steps, a partition that implies a better approximation for the potential near the continuum. In Figure 8, the energy correlation as a function of the number of steps using both HO-DVR and PT-DVR methods is displayed. In these plots, the left limit corresponds to the exact energy spectrum for the square well potential, while the right limit corresponds to the exact PT spectrum. The first feature we notice is that the anharmonic energy pattern is reached with relatively few steps. As a consequence the anharmonic spectrum is determined by the form of the potential near the continuum. This analysis resembles the notch test used to establish the confident limits of different regions of the potential [48].
Although the general correlation pattern displayed in Figure 8 by both methods is the same, there are some subtle differences concerned with the smoothness of the oscillations associated with the convergence previously discussed and depicted in Figure 5. Beyond this subtle differences, it is clear that there is an advantage of the PT-DVR approach over the HO-DVR method.
We have also carried out the corresponding analysis by approximating the Morse potential with a piecewise potential using the M-DVR approach. The corresponding energy correlation diagram is displayed in Figure 9. Again, the anharmonic pattern is reached with relatively few steps due to the fact that the approximation displayed in Figure 7 implies a better approximation near the continuum. It is worth stressing the fact that the description of the bound states for the half-space square well potential strongly depends on the basis dimension. In order to show this dependence, we present two diagrams corresponding to N = 150 and N = 2500 , respectively, where in the latter case the best description of the half-space square well potential becomes. The number of bound states for the half-space square well potential is much greater that the corresponding Morse spectrum. As the number of steps increases to approximate the Morse potential, most of the bound states move to the discretization description of the continuum.

4. 2D-Square Well Potential

The square 2D potential, although seldom discussed in detail in textbooks, represents a quite interesting problem from the symmetry point of view [9]. The solutions for the confined case with infinite walls are straightforwardly obtained in terms of the direct product of 1D solutions with energies.
E n 1 n 2 = 2 π 2 2 μ L 2 ( n 1 2 + n 2 2 ) .
This solution presents a degeneracy far from being trivial to be explained. In fact, from the point of view of the geometrical group C 4 v , the eigenstates present a systematic accidental degeneracy represented by the couple of degenerated levels carrying one-dimensional irreducible representations (irreps), as it is shown in the energy diagram depicted in Figure 10, where energy units were taken to be E ¯ n 1 n 2 = E n 1 n 2 / ( 2 π 2 / 2 μ L 2 ) . In order to explained this degeneracy, the new symmetry group G = T C 4 v was proposed where the accidental degeneracy with respect to the C 4 v group renders normal in the new group [9]. This peculiar property is particular to the confined system. When the walls become finite, degeneracy is either partially or completely removed depending of the depth and wide of the well. However, the description of the 2D system for a finite wall is a nontrivial task that will be considered in this work. The aim of this study is to provide physical insight into the symmetry breaking phenomenon induced by the depth of the potential. Starting with a confined particle with symmetry group G = T C 4 v , which is analyzed in detail in Appendix B, we shall proceed to present the general approach to tackle the 2D system in the framework of algebraic DVR approaches. We strongly recommend reading Appendix B in order to take full advantage of the discussion of this section.

4.1. Application of the Algebraic DVR Approaches to 2D Systems

Here, we consider a particle in a square well potential with potential depth V 0 . This problem will be solved in terms of the direct product of 1D bases corresponding to the x and y axes. In the following discussion, we shall use the HO-DVR method. The procedure using the PT-DVR approach is analogous, and we only have to establish the following mapping: | n x | n y | Φ n x σ | Φ n y σ .
The Hamiltonian associated with 2D square potential will be diagonalized in the direct product space of harmonic oscillators given by the following.
| n x n y = | n x | n y .
In order to generate the DVR basis, we first consider the discrete basis associated with the x coordinate, which takes the following form:
| x i = n x = 0 N 1 n x | x i | n x ,
| p x i = n x = 0 N 1 n x | p x i | n x ,
with a similar basis for the y degree of freedom.
| y i = n y = 0 N 1 n y | y i | n y ,
| p y i = n y = 0 N 1 n y | p y i | n y .
In terms of discrete bases, the direct product (29) takes the following form.
| n x n y = i = 1 N j = 1 N x i | n x y j | n y | x i | y j .
It is now convenient to introduce the following notation for the matrix elements involved in the following (34):
A i , n x = x i | n x ,
B j , n y = y j | n y ,
T i j , n x n y = x i | n x y j | n y ,
with the following relation.
T = A B .
In terms of the T matrix, the direct product (29) is recast in the following form:
| α = β = 1 N 2 T β , α | β ,
where | α = | n x | n y and | β = | x i | y j , with the following mapping.
α = ( n x 1 ) N + n y ,
β = ( i 1 ) N + j .
Following, the same procedure for the momentum representation, we define the following:
C i , n x = p x i | n x ,
D j , n y = p y j | n y ,
W i j , n x n y = p x i | n x p y j | n y ,
with the direct product basis:
| α = β = 1 N 2 W γ , α | γ ,
where | α = | n x | n y , | γ = | p x i | p y j , and γ = ( i 1 ) N + j with the following.
W = C D .
Using the direct product bases, the Hamiltonian matrix representation in the framework of the HO-DVR approach takes the following form:
H = W Λ ( p ) W + T Λ ( V ) T ,
where the matrix elements of Λ ( p ) in the momentum representations are as follows:
| | Λ ( p ) | | = 1 2 μ ( p x i 2 + p y j 2 ) δ i i δ j j ,
while, for the potential in the coordinate representation, we have the following.
| | Λ ( V ) | | = V ( x i , y j ) δ i i δ j j .
Hence, the diagonalization of the coordinate and momentum for the 1D systems in x and y provides the discrete localized bases to obtain the Hamiltonian matrix representation in terms of diagonal matrices. The direct product bases possess the relevant property of simplicity; however, the basis dimension increases as N 2 , and it may happen that a very large number of direct product functions are necessary for energy convergence. In order to reduce the cost of explicit diagonalization, the Lanczos algorithm may be used, which does not require storage of the Hamiltonian matrix [39]. In addition, it may be possible to reduce the size of the direct product by excluding functions that do not contribute to the associated wave functions involved in the energy levels of interest [41]. Instead, in this contribution, we are simplifying our analysis by introducing symmetry adapted functions. In the framework of this point of view, the projected basis takes the following form:
| q ϕ γ ( Γ ) = n x , n y S n x n y ; q Γ γ | n x n y ,
where Γ and γ are the irrep and component, respectively, of group C 4 v . Index q refers to the multiplicity of the irreps and it provides the dimension of the matrix representation to be diagonalized. In this context, the Hamiltonian matrix has the following form:
H s = S [ W Λ ( p ) W + T Λ ( V ) T ] S ,
which is in blocked form:
H s = Γ H s Γ ,
where each block is characterized by the one irrep of the group taking the first component for the degenerate irreps. The dimension of each block is given by q ( Γ ) 2 . We should remark that the use of symmetry, beyond the fact that simplifies the calculation, is a compulsory task because we are interested in identifying the degeneracy associated with symmetry and consequently the levels must carry the corresponding irrep. Symmetry projection, however, is not a trivial task due to the large amount of basis functions involved. Hence, we consider it to be convenient to present a highly efficient symmetry projection approach, which is called the eigenfunction method.

4.2. Symmetry Adapted Basis

In this section, we shall present the application of the eigenfunction method to project the direct product basis (29), which consists in diagonalizing a complete set of commuting operators (CSCO) given in terms of a linear combination of classes [54,55,56,57].
The CSCO used to project the basis is defined in accordance with a suitable canonical subgroup chain, which in our case has been chosen to be the following:
C 4 v C s d ,
where the following is the case:
C s d = { E , σ d a } ,
with symmetry elements depicted in panel (b) of Figure A1. Denoting the i-th class of group C 4 v as K i and the classes of the subgroup C s d as k j , we have, for the CSCO called C I I , the following linear combination:
C ^ I I = C ^ I C 4 v + C ^ I C s d ,
where the following is the case:
C ^ I C 4 v = K ^ 2 + 3 K ^ 5 ; C ^ I C s d = k ^ 2 ,
with the following classes:
K ^ 2 = C ^ 4 + C ^ 4 3 ; K ^ 5 = σ ^ d a + σ ^ d b ; k ^ 2 = σ ^ d a ,
following the standard class numbering for character tables. Operators C ^ I C 4 v and C ^ I C s d distinguish the irreps of groups C 4 v and C s d , respectively, while operator C ^ I I distinguishes the irreps and components of group C 4 v . Their simultaneous diagonalization leads to eigenkets | μ ν defined by the eigensystem:
C ^ I C 4 v | μ ν = μ | μ ν ; C ^ I C s d | μ ν = ν | μ ν , C ^ I I | μ ν = ( μ + ν ) | μ ν ,
where μ = μ ( Γ ) and ν = ν ( γ ) are functions of the group and the subgroup characters:
μ = 2 χ 2 ( Γ ) n Γ + 3 χ 5 ( Γ ) n Γ ; ν = χ 2 ( γ ) ,
where γ stands for the irrep of subgroup C s d , n Γ refers to the dimension of the Γ -th irrep and γ D ( Γ ) ( C 4 v ) C s d . Hence, the symmetry-adapted functions are equally expressed either with | Γ γ or | μ ν .
We now have to establish the effect of the operators involved in (55). Since C ^ 4 and σ ^ d a are generators of the group, it will not be necessary to consider operator σ ^ d b . Formally, we start establishing the effect of these operators over the following basis.
C ^ 4 | x = | C 4 x = | y ; σ ^ d a | x = | σ d a x = | y , C ^ 4 | y = | C 4 y = | x ; σ ^ d a | y = | σ d a y = | x .
Hence, for rotation C 4 , we have the following.
n x n y | C ^ 4 | n x n y = 4 d x d y d x d y n x n y | x y x y | C 4 | x y x y | n x n y = 4 d x d y d x d y n x n y | x y x y | n x n y δ ( y x ) δ ( y + x ) = ( 1 ) n x δ n y n x δ n x n y .
In the similar form, we obtain the following.
n x n y | C ^ 4 3 | n x n y = ( 1 ) n y δ n y n x δ n x n y ,
n x n y | σ ^ d a | n x n y = δ n y n x δ n x n y ,
n x n y | σ ^ d b | n x n y = ( 1 ) ( n x + n y ) δ n x n y δ n y n x .
Gathering all these matrix elements in operator C ^ I I , we obtain the following.
n x n y | C ^ I I | n x n y = [ ( 1 ) n x + ( 1 ) n y + 4 + 3 ( 1 ) n x + n y ] δ n x n y δ n y n x .
Although this matrix representation has dimension N 2 , in practice, it is a block diagonal with N one-dimensional matrices and a number of N ( N 1 ) / 2 two-dimensional matrices for even N, for instance. This fact significantly simplifies the projection. For example, taken N = 4 with dimension space N 2 = 16 , the symmetry projection reduces to the following:
L 12 = L 3 A 1 L 1 A 2 L 1 B 1 L 3 B 2 L 4 E 1 ,
where of course only the first component of the degenerate irrep E is necessary. Hence, the original Hamiltonian matrix of dimension N 2 given by (47) reduces to the diagonalization of five matrices of dimension a Γ in accordance to the reduction in representation Δ ( r e d ) ( C I I ) given by the following (65).
S 1 Δ ( r e d ) ( C I I ) S = μ a Γ D Γ ( C I I ) .
We now proceed to take advantage of this projection to obtain the description of the square well potential in the framework of the DVR approach through the diagonalization of the representation matrix (51).

4.3. Results and Discussion

In order to evaluate the validity of our approach, we start considering the case of a square well potential with infinite walls with exact energies given by (28). The parameters were chosen to be a = 4 Å and V 0 = 10 8 cm 1 , with the same mass previously considered. As in the 1D case, we describe the system using both HO-DVR and PT-DVR methods, fitting the parameters ω and α , respectively. The results are displayed in Table 3, where the criterion for the extension of the basis was given by obtaining errors of 1%. We have carried out calculations with the criterion of choosing the basis dimensions as close as possible but in accordance with an absolute minimum in rms. The basis dimensions are N = 52 and N = 60 for the HO-DVR and PT-DVR approaches respectively. The results are evaluated at the level of energies and wave functions. Concerning the energy spectrum, both methods present similar errors.
In order to show that the energy convergence showed in Table 3 provides also the correct wave functions, we have projected the eigenvectors to the position representation. Each eigenvector | ψ α Γ γ of Hamiltonian (51) is given in terms of a linear combination of the symmetry adapted basis (50).
| ψ α Γ , γ = q V q α | q ϕ γ ( Γ ) .
The wave functions are obtained through the following projection:
ψ α Γ , γ ( x , y ) = x y | ψ α Γ , γ = q V q α x y | q ϕ γ ( Γ )
with the following:
x y | q ϕ γ ( Γ ) = n 1 , n 2 S n 1 n 2 ; q Γ γ ψ n x ( x ) ψ n y ( y ) ,
where ψ n x ( x ) and ψ n y ( y ) are 1D harmonic oscillator functions. In Figure 11, a selected set of eigenstates has been depicted. These wave functions are identical to the ones obtained by projecting the exact functions (A52). For the first levels, we have the following correspondence.
| ψ 1 A 1 | ϕ 11 A 1 ; | ψ 1 B 2 | ϕ 22 B 2
| ψ 2 A 1 | ϕ 31 A 1 ; | ψ 2 B 2 | ϕ 42 B 2
| ψ 3 A 1 | ϕ 33 A 1 ; | ψ 1 E , A | ϕ 21 E , A
| ψ 1 A 2 | ϕ 42 A 2 ; | ψ 2 E , A | ϕ 32 E , A
| ψ 1 B 1 | ϕ 31 B 1 ; | ψ 3 E , A | ϕ 41 E , A .
It should be clear that the labels of these wave functions have their correspondence with the labeling of group G in accordance with Table A6.
In Figure 11, we show a set of representative wave functions illustrating the different types of degeneracy. The first row (a) shows the degenerate wave functions (3,1) with irreps ( A 1 , B 1 ) and energy E = 103.24 cm 1 . This degeneracy is accidental from the point of view of the C 4 v group but natural with respect to the G group with irrep k A 1 G . In the second row, (b) we show degenerate wave functions | ψ 3 E , A | ϕ 41 E , A , which in both groups render a natural degeneracy, a feature manifested by the fact that they are connected by rotation C 4 . In (c), we show the degenerate wave functions ( 4 , 2 ) corresponding to | ψ 1 A 2 , | ψ 2 B 2 with energy E = 206.48 cm 1 , which span the two dimensional irrep k A 2 G , but with accidental degeneracy with respect to the C 4 v group. Finally, we show in the last row wave function | ψ 3 A 1 with energy E = 185.84 cm 1 , which span the one dimensional irrep k 0 A 1 G A 1 . On the other hand, the wave functions are compared by plotting the absolute value for the difference between the exact and calculated wave functions displayed in Figure 12, a result that shows that the PT-DVR approach provides better wave functions.
It is interesting to show the N-dependence of the root mean square deviation, which is displayed in Figure 13. In both methods, the general pattern is similar but the detailed structure is quite different. The PT-DVR approach provides a somewhat more regular pattern. In this approach, potential V ( x i , y j ) is evaluated at the discrete points. From the total number of points N 2 , a number of points N i n fall inside the box (region with V 0 = 0 ). We may thus define the following:
% P = N i n N 2 × 100 ,
as the percentage of points falling inside the box. A significant feature of the N-dependence of convergence is that % P is correlated to rms: a minimum in % P corresponds to minimum in rms. This fact has been pointed out with dashed lines in Figure 13. We should observe, however, that several minima in rms are possible for choosing the basis dimension. We have chosen the criterion to select the first absolute minimum, selecting basis dimensions of the same order for both methods. This is a remarkable feature because it provides us with a criterion to select the appropriate basis dimension N to obtain an improvement. This information can be known before carrying out the calculation, since % P is obtained with the zeros of the polynomials [49]. This fact makes DVR methods efficient.
The spectrum depicted in Figure 10 is characteristic of the confined particle. When the potential depth becomes finite, the spectrum is expected to be modified in both energy and degeneracy. The latter feature will be proven by carrying out the calculation with the PT-DVR method.
In Figure 14, we show energy levels as a function of potential depth V 0 . It is worth stressing the splitting of the levels ( A 1 , B 1 ) and ( A 2 , B 2 ) as potential decreases. This means that the true symmetry group becomes geometrical group C 4 v . In order make the symmetry breaking G C 4 v in Figure 15 transparent, a zoom-in of levels ( A 1 , B 1 ) and ( A 2 , B 2 ) is shown. Notice that levels ( A 2 , B 2 ) split sooner than levels ( A 1 , B 1 ) , which is explained by the fact that states { | ϕ 42 A 2 , ϕ 42 B 2 } have higher energies. While the splitting of these levels manifests around V 0 = 450 cm 1 , the splitting of states { | ϕ 31 A 1 , ϕ 31 B 1 } are clearly manifested until V 0 = 250 cm 1 . This means that, at some values of potential V 0 , the systematic accidental degeneracy disappears, but near-degeneracy remains at the low lying region of the spectrum. In order to quantify this statement, we introduce parameter ζ defined as follows [58,59].
ζ = | 2 π arctan ( E 1 E 2 ) ( E 1 + E 2 ) / 2 | .
This parameter measures the relative splitting between a couple of levels E 1 and E 2 . For complete degeneracy ζ = 0 , while as splitting increases ζ 1 . In the same Figure 15, the parameters ζ A 1 B 1 and ζ A 2 B 2 associated with the degenerate levels for the square well potential are depicted as a function of the potential depth. As the potential depth decreases, the first parameter ζ A 2 B 2 starts increasing before ζ A 1 B 1 . It is until the value of V 0 = 350 cm 1 that ζ A 1 B 1 shows splitting, establishing the limit of accidental degeneracy. Using the PT-DVR approach, we have chosen four potential depths, namely V 0 = 10 8 , 350, 250 and 180 cm 1 , to exemplify the order of splitting induced by the finite potential depth. It is clear that, for potential depth V 0 = 350 cm 1 , the levels A 1 and A 2 are nearly degenerate, while the splitting of levels A 2 and B 2 becomes manifest. This degeneracy pattern was displayed with a given set of parameters V 0 and wide a, but with the appropriate parameters involving a larger number of bound states and the lower states A 1 , B 1 are expected to remain degenerate for finite potential depths. In such situations, an accidental degeneracy remains in the low lying are of the spectrum, with broken degeneracy in the rest of the spectrum. This correlation diagram is has been carried out without taking into account the criterion of % P , and consequently, it cannot be considered as a definitive result. However, it provides proof of symmetry breaking.

5. 2D-Rectangular Well Potential

In this section, we pay attention to the case of a rectangular square well potential with commensurate sides, as depicted in Figure 16. This case is interesting because systematic accidental degeneracy is also present. As in the square well potential case, the eigenstates —taking the origin in the left bottom corner of the rectangle with sides L 1 and L 2 — are given by the following:
Φ n 1 n 2 ( x , y ) = 2 L 1 L 2 sin n 1 π x L 1 sin n 2 π y L 2 ,
where n 1 and n 2 are positive integers. The energies levels are given by the following.
E n 1 n 2 = 2 π 2 2 μ n 1 2 L 1 2 + n 2 2 L 2 2 .
The commensurate sides may be imposed by considering the following:
n L 1 = m L 2 = L 0 ,
in such a manner that energy takes the following form.
E n 1 n 2 = 2 π 2 2 μ L 0 2 ( n 2 n 1 2 + m 2 n 2 2 ) .
Hence, states | ψ n 1 n 2 and | ψ n 1 n 2 are degenerate if the following is the case.
( n n 1 ) 2 + ( m n 2 ) 2 = ( n n 1 ) 2 + ( m n 2 ) 2 .
As an example, we shall consider the case where the length in the y direction is twice the length in the x direction, as depicted in Figure 16. In this case, we have n = 2 and m = 1 , which lresults in the following energy expression:
E ¯ n 1 n 2 = E n 1 n 2 E 0 = ( 4 n 1 2 + n 2 2 ) ,
where E 0 = 2 π 2 2 μ L 0 2 . The corresponding spectrum is shown in Figure 17, where the levels with accidental degeneracy have been remarked with broader lines. In this case, the geometrical symmetry is C 2 v with eigenfunctions (78) already spanning the irreps shown in Table 4 in accordance with the character in Table 5. The search of a dynamical symmetry leads to operator F ^ ( A 1 ) , spanning the A 1 irrep:
F ^ ( A 1 ) = 1 4 5 π 2 x 2 2 y 2 ,
which connects states B 2 . In contrast, it is not possible to find an operator D ^ ( B 1 ) from the set of spherical harmonics that preserve the boundary conditions. In Lemus et al. [12], the dynamical operation σ ^ d a was suggested and it connects states A 2 with B 2 . Here, we do not intend to obtain irreps based on this idea. Instead, we address the problem of describing accidental degeneracy using DVR approaches and prove that, for finite depth potential, a accidental degeneracy is removed.
To accomplish this task we proceed in similar form to the square well potential. The Hamiltonian associated with the 2D rectangular well potential will be diagonalized in the direct product basis of harmonic oscillators given by the following.
| n x n y = | n x | n y .
In terms of this basis, the Hamiltonian matrix representation takes the following form:
H = W Λ ( p ) W + T Λ ( V ) T ,
with the same meaning as in the square well potential. Again, the analysis will be simplified by introducing a symmetry adapted functions through (51), a procedure that should be followed in order to identify accidental degeneracy.

5.1. Symmetry Adapted Basis

We shall proceed to apply the eigenfunction method to project the direct product basis (29). Since group C 2 v presents only one-dimensional irreps, it is enough to consider the following combination of classes.
C ^ = 3 C ^ 2 + σ ^ v a .
The diagonalization of operators (87) leads to eigenkets | μ defined by the eigensystem:
C ^ | μ = μ | μ ,
where μ = μ ( Γ ) is a function of group characters:
μ = 3 χ 2 ( Γ ) + χ 3 ( Γ ) ,
where, here, Γ stands for the irrep of the group C 2 v . Hence, symmetry-adapted functions are equally expressed either with | Γ or | μ .
We now have to establish the effect of the operators involved in (87). Formally, we start establishing the effect of these operators over the following basis.
C ^ 2 | x = | C 2 x = | x ; σ ^ v a | x = | σ v a x = | x , C ^ 2 | y = | C 2 y = | y ; σ ^ v a | y = | σ v a y = | y .
Hence, we have the following.
n x n y | C ^ 4 | n x n y = ( 1 ) n x + n y δ n x n x δ n y n y ,
n x n y | σ ^ v a | n x , n y = ( 1 ) n y δ n x n x δ n y n y .
Gathering all these matrix elements in operator C ^ , we obtain the following.
n x n y | C ^ | n x n y = [ 3 ( 1 ) n x + n y + ( 1 ) n y ] δ n x n x δ n y n y .
Notice that this representation matrix is diagonal and, consequently, we use (93) to assign the corresponding irreps to construct the subspaces of symmetry adapted functions.
We now proceed to obtain the description of the rectangular well potential through the diagonalization of the representation matrix (51).

5.2. Results and Discussion

We first consider the case of a rectangular well potential with infinite walls where the exact energies are given by (83). The parameters were chosen to be as follows: a = 1 Å, b = 2 Å and V 0 = 10 8 cm 1 . Again, we described the system with both HO-DVR and PT-DVR methods, fitting parameters ω and α , respectively. In this system, because of the different dimensions of the box, two parameters corresponding to the different directions are free to be fitted. The results are displayed in Table 6, where the basis dimension was chosen in accordance with the criterion of % P to obtain a minimum error of 1%. It is important to remark that the absolute minimum of rms in HO-DVR corresponds to N = 97 , which is much larger than N = 73 . For the sake of comparison, we decided to take the values chosen in Table 6 in order to have similar dimensions for both methods.
From these results, we notice that the convergence of the expected degenerate states is difficult to be obtained although, it is quite clear that using the PT-DVR approach results in a better description: with a lower basis dimension, the expected degenerate values are pretty close. In the same venue, in Figure 18, the root mean square deviation is presented, where a similar behavior as in the square well potential is manifested: % P is correlated with the rms. Indeed the best description is obtained for a minimum % P , a fact that allows the search for the best dimension basis before carrying out the calculation.
In order to prove that the energy convergence shown in Table 6 provides the correct wave functions, we proceeded to project the eigenvectors to the position representation. In Figure 19, the functions associated with accidental degeneracy are shown. These functions coincide with the ones obtained by Cartesian solutions (78) and correspond to the following.
| Φ 14 B 2 ; | Φ 22 A 2 ,
| Φ 16 B 2 ; | Φ 32 B 2 .
It is clear that there is no geometrical transformation connecting these functions. The differential operator (84) connects the degenerate B 2 functions but not the states | ϕ 14 B 2 and | ϕ 22 A 2 . In Lemus et al. [12], it was suggested that operator σ ^ d a corresponding to the reflection with symmetry element depicted in Figure 20 with the previso that the space must be considered with periodic boundary conditions, which is explained by the need to have well defined functions in the original rectangular box.
The spectrum of Figure 17 corresponds to a confined particle. When the potential depth becomes finite, the spectrum is modified in both energy and degeneracy. In Figure 21, we show the energy levels as a function of potential depth V 0 . As the potential decreases, the accidental degeneracy associated with levels { | Φ 16 B 2 , Φ 32 B 2 } as well as with { | Φ 14 B 2 , Φ 22 A 2 } is broken. Again, the true symmetry group for the confined particle becomes geometrical group C 2 v and no degeneracy appears. The splitting of these couple of levels becomes more evident in the zooms depicted in Figure 22. It is clear that the splitting manifested through function ζ is grater for levels { B 2 , B 2 } than for levels { A 2 , B 2 } .
The analysis we have presented for the rectangular well has been limited to ratio L 1 / L 2 = m / n = 1 / 2 . In order to take into account other possibilities, we shall consider case m : n = 1 / p with p = 1 , 2 , 3 , 4 . In order to achieve this goal, we introduce continuous parameter τ defined as the following.
L 2 / L 1 = τ ; τ : [ 0.5 , 4 ] .
In this manner, we not only consider rectangular boxes with ratios L 2 / L 1 = 2 , 3 , 4 but also the case of a square well potential with τ = 1 . In Figure 23, we present the energy levels as a function of parameter τ . The diagram was generated using the PT-DVR approach with parameters N = 60 , α x = 9 × 10 9 m 1 and α y = 4.5 × 10 9 m 1 . The labeling scheme was chosen in accordance with the C 2 v group. For τ = 1 , we have the case of a square box. The degenerate states are pointed out with circles and correspond to the pair of levels E and { A 2 , B 2 } , which subduce to B 1 B 2 and { A 1 , A 1 } , respectively, in C 2 v . For relation τ = 2 , we identify degenerate levels { A 2 , B 2 } and { B 2 , B 2 } , which are previously analyzed and provided in Figure 17. For relation τ = 3 , two pairs of degenerate levels appear, corresponding to the following levels.
| Φ 16 B 2 ; | Φ 23 B 1 ,
| Φ 27 B 1 ; | Φ 32 B 2 .
Finally, for τ = 4 , we have a degeneracy associated to the following states:
| Φ 17 A 1 ; | Φ 21 B 1 ,
| Φ 18 B 2 ; | Φ 24 A 2 ,
which are pointed out in the correlation figure. Hence, the diagram of Figure 23 not only allows us to identify the accidental degeneracy due to the commensurate sides, but also shows symmetry breaking when ratio L 2 / L 1 stops being an integer. It should be clear that the accidental degeneracy with respect to group C 2 v renders it natural with respect to the true symmetry group.

6. Conclusions

In this contribution, we have presented a quantum mechanical analysis of one particle under 1D and 2D piecewise potentials. First, the benchmark problem of a particle confined in a 1D square well potential was studied in the framework of three methods: HO-DVR, PT-DVR and S U ( 2 ) -UGA. It was found that the PT-DVR method is the most appropriate approach, albeit the HO-DVR method also provides reasonable results. S U ( 2 ) -UGA turns out to be inappropriate for these type of potentials. Henceforth, the 1D square well potential with finite walls was studied using these methods. The general trend regarding the advantage of each method remained: A localized discrete basis is the best choice to describe these kind of systems. The next issue we considered was the energy correlation between a spectrum provided by a square and half-space well potentials and the anharmonic spectra associated with realistic interactions such as the case of PT and Morse potentials, respectively. To achieve this task, we have used a multistep piecewise potential. We have shown that the anharmonic spectrum is dominated by the form of the potential near the continuum but not by the presence of the continuum itself. In addition, we confirm the lack of exponential convergence as the basis dimension increases.
As a second issue in our analysis, we have considered a particle confined in a 2D square well potential, a system that is characterized by presenting systematic accidental degeneracy with respect to the geometrical point group C 4 v . The accidental degeneracy is explained by the identification of a new symmetry group, T C 4 v [9]. In contrast to the representation analysis provided by Leybraz et al. [9], a formal procedure to construct its irreducible representations is provided. The identification of group T C 4 v allows us to render systematic accidental degeneracy normal. We have proved that, when the walls of the square well potential become finite, symmetry breaking is present, a phenomenon formally described by the subduction T C 4 v C 4 v . A relevant aspect of this system with finite walls is that partial symmetry breaking may be manifested depending of the depth of the potential; the low lying levels may still keep accidental degeneracy for a given potential V 0 and width a. Regarding convergence, we have found that the oscillations of the rms displayed in Figure 13 are correlated with the number of discretized points x i falling inside the square box. Since this number of points corresponds to the zeroes of either the Hermite or Gegenbauer polynomials, the appropriate selection of the basis dimension can be determined before carrying out the calculation. This fact compensates for the lack of exponential convergence. The comparison between the exact results and the calculations obtained by HO-DVR and PT-DVR methods for the confined particle provides similar results in energy but the PT-DVR approach provides a better description for wave functions.
A square well potential is a quite simple albeit interesting system where all types of degeneracy may be present, a feature that can be used to exemplify the concept of symmetry group and the consequences of not having appropriately identified such group.
It should be stressed that a crucial aspect of this treatment is not only the symmetry projection of the basis to identify degeneracy but also to simplify calculations. The eigenfunction method turns out to be the best option that consists in diagonalizing a complete set of commuting operators defined in terms of a linear combination of classes associated with a chain of groups.
An additional system that we have included in our discussion is the case of a particle confined in a rectangular well potential with commensurate sides. This system also presents systematic accidental degeneracy. We have shown that this degeneracy is removed by the effect of the presence of finite walls. Choosing appropriately the potential depth V 0 and width a, a near-accidental degeneracy may remain for finite walls. We have also shown a correlation of energy levels when, for a given potential depth, the width of the well in the y direction is increased in order to obtain several multiples of the width in the x direction. The depicted correlation plot allows us to identify accidental degeneracy in the form of crossing of levels, a quite enlightening diagram showing the appearance of degeneracy due to commensurability.
Hence, the present analysis shows that algebraic DVR approaches are suitable for dealing with piecewise potentials. In particular, the PT-DVR approach seems to be more appropriate than the HO-DVR method due the fact that the PT-basis carries information of both anharmonicity and the continuum. An analogous situation occurs for the asymmetric square well potential described with the M-DVR approach.
The study we have presented paves the way to explore future applications in the field of quantum chemistry. In particular, the H 2 + molecule is a good candidate because the evaluation of electron–nuclei interactions is provided as an expansion in terms of spherical harmonics with coefficients depending on the radial coordinate, which behaves as a piecewise potential from the point of view of a DVR approach. A similar situation is present in the Helium atom, although in this case it will be necessary to combine this approach with the transformation brackets introduced by M. Moshinsky [51].

Author Contributions

E.S.: Software, formal analysis, investigation and writing—review and editing; R.D.S.-A.: investigation, validation, writing—review and editing; R.L.: writing—original draft, investigation and supervision. All authors have read and agreed to the published version of the manuscript.

Funding

This work is partially supported by DGAPA-UNAM, Mexico, under project IN-212020.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Alternative DVR Methods

In this appendix, we present the main features of the algebraic Pöschl-Teller-DVR (PT-DVR), Morse-DVR (M-DVR) and the S U ( 2 ) -UGA approaches. With the exception of M-DVR approach, these methods are appropriate for symmetric potentials. We start with the PT-DVR approach.

Appendix A.1. PT-DVR Method

The bound states of the Pöschl-Teller potential do not form a complete set of states in the Hilbert space [60,61,62]. The DVR approach needs a complete basis, which is provided by the following [63]:
Φ n σ ( u ) = A n σ ( 1 u 2 ) σ 2 C n σ 1 / 2 ( u ) ; n = 0 , 1 , 2 , ,
where C n σ 1 / 2 ( u ) are the Gegenbauer polynomials [64] with argument u = tanh ( α q ) , and σ is taken to be σ = 1 in order to decouple the bound from the continuum spectrum [63], while the normalization constant is as follows.
A n σ = α n ! ( σ + n 1 / 2 ) ( Γ [ σ 1 / 2 ] ) 2 π 2 2 2 σ Γ [ 2 σ + n 1 ] .
Using the factorization method, it is possible to obtain ladder operators { K ^ ± , K ^ 0 } with the following effect [63]:
K ^ + Φ n σ ( u ) = k + Φ n + 1 σ ( u ) ,
K ^ Φ n σ ( u ) = k Φ n 1 σ ( u ) ,
K ^ 0 Φ n σ ( u ) = k 0 Φ n σ ( u ) ,
with the following being the case.
k + = ( n + 1 ) ( 2 σ + n 1 ) ,
k = n ( 2 σ + n 2 ) ,
k 0 = ( n + σ 1 / 2 ) .
Based on these results, the momentum and natural variable u take the following form [63]:
p ^ = i α 2 [ B ^ + B ^ ] ; u ^ = 1 2 [ A ^ + + A ^ ] ;
where the new operators are defined by the following actions over the basis.
A ^ + Φ n σ ( u ) = ( n + 1 ) ( 2 σ + n 1 ) ( n + σ 1 / 2 ) ( n + σ + 1 / 2 ) Φ n + 1 σ ( u ) ,
A ^ Φ n σ ( u ) = n ( 2 σ + n 2 ) ( n + σ 1 / 2 ) ( n + σ 3 / 2 ) Φ n 1 σ ( u ) ,
B ^ + Φ n σ ( u ) = ( σ + n ) A ^ + Φ n σ ( u ) ,
B ^ Φ n σ ( u ) = ( σ + n 1 ) A ^ Φ n σ ( u ) .
Using these expressions, we can obtain the matrix elements of the coordinate and momentum, which induces diagonalization in the following subspace:
L N C P T = { | Φ n σ , n = 0 , 1 , , N 1 } ,
and it induces the following eigenvectors:
| u i = n = 0 N 1 Φ n σ | u i | Φ n σ ,
| p i = n = 0 N 1 Φ n σ | p i | Φ n σ ,
which provide momentum and coordinate representations characterized by the following.
u | u i = u i | u i , p | p i = p i | p i .
For the coordinate, we have relation u = tanh ( α q ) and, consequently, the following:
u j | V ( q ) | u i = V [ q i ( u i ) ] δ i j , p j | G ( p ) | p i = G ( p i ) δ i j ,
with the following being the case.
q i = 1 α arctanh ( u i ) .
Accordingly, the matrix representation of a Hamiltonian associated with a general potential takes the same form (13), albeit with new definitions T = | | Φ n σ | u i | | and W = | | Φ n σ | p i | | , with | | Λ ( p ) | | = ( p i 2 / 2 μ ) δ i j and | | Λ ( q ) | | = V [ q i ( u i ) ] δ i j in the momentum and coordinate representations, respectively.

Appendix A.2. SU(2)-UGA-DVR Method

The 1D unitary group approach has already been established by Lemus [43,44]. Here, we present salient features of the method. In this approach, a scalar boson s ( s ) is added to the physical space of a 1D harmonic oscillator with bosonic operators t ( t ) , called a ( a ) in Section 2. The bilinear products satisfy commutations relations associated with the generators of the U ( 2 ) group with the constraint that N ^ = t t + s s n ^ + n ^ s is preserved. A convenient form to express its generators is the following:
G U ( 2 ) = { N ^ , J ^ x , J ^ y , J ^ z } ,
with the following being the case.
J ^ x = 1 2 ( t s + s t ) ; J ^ y = i 2 ( t s s t ) ; J ^ z = 1 2 ( t t s s ) .
From this perspective, any dynamical variable of a 1D system may be expanded in terms of the generators (A20). The S U ( 2 ) group presents three group chains:
S U ( 2 ) U ( 1 ) ,
S U ( 2 ) S O ( 2 ) ,
S U ( 2 ) S O ¯ ( 2 ) ,
and each of them is characterized by operators J ^ z , J ^ x and J ^ y , respectively, with the following eigenkets:
J ^ z | [ N ] μ = μ | [ N ] μ = μ | [ N ] n ,
J ^ x | [ N ] ζ = ζ | [ N ] ζ ,
J ^ y | [ N ] ζ ¯ = ζ ¯ | [ N ] ζ ¯ ,
and the following eigenvalues:
μ , ζ , ζ ¯ = N 2 , N 2 + 1 , . N 2 ,
with the following relation between μ and the physical quantum number n:
n = j + μ ; N n = j μ ,
and j standing for the angular momentum associated with the S U ( 2 ) group. Since coordinate q and momentum p have the following realization:
q Q ^ = 1 2 ω μ 2 J ^ x N , p P ^ = 1 2 ω μ 2 J ^ y N ,
the bases (A26) and (A27) correspond to coordinate and momentum representations, respectively. We stress that the spectra of the coordinate and momentum are discrete, establishing a DVR approach.
The algebraic representation of any 1D system of the following:
H ^ a l g S U ( 2 ) = 1 2 μ P ^ 2 + V ( Q ^ )
is obtained through mapping q Q ^ ; p P ^ . Here, it is necessary to add and subtract a quadratic term in such a manner that the harmonic oscillator Hamiltonian is identified. We thus obtain the following:
H a l g S U ( 2 ) = ω n ^ + 1 2 n ^ 2 N + V ( Q ^ ) ,
where the following is the case.
V ( Q ^ ) = ω 2 μ 2 Q ^ 2 + V ( Q ^ ) .
Using transformation brackets of the following:
T = | | [ N ] ζ | [ N ] n | | ,
associated with the coordinates, the Hamiltonian in the energy representation takes the following form [43,44].
H = Λ ( E ) + T Λ ( Q ) T .
A remarkable fact of this approach is that the kets | [ N ] n are identified with the 1D harmonic oscillator, and consequently, it is possible to project kets | [ N ] ζ and | [ N ] ζ ¯ to position representations. It is important to take into account that, in the UGA, an accidental degeneracy appears, which must be removed [43,44]. This goal is achieved by introducing parameter ϵ in the matrix representation of the deformed oscillator contribution:
| | Λ ( E ) | | = ω n + 1 2 ϵ n 2 N δ n , n ,
where ϵ = 1 for n < N / 2 and ϵ = 0 for n N / 2 .

Appendix A.3. Morse DVR Method (M-DVR)

As in the case of the PT potential, the bound states of the Morse potential do not form a complete set of states in the Hilbert space [65,66]. The DVR approach needs a complete basis, which is provided by the tridiagonal Morse basis (TMB) given by the following [67,68,69,70]:
Φ n σ ( y ) = A n σ L n 2 σ 1 ( y ) y σ e y / 2 ,
with the following normalization constant: A n σ = ( β n ! ) / Γ ( 2 σ + 1 ) , where σ is a parameter that is chosen to be σ = 1 since, in this case, the bound states and the continuum part of the spectrum are decoupled [69]. Natural variable y is connected with the physical coordinate q through the following:
y = ( 2 j + 1 ) e β q ,
where j is related to the potential depth and β stands for the range of the Morse potential.
It is possible to construct ladder operators { K ^ ± , K ^ 0 } satisfying the s u ( 1 , 1 ) commutation relations with the following effect over the following basis (A37).
K ^ Φ n σ ( y ) = k ( n ) Φ n 1 σ ( y ) with k ( n ) = n ( 2 σ + n 1 ) ,
K ^ + Φ n σ ( y ) = k + ( n ) Φ n + 1 σ ( y ) with k + ( n ) = ( n + 1 ) ( 2 σ + n ) ,
K ^ 0 Φ n σ ( y ) = k 0 Φ n σ ( y ) with k 0 = σ + n .
On the other hand, natural variable y and momentum p = i d d q are related to the operators (A39)–(A41) in the following form:
p = i β 2 ( K ^ + K ^ ) ,
y = 2 K ^ 0 ( K ^ + + K ^ ) ,
with the following matrix elements.
Φ n σ | y | Φ n σ = 2 ( σ + n ) δ n , n ( k + ( n ) δ n , n + 1 + k ( n ) δ n , n 1 ) ,
Φ n σ | p | Φ n σ = i β 2 ( k + ( n ) δ n , n + 1 k ( n ) δ n , n 1 ) .
The diagonalization of these matrix representations in the following subspace:
L N T M B = { | Φ n σ , n = 0 , 1 , , N 1 } ,
produces the following eigenvectors:
| y i = n = 0 N 1 Φ n σ | y i | Φ n σ ,
| p i = n = 0 N 1 Φ n σ | p i | Φ n σ ,
where matrices T = | | Φ n σ | y i | | and W = | | Φ n σ | p i | | correspond to the coefficients that define basis transformations (A47) and (A48), which provide the coordinate and momentum representations characterized by the following.
y | y i = y i | y i , p | p i = p i | p i .
For the coordinate, we have relation (A38) and, consequently, the following:
y j | V ( q ) | y i = V [ q i ( y i ) ] δ i j , p j | G ( p ) | p i = G ( p i ) δ i j ,
with the following being the case.
q i = 1 β ln y i ( 2 j + 1 ) .
Equations (A50) establish the discrete variable representation and relation (A51) allows the description of systems beyond the Morse potential. Following the same procedure, the matrix representation of a Hamiltonian takes the same form (13) with the appropriate identification of the transformation coefficients and connection (A51).

Appendix B. Symmetry Group and Irreducible Representations

In this section, we present the symmetry group of particle in a 2D square potential with infinite walls together with the construction of its irreducible representations. Although this problem has already been presented in heuristic manner [9], it has not been formally presented in the framework of induced representations as it was for the cubic 3D square well potential [11].
Let us first start by presenting the solutions of one particle inside the potential depicted in Figure A1, where the origin is located at the left corner of the square. In this reference framework, the solutions take the following simple form:
ψ n 1 n 2 ( x , y ) = ψ n 1 ( x ) ψ n 2 ( y ) ,
where the following is the case.
ψ n 1 ( x ) = 2 L sin n 1 π x L ,
ψ n 2 ( y ) = 2 L sin n 2 π y L .
Figure A1. Square well potential and symmetry elements associated with geometrical symmetry C 4 v . Two reference systems are shown: (a) the origin is located at the left corner with solutions given by (A53) and (A54); and (b) the origin is located at the center of the square.
Figure A1. Square well potential and symmetry elements associated with geometrical symmetry C 4 v . Two reference systems are shown: (a) the origin is located at the left corner with solutions given by (A53) and (A54); and (b) the origin is located at the center of the square.
Symmetry 14 00445 g0a1
The symmetry projection of functions (A52) allows the assignment of irreducible representation to the energy spectrum depicted in Figure 10 [9]. In this system, we are able to identify two subspaces. The one-dimensional space of the following:
L 1 = { | ψ n n }
with the following reduction:
| ψ n n B 2 , A = n even
| ψ n n A 1 , A = n odd ,
and the following two dimensional subspaces:
L 2 = { | ψ n m } ; n m ,
with reductions provided in Table A1 is provided. The first 13 energy levels obtained with Equation (28) in dimensionless units E ¯ = E / ( 2 π 2 / 2 m L 2 ) are displayed in Figure 10. It is, thus, clear that double degeneracy A 1 B 1 and A 2 B 2 represent accidental degeneracy.
Table A1. Irreducible representations (irreps) of the geometrical group C 4 v contained in the subspaces (A58). The reduction depends of the parity of n and m with p and q integers.
Table A1. Irreducible representations (irreps) of the geometrical group C 4 v contained in the subspaces (A58). The reduction depends of the parity of n and m with p and q integers.
Irrepnm
E 2 p 2 q + 1
E 2 p + 1 2 q
A 1 B 1 2 p + 1 2 q + 1
A 2 B 2 2 p 2 q
However, since this degeneracy is systematic over the entire spectrum, we have the presence of systematic accidental degeneracy, which implies that symmetry group C 4 v is a subgroup of the true symmetry group. In order to obtain the true symmetry group, it is necessary to identify the operator connecting the accidental degenerate states. In order to achieve this goal, we start establishing that such operators F ^ ( ρ ) , carrying the ρ -th irreducible representation (irrep), must satisfy the following.
ψ Γ | F ^ ( ρ ) | ψ Γ 0 ; for ( Γ = A 1 , Γ = B 1 ) and ( Γ = A 2 , Γ = B 2 ) .
In other words, the operator must span the irrep ρ satisfying the following.
Γ ρ Γ .
In accordance with the character in Table A2, such operators should span the irreducible representation ρ = B 1 . Notice that degenerate states | ψ i E are left unaltered by the B 1 tensor.
Table A2. Character table of group C 4 v .
Table A2. Character table of group C 4 v .
C 4 v E 2 C 4 C 2 2 σ v 2 σ d
A 1 1 1 1 1 1 z x 2 + y 2 ; z 2
A 2 1 1 1 1 1 R z
B 1 1 1 1 1 1 x 2 y 2
B 2 1 1 1 1 1 x y
E2 0 2 0 0 ( x , y ) ; ( R x , R y ) ( x z , y z )
From the character in Table A2, we see that the Cartesian harmonic x 2 y 2 spans irrep B 1 , where the primes mean that the function is referred to the origin located in panel (b) of Figure A1. However, the same linear combination in terms of the square momenta transforms according to B 1 and, consequently, we have the following operator.
F ^ ( B 1 ) = 1 4 5 π 2 x 2 2 y 2 .
This operator commutes with the following Hamiltonian:
[ H ^ , F ^ ( B 1 ) ] = 0 ,
and conserves boundary conditions. Consequently, it is suitable to be added to the symmetry elements of the C 4 v group. Operator (A61) generates continuous group T with elements obtained by exponentiation in the following manner.
U ^ ( α ) = e i α F ^ ( B 1 ) .
In order to elucidate the structure of the new group it is convenient to observe whether subgroup T is invariant. Taking into account that elements R C 4 v are isomorphic to operators O ^ R that act over the space of functions, the conjugation of U ^ ( α ) with respect to the elements of C 4 v yields the following:
O ^ R U ^ ( α ) O ^ R 1 = U ^ ( α ) T ,
with the following being the case:
α = α χ B 1 ( R ) ,
where χ B 1 ( R ) denotes the character of the irrep B 1 in group C 4 v . This result means that group T is invariant, a fact that allows the new symmetry group G to be expressed as a semidirect product:
G = T C 4 v ,
and in terms of an expansion of left cosets.
G = λ = 1 | λ | s λ T ; | λ | = | G | | T | ; s λ C 4 v .
Hence, any element g G can be written in the following form:
g ^ = O ^ R U ^ ( α ) ; R C 4 v ,
where we have simplified the following notation g ^ = O ^ g .
We now proceed to construct irreducible representations of group G . To accomplish this task, we shall proceed through the induction method by constructing the irreps of group G from the irreps of subgroup T . We start with the construction of the irreps of the invariant subgroup, T .
First, let us consider the two dimensional representation spaces.
L 2 = { | ψ n 1 n 2 , | ψ n 2 n 1 } ; n 1 n 2 .
The action of the elements of T over this space is given by the following:
U ^ ( α ) | ψ n 1 n 2 = D ( k n ) ( α ) | ψ n 1 n 2
where the following is the case:
D ( k n ) ( α ) = e i α k n
with the following.
k n = ζ 0 ( n 1 2 n 2 2 ) ; ζ 0 = π 2 L 2 1 4 5 π .
We, thus, have states | ψ n 1 n 2 spanning representation k n . It is, thus, convenient to label the states according to this representation in the following form | ψ n 1 n 2 k n in such a manner that the following is the case.
U ^ ( α ) | ψ n 1 n 2 k n = e i α k n | ψ n 1 n 2 k n .
For the sake of convenience, we introduce the following definitions.
k 1 = k n 1 n 2 = ζ 0 ( n 1 2 n 2 2 ) ; k 2 = k n 2 n 1 = ζ 0 ( n 2 2 n 1 2 ) = k 1 .
Hence, the two-dimensional space (A69) is a representation space spanning the diagonal irreps of the subgroup T .
D ( k ) ( α ) = e i α k 1 0 0 e i α k 2
On the other hand, it is clear that for n 1 = n 2 = n , we have a one-dimensional space L 1 = { | ψ n n k 0 } spanning representation D ( k 0 ) ( α ) = 1 with k 0 = 0 .
U ^ ( α ) | ψ n n k 0 = D ( k 0 ) ( α ) | ψ n n k 0 = | ψ n n k 0 .
Once we count on the irreps of subgroup T , we are ready to proceed to construct the irreps of the group G . However, to accomplish this task, it is necessary to recall the concept of conjugate representation.
Let H G , with H invariant. Given D ( μ ) ( h ) , the μ -th irrep of h H , the matrix D ( μ ) ( g h g 1 ) μ D ( g ) ( h ) called conjugate representation is also a representation of H. In our case, H T is an invariant subgroup with representations D ( k n ) ( α ) = e i α k n . We, thus, have to identify the conjugate representation. In order to achieve this goal, consider the representation k 1 defined by the following.
U ^ ( α ) | ψ n 1 n 2 k 1 = e i α k 1 | ψ n 1 n 2 k 1 .
From (A64), we have the following.
O R U ^ ( α ) O R 1 = U ( α ) ; α = α χ B 1 ( R ) .
Hence, we have the following for the conjugate representation.
k 1 D ( R ) ( α ) = D ( k 1 ) ( α ) .
This result can be reinterpreted in the following form:
k 1 D ( R ) ( α ) = e i α R ^ F ^ ( B 1 ) R ^ 1 = e i α χ B 1 ( R ) F ^ ( B 1 ) = e i α k = D ( k ) ( α ) ,
where the following is the case:
k = k 1 χ B 1 ( R ) ,
with the proviso that the operators involved act over state | ψ n 1 n 2 k 1 . Here, we have used simplified notation, O R = R ^ . In addition, it is worth noticing that the following is the case:
ψ n 1 n 2 k 1 | R ^ F ^ ( B 1 ) R ^ 1 | ψ n 1 n 2 k 1 = R ^ 1 ψ n 1 n 2 k 1 | F ^ ( B 1 ) | R ^ 1 ψ n 1 n 2 k 1 ,
which means that the state spanning the conjugate representation k is obtained by the following.
R ^ 1 | ψ n 1 n 2 k 1 = | ψ n 1 n 2 k .
We notice that two irreps are identified, namely k 1 and k 2 . Starting with k 1 , we obtain conjugate representations k 1 D ( R ) ( α ) = D ( k 1 ) ( R U ( α ) R 1 ) , R C 4 v , although some of them are equivalent. The set of nonequivalent representations, denoted by k 1 S , is called the star of k 1 . In our case, from (A81), only two nonequivalent representations are obtained: k 1 itself and k = k 2 = k 1 , the latter generated with R C 4 v in (A81) satisfying χ B 1 ( R ) = 1 . The star is then given by the following:
k S = k 1 S = { k 1 , k 2 } ,
where each element is called the prong of the star. Starting with k 2 , the same star is obtained. This explains that the star will be referred as k S without specifying the particular k i from which it is generated. It should be clear that | C 4 v | = 8 conjugations are carried out to obtain (A84) of order | k S | = 2 . Consequently, there should be a set of elements that transform k 1 to either an equivalent or identical representation. Indeed this set of elements corresponds to the elements of T as well as the elements of subgroup C 2 v . This set of transformation form a group, the little group of k 1 , denoted by K ( k ) and is given by the following.
K ( k ) = T C 2 v .
Each prong in star (A84) has associated its own little group, but they are isomorphic, a fact that explains the notation of K ( k ) . The little group of k is infinite. To deal with a finite group, we consider the factor group called little co-group K ( k ) :
K ( k ) = K ( k ) T C 2 v .
in which its elements are provided by the left cosets in the following expansion.
K ( k ) = λ | λ | s λ T ; s λ C 2 v .
On the other hand, the factor group is isomorphic to group C 4 v :
G T C 4 v ,
with the following expansion in left cosets:
C 4 v = λ | λ | p λ C 2 v ; | λ | = | C 4 v | | C 2 v | ,
a fact that allows confrontation with finite groups. This is the basic background in the construction of the irreps of G . Expansion (A89) implies that every element g C 4 v can be expressed in terms of a product of the following form:
g = p λ h ; h C 2 v ,
and in the following explicit form:
C 4 v = C 2 v + σ d a C 2 v ,
with identification p 1 = E , p 2 = σ d a .
Let us now consider the ket | ψ n 1 n 2 k 1 . All the elements of C 2 v leave it invariant. On the other hand the action of the coset representatives in the expansion (A89) yields
O E | ψ n 1 n 2 k 1 =       | ψ n 1 n 2 k 1 ,
O σ d a | ψ n 1 n 2 k 1 =       | ψ n 2 n 1 k 2 .
Hence, the set of functions { | ψ n 1 n 2 k 1 , | ψ n 2 n 1 k 2 } form a representation space of G . However, for the construction of irreps, it is convenient to consider the basis in the following form: { p ^ λ | ψ n 1 n 2 k 1 ; λ = 1 , 2 } . The procedure to obtain the irreps of element g C 4 v consists in applying operator g ^ over the representation space. In this manner, we obtain the following:
g ^ p ^ λ | ψ n 1 n 2 k 1 = p ^ μ h ^ | ψ n 1 n 2 k 1 ; h C 2 v ,
where identity follows from expansion (A89). Here, h is called the sub-element of g in p λ denoted by h λ ( g ) . All sub-elements involved are listed in Table A3. We should remark that the action over states | ψ n 1 n 2 k 1 are well defined since they carry irreps of C 2 v .
Table A3. Sub-elements h λ ( g ) defined in (A94) through expansion (A89).
Table A3. Sub-elements h λ ( g ) defined in (A94) through expansion (A89).
g p λ p λ h h λ ( g )
EEE EE
C 4 E σ d a σ v a σ v a
C 4 3 E σ d a σ v b σ v b
σ v a E E σ v a σ v a
σ v b E E σ v b σ v b
σ d a E σ d a E E
σ d b E σ d a C 2 C 2
E σ d a σ d a E E
C 4 σ d a E σ v b σ v b
C 4 3 σ d a E σ v a σ v a
σ v a σ d a σ d a σ v b σ v b
σ v b σ d a σ d a σ v a σ v a
σ d a σ d a E E E
σ d b σ d a E C 2 C 2
We are now prepared to carry out the induction. However, before we accomplish this goal, it is convenient to present the general procedure. Let us consider the general situation of generating representations of a group G from the irreps D ( μ ) = μ H of H G [10]. First, we expand group G in cosets of H.
G = λ | λ | s λ H ; | λ | = | G | | H | .
The set of cosets { s λ H } span a representation of G:
g ^ ( s λ H ) = λ | λ | Δ H ( b ) ( g ) λ λ ( s λ H ) ,
where Δ ( b ) ( g ) is called the basal representation given by the following.
Δ H ( b ) ( g ) λ λ = δ λ τ ; g s λ = s τ h λ ( g ) .
Let us now consider a set of kets | ϕ i ( μ ) spanning the μ -th irrep μ H = D ( μ ) ( H ) of subgroup H.
O ^ R | ϕ i ( μ ) = j n μ D j i ( μ ) ( R ) | ϕ j ( μ ) .
Hence, set { s ^ λ | ϕ i ( μ ) } is a representation the space of group G:
g ^ ( s ^ λ | ϕ i ( μ ) ) = λ j Δ ( μ H G ) ( g ) λ j ; λ i ( s ^ λ | ϕ j ( μ ) )
where Δ ( μ H G ) ( g ) is called the induced representation with the following elements.
Δ ( μ H G ) ( g ) λ j ; λ i = Δ H ( b ) ( g ) λ λ D j i ( μ ) ( h λ ( g ) ) .
In general, the induced representation is reducible. However, when induction is carried out through the little co-groups associated with the irreps of the subgroup, the induced representations are complete and irreducible. We may, thus, sketch the general procedure as follows. In the frame work of our group G = T C 4 v , we begin these states, | ψ n 1 n 2 k 1 , by carrying irreps of the invariant subgroup T . Then, the little group K ( k ) is identified. This group in infinite; consequently, to deal with a finite group, it is considered the little co-group isomorphic to a subgroup of point group C 4 v . The states | ψ n 1 n 2 k 1 span irreps Γ of the little co-group, a fact that permits a labeling scheme, | ψ n 1 n 2 k 1 ; Γ , γ . Finally, induction is carried out by obtaining all irreducible representations of G .
k 1 Γ , γ G = Γ , γ K ( k 1 T ) G .
Here, there are two labels, the prong corresponding to the irrep k 1 T , and the irrep and component { Γ , γ } associated with the little co-group.
We now proceed to generate irreps (A101). In our previous analysis, we were dealing with two-dimensional space (A69) leading to the star (A84), which may be identified with vectors { k 1 , k 2 = k 1 } , as shown in Figure A2. In this figure, we denote k = k 1 in accordance with the notation for the little co-group. Each prong of the two dimensional star is called a general vector with little co-group C 2 v . On the other hand, the irrep k 0 T spanned by the kets | ψ n n k 0 leads to the star S 0 = { k = 0 } . Because this irrep k 0 = 0 is invariant under full point group C 4 v , it is called a special point.
Figure A2. Irreducible representations of the group T labeled by the corresponding little co-groups. Two little co-groups have been identified: K ( k ) corresponding to a general vector and K ( k 0 ) corresponding to special point k 0 = 0 .
Figure A2. Irreducible representations of the group T labeled by the corresponding little co-groups. Two little co-groups have been identified: K ( k ) corresponding to a general vector and K ( k 0 ) corresponding to special point k 0 = 0 .
Symmetry 14 00445 g0a2
We start considering general points associated with the star (A84). Since in this case the co-group is C 2 v , states { | ψ n 1 n 2 k 1 , | ψ n 2 n 1 k 2 } by itself carry irreps of C 2 v . Indeed, considering the generators { C 2 , σ v a } of the little co-group C 2 v , we obtain diagonal representations.
Δ ( C 2 ) = ( 1 ) n 1 + n 2 0 0 ( 1 ) n 1 + n 2 ; Δ ( σ v a ) = ( 1 ) n 2 + 1 0 0 ( 1 ) n 1 + 1 .
This representation is of course reducible and contains the two irreps given by the following:
μ : D ( μ ) ( C 2 ) = ( 1 ) n 1 + n 2 ; D ( μ ) ( σ v a ) = ( 1 ) n 2 + 1 ;
ν : D ( ν ) ( C 2 ) = ( 1 ) n 1 + n 2 ; D ( ν ) ( σ v a ) = ( 1 ) n 1 + 1 ,
which have to be identified with the irreps of the C 2 v group in accordance to the character in Table A4. Since the little co-groups for k 1 and k 2 are isomorphic, it is enough to consider the reduction associated with k 1 , which is displayed in Table A5.
Table A4. Character table for group C 2 v .
Table A4. Character table for group C 2 v .
C 2 v E C 2 σ v ( xz ) σ v ( yz )
A 1 1 1 1 1 z x 2 ; y 2 ; z 2
A 2 1 1 1 1 R z x y
B 1 1 1 1 1 x ; R y x z
B 2 1 1 1 1 y ; R x y z
Table A5. Identification of the irreps of C 2 v in accordance with parity.
Table A5. Identification of the irreps of C 2 v in accordance with parity.
n 1 = 2 m + 1 n 1 = 2 m n 1 = 2 m n 1 = 2 m + 1
n 2 = 2 n + 1 n 2 = 2 n n 2 = 2 n + 1 n 2 = 2 n
D ( μ ) ( C 2 ) 1 1 1 1
D ( μ ) ( σ v a ) 1 1 1 1
k 1 : μ A 1 A 2 B 1 B 2
We proceed to obtain the representations of G that are completely specified by the matrix representation of generators { U ( α ) , C 4 , σ d a } . Let us start considering element g = C 4 . In accordance with the coset expansion (A89), the representation space is { E | ψ k 1 , σ d a | ψ k 1 } . From Table A3 for the sub-elements, we obtain the following:
C 4 [ E | ψ k 1 , μ ] = σ d a σ v a | ψ k 1 , μ = χ ( μ ) ( σ v a ) [ σ d a | ψ k 1 , μ ] ,
C 4 [ σ d a | ψ k 1 , μ ] = E σ v b | ψ k 1 , μ = χ ( μ ) ( σ v b ) [ E | ψ k 1 , μ ] ,
from which we obtain the matrix representation for element C 4 .
D ( k μ G ) ( C 4 ) = 0 χ ( μ ) ( σ v b ) χ ( μ ) ( σ v a ) 0 .
In similar manner for element σ d a , we have the following:
σ d a [ E | ψ k 1 , μ ] = χ ( μ ) ( E ) [ σ d a | ψ k 1 , μ ] ,
σ d a [ σ d a | ψ k 1 , μ ] = χ ( μ ) ( E ) [ E | ψ k 1 , μ ] ,
and, consequently, we have the following.
D ( k μ G ) ( σ d a ) = 0 χ ( μ ) ( E ) χ ( μ ) ( E ) 0 .
Finally, we have to obtain the matrix representation of element U ( α ) . The action of an element of T over an element of the basis is given by the following:
U ^ ( α ) R ^ | ψ n 1 n 2 k 1 , μ = R ^ U ^ ( α ) | ψ n 1 n 2 k 1 , μ ; R = E , σ d a .
where we have taken into account the invariance of T . Nevertheless, from (A83), the new element U ( α ) is associated with the inverse of R.
U ^ ( α ) | ψ n 1 n 2 k 1 , μ = k 1 D ( R 1 ) ( α ) | ψ n 1 n 2 k 1 , μ = e i α k | ψ n 1 n 2 k 1 , μ .
However, in this case, R 1 = R and, consequently, the following is the case.
D ( k μ G ) ( α ) = e i α k 0 0 e i α k .
We now consider the special point k 0 = 0 . In this case, the little co-group is C 4 v and the basal representation is generated by element E C 4 v with basis | ψ n n k 0 , μ . The action of the generators is as follows.
U ^ ( α ) = 1 ,
C ^ 4 ( E | ψ n n k 0 , μ ) = ( 1 ) n + 1 ( E | ψ n n k 0 , μ ) ,
σ ^ d a ( E | ψ n n k 0 , μ ) = ( E | ψ n n k 0 , μ ) .
From the character in Table A2, we have identifications for the following representations:
n even D ( k 0 B 2 G ) ( C 4 ) = 1 ; D ( k 0 B 2 G ) ( σ d a ) = 1 ;
n odd D ( k 0 A 1 G ) ( C 4 ) = 1 ; D ( k 0 A 1 G ) ( σ d a ) = 1 ,
and in compact form in terms of characters.
D ( k 0 μ G ) ( α ) = 1 ;
D ( k 0 μ G ) ( C 4 ) = χ ( μ ) ( C 4 ) ; ; D ( k 0 μ G ) ( σ d a ) = χ ( μ ) ( σ d a ) .
Finally, we present subduction k μ G C 4 v in order to make clear that the accidental degeneracy previously identified in C 4 v renders a natural degeneracy in group G . In Table A6, we present subduction k μ G C 4 v where we can observe that the accidental degeneracy in accordance with geometrical group C 4 v is indeed a natural degeneracy in the new symmetry group, G .
Table A6. Subduction k μ G C 4 v .
Table A6. Subduction k μ G C 4 v .
C 4 v E C 4 , C 4 3 C 4 2 σ v a , σ v b σ d a , σ d b k μ G C 4 v
χ ( k o A 1 G ) 1 1 1 1 1 A 1
χ ( k o B 2 G ) 1 1 1 1 1 B 2
χ ( k A 1 G ) 2 0 2 2 0 A 1 B 1
χ ( k A 2 G ) 2 0 2 2 0 A 2 B 2
χ ( k B 1 G ) 2 0 2 0 0E
χ ( k B 2 G ) 2 0 2 0 0E
Remark A1.
In the construction of irreps of space groups, an induction process is followed from the irreps of the translational group to the space group. The procedure we have presented has a close analogy with space groups. Our group T corresponds to the translational group while irreps k corresponds to the momentum space involved in the Bloch functions. This is the reason we have used the language of solid state physics to developed the representation theory of the new group, T C 4 v .

References

  1. Cohen-Tannoudji, C.; Diu, B.; Laloë, F. Quantum Mechanics; John Wiley and Sons: New York, NY, USA, 1977; Volume I. [Google Scholar]
  2. Levin, F.S. An Introduction to Quantum Theory; Cambridge University Press: Cambridge, UK, 2002. [Google Scholar]
  3. Greiner, W. Quantum Mechanics. An Introduction; Springer: Frankfurt am Main, Germany, 1989. [Google Scholar]
  4. Harris, D.; Bertolucci, M.D. Symmetry and Spectroscopy. An Introduction to Vibrational and Electronic Spectroscopy; Dover Publications, Inc.: New York, NY, USA, 1989. [Google Scholar]
  5. MacIntosh, H. Symmetry and Degeneracy. Group Theory and Its Application; Loebel, E.M., Ed.; Academic Press: New York, NY, USA, 1971; Volume 2, pp. 75–144. [Google Scholar]
  6. Wybourne, B. Classical Groups for Physicists; Wyley: New York, NY, USA, 1974. [Google Scholar]
  7. Fock, V. Zur Theorie des Wasserstoffatoms. Z. Phys. 1935, 98, 145. [Google Scholar] [CrossRef]
  8. Bargmann, V. Zur Theorie des Wasserstoffatoms. Z. Phys. 1936, 99, 576. [Google Scholar] [CrossRef]
  9. Leybraz, F.; Frank, A.; Lemus, R.; Andrés, M.V. Accidental degeneracy in a simple quantum system: A new symmetry group for a particle in an impenetrable square-well potential. Am. J. Phys. 1997, 65, 1087. [Google Scholar] [CrossRef]
  10. Altmann, S.L. Induced Representations in Crystals and Molecules. Point, Space and Non-Rigid Molecular Groups; Academic: New York, NY, USA, 1977. [Google Scholar]
  11. Hernández-Castillo, A.O.; Lemus, R. Symmetry group of an impenetrable cubic well potential. J. Phys. A Math. Theor. 2013, 46, 465201. [Google Scholar] [CrossRef]
  12. Lemus, R.; Frank, A.; Andrés, M.V.; Leybraz, F. Accidental degeneracy and hidden symmetry: Rectangular wells with commensurate sides. Am. J. Phys. 1998, 66, 629. [Google Scholar] [CrossRef]
  13. Harris, D.O.; Engerholm, G.G.; Gwinn, W.D. Calculation of Matrix Elements for One-Dimensional Quantum-Mechanical Problems and the Application to Anharmonic Oscillators. J. Chem. Phys. 1965, 43, 1515. [Google Scholar] [CrossRef] [Green Version]
  14. Dickinson, A.S.; Certain, P.R. Calculation of Matrix Elements for One-Dimensional Quantum-Mechanical Problems. J. Chem. Phys. 1968, 49, 4209. [Google Scholar] [CrossRef]
  15. Lill, J.V.; Parker, G.A.; Light, J.C. Discrete variable representations and sudden models in quantum scattering theory. Chem. Phys. Lett. 1982, 89, 483. [Google Scholar] [CrossRef]
  16. Light, J.C.; Hamilton, I.P.; Lill, J.V. Generalized discrete variable approximation in quantum mechanics. J. Chem. Phys. 1985, 82, 1400. [Google Scholar] [CrossRef]
  17. Baĉić, Z.; Light, J.C. Highly excited vibrational levels of “floppy” triatomic molecules: A discrete variable representation—Distributed Gaussian basis approach. J. Chem. Phys. 1986, 85, 4594. [Google Scholar] [CrossRef]
  18. Lill, J.V.; Parker, G.A.; Light, J.C. The discrete variable–finite basis approach to quantum scattering. J. Chem. Phys. 1986, 85, 900. [Google Scholar] [CrossRef]
  19. Light, J.C.; Baĉić, Z. Adiabatic approximation and nonadiabatic corrections in the discrete variable representation: Highly excited vibrational states of triatomic molecules. J. Chem. Phys. 1987, 87, 4008. [Google Scholar] [CrossRef]
  20. Baĉić, Z.; Light, J.C. Accurate localized and delocalized vibrational states of HCN/HNC. J. Chem. Phys. 1987, 86, 3065. [Google Scholar] [CrossRef]
  21. Choi, S.E.; Light, J.C. Use of the discrete variable representation in the quantum dynamics by a wave packet propagation: Predissociation of NaI(1 Σ + 0 )→ NaI(0+)→ Na(2S) + I(2P). J. Chem. Phys. 1989, 90, 2593. [Google Scholar] [CrossRef]
  22. Tennison, J. Highly excited rovibrational states using a discrete variable representation: The H + 3 molecular ion. J. Chem. Phys. 1989, 91, 3815. [Google Scholar] [CrossRef]
  23. Light, J.C.; Carrington, T., Jr. Discrete-Variable Representations and their Utilization In Advances in Chemical Physics; Prigogine, I., Rice, S.A., Eds.; John Wiley & Sons: Hoboken, NJ, USA, 2000; Volume 114, p. 263. [Google Scholar]
  24. Shizgal, B.D.; Blackmore, R. A discrete ordinate method of solution of linear boundary value and eigenvalue problems. J. Comput. Phys. 1984, 55, 313. [Google Scholar] [CrossRef]
  25. Shizgal, B.D.; Chen, H. The quadrature discretization method (QDM) in the solution of the Schrödinger equation with nonclassical basis functions. J. Chem. Phys. 1996, 104, 4137. [Google Scholar] [CrossRef]
  26. Shizgal, B.D. Spectral Methods in Chemistry and Physics. Applications to Kinetic Theory and Quantum Mechanics; Springer: New York, NY, USA, 2015. [Google Scholar]
  27. Shizgal, B.D. Pseudospectral method of solution of the Schrödinger equation with non classical polynomials; the Morse and Pöschl–Teller (SUSY) potentials. Comp. Theor. Chem. 2016, 1084, 51. [Google Scholar] [CrossRef]
  28. Shizgal, B.D. A comparison of pseudospectral methods for the solution of the Schrödinger equation; the Lennard-Jones (n, 6) potential. Comp. Theor. Chem. 2017, 1114, 25. [Google Scholar] [CrossRef]
  29. Bao, J.; Shizgal, B.D. Pseudospectral method of solution of the Schrödinger equation for the Kratzer and pseudoharmonic potentials with nonclassical polynomials and applications to realistic diatom potentials. Comp. Theor. Chem. 2019, 1149, 49. [Google Scholar] [CrossRef]
  30. Baye, D.; Heenen, P.H. Generalised meshes for quantum mechanical problems. J. Phys. A Math. Gen. 1986, 19, 2041. [Google Scholar] [CrossRef]
  31. Vincke, M.; Malegat, L.; Baye, D. Regularization of singularities in Lagrange-mesh calculations. J. Phys. B At. Mol. Opt. Phys. 1993, 26, 811. [Google Scholar] [CrossRef]
  32. Baye, D. Lagrange-mesh method for quantum-mechanical problems. Phys. Stat. Sol. 2006, 243, 1095. [Google Scholar] [CrossRef]
  33. Baye, D. The Lagrange-mesh method. Phys. Rep. 2015, 565, 1. [Google Scholar] [CrossRef] [Green Version]
  34. Carvajal, M.; Arias, J.M.; Gómez-Camacho, J. Configuration localized Morse wave functions: Application to vibrational transitions in anharmonic diatomic molecules. Phys. Rev. A 1999, 59, 1852. [Google Scholar] [CrossRef]
  35. Carvajal, M.; Arias, J.M.; Gómez-Camacho, J. Analytic evaluation of Franck-Condon integrals for anharmonic vibrational wave functions. Phys. Rev. A 1999, 59, 3462. [Google Scholar] [CrossRef]
  36. Pérez-Bernal, F.; Arias, J.M.; Carvajal, M.; Gómez-Camacho, J. Configuration localized wave functions: General formalism and applications to vibrational spectroscopy of diatomic molecules. Phys. Rev. A 2000, 61, 042504. [Google Scholar] [CrossRef] [Green Version]
  37. Wei, H.; Carrington, T. The discrete variable representation of a triatomic Hamiltonian in bond length–bond angle coordinates. J. Chem. Phys. 1992, 97, 3029. [Google Scholar] [CrossRef]
  38. Wang, X.-G.; Carrington, T. A discrete variable representation method for studying the rovibrational quantum dynamics of molecules with more than three atoms. J. Chem. Phys. 2009, 130, 094101. [Google Scholar] [CrossRef]
  39. Bramley, M.J.; Carrington, T. A general discrete variable method to calculate vibrational energy levels of three- and four-atom molecules. J. Chem. Phys. 1993, 99, 8519. [Google Scholar] [CrossRef]
  40. Bramley, M.J.; Tromp, J.W.; Carrington, T.; Corey, G.C. Efficient calculation of highly excited vibrational energy levels of floppy molecules: The band origins of H + 3 up to 35,000 cm−1. J. Chem. Phys. 1994, 100, 6175. [Google Scholar] [CrossRef]
  41. Wang, X.-G.; Carrington, T. The Utility of Constraining Basis Function Indices When Using the Lanczos Algorithm to Calculate Vibrational Energy Levels. J. Phys. Chem. A 2001, 105, 2575. [Google Scholar] [CrossRef]
  42. Iachello, F.; Levine, R.D. Algebraic Theory of Molecules; Oxford University Press: Oxford, UK, 1995. [Google Scholar]
  43. Lemus, R. Unitary group approach for effective molecular potentials: 1D systems. Mol. Phys. 2019, 117, 167. [Google Scholar] [CrossRef]
  44. Lemus, R. A simple approach to solve the time independent Schröedinger equation for 1D systems. J. Phys. Comm. 2019, 3, 025012. [Google Scholar] [CrossRef]
  45. Estévez-Fregoso, M.M.; Arias, J.M.; Gómez-Camacho, J.; Lemus, R. An approach to establish a connection between algebraic and configuration spaces: su(ν + 1) algebraic model for vibrational excitations. Mol. Phys. 2018, 116, 2254. [Google Scholar] [CrossRef]
  46. Estévez-Fregoso, M.M.; Lemus, R. Connection between the su(3) algebraic and configuration spaces: Bending modes of linear molecules. Mol. Phys. 2018, 116, 2374. [Google Scholar] [CrossRef]
  47. Rodríguez-Arcos, M.; Lemus, R. Unitary group approach for effective potentials in 2D systems: Application to carbon suboxide C3O2. Chem. Phys. Lett. 2018, 713, 266. [Google Scholar] [CrossRef] [Green Version]
  48. Rodríguez-Arcos, M.; Lemus, R.; Arias, J.M.; Gómez-Camacho, J. Unitary group approach to describe interatomic potentials in 3D systems. Mol. Phys. 2019, 118, e1662959. [Google Scholar] [CrossRef]
  49. Rodríguez-Arcos, M.; Bermúdez-Montaña, M.; Arias, J.M.; Gómez-Camacho, J.; Orgaz, E.; Lemus, R. Algebraic discrete variable representation approaches: Application to interatomic effective potentials. Mol. Phys. 2021, 119, e1876264. [Google Scholar] [CrossRef]
  50. Bermúdez-Montaña, M.; Rodríguez-Arcos, M.; Lemus, R.; Arias, J.M.; Gómez-Camacho, J.; Orgaz, E. Algebraic DVR Approaches Applied to Describe the Stark Effect. Symmetry 2020, 12, 1719. [Google Scholar] [CrossRef]
  51. Moshinsky, M. The Harmonic Oscillator in Modern Physics: From Atoms to Quarks; Gordon and Breach: New York, NY, USA, 1969. [Google Scholar]
  52. Littlejohn, R.G.; Cargo, M. A general framework for discrete variable representation basis sets. J. Chem. Phys. 2002, 116, 8691. [Google Scholar] [CrossRef] [Green Version]
  53. Barletta, P.; Lombardi, M.; Kievsky, A. Variational DVR Calculations. Few-Body Syst. 2004, 34, 11. [Google Scholar] [CrossRef]
  54. Chen, J.-Q. Group Representation Theory for Physicists; World Scientific: New York, NY, USA, 1989. [Google Scholar]
  55. Lemus, R. A general method to obtain vibrational symmetry adapted bases in a local scheme. Mol. Phys. 2003, 101, 2511. [Google Scholar] [CrossRef]
  56. Álvarez-Bajo, O.; Lemus, R.; Carvajal, M.; Pérez-Bernal, F. Equivalent rotations associated with the permutation inversion group revisited: Symmetry projection of the rovibrational functions of methane. Mol. Phys. 2011, 109, 797. [Google Scholar] [CrossRef] [Green Version]
  57. Lemus, R. Quantum Numbers and the Eigenfunction Approach to Obtain Symmetry Adapted Functions for Discrete Symmetries. Symmetry 2012, 4, 667. [Google Scholar] [CrossRef]
  58. Bermúdez-Montaña, M.; Lemus, R.; Castaños, O. Fidelity, entropy, and Poincaré sections as tools to study the polyad breaking phenomenon. EPL 2016, 116, 13001. [Google Scholar] [CrossRef] [Green Version]
  59. Bermúdez-Montaña, M.; Lemus, R.; Castaños, O. Polyad breaking phenomenon associated with a local-to-normal mode transition and suitability to estimate force constants. Mol. Phys. 2017, 115, 3076. [Google Scholar] [CrossRef]
  60. Pöschl, G.; Teller, E. Bemerkungen zur Quantenmechanik des anharmonischen Oszillators. Z. Phys. 1933, 83, 143. [Google Scholar] [CrossRef]
  61. Rosen, N.; Morse, P.M. On the Vibrations of Polyatomic Molecules. Phys. Rev. 1932, 42, 210. [Google Scholar] [CrossRef]
  62. Landau, L.D.; Lifshitz, E.M. Quantum Mechanics, Non Relativistic Theory, 3rd ed.; Pergamon: Oxford, UK, 1997. [Google Scholar]
  63. Arias, J.M.; Gómez-Camacho, J.; Lemus, R. An su(1, 1) dynamical algebra for the Pöschl–Teller potential. J. Phys. A Math. Gen. 2004, 37, 877. [Google Scholar] [CrossRef]
  64. Gradshtein, I.S.; Rizhic, I.M. Tables of Integrals, Series and Product, 5th ed.; Academic Press: New York, NY, USA, 1994. [Google Scholar]
  65. Morse, P.M. Diatomic Molecules According to the Wave Mechanics. II. Vibrational Levels. Phys. Rev. 1929, 34, 57. [Google Scholar] [CrossRef]
  66. Matsumoto, A. Generalised matrix elements in discrete and continuum states for the Morse potential. J. Phys. B: At. Mol. Opt. Phys. 1988, 21, 2863. [Google Scholar] [CrossRef]
  67. Benedict, M.G.; Molnár, B. Algebraic construction of the coherent states of the Morse potential based on supersymmetric quantum mechanics. Phys. Rev. A 1999, 60, R1737. [Google Scholar] [CrossRef]
  68. Molnar, B.; Földi, P.; Benedict, M.G.; Bartha, F. Time evolution in the Morse potential using supersymmetry: Dissociation of the NO molecule. EPL 2003, 61, 445. [Google Scholar] [CrossRef]
  69. Lemus, R.; Arias, J.M.; Gómez-Camacho, J. An su(1, 1) dynamical algebra for the Morse potential. J. Phys. A Math. Gen. 2004, 37, 1805. [Google Scholar] [CrossRef]
  70. Bardoni, A.; Manini, N. Systematic calculation of molecular vibrational spectra through a complete Morse expansion. Int. J. Quantum Chem. 2007, 107, 782. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Symmetric square well potential. The origin was chosen in order to obtain wave functions with good parity.
Figure 1. Symmetric square well potential. The origin was chosen in order to obtain wave functions with good parity.
Symmetry 14 00445 g001
Figure 2. Comparison between the errors associated with the HO-DVR (solid circles) and PT-DVR (triangles) methods in accordance with the results provided in Table 2.
Figure 2. Comparison between the errors associated with the HO-DVR (solid circles) and PT-DVR (triangles) methods in accordance with the results provided in Table 2.
Symmetry 14 00445 g002
Figure 3. Fidelity obtained for each state displayed in Table 2 taking a step of Δ = 20 . The solid circles correspond to the HO-DVR method while the triangles stand for the PT-DVR approach.
Figure 3. Fidelity obtained for each state displayed in Table 2 taking a step of Δ = 20 . The solid circles correspond to the HO-DVR method while the triangles stand for the PT-DVR approach.
Symmetry 14 00445 g003
Figure 4. Particular wave functions provided by the HO-DVR and PT-DVR approaches taking equal basis dimension. The dashed lines correspond to the analytic wave functions. Values of the fidelity are included. At the right of the wave functions, the absolute value of the difference between the calculated and exact wave functions are also shown, reflecting a better description for PT-DVR.
Figure 4. Particular wave functions provided by the HO-DVR and PT-DVR approaches taking equal basis dimension. The dashed lines correspond to the analytic wave functions. Values of the fidelity are included. At the right of the wave functions, the absolute value of the difference between the calculated and exact wave functions are also shown, reflecting a better description for PT-DVR.
Symmetry 14 00445 g004
Figure 5. Errors associated with bound state q | Ψ 5 as a function of basis dimension N. The second row corresponds to a zoom of the results given in the first row.
Figure 5. Errors associated with bound state q | Ψ 5 as a function of basis dimension N. The second row corresponds to a zoom of the results given in the first row.
Symmetry 14 00445 g005
Figure 6. Comparison between the PT (left plot) and Morse (right plot) potentials with the simplest approximation in terms of piecewise potentials.
Figure 6. Comparison between the PT (left plot) and Morse (right plot) potentials with the simplest approximation in terms of piecewise potentials.
Symmetry 14 00445 g006
Figure 7. Representation of the PT and Morse potentials in terms of multistep-piecewise potentials.
Figure 7. Representation of the PT and Morse potentials in terms of multistep-piecewise potentials.
Symmetry 14 00445 g007
Figure 8. Correlation diagrams between the spectrum generated by the square well potential depicted in Figure 6 and the exact PT-potential as a function of the number of steps using the HO-DVR method (left figure) and the PT-DVR approach (right figure). In both cases, the basis dimension was taken to be N = 230 . The levels with gray color correspond to states in the square well potential that move to the continuum region as the number of steps are increased.
Figure 8. Correlation diagrams between the spectrum generated by the square well potential depicted in Figure 6 and the exact PT-potential as a function of the number of steps using the HO-DVR method (left figure) and the PT-DVR approach (right figure). In both cases, the basis dimension was taken to be N = 230 . The levels with gray color correspond to states in the square well potential that move to the continuum region as the number of steps are increased.
Symmetry 14 00445 g008
Figure 9. Correlation diagram between the spectrum generated by the half-space square well potential and the Morse potential as a function of the number of steps using the Morse-DVR approach. The left figure corresponds to the analysis using a basis dimension N = 150 , while in the right figure, N = 2500 is used. The levels with gray color correspond to states in the half square well potential that move to the continuum region as the number of steps increases.
Figure 9. Correlation diagram between the spectrum generated by the half-space square well potential and the Morse potential as a function of the number of steps using the Morse-DVR approach. The left figure corresponds to the analysis using a basis dimension N = 150 , while in the right figure, N = 2500 is used. The levels with gray color correspond to states in the half square well potential that move to the continuum region as the number of steps increases.
Symmetry 14 00445 g009
Figure 10. Irreducible representations associated with the energy levels of a particle confined in a square well potential. Systematic accidental degeneracy with irreps ( A 1 , B 1 ) and ( A 2 , B 2 ) has been remarked. See Leybraz et al. [9].
Figure 10. Irreducible representations associated with the energy levels of a particle confined in a square well potential. Systematic accidental degeneracy with irreps ( A 1 , B 1 ) and ( A 2 , B 2 ) has been remarked. See Leybraz et al. [9].
Symmetry 14 00445 g010
Figure 11. Selected symmetry projected wave functions associated with a particle in a square well potential. Both HO-DVR and PT-DVR provide the same plots at this level of resolution using the parameters associated with Table 3. The pair of wave functions are distributed in accordance with (a) accidental degeneracy with irreps ( A 1 , B 1 ), (b) natural degeneracy with components ( E A , E A ), (c) accidental degeneracy with irreps ( A 2 , B 2 ) and (d) totally symmetric single state.
Figure 11. Selected symmetry projected wave functions associated with a particle in a square well potential. Both HO-DVR and PT-DVR provide the same plots at this level of resolution using the parameters associated with Table 3. The pair of wave functions are distributed in accordance with (a) accidental degeneracy with irreps ( A 1 , B 1 ), (b) natural degeneracy with components ( E A , E A ), (c) accidental degeneracy with irreps ( A 2 , B 2 ) and (d) totally symmetric single state.
Symmetry 14 00445 g011
Figure 12. Absolute value of the difference between the exact and calculated wave functions for selected states. The differences show that the PT-DVR approach is more appropriate to describe this system. Hence, for HO-DVR, we have taken N = 52 , while for PT-DVR, we have N = 60 .
Figure 12. Absolute value of the difference between the exact and calculated wave functions for selected states. The differences show that the PT-DVR approach is more appropriate to describe this system. Hence, for HO-DVR, we have taken N = 52 , while for PT-DVR, we have N = 60 .
Symmetry 14 00445 g012
Figure 13. Root mean square deviation as a function of basis dimension N using both HO-DVR and PT-DVR approaches. In the upper panels, the basis dimension N 2 taken to obtain the results of Table 3 is pointed out. The average is included taking into account seven points to the right and to the left when possible. In the lower panels, the corresponding zooms are shown together with the percentage of points % P defined in Equation (76), where its correlation with rms is evident. Notice that several possibilities for the minimum in rms are possible. We have chosen the first absolute minimum in choosing the basis dimension.
Figure 13. Root mean square deviation as a function of basis dimension N using both HO-DVR and PT-DVR approaches. In the upper panels, the basis dimension N 2 taken to obtain the results of Table 3 is pointed out. The average is included taking into account seven points to the right and to the left when possible. In the lower panels, the corresponding zooms are shown together with the percentage of points % P defined in Equation (76), where its correlation with rms is evident. Notice that several possibilities for the minimum in rms are possible. We have chosen the first absolute minimum in choosing the basis dimension.
Symmetry 14 00445 g013
Figure 14. Correlation diagram of energy levels as a function of the potential depth.
Figure 14. Correlation diagram of energy levels as a function of the potential depth.
Symmetry 14 00445 g014
Figure 15. (a) A zoom-in of Figure 14 to show the correlation diagram of energy levels { A 2 , B 2 } and { A 1 , B 1 } as a function of the potential. In (b), the associated parameters ζ for these couple of levels are shown together with the corresponding energy splitting.
Figure 15. (a) A zoom-in of Figure 14 to show the correlation diagram of energy levels { A 2 , B 2 } and { A 1 , B 1 } as a function of the potential. In (b), the associated parameters ζ for these couple of levels are shown together with the corresponding energy splitting.
Symmetry 14 00445 g015
Figure 16. Rectangular well potential and symmetry elements associated with the geometrical symmetry C 4 v . Two reference systems are shown: (a) the origin is located at the left corner with solutions given by (78) and (b) the origin is located at the center of the rectangle.
Figure 16. Rectangular well potential and symmetry elements associated with the geometrical symmetry C 4 v . Two reference systems are shown: (a) the origin is located at the left corner with solutions given by (78) and (b) the origin is located at the center of the rectangle.
Symmetry 14 00445 g016
Figure 17. Irreducible representations associated with the energy levels of a particle confined in rectangular wells with n = 1 and m = 2 . Systematic accidental degeneracy with irreps ( A 2 , B 2 ) and ( B 2 , B 2 ) has been highlighted. The parenthesis ( n 1 n 2 ) corresponds to the labels associated with wave functions (78).
Figure 17. Irreducible representations associated with the energy levels of a particle confined in rectangular wells with n = 1 and m = 2 . Systematic accidental degeneracy with irreps ( A 2 , B 2 ) and ( B 2 , B 2 ) has been highlighted. The parenthesis ( n 1 n 2 ) corresponds to the labels associated with wave functions (78).
Symmetry 14 00445 g017
Figure 18. Root mean square deviation as a function of the basis dimension N using both HO-DVR and PT-DVR approaches. In the upper panels, basis dimension N 2 taken to obtain the results of Table 6 are pointed out. The average is included by taking into account seven points to the right and to the left when possible. In the lower panels, the corresponding zooms are shown together with the percentage of points % P defined in Equation (76), where its correlation with rms is evident.
Figure 18. Root mean square deviation as a function of the basis dimension N using both HO-DVR and PT-DVR approaches. In the upper panels, basis dimension N 2 taken to obtain the results of Table 6 are pointed out. The average is included by taking into account seven points to the right and to the left when possible. In the lower panels, the corresponding zooms are shown together with the percentage of points % P defined in Equation (76), where its correlation with rms is evident.
Symmetry 14 00445 g018
Figure 19. Wave functions associated with a particle in a rectangular well potential calculated with the PT-DVR approach with N = 60 , in accordance with the results of Table 6. In general, the HO-DVR approach does not provide wave functions of similar quality.
Figure 19. Wave functions associated with a particle in a rectangular well potential calculated with the PT-DVR approach with N = 60 , in accordance with the results of Table 6. In general, the HO-DVR approach does not provide wave functions of similar quality.
Symmetry 14 00445 g019
Figure 20. Action of reflection σ d a (with symmetry element in red) over the eigenfunctions with accidental degeneracy of the rectangular well potential. At the right of the wave functions the overlap between the original and the resulting wave functions are shown. Although they are coincidental, the space must be extended.
Figure 20. Action of reflection σ d a (with symmetry element in red) over the eigenfunctions with accidental degeneracy of the rectangular well potential. At the right of the wave functions the overlap between the original and the resulting wave functions are shown. Although they are coincidental, the space must be extended.
Symmetry 14 00445 g020
Figure 21. Correlation diagram of the energy levels as a function of the potential depth.
Figure 21. Correlation diagram of the energy levels as a function of the potential depth.
Symmetry 14 00445 g021
Figure 22. (a) A zoom of Figure 21 to show the correlation diagram of the energy levels { B 2 , B 2 } and { A 2 , B 2 } as a function of the potential. In (b), the associated parameters ζ for these couple of levels are shown together with the corresponding energy splitting.
Figure 22. (a) A zoom of Figure 21 to show the correlation diagram of the energy levels { B 2 , B 2 } and { A 2 , B 2 } as a function of the potential. In (b), the associated parameters ζ for these couple of levels are shown together with the corresponding energy splitting.
Symmetry 14 00445 g022
Figure 23. Correlation diagram of the energy levels for a given potential depth as a function of parameter τ defined in (96). The calculation was carried out with the PT-DVR method with parameters N = 60 , α x = 9 × 10 9 m 1 and α y = 4.5 × 10 9 m 1 . Degeneracy is pointed out with circles. However, for τ = 4 , degeneracy is not evident because of lack of convergence.
Figure 23. Correlation diagram of the energy levels for a given potential depth as a function of parameter τ defined in (96). The calculation was carried out with the PT-DVR method with parameters N = 60 , α x = 9 × 10 9 m 1 and α y = 4.5 × 10 9 m 1 . Degeneracy is pointed out with circles. However, for τ = 4 , degeneracy is not evident because of lack of convergence.
Symmetry 14 00445 g023
Table 1. Description of the 1D potential well with infinite walls with the algebraic DVR methods: HO-DVR, PT-DVR and S U ( 2 ) -UGA. The parameters for the potential were taken to be a = 4 Å and V 0 = 10 8 cm 1 , with μ = 1.6726219 × 10 27 kg. In addition, the fitted parameters associated with the algebraic DVR methods were found to be ω = 2.0 × 10 13 s 1 , α = 4.5 × 10 9 m 1 and ω = 1.6 × 10 12 s 1 for HO-DVR, PT-DVR and S U ( 2 ) -UGA, respectively. In each case, E c is the calculated energy in cm 1 . The zero of the energy was taken to be at the bottom of the well. The basis dimension N needed to obtain the displayed description was also included.
Table 1. Description of the 1D potential well with infinite walls with the algebraic DVR methods: HO-DVR, PT-DVR and S U ( 2 ) -UGA. The parameters for the potential were taken to be a = 4 Å and V 0 = 10 8 cm 1 , with μ = 1.6726219 × 10 27 kg. In addition, the fitted parameters associated with the algebraic DVR methods were found to be ω = 2.0 × 10 13 s 1 , α = 4.5 × 10 9 m 1 and ω = 1.6 × 10 12 s 1 for HO-DVR, PT-DVR and S U ( 2 ) -UGA, respectively. In each case, E c is the calculated energy in cm 1 . The zero of the energy was taken to be at the bottom of the well. The basis dimension N needed to obtain the displayed description was also included.
HO-DVRPT-DVR SU ( 2 ) -UGA
N = 58 N = 60 N = 2483
n Exact E c ( cm 1 ) % Error E c ( cm 1 ) % Error E c ( cm 1 ) % Error
110.3210.370.4510.320.0418.7381.41
241.3041.370.1641.160.3349.3519.50
392.9293.340.4592.880.04100.137.76
4165.19165.430.14164.620.35170.663.31
5258.11259.300.46258.020.03260.430.90
6371.67372.050.10370.280.37368.680.80
7505.89508.310.48505.750.03494.702.21
8660.75660.800.01657.970.42637.343.54
9836.26840.450.50836.110.02795.734.85
101032.431030.300.211027.250.50968.366.21
Table 2. Description of the 1D potential well using algebraic DVR methods: HO-DVR, PT-DVR and S U ( 2 ) -UGA. The parameters for the potential were taken to be a = 4 Å and V 0 = 1000 cm 1 . In this case the fitted parameters associated with the algebraic DVR methods were the same as in the previous case. The zero of the energy was located at the bottom of the well. The basis dimension N needed to obtained the displayed description was also included.
Table 2. Description of the 1D potential well using algebraic DVR methods: HO-DVR, PT-DVR and S U ( 2 ) -UGA. The parameters for the potential were taken to be a = 4 Å and V 0 = 1000 cm 1 . In this case the fitted parameters associated with the algebraic DVR methods were the same as in the previous case. The zero of the energy was located at the bottom of the well. The basis dimension N needed to obtained the displayed description was also included.
HO-DVRPT-DVR SU ( 2 ) -UGA
N = 58 N = 94 N = 2483
n Exact E c ( cm 1 ) % Error E c ( cm 1 ) % Error E c ( cm 1 ) % Error
19.119.020.999.070.4417.7995.28
236.4336.070.9936.250.4945.5825.12
381.8381.130.8681.510.3991.6411.99
4145.27144.030.85144.690.40155.557.08
5226.53224.900.72225.760.34236.754.51
6325.34323.160.67324.260.33334.442.80
7441.25439.070.49440.130.25447.651.45
8573.49571.190.40572.190.23575.080.28
9720.62719.820.11719.860.11714.900.79
10878.85878.850.00878.180.08863.811.71
Table 3. Results for one particle confined in a 2D square well potential using both HO-DVR and PT-DVR methods. The fitted parameters were chosen to be ω = 2.0 × 10 13 s 1 for the HO-DVR method and α = 4.5 × 10 9 m 1 for the PT-DVR method. The basis dimensions were chosen to obtain the absolute minimum of rms. Hence, for HO-DVR, we have taken N = 52 , while for PT-DVR, we have N = 60 .
Table 3. Results for one particle confined in a 2D square well potential using both HO-DVR and PT-DVR methods. The fitted parameters were chosen to be ω = 2.0 × 10 13 s 1 for the HO-DVR method and α = 4.5 × 10 9 m 1 for the PT-DVR method. The basis dimensions were chosen to obtain the absolute minimum of rms. Hence, for HO-DVR, we have taken N = 52 , while for PT-DVR, we have N = 60 .
HO-DVRPT-DVR
N 2 = 2704 N 2 = 3600
Irrep ( n 1 , n 2 ) Exact (cm 1 ) E c ( cm 1 ) % Error E c ( cm 1 ) % Error
A 1 (1,1)20.6520.620.1420.640.04
E(2,1)51.6251.570.1051.480.27
B 2 (2,2)82.5982.490.1282.320.33
A 1 (3,1)103.24103.250.01103.200.04
B 1 (3,1)103.24103.250.01103.200.04
E(3,2)134.22134.170.03134.040.13
E(4,1)175.51175.290.13175.940.24
A 1 (3,3)185.84186.850.54185.770.04
A 2 (4,2)206.48206.210.13205.780.34
B 2 (4,2)206.48206.210.13205.780.34
Table 4. Irreps spanned by the functions (78).
Table 4. Irreps spanned by the functions (78).
Irrep C 2 v n 1 n 2
A 1 oddodd
A 2 eveneven
B 1 evenodd
B 2 oddeven
Table 5. Character table for group C 2 v .
Table 5. Character table for group C 2 v .
C 2 v E C 2 σ v ( xz ) σ v ( yz )
A 1 1 1 1 1 z x 2 ; y 2 ; z 2
A 2 1 1 1 1 R z x y
B 1 1 1 1 1 x ; R y x z
B 2 1 1 1 1 y ; R x y z
Table 6. Results of the energies in cm 1 for one particle confined in a rectangular square well potential using both HO-DVR and PT-DVR methods. The dimensions of the rectangular well were chosen to be a = 1 Å and b = 2 Å with the fitted parameters ω x = 2.0 × 10 13 s 1 , ω y = 4.0 × 10 13 s 1 , and N = 73 for the HO-DVR method and α x = 9.0 × 10 9 , α x = 4.5 × 10 9 m 1 and N = 60 for the PT-DVR method.
Table 6. Results of the energies in cm 1 for one particle confined in a rectangular square well potential using both HO-DVR and PT-DVR methods. The dimensions of the rectangular well were chosen to be a = 1 Å and b = 2 Å with the fitted parameters ω x = 2.0 × 10 13 s 1 , ω y = 4.0 × 10 13 s 1 , and N = 73 for the HO-DVR method and α x = 9.0 × 10 9 , α x = 4.5 × 10 9 m 1 and N = 60 for the PT-DVR method.
HO-DVRPT-DVR
N 2 = 5329 N 2 = 3600
Irrep ( n 1 , n 2 ) Exact (cm 1 ) E c ( cm 1 ) % Error E c ( cm 1 ) % Error
A 1 ( 1 , 1 ) 51.6251.180.8551.600.04
B 2 ( 1 , 2 ) 82.5982.110.5882.440.19
A 1 ( 1 , 3 ) 134.22133.600.46134.160.04
B 1 ( 2 , 1 ) 175.51173.990.87174.960.31
B 2 ( 1 , 4 ) 206.49205.820.32205.900.28
A 2 ( 2 , 2 ) 206.49204.920.76205.800.33
B 1 ( 2 , 3 ) 258.11256.400.66257.520.23
A 1 ( 1 , 5 ) 299.40298.430.32299.300.03
A 2 ( 2 , 4 ) 330.38328.630.53329.250.34
A 1 ( 3 , 1 ) 382.00378.191.00381.850.04
B 2 ( 1 , 6 ) 413.00409.120.94411.560.34
B 2 ( 3 , 2 ) 413.00412.010.24412.690.07
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Suárez, E.; Santiago-Acosta, R.D.; Lemus, R. Algebraic DVR Approaches Applied to Piecewise Potentials: Symmetry and Degeneracy. Symmetry 2022, 14, 445. https://doi.org/10.3390/sym14030445

AMA Style

Suárez E, Santiago-Acosta RD, Lemus R. Algebraic DVR Approaches Applied to Piecewise Potentials: Symmetry and Degeneracy. Symmetry. 2022; 14(3):445. https://doi.org/10.3390/sym14030445

Chicago/Turabian Style

Suárez, Eduardo, Rubén Dario Santiago-Acosta, and Renato Lemus. 2022. "Algebraic DVR Approaches Applied to Piecewise Potentials: Symmetry and Degeneracy" Symmetry 14, no. 3: 445. https://doi.org/10.3390/sym14030445

APA Style

Suárez, E., Santiago-Acosta, R. D., & Lemus, R. (2022). Algebraic DVR Approaches Applied to Piecewise Potentials: Symmetry and Degeneracy. Symmetry, 14(3), 445. https://doi.org/10.3390/sym14030445

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop