Name: multiscale sampling of a heterogeneous water/metal catalyst interface using density functional theory and force-field molecular dynamics
Uploaded: 2019-04-12T14:45:00
Duration: 10 min 52 s
Description: Watch this Scientific Journal Video about Multiscale Sampling of a Heterogeneous Water/Metal Catalyst Interface using Density Functional Theory and Force-Field Molecular Dynamics at JoVE.com

Esta pesquisa foi financiada pela National Science Foundation através do número do prêmio CBET-1438325. Suporte de bolsa para CJB através da NASA formação Grant NX14AN43H é reconhecido com gratidão. Simulações foram realizadas no Cluster de supercomputador de Palmetto, que é mantido pelo grupo de tecnologia Cyberinfrastructure na Universidade de Clemson. Agradecemos o Dr. Paul J. Meza-Morales para o protocolo de teste.

1. gerar a estrutura de adsorbate
<ol><li>Crie um arquivo de VASP POSCAR compreendendo uma supercélula com condições de contorno periódicas, como você faria se você estava executando simulações de adsorbates nas superfícies metálicas sob vácuo. A supercélula deve incluir uma suposição inicial da estrutura adsorbate e a superfície de metal bem como espaço vazio acima da adsorbate para a adição de moléculas de H2O. Detalhes sobre como fazer isso são fornecidos no anterior trabalho35,44. Nota: É importante que a altura do espaço vácuo seja pelo menos 12 Å acima do topo da adsorbate.</li>
	<li>Relaxe a estrutura e minimizar sua energia utilizando o código da VASP. Detalhes sobre como fazer isso são fornecidos no anterior trabalho35,44. Isto produzirá um arquivo chamado CONTCAR, que será usado na próxima seção.</li>
</ol>2. Adicionar explícitas moléculas de H2O
<ol><li>Adicionar N moléculas de água para o espaço vazio no CONTCAR criado na etapa 1.2 utilizando o código de MCPliQ, onde N = ρV, ρ é a densidade da água, e V é o volume do espaço vácuo acima a adsorbate. Nota: ρ deve ser tomado como a densidade da água, conforme determinado pelo modelo TIP3P/CHARMM água à temperatura de simulação. V vai ser refinado na próxima etapa.<ol>
<li>Especifique as seguintes informações no arquivo de master_input.txt de MCPliQ: O número de moléculas de H2O para adicionar (N) alterando-se o primeiro argumento na linha 28, o caminho para o arquivo water.txt, alterando-se o segundo argumento na linha 28 e o mínimo e altura máxima de supercélula que pode ser ocupada por moléculas de água, alterando a coordenada z mínimo na linha 11 e a coordenada z máximo na linha 12.</li>
		<li>Execute o código MCPliQ, digitando mcpliq através da interface de linha de comando para inserir o arquivo CONTCAR com as moléculas de água. O código imprimirá um ou mais arquivos com a extensão de arquivo. POSCAR. Nota: Se mais de um POSCAR de arquivos são produzidos, eles irão ser nomeados POSCAR_n.POSCAR. Selecione o arquivo onde n é maior.</li>
	</ol>
</li>
	<li>Gerar os arquivos de entrada LAMMPS para uma simulação de NPT e equilibrar o volume da célula usando o FFMD do Ensemble de NPT em LAMMPS.<ol>
<li>Execute o script lmps_bond_angle.py na. Arquivo POSCAR gerado na etapa 2.1.2 digitando lmps_bond_angle.py $filename. POSCAR na interface de linha de comando, onde $filename é o nome da. Arquivo POSCAR gerado na etapa 2.1.2. Este script cria um arquivo chamado $filename. POSCAR.bond_angle_info.txt que lista títulos e ângulos que serão usados no arquivo de dados LAMMPS.</li>
		<li>Abra o VMD e selecione arquivo > Nova molécula na janela principal para abrir a janela do Navegador de arquivos da molécula . Selecione VASP_POSCAR no menu suspenso de Determinar o tipo de arquivo . Clique em procurar e navegue até o $filename. Arquivo POSCAR. Clique em Load para abrir o $filename. Arquivo POSCAR.</li>
		<li>Abra o Console de Tk dentro VMD selecionando extensões > Tk Console a partir da janela principal do VMD.</li>
		<li>Execute o seguinte comando no console do Tk: topo writelammpsdata $WDPATH/data.myadsorbate completo, onde $WDPATH é o diretório no computador onde VMD irá gravar o arquivo de dados LAMMPS e data.myadsorbate é o nome do arquivo de dados LAMMPS.</li>
		<li>Exclua a seção títulos e ângulos na parte inferior do arquivo data.myadsorbate. Em seguida, acrescente as listas de bond e ângulo no arquivo $filename. POSCAR.bond_angle_info.text em data.myadsorbate. Nota: Os índices para o tipo de ligação O-H e o tipo de ângulo H-O-H para as moléculas de água no $filename. Arquivo POSCAR.bond_angle_info.txt são definidas como 1. Assim, tipos de ligação e ângulo para adsorbates devem começar a contar no 2.</li>
		<li>Edite o arquivo data.myadsorbate adicionando os parâmetros Lennard-Jones para a seção de Par Coeffs e os parâmetros de Coulomb para a seção de átomos . Lennard-Jones e Coulomb parâmetros para moléculas de H2O, adsorbate átomos e átomos de superfície de metal precisam ser adicionados. Nota: Os parâmetros Lennard-Jones para Pt átomos, moléculas de água e átomos de adsorbate neste protocolo são obtidos da UFF38, TIP3P/CHARMM37e OPLS-AA39 campos de força, respectivamente. Parâmetros de Coulomb para átomos de moléculas e adsorbates de água são provenientes dos TIP3P/CHARMM37 e OPLS-AA39 campos de força, respectivamente. Parâmetros de Coulomb para Pt átomos estão definidos como 0 no presente protocolo. Alternativamente, cargas parciais calculadas podem ser usadas para os parâmetros de Coulomb para adsorbate átomos e os átomos do Pt.</li>
		<li>Copie o input.equil do arquivo de entrada LAMMPS no diretório $WDPATH. Edite a variável de grupo na linha 34 para indicar os índices do tipo de átomo para o oxigênio da água e átomos de hidrogênio da água e a variável de grupo na linha 35 para indicar os índices do tipo de átomo para os Pt e adsorbate átomos.</li>
		<li>Executar o software LAMMPS digitando mpiexec np - XX lmp_mpi < input.equil na interface de linha de comando, onde XX é o número de núcleos de CPU para usar, e lmp_mpi é o nome do executável LAMMPS. Fazendo isso irá executar uma minimização de energia para refinar a configuração de H2O, seguida de uma simulação de FFMD realizada em número constante de moléculas de H2O (N), volume (V) e temperatura (T) para trazer o água à temperatura de simulação, seguido de uma simulação de FFMD executada em constante N, pressão (P) e a temperatura (T) para determinar a altura fisicamente correta da caixa de simulação. Arquivos de saída que serão usados na seção 3 são chamados data.myadsorbate_npt e log.myadsorbate. Nota: A duração da simulação do TNP deve ser tempo suficiente para incluir um "equilíbrio" executá-los, onde o volume da supercélula trata de estado estacionário e uma "produção", que é usado para o conjunto das médias de amostra (aqui, a altura da supercélula). Durante a execução de equilibração, deve estabilizar o volume da supercélula quando plotado contra o tempo para um valor de estado estacionário. Uma vez que isso ocorre, a simulação de NPT pode ser dito para ser em sua produção. Verificar o equilibrio da simulação do TNP, garantindo que as flutuações na altura da supercélula (lz) são mínimas ou têm convergido para um valor constante. Se ocorrerem grandes flutuações, então re-gera uma configuração de2O H diminuindo o sapateado na linha 92 no arquivo input.equil e repetição etapa 2.2.8 ou começando novamente da etapa 2.1.1.</li>
	</ol>
