Red recíproca y difracción · Hoja de problemas

Hoja de problemas del módulo II.2

Seis problemas de dificultad creciente sobre la red recíproca y la difracción, con pistas y sin soluciones. La brújula del módulo: aquí los errores ya no son de recuento sino de convenio, y hay dos comprobaciones de diez segundos que cazan casi todos.

La hoja del módulo II.1 avisaba de que casi todos los errores de cristalografía son errores de recuento. Aquí cambia el enemigo: los errores de este módulo son de convenio. El 2π2\pi que sobra o que falta, el ángulo que se llama θ\theta cuando el aparato mide 2θ2\theta, el índice que se reduce a primos entre sí cuando en difracción no se reduce, la fórmula cúbica del espaciado aplicada a un cristal que no lo es. Ninguno de los cuatro produce un número absurdo: producen números creíbles, y por eso hay que cerrarles el paso con una comprobación mecánica en vez de con atención. Son dos, cuestan diez segundos y todos los problemas de abajo las piden. La primera: todo G|\mathbf{G}| que calcules, multiplicado por su dd, tiene que dar 2π2\pi — si da 1, estás en el otro convenio; si da 4π4\pi, te has comido un factor 2. La segunda: senθ=λ/2d\operatorname{sen}\theta = \lambda/2d nunca puede pasar de 1, y el número que sale del difractómetro es el doble del que entra en esa fórmula. Hay además un tercer aviso, y de él se encarga uno de los seis por las malas: las posiciones de los picos las decide la red y las intensidades las decide la base, así que diagnosticar una estructura mirando sólo dónde salen los picos no es difícil: con pocos picos es imposible, y el problema 5 está construido para que se vea con uno solo.

Necesitas: el módulo II.2 completo. En concreto: la construcción de los bi\mathbf{b}_i, el invariante VV=(2π)3V\cdot V^* = (2\pi)^3 y la fórmula general del espaciado del artículo 01 (problemas 1, 2 y 6); la primera zona de Brillouin y el recuento de vectores de onda del artículo 02 (problema 2); la ley de Bragg y la esfera límite de la construcción de Ewald del artículo 03 (problemas 1, 4 y 5); y el factor de estructura con el factor de Debye-Waller del artículo 04 (problemas 3 y 6). Del módulo II.1, los parámetros del magnesio y la distancia interplanar cúbica (problemas 2 y 3). Del módulo I.3, la temperatura de Debye y el criterio de Lindemann (problema 6). Ningún problema te pide resolver un material del que los artículos hayan dado cifras —el hierro aparece sólo como candidato a descartar en la tabla del problema 5—: todos los datos que hacen falta vienen en los enunciados, y las constantes y longitudes de onda están en el resumen en frío al final del artículo 04. Los problemas 2 y 3 son el mismo cristal en dos mitades: el 2 pone la geometría y el 3 decide qué se ve de ella.
Problema 1 ●○○ · Un difractograma sin etiqueta

Un polvo metálico sin identificar se mide con radiación de cobre, λ=1,5406\lambda = 1{,}5406 Å, y sus cinco primeras líneas salen a

2θ = 39,76° · 46,25° · 67,47° · 81,27° · 85,71°

Su densidad medida es 21,45 g/cm³. Es cúbico.

(a) Calcula dd para las cinco.

(b) Indexarlas es exactamente el procedimiento del ejercicio 1 del artículo 03, así que aquí se da hecho para no repetirlo: los cinco h2+k2+l2h^2+k^2+l^2 valen 3, 4, 8, 11 y 12. Di cuál de las tres redes cúbicas es y con qué regla lo decides. Y contesta después a lo que de verdad se pregunta aquí: esa lista se salta el 9 y el 10. Di exactamente qué ternas ocuparían esos dos sitios y por qué no están.

(c) Deduce aa de cada línea por separado y da la media. Las cinco no valen lo mismo: di cuál determina peor el parámetro de red y de dónde viene esa jerarquía. Cuidado con mezclar dos cosas que aquí no ordenan igual: lo que cada línea se desvía de la media y la incertidumbre que arrastra por el redondeo de su ángulo.

