stringtranslate.com

Factor Debye-Waller

El factor Debye-Waller (DWF), llamado así por Peter Debye e Ivar Waller , se utiliza en física de la materia condensada para describir la atenuación de la dispersión de rayos X o la dispersión coherente de neutrones causada por el movimiento térmico. [1] [2] También se le llama factor B , factor B atómico o factor de temperatura . A menudo, el "factor Debye-Waller" se utiliza como un término genérico que comprende el factor Lamb-Mössbauer de dispersión incoherente de neutrones y la espectroscopia Mössbauer .

La DWF depende del vector de dispersión q . Para un q dado , DWF( q ) da la fracción de dispersión elástica ; 1 – DWF( q ) da correspondientemente la fracción de dispersión inelástica (estrictamente hablando, esta interpretación de probabilidad no es verdadera en general [3] ). En estudios de difracción , solo la dispersión elástica es útil; en cristales, da lugar a picos de reflexión de Bragg distintos . Los eventos de dispersión inelástica son indeseables ya que causan un fondo difuso, a menos que se analicen las energías de las partículas dispersas, en cuyo caso llevan información valiosa (por ejemplo, en la dispersión inelástica de neutrones o en la espectroscopia de pérdida de energía de electrones ).

La expresión básica para el DWF está dada por

donde u es el desplazamiento de un centro de dispersión y denota promedio térmico o temporal.

Suponiendo la harmonicidad de los centros de dispersión en el material en estudio, la distribución de Boltzmann implica que se distribuye normalmente con media cero. Entonces, utilizando por ejemplo la expresión de la función característica correspondiente , la DWF toma la forma

Obsérvese que, si bien el razonamiento anterior es clásico, lo mismo se aplica a la mecánica cuántica.

Suponiendo también la isotropía del potencial armónico, se puede escribir

donde q , u son las magnitudes (o valores absolutos) de los vectores q , u respectivamente, y es el desplazamiento cuadrático medio . En las publicaciones cristalográficas, los valores de se dan a menudo donde . Nótese que si la onda incidente tiene una longitud de onda , y se dispersa elásticamente en un ángulo de , entonces

En el contexto de las estructuras de proteínas , se utiliza el término factor B. El factor B se define como

[4]

Se mide en unidades de Å 2 . Los factores B pueden tomarse como indicadores del movimiento vibratorio relativo de diferentes partes de la estructura. Los átomos con factores B bajos pertenecen a una parte de la estructura que está bien ordenada. Los átomos con factores B altos generalmente pertenecen a una parte de la estructura que es muy flexible. Cada registro ATOM ( formato de archivo PDB ) de una estructura cristalina depositada en el Protein Data Bank contiene un factor B para ese átomo.

Derivación

Introducción

Los experimentos de dispersión son un método común para aprender sobre los cristales . Estos experimentos suelen implicar una sonda (por ejemplo, rayos X o neutrones ) y un sólido cristalino. Una sonda bien caracterizada que se propaga hacia el cristal puede interactuar y dispersarse de una manera particular. Las expresiones matemáticas que relacionan el patrón de dispersión, las propiedades de la sonda, las propiedades del aparato experimental y las propiedades del cristal permiten entonces derivar las características deseadas de la muestra cristalina.

