Enlace y cohesión · Hoja de problemas

Hoja de problemas del módulo II.3

Seis problemas de dificultad creciente sobre enlace y cohesión, con pistas y sin soluciones. La brújula del módulo: aquí los errores ya no son de recuento ni de convenio, sino de EMPAREJAMIENTO — poner en una fórmula una energía y un volumen que no se refieren a la misma cosa. Hay dos comprobaciones de diez segundos que los cazan todos.

La hoja del II.1 avisaba de que casi todo error de cristalografía es un error de recuento; la del II.2, de que los de difracción son de convenio. Aquí el enemigo es otro y se llama emparejamiento: casi todos los fallos de este tema consisten en meter en la misma fórmula una energía y un volumen que no cuentan lo mismo. Energía de red frente a energía de cohesión. Por par de iones frente a por átomo. Volumen de la celda convencional frente a volumen por átomo. Ninguno produce un número absurdo — producen números creíbles, con las unidades correctas, y equivocados por factores 2, 4 u 8. Las dos comprobaciones que los cierran cuestan diez segundos y todos los problemas de abajo las piden. La primera: el número Kv0/E0K v_0/|E_0| tiene que caer en la banda de su familia —0,7 a 0,9 los iónicos, 2 a 3 los covalentes, 7 a 9 los moleculares y 1,2 a 4,3 los metálicos—, y si te sale ocho veces lo esperado has usado la celda convencional donde iba el átomo. Fíjate en la cuarta banda antes de usar la brújula: es la más ancha de las cuatro —un factor 3,5 de punta a punta— y se traga entera a la covalente, así que este número clasifica bien tres familias y no la cuarta. Eso es información y no un fallo, y los problemas 5 y 6 la explotan. La segunda: una energía de red y una energía de cohesión se diferencian en un ciclo de Born–Haber entero, que en un haluro alcalino son un factor 2,5, así que antes de comparar nada hay que decir cuál de las dos es. Y un tercer aviso que sólo aparece en el problema 1 y del que ningún artículo se ocupó del todo: hay sumas que dan lo que tú quieras según cómo las ordenes, y la comprobación no es sumar más términos sino mirar la carga del bloque.

Necesitas: el módulo II.3 completo. En concreto: el pozo E(r)E(r), K=Er02/(9v0)K = E''r_0^2/(9v_0) y el invariante Kv0/E0=mn/9K v_0/|E_0| = mn/9 del artículo 01 (problemas 1, 2, 4, 5 y 6); la constante de Madelung, el método de Evjen, el potencial de Born-Mayer y el ciclo de Born-Haber del artículo 02 (problemas 1, 2 y 6); el potencial de Lennard-Jones con sus sumas de red y la energía de punto cero del artículo 03 (problema 3); y la escala 2/med2\hbar^2/m_e d^2 del enlace covalente y el gas de electrones libres del artículo 04 (problemas 4 y 5). Del módulo II.1, la regla de los radios (problema 2) y las estructuras cúbicas (todos). Ningún problema te pide resolver un material del que los artículos hayan dado cifras: todos los datos que hacen falta vienen en los enunciados, y las constantes están en el resumen en frío al final del artículo 04.
Problema 1 ●○○ · La serie que se deja reordenar

Este problema baja a una y a dos dimensiones, donde todo se puede hacer con lápiz.

(a) El artículo 02 sumó la cadena alterna y obtuvo 2ln22\ln 2; eso se da por hecho. Lo que se pregunta aquí es cuánto cuesta obtenerlo por fuerza bruta: acota el error de truncar una serie alternada y di cuántos iones de la cadena hay que sumar para tener tres, cuatro y cinco cifras significativas. Compáralo con los 343 iones que le bastan al método de Evjen para dar cuatro cifras en tres dimensiones, y explica la paradoja aparente — que sumar un problema tridimensional salga más barato que el unidimensional.