(d) Con aa y la densidad, calcula la masa molar y di qué metal es. Cuidado: la masa que sale deja tres elementos dentro del 1,5 %, y hay que elegir con un argumento, no con el más próximo. Las masas molares y las densidades de los candidatos son datos de tabla periódica y no vienen en el resumen en frío: darlos aquí sería dar la respuesta.

(e) Da G|\mathbf{G}| de la primera línea en el convenio de este módulo y en el de la cristalografía, y comprueba Gd=2π|\mathbf{G}|\,d = 2\pi. Un catálogo publica ese punto recíproco con el segundo número: calcula el dd y el aa que obtendría quien lo tomara por el primero, y di en qué factor exacto se equivocaría.

Pista

Bragg en los dos sentidos: de 2θ2\theta a dd, y de dd a aa con la lista de enteros de por medio. La trampa está en el paso final: la masa molar que sale deja tres vecinos de la tabla periódica dentro del 1,5 %, y sólo uno de los tres tiene la densidad del enunciado. Comprobación: los cinco valores de aa tienen que coincidir en la cuarta cifra, y el que no lo haga delata un θ\theta donde iba 2θ2\theta.

Problema 2 ●○○ · El magnesio no es cúbico, y eso cambia la receta

El magnesio es hexagonal con a=3,2094a = 3{,}2094 Å y c=5,2107c = 5{,}2107 Å —los datos del resumen en frío del II.1—, y sus vectores primitivos son

a1=a(1,0,0),a2=a(12,32,0),a3=c(0,0,1).\mathbf{a}_1 = a(1,0,0), \qquad \mathbf{a}_2 = a\bigl(-\tfrac12,\tfrac{\sqrt3}{2},0\bigr), \qquad \mathbf{a}_3 = c(0,0,1).

(a) Calcula VV y los tres bi\mathbf{b}_i con sus módulos. Da el ángulo que forman b1\mathbf{b}_1 y b2\mathbf{b}_2 y compáralo con los 120° de los directos. Comprueba las nueve igualdades biaj=2πδij\mathbf{b}_i\cdot\mathbf{a}_j = 2\pi\delta_{ij} y el invariante de los volúmenes.

(b) Calcula dd y G|\mathbf{G}| de las familias (100), (002), (101) y (110) por los dos caminos —con los bi\mathbf{b}_i y con la fórmula hexagonal del artículo 01— y comprueba G=2π/d|\mathbf{G}| = 2\pi/d en las cuatro. Ordénalas por G|\mathbf{G}| creciente.

(c) Construye la primera zona de Brillouin: di qué poliedro es y de qué G\mathbf{G} sale cada cara. Calcula las distancias de Γ al centro de una cara lateral (punto M), a una arista vertical (punto K) y a una tapa (punto A), y el volumen de la zona por dos caminos que no se parezcan.

(d) La zona del magnesio es mucho más ancha que alta. Da el cociente ΓM/ΓA y escribe la desigualdad ΓM > ΓA como una condición sobre c/ac/a —quedará un número limpio—. Comprueba que la cumplen los cinco metales hexagonales cuyo c/ac/a da el ejemplo resuelto 1 del artículo 04 del II.1, y di qué forma tendría la zona de un cristal que estuviera justo en la frontera.

(e) ¿Cuántos vectores de onda distintos admite un milímetro cúbico de magnesio? El paso de la malla lo da el artículo 02 y no depende del material, así que la pregunta es otra: ¿de qué depende el número, y de qué no? Cambia la arista del cubo a 0,1 mm y di en qué factor cambia cada cosa.

Pista

Tres productos vectoriales, y el resto sale de biaj=2πδij\mathbf{b}_i\cdot\mathbf{a}_j = 2\pi\delta_{ij}. La trampa es el ángulo: los b\mathbf{b} no forman entre sí los 120° de los a\mathbf{a}, y arrastrarlos al producto escalar estropea todos los espaciados. Comprobación: cada G|\mathbf{G}| tiene que coincidir con 2π/d2\pi/d de la fórmula hexagonal, y el producto de los volúmenes con (2π)3(2\pi)^3.

