Ciclotrón

Laura Alejandra Ayala Cante y Santiago Steven Beltrán Romero

Universidad de los Andes.
Departamento de Fı́sica.
Aceleradores de partı́culas

(Dated: 14 de diciembre de 2020)

I. INTRODUCCIÓN

Los aceleradores de partı́culas son dispositivos que tie-

nen un espectro amplio de aplicaciones en cuanto a la
investigación cientı́fica. Por esta razón, es de gran im-
portancia estudiar con más detalle el mecanismo fı́sico
y todas las variables que influyen en su funcionamiento.
Particularmente, este trabajo se centrará en estudiar el
comportamiento de una partı́cula cargada positivamen-
te en un ciclotrón. Para tal fin, se profundizará en los
parámetros necesarios para modelar computacionalmen-
te la trayectoria que sigue un protón dentro de un ci-
clotrón teniendo en cuenta las fuerzas que actúan sobre
ella al estar dentro del acelerador. Con el fin de hacer una Figura 1. Funcionamiento general del ciclotrón. [1]
caracterización completa del comportamiento del protón,
además de obtener una gráfica de la posición del protón,
se analizarán gráficas de velocidad, tiempo, energı́a y mo- ca que se escribe como:
mento. Ası́ mismo, se hará una descripción de los paráme-
tros de entrada los cuales serán el determinante para ob- Fm = qvB (1)
tener resultados numéricos precisos.
Aplicando la segunda ley de Newton y que la aceleración
2
centrı́peta es ac = vR , se puede despejar el radio de la
II. MARCO TEÓRICO
trayectoria y la velocidad de la partı́cula

La base del funcionamiento del ciclotrón consiste en mv

R=
dos placas semicirculares huecas que se ubican en un cam- qB
(2)
po magnético paralelo al eje de las placas el cual se en- qBR
v=
carga de curvar la trayectoria del protón. Entre las dos m
placas, con un sistema de radiofrecuencia se genera un
campo eléctrico oscilante con frecuencia de resonancia, De allı́ se puede ver que la velocidad y el radio de la
de esta forma la partı́cula siempre estará en fase con el trayectoria aumentan linealmente el uno con el otro. Te-
campo E ~ y será acelerada. Finalmente, en un extremo de niendo la expresión de la velocidad se puede hallar la
la zona intermedia de las placas, sale el rayo final en don- energı́a cinética de la partı́cula que está dada por:
de está ubicado el detector. Como consecuencia, el campo
eléctrico acelera el protón hacia el interior de una de las q 2 B 2 R2
K= (3)
placas semicirculares donde realizan una trayectoria cir- 2m
cular debido al campo magnético que los lleva de nuevo Esta ecuación implica que al aumentar la energı́a, aumen-
a la región entre las placas. De esta forma se obtiene una ta el radio de la trayectoria y de esta forma se obtiene
trayectoria circular como lo muestra la Figura 1. El cam- la forma de espiral. Es decir, los arcos con un radio más
po magnético debe ser tal que el tiempo que demore la pequeño tienen menor energı́a que los arcos que descri-
partı́cula en recorrer un semicı́rculo sea igual a la mitad be la partı́cula con un radio mayor dado que la fuente
del periodo de las oscilaciones del campo eléctrico. Ası́, de partı́culas está en el centro y estas empiezan a ganar
cuando el protón llegue a la región intermedia el campo energı́a. Por lo tanto se debe fijar la magnitud del campo
eléctrico se habrá invertido de la forma que la partı́cula magnético en
se acelerada hacia la otra placa.
Para encontrar las relaciones correspondientes a la ci- 1 √
nemática del movimiento, partimos de la fuerza magnéti- B= 2mK (4)
qR
 2