(b) Ahora reordena esa misma serie tomando dos términos positivos por cada negativo, sin quitar ni añadir ninguno. Demuestra que converge a 32\tfrac32 del valor anterior. ¿Cuál de los dos números es el de un cristal, y con qué criterio se decide? Enuncia ese criterio de manera que se pueda aplicar en tres dimensiones.

(c) Ponle repulsión a la cadena: zλeR/ρz\lambda e^{-R/\rho} con ρ=0,32\rho = 0{,}32 Å y el mismo paso que el cloruro de sodio, R0=2,8201R_0 = 2{,}8201 Å. Cuidado con zz. Da la energía de red por par de iones y compárala con la del cristal tridimensional. ¿Qué fracción de la energía de Madelung del cloruro de sodio explica la sola fila de vecinos que pasa por el ion?

(d) Una red plana cuadrada de iones alternos, de paso RR. Aplica el método de Evjen con una sola capa —a mano, son ocho iones— y da el valor. Di si el bloque es neutro y cuánto te falta para el valor convergido, que es 1,615540. ¿Por qué en tres dimensiones una capa se queda un 17 % corta y aquí se queda un 20 %?

(e) Generaliza a dimensión dd cualquiera. El artículo 03 afirma que en tres dimensiones hace falta p>3p > 3; demuéstralo con el criterio de la integral, y para dd genérica y no sólo para 3. Con el resultado, rellena la tabla de d=1,2,3d = 1, 2, 3 —las tres dimensiones que este problema ha recorrido— por los cuatro exponentes que el interactivo del artículo 02 permite. ¿Cuál de las doce casillas es el peor caso de todos, y por cuántos órdenes se queda corto?

Pista

En (a) basta la cota clásica de una serie alternada —el error no supera el primer término omitido— y en (b) la reordenación se hace agrupando términos y comparando con la original. La trampa está en (c): en una cadena, cada ion tiene dos primeros vecinos y no seis, y ese número entra en la condición de equilibrio. Comprobación: el volumen por par de una cadena no es un volumen, así que la fórmula de KK del artículo 01 no se puede usar aquí — comprueba en su lugar que la profundidad sale con el mismo factor 1ρ/R01-\rho/R_0.

Problema 2 ●○○ · El yoduro de cesio y el volumen que cambia

El yoduro de cesio cristaliza en la estructura del cloruro de cesio con a=4,5667a = 4{,}5667 Å. Su módulo de compresibilidad medido es 12,5 GPa y su energía de red de Born–Haber, 585 kJ/mol. Los radios iónicos de Shannon son 1,67 Å para el Cs⁺ y 2,20 Å para el I⁻.

(a) Da la distancia al primer vecino y la energía de Madelung por par.

(b) Aquí está el problema. El artículo 02 dedujo K=(M/18R02)(1/ρ2/R0)K = (M/18R_0^2)(1/\rho - 2/R_0) para la estructura del cloruro de sodio, y ese 18 no es universal: es 9c9c con cc la constante que relaciona el volumen por par con R03R_0^3. Calcula cc para la estructura del cloruro de cesio y reescribe la fórmula. Comprueba tu cc de dos maneras que no se parezcan.

(c) Ajusta ρ\rho al módulo de compresibilidad medido y da λ\lambda; cuidado otra vez con el número de primeros vecinos. Comprueba que tu ρ\rho cae en la banda que el artículo 02 encontró para los seis haluros alcalinos.

(d) Predice la energía de red y compárala con el dato de Born–Haber. Di qué has ajustado y qué has predicho.

(e) La estructura del cloruro de cesio tiene la constante de Madelung más grande de las tres cúbicas. Explica entonces por qué el yoduro de sodio —mismo anión, catión más pequeño— no la adopta. Haz la comparación bien: los dos candidatos hay que evaluarlos a la distancia de contacto que cada estructura permite, no a la misma. Con r(Na+)=1,02r(\text{Na}^+) = 1{,}02 Å, ¿por qué margen gana la estructura del cloruro de sodio?

