The structure of a melt: The case of liquid bismuth

Este estudio emplea simulaciones de Dinámica Molecular y modelado de Monte Carlo inverso para caracterizar la estructura atómica del bismuto líquido a 573 K, revelando una disposición local dominada por triángulos y cuadrados deformados que se manifiestan como picos específicos en las Distribuciones de Pares y en las Distribuciones de Ángulos de Plano.

Lo esencial

Imagina una olla de metal fundido, específicamente Bismuto, un elemento plateado que parece un cristal de colores arcoíris cuando se enfría. Cuando está caliente y líquido, los átomos en su interior bailan caóticamente, chocando entre sí como en una multitud desordenada. La gran pregunta que los científicos han estado planteando es: Aunque se muevan aleatoriamente, ¿estos átomos aún conservan algunas formas o patrones ocultos de cuando estaban sólidos?

Este artículo es como una historia de detectives de alta tecnología donde los autores utilizaron computadoras potentes para congelar en el tiempo esta danza atómica y buscar esos patrones ocultos.

La Configuración: Una Pista de Baile Virtual

Los investigadores construyeron una "supercelda" virtual (una caja diminuta) que contenía 216 átomos de bismuto. Imagina esta caja como una pista de baile.

• Comenzaron con los átomos a una temperatura fresca (300 K) y los calentaron hasta que estaban derritiéndose y fluyendo (573 K).
• Ejecutaron esta simulación 100 veces seguidas, luego mantuvieron la "danza" durante otros 500 pasos para asegurar que el líquido fuera estable.
• Para asegurarse de que no estaban viendo simplemente un accidente, iniciaron cuatro simulaciones diferentes con "impulsos iniciales" ligeramente distintos (velocidades aleatorias) para los átomos, solo para ver si el resultado era el mismo cada vez.

Las Herramientas: Capturando Instantáneas del Caos

Para entender la estructura, utilizaron dos herramientas principales:

1. PDF (Función de Distribución de Pares): Imagina tomar una foto de la multitud y medir la distancia entre cada par de personas. Si ves a mucha gente parada exactamente a 3 pies de distancia, obtienes un "pico" en tu gráfico. Esto te dice a qué distancia suelen sentarse los átomos.
2. PAD (Distribución de Ángulos de Plano): Esto mide los ángulos formados por tres átomos. Si el Átomo A, el Átomo B y el Átomo C forman un triángulo, ¿cuál es el ángulo en el Átomo B? Esto te dice la forma de los grupos.

El Gran Descubrimiento: El Misterio del "Hombro"

En el estado líquido, el gráfico de distancias (el PDF) suele tener un pico principal grande, seguido de un segundo pico. Pero en el Bismuto, hay una extraña "protuberancia" o hombro justo después del primer pico.

• La Controversia: Algunos científicos pensaron que este hombro era solo un error causado por cómo medían los datos en la vida real (como una foto borrosa). Otros pensaron que era una característica real del líquido.
• El Veredicto: Dado que los autores crearon estos datos completamente dentro de una computadora (sin fotos borrosas involucradas), y el hombro apareció cada vez, es real. Es una característica genuina del Bismuto líquido.

¿Qué Formas se Están Escondiendo en el Líquido?

Al analizar los ángulos y las distancias, los autores descubrieron que incluso en el líquido caótico, los átomos no son totalmente aleatorios. Están formando formas específicas, ligeramente aplastadas:

1. Triángulos Deformados: A los átomos les gusta agruparse de tres en tres, formando triángulos. Sin embargo, no son triángulos equiláteros perfectos; están aplastados o estirados. Esto corresponde a un ángulo específico de aproximadamente 53° a 58°.
2. Cuadrados Deformados: Los átomos también forman grupos de cuatro que parecen cuadrados o rombos, pero nuevamente, están distorsionados. Esto corresponde a ángulos alrededor de 85° a 90°.

El "Hombro" Explicado:
La misteriosa protuberancia de "hombro" en el gráfico de distancias es en realidad causada por las líneas diagonales de estos cuadrados aplastados. Cuando miras un cuadrado, la distancia de esquina a esquina (la diagonal) es más larga que la distancia de lado a lado. En el líquido, estas distancias diagonales crean ese extra de "protuberancia" en los datos.