Adicionalmente, la energı́a cinética que tendrá el protón y el momento en ambas variables.También debido al he-
puede expresarse como el producto de la diferencia de cho de que el campo magnético perpendicular (en el eje
potencial existente entre las D’s y el número de veces z) B acopla las ecuaciones diferenciales que describen la
que los iones son acelerados por ésta. Sin embargo, el cinemática de la partı́cula.
tiempo que demora la partı́cula en recorrer cada arco es
Inicialmente, se definen los campos magnéticos y
una constante; entonces el tiempo se puede expresar en
eléctricos que se van a usar. Para el campo magnético
términos de la frecuencia del ciclotrón como:
B se emplea la ecuación 4 con lo cual el protón adquiere
la energı́a necesaria al momento de salir del Ciclotrón.
2πR 2πm π Por otro lado, el campo eléctrico E se aplica cuando la
T = = = . (5) partı́cula tiene una posición en y entre − d2 y d2 , con di-
v qB ωcicl
rección en el eje y y viene dado por la ecuación 10. En
Debido a la condición de resonancia de fase en la que ella se va a considerar una magnitud 107 veces mayor que
se debe encontrar la partı́cula, se obtiene la relación pa- el campo magnético para tener un número grande pero
ra la frecuencia del ciclotrón que debe ser la misma del perceptible de vueltas y se tomará como fase φ = 0 con
oscilador el fin de que se aplique en esa región, los valores máxi-
mos del campo eléctrico para acelerar el protón. Por esta
razón, se define la función para este campo dada por
v qB
ωosc = ωcicl = = (6)
R m

E0 cos (ωosc t), si − d2 ≤ y ≤ d2 ,
Debido al principio del funcionamiento del ciclotrón E(y, t) = , (10)
0 de lo contrario.
(campo magnético y eléctrico) se observa una propiedad
de enfoque vertical en la trayectoria de la partı́cula. Esto
quiere decir que gracias al campo eléctrico las partı́culas
la cual será de utilidad en la definición de la aceleración
se enfocan en el centro corrigiendo la desviación con la
en y.
que emergen de la fuente. Por esta razón, es más eficien-
te para partı́culas de baja energı́a dada la fuerza de las
lı́neas de campo eléctrico. También hay un enfoque débil
debido al cambio en el campo magnético en la compo-
nente z.
A partir de estas relaciones que modelan el comporta- 1.00
miento de la partı́cula, se puede escribir la ecuación que
modela el movimiento del protón entre la D de la forma: 0.75

ma = qv × Bo ẑ (7) 0.50

Asimismo, se obtiene la siguiente ecuación para la región 0.25

Posición en y (m)

entre las placas:

0.00
ma = qv × Bo ẑ + qE (8)
−0.25

III. MARCO COMPUTACIONAL −0.50

Con el fin de comprender el comportamiento de una −0.75

partı́cula cargada que se acelera en un Ciclotrón se reali-

−1.00
zará un modelo computacional del paso a paso que surge
−1.00 −0.75 −0.50 −0.25 0.00 0.25 0.50 0.75 1.00
dicha partı́cula. De esta manera se empleará el método Posición en x (m)
explicado por Feynmann [2] en el que la evolución de po-
sición y velocidad de la partı́cula se da a primer orden.
Esto quiere decir que en el caso de la velocidad para una
Figura 2. Trayectoria del protón en un Ciclotrón desde el
aceleración a dada en un determinado tiempo se cumple
origen coordenado hasta obtener un radio R = 1m con
que campo magnético B = 2, 044ẐT y la aplicación del cam-
po eléctrico dada por la ecuación 10 de valor máximo E0 =
dv vf − vo
a= ≈ . (9) 2, 044 × 107 V m.
dt dt
Este proceso se observa en dimensión dos x y y por lo
que simultáneamente se debe ir actualizando la posición El paso a paso se inicia con la partı́cula con las siguien-
 3

1e−19 1e−19

3 3

2 2

1 1
Momento en x (kg m/s)

Momento en y (kg m/s)

0 0

−1 −1

−2 −2

−3 −3

−1.00 −0.75 −0.50 −0.25 0.00 0.25 0.50 0.75 1.00 −1.00 −0.75 −0.50 −0.25 0.00 0.25 0.50 0.75 1.00
Posición en x (m) Posición en y (m)

Figura 3. Diagramas de fase en x (izquierda) y y (derecha) para la trayectoria de un protón en un Ciclotrón.

