Name: analyzing melts and fluids from ab initio molecular dynamics simulations with the umd package
Uploaded: 2021-09-17T13:00:00
Duration: 6 min 37 s
Description: Watch this Scientific Journal Video about Analyzing Melts and Fluids from Ab Initio Molecular Dynamics Simulations with the UMD Package at JoVE.com

Este trabajo fue apoyado por el Consejo Europeo de Investigación (ERC) en el marco del programa de investigación e innovación Horizonte 2020 de la Unión Europea (acuerdo de subvención número 681818 IMPACT a RC), por la Dirección de Física y Química Extrema del Observatorio de Carbono Profundo, y por el Consejo de Investigación de Noruega a través de su esquema de financiación de Centros de Excelencia, proyecto número 223272. Reconocemos el acceso a los superordenadores GENCI a través de la serie stl2816 de subvenciones informáticas eDARI, al superordenador Irene AMD a través del proyecto PRACE RA4947 y al superordenador Fram a través del UNINETT Sigma2 NN9697K. FS fue apoyado por un proyecto Marie Skłodowska-Curie (acuerdo de subvención ABISSE No.750901).

<ol>
	<li>Frenkel, D., Smit, B. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=.">.</a> Understanding Molecular Simulation. From Algorithms to Applications. , (2001).</li><li>Allen, M. P., Tildesley, D. J., Allen, T. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=.">.</a> Computer Simulation of Liquids. , (1989).</li><li>Zepeda-Ruiz, L. A., Stukowski, A., Oppelstrup, T., Bulatov, V. V. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Probing+the+limits+of+metal+plasticity+with+molecular-dynamics+simulations.">Probing the limits of metal plasticity with molecular-dynamics simulations.</a> Nature Publishing Group. 550 (7677), 492-495 (2017).</li><li>Humphrey, W., Dalke, A., Schulten, K. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=VMD:+Visual+molecular+dynamics.">VMD: Visual molecular dynamics.</a> Journal of Molecular Graphics &amp; Modeling. 14 (1), 33-38 (1996).</li><li>Momma, K., Izumi, F. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=VESTA3+for+three-dimensional+visualization+of+crystal,+volumetric+and+morphology+data.">VESTA3 for three-dimensional visualization of crystal, volumetric and morphology data.</a> Journal of Applied Crystallography. 44 (6), 1272-1276 (2011).</li><li>Brehm, M., Kirchner, B. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=TRAVIS+-+A+free+Analyzer+and+Visualizer+for+Monte+Carlo+and+Molecular+Dynamics+Trajectories.">TRAVIS - A free Analyzer and Visualizer for Monte Carlo and Molecular Dynamics Trajectories.</a> Journal of Chemical Information and Modeling. 51 (8), 2007-2023 (2011).</li><li>Stixrude, L. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Visualization-based+analysis+of+structural+and+dynamical+properties+of+simulated+hydrous+silicate+melt.">Visualization-based analysis of structural and dynamical properties of simulated hydrous silicate melt.</a> Physics and Chemistry of Minerals. 37 (2), 103-117 (2009).</li><li>Kresse, G., Hafner, J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Ab+initio+Molecular-Dynamics+for+Liquid-Metals.">Ab initio Molecular-Dynamics for Liquid-Metals.</a> Physical Review B. 47 (1), 558-561 (1993).</li><li>Gygi, F. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Architecture+of+Qbox:+A+scalable+first-principles+molecular+dynamics+code.">Architecture of Qbox: A scalable first-principles molecular dynamics code.</a> IBM Journal of Research and Development. 52 (1-2), 137-144 (2008).</li><li>Harvey, J. P., Asimow, P. D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Current+limitations+of+molecular+dynamic+simulations+as+probes+of+thermo-physical+behavior+of+silicate+melts.">Current limitations of molecular dynamic simulations as probes of thermo-physical behavior of silicate melts.</a> American Mineralogist. 100 (8-9), 1866-1882 (2015).</li><li>Caracas, R., Hirose, K., Nomura, R., Ballmer, M. D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Melt-crystal+density+crossover+in+a+deep+magma+ocean.">Melt-crystal density crossover in a deep magma ocean.</a> Earth and Planetary Science Letters. 516, 202-211 (2019).</li><li>Green, M. S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Markoff+Random+Processes+and+the+Statistical+Mechanics+of+Time-Dependent+Phenomena.+II.+Irreversible+Processes+in+Fluids.">Markoff Random Processes and the Statistical Mechanics of Time-Dependent Phenomena. II. Irreversible Processes in Fluids.</a> The Journal of Chemical Physics. 22 (3), 398-413 (1954).</li><li>Kubo, R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Statistical-Mechanical+Theory+of+Irreversible+Processes.+I.+General+Theory+and+Simple+Applications+to+Magnetic+and+Conduction+Problems.">Statistical-Mechanical Theory of Irreversible Processes. I. General Theory and Simple Applications to Magnetic and Conduction Problems.</a> Journal of the Physical Society of Japan. 12 (6), 570-586 (1957).</li><li>Lin, S. T., Blanco, M., Goddard, W. A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=The+two-phase+model+for+calculating+thermodynamic+properties+of+liquids+from+molecular+dynamics:+Validation+for+the+phase+diagram+of+Lennard-Jones+fluids.">The two-phase model for calculating thermodynamic properties of liquids from molecular dynamics: Validation for the phase diagram of Lennard-Jones fluids.</a> The Journal of Chemical Physics. 119 (22), 11792-11805 (2003).</li><li>Meyer, E. R., Kress, J. D., Collins, L. A., Ticknor, C. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Effect+of+correlation+on+viscosity+and+diffusion+in+molecular-dynamics+simulations.">Effect of correlation on viscosity and diffusion in molecular-dynamics simulations.</a> Physical Review E. 90 (4), 1198-1212 (2014).</li><li>Soubiran, F., Militzer, B., Driver, K. P., Zhang, S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Properties+of+hydrogen,+helium,+and+silicon+dioxide+mixtures+in+giant+planet+interiors.">Properties of hydrogen, helium, and silicon dioxide mixtures in giant planet interiors.</a> Physics of Plasmas. 24 (4), 041401-041407 (2017).</li><li>Flyvbjerg, H., Petersen, H. G. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Error+estimates+on+averages+of+correlated+data.">Error estimates on averages of correlated data.</a> The Journal of Chemical Physics. 91 (1), 461-466 (1989).</li><li>Tuckerman, M. E. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=.">.</a> Statistical mechanics: theory and molecular simulation. , (2010).</li><li>McDonough, W. F., Sun, S. S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=The+composition+of+the+Earth.">The composition of the Earth.</a> Chemical Geology. 120, 223-253 (1995).</li><li>Elkins-Taton, L. T. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Magma+oceans+in+the+inner+solar+system.">Magma oceans in the inner solar system.</a> Annual Review of Earth and Planetary Sciences. 40, 113-139 (2012).</li><li>Lock, S. J., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=The+origin+of+the+Moon+within+a+terrestrial+synestia.">The origin of the Moon within a terrestrial synestia.</a> J. Geophysical Research: Planets. 123, 910-951 (2018).</li><li>Solomatova, N. V., Caracas, R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Pressure-induced+coordination+changes+in+a+pyrolitic+silicate+melt+from+ab+initio+molecular+dynamics+simulations.">Pressure-induced coordination changes in a pyrolitic silicate melt from ab initio molecular dynamics simulations.</a> Journal of Geophysical Research: Solid Earth. 124, 11232-11250 (2019).</li></ol>

