Movimiento bajo una fuerza inversamente proporcional al cubo de la distancia

Posición, velocidad y aceleración en coordenadas polares

La posición de una partícula es (x, y) en coordenadas rectangulares y (r, θ) en coordenadas polares. La relación es

x=r·cosθ, x=r·sinθ

Expresamos la velocidad de la partícula en coordenadas polares

v = d r dt = r ^ dr dt +r d r ^ dt

Calculamos las componentes rectangulares de los vectores unitarios

r ^ = i ^ cosθ+ j ^ sinθ θ ^ = i ^ sinθ+ j ^ cosθ

Derivamos con respecto del tiempo

d r ^ dt =( i ^ sinθ+ j ^ cosθ) dθ dt = θ ^ dθ dt d θ ^ dt =( i ^ cosθ j ^ sinθ) dθ dt = r ^ dθ dt

La expresión del vector velocidad en coordenadas polares es

v = dr dt r ^ +r dθ dt θ ^

Derivando el vector velocidad con respecto del tiempo, obtenemos el vector aceleración

a = d v dt = d 2 r d t 2 r ^ + dr dt d r ^ dt +( r d 2 θ d t 2 + dr dt dθ dt ) θ ^ +r dθ dt d θ ^ dt = d 2 r d t 2 r ^ + dr dt dθ dt θ ^ +( r d 2 θ d t 2 + dr dt dθ dt ) θ ^ r ( dθ dt ) 2 r ^ = ( d 2 r d t 2 r ( dθ dt ) 2 ) r ^ +( r d 2 θ d t 2 +2 dr dt dθ dt ) θ ^

Fuerza repulsiva inversamente proporcional al cubo de la distancia

Para deducir la ecuación de la trayectoria seguiremos los mismos pasos que para una fuerza inversamente proporcional al cuadrado de la distancia (última sección al final de la página)

Supongamos una partícula de masa m que se mueve bajo la acción de una fuerza central atractiva F(r)=k/r3. Las ecuaciones del movimiento son

m( d 2 r d t 2 r ( dθ dt ) 2 )= k r 3 m( r d 2 θ d t 2 +2 dθ dt dr dt )=0

La segunda ecuación nos indica que el momento angular L (expresado en coordenadas polares es constante)

d dt ( m r 2 dθ dt )=0L=m r 2 dθ dt =cte

Teniendo en cuenta la constancia del momento angular L, escribimos la ecuación del movimiento en la dirección radial, de la forma

d 2 r d t 2 1 r 3 ( k m + L 2 m 2 )=0

Para calcular la ecuación de la trayectoria, r=r(θ), eliminamos el tiempo t del siguiente modo

dr dt = dr dθ dθ dt = L m r 2 dr dθ d 2 r d t 2 = d dt ( L m r 2 dr dθ )= d dθ ( L m r 2 dr dθ ) dθ dt = L m r 2 d dθ ( L m r 2 dr dθ ) d 2 r d t 2 = L 2 m 2 r 4 ( d 2 r d θ 2 2 r ( dr dθ ) 2 )

Obtenemos una ecuación en términos de r y de sus derivadas respecto del ángulo θ

d 2 r d θ 2 2 r ( dr dθ ) 2 r( 1+ mk L 2 )=0

Hacemos el cambio de variable u=1/r

dr dθ = 1 u 2 du dθ d 2 r d θ 2 = 2 u 3 ( du dθ ) 2 1 u 2 d 2 u d θ 2

La ecuación del movimiento se convierte en

d 2 u d θ 2 +( 1+ mk L 2 )u=0

La solución de la ecuación diferencial es

u=Asin(βθ)+Bcos(βθ) β= 1+ mk L 2

Como F(r)=k/r3 es una fuerza conservativa, su energía potencial es Ep(r)=k/(2r2). Recuérdese que

F(r)= d dr E p (r)

El momento angular y la energía de la partícula es