tes posiciones, velocidades y aceleraciones iniciales Este proceso se repite cambiando t por t + dt, lo cual
genera la evolución deseada al repetir sucesivamente es-
x(0) = 0.0m y y(0) = 0.0m, te algoritmo. En cada repetición es necesario calcular el
dt m dt m radio como
vx ( ) = 0.0 y vy ( ) = 0.0 , (11)
2 s 2 s
p
m m r(t) = x(t)2 + y(t)2 (15)
ax (0) = 0.0 2 y ay (0) = 0.0 2
s s para establecer el punto en que se finalizan las iteraciones,
T 2π
esta finalización se da cuando r(t) ≥ R con R el radio
y se evoluciona temporalmente con dt = = ωosc N para N máximo del Ciclotrón.
N arbitrario y grande, de la siguiente manera. Primero En lo que respecta al análisis experimental de la itera-
se calcula la aceleración con los datos iniciales, de allı́ se ción en la trayectoria del protón se tiene como salida las
obtiene en general que posiciones x y y, las velocidades vx y vy , y el tiempo t.
De esta forma, se representa con las posiciones, la trayec-
dt toria de la partı́cula parametrizada por el tiempo t en el
ax (t + dt) =ωosc vy (t + )
2 (12) que se desarrolla el movimiento. Similarmente, es posible
dt q representar los diagramas de fase en x y y, calculando los
ay (t + dt) = − ωosc vx (t + ) + E(y, t).
2 m momentos px y py a partir de las respectivas velocidades
y la masa del protón. De este análisis se espera obtener
Luego, con los datos previos de velocidad se hallan las el comportamiento en espiral caracterı́stico del ciclotrón,
nuevas posiciones como sigue donde se observe la influencia del campo eléctrico en la
aceleración de la carga. Posteriormente, se realizó del ra-
dt
x(t + dt) =x(t) + vx (t + ) dio de la trayectoria cuando x es cercano a cero o está en
2 (13) el rango − ≤ x ≤ , con  = 7 × 10−4 m. En este rango
dt se calcularon distintos radios a través de dos métodos:
y(t + dt) =y(t) + vy (t + ).
2 el primero según la ecuación 15 en el que se emplean las
posiciones, y el segundo a partir de la ecuación
Para después determinar ambas velocidades, a partir de
las aceleraciones ya calculadas ası́ 1 q 2
rp = px + p2y , (16)
qB
3 dt
vx (t + dt) = vx (t + ) + ax (t + dt) la cual hace uso de los momentos en los puntos escogidos.
2 2 (14)
3 dt A partir de esta comparación se espera obtener que r =
vy (t + dt) = vy (t + ) + ay (t + dt). rp en los distintos puntos de la trayectoria del protón.
2 2
 4

Por otro lado, también se revisa el comportamiento

del periodo para algunas de las vueltas descritas por el
protón en el acelerador de Ciclotrón. Para ello, se con-
Fit lineal
sideran los puntos con x > 0 y − ≤ y ≤  en especial
0.7
los tiempos de dichos puntos. Seguidamente, se restan los
tiempos de puntos consecutivos y se ignoran las diferen-
0.6
cias de tiempos T < 1.0 × 10−10 s, ya que representan
puntos en una misma vuelta. Respecto a las vueltas ele-

Radio según el momento (m)

0.5
gidas se comparan los periodos consecutivos con el valor
del periodo del Ciclotrón T de la ecuación 5. Un último
0.4
análisis se realiza acerca del incremento de la energı́a del
protón. Esta energı́a se determina a partir de
0.3
1 1 2
E = m(vx2 + vy2 ) = (p + p2y ) (17)
2 2m x 0.2

y se espera que la partı́cula alcance sucesivamente la 0.1

energı́a máxima fijada, donde en cada vuelta la energı́a
sea constante. Esta serie de conclusiones esperables serán 0.0
estudiadas a partir del modelo planteado para el Ci-
0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7
clotrón y sus resultados se encuentran condensados y ana- Radio según la posición (m)
lizados en la siguiente Sección.

Figura 4. Gráfica del radio obtenido por medio del momento