</li>
</ol>3. extrair a altura apropriada da supercélula
<ol><li>Execute o script get_npt_lz.py no arquivo log.myadsorbate, digitando get_npt_lz.py log.myadsorbate na interface de linha de comando. Este script produz a altura média supercélula da porção de "produção" da simulação NPT no arquivo avg_lz.txt. Nota: O script get_npt_lz.py assume que LAMMPS grava o comprimento da célula z-dimensão (lz) para o arquivo de log.myadsorbate cada 1000 fs (personalizável na linha 20 do script get_npt_lz.py), que é o padrão no arquivo de entrada fornecido input.equil LAMMPS. O script get_npt_lz.py detecta e descarta os 2 primeiros no valor de lz ns (personalizável na linha 19 do script get_npt_lz.py) valores no arquivo log.myadsorbate, como eles compreendem a porção de equilibração da simulação, enquanto os restantes 3 ns compõem o " porção de produção"e são, portanto, utilizadas pelo script get_npt_lz.py para calcular o comprimento de média dimensão-z. Além do arquivo de avg_lz.txt, o script de get_npt_lz.py gera um arquivo chamado npt_data.txt, que fornece valores de lz em função do temporal, bem como um arquivo chamado npt_plot.png, que plota os mesmos dados. O enredo pode ser usado para verificar o equilíbrio da simulação do TNP.</li>
	<li>Reconstrua a supercélula usando a altura média determinada em NPT.<ol>
<li>Copie o arquivo data.myadsorbate_npt em um novo diretório, aqui referido como $WD2PATH e renomeá-lo data.myadsorbate.</li>
		<li>Editar o novo arquivo data.myadsorbate para que a altura de lz é igual à valor médio de saída do script get_npt_lz.py alterando os argumentos zlo e zhi, no arquivo data.myadsorbate zlo 0.0 e zhi é o valor de lz do arquivo avg_lz.txt produzido em s 3.1 de tep.</li>
	</ol>
</li>
</ol>4. gerar configurações de moléculas de H2O
<ol><li>Copie o input.prod de arquivo de entrada LAMMPS $WD2PATH. Edite a variável de grupo na linha 32 para indicar os índices do tipo de átomo para o oxigênio da água e átomos de hidrogênio da água e a variável de grupo na linha 33 para indicar os índices do tipo de átomo para os Pt e adsorbate átomos.</li>
	<li>Executar o software LAMMPS digitando mpiexec np - XX lmp_mpi < input.prod na interface de linha de comando, onde XX é o número de núcleos de CPU para usar, e lmp_mpi é o nome do executável LAMMPS. Fazendo isso será executado uma simulação de NVT constante sobre as moléculas de H2O. A chave de saída desta simulação é o arquivo dump.myadsorbate.lammpstrj. Nota: A duração da simulação NVT deve ser tempo suficiente para compõem um equilibrio executá-los, onde a energia do sistema vem para firmar o estado, e uma produção, do qual o conjunto das médias (aqui, as posições espaciais das moléculas de água) são amostrados. Durante a execução de equilibração, deve estabilizar a energia do sistema quando plotado contra o tempo para um valor de estado estacionário. Uma vez que isso ocorre, a simulação de NVT pode ser dito para ser em sua produção.</li>
