Reactive Force Field for P/Sn/I System: Atomistic Insight into the Early Stage of Black Phosphorus and Phosphorene Synthesis Process

Resumo Técnico: Campo de Força Reativo para o Sistema P/Sn/I e Visão Atomística da Síntese de Fósforo Negro

Definição do Problema
O fósforo negro (BP) e sua contraparte bidimensional, o fosforeno, são materiais promissores para a próxima geração de eletrônicos devido à sua alta mobilidade de portadores e band gaps ajustáveis. No entanto, a síntese de filmes de BP de alta qualidade e grande área via Transporte Químico de Vapor (CVT) continua sendo um desafio. O processo, que tipicamente utiliza fósforo vermelho com aditivos de estanho (Sn) e iodo (I₂), não é totalmente compreendido em nível atomístico. Persistem incertezas fundamentais em relação aos papéis específicos dos aditivos, a identidade das espécies estáveis em fase gasosa, os mecanismos de transporte de fósforo e a natureza dos intermediários de pré-nucleação. A observação experimental dessas dinâmicas iniciais é difícil devido às altas temperaturas e ambientes reativos, enquanto a Teoria do Funcional da Densidade (DFT) padrão é limitada a sistemas pequenos e configurações estáticas a 0K, falhando em capturar as redes de reação dinâmicas e de temperatura finita necessárias para modelar os mecanismos de crescimento.

Metodologia
Para abordar essas limitações, os autores desenvolveram um novo campo de força reativo (ReaxFF) especificamente para o sistema P/Sn/I, designado como P/Sn/I-2512. Este campo de força foi treinado contra um conjunto de dados abrangente de mecânica quântica (QM) gerado por DFT (funcional PBE com base TZ2P e correção de dispersão D3 de Grimme). O conjunto de treinamento incluiu:

• Energias de dissociação de ligação para ligações Sn–Sn, Sn–I, Sn–P, I–P e I–I.
• Perfis de ângulos de valência e torção para várias configurações moleculares.
• Calores de formação para espécies-chave (SnI₂, SnI₄, PI₃).
• Equações de estado para o estanho em fase condensada ($\alpha$-Sn e $\beta$-Sn) e energias de formação de clusters de estanho (Sn$_x$, $x=2-6$).

Os parâmetros do campo de força foram otimizados usando a Estratégia Evolutiva de Adaptação de Matriz de Covariância (CMA-ES) dentro do Amsterdam Modeling Suite. A validação foi realizada comparando os resultados do ReaxFF com os dados de DFT para dissociação de ligação, equações de estado e energias de reação, mostrando forte concordância.

Subsequentemente, simulações de Dinâmica Molecular Reativa (ReaxFF-MD) foram conduzidas usando LAMMPS para investigar os processos de pré-nucleação em estágio inicial a 923 K. O estudo empregou o ensemble NVT com um termostato de Nosé-Hoover. As simulações foram estruturadas para isolar variáveis específicas:

1. Validação do Sistema: Simulações de MD dos sistemas Sn/I₂ e P/I₂ em diferentes densidades (pressões efetivas de 0,01, 0,1 e 1 atm) para verificar a formação de espécies-chave como SnI₂, SnI₄, PI₂ e PI₃.
2. Papel do Iodo: Três sistemas foram comparados (apenas P; P + Sn; P + Sn + I₂) para isolar o efeito do iodo na recombinação do fósforo.
3. Efeitos de Densidade: O sistema P/Sn/I₂ foi simulado sob três regimes de densidade (pressões efetivas de 5, 10 e 100 atm) para modelar a transição de vapor diluído para fases condensadas, mimetizando as condições onde a formação de BP é experimentalmente observada para ocorrer.

Principais Resultados

• Validação do Campo de Força: O campo de força P/Sn/I-2512 reproduziu com precisão os dados de referência de DFT para dissociação de ligação, ângulos de valência, barreiras de torção e as equações de estado para cristais e clusters de estanho. Ele previu com sucesso a formação de SnI₂, SnI₄ e PI₃ em altas temperaturas.
• Dependência de Densidade em Sn/I₂ e P/I₂: No sistema Sn/I₂, o aumento da densidade deslocou a composição de espécies elementares para ambientes de coordenação complexos, promovendo a formação de SnI₄ e Sn₂I₄ (dímeros de SnI₂) além de SnI₂. No sistema P/I₂, densidades mais altas facilitaram a formação de PI₃ e espécies poliméricas de PI$_x$, consistente com as dificuldades experimentais em isolar o PI₃ puro.
• Papel do Iodo na Recombinação: As simulações revelaram que, na ausência de iodo, os átomos de fósforo permaneciam majoritariamente atômicos ou formavam clusters mínimos. A adição de I₂ aumentou significativamente a taxa de recombinação P–P, levando à formação de motivos binários de P$_x$I$_y$ e compostos ternários transitórios de Sn$_x$P$_y$I$_z$. Isso confirma o papel do iodo na solubilização do fósforo e na facilitação do transporte.
• Via de Pré-nucleação em Alta Densidade: Em baixas densidades (5 atm), o sistema consistia principalmente de fósforo atômico e pequenos clusters (P₂, P₃). Em densidade média (10 atm), clusters maiores (P₄, P₅) começaram a se formar. Em alta densidade (100 atm), o sistema condensou-se rapidamente em clusters de ordem superior de P$_x$ (x = 8–12) e grandes compostos ternários de Sn$_x$P$_y$I$_z$.
• Identificação de Intermediários: Na densidade mais alta, as simulações produziram estruturas de fósforo do tipo Hittorf condensadas nas bordas dos clusters de Sn$_x$P$_y$I$_z$. Além disso, motivos hexagonais do tipo BP estabilizados pelo iodo foram observados, identificados como potenciais sementes de nucleação.

Significância e Alegações
O artigo alega fornecer um modelo atomístico prático para estudar as condições de crescimento e os efeitos de aditivos na síntese de vapor de BP/fosforeno. As principais contribuições são:

1. Desenvolvimento de um Campo de Força Validado: O campo de força P/Sn/I-2512 oferece uma ferramenta confiável para simular sistemas reativos complexos envolvendo fósforo, estanho e iodo, superando as limitações de tamanho e escala de tempo da DFT.
2. Visão Mecanística: O estudo elucida que o iodo e a densidade do sistema são os fatores de controle para a recombinação do fósforo. Demonstra que o iodo não é apenas um agente de transporte, mas promove ativamente a recombinação P–P através da formação de intermediários P$_x$I$_y$ e compostos ternários de Sn$_x$P$_y$I$_z$.
3. Modelo de Pré-nucleação: Os autores propõem uma via específica de pré-nucleação onde condições de alta densidade, combinadas com o iodo, facilitam a condensação do fósforo em intermediários do tipo Hittorf e hexágonos do tipo BP. Essas espécies são sugeridas como precursores para o subsequente crescimento do fósforo negro durante a etapa de resfriamento.

O trabalho não propõe novos protocolos experimentais, mas oferece um arcabouço teórico para interpretar observações experimentais existentes, especificamente a necessidade de alta pressão (densidade) e o papel crítico do iodo na síntese por CVT de fósforo negro.