Research Article  Open Access
Cornell Potential: A Neural Network Approach
Abstract
We solved Schrödinger equation with Cornell potential (Coulombpluslinear potential) by using neural network approach. Four different cases of Cornell potential for different potential parameters were used without a physical relevance. Besides that charmonium, bottomonium and bottomcharmed spinaveraged spectra were also calculated. Obtained results are in good agreement with the reference studies and available experimental data.
1. Introduction
The Schrödinger equation with Coulombpluslinear potential (Cornell potential) has received a great deal of attention as an important nonrelativistic model for the study of quarkantiquark systems, namely, mesons [1–18]. Especially Coulombpluslinear potential was the first potential model to study heavy quarkonium systems and inspired other phenomenological models. The quarkonium spectrum is a substantial field to improve our understanding about nature of QCD. The Cornell potential reads aswhere and are some constants. At small quarkantiquark distances, the Cornell potential goes like which is known as Coulomblike term and at large distances it goes as , which is known as linear term. The Coulomblike term arises from onegluon exchange and linear term presumably arises from higher order effects [19]. Up to now, the linear term has not been calculated from the first principles of QCD.
Aside from its physical relevance, the Schrödinger equation with Coulombpluslinear potential has been studied with pure mathematical techniques. Hall used the method of potential envelopes to construct general upper and lower bounds on the eigenvalues of the Hamiltonian representing a single particle in Coulombpluslinear potential [20]. In [21], Chaudhuri et al. used Hill determinant method to bound state eigenvalue problem with Coulombpluslinear potential. Plante and Antippa solved the Schrödinger equation for a quarkantiquark system interacting via a Coulombpluslinear potential and obtained the wave functions as power series, with their coefficients given in terms of the combinatorics functions [22]. In [23], the authors presented a numerically precise treatment of the Crank–Nicolson method with an imaginary time evolution operator in order to solve the Schrödinger equation. Although such models have been studied to solve Schrödinger equations numerically or analytically, exact solutions of Schrödinger equation with Cornell potential are still unknown.
Artificial neural networks (ANNs) can be used as an alternative method to solve both ordinary and partial differential equations. ANN is a parallel distributed processor in which numerous numbers of simply designed computing units exist. These units are called neurons. Its massive connections between neurons can be used for storing various types of informations, particularly the ones classified as knowledge or experience [24, 25]. This information is acquired as interneuron connection strengths, or synaptic weights.
Especially with the work of [26, 27], ANNs are being used widely for solving ordinary and partial differential equations. ANNs have many advantages compared to existing semianalytic and numerical methods. The main advantages of using neural networks to solve differential equations can be stated as follows [28]:(i)Neural networks (NNs) are universal functions approximators which are useful in solving differential equations.(ii)Trial solutions of artificial neural networks involve a single independent variable regardless of the dimension of the problem.(iii)The approximate solutions are continuous over all the domain of integration. In contrast, the numerical methods provide solutions only over discrete points and the solution between these points must be interpolated.(iv)The computational complexity does not increase considerably with the number of sampling points and with the number of dimensions in problem.
Being an eigenvalue problem, Schrödinger equation had got attention by ANNs at the beginning of 90s. In [29], the authors taught a neural network the solution of the two dimensional Schrödinger equation for some model potential energy functions. In work of Androsiuk et al. [30], they presented computer simulations of a neural network capable of learning to perform transformations generated by the Schrödinger equation required to find eigenenergies of two dimensional harmonic oscillator. Sugawara presented a new approach for solving the Schrödinger equation based on genetic algorithm and artificial neural network. The method was tested in the calculation of onedimensional harmonic oscillator and other model potentials such as Morse potential [31].
In the present paper, we solve Schrödinger equation with Coulombpluslinear potential via ANN system. In Section 2, we introduce the formalism of ANN method and describe how it can be applied to the quantum mechanical calculations. In Section 3, we give and discuss numerical results for the eigenvalues of Coulombpluslinear potential and in Section 4 we summarize our findings.
2. Formalism
The neural network (NN) is constructed as a model of simply designed computing unit, called neuron. Figure 1 illustrates a simple model of a single neuron with multiple inputs and one output.
All the input signals are summed up as and the output signal is determined by the nonlinear activation function . In this work we use a sigmoid functionas an activation function since it is possible to derive all the derivatives of in terms of itself. This differentiability is an important aspect for the Schrödinger equation. Here we use a multilayer perceptron neural network (MLPN). The architecture of MLPN is shown in Figure 2.
This is an example of feed forward neural network. Feed forward neural networks are the most popular architectures due to their structural flexibility, good representational capabilities, and availability of a large number of training algorithm [24]. This network consists of neurons arranged in layers in which every neuron is connected to all neurons of the next layer and the information process can be in only one way, from input to the output. The inputoutput properties of the neurons (perceptrons in computerized systems) can be written as where is for input, is for hidden, and is for output layers. Input to the perceptrons is given as Here, and represent the numbers of the units belonging to input and hidden layers, is the synaptic weight parameter which connects the neurons and , and represents threshold parameter for the neuron [30]. The overall response of the network is given asThe derivatives of with respect to network parameters (weights and thresholds) can be done by differentiating (9) as Neural network systems are selftraining systems. Therefore parameters of training need updating in the transfer process and these equations play important roles on the learning processes of the neural networks.
2.1. Implementation of ANN to Quantum System
We consider the application of the ANN to a quantum mechanical calculation. We will follow the formalism which was developed in [26]. Consider the following differential equation:where is a linear operator, is a known function, and at the boundaries. In order to solve this differential equation a trial function can be written as which employs a feed forward neural network with parameter vectors and to be adjusted. The parameter vector corresponds to the weights and biases of the neural architecture. The functions and should be specified in a convenient way so that satisfies the boundary conditions regardless of the and values. To obtain a solution for (14), the collocation method can be used and the differential equation can be transformed into a minimization problemFor Schrödinger equation, (14) takes the form with the boundary condition . In this case, the trial solution can be written as where at boundary conditions for a range of values. By discretizing the domain of the problem, it is transformed into a minimization problem with respect to the parameters and where is the error function and can be computed as
Let us consider a multilayer perceptron with input units, one hidden layer with units and one output. This multilayer perceptron looks like Figure 2 but with just one output. For a given input vector the output of the network is where Here, denotes the weight from the input unit to the hidden unit , is the weight from the hidden unit to the output, is the bias of hidden unit , and is the sigmoid function, given in (2). The derivatives of network output can be defined as where and is the th order derivative of the sigmoid.
Once the derivatives of the error with respect to the network parameters have been defined, any minimization technique can be carried out. In this work, we used a feed forward neural network with a back propagation algorithm.
By employing this approach it is possible to obtain energy eigenvalues of the Cornell potential. Before obtaining Cornell potential eigenvalues, it will be useful to give an example.
2.1.1. A Warm up Example: Yukawa Potential
The Yukawa potential has an important role in various branches of physics. For example, in plasma physics, it is known as the DebyeHückel potential; in solidstate physics and atomic physics, it is called the ThomasFermi or the screened Coulomb potential and has a role in the nucleonnucleon interaction arising out of the onepionexchange mechanism in nuclear physics [34]. The Yukawa potential reads as where the parameters and are given by different expressions, depending on the problem under consideration. Yukawa potential does not admit an exact solution and therefore various approximate analytic and numerical methods. The ground state energy for , , and in reported by [34] is , by [35] is , and by [36] is .
We parametrize trial function as where denotes a feed forward artificial neural network with one hidden layer and sigmoid hidden units The minimization problem is The energy eigenvalue obtained with this scheme is which is in good agreement with the numerical results in first five digits.
3. Numerical Results and Discussion
The radial part of the Schrödinger equation for Coulombpluslinear potential iswhere and . We trained the network with 200 equidistance points with in the interval and solved the Schrödinger equation in four cases.
3.1. Case 1:
Putting into (29) we get The Hamiltonian of this equation is By applying the procedure which was mentioned in the previous section, we get eigenvalues of this potential. In Tables 1 and 2, we present eigenvalues of the potential for and states, respectively.
It can be seen from Tables 1 and 2 that eigenvalues of the radially and orbitally excited states by ANN method are in good agreement with the reference studies.
3.2. Case 2:
The Hamiltonian for this potential is In Tables 3 and 4, we present eigenvalues of the potential for and states, respectively.
It can be seen from Tables 3 and 4 that eigenvalues of the radially and orbitally excited states by ANN method are in good agreement with the reference studies.
3.3. Case 3:
The Hamiltonian for this potential is In Tables 5 and 6, we present eigenvalues of the potential for and states, respectively.
It can be seen from Tables 5 and 6 that eigenvalues of the radially and orbitally excited states by ANN method are in good agreement with the reference studies. In this case, numerical results of orbitally excited sates can differ in first two or three digits. The reason for this difference can be the dominance of the linear part of the potential compared to the Coulomb part.
3.4. Case 4:
The Hamiltonian for this potential is In Table 7, we present ground state eigenvalues of the potential for varying values. It can be seen from Table 7 that obtained results are in good agreement with the given studies.