</ol>5. determinar o tempo de vida de ligação de hidrogênio para amostragem de tempo próprio
<ol><li>Edite o script hb_lifetime_dist.py para especificar: a iteração do primeiro quadro do arquivo dump.myadsorbate.lammpstrj, alterando a variável actualStart na linha 22, quantas vezes quadros são gravados no arquivo de trajetória LAMMPS alterando a variável de iteração no linha 23, recomendando-se o primeiro e o último timesteps o script (ou seja, a parte de produção da trajetória), alterando as variáveis N_first e N_last nas linhas 24 e 25, se frames consecutivos são considerados ou quadros são ignorados, alterando o nevery variável na linha 26 e o número de linhas por seção de quadros do arquivo trajetória alterando a variável frameLine na linha 27. Além disso, edite linhas 31 a 35 para especificar quais tipos de átomo dentro do arquivo data.myadsorbate pertencem a adsorbate e quais tipos de átomo pertencem as moléculas de H2O. Nota: O script hb_lifetime_dist.py analisa as configurações de2O H o prazo de produção e determina se qualquer moléculas de O2H são hidrogênio ligado ao adsorbate. Ele então conta o tempo de simulação que cada hidrogênio ligação permanece intacta e relata esta informação como uma distribuição de vidas de ligação de hidrogênio em unidades do ps. A versão específica do script que é fornecido com este protocolo assume que LAMMPS grava a configuração das moléculas de H2O arquivo de dump.myadsorbate.lammpstrj todos os fs 1000, que é o padrão no arquivo de entrada fornecido input.prod LAMMPS. Ele detecta e descarta as 2 primeiras vale ns de configurações no o dump.myadsorbate.lammpstrj de arquivo, como eles compreendem a porção de equilibração da simulação e utiliza o restante 3 ns para calcular a vida útil da ligação de hidrogênio.</li>
	<li>Execute o script hb_lifetime_dist.py no arquivo dump.myadsorbate.lammpstrj, digitando hb_lifetime_dist.py na interface de linha de comando. Fazendo isto irá produzir um arquivo chamado distribution_HB_lifetime.dat.</li>
	<li>Plotar os dados no arquivo distribution_HB_lifetime.dat para visualizar a distribuição das vidas de ligação de hidrogênio que ocorreram durante a simulação de NVT.</li>
	<li>Determine o incremento de tempo para usar o intervalo de amostragem baseado em vidas a ligação de hidrogênio calculado. A melhor escolha é o tempo de vida máximo de ligação de hidrogênio; Alternativamente, pode ser usado um valor que capturasse o intervalo de confiança de 95%.</li>
</ol>6. configurações de amostra de líquido moléculas de H2O
<ol><li>Determine o número de configurações a partir da execução de produção da trajetória do NVT FFMD para cálculos ainda mais. O número de configurações deve ser selecionado para que o tempo mínimo entre configurações é igual ou maior que o intervalo de amostragem de tempo identificado na seção 5.<ol>
<li>Edite o valor padrão para a variável num_frames na linha 21 do script lammps_frames.py para especificar o número de configurações para extrair.</li>
		<li>Execute o script lammps_frames.py dump.myadsorbate.lammpstrj o arquivo digitando lammps_frames.py na interface de linha de comando. Fazendo isto produzirá uma lista de simulação vezes correspondente às configurações que devem ser extraídas do arquivo dump.myadsorbate.lammpstrj. Essas configurações podem ser usadas como estruturas a partir de simulações AIMD ou QM. Nota: 1) o script lammps_frames.py detecta automaticamente o LAMMPS log e arquivos de despejo, bem como a parte de produção da trajetória dentro do arquivo de despejo e divide o número de configurações dentro do arquivo de despejo em 10 grupos. Como alternativa, o usuário pode especificar o arquivo de log, o arquivo de despejo e o número de configurações da interface de linha de comando usando o -l, -d e opções - n, respectivamente. Para fazer isso, o usuário deve digitar lammps_frames.py - n XX -l $logfilename -d $dumpfilename na interface de linha de comando, onde XX é o número desejado de configurações, $logfilename é o nome do arquivo de log LAMMPS e $dumpfilename é o nome do Arquivo LAMMPS trajetória (despejo). Os tempos de simulação que são saída referem-se às vezes medianos em cada grupo. 2) se as configurações serão calculadas na VASP com a bandeira LDIPOLE ativada, uma pequena camada de espaço vazio deve ser adicionada ao topo da supercélula acima da camada de água. Isto irá facilitar a convergência da estrutura eletrônica para o cálculo da VASP. Adicionar um Å 3 adicionais de espaço vazio acima das moléculas de O2H tem sido bem sucedido nas simulações discutidas abaixo.</li>
	</ol>
</li>
</ol>

Os autores divulgar sem conflitos de interesse.

O método apresentado foi selecionado por sua facilidade de implementação, mas várias personalizações poderiam ser feitas. Por um lado, os campos de força usados nas simulações FFMD pode ser modificados. Alterar os parâmetros de campo de força e/ou potenciais pode ser feito editando os arquivos de dados e entrada LAMMPS. Da mesma forma, poderiam ser empregados solventes que não sejam de H2O. Para fazer essa modificação, a molécula do solvente desejada precisa ser inserido a partir de passo 2.1.1, e os arquivos de entrada LAMMPS precisaria ser editado para incorporar as potencialidades apropriadas e parâmetros. Inserindo a nova molécula de solvente também exigiria fornecendo as coordenadas internas da molécula do solvente em um arquivo txt. análogo ao arquivo water.txt.
Outra modificação que poderia ser feita é modificar a área da superfície laje. Os resultados discutidos neste manuscrito empregadas lajes superfície 3 Pt x 3 Pt ou Pt 4 x 4 Pt, que têm áreas de superfície menos de 120 Å2. À medida que aumenta a área de superfície da laje, o gasto computacional também aumenta. Gasto computacional tem o maior impacto na seção 5 do presente protocolo. Se as etapas de processamento de dados na secção 5 tornam-se computacionalmente proibitivas, grande volume de dados post estratégias de processamento tais como aqueles discutidos em Li et al 201845 podem ser empregadas.
Possíveis fontes de incerteza para este procedimento incluem o campo de força empregado, o método de amostragem e a frequência de amostragem. A estrutura da água é determinada pelo campo de força que é usado, significa que a escolha do campo de força pode influenciar as configurações específicas de moléculas de H2O. Nosso grupo avaliou como as escolhas de campo de força para H2O moléculas e átomos Pt influenciam as energias de interação, calculado em FFMD e achei que a escolha do campo de força contribui menos de 0,1 eV a esta energia de interação. Outra fonte de incerteza é o método de amostragem, que influencia as configurações específicas que são usadas para calcular uma quantidade de interesse. Nosso grupo tem comparado o desempenho do método "amostragem de tempo" apresentado no presente protocolo, com um método de "amostragem de energia", que é tendencioso para configurações de energia menores de moléculas de H2O, na interação energias calculadas em DFT e encontraram os dois Esses métodos de amostragem estatisticamente igual dar valores35,46. A frequência de amostragem também pode influenciar os resultados. Temos encontrado que a frequência de amostragem contribui menos e avaliados como aumentar o número de configurações de 10 a 30.000 influencia as energias de interação média calculadas em FFMD para 40 diferentes C3HxO3 adsorbates do que 0,1 eV para a energia de interação média44.
A principal limitação para este método é que os adsorbates são aproximados pelas estruturas sob vácuo durante as simulações de FFMD. Na realidade, os adsorbates que exibem mudanças conformacionais (trechos de ligação, curvas de ângulo, movimentos de torção, etc.), devido a movimentos térmicos normais, incluindo interações com moléculas de solventes. Tentativas de incluir mudanças conformacionais de adsorbates para as simulações de FFMD exigiria desenvolvimento detalhado de campos de força para adsorbates superfície catalíticas, ou seja, que compreendem os termos que descrevem os trechos de ligação, curvas de ângulo e termos de torção, entre outros. Como uma direção futura do presente protocolo, estamos a desenvolver tais campos de força para adsorbates em superfícies sólidas, que será usada para determinar a extensão a que usar adsorbates rígidas influencia os resultados.

