stringtranslate.com

Poder de frenado (radiación de partículas)

En física nuclear y de materiales , el poder de frenado es la fuerza retardadora 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 las partículas . [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]

Gráfico que muestra las relaciones entre la radiactividad y la radiación ionizante detectada.

Definición y curva de Bragg

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 iónico (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 camino 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 del recorrido, x :

El signo menos hace que S sea positivo.

Curva de Bragg de partículas alfa de 5,49 MeV en el aire

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 llama curva de Bragg . Esto es de gran importancia práctica para la radioterapia .

La ecuación anterior define la potencia de frenado lineal que en el sistema internacional se expresa en N pero 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 detención lineales de los dos estados son muy diferentes simplemente debido a la diferente densidad. Por lo tanto, a menudo se divide la fuerza por la densidad del material para obtener la potencia de frenado de la masa 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. La capacidad de frenado de la masa 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 partículas alfa del gas naturalmente radiactivo radón ( 222 Rn), que está presente en el aire en cantidades mínimas.

El alcance medio se puede calcular integrando la potencia de frenado recíproca sobre la energía: [5]

dónde:

E 0 es la energía cinética inicialde la partícula.
Δ x es el rango de "aproximación de desaceleración continua (CSDA)" y
S ( E ) es la potencia de frenado lineal.

La energía depositada se puede obtener integrando el poder de frenado en todo el recorrido del ion mientras se mueve en el material.

Parada electrónica, nuclear y radiativa.

La detención electrónica 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 resultar tanto en excitaciones de los electrones ligados al medio como también en excitaciones de la nube de electrones del ion). La potencia de frenado electrónica lineal es idéntica a la transferencia de energía lineal sin restricciones .

En lugar de la transferencia de energía, algunos modelos consideran el poder de frenado electrónico como una transferencia de impulso entre el gas de electrones y el ion energético. Esto es consistente con el resultado de Bethe en el rango de alta energía. [6]

Dado que el número de colisiones que un ion experimenta con los 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 posibles de los iones. En cambio, la potencia de frenado electrónica a menudo se expresa como una función simple de la energía, que es un promedio de todos los procesos de pérdida de energía para diferentes estados de carga. Puede determinarse 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, resulta más difícil determinar la parada electrónica 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 la detención electrónica de varios sistemas de objetivos iónicos en una amplia gama de energías, incluido el régimen de baja energía. [8] [9]

Poder de frenado electrónico y nuclear para iones de aluminio en aluminio, versus energía de partículas por nucleón. El máximo de la curva de parada nuclear ocurre típicamente a energías del orden de 1 keV por nucleón .

Paul ha realizado presentaciones gráficas de valores experimentales del poder de parada electrónico de muchos iones en muchas sustancias. [10] La precisión de varias tablas de parada se ha determinado mediante comparaciones estadísticas. [11]

El poder de detención nuclear se refiere a las colisiones elásticas entre el ion del proyectil y los átomos de la muestra (la designación establecida "nuclear" puede resultar confusa ya que la detención nuclear no se debe a fuerzas nucleares, [12] pero se pretende señalar que este tipo de la parada implica la interacción del ion con los núcleos en el objetivo). Si se conoce la forma de la energía potencial repulsiva entre dos átomos (ver más abajo), es posible calcular el poder de frenado nuclear . En la figura del poder de frenado que se muestra arriba para los iones de aluminio en aluminio, el frenado nuclear es insignificante excepto en la energía más baja. La parada nuclear aumenta cuando aumenta la masa del ion. En la figura que se muestra a la derecha, la parada nuclear es mayor que la parada electrónica a baja energía. Para iones muy ligeros que se desaceleran en materiales pesados, la detención nuclear es más débil que la electrónica en todas las energías.

Especialmente en el campo de los daños por radiación en detectores, el término " pérdida de energía no ionizante " (NIEL) se utiliza como término opuesto a la transferencia lineal de energía (LET), véanse, por ejemplo, las Refs. [13] [14] [15] Dado que, por definición, la potencia de parada nuclear no implica excitaciones electrónicas, se puede considerar que NIEL y la parada nuclear son la misma cantidad en ausencia de reacciones nucleares.

