Anharmonic oscillators with closed-form solutions

  • #1
hilbert2
Science Advisor
Insights Author
Gold Member
1,600
607
TL;DR Summary
Some quantum oscillators with higher than quadratic V(x) seem to have an exact solution. Does this have any applications in QFT?
There are some articles from the 1980s where the authors discuss 1D quantum oscillators where has higher than quadratic terms in it but an exact solution can still be found. One example is in this link:

https://iopscience.iop.org/article/10.1088/0305-4470/14/9/001

Has anyone tried to form a quantum field where the normal modes have a behavior similar to this kind of oscillator, and write the ground state of this interacting field in terms of the ground state and creation operators of the non-interacting case? It would likely not have any practical use, but it would just be a QFT system where something can be solved non-perturbatively.
 
  • Like
Likes vanhees71 and Demystifier
Physics news on Phys.org
  • #3
hilbert2 said:
Has anyone tried to form a quantum field where the normal modes have a behavior similar to this kind of oscillator, and write the ground state of this interacting field in terms of the ground state and creation operators of the non-interacting case? It would likely not have any practical use, but it would just be a QFT system where something can be solved non-perturbatively.
There are quite a number of exactly solvable relativistic QFTs in 2 spacetime dimensions. The most important one is perhaps the Schwinger model, which is QED in 2 instead of 4 dimensions.

Quantum wires are described by nonrelativistic QFTs, some of which are exactly solvable and explain (in simplified models) the fractional quantum Hall effect.
 
  • Like
Likes vanhees71 and hilbert2
  • #4
Keith_McClary said:
That one seems to be on the half line.

Here is one in 2D:
https://arxiv.org/pdf/quant-ph/0003100.pdf
Oh, I didn't notice that right away. A system defined only on one half-line of the position axis wouldn't be very useful, because the degrees of freedom of the corresponding field would have only positive values then.

A. Neumaier said:
There are quite a number of exactly solvable relativistic QFTs in 2 spacetime dimensions. The most important one is perhaps the Schwinger model, which is QED in 2 instead of 4 dimensions.

Quantum wires are described by nonrelativistic QFTs, some of which are exactly solvable and explain (in simplified models) the fractional quantum Hall effect.
Thanks for the reply, I didn't remember the Schwinger model. A solvable scalar field system where the Lagrangian density contains a 6th or higher power of the field function would be in some way interesting because higher than phi-fourth interactions are not renormalizable (if I understood D. Tong's lecture notes correctly) and with an exact solution you could still see how the system behaves.
 
  • #5
hilbert2 said:
A solvable scalar field system where the Lagrangian density contains a 6th or higher power of the field function would be in some way interesting because higher than phi-fourth interactions are not renormalizable (if I understood D. Tong's lecture notes correctly) and with an exact solution you could still see how the system behaves.
In 2D all bosonic fields with polynomial interactions are renormalizable, and existence has been proved rigorously if the interactions are bounded below. See the book by Glimm and Jaffe.

But for field theories, sextic interactions are not exactly solvable; at least I don't know of an example. Typically, the exactly solvable theories have very special quartic interactions.

The 2D Gross-Neveu Model is nonrenormalizable and exactly solvable, I believe.
 
Last edited:
  • Like
Likes vanhees71
  • #6
A. Neumaier said:
The 2D Gross-Neveu Model is nonrenormalizable and exactly solvable, I believe.
How can a theory be both nonrenormalizable and exactly solvable at the same time? Does it mean that it is ill defined perturbatively but well defined non-perturbatively?
 
  • Like
Likes vanhees71
  • #7
Demystifier said:
How can a theory be both nonrenormalizable and exactly solvable at the same time? Does it mean that it is ill defined perturbatively but well defined non-perturbatively?
Yes if you adhere to the old belief that nonrenormalizable theories are perturbatively undefined.

But the modern view is that nonrenormalizable theories are as well-defined perturbatively as renormalizable ones, except that they depend on a countable (rather than finite) number of parameters. This is not a real defect since analytic functions also depend on in finitely many parameters but as long as the power series converge fast enough a few known parameters produce accurate results.

Thus the most correct statement is that nonrenormalizable and exactly solvable theories are underspecified perturbatively but well defined non-perturbatively.
 
  • Like
Likes vanhees71 and Demystifier
  • #8
I believe that the way how the authors of the linked articles have found the exact solutions for those oscillators has been by rearranging the nondimensional () Schrödinger equation



to the inverse problem equation

(*)

and then by noting that for instance in the case the resulting potential energy is

,

which is a sixth-degree polynomial. So, the ground state of the system with having the form of that polynomial is like a Gaussian but with higher than second powers of in the exponent. However, if the ground state is like that, then none of the excited state wave functions have the simple form of the same exponential multiplied with polynomials (like in the harmonic oscillator case, where those multipliers are Hermite polynomials). Actually, I'm able to prove this in the following way (tell me if there's a logical error somewhere):