La Conclusión

El artículo concluye que el Bismuto líquido no es simplemente una sopa aleatoria de átomos. Retiene una "memoria" de su estructura sólida. Incluso cuando se derrite, los átomos prefieren organizarse en triángulos aplastados y cuadrados aplastados.

Esto explica el "hombro" en los datos: es la huella dactilar de esas formas similares a cuadrados. Los autores también notaron que podría haber formas aún más complejas (como pentágonos o hexágonos) acechando en los datos, pero eso es un misterio para otro día.

En resumen: El líquido es caótico, pero es un caos estructurado, lleno de triángulos y cuadrados aplastados que dejan una marca distintiva en los datos, demostrando que la estructura líquida está más organizada de lo que pensábamos.

Resumen técnico

Resumen Técnico: La Estructura de un Fundido: El Caso del Bismuto Líquido

Planteamiento del Problema
A pesar de que el bismuto (Bi) es un material ampliamente estudiado con aplicaciones industriales significativas y transiciones de fase únicas (incluyendo superconductividad en estados amorfo y líquido), la estructura atómica precisa de su fase líquida sigue siendo objeto de debate. Específicamente, existen discrepancias en la literatura respecto a las Funciones de Distribución de Pares (PDF) y las Distribuciones de Ángulos de Plano (PAD) del bismuto líquido. Un punto central de controversia es la existencia y el origen de un "pseudo-pico conspicuo" (o hombro) observado en las PDF a aproximadamente 4.6 Å. Algunos investigadores atribuyen esta característica a arreglos atómicos específicos (como cadenas triangulares o estructuras cuadradas), mientras que otros argumentan que es un artefacto de los procedimientos de medición experimental o de la transformación de Fourier de los factores de estructura. Además, la relación entre estas características estructurales y las propiedades electrónicas mejoradas (como la superconductividad de alta-$T_c$) de las fases desordenadas del bismuto permanece poco clara. Este estudio tiene como objetivo resolver estas discrepancias investigando la topología atómica del bismuto líquido a 573 K para determinar si existen estructuras geométricas específicas que puedan explicar las características observadas en las PDF y potencialmente influir en las propiedades electrónicas.

Metodología
Los autores emplearon un enfoque computacional de múltiples pasos utilizando simulaciones de Dinámica Molecular (DM) y Monte Carlo Inverso (RMC):

1. Configuración del Sistema: Se construyeron cuatro superceldas, cada una conteniendo 216 átomos de bismuto, basadas en una multiplicación 3x3x3 de una celda unitaria de tipo diamante. Esta estructura inicialmente inestable fue elegida para facilitar la evolución hacia una topología desordenada. El sistema mantuvo la densidad líquida experimental de 10.029 g/cm³ a 573 K.
2. Dinámica Molecular (DM): Las simulaciones se realizaron utilizando el paquete Materials Studio (código DMol3) con Teoría del Funcional de la Densidad (DFT). El proceso involucró:
   • Calentamiento desde 300 K hasta 573 K en 100 pasos (de 19.4 fs cada uno).
   • Mantenimiento de la temperatura a 573 K durante 500 pasos utilizando un termostato Nosé-Hoover NVT.
   • Se utilizaron cuatro conjuntos iniciales distintos de velocidades aleatorias ($v_0$ a $v_3$) para probar la independencia de las estructuras finales respecto a las condiciones iniciales.
   • Los cálculos electrónicos utilizaron el funcional LDA-VWN, conjuntos de base numérica doble con polarización de función d (dnd) y Potenciales Pseudo de Núcleo Semi-Interior (dspp).