Problema 3 ●●○ · La extinción que este módulo no ha deducido

Continúa el problema 2 y usa sus espaciados.

El magnesio tiene dos átomos por celda, en (0,0,0) y en (13,23,12)(\tfrac13,\tfrac23,\tfrac12). El artículo 04 dedujo las extinciones de la centrada en el cuerpo, la centrada en las caras y el diamante; ésta no, y es la de la mitad de los metales de la tabla periódica.

(a) Escribe ShklS_{hkl} para esa base y demuestra que S2/f2|S|^2/f^2 sólo puede tomar cuatro valores. Da los cuatro y la condición sobre (h,k,l)(h,k,l) de cada uno. Sugerencia: agrupa el exponente en (h+2k)/3+l/2(h+2k)/3 + l/2.

(b) Enuncia la regla de extinción en una línea y da las cinco primeras (hkl)(hkl) extinguidas, ordenadas por G|\mathbf{G}| creciente.

(c) Calcula 2θ2\theta y las intensidades relativas de las tres primeras líneas del magnesio con radiación de cobre, con IS2×m×LPI \propto |S|^2 \times m \times \text{LP}mm es la multiplicidad y LP el factor de Lorentz-polarización (1+cos22θ)/(sen2θcosθ)(1+\cos^2 2\theta)/(\operatorname{sen}^2\theta\cos\theta)—. En el patrón publicado del magnesio la línea más intensa es con diferencia la (101) y las dos anteriores se quedan alrededor de una cuarta parte de ella: ¿lo reproduces? Fíjate en que la que gana no es la de mayor S2|S|^2 y di qué la hace ganar.

(d) Contesta en el espacio real, sin factores de estructura: la (001) falta y la (002) no; y la (100), con l=0l = 0, tampoco falta. Di dónde están exactamente los planos de átomos de la primera familia y por qué con la tercera no pasa lo mismo. El mecanismo lo tienes en el artículo 01 para el grafito; lo que se pide aquí es la cuenta: qué fracción de la periodicidad separa las dos capas del compacto hexagonal, y qué desfase introduce eso en la (001) y en la (100).

(e) El carburo de wolframio tiene esa misma geometría con dos átomos distintos: W en (0,0,0) y C en (13,23,12)(\tfrac13,\tfrac23,\tfrac12). Con fZf \approx Z, ¿qué les pasa a las líneas extinguidas? Da el cociente de intensidades entre una de ellas y una de la clase más fuerte, y di qué habría que poner en lugar del carbono para que la extinción volviera a ser exacta.

Pista

Dos átomos, dos exponenciales: agrupa el exponente y mira cuándo el corchete se anula, que es cuando vale un semientero. La trampa es dar por buena aquí la regla de la misma paridad, que es de otra red y en un hexagonal no significa nada. Comprobación: el promedio de S2|S|^2 sobre muchas ternas, pesando cada clase por su fracción de (hkl)(hkl), tiene que dar 2f22f^2.

Problema 4 ●●○ · Tres límites que se confunden todo el rato

El germanio tiene la estructura del diamante con a=5,6575a = 5{,}6575 Å. Se va a medir en un difractómetro de sobremesa que sólo barre 2θ2\theta entre 10° y 120°, y que admite tubo de cobre (λ=1,5406\lambda = 1{,}5406 Å) o de molibdeno (λ=0,7107\lambda = 0{,}7107 Å).

(a) Escribe la lista de las (hkl)(hkl) permitidas de la estructura del diamante hasta h2+k2+l2=51h^2+k^2+l^2 = 51 con sus dd, y comprueba que ese 51 no es un corte arbitrario: es exactamente donde se acaba lo accesible con cobre por el límite físico dλ/2d \ge \lambda/2. ¿Cuántas líneas son? ¿Y cuántas permitiría el molibdeno —basta el número—?