IV. RESULTADOS Y ANÁLISIS
de la partı́cula en función del radio obtenido a través de la
posición, de un protón en un Ciclotrón.
Para obtener las gráficas que permitan caracterizar
completamente la trayectoria del protón dentro del ci-
clotrón, se usaron los siguientes datos de entrada: 2, especialmente en el diagrama en x. En efecto, este
comportamiento se debe a que la partı́cula describe un
R = 1m comportamiento circular y con cada vuelta aumenta su
velocidad. Por lo tanto es clara la influencia de ambos
E = 200M eV campos eléctrico y magnético. En particular, la influen-
cia del campo eléctrico es más evidente en la distorsión
q = 1.602 × 10−19 C
inicial del diagrama de fase de y para velocidades bajas,
m = 1, 673 × 10−27 kg ya que el protón tenderá inicialmente a tener una tra-
yectoria rectilı́nea y a medida que aumenta que aumenta
B = 2, 044ẑT la velocidad la influencia del campo magnético es ma-
yor y en tanto el comportamiento circular pierde dicha
E0 = 2, 044 × 107 V m distorsión. Adicionalmente, se observó, como es intuiti-
d = R/10 vo notarlo, que la magnitud del campo eléctrico define el
número de vueltas efectuadas por la partı́cula cargada.
La Figura 2 muestra la gráfica de posición en x vs Respecto a la comparación de radios obtenidos a partir
la posición en y del protón. Se puede observar que la del momento y de la posición para algunos puntos en par-
partı́cula parte del reposo y empieza a acelerarse debido ticular, se obtiene la gráfica de la Figura 4 la cual muestra
al campo eléctrico y a describir trayectorias circulares un comportamiento lineal entre los valores para el radio
debido al campo magnético. En esta gráfica se observa el según el momento y según la posición. La ecuación de la
espiral que describe la partı́cula y el aumento gradual el recta es:
radio de la misma. Al final de la trayectoria, al alcanzar el
radio máximo, se observa que las circunferencias son más
rp = 0.983r + 0.00894m (18)
pegadas que al inicio; sin embargo, el tiempo de cada una
es el mismo. En las siguientes gráficas se podrá analizar A partir de este resultado se comprueba que la rela-
más a fondo como cambian estos factores como el tiempo ción es lineal y por tanto los valores para el radio dan
y el radio de la trayectoria. Ası́, de la gráfica de posición resultados similares cumpliendo la igualdad r = rp . Ası́
se puede verificar el comportamiento circular del protón mismo, al utilizar dos métodos diferentes para hallar el
y la influencia del campo magnético y eléctrico. radio, comprobamos que la trayectoria que describe el
Similarmente, se representan los diagramas de fase en protón si es completamente circular como se esperaba.
x y y en la Figura 3. En ellos se observa un comporta- Por otra parte, del mecanismo del Ciclotrón se desea
miento similar de espiral a la trayectoria de la Figura que en cada una de sus vueltas tarde un mismo tiempo
 5

1e−8 1e−11

Periodo teórico T
6 Datos Ciclotrón 3.0

2.5
5
2.0
Periodo (s)

Energía (J)
4
1.5

3 1.0

0.5
2
0.0

0 2 4 6 8 10 0 2 4
Tiempo (s)
6 8
1e−7
Vueltas

Figura 5. Gráfica del periodo de distintas vueltas realizadas Figura 6. Energı́a de la partı́cula en función del tiempo du-
por el protón en el Ciclotrón. rante su movimiento al interior del Ciclotrón.

partı́cula dentro del ciclotrón. Se observa que la energı́a