La siguiente derivación se basa en el capítulo 14 de The Oxford Solid State Basics [5] de Simon y en el informe Atomic Displacement Parameter Nomenclature de Trueblood et al. [6] (disponible en #Enlaces externos). Se recomienda consultar estas fuentes para una discusión más explícita. Se puede encontrar información de fondo sobre la mecánica cuántica involucrada en Modern Quantum Mechanics de Sakurai y Napolitano . [7]

Los experimentos de dispersión a menudo consisten en una partícula con un momento cristalino inicial que incide sobre un sólido. La partícula pasa a través de un potencial distribuido en el espacio, , y sale con un momento cristalino . Esta situación se describe mediante la regla de oro de Fermi , que da la probabilidad de transición por unidad de tiempo, , al estado propio de energía desde el estado propio de energía debido a la perturbación débil causada por nuestro potencial .

. (1)

Al insertar un conjunto completo de estados de posición y luego utilizar la expresión de onda plana que relaciona la posición y el momento, encontramos que el elemento de la matriz es simplemente una transformada de Fourier del potencial.

. (2)

Arriba, la longitud de la muestra se denota por . Ahora asumimos que nuestro sólido es un cristal periódico con cada celda unitaria etiquetada por un vector de posición reticular . La posición dentro de una celda unitaria está dada por un vector tal que la posición general en el cristal puede expresarse como . Debido a la invariancia traslacional de nuestras celdas unitarias, la distribución de potencial de cada celda es idéntica y .

. (3)

Ecuación de Laue

Según la fórmula de suma de Poisson :

. (4)

es un vector reticular recíproco del potencial periódico y es el volumen de su celda unitaria . Comparando (3) y (4), encontramos que la ecuación de Laue debe cumplirse para que se produzca dispersión:

. (5)

(5) es una afirmación de la conservación del momento cristalino. Las partículas dispersas en un cristal experimentan un cambio en el vector de onda igual a un vector reticular recíproco del cristal. Cuando lo hacen, la contribución al elemento de la matriz es simplemente una constante finita. Por lo tanto, encontramos un vínculo importante entre las partículas dispersas y el cristal dispersante. La condición de Laue, que establece que el momento cristalino debe conservarse, es equivalente a la condición de Bragg , que exige una interferencia constructiva para las partículas dispersas. Ahora que vemos cómo el primer factor de (3) determina si las partículas incidentes se dispersan o no, consideramos cómo el segundo factor influye en la dispersión.

Factor de estructura

El segundo término en el lado derecho de (3) es el factor de estructura .

. (6)

Para un vector reticular recíproco dado (correspondiente a una familia de planos reticulares etiquetados por índices de Miller ), la intensidad de las partículas dispersas es proporcional al cuadrado del factor de estructura.

. (7)

Enterrados en (6) hay aspectos detallados de la estructura cristalina que vale la pena distinguir y discutir.

Factor Debye-Waller

La consideración del factor de estructura (y nuestra suposición sobre la invariancia traslacional) se complica por el hecho de que los átomos en el cristal pueden desplazarse de sus respectivos sitios reticulares. Suponiendo que el potencial de dispersión es proporcional a la densidad de la materia dispersa, reescribimos el factor de estructura.

. (8)

A partir de aquí, se entiende que la integral se toma sobre la celda unitaria. es la densidad de la materia dispersa. Los corchetes angulares indican un promedio temporal de cada celda unitaria seguido de un promedio espacial sobre cada celda unitaria. Suponemos además que cada átomo se desplaza independientemente de los demás átomos.

. (9)

El número de átomos en la celda unitaria es y el factor de ocupación para el átomo es . representa el punto en la celda unitaria para el cual nos gustaría saber la densidad de materia dispersa. es la densidad de materia dispersa del átomo en una posición separada de la posición nuclear por un vector . es la función de densidad de probabilidad para el desplazamiento. es el sitio reticular de referencia desde el cual el átomo puede ser desplazado a una nueva posición . Si es suficientemente simétrico (por ejemplo, esféricamente simétrico), es simplemente la posición nuclear media. Al considerar la dispersión de rayos X, la densidad de materia dispersa consiste en la densidad de electrones alrededor del núcleo. Para la dispersión de neutrones, tenemos funciones ponderadas por una longitud de dispersión para el núcleo respectivo (ver pseudopotencial de Fermi ). Nótese que en la discusión anterior, asumimos que los átomos no eran deformables. Con esto en mente, (9) puede insertarse en la expresión (8) para el factor de estructura.

; . (10)

Ahora vemos que el factor de estructura general puede representarse como una suma ponderada de los factores de estructura correspondientes a cada átomo. Establezcamos el desplazamiento entre la ubicación en el espacio para la que nos gustaría conocer la densidad de dispersión y la posición de referencia para el núcleo igual a una nueva variable . Hagamos lo mismo para el desplazamiento entre las posiciones nucleares desplazadas y de referencia . Sustituyamos en (10).

. (11)

Dentro de los corchetes de (11), convolucionamos la densidad de materia dispersa del átomo con la función de densidad de probabilidad para algún desplazamiento nuclear. Luego, entre corchetes, realizamos la transformada de Fourier de la convolución resultante. El paso final es multiplicar por una fase que depende de la posición de referencia (por ejemplo, la media) del átomo . Pero, según el teorema de convolución , la transformada de Fourier de una convolución es lo mismo que multiplicar las dos funciones transformadas de Fourier. Establezcamos el desplazamiento entre la ubicación en el espacio para la que nos gustaría saber la densidad de dispersión y la posición del núcleo igual a una nueva variable .

. (12)

Sustituya (12) en (10).

. (13)

Eso es:

; , . (14)

es el factor de forma atómica del átomo ; determina cómo la distribución de materia dispersa alrededor de la posición nuclear influye en la dispersión. es el factor atómico de Debye-Waller; determina cómo la propensión al desplazamiento nuclear desde la posición de red de referencia influye en la dispersión. La expresión dada para en la introducción del artículo es diferente debido a 1) la decisión de tomar el promedio térmico o temporal, 2) la elección arbitraria del signo negativo en el exponencial y 3) la decisión de elevar al cuadrado el factor (que lo conecta más directamente con la intensidad observada).