Suppose the unnormalized ground state wave function of an anharmonic oscillator is , where is some polynomial of containing only terms of order divisible by 2. Then, according to (*), the 1D potential energy for which this is one of the energy eigenfunctions is

.

The constant is the ground state energy. Now, if another eigenfunction of the same Hamiltonian is with a polynomial, the function should obey the differential equation (derived from (*))



where is a constant that doesn't depend on . Suppose the polynomial is of :th order and is of :th order. Then the term on the LHS that contains highest powers of is the , where the order is . The highest power of on the RHS is of :th order. But the terms and can't be of same order in unless , as in the harmonic oscillator. Therefore, none of the excited states of this kind of anharmonic oscillators can't have the form of a simple finite-degree polynomial multiplying the ground state wave function.

Edit: The term is of order , not as in the unedited post.
 
Last edited:
  • Like
Likes vanhees71
  • #9
hilbert2 said:
I believe that the way how the authors of the linked articles have found the exact solutions for those oscillators has been by rearranging the nondimensional () Schrödinger equation



to the inverse problem equation

(*)

and then by noting that for instance in the case the resulting potential energy is

,

which is a sixth-degree polynomial.
This approach only works for a single oscillator. But you asked about quantum fields, where one needs infinitely many of them!
 
  • Like
Likes vanhees71
  • #10
A. Neumaier said:
This approach only works for a single oscillator. But you asked about quantum fields, where one needs infinitely many of them!
Yes, if the Fourier modes of the self-interacting field behave like anharmonic oscillators with different "spring constants", not all of them necessarily have an exact solution for the ground state energy and wave function.

Edit: It's possible that I'm considering field systems in a too simplified way here; if some scalar field equation with a nonlinear interaction has a Fourier transform

,

where is some parameter, then the original untransformed field equation is not necessarily an acceptable local equation of motion. Or is it?
 
Last edited:
  • Like
Likes vanhees71
  • #11
Using the shooting method to calculate first 10 energy eigenvalues of some pure anharmonic oscillators

,

with particle mass and set to value 1, I found out that a function like fits to the data points (th energy eigenvalue) if the quantum number is seen as a continuum variable. For multiplier and exponents to , the curve fits look like in this graph (note that the ground state of the harmonic oscillator case is also labelled , not the usual ):

power_law_fits_for_anharmonic_eigenvalues.jpg


The fitting parameters , and are shown in this table:

power_law_fit_to_anharmonic_energies_table.jpg


So this seems to be kind of an exact result, but not completely as you have to determine the parameters empirically for each oscillator.

I initially though I was the first to notice this feature of those energy levels, but yesterday I found this article preprint made by researchers in Iran in 2019 and describing the same thing:

https://arxiv.org/pdf/1910.02746.pdf

Is anyone else familiar with this power-law property of the pure anharmonic oscillator? The exponent clearly approaches when , because an anharmonic oscillator of infinite order is practically the same as a 1D particle-in-box system with sudden infinite walls around it.

Actually, that same trial function also fits to data points describing the spectrum of an oscillator with both a quadratic and quartic term in , but then the large variety of combinations of parameter values makes it difficult to find a general logic behind the results.
 
  • Like
Likes Keith_McClary and vanhees71
Back
Top