(b) La ventana del aparato es otro límite y no es el mismo. ¿A qué intervalo de dd corresponde con cada radiación? Da los dos extremos.

(c) ¿Cuántas líneas caen dentro de la ventana con cada radiación? Ojo con el conjunto sobre el que se cuenta, que no es el mismo: con cobre las candidatas son las 13 de (a) y no hay más; con molibdeno hay que contar sobre todas las que él permite. De las que se pierden, di por qué extremo se pierde cada una — las dos radiaciones las pierden por el mismo lado, y eso no es lo interesante. Lo interesante es el riesgo del otro extremo, que sólo corre el molibdeno: di cuál es y por cuánto margen se salva aquí.

(d) Calcula la separación media entre líneas consecutivas dentro de la ventana con cada radiación. El ejemplo resuelto 1 del artículo 03 ya dijo que elegir λ\lambda es elegir entre número de datos y separación angular: ponle cifras a esa frase para este material y di cuántas veces peor separadas quedan las líneas con molibdeno. Después decide algo que el artículo no decide: con una anchura instrumental de 0,10°, ¿es el difractograma de molibdeno inutilizable para identificar la fase, o el problema es otro?

(e) Cuantifica la segunda mitad de (d). Para cada línea de la ventana con cobre, ¿qué error en 2θ2\theta admite si aa ha de salir con una parte en 10⁴? Parte de Δa/a=cotgθ  Δθ\Delta a/a = \operatorname{cotg}\theta\;\Delta\theta y dedúcela antes de usarla.

(f) Exactamente una de las líneas de la ventana con cobre no es una familia sino dos, que caen en el mismo sitio. Di cuál es, cuáles son las dos ternas, cuáles sus multiplicidades y qué le pasa a la altura de esa línea. ¿Podría separarlas girando la muestra, o hace falta otra cosa?

Pista

Todo sale de senθ=λ/2d\operatorname{sen}\theta = \lambda/2d en las dos direcciones: una vez para el límite físico y otra para los extremos de la ventana. La trampa es contar líneas creyendo que cuentas familias: hay una posición a la que llegan dos ternas y sus intensidades se suman. Comprobación: el intervalo de dd de la ventana tiene que quedar dentro del que permite λ/2\lambda/2, nunca al revés.

Problema 5 ●●● · Un pico, un parámetro de red y una pregunta imposible

Llega este correo:

«He medido un metal cúbico desconocido con Cu-Kα. El primer pico sale a 2θ=44,45±0,102\theta = 44{,}45 \pm 0{,}10^\circ, así que d=2,0365d = 2{,}0365 Å. Como la primera línea de un cúbico es la (111), el parámetro de red es a=d3a = d\sqrt3. ¿Me confirmas el número?»

(a) Recalcula dd y propaga la incertidumbre del ángulo. ¿Está bien ese primer paso?

(b) Audita la frase «la primera línea de un cúbico es la (111)». Enumera todas las hipótesis compatibles con lo medido y da el aa —con su incertidumbre— que resulta de cada una.

(c) Con esta lista de candidatos, di cuáles son compatibles con la medida y cuáles quedan fuera, midiendo la discrepancia en unidades de la propia incertidumbre.

MetalRedaa (Å)MM (g/mol)ρ\rho (g/cm³)
Vbcc3,024050,9426,11
Crbcc2,883951,9967,19
Fe αbcc2,866555,8457,874
Nifcc3,524058,6938,908
Gediamante5,657572,6305,323

(d) El compañero contesta que también ha pesado la muestra y le sale ρ=8,90\rho = 8{,}90 g/cm³, así que «ya está resuelto». Escribe ρ=nM/NAa3\rho = nM/N_A a^3, cuenta incógnitas y ecuaciones, y da la masa molar que saldría con esa densidad para cada una de las hipótesis de (b). ¿Cierra la balanza la pregunta o la traslada?