Parámetro de desplazamiento anisotrópico, U

Una simplificación común de (14) es la aproximación armónica, en la que la función de densidad de probabilidad se modela como una gaussiana . Bajo esta aproximación, se ignora el desorden desplazativo estático y se supone que los desplazamientos atómicos están determinados completamente por el movimiento (se han considerado modelos alternativos en los que la aproximación gaussiana no es válida en otro lugar [8] ).

; ; . (15)

Hemos eliminado el índice atómico. pertenece a la red directa mientras que pertenecería a la red recíproca. Al elegir la base adimensional conveniente , garantizamos que tendrá unidades de longitud y describirá el desplazamiento. El tensor en (15) es el parámetro de desplazamiento anisotrópico. Con dimensión (longitud) , está asociado con los desplazamientos cuadráticos medios. Para el desplazamiento cuadrático medio a lo largo del vector unitario , simplemente tome . Los esquemas relacionados utilizan los parámetros o B en lugar de (consulte Trueblood et al. [6] para una discusión más completa). Finalmente, podemos encontrar la relación entre el factor de Debye-Waller y el parámetro de desplazamiento anisotrópico.

. (16)

A partir de las ecuaciones (7) y (14), el factor Debye-Waller contribuye a la intensidad observada de un experimento de difracción. Y con base en (16), vemos que nuestro factor de desplazamiento anisotrópico es responsable de determinar . Además, (15) muestra que puede estar directamente relacionado con la función de densidad de probabilidad para un desplazamiento nuclear desde la posición media. Como resultado, es posible realizar un experimento de dispersión en un cristal, ajustar el espectro resultante para los diversos valores atómicos y derivar la tendencia de cada átomo para el desplazamiento nuclear a partir de .

Aplicaciones

Un modelo de elipsoide térmico de probabilidad del 50% de H 8 Si 8 O 12 construido con ORTEP-3 [9] a partir de un archivo .cif en ICSD. [10] El análisis posterior a un experimento de difracción consiste en ajustar el espectro observado de partículas dispersas. U puede refinarse para cada átomo distinto durante el proceso. Para el modelo de probabilidad del 50% anterior, en la ecuación (15). Esto define una superficie de desplazamientos nucleares para cada U. Por lo tanto, esperamos que cada elipsoide varíe dependiendo del tipo y entorno de su átomo. Nótese que las superficies representan desplazamientos nucleares; los modelos de elipsoide térmico no deben confundirse con otros modelos (por ejemplo, densidad electrónica, radios de Van der Waals). Se muestran menos de 28 átomos debido a la redundancia de consideraciones de simetría.