Pista

Todo sale de rehacer la deducción del artículo 02 sin dar por bueno ningún número que dependa de la estructura. La trampa son los dos números que cambian a la vez y no se ven: el volumen por par y el número de primeros vecinos. Comprobación: si tu ρ\rho se sale de la banda de los haluros o si Kv0/UK v_0/U no cae en la banda iónica del artículo 01, has arrastrado una geometría que no era.

Problema 3 ●●○ · El argón en pareja y en dos dimensiones

Con los parámetros del artículo 03, ε=0,0104\varepsilon = 0{,}0104 eV y σ=3,40\sigma = 3{,}40 Å.

(a) El dímero. El ejercicio 1 del artículo 03 ya dio la distancia de equilibrio del pozo de dos átomos, 21/6σ2^{1/6}\sigma, y su curvatura, 72ε/rmin272\varepsilon/r_{\min}^2; parte de ahí. Calcula la frecuencia de vibración del dímero Ar₂ —cuidado con la masa que entra— y su energía de punto cero. ¿Qué fracción de la profundidad se lleva? ¿Cuántos números de onda vale esa vibración, y en qué región del espectro cae? El dímero medido tiene De=12,34D_e = 12{,}34 meV y re=3,757r_e = 3{,}757 Å: compara y explica el signo de la discrepancia con lo que el callout del artículo 03 dice sobre el C6C_6.

(b) La monocapa. Una capa de argón adsorbida sobre grafito forma una red triangular plana. Calcula A6A_6 y A12A_{12} quedándote con tres capas de vecinos —6 a R0R_0, 6 a 3R0\sqrt3 R_0 y 6 a 2R02R_0—, y estima el error de cada una por separado; no es el mismo ni de lejos. Con ellas, da R0/σR_0/\sigma y Ecoh/εE_{\text{coh}}/\varepsilon de la monocapa, y compáralos con los del cristal. ¿Qué fracción de la cohesión tridimensional conserva una monocapa, y por qué no es la mitad si tiene la mitad de vecinos próximos?

(c) El invariante del artículo 01 se dedujo con v=cr3v = c\,r^3. Rehazlo en dos dimensiones, donde el «volumen» es un área y va como r2r^2, y di qué sustituye al 9.

(d) Estima la temperatura de fusión de la monocapa escalando desde los 83,8 K del argón másico. El valor experimental está entre 47 y 50 K: ¿te pasas o te quedas corto, y qué efecto que este cálculo ignora explica el signo?

Pista

El apartado (a) es el pozo de pares y no el del cristal: no lleva sumas de red por ninguna parte. La trampa de (b) es truncar las dos sumas y creer que el error es parecido — mira los exponentes antes de decidir cuántas capas te hacen falta de cada una. Comprobación: en (b) el cociente R0/σR_0/\sigma tiene que salir mayor que el del cristal, y si te sale menor has intercambiado las dos sumas.

Problema 4 ●●○ · El carburo de silicio, entre dos mundos

El carburo de silicio cúbico tiene la estructura de la blenda de cinc con a=4,3596a = 4{,}3596 Å y un módulo de compresibilidad medido de 220 GPa. Su entalpía de formación es −73 kJ/mol. Las entalpías de enlace sencillo C—C y Si—Si son 348 y 222 kJ/mol, y las electronegatividades de Pauling del silicio y el carbono, 1,90 y 2,55.

(a) Da la distancia Si—C y compárala con la media de las distancias Si—Si y C—C de los elementos. No coinciden: di en qué sentido falla y relaciónalo con el apartado (c).

(b) Calcula la energía de cohesión del carburo por fórmula, por átomo y por enlace, partiendo de las cohesiones del silicio y del carbono y de la entalpía de formación. Cuenta los enlaces con cuidado: no son dos por fórmula.