(e) Di cuál es la medida adicional más barata que sí la cierra, qué valor tomaría para cada tipo de red, y por qué esa medida no necesita conocer ni aa ni λ\lambda.

(f) Escribe la respuesta al correo. Incluye lo único que está determinado por el dato que te han dado, y una hipótesis que el compañero está usando sin haberla escrito y que ni la segunda línea ni la balanza detectarían: piensa en qué pasaría si la muestra fuera una disolución sólida de dos de los metales de la tabla, cuyo parámetro de red interpola entre los de sus componentes.

Pista

Antes de calcular nada, cuenta las incógnitas del problema y las ecuaciones que te han dado. La trampa es la frase «como es cúbico, la primera línea es la (111)»: eso no es un dato sino una hipótesis sobre la red, y hay otras tres igual de legítimas. Comprobación: la medida que cierra la pregunta tiene que ser un cociente, porque un cociente no lleva dentro ni el parámetro ni la longitud de onda.

Problema 6 ●●● · El plomo vibra demasiado

El plomo es cúbico centrado en las caras con a=4,9508a = 4{,}9508 Å y masa molar 207,2 g/mol. Su temperatura de Debye es 105 K —la más baja de los metales corrientes— y funde a 600,6 K. Se mide con radiación de cobre.

(a) Calcula u2\langle u^2\rangle y urmsu_{\text{rms}} a 293 K con el modelo de Debye del artículo 04 y compáralos con los del cobre. Antes de hacer la integral, di cuánto vale T/θDT/\theta_D y qué esperas de ese cociente.

(b) Deduce BB y calcula la pérdida de intensidad de todas las líneas del plomo accesibles con Cu-Kα. ¿Cuál es la primera que pierde más de la mitad, y a qué 2θ2\theta está? Y una advertencia que el ejemplo resuelto 2 del artículo 04 ya hizo con el cobre: la θD\theta_D que ajusta el calor específico no es la que ajusta BB, así que el BB que salga de los 105 K del enunciado es una cota inferior. El valor medido del plomo a 293 K está entre 2,2 y 2,6 Ų: repite la pregunta con 2,4 y di cuánto se adelanta la respuesta.

(c) Repite el cálculo en el punto de fusión y da el cociente entre urmsu_{\text{rms}} y la distancia al primer vecino. Compáralo con el criterio de Lindemann. ¿Sale por encima o por debajo del 10 %, y qué dirías de la regla a la vista de ese número?

(d) El callout del artículo 01 dejó dicho que un cristal finito de tamaño LL no tiene puntos recíprocos sino manchas de anchura ΔG2π/L\Delta G \approx 2\pi/L, y ahí lo dejó. Conviértelo en algo medible: deduce la anchura angular β\beta que eso produce en un difractograma —diferencia la ley de Bragg— y comprueba que sale β=λ/(Lcosθ)\beta = \lambda/(L\cos\theta). Aplícalo a un polvo de plomo de 12 nm y da β\beta de la primera línea y de la última.

(e) Con esa anchura, ¿se solapan dos líneas vecinas del plomo a 12 nm? Busca el par más apretado de líneas a ángulos distintos y di por debajo de qué tamaño de grano se funden en una. Lo de «ángulos distintos» no es una precaución retórica: en el difractograma del plomo hay dos posiciones a las que llegan dos familias cada una, separadas exactamente 0°, igual que la del apartado (f) del problema 4. Localízalas y di por qué no cuentan aquí. Y por el otro lado: ¿a partir de qué tamaño β\beta queda por debajo de la anchura instrumental de un aparato corriente, 0,10°?

(f) Los dos efectos borran la parte derecha del difractograma y no son lo mismo. El artículo 04 ya dijo que el Debye–Waller rebaja los picos sin ensancharlos y que la intensidad que se va reaparece como fondo difuso: completa el cuadro para el otro efecto —qué le pasa a su anchura, qué le pasa a su área y adónde va lo que falta— y usa esa diferencia para diseñar la medida más simple que los separe sin tocar la muestra.

Pista