3.5. Heavy Quarkonium Spectra
Quarkonium systems are an ideal area for clarifying our understanding of QCD. They probe nearly all the energy regimes of QCD from high energy region to low energy region. In the high energy region, an expansion of coupling constant is possible and perturbative QCD is applicable. In the low energy region, such an expansion is not possible in the coupling constant and therefore nonperturbative methods need to be used. Besides that nonrelativistic QCD (NRQCD) approximation is also used for spectroscopy, decay, and production of heavy quarkonium [37, 38]. For an overview of NRQCD, see [39].
In this section we obtained spinaveraged mass spectra charmonium, bottomonium, and bottomcharmed system by solving nonrelativistic Schrödinger equation. It is possible to obtain full spectra for quarkonium systems including relativistic effects, spinspin, and spinorbit interactions. Since most of these contributions are really small compared to the given potential, even by neglecting those effects one can find results that are close to the experimental data.
The related Cornell potential is [40] with
In Tables 8, 9, and 10, charmonium, bottomonium, and bottomcharmed spectra are presented, respectively.
As can be seen in Tables 8, 9, and 10 charmonium, bottomonium, and bottomcharmed spectra of ANN agree well with given reference and available experimental data.
4. Summary
In this paper, we applied ANN method to deal with the solution of the Schrödinger equation with Coulombpluslinear potential. This potential belongs to the nonsolvable potentials class which has an exactly/analytically solvable part, together with a modifying term. We obtained the eigenvalues of Schrödinger equation and spinaveraged mass spectra of charmonium, bottomonium, and bottomcharmed systems. The obtained eigenvalues and heavy quarkonium spectra are in agreement with the theoretical studies and available experimental data.
The feed forward ANNs method has a good property of function approximation. A function approximation problem is to select or find a function among a welldefined functions set that closely matches a target function in a task specific way. This form employs a feed forward neural network as the basic approximation element, whose parameters (weights and biases) are adjusted to minimize an appropriate error function. In this study, the wave function is represented by the feed forward ANN and its inputs are taken as coordinate values. A trial solution is written as a feed forward neural network which contains adjustable parameters (weights and biases) and eigenvalue is refined to the known solutions by training the neural network.
This study would be useful for the exact or quasiexact spectra of a few body systems. It is also possible in principle to handle many body problems but such problems will impose much more heavier computational load and other difficulties such as convergence of sigmoid functions [41] and availability of hardware.
Data Availability
No data were used to support this study.
Conflicts of Interest
The author declares that they have no conflicts of interest.
References
 E. Eichten, K. Gottfried, T. Kinoshita, J. Kogut, K. D. Lane, and T.M. Yan, “Spectrum of charmed quarkantiquark bound states,” Physical Review Letters, vol. 34, no. 6, pp. 369–372, 1975. View at: Publisher Site  Google Scholar
 E. Eichten, K. Gottfried, T. Kinoshita, K. D. Lane, and T. M. Yan, “Interplay of confinement and decay in the spectrum of charmonium,” Physical Review Letters, vol. 36, p. 500, 1976. View at: Publisher Site  Google Scholar
 E. Eichten, K. Gottfried, T. Kinoshita, K. D. Lane, and T.M. Yan, “Charmonium: The model,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 17, no. 11, pp. 3090–3117, 1978. View at: Publisher Site  Google Scholar
 C. Quigg and J. L. Rosner, “Quantum mechanics with applications to quarkonium,” Physics Reports, vol. 56, no. 4, pp. 167–235, 1979. View at: Publisher Site  Google Scholar  MathSciNet
 S. Godfrey and N. Isgur, “Mesons in a relativized quark model with chromodynamics,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 32, no. 1, pp. 189–231, 1985. View at: Publisher Site  Google Scholar
 R. Roychoudhury, Y. P. Varshni, and M. Sengupta, “Family of exact solutions for the Coulomb potential perturbed by a polynomial in ,” Physical Review A: Atomic, Molecular and Optical Physics, vol. 42, no. 1, pp. 184–192, 1990. View at: Publisher Site  Google Scholar  MathSciNet
 L. P. Fulcher, “Matrix representation of the nonlocal kinetic energy operator, the spinless Salpeter equation and the Cornell potential,” Physical Review D, vol. 50, p. 447, 1994. View at: Publisher Site  Google Scholar
 T. Barnes, S. Godfrey, and E. S. Swanson, “Higher charmonia,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 72, article 054026, 2005. View at: Publisher Site  Google Scholar
 H. S. Chung, J. Lee, and D. Kang, “Cornell potential parameters for Swave heavy quarkonia,” Journal of the Korean Physical Society, vol. 52, pp. 1151–1154, 2008. View at: Publisher Site  Google Scholar
 C. O. Dib and N. A. Neill, “ splitting predictions in potential models,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 86, Article ID 094011, 2012. View at: Google Scholar
 M. Hamzavi and A. A. Rajabi, “Scalar–vector–pseudoscalar Cornell potential for a spin1/2 particle under spin and pseudospin symmetries: 1 + 1 dimensions,” Annals of Physics, vol. 334, pp. 316–320, 2013. View at: Publisher Site  Google Scholar
 J.K. Chen, “Spectral method for the Cornell and screened Cornell potentials in momentum space,” Physical Review D, vol. 88, Article ID 076006, 2013. View at: Google Scholar
 J. Alford and M. Strickland, “Charmonia and bottomonia in a magnetic field,” Physical Review D, vol. 88, Article ID 105017, 2013. View at: Google Scholar
 P. W. M. Evans, C. R. Allton, and J.I. Skullerud, “Ab initio calculation of finitetemperature charmonium potentials,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 89, Article ID 071502, 2014. View at: Publisher Site  Google Scholar
 R. L. Hall and N. Saad, “Schrödinger spectrum generated by the Cornell potential,” Open Physics Journal, vol. 13, pp. 83–89, 2015. View at: Google Scholar
 A. Vega and J. Flores, “Heavy quarkonium properties from Cornell potential using variational method and supersymmetric quantum mechanics,” Pramana, vol. 87, p. 73, 2016. View at: Publisher Site  Google Scholar  MathSciNet
 B. Gönül and Y. Cançelik, “Remarks on the treatments of nonsolvable potentials,” Turkish Journal of Physics, vol. 41, pp. 326–336, 2017. View at: Google Scholar
 N. R. Soni, B. R. Joshi, R. P. Shah, H. R. Chauhan, and J. N. Pandya, “ spectroscopy using the Cornell potential,” The European Physical Journal C, vol. 78, no. 7, 2018. View at: Publisher Site  Google Scholar
 D. B. Lichtenberg, “Energy levels of quarkonia in potential models,” International Journal of Modern Physics A, vol. 2, no. 6, pp. 1669–1705, 1987. View at: Publisher Site  Google Scholar
 R. L. Hall, “Simple eigenvalue formula for the Coulombpluslinear potential,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 30, no. 2, pp. 433–436, 1984. View at: Publisher Site  Google Scholar
 R. N. Chaudhuri, M. Tater, and M. Znojil, “The Hill determinant approach to the Coulomb plus linear confinement,” Journal of Physics A: Mathematical and General, vol. 20, no. 6, pp. 1401–1410, 1987. View at: Publisher Site  Google Scholar  MathSciNet
 G. Plante and A. F. Antippa, “Analytic solution of the Schrödinger equation for the Coulombpluslinear potential. I. The wave functions,” Journal of Mathematical Physics, vol. 46, no. 6, Article ID 062108, 2005. View at: Publisher Site  Google Scholar  MathSciNet
 D. Kang and E. Won, “Precise numerical solutions of potential problems using the Crank–Nicolson method,” Journal of Computational Physics, vol. 20, Article ID 2970, 2008. View at: Google Scholar
 S. Haykin, Neural Networks: A Comprehensive Foundation, PrenticeHall, New York, NY, USA, 2nd edition, 1999, Neural Networks: A Comprehensive Foundation.
 M. T. Spinning, J. A. Darsey, B. G. Sumpter, and D. W. Noid, “Opening up the black box of artificial neural networks,” Journal of Chemical Education, vol. 71, Article ID 406, 1994. View at: Google Scholar
 I. E. Lagaris, A. Likas, and D. I. Fotiadis, “Artificial neural network methods in quantum mechanics,” Computer Physics Communications, vol. 104, no. 13, pp. 1–14, 1997. View at: Publisher Site  Google Scholar
 I. E. Lagaris, A. Likas, and D. I. Fotiadis, “Artificial neural networks for solving ordinary and partial differential equations,” IEEE Transactions on Neural Networks and Learning Systems, vol. 9, no. 5, pp. 987–1000, 1998. View at: Publisher Site  Google Scholar
 N. Yadav, A. Yadav, and M. Kumar, An Introduction to Neural Network Methods for Differential Equations, SpringerBriefs in Applied Sciences and Technology, Springer, 2015. View at: Publisher Site  MathSciNet
 J. A. Darsey, D. W. Noid, and B. R. Upadhyaya, “Application of neural network computing to the solution for the groundstate eigenenergy of twodimensional harmonic oscillators,” Chemical Physics Letters, vol. 177, pp. 189–194, 1991. View at: Publisher Site  Google Scholar
 J. Androsiuk, L. Kutak, and K. Sienicki, “Neural network solution of the Schrödinger equation for a twodimensional harmonic oscillator,” Chemical Physics, vol. 173, pp. 377–383, 1993. View at: Publisher Site  Google Scholar
 M. Sugawara, “Numerical solution of the Schrödinger equation by neural network and genetic algorithm,” Computer Physics Communications, vol. 140, no. 3, pp. 366–380, 2001. View at: Google Scholar
 H. Ciftci and H. F. Kisoglu, “Nonrelativistic arbitrary lstates of quarkonium through asymptotic iteration method,” Advances in High Energy Physics, vol. 2018, Article ID 4549705, 7 pages, 2018. View at: Google Scholar
 C. Patrignani et al., “Review of particle physics,” Chinese Physics C, vol. 40, Article ID 100001, 2016. View at: Google Scholar
 M. Karakoc and I. Boztosun, “Accurate iterative and perturbative solutions of the yukawa potential,” International Journal of Modern Physics E, vol. 15, no. 6, pp. 1253–1262, 2006. View at: Google Scholar
 B. Gönül and K. Köksal, “Equivalence of two alternative approaches to Schrödinger equations,” Physica Scripta. An International Journal for Experimental and Theoretical Physics, vol. 73, no. 6, pp. 629–631, 2006. View at: Publisher Site  Google Scholar  MathSciNet
 C. Lee, “Equivalence of logarithmic perturbation theory and expansion of the superpotential in supersymmetric quantum mechanics,” Physics Letters A, vol. 267, no. 23, pp. 101–108, 2000. View at: Publisher Site  Google Scholar  MathSciNet
 B. A. Thacker and G. P. Lepage, “Heavyquark bound states in lattice QCD,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 43, no. 1, pp. 196–208, 1991. View at: Publisher Site  Google Scholar
 G. P. Lepage, L. Magnea, C. Nakhleh, U. Magnea, and K. Hornbostel, “Improved nonrelativistic QCD for heavyquark physics,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 46, no. 9, pp. 4052–4067, 1992. View at: Publisher Site  Google Scholar
 J. Soto, “Overview of nonrelativistic QCD,” The European Physical Journal A, vol. 31, pp. 705–710, 2007. View at: Google Scholar
 L. P. Fulcher, Z. Chen, and K. C. Yeong, “Energies of quarkantiquark systems, the Cornell potential, and the spinless Salpeter equation,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 47, no. 9, pp. 4122–4132, 1993. View at: Publisher Site  Google Scholar
 A. R. Barron, “Universal approximation bounds for superpositions of a sigmoidal function,” IEEE Transactions on Information Theory, vol. 39, no. 3, pp. 930–945, 1993. View at: Publisher Site  Google Scholar  MathSciNet
Copyright
Copyright © 2019 Halil Mutuk. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The publication of this article was funded by SCOAP^{3}.