(c) Pauling propuso que el enlace A—B vale la media geométrica de A—A y B—B más un incremento iónico 96,5(Δχ)296{,}5\,(\Delta\chi)^2 kJ/mol. Aplícalo y compara con lo que has obtenido en (b). Después usa las dos ionicidades del carburo que el módulo ya da —la de Pauling en el resumen en frío y la de Phillips, 0,177— para decidir dos cosas: qué coordinación predice cada escala con la frontera 0,785, y si esa predicción habría sido igual de fiable en un compuesto que cayera cerca de la frontera. Contéstalo con números y no de palabra: calcula tú la ionicidad de Pauling del sulfuro y del seleniuro de cinc con la fórmula del artículo 01 —χ(Zn) = 1,65, χ(S) = 2,58, χ(Se) = 2,55— y ponlas al lado de las de Phillips, 0,623 y 0,676, que están en el resumen en frío. Uno de los dos está a una décima de la frontera: ¿qué escala te fiarías de usar ahí y por qué? Relaciona además el incremento iónico que has calculado con el acortamiento del enlace del apartado (a): ¿son el mismo efecto?

(d) El artículo 04 ajustó los módulos de compresibilidad del grupo 14 a K=2104d3,62K = 2104\,d^{-3{,}62} GPa con dd en ångströms. Aplícalo al carburo y compara con los 220 GPa medidos. Di por qué esa extrapolación es legítima y a qué compuesto del artículo 01 no se la aplicarías.

(e) Calcula Kv0/EcohK v_0/E_{\text{coh}} y el mnmn equivalente del carburo, y di en qué banda de la tabla del artículo 01 cae. Comprueba de paso que la escala 2/med2\hbar^2/m_e d^2 sigue dando el orden de magnitud del enlace.

Pista

La cohesión de un compuesto sale de las de sus elementos más la entalpía de formación, con el signo cuidado. La trampa es el recuento de enlaces: una fórmula SiC tiene dos átomos con cuatro enlaces cada uno, y cada enlace pertenece a dos átomos. Comprobación: la energía por enlace tiene que quedar entre la del Si—Si y la del C—C, y si se sale de ese intervalo has contado mal.

Problema 5 ●●● · El gas de electrones fuera de su terreno

El artículo 04 llevó el gas de electrones libres a los cinco alcalinos, todos monovalentes, y encontró cocientes entre 2,07 y 1,07. Aquí se le lleva a dos metales polivalentes: el aluminio, cúbico centrado en las caras con a=4,0495a = 4{,}0495 Å, tres electrones de valencia y K=76K = 76 GPa; y el plomo, también centrado en las caras con a=4,9508a = 4{,}9508 Å, cuatro electrones de valencia y K=46K = 46 GPa.

(a) Calcula para los dos la densidad de electrones, kFk_F, EFE_F y la temperatura de Fermi. Cuidado con los átomos por celda y con los electrones por átomo, que son dos recuentos distintos y ninguno es uno.

(b) Da el módulo de compresibilidad del gas de electrones y el cociente frente al medido. Hazlo por dos caminos que no se parezcan: desde la energía y desde la presión.

(c) Compara los dos cocientes con los de los cinco alcalinos. La tendencia que allí valía —el modelo mejora al bajar por la columna— no explica lo que sale aquí. Formula la regularidad que sí lo explica y dila en función de una sola variable.

(d) Dale la vuelta al modelo: como Kn5/3K \propto n^{5/3}, cada metal tiene una densidad de electrones «efectiva» que reproduciría su KK medido. Calcúlala para los dos y exprésala como fracción de la valencia nominal —ojo, en electrones por átomo los dos números no se parecen en nada—. Así sí sale un resultado sorprendentemente parecido en ambos: di qué crees que significa y qué comprobarías para saber si es casualidad.