{ L=m r 2 dθ dt E= 1 2 m ( dr dt ) 2 + L 2 2m r 2 + k 2 r 2 = 1 2 m ( dr dt ) 2 + L 2 +mk 2m r 2

El segundo término es la energía potencial efectiva, inversamente proporcional al cuadrado de la distancia

fplot(@(x) 1./x.^2, [0.5,4])
E=1.5; %energía 
line([0,4],[E,E],'color','k')
rm=sqrt(1/E); % distancia mínima
line([rm,rm],[0,E],'color','k', 'lineStyle','--')
grid on
xlabel('r')
ylabel('V_{ef}(r)')
title('Energía potencial efectiva')

La partícula se acerca hasta una distancia rm del origen, su componente radial de la velocidad dr/dt=0.

E= L 2 +mk 2m r m 2 r m = L 2 +mk 2mE

La partícula de energía E, parte del infinito r→∞, u=0, su posición angular θ=0. El coeficiente B deberá ser nulo.

La ecuación de la trayectoria es 1/r=Asin(βθ). El coeficiente A se determina sabiendo que para r→∞, la energía de la partícula es E

Teniendo en cuenta que

dr dt = dr dθ dθ dt = β A cos( βθ ) sin 2 ( βθ ) L m r 2 = L m βAcos( βθ )

La energía de la partícula E está relacionada con el coeficiente A

E= 1 2 m ( dr dt ) 2 E= 1 2 L 2 β 2 m A 2

La ecuación de la trayectoria es

1 r = 2mE L 2 +mk sin( 1+ mk L 2 θ )

Dispersión

Una introducción al fenómeno de la dispersión se estudia en la página titulada 'Parámetro de impacto y ángulo de dispersión'

Si v es la velocidad de la partícula cuando r→∞, y b es el parámetro de impacto, la energía E y el momento angular L son

E= 1 2 m v 2 ,L=mvb 1 r =v m m b 2 v 2 +k sin( 1 bv m b 2 v 2 +k m θ )

Dibujamos la trayectoria, tomando el término común bajo la raiz igual a la unidad y tomando una velocidad v=1

b=0.8; %parámero de impacto
th=linspace(10,170,40)*pi*b/180;
r=1./sin(th/b);
plot(r.*cos(th),r.*sin(th))
line([0,1],[0,tan(pi*b/2)], 'lineStyle','--', 'color','k') %vectores
line([0,-1],[0,-tan(pi*b)],'color','k')
grid on
axis equal
xlim([-5,5])
xlabel('x')
ylabel('y')
title('Dispersión')

La velocidad inicial y final de la partícula no cambia, cuando r→∞, cambia su dirección. Se define el ángulo de dispersión Φ el que forma la dirección de la velocidad inicial y final de la partícula, tal como se muestra en la figura

Cuando θ→0 y θ→π/β, entonces r→∞. La posición angular de máximo acercamiento es θ=π /(2β)

1 r = 2mE L 2 +mk sin( βθ ) 1 r m = 2mE L 2 +mk sin( β π 2β ) r m = L 2 +mk 2mE

coincide con la expresión obtenida anteriormente

La relación entre parámetro de impacto b y el ángulo de dispersión Φ es

Φ=π π β ( 1 Φ π ) 2 = m b 2 v 2 m b 2 v 2 +k b 2 = ( πΦ ) 2 m v 2 Φ(2πΦ)

Representamos las trayectorias de las partículas para varios valores del parámetros de impacto b

hold on
for b=0.2:0.2:1 %parámero de impacto
    th=linspace(10,170,100)*pi*b/180;
    r=1./sin(th/b);
    plot(r.*cos(th),r.*sin(th))
end
hold off
grid on
axis equal
xlim([-5,5])
xlabel('x')
ylabel('y')
title('Dispersión')

Fuerza atractiva inversamente proporcional al cubo de la distancia

Como F(r)=-k/r3 es una fuerza conservativa, su energía potencial es Ep(r)=-k/(2r2).

{ L=m r 2 dθ dt E= 1 2 m ( dr dt ) 2 + L 2 2m r 2 k 2 r 2 = 1 2 m ( dr dt ) 2 + L 2 mk 2m r 2

El segundo término es la energía potencial efectiva, inversamente proporcional al cuadrado de la distancia

Se pueden dar tres casos

hold on
fplot(@(x) 1./x.^2, [0.5,4])
fplot(@(x) -1./x.^2, [0.5,4])
line([0,4],[0,0])
E=1.5; %energía 
line([0,4],[E,E],'color','k')
rm=sqrt(1/E); % distancia mínima
line([rm,rm],[0,E],'color','k', 'lineStyle','--')
line([0,4],[-E,-E],'color','k')
line([rm,rm],[0,-E],'color','k', 'lineStyle','--')
hold off
grid on
xlabel('r')
ylabel('V_{ef}(r)')
title('Energía potencial efectiva')

La ecuación del movimientos es la misma que en el apartado anterior sustituyendo k por -k

d 2 u d θ 2 +( 1 mk L 2 )u=0

Examinamos los distintos casos

L2>mk

La solución de la ecuación diferencial es la misma que en el apartado anterior sustituyendo k por -k

u=Asin(βθ)+Bcos(βθ) β= 1 mk L 2

Otra forma alternativa es

u=Ccos( β(θ θ 0 ) )

Calculamos la constante C a partir de la conservación de la energía

La energía de la partícula en coordenadas polares es

E= 1 2 m{ ( dr dt ) 2 + r 2 ( dθ dt ) 2 } k 2 r 2 E= 1 2 L 2 m r 4 { ( dr dθ ) 2 + r 2 } k 2 r 2 E= 1 2 L 2 m u 4 { 1 u 4 ( du dθ ) 2 + 1 u 2 } k 2 u 2 E= 1 2 { L 2 m ( du dθ ) 2 +( L 2 m k ) u 2 }

Conocido u y du/dθ, obtenemos la energía E y despejamos la constante C

C 2 = 2mE β 2 L 2

Para que C sea real, E>0

Cuando r→∞, u→0, β(θ-θ0)=π/2, θ=π/(2β)+θ0

beta=1/7;
th_0=pi/3;
th=th_0:pi/18:th_0+pi/(2*beta)-pi/18;
r=1./cos(beta*(th-th_0));
plot(r.*cos(th),r.*sin(th))
grid on
xlabel('x')
ylabel('y')
title('L^2>mk')

L2<mk

La solución de la ecuación diferencial en u es

u=Aexp(βθ)+Bexp(βθ) β= mk L 2 1

Conocido u y du/dθ, obtenemos la energía E de la partícula

E= 1 2m ( mk L 2 )( 4AB ) AB= mE 2 β 2 L 2

E<0

De forma alternativa

u=Ccoshβ(θ θ 0 )

Conocido u y du/dθ, obtenemos la energía E de la partícula y despejamos la constante C como en el apartado anterior

C= 2mE β 2 L 2

Para que C sea real, E<0. Como la función |cosh(x)|≥1, la trayectoria r=1/u está limitada

beta=1/7;
th_0=pi/6;
th=th_0:pi/18:10*pi;
r=1./cosh(beta*(th-th_0));
plot(r.*cos(th),r.*sin(th))
grid on
axis equal
xlabel('x')
ylabel('y')
title('L^2<mk, E<0')

E>0

De forma alternativa

u=Csinhβ(θ θ 0 )

Conocido u y du/dθ, obtenemos la energía E de la partícula y despejamos la constante C como en el apartado anterior

C= 2mE β 2 L 2

Para que C sea real E>0. r tiende a infinito cuando sinh se hace cero, una trayectoria ilimitada

beta=1/100;
th_0=pi/6;
th=3*pi:-pi/18:th_0-pi/18;
r=1./sinh(beta*(th-th_0));
plot(r.*cos(th),r.*sin(th))
grid on
axis equal
xlabel('x')
ylabel('y')
title('L^2<mk, E>0')

E=0

Cuando E=0, A o B es cero

u=Aexp(βθ) o bien, u=Aexp(-βθ)

beta=1/100;
th_0=pi/6;
th=th_0:pi/36:5*pi;
r=1./exp(beta*th);
plot(r.*cos(th),r.*sin(th))
grid on
axis equal
xlabel('x')
ylabel('y')
title('L^2<mk, E=0')

L2=mk

La ecuación diferencial en u se reduce a

d 2 u d θ 2 =0 u=Aθ+B

E>0

Conocido u y du/dθ, obtenemos la energía E de la partícula y despejamos la constante A como en el apartado anterior

E= 1 2 L 2 m A 2 A= 2mE L 2

th_0=pi/6;
th=th_0:pi/18:5*pi;
r=1./(th+pi);
plot(r.*cos(th),r.*sin(th))
grid on
axis equal
xlabel('x')
ylabel('y')
title('L^2=mk, E>0')

u=0, r=∞, para el ángulo θ=-B/A

E=0

Cuando E=0, el coeficiente A=0, y u=B o r=1/B. La trayectoria es una circunferencia

Referencias

O.L. de Lange, J. Pierrus. Solved Problems in Classical Mechanics. Analytical and numerical solutions with comments. Oxford University Press (2010). Questions 8.11

Douglas Cline. Variational Principles in Classical Mechanics. Second edition, 2019, pp. 273