Modelagem moleculares fenômenos envolvidos na catálise heterogênea sob condições de líquido é necessária para a função catalítica de compreensão; no entanto, esta continua a ser desafiador porque ele requer um delicado equilíbrio entre precisão química e custo computacional. Em geral, como catálise envolve a quebra e formação de ligações químicas, mecânica quântica deve ser usado para pelo menos algum grau; no entanto, longas simulações estão desafiando na mecânica quântica, como eles exigem recursos significativos do computador. Desde moléculas na fase líquida estão em constante movimento térmico, simulações também devem incluir a amostragem configuracional, ou seja, eles devem incorporar vários arranjos espaciais das moléculas do líquido, como cada arranjo espacial diferente (ou seja, cada um configuração) tem uma energia diferente. Isto significa que várias configurações de moléculas do líquido devem ser simuladas para cada espécie catalítica de interesse. Estas necessidades – para usar a mecânica quântica e realizar vários cálculos por espécie catalítica – podem processar modelagem em catálise heterogênea em fase líquida intratáveis computacionalmente. O objetivo do método aqui descrito é permitir computacionalmente tractable simulações de fenômenos em catálise heterogênea em fase líquida.
Estamos particularmente interessados em reações catalisadas forma heterogénea que sejam efectuadas sob a água líquida. As moléculas de água têm uma influência significativa sobre fenômenos catalíticos, tais como interação com espécies catalítico (por exemplo, via forças de dispersão e ligações de hidrogênio)1,2,3,4,5 ,6,7,8,9,10,11,12,13,14 , 15 , 16 , 17 , 18 , 19 , 20 , 21 , 22 , 23, participando de reações catalíticas1,7,8,9,15,21,22,24 ,25,26,27e influenciando os caminhos de reação e/ou taxas catalítica1,11,12,15, 18,23,25,,27,28,29,30,31. Modelagem desses fenômenos foi executada usando QM e/ou ab initio a dinâmica molecular (AIMD)1,2,6,7,14,22 ,25,27,28,32,33,34, força dinâmica molecular de campo (FFMD)35 e a mecânica quântica/molecular mecânica (QM/MM)10. Em AIMD e FFMD, os átomos no sistema são movidos em conformidade com as equações de Newton de movimento de acordo com as forças que atuam sobre eles. Em AIMD, a energia do sistema e as forças são calculadas com a mecânica quântica, Considerando que em FFMD, a energia do sistema e as forças são calculadas usando força campos, que são expressões algébricas que são parametrizadas com base em experimental ou dados de QM. Em QM/MM, a parte do sistema onde ocorre o vínculo quebrando e formando é calculada com QM, e o restante do sistema é calculado com MM, que emprega campos de força. Porque eles empregam diretamente QM, AIMD e QM/MM são mais adequados para capturar a quebra de vínculo e formar que ocorre na catálise heterogênea da fase aquosa; no entanto, FFMD é significativamente mais computacionalmente tractable e, portanto, mais adequado para gerar as configurações das moléculas de líquido H2O. O método apresentado neste protocolo equilibra precisão química e despesa computacional empregando uma combinação de QM e FFMD.
Especificamente, esse método usa simulações de FFMD para gerar configurações de líquido H2O e QM para calcular as energias do sistema. FFMD é realizada usando LAMMPS. 36 os campos de força usados em FFMD neste trabalho empregar Lennard-Jones + potenciais de Coulomb (LJ + C), onde os parâmetros LJ foram retirados do modelo TIP3P/CHARMM37 para H2O, o campo de força universal38 (UFF) para Pt e o Campo de força de OPLS-AA39 espécies catalítico e os parâmetros de Coulomb foram retirados do modelo de37 TIP3P/CHARMM para H2O e o campo de força OPLS-AA39 espécies catalítico. Os parâmetros de Coulomb para Pt átomos tem sido definidos como 0. Cálculos de QM são executados usando a VASP código40,41,,42, que é um código de teoria funcional (DFT) densidade. Inserções de molécula de água são executadas com um código desenvolvido in-house chamado Monte Carlo plug-in para métodos de Quantum (MCPliQ). Arquivo conversões de VASP LAMMPS neste protocolo são executadas com o software de dinâmica Molecular Visual (VMD)43.
O protocolo destina-se para gerar configurações de moléculas de água líquida em torno de espécie catalítica em superfícies de metal de transição plana na cobertura de baixa. Cobertura é denotado por θ e definida como o número de adsorbates por superfície metal átomo (ou seja, o número de adsorbates de superfície normalizado pelo número de átomos de metal na camada superior da laje metal no modelo de catalisador). Neste manuscrito, baixa cobertura é definida como θ ≤ 1/9 monocamada (ML), onde 1 ML significa uma espécie catalítica por superfície átomo de metal. Os modelos de catalisador devem ser colocados em caixas de simulação periódica. As caixas de simulação não tem que ser cubos. Este manuscrito demonstra o uso do protocolo para gerar configurações de líquido H2O que pode ser usado para calcular as quantidades de interesse em catálise heterogênea de fase aquosa.
Este protocolo requer que o usuário tem acesso às versões instalados e funcionando do software VASP, MCPliQ, LAMMPS e VMD. Obter mais informações sobre a VASP (https://www.vasp.at/), LAMMPS (https://Lammps.sandia.gov/) e VMD (https://www.ks.uiuc.edu/Research/vmd/) estão disponíveis em seus respectivos sites. O software MCPliQ está documentado em https://github.com/getman-research-group/JoVE_article, juntamente com todos os arquivos de entrada e os scripts de Python mencionados no presente protocolo. Este protocolo assume que os arquivos executáveis e scripts mencionadas dentro serão executados em um computador de alto desempenho de pesquisa e são instalados em um diretório que está na variável $PATH do usuário. Se um arquivo executável ou script é colocado em um local que não seja no do usuário $PATH, em seguida, o caminho para o executável deve ser incluído para executá-lo. Executáveis e scripts são executados em etapas, 2.1.2, 2.2.1, 2.2.8, 3.1, 4.2, 5.2 e 6.1.2. Por exemplo executar o código MCPliQ na etapa 2.1.2 de um diretório que não está no usuário do $PATH, o usuário deve digitar $PATHTOMCPLIQ/mcpliq na interface de linha de comando em vez de mcpliq, onde $PATHTOMCPLIQ é o local onde o mcpliq executável foi armazenado (por exemplo, $PATHTOMCPLIQ pode ser ~ / bin). Antes de iniciar este protocolo, todos os arquivos executáveis e scripts devem dar permissões executáveis (por exemplo, no Linux, isso poderia ser feito digitando chmod + x mcpliq na interface de linha de comando do diretório onde o executável do mcpliq é armazenado). Além disso, quaisquer módulos requeridos por qualquer um do software ou scripts devem ser carregados (essas dependências serão específicas para cada uma das instalações de vários software e o computador onde serão executadas as simulações).

<table>
	<tbody>
		<tr>
			<td>VASP software</td>
			<td>Computational Materials Physics, Dept. of Physics, University of Vienna</td>
			<td>vasp.5.4.4</td>
			<td>Standard parallel VASP executable in the newest version.</td>
		</tr>
		<tr>
			<td>LAMMPS software</td>
			<td>Sandia National Laboratory</td>
			<td>31Mar17-dp</td>
			<td>Double-precision, parallel LAMMPS executable from 31 March 2017.</td>
		</tr>
		<tr>
			<td>VMD software</td>
			<td>Theoretical and Computational Biophysics Group, University of Illinois at Urbana-Champaign</td>
			<td>1.9.3</td>
			<td>Standard VMD executable in the newest version.</td>
		</tr>
		<tr>
			<td>MCPliQ software</td>
			<td>Getman Research Group, Dept. of Chemical and Biomolecular Engineering, Clemson University</td>
			<td></td>
			<td>Executable and input files for the MCPliQ software availabe from the Getman Research Group GitHub page.</td>
		</tr>
		<tr>
			<td>JoVE article scripts</td>
			<td>Getman Research Group, Dept. of Chemical and Biomolecular Engineering, Clemson University</td>
			<td></td>
			<td>Python scripts for this JoVE manuscript available from the Getman Research Group GitHub page.</td>
		</tr>
		<tr>
			<td>H2O PDB file</td>
			<td>Getman Research Group, Dept. of Chemical and Biomolecular Engineering, Clemson University or RCSB Protein Data Bank</td>
			<td></td>
			<td>PDB file for a water molecule, available from the Getman Research Group GitHub page or at http://www.rcsb.org/ligand/HOH.</td>
		</tr>
	</tbody>
</table>

multiscale sampling of a heterogeneous water/metal catalyst interface using density functional theory and force-field molecular dynamics

Um número significativo de forma heterogénea catalisada por processos químicos ocorre sob condições de líquido, mas simulando a função de catalisador sob tais condições é desafiador quando é necessário incluir as moléculas de solventes. A quebra de vínculo e formando processos modelados nestes sistemas exigem o uso de métodos de química quântica. Desde que as moléculas na fase líquida estão em constante movimento térmico, simulações também devem incluir a amostragem configuracional. Isto significa que várias configurações de moléculas do líquido devem ser simuladas para cada espécie catalítica de interesse. O objetivo do protocolo aqui apresentado é gerar e trajetórias de configurações de moléculas de água líquida em torno de espécie catalítica em superfícies planas de metal de transição de uma forma que equilibra precisão química com gasto computacional da amostra. Especificamente, simulações de dinâmica molecular (FFMD) de campo de força são usadas para gerar configurações de moléculas de líquido que posteriormente podem ser usadas em métodos baseados em mecânica quântica, tais como a teoria do funcional da densidade ou ab initio molecular dinâmica. Para ilustrar isso, neste manuscrito, o protocolo é usado para intermediários catalíticos que poderiam estar envolvidos na via para a decomposição de glicerol (C3H8O3). As estruturas que são geradas usando FFMD são modeladas em DFT para estimar os enthalpies de solvatação da espécie catalítica e identificar como moléculas de H2O participarem decomposições catalíticas.

Um uso do presente protocolo é calcular as energias de interação entre a água líquida e a espécie catalítica, ou seja, ΔEint35:
∆Eint=Eespécie catalítica + H2O+Esuperfície limpa catalisador-Eespécie catalítica-Ecatalisador limpa superfície + H2O
onde Eespécie catalítica + H2Oé a energia de uma configuração de moléculas de H2O em torno de uma espécie catalítica sobre uma superfície de metal, Esuperfície do catalisador limpa é a energia da superfície do catalisador limpo no vácuo, E Espécie catalítica é a energia da espécie catalítica sobre uma superfície de metal no vácuo, e Elimpar o catalisador superfície + H2O é a energia da configuração de H2O sobre a superfície do catalisador com a espécie catalítica removida. As posições das moléculas de H2O usado para calcular Eespécie catalítica + H2O e Elimpar o catalisador superfície + H2O devem ser idênticos. Todos os valores de E são calculados usando o código da VASP. A quantidade ΔEint inclui todas as interações físicas e químicas entre todas as moléculas na estrutura de água líquida e a espécie catalítica e dá uma estimativa razoável da entalpia de solvatação dos catalítico espécie, que é necessário para calcular sua energia livre de solvatação e total de energia livre. A tabela 1 fornece valores para ΔEint calculado para a espécie em um catalisador Pt(111) de superfície com fórmulas químicas iguais a CxHyOz em unidades de eV (1 eV = 96.485 kJ/mol). Os valores foram calculados em coberturas ≤ 1/9 ML.35,46 os valores relatados são as médias assumidas 10 configurações de líquido H2O, e as incertezas são relatadas como desvios-padrão. Todos os valores são negativos, indicando interações favoráveis com água.
Outra aplicação do presente protocolo é gerar estruturas partidas para AIMD. Filme 1 é um filme de uma trajetória AIMD que foi iniciado a partir de uma configuração gerada pelo presente protocolo. No início do filme, um adsorbate COH é mostrado em uma superfície de Pt(111) sob uma estrutura de líquido H2O. Molécula de um H2O é enfatizada, que formou uma ligação de hidrogênio com COH. Ao longo do filme, esta molécula de2O H abstrai o próton do adsorbate COH e deposita um segundo átomo de hidrogênio na superfície Pt(111). A molécula de H2O, assim, ajuda a catalisar a reação COH * + * → CO * + H *, onde o * s indicar sites catalíticos. Esta simulação destaca a força principal e o objetivo principal do método de amostragem múltipla escala aqui descrito. Inúmeras configurações de moléculas de H2O são geradas com FFMD, devido a sua força na rastreabilidade computacional. No entanto, uma limitação do FFMD é que ele não pode capturar quebra de vínculo e formando a menos que um campo de força reativo é implementado. AIMD usa mecânica quântica para calcular as energias e assim pode capturar o vínculo quebrando e formando. No entanto, AIMD demasiado computacionalmente é exigente gerar todas as configurações de moléculas de H2O necessárias para garantir a amostragem suficiente foi alcançada. Assim, o presente protocolo combina os dois métodos.
As estruturas do líquido H2O as moléculas geradas por esse procedimento dependem de configurações de entrada. Estas configuração incorretamente pode ter influências involuntárias nas estruturas de água. Por exemplo, quando as distâncias intermoleculares se tornam muito pequenas ou quando outros parâmetros nos arquivos de entrada de dinâmica molecular são definidos incorretamente ou nos valores não físicos, a estrutura da água pode tornar-se irracional. Nestas circunstâncias, a estrutura da água irá "explodir" unintendedly durante a trajetória FFMD. A Figura 1 mostra um exemplo disso. O instantâneo do lado esquerdo é a estrutura inicial para uma corrida FFMD, e o instantâneo do lado direito é um instantâneo tomado no prazo de 1 ps de iniciar a simulação. Como pode ser visto, as moléculas de O2H moveram longe da superfície. Isto é causado pelo indevido configurações feitas em arquivos de simulação de entrada e não é uma estrutura que é provável de ocorrer na realidade.
<img alt="Figure 1" class="xfigimg" src="/files/ftp_upload/59284/59284fig1v2.jpg"> Figura 1: Exemplo de um resultado negativo. A simulação de dinâmica molecular de campo de força "explode" devido a uma configuração física ou valor. Imagem da mão esquerda: A geometria inicial da superfície do Pt(111), adsorbate e estrutura de água líquida. Imagem direita: A geometria da superfície de Pt(111), adsorbate, água líquida e estrutura inferior a 1 ps mais tarde. Na imagem da direita, as moléculas de H2O tem separado da superfície devido às forças de unphysically grandes. <a href="https://www.jove.com/files/ftp_upload/59284/59284fig1v2large.jpg" target="_blank">Clique aqui para ver uma versão maior desta figura.</a>
<img alt="Movie 1" src="/files/ftp_upload/59284/59284mov1.jpg"> Filme 1: Ab initio a dinâmica molecular (AIMD) simulação iniciada a partir de uma configuração gerada no amostragem múltipla escala. Uma molécula de2O H que é originalmente hidrogênio ligado a um adsorbate COH sobre uma superfície de Pt(111) abstrai o próton do COH e deposita um segundo hidrogênio na superfície Pt(111). Este vínculo quebrando e formando o evento pode ser capturado por AIMD, mas não com campo de força dinâmica molecular (FFMD) a menos que um campo de força reactivo é usado. A configuração inicial de moléculas de H2O usado nesta simulação AIMD foi gerada usando o FFMD, conforme descrito neste manuscrito. <a href="https://www.jove.com/files/ftp_upload/59284/movie1.mp4" target="_blank">Por favor clique aqui para ver este vídeo. (Botão direito do mouse para fazer o download.)</a>
<table border="1" fo:keep-together.within-page="1" fo:keep-with-next.within-page="always"><tbody>
<tr>
<td>Espécie catalítica</td>
			<td>∆ eint (eV)</td>
		</tr>
<tr>
<td>COH</td>
			<td>-0.70 ± 0.07</td>
		</tr>
<tr>
<td>CO</td>
			<td>-0.03 ± 0,03</td>
		</tr>
<tr>
<td>CH2OH</td>
			<td>-0.64 ± 0,12</td>
		</tr>
<tr>
<td>CHO-CHOH-CH2OH</td>
			<td>-0.93 ± 0.22</td>
		</tr>
<tr>
<td>COH-COH-CH2OH</td>
			<td>-0.87 ± 0.23</td>
		</tr>
<tr>
<td>COH-CHOH-COH</td>
			<td>-1.72 ± 0,26</td>
		</tr>
<tr>
<td>CHOH-COH-CO</td>
			<td>-1.57 ± 0,25</td>
		</tr>
<tr>
<td>CHO-CO-CO</td>
			<td>-0.31 ± 0,19</td>
		</tr>
</tbody></table>Tabela 1: Água-catalítica resultados de energia de interação espécie. Energias de interação em eV calculado para oito CxHyOz adsorbates na Pt(111). Valores relatados são as médias assumiu várias configurações de líquido H2O. As incertezas são os desvios-padrão das médias. 1 eV = 96.485 kJ/mol.

Watch this Scientific Journal Video about Multiscale Sampling of a Heterogeneous Water/Metal Catalyst Interface using Density Functional Theory and Force-Field Molecular Dynamics at JoVE.com

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

O objetivo do protocolo aqui apresentado é gerar e trajetórias de configurações de moléculas de água líquida em torno de espécie catalítica sobre uma superfície plana de metal de transição da amostra. As configurações incluídas na amostra podem ser usadas como estruturas a partir de métodos baseados na mecânica quântica.

Amostragem múltipla escala de uma Interface de catalisador heterogêneo água/Metal usando a teoria do funcional da densidade e do campo de força dinâmica Molecular

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

Os modelos de protocolo, fenômenos térmicos e quânticos em fase líquida de catalise heterogênea. É o primeiro, ao nosso conhecimento, a incorporar a mecânica quântica com amostragem completa de um ambiente líquido explícito. As configurações de moléculas líquidas geradas a partir deste protocolo representam as esperadas em condições reais de reação e podem ser usadas para explorar fenômenos de nível molecular que dependem dos arranjos espaciais das moléculas.As configurações de moléculas líquidas geradas por este protocolo fornecem insights sobre os papéis que o solvente desempenha na fase líquida catálise heterogênea. Se você está testing este protocolo pela primeira vez, recomendo que você primeiro garanta que você tenha acesso e possa executar o software VASP, MCPliQ, VMD e LAMMPS. Demonstrando o procedimento será Tianjun Xie, um estudante de pós-graduação do meu laboratório.Depois de gerar a estrutura adsorbate, gere os arquivos de entrada LAMMPS para uma simulação de NPT e equilibre o volume da célula usando FFMD. Copie o arquivo de entrada LAMMPS no diretório de trabalho. Edite a variável de grupo na linha 34 para indicar os índices do tipo átomo para os átomos de oxigênio da água e hidrogênio da água.Na variável grupo na linha 35, para indicar os índices do tipo átomo para os átomos de platina e adsorbato. Na linha 17 do arquivo de entrada, edite a variável etapa de execução para definir a duração da simulação NPT tempo suficiente para incluir uma execução de equilíbrio e uma execução de produção. Execute o software LAMMPS, digitando o comando na interface de linha de comando, que inclui as informações do número de núcleos de CPU para usar e o nome do LAMMPS executável.A minimização de energia refina a configuração da molécula de água e é seguida por uma simulação de FFMD realizada em número constante de moléculas de água, volume e temperatura para levar a água à temperatura da simulação. Outra simulação FFMD é então executada em número constante de moléculas de água, pressão e temperatura para determinar a altura fisicamente correta da caixa de simulação. Os arquivos de saída serão usados mais tarde.Após a simulação do NPT, plote a altura da supercélula contra o tempo. O ponto em que ele se nivela para um valor de estado estável é o ponto na simulação NPT quando a produção pode começar. Verifique o equilíbrio da simulação do NPT, garantindo que as flutuações na altura da supercélula sejam mínimas ou tenham convergido para um valor constante.Se ocorrerm grandes flutuações, abra a entrada. equil arquivo e diminuir o passo de tempo na linha 92 para regenerar a configuração da molécula de água e executar o software LAMMPS como antes. Para começar, digite a interface de linha de comando para executar o script.Este script produz a altura média da supercélula da parte de execução de produção da simulação NPT para um arquivo TXT. O script lê o comprimento da dimensão z da célula em intervalos de 1.000 femtosegundos, que é o intervalo padrão para impressão de informações em LAMMPS. Se um intervalo de impressão diferente for desejado, ele pode ser alterado pela linha de edição 20 da get_npt_lz.py script e linha 16 de entrada LAMMPS. arquivo equil. O script detecta e descarta os dois primeiros nanossegundos do valor dos valores lz, pois compreendem a parte de equilíbrio da simulação.A duração da execução de equilíbrio pode ser alterada pela linha de edição 19 do arquivo. Os três nanossegundos restantes compõem a porção de produção e, portanto, são usados para calcular o comprimento médio da dimensão Z. Além disso, o script produz outro arquivo TXT que fornece valores de lz em função do passo de tempo, bem como um arquivo PNG, que plota os mesmos dados.O enredo pode ser usado para verificar o equilíbrio da simulação do NPT. Para reconstruir a supercélula usando a altura média determinada no NPT, copie o arquivo de dados gerado anteriormente em um novo diretório de trabalho e renomeie-o como data.myadsorbate. Em seguida, edite o novo arquivo de dados para alterar zlo para 0.0 e zhi para o valor lz a partir da saída de valor médio no arquivo TXT.Copie o arquivo de entrada LAMMPS no novo diretório de trabalho. Edite a variável de grupo na linha 32 para indicar os índices do tipo átomo para os átomos de oxigênio da água e hidrogênio da água e a variável grupo na linha 33 para indicar os índices do tipo átomo para os átomos de platina e adsorbate. Em seguida, na linha 16, edite a variável runStep para que seja longa o suficiente para incluir uma corrida de equilíbrio e uma execução de produção.Digite o comando para executar LAMMPS na interface de linha de comando para executar o software LAMMPS. Isso executará uma simulação NVT constante nas moléculas de água, e o arquivo de saída chave é gerado. A simulação de NVT compreende uma parte de equilíbrio e uma parte de produção.A parte de produção começa quando a energia do sistema plotada contra os níveis de tempo desligados. Para começar, abra o script vitalício da ligação de hidrogênio. Altere a variável inicial real na linha 22 para definir o passo do tempo do primeiro período de tempo.Altere a variável timestep na linha 23 para definir com que frequência os quadros são escritos no arquivo de trajetória LAMMPS. Altere N_first e N-última variáveis na linha 24 e 25 para definir os primeiros e últimos passos do tempo, e altere a variável nevery na linha 26 para definir se quadros consecutivos são considerados ou ignorados. Defina o número de linhas por seção de quadro do arquivo de trajetória alterando a variável da linha de quadro na linha 27.Além disso, edite as linhas 31 a 35, para especificar quais tipos de átomos dentro dos dados. o arquivo myadsorbate pertence ao adsorbate, e quais tipos de átomos pertencem às moléculas de água. O script analisa as configurações de água na execução da produção e determina se alguma molécula de água está ligada a hidrogênio ao adsorbato.Em seguida, conta o tempo de simulação de que cada ligação de hidrogênio permanece intacta e relata essas informações como uma distribuição de vida útil de ligação de hidrogênio em unidades de picosegundos. O LAMMPS grava a configuração de moléculas de água no arquivo a cada 1.000 femtosegundos, que é o padrão no arquivo de entrada LAMMPS fornecido. Ele detecta e descarta os dois primeiros nanossegundos de configurações no arquivo, pois compreendem a parte de equilíbrio da simulação e usa os três nanossegundos restantes para calcular as vidas de ligação de hidrogênio.Para executar o script, digite na interface de linha de comando. Em seguida, gera um arquivo DAT. Plote os dados no arquivo para visualizar a distribuição de vida útil de títulos de hidrogênio que ocorreu durante a simulação NVT.Para determinar o incremento de tempo a ser usado para o intervalo de amostragem de tempo, use um incremento de tempo maior ou igual à vida útil máxima da ligação de hidrogênio. Determine o número de configurações a partir da execução de produção da trajetória NVT FFMD para amostrar de tal forma que o tempo mínimo entre as configurações seja igual ou maior do que o intervalo de amostragem de tempo identificado anteriormente. No script de extração de quadros escrito anteriormente, edite o valor padrão para a variável número de quadros na linha 21 para especificar o número de configurações a serem extraídas.Para executar o script, digite o nome do script na interface de linha de comando. Isso irá produzir uma lista de tempos de simulação correspondentes às configurações que devem ser extraídas do arquivo de simulação NVT. Essas configurações podem ser usadas como estruturas de partida em simulações AIMD ou QM.Neste procedimento, o FFMD foi utilizado para gerar a configuração inicial das moléculas de água. A simulação AIMD mostra que uma molécula de água que é originalmente ligada a hidrogênio a um adsorbate de álcool açúcar em uma superfície de platina-111 abstrai o hidrogênio do adsorbate do álcool e deposita um segundo hidrogênio na superfície de platina-111. As estruturas das moléculas de água líquida dependem das configurações de entrada.Definir essas influências inadequadas pode ter influências não intencionais nas estruturas da água. Nesta figura, o lado esquerdo é a estrutura inicial para uma corrida FFMD. E o lado direito está a um segundo de iniciar a simulação.A simulação FFMD explode devido a configurações de força unphysicamente grandes, fazendo com que as moléculas de água se movam para longe da superfície. As configurações podem ser usadas na mecânica quântica, ou simulação QMM, ou podem ser usadas para analisar estatísticas relacionadas às posições espaciais das moléculas. Essa técnica abre caminho para os pesquisadores explorarem os papéis que os ambientes de reação líquida têm sobre a catálise, gerando configurações reais de moléculas líquidas em interfaces catalíticas.

Research

JoVE Journal

Chemistry

Layer-by-layer Synthesis and Transfer of Freestanding Conjugated Microporous Polymer Nanomembranes

Síntese camada por camada e Transferência de autônomos conjugados microporosas Polymer Nanomembranes

CMP as large surface area materials have attracted growing interest recently, due to their high variability in the incorporation of functional groups in ...

Química

HKUST-1 as a Heterogeneous Catalyst for the Synthesis of Vanillin

HKUST-1 como um catalisador heterogéneo para a Síntese de vanilina

Vanillin (4-hydoxy-3-methoxybenzaldehyde) is the main component of the extract of vanilla bean. The natural vanilla scent is a mixture of approximately ...

Facile Synthesis of Worm-like Micelles by Visible Light Mediated Dispersion Polymerization Using Photoredox Catalyst

Síntese Facile de micelas Worm-like por luz visível Mediated Dispersão Polimerização Usando Photoredox Catalisador

Presented herein is a protocol for the facile synthesis of worm-like micelles by visible light mediated dispersion polymerization. This approach begins ...

Synthesis and Performance Characterizations of Transition Metal Single Atom Catalyst for Electrochemical CO2 Reduction

Synthesis and Performance Characterizations of Transition Metal Single Atom Catalyst for Electrochemical CO2 Reduction

Síntese e caracterizações de desempenho de catalisador de átomo de Metal de transição para redução de eletroquímica CO2

This protocol presents both the synthesis method of the Ni single atom catalyst, and the electrochemical testing of its catalytic activity and selectivity ...

Preparation of Functional Silica Using a Bioinspired Method

Preparação de sílica funcional usando um método Bioinspirada

The goal of the protocols described herein is to synthesize bioinspired silica materials, perform enzyme encapsulation therein, and partially or totally ...

Heterogeneous Removal of Water-Soluble Ruthenium Olefin Metathesis Catalyst from Aqueous Media Via Host-Guest Interaction

Remoção heterogênea de rutênio Water-soluble Olefin Metathesis catalisador de meios aquosos através da interação do anfitrião-convidado

A highly efficient transition metal catalyst removal method is developed. The water-soluble catalyst contains a newly-designed NHC ligand for the catalyst ...

Imaging Corrosion at the Metal-Paint Interface Using Time-of-Flight Secondary Ion Mass Spectrometry

Corrosão por imagem na interface metal-Paint usando espectrometria de massa de íons secundários de tempo de voo

Corrosion developed at the paint and aluminum (Al) metal-paint interface of an aluminum alloy is analyzed using time-of-flight secondary ion mass ...

[(DPEPhos)(bcp)Cu]PF6: A General and Broadly Applicable Copper-Based Photoredox Catalyst

[DPEPhosbcpCu]PF6: A General and Broadly Applicable Copper-Based Photoredox Catalyst

[(DPEPhos) (BCP) UC] PF6: um catalisador Photoredox cobre-baseado geral e amplamente aplicável

Our group recently reported the use of [(DPEPhos)(bcp)Cu]PF6 as a general copper-based photoredox catalyst which proved efficient to promote the ...

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

Espectro vibracional de uma interface N719-Cromosforre/Titania da Simulação Empírica-Potencial de Dinâmica Molecular, solvated por um 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 ...

Development of Heterogeneous Enantioselective Catalysts using Chiral Metal-Organic Frameworks (MOFs)

Development of Heterogeneous Enantioselective Catalysts using Chiral Metal-Organic Frameworks MOFs

Desenvolvimento de catalisadores enantioseseseseses heterogêneos usando estruturas quirais metal-orgânicas (MOFs)

Substrate size discrimination by the pore size and homogeneity of the chiral environment at the reaction sites are important issues in the validation of ...