Dos efectos independientes sobre el mismo patrón: uno rebaja las alturas y el otro reparte cada pico a lo ancho. La trampa es el factor 3 que separa u2\langle u^2\rangle del ux2\langle u_x^2\rangle con el que se define BB, y el criterio de Lindemann usa el primero. Comprobación: repite el cálculo de (a) con el cobre y tiene que devolverte el u2\langle u^2\rangle del artículo 04.

Fin del módulo II.2. Con esto ya sabes construir la red donde vive todo lo que un cristal le hace a una onda: los tres bi\mathbf{b}_i a partir de tres productos vectoriales, con Ghkl\mathbf{G}_{hkl} perpendicular a los planos que nombra y de longitud 2π/d2\pi/d en los siete sistemas cristalinos, no sólo en el cúbico. Sabes que la celda de Wigner-Seitz de esa red es la primera zona de Brillouin, que dentro caben tantos vectores de onda como celdas primitivas tiene el cristal, y que su borde es exactamente el sitio donde el cristal difracta. Sabes que la condición de Laue y la ley de Bragg son la misma condición en dos idiomas, leerlas las dos sobre la esfera de Ewald, y por qué con una longitud de onda λ\lambda no hay manera de ver nada menor que λ/2\lambda/2. Y sabes separar las dos preguntas que la gente mezcla: dónde salen los picos lo decide la red y cuánto valen lo decide la base, con el factor de estructura dando las extinciones, el factor de forma atómico explicando que el hidrógeno no se vea, los neutrones arreglándolo con una longitud de dispersión que ni cae con el ángulo ni tiene por qué ser positiva, y Debye–Waller diciendo cuánta señal se lleva la temperatura. Sobre todo, está pagada la deuda con la que se cerró el II.1: aquel módulo citó una docena de veces «el cobre mide 3,6149 Å» y ahora sabes de dónde sale ese número, con qué aparato, con qué error y por qué hay que medirlo en las líneas de ángulo alto. En términos de Simon, el capítulo 13 queda cubierto entero y el 14 también.

Lo que falta es la energía. Mira hacia atrás y compruébalo: aquí han salido kilovoltios de un tubo, temperaturas de Debye y desplazamientos cuadráticos medios, pero ni una sola vez la energía que mantiene unidos los átomos. Todo ha sido geometría y fases —dónde están los átomos y cómo se suman sus ondas—, y en ningún momento se ha preguntado por qué están ahí. El parámetro de red ha entrado siempre como dato medido, incluso en la hoja que acabas de terminar: al plomo del problema 6 le hemos dado 4,9508 Å y al germanio del 4, 5,6575 Å, sin ninguna idea de por qué no son otros. Y la temperatura de Debye del plomo, que en el problema 6 decide cuántas líneas sobreviven, ha entrado igual que entró en el módulo I.3: como un número de tabla. Un cristal es una configuración de equilibrio de un montón de átomos, y de la función que se minimiza para llegar a ella no hemos dicho nada.

El módulo II.3, «Enlace y cohesión», escribe esa función. Un término por cada manera que tiene la materia de sujetarse —iónica, covalente, metálica y de van der Waals— con sus energías en electronvoltios: la constante de Madelung, que convierte una suma condicionalmente convergente sobre toda la red en un número; el potencial de Lennard-Jones; y el módulo de compresibilidad deducido de la curvatura de E(V)E(V) en el mínimo, que es la primera vez en el sitio en que una propiedad mecánica medible sale de un cálculo y no de una tabla. Con eso se contesta por fin por qué el cloruro de sodio tiene el parámetro que tiene y por qué el plomo es blando y el diamante duro. Y abre lo siguiente solo: si hay una energía con un mínimo, hay una curvatura alrededor de él, y apartar un átomo de su sitio cuesta un trabajo que se devuelve oscilando. Ésos son los fonones del módulo II.4 — que se moverán, uno por uno, dentro de la zona de Brillouin construida aquí. El mapa ya está dibujado; el II.3 pone la fuerza y el II.4, el movimiento.