Los autores no tienen nada que revelar.

El paquete UMD ha sido diseñado para funcionar mejor con simulaciones ab initio, donde el número de instantáneas generalmente se limita a decenas a cientos de miles de instantáneas, con unos pocos cientos de átomos por unidad de celda. Las simulaciones más grandes también son manejables siempre que la máquina en la que se ejecuta el posprocesamiento tenga suficientes recursos de memoria activos. El código se distingue por la variedad de propiedades que puede calcular y por su licencia de código abierto.
Los archivos umd.dat son apropiados para los conjuntos que conservan el número de partículas sin cambios a lo largo de la simulación. El paquete UMD puede leer archivos derivados de cálculos donde varía la forma y el volumen del cuadro de simulación. Estos cubren los cálculos más comunes, como NVT y NPT, donde el número de partículas, N, temperatura T, volumen, V y / o presión, P, se mantienen constantes.
Para el momento comienzan la función de distribución de pares, así como todas las escrituras que necesitan estimar las distancias interatómicas, como las escrituras de especiación, funcionan solo para células unitarias ortogonales, es decir, para células cúbicas, tetragonales y ortorrómbicas, donde los ángulos entre los ejes son de 90 °.
Las principales líneas de desarrollo para la versión 2.0 son la eliminación de la restricción de ortogonalidad para distancias y la adición de más características para los scripts de especiación: analizar enlaces químicos individuales, analizar los ángulos interatómicos e implementar la segunda esfera de coordinación. Con la ayuda de la colaboración externa, estamos trabajando en la migración del código a una GPU para un análisis más rápido en sistemas más grandes.