Los parámetros de desplazamiento anisotrópico suelen ser útiles para visualizar la materia. A partir de (15), podemos definir elipsoides de probabilidad constante para los cuales , donde es alguna constante. Dichos " elipsoides de vibración " se han utilizado para ilustrar estructuras cristalinas. [9] Alternativamente, las superficies de desplazamiento cuadrático medio a lo largo de pueden definirse por . Consulte los enlaces externos "Galería de ORTEP trazados por rayos", "Artículo de 2005 de Rowsell et al ." y "Artículo de 2009 de Korostelev y Noller" para obtener más imágenes. Los parámetros de desplazamiento anisotrópico también se refinan en programas (por ejemplo, GSAS-II [11] ) para resolver espectros de dispersión durante el refinamiento de Rietveld .

Referencias

  1. ^ Adiós, Peter (1913). "Interferencia de Röntgenstrahlen und Wärmebewegung". Annalen der Physik (en alemán). 348 (1): 49–92. Código bibliográfico : 1913AnP...348...49D. doi : 10.1002/andp.19133480105.
  2. ^ Waller, Ivar (1923). "Zur Frage der Einwirkung der Wärmebewegung auf die Interferenz von Röntgenstrahlen". Zeitschrift für Physik A (en alemán). 17 (1): 398–408. Código Bib : 1923ZPhy...17..398W. doi :10.1007/BF01328696. S2CID  121962265.
  3. ^ Lipkin, Harry (2004). "Física de los factores de Debye-Waller". arXiv : cond-mat/0405023v1 .
  4. ^ Bahar, Ivet; Atilgan, Ali Rana; Erman, Burak (1997). "Evaluación directa de fluctuaciones térmicas en proteínas usando un potencial armónico de un solo parámetro". Folding and Design . 2 (3): 173–181. doi :10.1016/S1359-0278(97)00024-2. PMID  9218955.
  5. ^ Simon, Steven H. (20 de junio de 2013). Fundamentos del estado sólido de Oxford . OUP Oxford. ISBN 9780199680771.OCLC 1038069097  .
  6. ^ ab Trueblood, KN; Bürgi, HB; Burzlaff, H.; Dunitz, JD; Gramaccioli, CM; Schulz, HH; Shmueli, U.; Abrahams, SC (1996-09-01). "Nomenclatura de parámetros de desplazamiento atómico. Informe de un subcomité sobre nomenclatura de parámetros de desplazamiento atómico". Acta Crystallographica Sección A . 52 (5): 770–781. doi : 10.1107/s0108767396005697 . ISSN  0108-7673.
  7. ^ Sakurai, JJ; Napolitano, Jim (21 de septiembre de 2017). Mecánica cuántica moderna . Cambridge University Press. doi :10.1017/9781108499996. ISBN 9781108499996.
  8. ^ "3. Más allá de la aproximación gaussiana". ww1.iucr.org . Consultado el 15 de mayo de 2019 .
  9. ^ ab Burnett, MN; Johnson, CK (1 de julio de 1996). ORTEP-III: Programa de gráficos de elipsoides térmicos de Oak Ridge para ilustraciones de estructuras cristalinas (informe). doi :10.2172/369685. OSTI  369685.
  10. ^ Törnroos, KW (15 de noviembre de 1994). "Octahidridosilasesquioxano determinado por difracción de neutrones". Acta Crystallographica Sección C . 50 (11): 1646–1648. doi :10.1107/S0108270194005342.
  11. ^ "Ayuda para GSAS-II". subversion.xray.aps.anl.gov . Consultado el 30 de abril de 2019 .

Enlaces externos