Por tanto, el poder de frenado total no relativista es la suma de dos términos: . Se han ideado varias fórmulas semiempíricas de poder de frenado. El modelo dado por Ziegler, Biersack y Littmark (la llamada parada "ZBL", ver el próximo capítulo), [16] [17] implementado en diferentes versiones de los códigos TRIM/SRIM , [18] se usa con mayor frecuencia en la actualidad.

En el caso de energías iónicas extremadamente altas [3] , también hay que tener en cuenta el poder de frenado de la radiación , que se debe a la emisión de bremsstrahlung en los campos eléctricos de las partículas del material atravesado. [12] Para los proyectiles de electrones, la detención radiativa siempre es importante. A altas energías de iones, también puede haber pérdidas de energía debido a reacciones nucleares, pero estos procesos normalmente no se describen mediante energía de parada. [12]

Cerca de la superficie de un material objetivo sólido, la detención tanto nuclear como electrónica puede provocar chisporroteo .

El proceso de desaceleración en sólidos.

Ilustración de la desaceleración de un solo ion en un material sólido.

Al comienzo del proceso de desaceleración a altas energías, el ion se desacelera principalmente mediante detención electrónica y se mueve casi en línea recta. Cuando el ion se ha desacelerado lo suficiente, las colisiones con los núcleos (la parada 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, serán removidos de sus posiciones en la red y producirán 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 parada nuclear ya no tiene sentido. La cantidad total de energía depositada por las colisiones nucleares en los átomos de los materiales se denomina energía nuclear depositada.

El recuadro de la figura muestra una distribución típica de 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 suele estar en el rango micrométrico .

Potenciales interatómicos repulsivos

A distancias muy pequeñas entre los núcleos, la interacción repulsiva puede considerarse esencialmente coulómbica. A mayores distancias, las nubes de electrones protegen los núcleos entre sí. Por tanto, el potencial de repulsión se puede describir multiplicando la repulsión de Coulomb entre núcleos con una función de detección φ(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 detección.

A lo largo de los años se han propuesto una gran cantidad de potenciales repulsivos y funciones de detección 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 llamado potencial repulsivo ZBL. Se ha construido ajustando una función de detección universal a potenciales obtenidos teóricamente y calculados para una gran variedad de pares de átomos. [16] El parámetro y la función de detección ZBL tienen la forma

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 del par 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 autoconsistentes de energía total utilizando la teoría del funcional de densidad y la aproximación de densidad local (LDA) para el intercambio y la correlación electrónicos. [19]

canalización

En materiales cristalinos, el ion puede en algunos casos "canalizarse", es decir, concentrarse en un canal entre planos cristalinos donde casi no experimenta colisiones con los núcleos. Además, la potencia de frenado electrónica puede ser más débil en el canal. Por tanto, la parada nuclear y electrónica no sólo depende del tipo y densidad del material sino también de su estructura microscópica y sección transversal.

Simulaciones por computadora de iones desacelerando

Los métodos de simulación por computadora para calcular el movimiento de los 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 en ellos es seguir el movimiento del ion en el medio simulando las colisiones con los núcleos en el medio. El poder de frenado electrónico generalmente se tiene 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). [20] 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 de 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 sólo es adecuado 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 potencial interatómico y de parada electrónica ZBL . [16] [18] [21] 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 gravemente su utilidad en muchos casos. Varios programas BCA superan esta dificultad; algunos bastante conocidos son MARLOWE, [22] BCCRYS y crystal-TRIM.

Aunque los métodos BCA se han utilizado con éxito para describir muchos procesos físicos, tienen algunos obstáculos para describir de manera realista el proceso de desaceleración de los iones energéticos. La suposición básica de que las colisiones son binarias genera graves problemas al intentar tener en cuenta múltiples interacciones. Además, al simular materiales cristalinos, el proceso de selección del siguiente átomo de la red que colisiona y el parámetro de impacto p siempre involucran varios parámetros que pueden no tener valores perfectamente definidos, lo que puede afectar los resultados entre un 10% y un 20% incluso para elecciones que parecen bastante razonables. los valores de los parámetros. La mejor confiabilidad en BCA se obtiene incluyendo múltiples colisiones en los cálculos, lo cual no es fácil de hacer correctamente. Sin embargo, al menos MARLOWE hace esto.

Una forma fundamentalmente más sencilla de modelar múltiples colisiones atómicas la proporcionan las simulaciones de dinámica molecular (MD), en las que la evolución temporal de un sistema de átomos se calcula 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 involucrados en las simulaciones MD para hacerlos lo suficientemente eficientes para calcular rangos de iones. [23] [24] Las simulaciones MD describen automáticamente la potencia de frenado nuclear. El poder de frenado electrónico se puede incluir fácilmente en simulaciones de dinámica molecular, ya sea como una fuerza de fricción [23] [25] [26] [27] [24] [28] [29] [30] o de una manera más avanzada también siguiendo el calentamiento de los sistemas electrónicos y acoplando los grados de libertad electrónicos y atómicos. [31] [32] [33]

Partícula ionizante mínima

Más allá del máximo, el poder de frenado disminuye aproximadamente como 1/v 2 al aumentar la velocidad de la partícula v , pero después de un mínimo, aumenta nuevamente. [34] 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 es cercana al mínimo. En muchos casos prácticos, las partículas relativistas (por ejemplo, los muones de rayos cósmicos ) son partículas mínimamente ionizantes. Una propiedad importante de todas las partículas ionizantes mínimas es que es aproximadamente cierta donde y son las cantidades cinemáticas relativistas habituales. Además, todos los MIP tienen casi la misma pérdida de energía en el material cuyo valor es: . [34]

Ver también

Referencias

  1. ^ Bragg, WH (1905). "Sobre las partículas α de radio y su pérdida de alcance al pasar a través de varios átomos y moléculas". Fil. Mag . 10 (57): 318. doi : 10.1080/14786440509463378.
  2. ^ Bohr, N. (1913). "Sobre la teoría de la disminución de la velocidad de las partículas electrificadas en movimiento al atravesar la materia". Fil. Mag . 25 (145): 10. doi :10.1080/14786440108634305.
  3. ^ ab Informe ICRU 73: Detención de iones más pesados ​​que el helio, Journal of the ICRU, 5 No. 1 (2005), Oxford Univ. Presione ISBN 0-19-857012-0 
  4. ^ Podgorsak, EB, ed. (2005). Física de oncología radioterápica: manual para profesores y estudiantes (PDF) . Viena: Organismo Internacional de Energía Atómica. ISBN 978-92-0-107304-4. Consultado el 25 de noviembre de 2012 .
  5. ^ Comisión Internacional de Unidades y Medidas de Radiación (1970). Transferencia de energía lineal (PDF) . Washington DC ISBN 978-0913394090. Informe de la ICRU 16 . Consultado el 1 de diciembre de 2012 .{{cite book}}: CS1 maint: location missing publisher (link)
  6. ^ Yang, C.; Di Li, Di Li; Geng Wang, Geng Wang; Li Lin, Li Lin; Tasch, AF; Banerjee, S. (2002). "Modelo mecánico cuántico de poder de parada electrónico para iones en un gas de electrones libres". Tecnología de implantación de iones. 2002. Actas de la 14ª Conferencia Internacional sobre . págs. 556–559. doi :10.1109/IIT.2002.1258065. ISBN 0-7803-7155-0. S2CID  117936302.
  7. ^ P. Sigmund: Detención de iones pesados. Springer Tracts en física moderna vol. 204 (2004) ISBN 3-540-22273-1 
  8. ^ Zeb, M. Ahsan; Kohanoff, J.; Sánchez-Portal, D.; Arnau, A.; Juaristi, JI; Artacho, Emilio (31 de mayo de 2012). "Poder de frenado electrónico en oro: el papel de los electrones d y la anomalía H / He". Cartas de revisión física . 108 (22): 225504. arXiv : 1205.1728 . Código Bib : 2012PhRvL.108v5504Z. doi : 10.1103/PhysRevLett.108.225504. PMID  23003620. S2CID  2682651.
  9. ^ Ullah, Rafi; Corsetti, Fabiano; Sánchez Portal, Daniel; Artacho, Emilio (11 de marzo de 2015). "Poder de frenado electrónico en un semiconductor de banda prohibida estrecha desde los primeros principios". Revisión física B. 91 (12): 125203. arXiv : 1410.6642 . Código Bib : 2015PhRvB..91l5203U. doi : 10.1103/PhysRevB.91.125203. S2CID  31233417.
  10. ^ "Poder de frenado de iones ligeros". Archivado desde el original el 6 de febrero de 2012 . Consultado el 19 de febrero de 2014 .
  11. ^ Pablo, H (2006). "Una comparación de tablas de poder de parada recientes para iones ligeros y medio pesados ​​con datos experimentales y aplicaciones a la dosimetría de radioterapia". Instrumentos y métodos nucleares en la investigación en física B. 247 (2): 166-172. Código Bib : 2006NIMPB.247..166P. doi :10.1016/j.nimb.2006.01.059.
  12. ^ abc Comisión Internacional de Unidades y Medidas de Radiación (octubre de 2011). Seltzer, Stephen M. (ed.). «Cantidades y unidades fundamentales de las radiaciones ionizantes» (PDF) . Revista de la ICRU (edición revisada). 11 (1): NP.2–NP. doi :10.1093/jicru/ndr012. PMID  24174259. Informe ICRU 85a . Consultado el 14 de diciembre de 2012 .[ enlace muerto permanente ]
  13. ^ Huhtinen, Mika (2002). "Simulación de pérdida de energía no ionizante y formación de defectos en silicio". Instrumentos y métodos nucleares en la investigación en física B. 491 (1–2): 194–215. Código Bib : 2002NIMPA.491..194H. doi :10.1016/s0168-9002(02)01227-5.
  14. ^ Barry, Alabama; Houdayer, AJ; Hinrichsen, PF; Letourneau, WG; Vicente, J (1995). "La dependencia energética de las constantes de daño de vida útil en LED de GaAs para protones de 1 a 500 MeV". Transacciones IEEE sobre ciencia nuclear . 42 (6): 2104–2107. Código bibliográfico : 1995ITNS...42.2104B. doi : 10.1109/23.489259.
  15. ^ Lindström, G (2001). "Detectores de radiación de silicio duro: desarrollos de la colaboración RD48 (ROSE)". Instrumentos y métodos nucleares en la investigación en física A. 466 (2): 308–326. Código Bib : 2001NIMPA.466..308L. doi :10.1016/S0168-9002(01)00560-5. hdl : 11568/67464 .
  16. ^ abcd JF Ziegler, JP Biersack y U. Littmark. En The Stopping and Range of Ions in Matter, volumen 1, Nueva York, 1985. Pergamon. ISBN 0-08-022053-3 
  17. ^ JF Ziegler, JP Biersack y MD Ziegler: SRIM: la detención y el alcance de los iones en la materia, SRIM Co., 2008. ISBN 0-9654207-1-X 
  18. ^ ab sitio web de SRIM
  19. ^ Nordlund, K; Runeberg, N; Sundholm, D (1997). "Potenciales interatómicos repulsivos calculados utilizando métodos de Hartree-Fock y teoría de densidad funcional". Instrumentos y métodos nucleares en la investigación en física B. 132 (1): 45. Código bibliográfico : 1997NIMPB.132...45N. doi :10.1016/S0168-583X(97)00447-3.
  20. ^ Robinson, marca; Torrens, Ian (1974). "Simulación por computadora de cascadas de desplazamiento atómico en sólidos en la aproximación de colisión binaria". Revisión física B. 9 (12): 5008. Código bibliográfico : 1974PhRvB...9.5008R. doi : 10.1103/PhysRevB.9.5008.
  21. ^ Biersack, J; Haggmark, L (1980). "Un programa informático Monte Carlo para el transporte de iones energéticos en objetivos amorfos ☆". Instrumentos y métodos nucleares . 174 (1): 257. Código bibliográfico : 1980NucIM.174..257B. doi :10.1016/0029-554X(80)90440-1.
  22. ^ Robinson, M (1992). "Estudios de simulación por computadora de cascadas de colisiones de alta energía1". Instrumentos y métodos nucleares en la investigación en física B. 67 (1–4): 396–400. Código Bib : 1992NIMPB..67..396R. doi :10.1016/0168-583X(92)95839-J.
  23. ^ ab Nordlund, K (1995). "Simulación de dinámica molecular de rangos de iones en el rango de energía de 1 a 100 keV". Ciencia de Materiales Computacionales . 3 (4): 448–456. doi :10.1016/0927-0256(94)00085-Q.
  24. ^ ab Beardmore, Keith; Grønbech-Jensen, Niels (1998). "Esquema de dinámica molecular eficiente para el cálculo de perfiles dopantes por implantación de iones". Revisión física E. 57 (6): 7278–7287. arXiv : física/9901054 . Código bibliográfico : 1998PhRvE..57.7278B. CiteSeerX 10.1.1.285.6727 . doi : 10.1103/PhysRevE.57.7278. S2CID  13994369. 
  25. ^ Hobler, G. (2001). "Sobre el rango útil de aplicación de simulaciones de dinámica molecular en la aproximación de la interacción de retroceso". Instrumentos y métodos nucleares en la investigación en física B. 180 (1–4): 203. Código bibliográfico : 2001NIMPB.180..203H. doi :10.1016/s0168-583x(01)00418-9.
  26. ^ Caturla, M. (1996). "Procesamiento de silicio por haz de iones a energías keV: un estudio de dinámica molecular". Revisión física B. 54 (23): 16683–16695. Código bibliográfico : 1996PhRvB..5416683C. doi : 10.1103/PhysRevB.54.16683. PMID  9985796.
  27. ^ Smith, R. (1997). "Simulación de dinámica molecular de bombardeo de iones de Ni {100} de 0,1 a 2 keV". Efectos de la radiación y defectos en los sólidos . 141 : 425. doi : 10.1080/10420159708211586.
  28. ^ Duvenbeck, A. (2007). "Promoción de electrones y fricción electrónica en cascadas de colisiones atómicas". Nueva Revista de Física . 9 (2): 38. Código bibliográfico : 2007NJPh....9...38D. doi : 10.1088/1367-2630/9/2/038 .
  29. ^ Hou, M. (2000). "Deposición de grupos de AuN en superficies de Au (111). I. Modelado a escala atómica". Revisión física B. 62 (4): 2825. Código bibliográfico : 2000PhRvB..62.2825H. doi : 10.1103/PhysRevB.62.2825.
  30. ^ Björkas, C. (2009). "Evaluación de la relación entre la mezcla de haces de iones, el acoplamiento de electrones y fonones y la producción de daños en Fe". Instrumentos y métodos nucleares en la investigación en física B. 267 (10): 1830. Código bibliográfico : 2009NIMPB.267.1830B. doi :10.1016/j.nimb.2009.03.080.
  31. ^ Pronnecke, S. (1991). "El efecto de la pérdida de energía electrónica sobre la dinámica de los picos térmicos en Cu" (PDF) . Revista de investigación de materiales . 6 (3): 483. Código bibliográfico : 1991JMatR...6..483P. doi :10.1557/jmr.1991.0483. S2CID  96618269.
  32. ^ Duffy, DM (2007). "Incluidos los efectos de la detención electrónica y las interacciones electrón-ion en simulaciones de daños por radiación". Revista de Física: Materia Condensada . 17 (1): 016207. Código bibliográfico : 2007JPCM...19a6207D. doi :10.1088/0953-8984/19/1/016207. S2CID  122777435.
  33. ^ Tamm, A. (2016). "Interacción electrón-fonón dentro de la dinámica molecular clásica". Revisión física B. 94 (1): 024305. Código bibliográfico : 2016PhRvB..94a4305L. doi : 10.1103/PhysRevB.94.014305.
  34. ^ ab http://pdg.lbl.gov/2005/reviews/passagerpp.pdf [ URL básica PDF ]

Otras lecturas

enlaces externos