En física nuclear y de materiales , el poder de frenado es la fuerza retardante que actúa sobre partículas cargadas , típicamente partículas alfa y beta , debido a la interacción con la materia , lo que resulta en la pérdida de energía cinética de la partícula . [1] [2] El poder de frenado también se interpreta como la velocidad a la que un material absorbe la energía cinética de una partícula cargada . Su aplicación es importante en una amplia gama de áreas termodinámicas, como la protección radiológica , la implantación de iones y la medicina nuclear . [3]
Tanto las partículas cargadas como las no cargadas pierden energía al atravesar la materia. Los iones positivos se consideran en la mayoría de los casos a continuación. El poder de frenado depende del tipo y la energía de la radiación y de las propiedades del material por el que pasa. Dado que la producción de un par de iones (normalmente un ion positivo y un electrón (negativo)) requiere una cantidad fija de energía (por ejemplo, 33,97 eV en aire seco [4] : 305 ), el número de ionizaciones por longitud de trayectoria es proporcional al poder de frenado. El poder de frenado del material es numéricamente igual a la pérdida de energía E por unidad de longitud de trayectoria, x :
El signo menos hace que S sea positivo.
La fuerza suele aumentar hacia el final del rango y alcanza un máximo, el pico de Bragg , poco antes de que la energía caiga a cero. La curva que describe la fuerza en función de la profundidad del material se denomina curva de Bragg . Esto tiene una gran importancia práctica para la radioterapia .
La ecuación anterior define el poder de frenado lineal que en el sistema internacional se expresa en N pero que normalmente se indica en otras unidades como M eV /mm o similares. Si se compara una sustancia en forma gaseosa y sólida, entonces los poderes de frenado lineal de los dos estados son muy diferentes simplemente por la diferente densidad. Por lo tanto, a menudo se divide la fuerza por la densidad del material para obtener el poder de frenado másico que en el sistema internacional se expresa en m 4 / s 2 pero que normalmente se encuentra en unidades como MeV/(mg/cm 2 ) o similares. El poder de frenado másico depende entonces muy poco de la densidad del material.
La imagen muestra cómo el poder de frenado de las partículas alfa de 5,49 MeV aumenta a medida que la partícula atraviesa el aire, hasta alcanzar el máximo. Esta energía particular corresponde a la radiación de las partículas alfa del gas radiactivo natural radón ( 222 Rn), que está presente en el aire en cantidades mínimas.
El rango medio se puede calcular integrando la potencia de frenado recíproca sobre la energía: [5]
dónde:
La energía depositada se puede obtener integrando el poder de frenado a lo largo de toda la longitud del recorrido del ion mientras se mueve en el material.
El frenado electrónico se refiere a la desaceleración de un ion proyectil debido a las colisiones inelásticas entre los electrones ligados en el medio y el ion que se mueve a través de él. El término inelástico se utiliza para indicar que se pierde energía durante el proceso (las colisiones pueden dar lugar tanto a excitaciones de los electrones ligados del medio como a excitaciones de la nube de electrones del ion). La potencia de frenado electrónico lineal es idéntica a la transferencia de energía lineal sin restricciones .
En lugar de transferencia de energía, algunos modelos consideran la potencia de frenado electrónica como transferencia de momento entre el gas de electrones y el ion energético. Esto es coherente con el resultado de Bethe en el rango de alta energía. [6]
Dado que el número de colisiones que experimenta un ion con electrones es grande, y dado que el estado de carga del ion mientras atraviesa el medio puede cambiar con frecuencia, es muy difícil describir todas las interacciones posibles para todos los estados de carga de iones posibles. En cambio, la potencia de frenado electrónico a menudo se da como una función simple de la energía que es un promedio tomado de todos los procesos de pérdida de energía para diferentes estados de carga. Se puede determinar teóricamente con una precisión de unos pocos % en el rango de energía por encima de varios cientos de keV por nucleón a partir de tratamientos teóricos, siendo el más conocido la fórmula de Bethe . A energías inferiores a aproximadamente 100 keV por nucleón, se vuelve más difícil determinar el frenado electrónico utilizando modelos analíticos. [7] Recientemente , la teoría funcional de la densidad dependiente del tiempo en tiempo real se ha utilizado con éxito para determinar con precisión el frenado electrónico para varios sistemas ion-objetivo en un amplio rango de energías, incluido el régimen de baja energía. [8] [9]
Paul ha realizado representaciones gráficas de valores experimentales del poder de frenado electrónico para muchos iones en muchas sustancias. [10] La precisión de varias tablas de frenado se ha determinado mediante comparaciones estadísticas. [11]
El poder de frenado nuclear se refiere a las colisiones elásticas entre el ion del proyectil y los átomos en la muestra (la designación establecida "nuclear" puede ser confusa ya que el frenado nuclear no se debe a fuerzas nucleares, [12] pero tiene la intención de señalar que este tipo de frenado implica la interacción del ion con los núcleos en el objetivo). Si uno conoce la forma de la energía potencial repulsiva entre dos átomos (ver a continuación), es posible calcular el poder de frenado nuclear . En la figura de poder de frenado que se muestra arriba para iones de aluminio en aluminio, el frenado nuclear es insignificante excepto en la energía más baja. El frenado nuclear aumenta cuando aumenta la masa del ion. En la figura que se muestra a la derecha, el frenado nuclear es mayor que el frenado electrónico a baja energía. Para iones muy ligeros que se desaceleran en materiales pesados, el frenado nuclear es más débil que el electrónico en todas las energías.
Especialmente en el campo de daños por radiación en detectores, el término " pérdida de energía no ionizante " (NIEL) se utiliza como un término opuesto a la transferencia de energía lineal (LET), véanse, por ejemplo, las referencias. [13] [14] [15] Dado que, por definición, la potencia de frenado nuclear no implica excitaciones electrónicas, la NIEL y el frenado nuclear pueden considerarse la misma cantidad en ausencia de reacciones nucleares.
Por lo tanto, la potencia de frenado total no relativista es la suma de dos términos: . Se han ideado varias fórmulas de potencia de frenado semiempíricas. El modelo propuesto por Ziegler, Biersack y Littmark (el llamado frenado "ZBL", véase el capítulo siguiente), [16] [17] implementado en diferentes versiones de los códigos TRIM/SRIM , [18] es el que se utiliza con más frecuencia en la actualidad.
El poder de frenado radiativo , que se debe a la emisión de radiación de frenado en los campos eléctricos de las partículas en el material atravesado, debe considerarse a energías iónicas extremadamente altas. [12] Para los proyectiles de electrones, el frenado radiativo siempre es importante. A energías iónicas altas, también puede haber pérdidas de energía debido a reacciones nucleares, pero tales procesos normalmente no se describen mediante el poder de frenado. [12]
Cerca de la superficie de un material objetivo sólido, tanto la detención nuclear como la electrónica pueden provocar pulverización catódica .
Al principio del proceso de desaceleración a altas energías, el ion se desacelera principalmente por detención electrónica y se mueve casi en línea recta. Cuando el ion se ha desacelerado lo suficiente, las colisiones con núcleos (detención nuclear) se vuelven cada vez más probables, dominando finalmente la desaceleración. Cuando los átomos del sólido reciben energías de retroceso significativas cuando son golpeados por el ion, se retiran de sus posiciones reticulares y producen una cascada de colisiones adicionales en el material. Estas cascadas de colisiones son la principal causa de producción de daños durante la implantación de iones en metales y semiconductores.
Cuando las energías de todos los átomos del sistema han caído por debajo del umbral de energía de desplazamiento , cesa la producción de nuevos daños y el concepto de detención nuclear ya no tiene sentido. La cantidad total de energía depositada por las colisiones nucleares sobre los átomos de los materiales se denomina energía nuclear depositada.
El recuadro de la figura muestra una distribución típica de los iones depositados en el sólido. El caso que se muestra aquí podría ser, por ejemplo, la desaceleración de un ion de silicio de 1 MeV en silicio. El rango medio para un ion de 1 MeV está típicamente en el rango micrométrico .
A distancias muy pequeñas entre los núcleos, la interacción repulsiva puede considerarse esencialmente coulombiana. A distancias mayores, las nubes de electrones apantallan los núcleos entre sí. Por lo tanto, el potencial repulsivo puede describirse multiplicando la repulsión coulombiana entre núcleos por una función de apantallamiento φ(r/a),
donde φ(r/a) → 1 cuando r → 0. Aquí y son las cargas de los núcleos que interactúan, y r la distancia entre ellos; a es el llamado parámetro de apantallamiento.
A lo largo de los años se han propuesto una gran cantidad de potenciales repulsivos y funciones de apantallamiento diferentes, algunos determinados de forma semiempírica y otros a partir de cálculos teóricos. Un potencial repulsivo muy utilizado es el propuesto por Ziegler, Biersack y Littmark, el denominado potencial repulsivo ZBL. Se ha construido ajustando una función de apantallamiento universal a potenciales obtenidos teóricamente calculados para una gran variedad de pares de átomos. [16] El parámetro y la función de apantallamiento ZBL tienen las formas
y
donde x = r/a u , y a 0 es el radio atómico de Bohr = 0,529 Å.
La desviación estándar del ajuste del potencial repulsivo universal ZBL a los potenciales específicos de pares calculados teóricamente a los que se ajusta es del 18 % por encima de 2 eV. [16] Se pueden obtener potenciales repulsivos aún más precisos a partir de cálculos de energía total autoconsistentes utilizando la teoría funcional de la densidad y la aproximación de densidad local (LDA) para el intercambio y correlación electrónicos. [19]
En los materiales cristalinos, en algunos casos, el ion puede "canalizarse", es decir, enfocarse en un canal entre los planos cristalinos donde casi no experimenta colisiones con los núcleos. [20] Además, el poder de frenado electrónico puede ser más débil en el canal. [21] Por lo tanto, el frenado nuclear y electrónico no solo dependen del tipo y la densidad del material, sino también de su estructura microscópica y sección transversal.
Los métodos de simulación por ordenador para calcular el movimiento de iones en un medio se han desarrollado desde la década de 1960 y ahora son la forma dominante de tratar teóricamente el poder de frenado. La idea básica de estos métodos es seguir el movimiento del ion en el medio simulando las colisiones con los núcleos del medio. El poder de frenado electrónico suele tenerse en cuenta como una fuerza de fricción que frena el ion.
Los métodos convencionales utilizados para calcular los rangos de iones se basan en la aproximación de colisión binaria (BCA). [22] En estos métodos, el movimiento de iones en la muestra implantada se trata como una sucesión de colisiones individuales entre el ion de retroceso y los átomos en la muestra. Para cada colisión individual, la integral de dispersión clásica se resuelve mediante integración numérica.
El parámetro de impacto p en la integral de dispersión se determina a partir de una distribución estocástica o de una manera que tenga en cuenta la estructura cristalina de la muestra. El primer método es adecuado únicamente en simulaciones de implantación en materiales amorfos, ya que no tiene en cuenta la canalización.
El programa de simulación BCA más conocido es TRIM/SRIM ( acrónimo de TRansport of Ions in Matter, en versiones más recientes llamado Stopping and Range of Ions in Matter), que se basa en el frenado electrónico ZBL y el potencial interatómico . [16] [18] [23] Tiene una interfaz de usuario muy fácil de usar, y tiene parámetros predeterminados para todos los iones en todos los materiales hasta una energía iónica de 1 GeV, lo que lo ha hecho inmensamente popular. Sin embargo, no tiene en cuenta la estructura cristalina, lo que limita severamente su utilidad en muchos casos. Varios programas BCA superan esta dificultad; algunos bastante conocidos son MARLOWE, [24] BCCRYS y crystal-TRIM.
Aunque los métodos BCA se han utilizado con éxito para describir muchos procesos físicos, presentan algunos obstáculos para describir el proceso de desaceleración de iones energéticos de forma realista. La suposición básica de que las colisiones son binarias da lugar a graves problemas cuando se intenta tener en cuenta múltiples interacciones. Además, al simular materiales cristalinos, el proceso de selección del siguiente átomo reticular que colisiona y el parámetro de impacto p siempre implican varios parámetros que pueden no tener valores perfectamente bien definidos, lo que puede afectar los resultados en un 10-20% incluso para elecciones aparentemente bastante razonables de los valores de los parámetros. La mejor fiabilidad en BCA se obtiene incluyendo múltiples colisiones en los cálculos, lo que no es fácil de hacer correctamente. Sin embargo, al menos MARLOWE lo hace.
Una forma fundamentalmente más sencilla de modelar colisiones atómicas múltiples es la que proporcionan las simulaciones de dinámica molecular (MD), en las que se calcula la evolución temporal de un sistema de átomos resolviendo numéricamente las ecuaciones de movimiento. Se han ideado métodos MD especiales en los que se ha reducido el número de interacciones y átomos implicados en las simulaciones MD para que sean lo suficientemente eficientes para calcular rangos de iones. [25] [26] Las simulaciones MD describen automáticamente el poder de frenado nuclear. El poder de frenado electrónico se puede incluir fácilmente en las simulaciones de dinámica molecular, ya sea como una fuerza de fricción [25] [27] [ 28] [29] [26] [30] [31] [32] o de una manera más avanzada siguiendo también el calentamiento de los sistemas electrónicos y acoplando los grados de libertad electrónicos y atómicos. [33] [34] [35]
Más allá del máximo, el poder de frenado disminuye aproximadamente como 1/v 2 con el aumento de la velocidad de la partícula v , pero después de un mínimo, aumenta de nuevo. [36] Una partícula ionizante mínima (MIP) es una partícula cuya tasa media de pérdida de energía a través de la materia está cerca del mínimo. En muchos casos prácticos, las partículas relativistas (por ejemplo, los muones de rayos cósmicos ) son partículas ionizantes mínimas. Una propiedad importante de todas las partículas ionizantes mínimas es que es aproximadamente verdadera donde y son las cantidades cinemáticas relativistas habituales. Además, todas las MIP tienen casi la misma pérdida de energía en el material cuyo valor es: . [36]
{{cite book}}
: CS1 maint: location missing publisher (link)