Summary
El objetivo del protocolo presentado aquí es generar y trayectorias de configuraciones de moléculas de agua líquida en especie catalítica sobre una superficie plana de metal de transición de la muestra. Las configuraciones de muestras pueden utilizarse como estructuras a partir de métodos basados en la mecánica cuántica.
Abstract
Un número significativo de procesos químicos catalizados heterogéneo se producen bajo condiciones de líquido, pero que simula la función de catalizador en estas condiciones es difícil cuando es necesario incluir las moléculas solventes. La ruptura del vínculo y la formación de procesos modelados en estos sistemas requieren el uso de métodos de química cuántica. Puesto que las moléculas en la fase líquida están en constante movimiento térmico, simulaciones también deben incluir muestreo configuracional. Esto significa que múltiples configuraciones de las moléculas de líquido deben ser simuladas para cada especie catalítica de interés. El objetivo del protocolo presentado aquí es generar y trayectorias de configuraciones de moléculas de agua líquida alrededor de especies catalíticas en superficies planas de metales de transición en una forma que equilibra la exactitud Química computacional gasto de la muestra. En concreto, se utilizan simulaciones de dinámica molecular (FFMD) del campo de fuerza para generar configuraciones de moléculas de líquido que posteriormente pueden ser empleadas en métodos basados en la mecánica cuántica como teoría de densidad funcional o ab initio molecular dinámica. Para ilustrar esto, en este manuscrito, el protocolo se utiliza para productos intermedios catalizadoras que podrían estar implicados en el camino para la descomposición del glicerol (C3H8O3). Las estructuras que se generan utilizando FFMD se modelan en DFT para estimar entalpías de solvatación de la especie catalítica e identificar cómo participan moléculas de H2O en descomposición catalítica.
Introduction
Modelado de fenómenos moleculares implicados en la catálisis heterogénea bajo condiciones de líquido es necesario para la función catalítica de la comprensión; sin embargo, esto sigue siendo desafiante porque requiere un fino equilibrio entre la precisión de químico y costo computacional. En general, desde catálisis implica la ruptura y formación de enlaces químicos, la mecánica cuántica debe utilizarse al menos en cierto grado; sin embargo, simulaciones de largo están desafiando en mecánicos del quántum, que requieren recursos informáticos importantes. Puesto que las moléculas en la fase líquida están en constante movimiento térmico, simulaciones también deben incluir muestreo configuracional, es decir, deben incorporar múltiples disposiciones espaciales de las moléculas del líquido, como cada arreglo espacial diferente (es decir, cada uno configuración) tiene una energía diferente. Esto significa que múltiples configuraciones de las moléculas de líquido deben ser simuladas para cada especie catalítica de interés. Estas necesidades – para usar la mecánica cuántica y para realizar múltiples cálculos por especie catalítica – pueden hacer modelado en catálisis heterogénea en fase líquida computacionalmente intratable. El propósito del método descrito en este documento es permitir computacionalmente manejables simulaciones de fenómenos de catálisis heterogénea en fase líquida.
Nos interesa particularmente heterogéneo catalizadas reacciones que se llevan a cabo bajo el agua líquido. Las moléculas de agua tienen una influencia significativa en los fenómenos catalíticos, como interactuar con la especie catalítica (por ejemplo, a través de las fuerzas de dispersión e hidrógeno)1,2,3,4,5 ,6,7,8,9,10,11,12,13,14 , 15 , 16 , 17 , 18 , 19 , 20 , 21 , 22 , 23, participan en reacciones catalíticas1,7,8,9,15,21,22,24 ,25,26,27y que influyen en las vías de reacción o tasa catalítica1,11,12,15, 18,23,25,27,28,29,30,31. Modelado de estos fenómenos ha llevado a cabo utilizando QM o ab initio dinámica molecular (AIMD)1,2,6,7,14,22 ,25,27,28,32,33,34, fuerza dinámica molecular de campo (FFMD)35 y la mecánica cuántica/molecular mecánica (QM/MM)10. En AIMD y FFMD, se mueven los átomos en el sistema en virtud de las ecuaciones de Newton de movimiento según las fuerzas que actúan sobre ellos. En AIMD, las fuerzas y energía del sistema se calculan con la mecánica cuántica, mientras que en FFMD, la energía del sistema y las fuerzas se calculan utilizando datos QM o campos, que son expresiones algebraicas que son parametrizadas en basan experimental de la fuerza. En QM/MM, se calcula la parte del sistema donde se produce el vínculo rompiendo y formando con QM, y el resto del sistema se calcula con m, que emplea campos de fuerza. Porque emplean directamente QM, AIMD y QM/MM son más apropiados para la captura de la ruptura del vínculo y formación se produce en catálisis heterogénea de la fase acuosa; sin embargo, FFMD es significativamente más cómputo manejable y por lo tanto más adecuado para la generación de las configuraciones de moléculas de líquido H2O. El método presentado en este protocolo equilibrio químico exactitud y costo computacional empleando una combinación de QM y FFMD.
Específicamente, este método utiliza simulaciones de FFMD para generar configuraciones de líquido H2O y QM para calcular energías de sistema. FFMD se lleva a cabo utilizando LAMMPS. 36 los campos de fuerza utilizados en FFMD en este trabajo emplean Lennard-Jones + potencial de Coulomb (LJ + C), donde los parámetros LJ se han tomado desde el TIP3P/CHARMM modelo37 de H2O, el campo de fuerza universal38 (UFF) para el Pt y el Campo de fuerza de OPLS AA39 especies catalíticas y los parámetros de Coulomb se han adoptado desde el TIP3P/CHARMM37 modelo H2O y el campo de fuerza de OPLS AA39 especies catalíticas. Se han creado los parámetros Coulomb para átomos de Pt a 0. Cálculos de QM se realizan utilizando la VASP código40,41,42, que es un código (DFT) la teoría del funcional de densidad. Inserciones de la molécula de agua se realizan con un código desarrollado in-House llamada Monte Carlo plug-in para métodos de cuántica (MCPliQ). Conversiones de archivo de VASP a LAMMPS en este protocolo se realizan con el software de dinámica Molecular Visual (VMD)43.
El protocolo pretende generar configuraciones de moléculas de agua líquida alrededor de especies catalíticas en superficies de metal de transición plana en baja cobertura. Cobertura es θ denota y se define como el número de adsorbates por átomo metálico superficial (es decir, el número de superficie adsorbates normalizado por el número de átomos del metal en la capa superior de la placa metálica en el modelo de catalizador). En este manuscrito, baja cobertura se define como θ ≤ 1/9 monocapa (ML), donde 1 ML significa una especie catalítica por átomo metálico superficial. Los modelos de catalizador deben colocarse en cajas de simulación periódica. Las cajas de simulación debe ser cubos. Este manuscrito muestra el uso del protocolo para la generación de configuraciones de líquido H2O que pueden utilizarse para calcular cantidades de interés en catálisis heterogénea de la fase acuosa.
Este protocolo requiere que el usuario tiene acceso a versiones instaladas y funcionando de la VASP, MCPliQ, LAMMPS y VMD. Más información sobre VMD (https://www.ks.uiuc.edu/Research/vmd/), VASP (https://www.vasp.at/) y LAMMPS (https://Lammps.sandia.gov/) están disponibles en sus sitios Web. El software MCPliQ está documentado en https://github.com/getman-research-group/JoVE_article, junto con todos los archivos y scripts de Python, mencionados en el presente Protocolo. Este protocolo asume que los archivos ejecutables y scripts mencionados dentro de se ejecutará en un equipo de investigación de alto rendimiento y se instalan en un directorio que se encuentra en la variable $PATH del usuario. Si un ejecutable o una secuencia de comandos se coloca en una ubicación que no está en el usuario de $PATH, entonces debe incluirse la ruta al ejecutable para ejecutarlo. Ejecutables y scripts se ejecutan en pasos 2.1.2 2.2.1, 2.2.8, 3.1, 4.2, 5.2 y 6.1.2. Por ejemplo ejecutar el código MCPliQ en el paso 2.1.2 desde un directorio que no está en el usuario de $PATH, el usuario teclea $PATHTOMCPLIQ/mcpliq en la interfaz de línea de comandos en lugar de mcpliq, donde $PATHTOMCPLIQ es el lugar donde la mcpliq se ha almacenado el archivo ejecutable (por ejemplo, podría ser $PATHTOMCPLIQ ~ / bin). Antes de iniciar este protocolo, todos los archivos ejecutables y scripts se conceda permisos de ejecutables (por ejemplo, en Linux, esto puede hacerse escribiendo chmod + x mcpliq en la interfaz de línea de comandos desde el directorio donde se guarda el ejecutable mcpliq). Además, se deben cargar los módulos requeridos por cualquiera de los software o scripts (estas dependencias serán específicas para instalaciones individuales de software y el equipo donde se ejecutarán las simulaciones).
Protocol
1. generar la estructura de adsorbato
- Crear un archivo de VASP POSCAR que comprende un supercell con condiciones de límite periódicas, como lo haría si estaban realizando simulaciones de adsorbatos en superficies de metal bajo vacío. El supercell debe incluir una suposición inicial de la estructura del adsorbato y la superficie de metal como espacio vacío por encima del adsorbate para agregar moléculas de H2O. Más detalles sobre cómo hacer esto se encuentran en anterior trabajo35,44.
Nota: Es importante que la altura del espacio vacío por lo menos 12 Å arriba del adsorbato. - La estructura de relajarse y reducir al mínimo su energía utilizando el código de la VASP. Más detalles sobre cómo hacer esto se encuentran en anterior trabajo35,44. Esto producirá un archivo llamado CONTCAR, que se utilizará en la siguiente sección.
2. Añadir expresas moléculas de2O de H
- Añadir N moléculas de agua al espacio vacío en la CONTCAR creado en el paso 1.2 usando el código MCPliQ, donde N = ρV, ρ es la densidad del agua, y V es el volumen del espacio vacío por encima del adsorbate.
Nota: se deben tomar ρ la densidad del agua determinada por el modelo de TIP3P/CHARMM agua a la temperatura de la simulación. V será refinado en el paso siguiente.- Especifique la siguiente información en el archivo de master_input.txt de MCPliQ: el número de moléculas de H2O para agregar (N) al cambiar el primer argumento en la línea 28, la ruta para el archivo water.txt cambiando el segundo argumento en la línea 28 y el mínimo y altura máxima de la supercélula que puede ser ocupado por las moléculas de agua cambiando la coordenada z de mínimo en la línea 11 y la máxima coordenada z en la línea 12.
- Ejecutar el código de MCPliQ escribiendo mcpliq desde la interfaz de línea de comandos para introducir las moléculas de agua en el archivo CONTCAR. El código de salida uno o más archivos con la extensión de archivo. POSCAR.
Nota: Si se producen más de un archivos POSCAR, ellos serán nombrados POSCAR_n.POSCAR. Seleccione el archivo donde n es el más grande.
- Generar los archivos de entrada LAMMPS para una simulación de NPT y equilibrar el volumen de la célula usando FFMD en el conjunto de la NPT en LAMMPS.
- Ejecutar el script lmps_bond_angle.py en el. Archivo POSCAR generado en el paso 2.1.2 escribiendo lmps_bond_angle.py $filename. POSCAR en la interfaz de línea de comandos, donde $filename es el nombre del. Archivo POSCAR generado en el paso 2.1.2. Este script crea un archivo llamado $filename. POSCAR.bond_angle_info.txt lista de bonos y los ángulos que se utilizarán en el archivo de datos LAMMPS.
- Abra VMD y seleccione archivo > Nueva molécula en la ventana principal para abrir la ventana de Explorador de archivos de la molécula . Seleccione VASP_POSCAR en el menú desplegable de Determinar el tipo de archivo . Haga clic en examinar y desplácese hasta el $filename. Archivo de la POSCAR. Haga clic en carga para abrir la $filename. Archivo de la POSCAR.
- Abra la consola de Tk en VMD seleccionando extensiones > Tk consola de la ventana principal de VMD.
- Ejecute el siguiente comando en la consola de Tk: topo writelammpsdata $WDPATH/data.myadsorbate completo, donde $WDPATH es el directorio en el equipo donde VMD escribirá el archivo de datos LAMMPS y data.myadsorbate es el nombre del archivo de datos LAMMPS.
- Eliminar la sección de enlaces y ángulos en la parte inferior del archivo data.myadsorbate. A continuación, anexar las listas de enlace y ángulo en el archivo $filename. POSCAR.bond_angle_info.text en data.myadsorbate.
Nota: Los índices para el tipo de enlace O-H y tipo del ángulo H-O-H de las moléculas de agua en el $filename. Ambos POSCAR.bond_angle_info.txt archivo se establece en 1. Así, tipos de enlace y ángulo para adsorbatos deben empezar a contar en 2. - Edite el archivo data.myadsorbate añadiendo los parámetros de Lennard-Jones a la sección de Coeffs par y los parámetros de Coulomb a la sección de átomos . Lennard-Jones y parámetros de Coulomb para las moléculas de H2O átomos del adsorbato y átomos metálicos de superficie hay que agregar.
Nota: Los parámetros de Lennard-Jones para Pt átomos, las moléculas de agua y átomos del adsorbato en el presente Protocolo se obtienen de la UFF38TIP3P/CHARMM37y OPLS AA39 campos de fuerza, respectivamente. Parámetros de Coulomb para los átomos de las moléculas y adsorbatos de agua se obtienen de TIP3P/CHARMM37 los OPLS AA39 campos de fuerza, respectivamente. Parámetros de Coulomb para átomos de Pt estén configurados en 0 en este protocolo. Alternativamente, podrían utilizarse cargas parciales calculados para los parámetros de Coulomb para el adsorbato átomos y átomos de Pt. - Copie el LAMMPS archivo de entrada input.equil en el directorio $WDPATH. Editar la variable de grupo en la línea 34 para indicar los índices de tipo de átomo agua oxígeno y átomos de hidrógeno del agua y la variable de grupo en la línea 35 para indicar los índices de tipo de átomo para los átomos de Pt y adsorbato.
- Ejecutar el software LAMMPS escribiendo mpiexec np - XX lmp_mpi < input.equil en la interfaz de línea de comandos, donde XX es el número de núcleos de CPU a usar, y lmp_mpi es el nombre del ejecutable LAMMPS. Haciendo esto se ejecutará una minimización de la energía para refinar la configuración de H2O, seguida de una simulación de FFMD en número constante de moléculas de H2O (N), volumen (V) y temperatura (T) para el agua a la temperatura de la simulación, seguida por una simulación de FFMD corren constante N, presión (P) y la temperatura (T) para determinar la altura físicamente correcta de la caja de simulación. Archivos de salida que se utilizará en la sección 3 se llaman data.myadsorbate_npt y log.myadsorbate.
Nota: La duración de la simulación de NPT debe tiempo suficiente para formar un "equilibrado" ejecutar, donde el volumen de la supercélula trata de estado estable y una «producción» ejecutar, que se utiliza para la muestra que el conjunto de promedios (aquí, la altura del supercell). Durante el funcionamiento equilibrado, debe Nivele el volumen de supercell cuando graficados contra el tiempo a un valor de estado estacionario. Una vez que esto ocurre, la simulación de la NPT puede decirse que en su funcionamiento de producción. Comprobar equilibrado de la simulación de NPT por asegurar que las fluctuaciones en la altura de supercell (lz) son mínimas o se han acercado a un valor constante. Si se producen grandes fluctuaciones, luego volver a generar una configuración de2O H disminuyendo las incógnitas en la línea 92 en el archivo input.equil y repetición paso 2.2.8 o volver a empezar desde el paso 2.1.1.
3. extraer la altura apropiada de la supercélula
- Ejecutar el script get_npt_lz.py en el archivo log.myadsorbate si log.myadsorbate de get_npt_lz.py en la interfaz de línea de comandos. Este script genera la altura promedio de supercell de la porción de "producción" de la simulación de la NPT en el archivo avg_lz.txt.
Nota: El script de get_npt_lz.py asume que LAMMPS escribe la longitud de la célula dimensión z (lz) en el archivo log.myadsorbate cada 1000 fs (personalizables en línea 20 de la escritura de la get_npt_lz.py), que es el valor predeterminado en el archivo de entrada input.equil proporcionado LAMMPS. El script de get_npt_lz.py detecta y descarta los primeros 2 valor ns (personalizables en línea 19 de la secuencia de comandos de get_npt_lz.py) de lz valores en el archivo log.myadsorbate, que comprenden la porción del equilibrado de la simulación, mientras que los restantes 3 ns comprenden el " parte de la producción"y se utilizan así el script get_npt_lz.py para calcular la longitud promedio z-dimensión. Además del archivo avg_lz.txt, el script de get_npt_lz.py salidas de un archivo llamado npt_data.txt, que proporciona valores de lz en función de los rebotes, así como un archivo llamado npt_plot.png, que los mismos datos de parcelas. La trama puede utilizarse para comprobar el equilibrado de la simulación de la NPT. - Reconstruir el supercell utilizando la altura media en NPT.
- Copie el archivo data.myadsorbate_npt en un nuevo directorio, que se refiere aquí como $WD2PATH y cámbiele el nombre data.myadsorbate.
- Edite el nuevo archivo de data.myadsorbate para que la altura de lz es igual a la salida de valor medio de la escritura de la get_npt_lz.py cambiando los argumentos zlo y zhi en el archivo data.myadsorbate zlo es 0.0 y zhi es el valor de lz desde el archivo avg_lz.txt en s TEP 3.1.
4. generar configuraciones de moléculas de H2O
- Copie el LAMMPS archivo de entrada input.prod en $WD2PATH. Editar la variable de grupo en la línea 32 para indicar los índices de tipo de átomo agua oxígeno y átomos de hidrógeno del agua y la variable de grupo en la línea 33 para indicar los índices de tipo de átomo para los átomos de Pt y adsorbato.
- Ejecutar el software LAMMPS escribiendo mpiexec np - XX lmp_mpi < input.prod en la interfaz de línea de comandos, donde XX es el número de núcleos de CPU a usar, y lmp_mpi es el nombre del ejecutable LAMMPS. Haciendo esto se ejecutará una constante simulación de NVT en las moléculas de H2O. La clave de la salida de esta simulación es el archivo dump.myadsorbate.lammpstrj.
Nota: La duración de la simulación de NVT debe tiempo suficiente para abarcar un equilibrado ejecutar, donde la energía del sistema viene a constante del estado, y una funcionamiento de producción, de la cual el conjunto promedio (aquí, las posiciones espaciales de las moléculas de agua) son muestreadas. Durante el funcionamiento equilibrado, la energía del sistema cuando se enfrenta a tiempo debe nivele a un valor de estado estacionario. Una vez que esto ocurre, la simulación de NVT puede decirse que en su funcionamiento de producción.
5. determinar la duración del enlace del hidrógeno para el muestreo de tiempo
- Editar el script hb_lifetime_dist.py para especificar: las incógnitas del primer fotograma del archivo dump.myadsorbate.lammpstrj cambiando la variable ///ActualStart en la línea 22, con qué frecuencia marcos se escriben en el archivo de trayectoria LAMMPS cambiando la variable de incógnitas en línea 23, debería considerar la primera y la última timesteps el script (es decir, la porción de la producción de la trayectoria) cambiando las variables N_first y N_last en las líneas 24 y 25, si fotogramas consecutivos se consideran o se saltan frames cambiando el todos variable en la línea 26 y el número de líneas por el marco de la sección del trayectoria archivo cambiando la variable del frameLine en línea 27. Además, editar líneas del 31 al 35 para especificar que tipos de átomo dentro del archivo data.myadsorbate pertenecen al adsorbato y que tipos de átomo pertenecen a las moléculas de H2O.
Nota: El script de hb_lifetime_dist.py analiza las configuraciones de H2O en el plazo de producción y determina si cualquier moléculas de H2O hidrógeno enlazado al adsorbato. Entonces cuenta el tiempo de simulación que cada hidrógeno bonos sigue intacta y esta información como una distribución de tiempos de vida del enlace de hidrógeno en unidades de ps. La versión específica de la secuencia de comandos que se proporciona con este protocolo asume que LAMMPS escribe la configuración de moléculas de H2O en el archivo dump.myadsorbate.lammpstrj cada 1000 fs, que es el valor predeterminado en el archivo de entrada input.prod proporcionado LAMMPS. Detecta y descarta los primeros 2 valor ns de configuraciones en el archivo el dump.myadsorbate.lammpstrj, que comprenden la porción del equilibrado de la simulación y ns utiliza los 3 restantes para calcular la vida del enlace de hidrógeno. - Ejecutar el script hb_lifetime_dist.py en el archivo de dump.myadsorbate.lammpstrj escribiendo hb_lifetime_dist.py en la interfaz de línea de comandos. Hacer esto producirá un archivo llamado distribution_HB_lifetime.dat.
- Representar los datos en el archivo distribution_HB_lifetime.dat para ver la distribución de hidrógeno vidas que ocurrió durante la simulación de NVT.
- Determinar el incremento del tiempo a utilizar para el intervalo de muestreo basado en la vida de hidrógeno calculado. La mejor opción es la vida útil máxima de hidrógeno; Alternativamente, puede utilizarse un valor que captaría el intervalo de confianza de 95%.
6. configuraciones de ejemplo de moléculas de2O de líquido H
- Determinar el número de configuraciones de la producción de la trayectoria de NVT FFMD para cálculos posteriores. El número de configuraciones debe seleccionarse para que el tiempo mínimo entre configuraciones es igual o mayor que el intervalo de muestreo identificado en la sección 5.
- Editar el valor predeterminado para la variable num_frames en la línea 21 de la lammps_frames.py secuencia de comandos para especificar el número de configuraciones para extraer.
- Ejecutar el script lammps_frames.py en el dump.myadsorbate.lammpstrj archivo escribiendo lammps_frames.py en la interfaz de línea de comandos. Esto generará una lista de simulación de tiempos correspondientes a las configuraciones que deben ser extraídas del archivo dump.myadsorbate.lammpstrj. Estas configuraciones se pueden utilizar como a partir de estructuras AIMD o QM simulaciones.
Nota: 1) la secuencia de comandos de lammps_frames.py automáticamente detecta el registro LAMMPS y volcado de archivos así como la porción de la producción de la trayectoria dentro del archivo de volcado y divide el número de configuraciones en el archivo de volcado en 10 grupos. Alternativamente, el usuario puede especificar el archivo de registro, el archivo de volcado y el número de configuraciones de la interfaz de línea de comandos utilizando la -l -d y opciones - n, respectivamente. Para ello, el usuario debe escribir lammps_frames.py - n XX -l $logfilename -d $dumpfilename en la interfaz de línea de comandos, donde XX es el número de configuraciones, $logfilename es el nombre del archivo de registro LAMMPS y $dumpfilename es el nombre de la Archivo de la trayectoria (descarga) de LAMMPS. Los tiempos de simulación que son salida se refieren a los tiempos promedio en cada grupo. 2) si las configuraciones se calculará en VASP con la bandera LDIPOLE encendida, debe añadirse una capa pequeña de espacio vacío en la parte superior de la supercélula sobre la capa de agua. Esto facilitará la convergencia de la estructura electrónica en el cálculo de la VASP. Agregar un adicional 3 Å de espacio vacío por encima de las moléculas de H2O ha tenido éxito en las simulaciones que se examinan a continuación.
Representative Results
Un uso de este protocolo es calcular las energías de interacción entre el agua líquida y especies catalíticas, es decir, ΔEint35:
∆Eint=Eespecies catalíticas + H2O+Esuperficie del catalizador limpio-Eespecie catalítica-Ecatalizador limpio superficie + H2O
Eespecies catalíticas + H2Oes la energía de una configuración de moléculas de H2O alrededor de una especie catalítica en una superficie de metal, Esuperficie del catalizador limpia es la energía de la superficie del catalizador limpio en vacío, E Especies catalíticas es la energía de la especie catalítica en una superficie de metal en el vacío, y Elimpiar catalizador superficie + H2O es la energía de la configuración de H2O sobre la superficie del catalizador con la especie catalítica quitada. Las posiciones de las moléculas de H2O para calcular Eespecies catalíticas + H2O y Elimpiar catalizador superficie + H2O deben ser idénticos. Todos los valores de E se calculan utilizando el código de la VASP. La cantidad ΔEint incluye todas las interacciones físicas y químicas entre todas las moléculas en la estructura del agua líquida y la especie catalítica y da una estimación razonable de la entalpia de solvatación de la catalítica especie, que es necesario para calcular la energía libre de solvatación y total energía libre. Tabla 1 proporciona valores para ΔEint calculado para las especies en un catalizador de Pt(111) la superficie con fórmulas químicas iguales a CxHyOz en unidades de eV (1 eV = 96.485 kJ/mol). Los valores se calcularon en coberturas ≤1/9 ML.35,46 los valores reportados son los promedios tomados en 10 configuraciones de líquido H2O, y las incertidumbres se divulgan como desviaciones estándar. Todos los valores son negativos, indicando las interacciones favorables con agua.
Otra aplicación de este protocolo es generar estructuras partidas para AIMD. Película 1 es una película de una trayectoria AIMD que se inició desde una configuración generada por este protocolo. Al comienzo de esta película, un adsorbato COH se muestra sobre una superficie de Pt(111) bajo una estructura de líquido H2O. Molécula de un H2O se acentúa, que formó un enlace del hidrógeno con COH. En el transcurso de la película, esta molécula de H2O abstrae el protón de adsorbato COH y deposita un segundo átomo de hidrógeno en la superficie de Pt(111). La molécula de H2O, por lo tanto ayuda a catalizar la reacción COH * + * → CO * + H *, donde el * s indican sitios catalíticos. Esta simulación pone de relieve la fuerza principal y el propósito principal del método de muestreo multiescala descrito. Numerosas configuraciones de moléculas de H2O son generadas con FFMD, debido a su fuerza en tratabilidad computacional. Sin embargo, una limitación de FFMD es que no puede capturar bonos ruptura y formando a menos que se aplica un campo de fuerza reactivo. AIMD utiliza la mecánica cuántica para calcular las energías y por lo tanto puede capturar bonos rompiendo y formando. Sin embargo, AIMD está exigiendo demasiado de cómputo generar todas las configuraciones de las moléculas de H2O necesarias para muestreo suficiente se ha logrado. Así, este protocolo combina los dos métodos.
Las estructuras de moléculas de O2de H líquidas generadas por este procedimiento dependen de la configuración de entrada. Ajuste estas mal puede tener influencias no deseadas en las estructuras del agua. Por ejemplo, cuando las distancias intermoleculares son muy pequeñas o cuando otros parámetros en los archivos de entrada dinámica molecular son mal ajustados o asumir valores incontrolados, la estructura del agua puede llegar a ser irracional. En estas circunstancias, la estructura del agua será "blow-up" unintendedly durante la trayectoria FFMD. La figura 1 muestra un ejemplo de esto. La instantánea en el lado izquierdo es la estructura de partida de una carrera FFMD, y la instantánea a la derecha es una instantánea tomada en ps 1 de comenzar la simulación. Como puede verse, las moléculas de H2O han movido lejos de la superficie. Esto es causado por la configuración incorrecta en los archivos de entrada de simulación y no es una estructura que es probable que ocurra en la realidad.
Figura 1: Ejemplo de un resultado negativo. La simulación de dinámica molecular de campo de fuerza "saltar" debido a un ajuste incontrolado o valor. Imagen izquierda: la geometría a partir de la superficie Pt(111), el adsorbato y la estructura del agua líquida. Imagen derecha: la geometría de la superficie de Pt(111), el adsorbato y el agua líquida estructura ps menos 1 más adelante. En la imagen derecha, las moléculas de H2O se separan de la superficie debido a las fuerzas unphysically grandes. Haga clic aquí para ver una versión más grande de esta figura.
Película 1: Ab initio simulación dinámica molecular (AIMD) iniciada de una configuración en multiescala muestreo. Una molécula de H2O que originalmente hidrógeno enlazado a un adsorbato COH en una superficie de Pt(111) abstrae el protón de COH y deposita un segundo hidrógeno sobre la superficie Pt(111). Este bono rompiendo y formando el evento puede ser capturado por AIMD pero no con el campo de fuerza dinámica molecular (FFMD) salvo que se utilice un campo de fuerza reactivo. La configuración inicial de las moléculas de2O H utilizada en esta simulación AIMD fue generada usando FFMD como se describe en este manuscrito. Por favor haga clic aquí para ver este video. (Clic derecho para descargar)
| Especies catalíticas | ∆Eint (eV) |
| COH | -0.70 ± 0.07 |
| CO | -0,03 ± 0.03 |
| CH2OH | -0.64 ± 0,12 |
| CHO-CHOH-CH2OH | -0.93 ± 0.22 |
| COH-COH-CH2OH | -0.87 ± 0.23 |
| COH-CHOH-COH | -1.72 ± 0.26 |
| CHOH-COH-CO | -1.57 ± 0.25 |
| CHO-CO-CO | -0.31 ± 0.19 |
Tabla 1: Resultados de interacción de la energía de la especie catalítica de agua. Las energías de interacción en eV calculan para ocho CxHyOz adsorbatos en Pt(111). Valores registrados se toman los promedios sobre múltiples configuraciones de líquido H2O. Las incertidumbres son las desviaciones estándar de los promedios. 1 eV = 96.485 kJ/mol.
Discussion
El método presentado fue seleccionado por su facilidad de aplicación, pero podrían realizarse múltiples personalizaciones. Uno, pueden modificar los campos de fuerza utilizados en las simulaciones FFMD. Cambiar los parámetros de campo de fuerza o potenciales puede hacerse editando los archivos de entrada y datos LAMMPS. Del mismo modo, podrían emplearse solventes distintos de H2O. Para hacer esta modificación, la molécula de disolvente deseada tendría que insertarse a partir de paso 2.1.1, y los archivos de entrada LAMMPS tendrían que editarse para incorporar los potenciales apropiados y los parámetros. Insertar la nueva molécula solvente también requeriría suministrando las coordenadas internas de la molécula del solvente en un archivo .txt análogo al archivo water.txt.
Otra modificación que pudiera hacerse es modificar la zona de la losa superficial. Los resultados discutidos en este manuscrito emplean losas superficiales 3 Pt x 3 Pt o Pt 4 x 4 Pt, que tienen áreas de superficie inferior a 120 Å2. A medida que aumenta la superficie de la losa, también aumenta el costo computacional. Costo computacional tiene el impacto más grande en la sección 5 del presente Protocolo. Si los pasos de procesamiento de datos en la sección 5 es computacionalmente prohibitivos, grandes datos post procesamiento estrategias tales como ésos discutidos en Li et al 201845 pueden ser empleados.
Posibles fuentes de incertidumbre para este procedimiento incluyen el campo de fuerza empleado, el método de muestreo y la frecuencia de muestreo. La estructura del agua está determinada por el campo de fuerza que se utiliza, lo que significa que la elección del campo de fuerza podría influir en las configuraciones específicas de moléculas de H2O. Nuestro grupo ha evaluado cómo influyen las energías de interacción calculados en FFMD y encuentra que la elección del campo de fuerza contribuye menos de 0,1 eV a esta energía de interacción en las opciones de campo de fuerza para moléculas de H2O y los átomos de Pt. Otra fuente de incertidumbre es el método de muestreo, que influye en las configuraciones específicas que se utilizan para calcular la cantidad de interés. Nuestro grupo ha comparado el rendimiento del "tiempo" método presentado en este protocolo con un método de "muestreo de energía", que está sesgado a configuraciones de energía más baja de moléculas de H2O, en la interacción energías calculan en DFT y encontraron a ambos de estos métodos de muestreo estadísticamente igual dan valores de35,46. La frecuencia de muestreo también puede influir en los resultados. Hemos evaluado cómo aumentar el número de configuraciones de 10 a 30.000 influye en las energías de interacción promedio calculadas en FFMD para 40 diferentes C3HxO3 adsorbatos y encontró que la frecuencia de muestreo contribuye menos de 0,1 eV a la energía de interacción promedio44.
La principal limitación de este método es que los adsorbates se aproximan por las estructuras bajo vacío durante las simulaciones de FFMD. En realidad, los adsorbates exhibirían cambios conformacionales (tramos de enlace, ángulo curvas, movimientos de torsión, etc.) debido a movimientos térmicos normales, incluyendo las interacciones con las moléculas de solvente. Intentos de incluir los cambios conformacionales de adsorbatos en las simulaciones FFMD requeriría detallado desarrollo de campos de fuerza para los adsorbates superficie catalíticas, es decir, que comprenden términos que describen tramos de enlace, ángulo de codos y términos torsionales, entre otros. Como una dirección futura de este protocolo, estamos desarrollando dichos campos de fuerza para adsorbatos en superficies sólidas, que vamos a utilizar para determinar el grado en que usando adsorbates rígidos influye en los resultados.
Disclosures
Los autores no divulgar conflictos de interés.
Acknowledgments
Esta investigación fue financiada por la National Science Foundation a través del número de concesión 1438325 CBET. Se agradece el apoyo de beca a CJB a través de NASA formación Grant NX14AN43H. Se realizaron simulaciones en el clúster de supercomputación Palmetto, mantenido por el grupo de tecnología de infraestructura cibernética en la Universidad de Clemson. Agradecemos al Dr. Paul J. Morales Meza para el protocolo de la prueba.
Materials
| Name | Company | Catalog Number | Comments |
| VASP software | Computational Materials Physics, Dept. of Physics, University of Vienna | vasp.5.4.4 | Standard parallel VASP executable in the newest version. |
| LAMMPS software | Sandia National Laboratory | 31Mar17-dp | Double-precision, parallel LAMMPS executable from 31 March 2017. |
| VMD software | Theoretical and Computational Biophysics Group, University of Illinois at Urbana-Champaign | 1.9.3 | Standard VMD executable in the newest version. |
| MCPliQ software | Getman Research Group, Dept. of Chemical and Biomolecular Engineering, Clemson University | Executable and input files for the MCPliQ software availabe from the Getman Research Group GitHub page. | |
| JoVE article scripts | Getman Research Group, Dept. of Chemical and Biomolecular Engineering, Clemson University | Python scripts for this JoVE manuscript available from the Getman Research Group GitHub page. | |
| H2O PDB file | Getman Research Group, Dept. of Chemical and Biomolecular Engineering, Clemson University or RCSB Protein Data Bank | PDB file for a water molecule, available from the Getman Research Group GitHub page or at http://www.rcsb.org/ligand/HOH. |
References
- Liu, J. L., Cao, X. M., Hu, P. Density functional theory study on the activation of molecular oxygen on a stepped gold surface in an aqueous environment: A new approach for simulating reactions in solution. Physical Chemistry Chemical Physics. 16 (9), 4176-4185 (2014).
- Okamoto, Y., Sugino, O., Mochizuki, Y., Ikeshoji, T., Morikawa, Y. Comparative study of dehydrogenation at Pt(111)/water and Pt(111)/vacuum of methanol interfaces. Chemical Physics Letters. 377 (1-2), 236-242 (2003).
- Santana, J. A., Cabrera, C. R., Ishikawa, Y. A density-functional theory study of electrochemical adsorption of sulfuric acid anions on Pt(111). Physical Chemistry Chemical Physics. 12 (32), 9526-9534 (2010).
- Artrith, N., Kolpak, A. M. Understanding the composition and activity of electrocatalytic nanoalloys in aqueous solvents: A combination of DFT and accurate neural network potentials. Nano Letters. 14 (5), 2670-2676 (2014).
- Jinnouchi, R., Kodama, K., Morimoto, Y. DFT calculations on H, OH and O adsorbate formations on Pt(111) and Pt(332) electrodes. Journal of Electroanalytical Chemistry. 716, 31-44 (2014).
- Yoon, Y., Rousseau, R., Weber, R. S., Mei, D. H., Lercher, J. A. First-principles study of phenol hydrogenation on Pt and Ni catalysts in aqueous phase. Journal of the American Chemical Society. 136 (29), 10287-10298 (2014).
- Desai, S. K., Pallassana, V., Neurock, M. A periodic density functional theory analysis of the effect of water molecules on deprotonation of acetic acid over Pd(III). Journal of Physical Chemistry B. 105 (38), 9171-9182 (2001).
- Huang, Z. Q., Long, B., Chang, C. R. A theoretical study on the catalytic role of water in methanol steam reforming on PdZn(111). Catalysis Science & Technology. 5 (5), 2935-2944 (2015).
- Chang, C. R., Huang, Z. Q., Li, J. Hydrogenation of molecular oxygen to hydroperoxyl: An alternative pathway for O2 activation on nanogold catalysts. Nano Research. 8 (11), 3737-3748 (2015).
- Faheem, M., Heyden, A. hybrid quantum mechanics/molecular mechanics solvation scheme for computing free energies of reactions at metal-water interfaces. Journal of Chemical Theory and Computation. 10 (8), 3354-3368 (2014).
- Behtash, S., et al. Solvation effects in the hydrodeoxygenation of propanoic acid over a model Pd(211) catalyst. Journal of Physical Chemistry C. 120 (5), 2724-2736 (2016).
- Behtash, S., Lu, J. M., Walker, E., Mamun, O., Heyden, A. Solvent effects in the liquid phase hydrodeoxygenation of methyl propionate over a Pd(111) catalyst model. Journal of Catalysis. 333, 171-183 (2016).
- Norskov, J. K., et al. Origin of the overpotential for oxygen reduction at a fuel-cell cathode. Journal of Physical Chemistry B. 108 (46), 17886-17892 (2004).
- Skachkov, D., Rao, C. V., Ishikawa, Y. Combined first-principles molecular dynamics/density functional theory study of ammonia electrooxidation on Pt(100) electrode. Journal of Physical Chemistry C. 117 (48), 25451-25466 (2013).
- Hibbitts, D. D., Loveless, B. T., Neurock, M., Iglesia, E. Mechanistic role of water on the rate and selectivity of Fischer-Tropsch synthesis on ruthenium catalysts. Angewandte Chemie International Edition. 52 (47), 12273-12278 (2013).
- Abdelrahman, O. A., Heyden, A., Bond, J. Q. Analysis of kinetics and reaction pathways in the aqueous-phase hydrogenation of levulinic acid to form γ-valerolactone over Ru/C. ACS Catalysis. 4 (4), 1171-1181 (2014).
- Wang, H. F., Liu, Z. P. Formic acid xxidation at Pt/H2O interface from periodic DFT calculations integratd with a continuum solvation model. Journal of Physical Chemistry C. 113, 17502-17508 (2009).
- Behtash, S., Lu, J., Faheem, M., Heyden, A. Solvent effects on the hydrodeoxygenation of propanoic acid over Pd(111) model surfaces. Green Chemistry. 16, 605-616 (2014).
- Montemore, M. M., Andreussi, O., Medlin, J. W. Hydrocarbon adsorption in an aqueous environment: A computational study of alkyls on Cu(111). The Journal of Chemical Physics. 145, 074702 (2016).
- Hartnig, C., Grimminger, J., Spohr, E. Adsorption of formic acid on Pt(111) in the presence of water. Journal of Electroanalytical Chemistry. 607, 133-139 (2007).
- Hartnig, C., Grimminger, J., Spohr, E. The role of water in the initial steps of methanol xxidation on Pt (211). Electrochimica Acta. 52 (6), 2236-2243 (2007).
- Hartnig, C., Spohr, E. The role of water in the initial steps of methanol xxidation on Pt (111). Chemical Physics. 319, 185-191 (2005).
- Michel, C., et al. Role of water in metal catalyst performance for ketone hydrogenation: A joint experimental and theoretical study on levulinic acid conversion into gamma-valerolactone. Chemical Communications. 50 (83), 12450-12453 (2014).
- Zope, B. N., Hibbitts, D. D., Neurock, M., Davis, R. J. Reactivity of the gold/water interface during selective oxidation catalysis. Science. 330 (6000), 74-78 (2010).
- Pavlova, A., Meijer, E. J. Understanding the role of water in aqueous ruthenium-catalyzed transfer hydrogenation of ketones. ChemPhysChem. 13 (15), 3492-3496 (2012).
- Saavedra, J., Doan, H. A., Pursell, C. J., Grabow, L. C., Chandler, B. D. The critical role of water at the gold-titania interface in catalytic CO oxidation. Science. 345 (6204), 1599-1602 (2014).
- Desai, S., Neurock, M. A first principles analysis of CO oxidation over Pt and Pt66.7%Ru33.3%(111) surfaces. Electrochimica Acta. 48 (25-26), 3759-3773 (2003).
- Gohda, Y., Schnur, S., Gross, A. Influence of water on elementary reaction steps in electrocatalysis. Faraday Discussions. 140, 233-244 (2008).
- Nie, X. W., Luo, W. J., Janik, M. J., Asthagiri, A. Reaction mechanisms of CO2 electrochemical reduction on Cu(111) determined with density functional theory. Journal of Catalysis. 312, 108-122 (2014).
- Michel, C., Auneau, F., Delbecq, F., Sautet, P. C-H versus O-H bond dissociation for alcohols on a Rh(111) surface: A strong assistance from hydrogen bonded neighbors. ACS Catalysis. 1 (10), 1430-1440 (2011).
- Neurock, M., Wasileski, S. A., Mei, D. From first principles to catalytic performance: tracking molecular transformations. Chemical Engineering Science. 59 (22-23), 4703-4714 (2004).
- Camellone, M. F., Marx, D. On the impact of solvation on a Au/TiO2 nanocatalyst in contact with water. Journal of Physical Chemistry Letters. 4 (3), 514-518 (2013).
- Santana, J. A., Mateo, J. J., Ishikawa, Y. Electrochemical hydrogen oxidation on Pt(110): A combined direct molecular dynamics/density functional theory study. Journal of Physical Chemistry C. 114 (11), 4995-5002 (2010).
- Santana, J. A., Saavedra-Arias, J. J., Ishikawa, Y. Electrochemical hydrogen xxidation on Pt(100): A combined direct molecular dynamics/density functional theory study. Electrocatalysis-US. 6 (6), 534-543 (2015).
- Bodenschatz, C. J., Sarupria, S., Getman, R. B. Molecular-level details about liquid H2O interactions with CO and sugar alcohol adsorbates on Pt(111) calculated using density functional theory and molecular dynamics. Journal of Physical Chemistry C. 119 (24), 13642-13651 (2015).
- Plimpton, S. Fast parallel algorithms for short-range molecular dynamics. Journal of Computational Physics. 117, 1-19 (1995).
- MacKerell, A. D. Jr, et al. All-atom empirical potential for molecular modeling and dynamics studies of proteins. Journal of Physical Chemistry B. 102, 3586-3616 (1998).
- Rappe, A. K., Casewit, C. J., Colwell, K. S., Goddard, W. A. III, Skiff, W. M. UFF, A full periodic table force field for molecular mechanics and molecular dynamics simulations. Journal of the American Chemical Society. 114, 10024-10035 (1992).
- Kahn, K., Bruice, T. C. Parameterization of OPLS-AA force field for the conformational analysis of macrocyclic polyketides. Journal of Computational Chemistry. 23, 977-996 (2002).
- Kresse, G., Furthmuller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Computational Materials Science. 6, 15-50 (1996).
- Kresse, G., Furthmuller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Physical Review B. 54, 11169-11186 (1996).
- Kresse, G., Hafner, J. Ab initio molecular dynamics for liquid metals. Physical Review B. 47, 558-561 (1993).
- Humphrey, W., Dalke, A., Schulten, K. VMD - visual molecular dynamics. Journal of Molecular Graphics and Modelling. 14, 33-38 (1996).
- Xie, T., Sarupria, S., Getman, R. B. A DFT and MD study of aqueous phase dehydrogenation of glycerol on Pt(111): Comparing chemical accuracy versus computational expense in different methods for calculating aqueous phase system energies. Molecular Simulation. 43, 370-378 (2017).
- Li, Y., Zhang, X., Srinath, A., Getman, R. B., Ngo, L. B. Combining HPC and big data infrastructures in large-scale post-processing of simulation data: A case3 study in PEARC ’18. Proceedings of the Practice and Experience on Advanced Research Computing. , Article No. 41 (2018).
- Bodenschatz, C. J., Sarupria, S., Getman, R. B. Correction to "Molecular-level details about liquid H2O interactions with CO and sugar alcohol adsorbates on Pt(111) calculated using density functional theory and molecular dynamics. Journal of Physical Chemistry C. 120, 801 (2016).