(e) El artículo 04 dice que este modelo no da la cohesión porque su energía cinética es positiva. Ponle número al agujero: calcula la energía cinética total por átomo del gas de electrones del aluminio y di cuánta atracción por átomo haría falta añadir para cerrar los 3,39 eV de cohesión medida. Compara esa cifra con la propia cohesión y con la energía de Madelung de un cristal iónico del módulo, y di qué habría que añadirle al modelo — sin añadirlo.

Pista

Dos recuentos independientes —átomos por celda y electrones por átomo— que hay que multiplicar, no elegir. La trampa de (c) es buscar una tendencia con el tamaño del ion, que aquí no es la variable: mira qué tienen en común el aluminio y el plomo y que ninguno de los alcalinos tiene. Comprobación: EFE_F del aluminio es de las pocas que están tabuladas en cualquier manual, y si no te sale rodeando los doce electronvoltios, revisa el recuento antes de seguir.

Problema 6 ●●● · Tres sólidos sin etiqueta

Llegan tres muestras con estos datos medidos y sin decir de qué familia es cada una:

MuestraRedaa (Å)KK (GPa)Energía
Acúbica centrada en las caras4,08531002,95 eV/átomo (cohesión)
Bcúbica centrada en las caras3,52401804,44 eV/átomo (cohesión)
Cestructura del cloruro de sodio5,973219,9751 kJ/mol (energía de red)

(a) Calcula Kv0/EK v_0/E y el mnmn equivalente de las tres, y clasifícalas con la tabla del artículo 01.

(b) La muestra C tiene una trampa que se cobra a casi todo el mundo. Su energía viene dada como energía de red y no como cohesión. Calcula también su cohesión por átomo con I=5,139I = 5{,}139 eV y A=3,364A = 3{,}364 eV, repite el apartado (a) con esa energía y con el volumen por átomo, y di en qué familia la habrías metido. Explica por qué la versión equivocada también parece razonable.

(c) Confirma tu clasificación de C por un camino independiente: ajústale un ρ\rho de Born–Mayer y comprueba si su energía de red predicha cuadra con la medida y si ρ\rho cae en la banda de los haluros.

(d) Deduce que la derivada del módulo de compresibilidad respecto de la presión, evaluada a presión nula, vale

K0dKdPP=0=m+n+63K'_0 \equiv \left.\frac{dK}{dP}\right|_{P=0} = \frac{m+n+6}{3}

para un pozo A/rm+B/rn-A/r^m + B/r^n, y explica por qué es otra medida de la forma del pozo, independiente de Kv0/E0K v_0/|E_0|. Calcula qué predice para un cristal molecular (6, 12), para uno iónico (1, 9) y para tu muestra C.

(e) Los valores medidos de dK/dPdK/dP son 5,35 en el NaCl, 5,50 en el KCl, 7,40 en el argón sólido, 5,48 en el cobre y 4,15 en el MgO. Ahora aprovecha (d) del todo: con mnm\,n del apartado (a) y m+n=3K06m+n = 3K'_0 - 6 de éste, los dos exponentes quedan determinados como raíces de una cuadrática. Hazlo para los cinco —para el MgO, a=4,2117a = 4{,}2117 Å, K=160K = 160 GPa y U=3795U = 3795 kJ/mol; para el argón, R0=3,76R_0 = 3{,}76 Å, K=2,7K = 2{,}7 GPa y 0,080 eV/átomo; para el cobre, los datos del artículo 01—. Con las dos sales tienes una comprobación fuerte: mm tiene que salir 1, porque la atracción es Coulomb exacto. Uno de los cinco no da solución real: localízalo, di qué desigualdad se está violando y qué significa eso — que sus dos medidas no caben en ningún pozo de dos potencias. Y de los que sí la dan, di cuál es el que peor encaja con la pareja típica de su familia y por qué no es una sorpresa.