3. Análisis de Datos:
   • PDF: Calculadas para los últimos 1, 10, 25 y 100 pasos de la trayectoria de DM.
   • Generación RMC: Para obtener estructuras atómicas representativas, se aplicó Monte Carlo Inverso a las PDF promediadas de los últimos 100 pasos para cada conjunto de velocidades, generando cuatro configuraciones atómicas distintas.
   • PAD: Las Distribuciones de Ángulos de Plano se calcularon tanto para las estructuras promediadas por DM como para las generadas por RMC utilizando un código personalizado ("Correlation"). Se definió un corte de enlace de 4.185 Å (el mínimo entre el primer pico y el hombro) para identificar átomos enlazados.
4. Validación: La estabilidad de la simulación se verificó monitoreando la variación de energía, la cual resultó ser despreciable. Los resultados se compararon con datos experimentales existentes y otros estudios computacionales.

Resultados Clave

1. Convergencia y Consistencia: Las PDF obtenidas de los cuatro conjuntos iniciales diferentes de velocidades se volvieron indistinguibles al promediarlas sobre los últimos 100 pasos, confirmando que el sistema alcanzó un estado líquido representativo independiente de las condiciones iniciales. Las estructuras generadas por RMC produjeron PDF idénticas a los promedios de DM, validando el enfoque RMC para describir la fase líquida.
2. Funciones de Distribución de Pares (PDF): Las PDF calculadas para el bismuto líquido a 573 K exhiben dos picos prominentes a 3.25 Å y 6.55 Å, junto con un hombro distinto (pseudo-pico) a 4.6 Å.
   • El primer pico (3.25 Å) resulta de la coalescencia de los dos primeros picos cristalinos (3.12 Å y 3.47 Å).
   • El hombro a 4.6 Å corresponde a la coalescencia de los tercer y cuarto picos de vecinos cristalinos (4.52 Å y 4.72 Å).
   • Estos resultados se alinean con muchos estudios experimentales y computacionales que reportan el hombro, aunque existen algunas variaciones experimentales respecto a su prominencia.
3. Distribuciones de Ángulos de Plano (PAD): Las PAD revelan preferencias angulares específicas indicativas de orden geométrico local:
   • ~53° (DM) y ~58° (RMC): Estos picos sugieren la presencia de triángulos equiláteros deformados (tríadas).
   • ~85° (DM) y ~90° (RMC): Estos picos indican la existencia de cuadrados o rombos deformados.
   • 120°–140°: Una región amplia y menos pronunciada sugiere la presencia de estructuras geométricas más complejas y de orden superior.
4. Origen del Hombro: El estudio establece una correlación directa entre el hombro en la PDF y la distribución angular. La distancia interatómica correspondiente a la diagonal de las estructuras cuadradas deformadas (identificadas por los ángulos de 90°) coincide con la posición del hombro (~4.6 Å).

Significado y Afirmaciones
El artículo afirma resolver el debate sobre el hombro en la PDF del bismuto líquido. Los autores afirman que:

• Realidad Estructural: El hombro es una característica estructural genuina del bismuto líquido, no un artefacto de la medición experimental o errores en la transformación de Fourier. Surge del agrupamiento estadístico de átomos en geometrías específicas.
• Origen Geométrico: El hombro es causado por la coalescencia de los picos de segundos vecinos cristalinos y, más específicamente, por la existencia de triángulos equiláteros deformados y cuadrados deformados dentro de la estructura líquida. Las distancias diagonales en estos cuadrados deformados corresponden al hombro de 4.6 Å.
• Contribución Metodológica: Al combinar DM con RMC, el estudio proporciona un conjunto consistente de estructuras atómicas que reproducen tanto las características de PDF como de PAD, ofreciendo una imagen más clara de la topología local que los informes previos contradictorios.
• Implicaciones para las Propiedades: Aunque el artículo no demuestra definitivamente un vínculo causal con la superconductividad, postula que la identificación de estos arreglos atómicos específicos (desorden, triángulos y cuadrados) en la fase líquida proporciona el contexto estructural necesario para investigar por qué las fases desordenadas del bismuto exhiben propiedades electrónicas elevadas y superconductividad.

En conclusión, el trabajo identifica la topología atómica local del bismuto líquido como dominada por motivos triangulares y cuadrados deformados, los cuales explican directamente el controvertido hombro en la Función de Distribución de Pares.