Los fluidos y los fundidos son vectores activos de transporte químico y físico en entornos naturales. Las elevadas tasas de difusión atómica favorecen los intercambios y reacciones químicas, la baja viscosidad junto con la flotabilidad variable favorecen la gran transferencia de masa y las relaciones de densidad de fusión cristalina favorecen la estratificación dentro de los cuerpos planetarios. La ausencia de una red periódica, las altas temperaturas típicas requeridas para alcanzar el estado fundido y la dificultad para el enfriamiento hacen que la determinación experimental de una serie de propiedades obvias, como la densidad, la difusión y la viscosidad, sea extremadamente desafiante. Estas dificultades hacen que los métodos computacionales alternativos sean herramientas fuertes y útiles para investigar esta clase de materiales.
Con el advenimiento de la potencia de cálculo y la disponibilidad de superordenadores, actualmente se emplean dos grandes técnicas de simulación atomística numérica para estudiar el estado dinámico de un sistema atomístico no cristalino, Monte Carlo1 y la dinámica molecular (MD)1,2. En las simulaciones de Monte Carlo, el espacio de configuración se muestrea aleatoriamente; Los métodos de Monte Carlo muestran una escala lineal en paralelización si todas las observaciones de muestreo son independientes entre sí. La calidad de los resultados depende de la calidad del generador de números aleatorios y de la representatividad del muestreo. Los métodos de Monte Carlo muestran escala lineal en paralelización si el muestreo es independiente entre sí. En dinámica molecular (DM) el espacio configuracional es muestreado por trayectorias atómicas dependientes del tiempo. A partir de una configuración dada, las trayectorias atómicas se calculan integrando las ecuaciones newtonianas de movimiento. Las fuerzas interatómicas se pueden calcular utilizando potenciales interatómicos modelo (en la DM clásica) o utilizando métodos de primeros principios (in ab initio, o primeros principios, MD). La calidad de los resultados depende de la longitud de la trayectoria y su capacidad para no sentirse atraído por los mínimos locales.
Las simulaciones de dinámica molecular contienen una gran cantidad de información, toda relacionada con el comportamiento dinámico del sistema. Las propiedades medias termodinámicas, como la energía interna, la temperatura y la presión, son bastante estándar para calcular. Pueden extraerse de los archivos de salida de las simulaciones y promediarse, mientras que las cantidades relacionadas directamente con el movimiento de los átomos, así como con su relación mutua, deben calcularse después de la extracción de las posiciones y velocidades atómicas.
En consecuencia, se ha dedicado mucho esfuerzo a visualizar los resultados, y varios paquetes están disponibles hoy en día, en diferentes plataformas, de código abierto o no [Ovito3, VMD4, Vesta5, Travis6, etc.]. Todas estas herramientas de visualización tratan eficientemente las distancias interatómicas, y como tales, permiten el cálculo eficiente de las funciones de distribución de pares y los coeficientes de difusión. Varios grupos que realizan simulaciones de dinámica molecular a gran escala tienen software patentado para analizar varias otras propiedades resultantes de las simulaciones, a veces en shareware u otras formas de acceso limitado a la comunidad, y a veces limitado en alcance y uso a algunos paquetes específicos. Algoritmos sofisticados para extraer información sobre enlaces interatómicos, patrones geométricos y termodinámica se desarrollan e implementan en algunos de estos paquetes3,4,5,6,7, etc.
Aquí proponemos el paquete UMD, un paquete de código abierto escrito en Python para analizar la salida de simulaciones de dinámica molecular. El paquete UMD permite el cálculo de una amplia gama de propiedades estructurales, dinámicas y termodinámicas (Figura 1). El paquete está disponible a través del sitio web de GitHub (https://github.com/rcaracas/UMD_package) y a través de una página dedicada (http://moonimpact.eu/umd-package/) del proyecto ERC IMPACT como un paquete de acceso abierto.
Para que sea universal y más fácil de manejar, nuestro enfoque es extraer primero toda la información relacionada con el estado termodinámico y las trayectorias atómicas del archivo de salida de la ejecución de dinámica molecular real. Esta información se almacena en un archivo dedicado, cuyo formato es independiente del paquete MD original donde se ejecutó la simulación. Llamamos a estos archivos archivos "umd", que significa Dinámica Molecular Universal. De esta manera, nuestro paquete UMD puede ser utilizado fácilmente por cualquier grupo ab initio con cualquier software, todo con un mínimo esfuerzo de adaptación. El único requisito para utilizar el paquete actual es escribir el analizador apropiado desde la salida del software MD en particular en el formato de archivo umd, si aún no existe. Por el momento, proporcionamos dichos analizadores para los paquetes VASP8 y QBox9 .
<img alt="Figure 1" class="xfigimg" src="/files/ftp_upload/61534/61534fig01.jpg"> Figura 1: Diagrama de flujo de la biblioteca UMD. Las propiedades físicas están en azul, y los principales scripts de Python y sus opciones están en rojo. <a href="https://www.jove.com/files/ftp_upload/61534/61534fig01large.jpg" target="_blank">Haga clic aquí para ver una versión más grande de esta figura.</a>
Los archivos umd son archivos ASCII; la extensión típica es "umd.dat" pero no obligatoria. Todos los componentes de análisis pueden leer archivos ASCII del formato umd, independientemente de la extensión del nombre real. Sin embargo, algunos de los scripts automáticos diseñados para realizar estadísticas rápidas a gran escala en varias simulaciones buscan específicamente archivos con la extensión umd.dat. Cada propiedad física se expresa en una línea. Cada línea comienza con una palabra clave. De esta manera, el formato es altamente adaptable y permite agregar nuevas propiedades al archivo umd, al tiempo que conserva su legibilidad en todas las versiones. Las primeras 30 líneas del archivo umd de la simulación de pirolita a 4.6 GPa y 3000 K, utilizadas a continuación en la discusión, se muestran en la Figura 2.
<img alt="Figure 2" class="xfigimg" src="/files/ftp_upload/61534/61534fig02.jpg"> Figura 2: El comienzo del archivo umd que describe la simulación de pirolita líquida a 4.6 GPa y 3000 K. El encabezado va seguido de la descripción de cada instantánea. Cada propiedad se escribe en una línea, que contiene el nombre de la propiedad física, los valores y las unidades, todas separadas por espacios. <a href="https://www.jove.com/files/ftp_upload/61534/61534fig02large.jpg" target="_blank">Haga clic aquí para ver una versión más grande de esta figura.</a>
Todos los archivos umd contienen un encabezado que describe el contenido de la celda de simulación: el número de átomos, electrones y tipos atómicos, así como detalles para cada átomo, como su tipo, símbolo químico, número de electrones de valencia y su masa. Una línea vacía marca el final del encabezado y lo separa de la parte principal del archivo umd.
Luego se detalla cada paso de la simulación. Primero, se dan los parámetros termodinámicos instantáneos, cada uno en una línea diferente, especificando (i) el nombre del parámetro, como energía, tensiones, presión hidrostática equivalente, densidad, volumen, parámetros de red, etc., (ii) su(s) valor(es), y (iii) sus unidades. A continuación se presenta una tabla que describe los átomos. Una línea de cabecera da las diferentes medidas, como posiciones cartesianas, velocidades, cargas, etc., y sus unidades. Luego, cada átomo se detalla en una línea. Por grupos de tres, correspondientes a los tres ejes x, y, z , las entradas son: las posiciones reducidas, las posiciones cartesianas plegadas en la celda de simulación, las posiciones cartesianas (que tienen debidamente en cuenta el hecho de que los átomos pueden atravesar varias celdas unitarias durante una simulación), las velocidades atómicas y las fuerzas atómicas. Las dos últimas entradas son escalares: carga y momento magnético.
Dos bibliotecas principales garantizan el correcto funcionamiento de todo el paquete. La biblioteca umd_process.py se ocupa de los archivos umd, como la lectura y la impresión. La biblioteca crystallography.py se ocupa de toda la información relacionada con la estructura atómica real. La filosofía subyacente de la biblioteca crystallography.py es tratar la red como un espacio vectorial. Los parámetros de la celda unitaria junto con su orientación representan los vectores base. El "Espacio" tiene una serie de atributos escalares (volumen específico, densidad, temperatura y número específico de átomos), propiedades termodinámicas (energía interna, presión, capacidad calorífica, etc.) y una serie de propiedades tensoriales (tensión y elasticidad). Los átomos pueblan este espacio. La clase "Lattice" define este conjunto, junto con varios cálculos cortos, como volumen específico, densidad, obtención de la red recíproca del directo, etc. La clase "Átomos" define los átomos. Se caracterizan por una serie de propiedades escalares (nombre, símbolo, masa, número de electrones, etc.) y una serie de propiedades vectoriales (la posición en el espacio, ya sea relativa a la base vectorial descrita en la clase Lattice, o relativa a coordenadas cartesianas universales, velocidades, fuerzas, etc.). Aparte de estas dos clases, la biblioteca de crystallography.py contiene una serie de funciones para realizar una variedad de pruebas y cálculos, como distancias atómicas o multiplicación de celdas. La tabla periódica de los elementos también se incluye como diccionario.
Los diversos componentes del paquete umd escriben varios archivos de salida. Como regla general, todos son archivos ASCII, todas sus entradas están separadas por pestañas y se hacen lo más autoexplicativas posible. Por ejemplo, siempre indican claramente la propiedad física y sus unidades. Los archivos umd.dat cumplen plenamente con esta regla.

<table>
	<tbody>
		<tr>
			<td>getopt library</td>
			<td></td>
			<td></td>
			<td>open-source</td>
		</tr>
		<tr>
			<td>glob library</td>
			<td></td>
			<td></td>
			<td>open-source</td>
		</tr>
		<tr>
			<td>matplotlib library</td>
			<td></td>
			<td></td>
			<td>open-source</td>
		</tr>
		<tr>
			<td>numpy library</td>
			<td></td>
			<td></td>
			<td>open-source</td>
		</tr>
		<tr>
			<td>os library</td>
			<td></td>
			<td></td>
			<td>open-source</td>
		</tr>
		<tr>
			<td>Python software</td>
			<td>The Python Software Foundation</td>
			<td>Version 2 and 3</td>
			<td>open-source</td>
		</tr>
		<tr>
			<td>random library</td>
			<td></td>
			<td></td>
			<td>open-source</td>
		</tr>
		<tr>
			<td>re library</td>
			<td></td>
			<td></td>
			<td>open-source</td>
		</tr>
		<tr>
			<td>scipy library</td>
			<td></td>
			<td></td>
			<td>open-source</td>
		</tr>
		<tr>
			<td>subprocess library</td>
			<td></td>
			<td></td>
			<td>open-source</td>
		</tr>
		<tr>
			<td>sys library</td>
			<td></td>
			<td></td>
			<td>open-source</td>
		</tr>
	</tbody>
</table>

analyzing melts and fluids from ab initio molecular dynamics simulations with the umd package

Hemos desarrollado un paquete de código abierto basado en Python para analizar los resultados derivados de simulaciones de dinámica molecular ab initio de fluidos. El paquete es el más adecuado para aplicaciones en sistemas naturales, como fundidos de silicato y óxido, fluidos a base de agua y varios fluidos supercríticos. El paquete es una colección de scripts de Python que incluyen dos bibliotecas principales que se ocupan de los formatos de archivo y la cristalografía. Todos los scripts se ejecutan en la línea de comandos. Proponemos un formato simplificado para almacenar las trayectorias atómicas y la información termodinámica relevante de las simulaciones, que se guarda en archivos UMD, que significa Dinámica Molecular Universal. El paquete UMD permite el cálculo de una serie de propiedades estructurales, de transporte y termodinámicas. Comenzando con la función de distribución de pares, define las longitudes de enlace, construye una matriz de conectividad interatómica y, finalmente, determina la especiación química. Determinar la vida útil de las especies químicas permite realizar un análisis estadístico completo. Luego, los scripts dedicados calculan los desplazamientos cuadráticos medios para los átomos, así como para las especies químicas. El análisis de autocorrelación implementado de las velocidades atómicas arroja los coeficientes de difusión y el espectro vibratorio. El mismo análisis aplicado sobre las tensiones produce la viscosidad. El paquete está disponible a través del sitio web de GitHub y a través de su propia página dedicada del proyecto ERC IMPACT como paquete de acceso abierto.

La pirolita es un modelo de fusión de silicato multicomponente (0.5Na2O 2CaO 1.5Al2O3 4FeO 30MgO 24SiO2) que mejor se aproxima a la composición del silicato a granel de la Tierra, el promedio geoquímico o de todo nuestro planeta, excepto su núcleo a base de hierro19. La Tierra primitiva estuvo dominada por una serie de eventos de fusión a gran escala20, el último podría haber envuelto a todo el planeta, después de su condensación para el disco protolunar21. La pirolita representa la mejor aproximación a la composición de tales océanos de magma a escala planetaria. En consecuencia, estudiamos ampliamente las propiedades físicas de la fusión de pirolita en el rango de temperatura de 3.000\u20125.000 K y el rango de presión de GPa de 0\u2012150 a partir de simulaciones de dinámica molecular ab initio en la implementación de VASP. Estas condiciones termodinámicas caracterizan completamente las condiciones oceánicas de magma más extremas de la Tierra. Nuestro estudio es un excelente ejemplo de un uso exitoso del paquete UMD para todo el análisis en profundidad de los fundidos22. Calculamos la distribución y las longitudes medias de los enlaces, trazamos los cambios en la coordinación catiónico-oxígeno y comparamos nuestros resultados con estudios experimentales y computacionales previos sobre silicatos amorfos de diversas composiciones. Nuestro análisis en profundidad ayudó a descomponer los números de coordinación estándar en sus constituyentes básicos, describir la presencia de poliedros de coordinación exóticos en la masa fundida y extraer vidas útiles para todos los poliedros de coordinación. También destacó la importancia del muestreo en simulaciones en términos de longitud de la trayectoria y también del número de átomos presentes en el sistema que se modela. En cuanto al post-procesamiento, el análisis UMD es independiente de estos factores, sin embargo, deben tenerse en cuenta al interpretar los resultados proporcionados por el paquete UMD. Aquí, mostramos algunos ejemplos de cómo el paquete UMD se puede utilizar para extraer varios rasgos característicos de las fundiciones, con una aplicación a la pirolita fundida.
La función de distribución de pares Si-O obtenida del script gofrs_umd.py muestra que el radio de la primera esfera de coordinación, que es el primer mínimo de la función g(r), se encuentra alrededor de 2,5 angstroms en T = 3000 K y P = 4,6 GPa. El máximo de g(r) se encuentra en 1.635 Å, esta es la mejor aproximación a la longitud de la curva. La cola larga se debe a la temperatura. Usando este límite como la distancia de enlace si-O, el análisis de especiación muestra que las unidades de SiO4, que pueden durar hasta unos pocos picosegundos, dominan la fusión (Figura 5). Hay una parte importante de la masa fundida que muestra polimerización parcial, como lo refleja la presencia de dímeros como Si2O7 y trímeros como las unidades Si3Ox. Su vida útil correspondiente es del orden del picosegundo. Todos los polímeros de orden superior tienen una vida útil considerablemente más corta.
<img alt="Figure 5" class="xfigimg" src="/files/ftp_upload/61534/61534fig05.jpg"> Figura 5: Vida útil de la especie química Si-O. La especiación se identifica en una fusión multicomponente a 4,6 GPa y 3000 K. Las etiquetas marcan los monómeros siO3, SiO4 y SiO5 y los diversos polímeros SixOy. <a href="https://www.jove.com/files/ftp_upload/61534/61534fig05large.jpg" target="_blank">Haga clic aquí para ver una versión más grande de esta figura.</a>
Los diferentes valores de los pasos verticales y horizontales, definidos por las banderas –z y –v anteriores, producen varios muestreos del MSD (Figura 6). Incluso grandes valores de z y v son suficientes para definir las pendientes y, por lo tanto, los coeficientes de difusión de los diferentes átomos. La ganancia en el tiempo para el post-procesamiento es notable cuando se va a valores grandes de z y v. El MSD ofrece un criterio de validación muy fuerte para la calidad de las simulaciones. Si la parte de difusión del MSD no es lo suficientemente larga, eso es una señal de que la simulación es demasiado corta y no alcanza el estado fluido en un sentido estadístico. El requisito mínimo para la parte difusiva del MSD depende en gran medida del sistema. Se puede requerir que todos los átomos cambien su sitio al menos una vez en la estructura de la masa fundida para que sea considerada como un fluido10. Un excelente ejemplo con aplicaciones en ciencias planetarias son las complejas fundiciones de silicato a altas presiones cercanas o incluso por debajo de su línea liquidus11. Los átomos de Si, los principales cationes formadores de redes, cambian de sitio después de más de dos docenas de picosegundos. Las simulaciones más cortas que este umbral estarían considerablemente inframuestreando el posible espacio de configuración. Sin embargo, como los aniones coordinadores, es decir, los átomos de O, se mueven más rápido que los átomos centrales de Si, pueden compensar parte de la lenta movilidad del Si. Como tal, todo el sistema podría cubrir un mejor muestreo del espacio de configuración que el asumido solo a partir de los desplazamientos de Si.
<img alt="Figure 6" class="xfigimg" src="/files/ftp_upload/61534/61534fig06.jpg"> Figura 6: Desplazamientos cuadráticos medios (TME). Los MSD se ilustran para algunos tipos atómicos de una fusión de silicato de múltiples componentes. El muestreo con varios pasos horizontales y verticales, z y v, produce resultados consistentes. Círculos sólidos: -z 50 –v 50. Círculos abiertos: -z 250 –v 500. <a href="https://www.jove.com/files/ftp_upload/61534/61534fig06large.jpg" target="_blank">Haga clic aquí para ver una versión más grande de esta figura.</a>
Finalmente, las funciones de VCA atómica producen el espectro vibratorio de la masa fundida. La Figura 7 muestra el espectro en las mismas condiciones de presión y temperatura que las anteriores. Representamos las contribuciones de los átomos de Mg, Si y O, así como el valor total. A la frecuencia cero hay un valor finito del espectro, que corresponde al carácter de difusión de la masa fundida. La extracción de las propiedades termodinámicas del espectro vibratorio debe eliminar este carácter difusivo similar al gas de cero, pero también tener en cuenta adecuadamente su desintegración a frecuencias más altas.
<img alt="Figure 7" class="xfigimg" src="/files/ftp_upload/61534/61534fig07.jpg"> Figura 7: El espectro vibratorio de la fundición de pirolita. La parte real de la transformada de Fourier de la función de autocorrelación velocidad-velocidad atómica produce el espectro vibratorio. Aquí el espectro se calcula para una fusión de silicato de múltiples componentes. Los fluidos tienen un carácter difusivo similar al gas distinto de cero a una frecuencia cero. <a href="https://www.jove.com/files/ftp_upload/61534/61534fig07large.jpg" target="_blank">Haga clic aquí para ver una versión más grande de esta figura.</a>

Watch this Scientific Journal Video about Analyzing Melts and Fluids from Ab Initio Molecular Dynamics Simulations with the UMD Package at JoVE.com

Analyzing Melts and Fluids from Ab Initio Molecular Dynamics Simulations with the UMD Package

Los fundidos y fluidos son vectores ubicuos de transporte de masa en sistemas naturales. Hemos desarrollado un paquete de código abierto para analizar simulaciones de dinámica molecular ab initio de dichos sistemas. Calculamos las propiedades estructurales (unión, clusterización, especiación química), transporte (difusión, viscosidad) y termodinámicas (espectro vibratorio).

Análisis de fundiciones y fluidos de ab initio Molecular Dynamics Simulations con el paquete UMD

We have developed a Python-based open-source package to analyze the results stemming from ab initio molecular-dynamics simulations of fluids. The package ...

El análisis de simulaciones de dinámica molecular de primer principio nos permite predecir con precisión las propiedades físicas y químicas de los fluidos. Este método se puede aplicar a cualquier análisis atómico en física y química. Para comenzar, extraiga cada conjunto específico de propiedades físicas utilizando uno o más scripts de Python dedicados del paquete.Ejecute todos los scripts en la línea de comandos. Todos ellos emplean una serie de banderas, que son lo más consistentes posible de un guión a otro. Transforme la salida de la simulación MD realizada en un código de primeros principios en un archivo UMD, luego transfiera los archivos umd a archivos xyz para facilitar la visualización en varios otros paquetes, como VMD o VESTA.Invierta el archivo UMD en archivos POSCAR de tipo VASP utilizando umd2poscar. py script, seleccionando instantáneas de las simulaciones con la frecuencia predefinida. Ejecute el gofrs_umd.py script para calcular la función de distribución de pares para todos los pares de tipos atómicos A y B.La salida se escribe en una pestaña de archivo ASCII separada con la extensión gofrs.dat. Extraer las distancias medias de enlace interatómico como los radios de las primeras esferas de coordinación. Para ello, identifique la posición del primer máximo de las funciones de distribución de pares trazando los gofrs.dat en una aplicación de hoja de cálculo y buscando los máximos y mínimos para cada par de átomos, luego identifique el radio de la primera esfera de coordinación como el primer mínimo del PDF utilizando el software de hoja de cálculo. Ejecute el script de especiación para obtener la matriz de conectividad y obtener los poliedros de coordinación o la polimerización. Ejecute el speciation_umd.py script con el indicador r0, que muestrea el gráfico de conectividad en el primer nivel para identificar los poliedros de coordinación. Ejecute el speciation_umd. py script con el flag r1, que muestrea el gráfico de conectividad en todos los niveles de profundidad para obtener la polimerización.Trazar la vida útil de cada cúmulo atómico de todas las especies químicas encontradas en la simulación tal como se encuentran en la pápula. archivos dat. Extraiga los desplazamientos cuadráticos medios o MSD de los átomos en función del tiempo para obtener la autodifusividad, luego calcule el MSD utilizando la serie de msd_umd.py scripts y calcular el MSD promedio de cada tipo atómico. Calcular el MSD de cada átomo y de la especie química. Traza el MSD usando un software basado en hojas de cálculo y calcula los coeficientes de difusión a partir de la pendiente del MSD.Ejecute el vibr_spectrum_umd. py script para calcular la autocorrelación de velocidad atómica o función VAC para cada tipo atómico y realizar su rápida transformada de Fourier. Traza el espectro vibratorio a partir del vibr.dat utilizando un software similar a una hoja de cálculo. Identificar el valor finito en omega es igual a cero que corresponde al carácter difusivo del fluido en los diversos picos del espectro a frecuencia finita. Ejecutar promedios.py para extraer los valores promedio y la dispersión de presión, temperatura, densidad y energía interna de los archivos UMD. Finalmente, ejecute los promedios completos. py script para realizar el análisis estadístico completo incluyendo el error de la media.La pirolita es un modelo de fusión de silicato multicomponente que se aproxima mejor a la composición de la Tierra de silicato a granel. El paquete UMD se utilizó para extraer varios rasgos característicos de la pirolita fundida. El máximo de la función de distribución del par silicio-oxígeno se encuentra en 1.635 angstrom, que es la mejor aproximación a la longitud de la curvatura.Usando este límite como la distancia de enlace silicio-oxígeno, el análisis de especiación muestra que las unidades de ortosilicato que pueden durar hasta unos pocos picosegundos dominan la fusión. Hay una parte importante de la masa fundida que muestra una polimerización parcial reflejada por la presencia de dímeros como el disilicato y trímeros como las unidades Si3Ox. Su vida útil correspondiente es del orden del picosegundo.Todos los polímeros de orden superior tienen una vida útil considerablemente más corta. Los diferentes valores de los pasos verticales y horizontales producen varios muestreos del MSD. Incluso los valores grandes de Z y V son suficientes para definir las pendientes y, por lo tanto, los coeficientes de difusión de los diferentes átomos.El tiempo de post-procesamiento aumenta dramáticamente para grandes valores de Z y V.Finalmente, las funciones de autocorrelación de velocidad atómica producen el espectro vibratorio de la masa fundida. Aquí se muestran las contribuciones de los átomos de magnesio, silicio y oxígeno, así como el valor total. Al intentar este protocolo, compruebe siempre la convergencia.Asegúrese de que las trayectorias atómicas sean lo suficientemente largas como para capturar adecuadamente el fenómeno que le interesa. Esta técnica cubre el post-procesamiento de los resultados de la simulación. Las simulaciones y su análisis deben hacerse en paralelo.

Research

JoVE Journal

Chemistry

Isolation and Chemical Characterization of Lipid A from Gram-negative Bacteria

Aislamiento y Caracterización química del lípido A de bacterias gram-negativas

Lipopolysaccharide (LPS) is the major cell surface molecule of gram-negative bacteria, deposited on the outer leaflet of the outer membrane bilayer. LPS ...

Investigación

Química

Analyzing Protein Dynamics Using Hydrogen Exchange Mass Spectrometry

El análisis de dinámica de proteínas utilizando hidrógeno Intercambio de Espectrometría de Masas

All cellular processes depend on the functionality of proteins. Although the functionality of a given protein is the direct consequence of its unique ...

Visualizing Lignification Dynamics in Plants with Click Chemistry: Dual Labeling is BLISS!

Visualización dinámica de lignificación en plantas con la química Click: doble etiquetado es felicidad!

Lignin is one of the most prevalent biopolymers on the planet and a major component of lignocellulosic biomass. This phenolic polymer plays a vital ...

Combustion Chemistry of Fuels: Quantitative Speciation Data Obtained from an Atmospheric High-temperature Flow Reactor with Coupled Molecular-beam Mass Spectrometer

Química de la combustión de combustibles: especiación cuantitativa datos obtenidos de un Reactor de flujo de alta temperatura atmosférica con junto haces moleculares de espectrómetro de masas

This manuscript describes a high-temperature flow reactor experiment coupled to the powerful molecular beam mass spectrometry (MBMS) technique. This ...

Probing the Structure and Dynamics of Interfacial Water with Scanning Tunneling Microscopy and Spectroscopy

Sondear la estructura y dinámica del agua superficial con túnel microscopía de barrido y espectroscopia

Water/solid interfaces are ubiquitous and play a key role in many environmental, biophysical, and technological processes. Resolving the internal ...

A Closed-Type Wireless Nanopore Electrode for Analyzing Single Nanoparticles

Un electrodo Nanopore inalámbricas cerradas para el análisis de las nanopartículas individuales

Measuring the intrinsic features of single nanoparticles by nanoelectrochemistry holds deep fundamental importance and has potential impacts in ...

Multiscale Sampling of a Heterogeneous Water/Metal Catalyst Interface using Density Functional Theory and Force-Field Molecular Dynamics

Muestreo multiescala de una heterogénea agua Metal catalizador interfaz utiliza teoría del funcional densidad y campo de fuerza dinámica Molecular

A significant number of heterogeneously-catalyzed chemical processes occur under liquid conditions, but simulating catalyst function under such conditions ...

Vibrational Spectra of a N719-Chromophore/Titania Interface from Empirical-Potential Molecular-Dynamics Simulation, Solvated by a Room Temperature Ionic Liquid

Espectros vibratorios de una interfaz N719-Chromophore/Titania de la simulación molecular-dinámica empírica-potencial, Solvada por un líquido iónico de temperatura ambiente

The accurate molecular simulation prediction of vibrational spectra, and other structural, energetic and spectral characteristics, of photo-active ...

Computation of Atmospheric Concentrations of Molecular Clusters from ab initio Thermochemistry

Computation of Atmospheric Concentrations of Molecular Clusters from ab initio Thermochemistry

Cálculo de las concentraciones atmosféricas de los cúmulos moleculares a partir de la termoquímica ab initio

The computational study of the formation and growth of atmospheric aerosols requires an accurate Gibbs free energy surface, which can be obtained from gas ...

Interactive Molecular Model Assembly with 3D Printing

Ensamblaje de modelo molecular interactivo con impresión 3D

With the growth in accessibility of 3D printing, there has been a growing application of and interest in additive manufacturing processes in chemical ...