Modeling Compressive Instability in Two-Dimensional Ti2COx MXenes

Resumo Técnico: Modelagem de Instabilidade de Compressão em MXenes Ti2COx Bidimensionais

Definição do Problema
Embora as propriedades mecânicas de MXenes sob carga de tração tenham sido extensivamente caracterizadas, seu comportamento sob cargas de compressão e cisalhamento permanece amplamente inexplorado. Em aplicações práticas, como eletrônica flexível, dispositivos vestíveis e compósitos de MXene-polímero, as nanofolhas encontram frequentemente tensões de compressão devido à flexão, retração da matriz ou carga externa. Este ambiente compressivo pode levar ao enrugamento (buckling), ondulação ou falha estrutural. A literatura atual carece de uma compreensão fundamental dos mecanismos de instabilidade compressiva de MXenes, particularmente no que diz respeito aos efeitos de terminações superficiais, defeitos atômicos e fenômenos dependentes do tamanho que a mecânica do contínuo clássica pode não ser capaz de capturar.

Metodologia
O estudo emprega simulações de dinâmica molecular (MD) reativa usando o Large Scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) e o campo de força ReaxFF (parametrizado em [36]) para investigar o comportamento de compressão e pós-enrugamento de nanofolhas de Ti2C e Ti2CO2.

• Configuração da Simulação: Nanofolhas com dimensões intraplanares de 14×14 nm² foram modeladas a 1 K para isolar respostas mecânicas intrínsecas. As simulações utilizaram um passo de tempo (timestep) de 0,2 fs.
• Condições de Carga: O estudo examinou compressão uniaxial e biaxial, bem como carga de cisalhamento, ao longo das direções cristalográficas armchair e zigzag.
• Variáveis Investigadas: A pesquisa analisou sistematicamente os efeitos de:
  • Taxas de deformação (variando de $0,725 \times 10^6$ a $2,85 \times 10^8$ s⁻¹).
  • Tipos de defeitos: Defeitos pontuais isolados (vagas de Ti e C), aglomerados de vagas aleatórias e aglomerados específicos de vagas de Ti.
  • Pressão de confinamento lateral (modelada via molas harmônicas/fundação de Winkler).
  • Terminação superficial (comparando Ti2C prístino com Ti2CO2 terminado em oxigênio).
• Validação e Modelagem: O potencial ReaxFF foi validado contra Teoria do Funcional da Densidade (DFT) e dados experimentais para parâmetros estruturais e módulo de Young. Para abordar as discrepâncias entre resultados atomísticos e a teoria clássica, um modelo de contínuo não local impulsionado por tensão foi implementado e calibrado com os dados de MD.

Principais Resultados

1. Taxa de Deformação e Anisotropia: Taxas de deformação mais altas resultaram em maiores tensões e deformações de enrugamento, acompanhadas por um comportamento de pós-enrugamento mais instável, caracterizado por maiores oscilações de tensão. As nanofolhas de Ti2C exibiram maior resistência ao enrugamento ao longo da direção armchair em comparação com a direção zigzag.
2. Influência de Defeitos:
   • Defeitos Pontuais Isolados: Com até 3% de fração de defeitos, as vagas isoladas reduziram primariamente a tensão crítica de enrugamento e aumentaram as oscilações de tensão de pós-enrugamento, mas não alteraram significativamente a deformação de enrugamento ou os modos globais de forma. As vagas de Ti tiveram um efeito de enfraquecimento mais pronunciado do que as de C.
   • Aglomerados de Vagas: Aglomerados de vagas de Ti (especificamente grupos de 4–5 vagas) alteraram fundamentalmente o mecanismo de deformação. Em vez do enrugamento clássico de placa fina, as nanofolhas exibiram deflexões fora do plano, onduladas e localizadas, impulsionadas por uma severa complacência localizada ao redor dos aglomerados de defeitos.
3. Terminação Superficial e Confinamento:
   • Terminação de Oxigênio: A presença de grupos superficiais –O aumentou substancialmente a tensão de enrugamento (de ~1 GPa para ~3,5 GPa) e reduziu a anisotropia entre as direções de carga. No entanto, as nanofolhas de Ti2CO2 apresentaram fratura em grandes deformações compressivas, enquanto o Ti2C prístino manteve a integridade estrutural até deformações superiores a 0,35.
   • Confinamento Lateral: A aplicação de pressão de confinamento lateral (simulando uma matriz elástica) aumentou significativamente tanto a tensão quanto a deformação de enrugamento, além de amortecer as oscilações de pós-enrugamento.
4. Modos de Enrugamento: Sob compressão uniaxial, as nanofolhas exibiram primeiros modos de enrugamento do tipo placa. Sob compressão biaxial, observou-se um formato de enrugamento simétrico, do tipo domo. A carga de cisalhamento produziu modos de deflexão elípticos.
5. Contínuo vs. Atomístico: A teoria clássica de vigas de Euler-Bernoulli subestimou significadamente as deformações críticas de enrugamento das nanofolhas. Uma formulação de contínuo não local impulsionada por tensão, calibrada com um parâmetro de comprimento característico de 5 Å, reproduziu com sucesso as deformações de enrugamento dependentes do tamanho observadas nas simulações de MD.

Significância e Alegações
O artigo afirma preencher uma lacuna crítica na literatura ao fornecer a primeira investigação atomística sistemática da instabilidade de compressão e cisalhamento de nanofolhas de MXene. O estudo demonstra que:

• A mecânica do contínuo clássica é insuficiente para prever o comportamento de enrugamento de MXenes 2D devido aos efeitos de tamanho, necessitando de formulações não locais para uma modelagem precisa.
• A estabilidade mecânica de MXenes é altamente sensível à natureza específica dos defeitos (isolados vs. aglomerados) e à química de superfície.
• O confinamento lateral e a terminação superficial são alavancas eficazes para aumentar a estabilidade compressiva.

Os autores concluem que esses achados fornecem uma compreensão fundamental das transformações morfológicas de MXenes sob compressão. Esse conhecimento é apresentado como essencial para o design de compósitos de MXene robustos e eletrônica flexível, e sugere-se que esses insights possam estimular pesquisas futuras em nanotubos, nanoscrolls e arquiteturas dobradas baseadas em MXene. O estudo enfatiza que as tendências observadas são insights mecanísticos esperados de serem qualitativamente transferíveis para condições de menor taxa de deformação, apesar das altas taxas de deformação inerentes às simulações de MD.