Harmonic Oscillator with Additional Repulsive Cubic Force: Solutions and Study

AI Thread Summary
The discussion centers on solving a harmonic oscillator problem with an additional repulsive cubic force, represented by the potential U(q1)=k/2*q1² - k'q1⁴. The original poster is seeking help with the integral solution for the system, having attempted trigonometric substitution without success. A response suggests transforming the variables to simplify the equation into a form resembling the Jacobi elliptic integral of the first kind. This approach may provide a pathway to find both periodic and non-periodic solutions. The conversation highlights the challenges of tackling complex integrals in classical mechanics.
Miesvama
Messages
2
Reaction score
0
Hi all,

this is my first time on PF.

I do not know English but I have a problem of a harmonic oscillator.
I have rather large head, help me please , I do not know what else to do ...
I have this problem:

Consider the harmonic oscillator with an additional repulsive
cubic force, whose potential is U(q1)=\frac{k}{2}*q1^{2} - k'q1^{4}, (k, k > 0), and study all
possible solutions, periodic and non-periodic.

I do know the Hamiltonian and the equation solution of the system, giving

q1*=\sqrt{1/2}\int(dq1/\sqrt{(E/m)-U(q1)})

I tried to do it by trigonometric substitution but does not work, i do not know if anyone could give me some idea of how I can solve, I'll be very grateful.
 
Physics news on Phys.org
Miesvama said:
Hi all,

this is my first time on PF.

I do not know English but I have a problem of a harmonic oscillator.
I have rather large head, help me please , I do not know what else to do ...
I have this problem:

Consider the harmonic oscillator with an additional repulsive
cubic force, whose potential is U(q1)=\frac{k}{2}*q1^{2} - k'q1^{4}, (k, k > 0), and study all
possible solutions, periodic and non-periodic.

I do know the Hamiltonian and the equation solution of the system, giving

q1*=\sqrt{1/2}\int(dq1/\sqrt{(E/m)-U(q1)})

I tried to do it by trigonometric substitution but does not work, i do not know if anyone could give me some idea of how I can solve, I'll be very grateful.

Hi. Nobody else has responded to your question, so I guess I'll give it a try.

The equation that I think you meant is this:

t = \sqrt{m/2}\int(dq1/\sqrt{E-U(q1)})

If we change to new variables x and s where
x = A q1 and s = B t, where A and
B are constants, we can choose the constants to make the
equation look like this:

t = \int(dx/\sqrt{(1-x^2)(1- \lambda^2 x^2)}

where \lambda is yet another constant.

That's called the "Jacobi form of the incomplete elliptic integral of the first kind", F(s;\lambda).

http://en.wikipedia.org/wiki/Jacobi..._of_nonlinear_ordinary_differential_equations
 
Thank you.. :D
 
Last edited:
Thread ''splain this hydrostatic paradox in tiny words'
This is (ostensibly) not a trick shot or video*. The scale was balanced before any blue water was added. 550mL of blue water was added to the left side. only 60mL of water needed to be added to the right side to re-balance the scale. Apparently, the scale will balance when the height of the two columns is equal. The left side of the scale only feels the weight of the column above the lower "tail" of the funnel (i.e. 60mL). So where does the weight of the remaining (550-60=) 490mL go...
Consider an extremely long and perfectly calibrated scale. A car with a mass of 1000 kg is placed on it, and the scale registers this weight accurately. Now, suppose the car begins to move, reaching very high speeds. Neglecting air resistance and rolling friction, if the car attains, for example, a velocity of 500 km/h, will the scale still indicate a weight corresponding to 1000 kg, or will the measured value decrease as a result of the motion? In a second scenario, imagine a person with a...
Scalar and vector potentials in Coulomb gauge Assume Coulomb gauge so that $$\nabla \cdot \mathbf{A}=0.\tag{1}$$ The scalar potential ##\phi## is described by Poisson's equation $$\nabla^2 \phi = -\frac{\rho}{\varepsilon_0}\tag{2}$$ which has the instantaneous general solution given by $$\phi(\mathbf{r},t)=\frac{1}{4\pi\varepsilon_0}\int \frac{\rho(\mathbf{r}',t)}{|\mathbf{r}-\mathbf{r}'|}d^3r'.\tag{3}$$ In Coulomb gauge the vector potential ##\mathbf{A}## is given by...
Back
Top