en llegar a un mismo punto. En particular, por el he- aumenta conforme pasa el tiempo y que en cada recorrido
cho de considerar el campo magnético perpendicular y de las D’s es constante. Ası́ mismo, se observan pequeños
constante, aunado a las caracterı́sticas del protón se ob- saltos en los valores de la energı́a los cuales correspon-
tuvo que ese tiempo entre vueltas debe ser teóricamente den al cambio de orbita que describe la partı́cula en su
Tteo = 3, 2105 × 10−8 s gracias a usar la ecuación 5. En trayectoria. Es decir, este comportamiento describe cla-
este sentido, en el desarrollo computacional se intentaron ramente la energı́a del protón al pasar a una orbita más
tomar la mayor cantidad de vueltas consecutivas con el grande debido a la energı́a que se transfiere por el campo
obstáculo de que definir los puntos a considerar como un eléctrico presente en la región intermedia. En la gráfi-
rango lleva a tomar puntos en una misma vuelta. Por lo ca también se observa que efectivamente la energı́a llega
tanto, como se muestra en la Figura 5, estos puntos de al valor esperado de energı́a final, empero se puede ver
una misma vuelta representan diferencias de tiempos que que sobrepasa mı́nimamente el lı́mite. Los siguientes va-
son descartadas al ser menores a 1.0 × 10−10 s. En conse- lores muestran la diferencia entre el valor de energı́a final
cuencia, se estudiaron 11 vueltas, un número menor que teórico y el que se obtuvo en la simulación.
las realizadas por el protón, pero que muestran la ten-
dencia de la mayorı́a de ellas a tener un periodo cercano Eteo = 3, 204 × 10−11 J
al valor teórico. (19)
Eexp = 3, 232 × 10−11 J,
Como consecuencia de analizar la discrepancia con el
valor real se obtuvo una desviación de 1, 095 × 10−8 s. lo cual nos da un error del 0, 87 %. Este mı́nimo error se
Este valor se debe en gran medida a la diferencia gene- debe a que la trayectoria no está completamente centrada
rada por la vuelta 4, la cual muy probablemente se trate en el origen ya que en este punto ya hay una velocidad
de tres vueltas consecutivas. Por lo tanto, este resulta- inicial en la partı́cula.
do hace referencia en mayor medida a la efectividad del
método de elección de un punto por cada vuelta que de
la trayectoria en sı́, ya que en bastantes vueltas los datos V. CONCLUSIONES
coincidieron con el valor real. De esta manera, se hace
necesario proponer algoritmos basados, por ejemplo, en En conclusión, el Ciclotrón es un acelerador de partı́cu-
el cambio de signo en la trayectoria guarde los datos, pa- las que emplea campos magnéticos para describir las tra-
ra que se consideren más vueltas y el cálculo del periodo yectorias circulares y simultáneamente contribuir al enfo-
tenga una menor desviación. que, a medida que un campo eléctrico aplicado en medio
El resultado de la Figura 6 muestra el comportamien- de las D’s ayuda al aumento de la energı́a de las partı́cu-
to de la energı́a del protón respecto al tiempo que dura la las. La aplicación de este campo eléctrico a una frecuencia
 6

cercana a la frecuencia de ciclotrón es lo que contribuye a de 1, 095 × 10−8 s con respecto al valor teórico. En este
que por resonancia la partı́cula cargada se acelere y des- proceso las discrepancias con respecto a la teorı́a se deben
criba trayectorias circulares de mayor radio. Mediante el a las dificultades del método para seleccionar un punto
proceso computacional de Feynmann se pudo modelar las por cada vuelta consecutiva, de manera que con estas
espirales caracterı́sticas en la trayectoria y diagramas de caracterı́sticas se identificaron 12 vueltas, una cantidad
fase de las partı́culas en un Ciclotrón, donde se represen- menor a las vueltas totales. En cuanto a la energı́a se
tan los efectos de ambos campos. Además, se estableció la observó el aumentó progresivo de la misma a medida que
igualdad en los radios calculados empleando la posición se pasaba a una orbita de radio mayor. Después de cier-
y el momento, lo cual refleja la relación de estas variables to tiempo, se notó que el Ciclotrón llegó a una energı́a
conjugadas en el movimiento. bastante próxima a la deseada con un error de 0, 87 %
Asimismo, se comprobó la caracterı́stica de periodos causado por el hecho de que la trayectoria no está del
iguales en las vueltas del Ciclotrón con una desviación todo centrada en el origen.

[1] Sc.ehu.es, “ Acelerador de partı́culas cargadas. El ci- [2] R. F. M. Sands and R. Leighton, The Feynman Lectures
clotrón,” (2007), [Web; accedido el 13-12-2020]. on physics, Vol. 2 (Addison–Wesley, 1964).