(f) Con todo lo anterior, escribe en cinco líneas el procedimiento que seguirías para clasificar el enlace de un sólido desconocido midiendo lo menos posible, y di qué medida no pedirías porque no discrimina.

Pista

El apartado (a) es aritmética con la tabla del artículo 01; todo el problema está en emparejar bien la energía con el volumen. La trampa es exactamente ésa y el apartado (b) la enseña en vivo: con la otra energía y el otro volumen, la muestra C cambia de familia sin que nada chirríe. Comprobación: el K0K'_0 de (d) y el Kv0/E0K v_0/|E_0| de (a) tienen que ser compatibles con la misma pareja (m, n) — si no lo son, uno de los dos está mal emparejado.

Fin del módulo II.3. Con esto ya sabes escribir la función que un cristal minimiza y leerle las tres cosas que se miden: dónde está su mínimo, cuánto hondo es y cuánto se curva. Sabes que el módulo de compresibilidad es la curvatura de E(V)E(V) y sabes calcularlo, que es la primera vez en el sitio que una propiedad mecánica sale de una cuenta y no de una tabla. Sabes que Kv0/E0=mn/9K v_0/|E_0| = mn/9 no depende de nada salvo de los dos exponentes, y que ese número separa las familias de enlace mejor que la energía. Sabes calcular una constante de Madelung —y, lo que es más difícil, sabes por qué no se puede sumar a lo bruto y qué hay que exigirle a un bloque para que su suma signifique algo—, ajustar un potencial de Born-Mayer y cerrar el ciclo de Born-Haber hasta el 3,3 eV por ion que el Nivel I dio como dato. Sabes de dónde sale el 1/r61/r^6 de van der Waals, y has visto el único caso del módulo en que el modelo predice tres propiedades de un sólido con dos parámetros medidos en un gas. Y sabes que el enlace covalente es 2t2t con t2/med2t \sim \hbar^2/m_ed^2, que sus 109,471° son arccos(1/3)\arccos(-1/3) y no un dato cristalográfico, y que la cohesión de un cristal covalente es la suma de sus enlaces contados uno a uno. En términos de Simon, el capítulo del enlace químico queda cubierto entero, con Kittel puesto para las tablas de haluros alcalinos.

Lo que falta es el movimiento. Todo lo de este módulo ha sido energía en el mínimo, con los átomos quietos en su sitio. Y ya se ha colado dos veces la señal de que eso no basta: la energía de punto cero, que en el neón se lleva el 27 % del pozo y sin la cual los cuatro gases nobles no cuadran, y la temperatura de Debye, que este módulo ha tenido que pedir prestada cinco veces —los cuatro gases nobles y el cobre— para calcular ese punto cero, sin poder deducirla de nada. Un pozo con una curvatura es un muelle, y un muelle no se está quieto — ni siquiera a cero kelvin. Nos falta además la mitad del pozo: aquí sólo se ha usado E(r0)E''(r_0), y el pozo es asimétrico, lo que el módulo I.3 llamó anarmonicidad para explicar la dilatación térmica sin poder deducirla.

El módulo II.4, «Vibraciones de la red: fonones», pone en movimiento exactamente esta curvatura. La constante de fuerza que sale de E(r0)E''(r_0) es la que sujeta cada átomo a su sitio; con ella se monta la cadena de osciladores acoplados, sale una relación de dispersión ω(k)\omega(\mathbf{k}), y esa k\mathbf{k} vive —por fin, y era el otro cabo suelto— dentro de la zona de Brillouin que construyó el II.2. Con una base de dos átomos aparecen dos ramas, y una de ellas explica por qué el cloruro de sodio absorbe en el infrarrojo. De ahí sale, ya en el II.5, el calor específico que el módulo I.3 midió sin poder deducir, con la ley T3T^3 y la temperatura de Debye que este módulo ha estado usando como dato prestado. El mapa lo dibujó el II.2, la fuerza la ha puesto el II.3, y el II.4 la pone a vibrar.