Molecular Hyperpolarisability as a Screening Descriptor for Second-Order Nonlinear Optics in Ferroelectric Nematic Liquid Crystals

Resumo Técnico: Hiperpolarizabilidade Molecular como Descritor de Triagem para Óptica Não Linear de Segunda Ordem em Cristais Líquidos Ferroelétricos Nemáticos

Definição do Problema
Os cristais líquidos ferroelétricos nemáticos ($N_F$) oferecem uma combinação única de fluidez e ordem polar macroscópica, tornando-os candidatos atraentes para aplicações ópticas não lineares (NLO) de segunda ordem, como a geração de segunda harmônica (SHG). Embora o material nemático $N_F$ arquetípico, RM734, exiba grandes coeficientes de óptica não linear, a maioria dos materiais nemáticos ferroelétricos conhecidos foi projetada para ordenamento polar, e não para desempenho NLO. Consequentemente, as características estruturais que maximizam a hiperpolarizabilidade molecular ($\beta$) dentro das restrições da organização de cristal líquido polar ainda não estão bem estabelecidas. Este hiato limita a descoberta racional de novos materiais $N_F$ otimizados para NLO. O desafio central abordado é se cálculos de estrutura eletrônica podem rastrear de forma confiável materiais de cristal líquido polar conhecidos para identificar candidatos promissores de NLO e estabelecer regras de design, apesar das complexidades de converter propriedades moleculares microscópicas em coeficientes macroscópicos.

Metodologia
Os autores desenvolveram um fluxo de trabalho computacional para fazer a ponte entre a estrutura eletrônica molecular e os coeficientes NLO macroscópicos ($d_{ij}$):

1. Cálculos de Estrutura Eletrônica:

   • Geometria: As geometrias iniciais foram otimizadas usando o método semiemperical GFN2-xTB. Conformeros foram gerados para capturar 80% da população cumulativa de Boltzmann.
   • Hiperpolarizabilidade ($\beta$): Tensores de hiperpolarizabilidade molecular dependentes da frequência foram calculados a 800 nm utilizando uma gama de métodos de Teoria do Funcional da Densidade (DFT) (HF, híbrido GGA, híbrido meta-GGA e híbridos de intervalo separado) combinados com várias bases (6-31+G(d), 6-311++G(d,p), aug-cc-pVDZ, aug-cc-pVTZ).
   • Média de Conformeros: O estudo comparou cálculos em geometrias de mínimo global único contra conjuntos de conformeros ponderados por Boltzmann para determinar se o escalonamento explícito de médias melhorava a precisão.

2. Conversão Macroscópica (Modelo de Gás Orientado):

   • Os tensores moleculares foram convertidos em coeficientes $d$ macroscópicos usando um modelo de gás orientado.
   • Isso envolveu a média térmica dos componentes do tensor usando parâmetros de ordem orientacional ($\langle P_1 \rangle$ e $\langle P_3 \rangle$) derivados de dados experimentais de polarização e birrefringência para materiais de referência (RM734, DIO, C1).
   • Correções de campo local foram aplicadas usando a aproximação do índice de refração efetivo de Vuks ($n_{eff}$), e a densidade numérica foi estimada com base em uma densidade representativa de 1,3 g cm⁻³.

3. Benchmarking e Triagem:

   • O fluxo de trabalho foi testado contra os coeficientes experimentais $d_{22}$, $d_{15}$ e $d_{31}/d_{32}$ para RM734, DIO e C1.
   • As métricas de desempenho incluíram o viés (erro percentual médio assinado) e a dispersão (coeficiente de variação) das razões entre o calculado e o experimental.
   • O melhor protocolo de desempenho foi aplicado a um conjunto de dados curados de cristais líquidos polares conhecidos ($N_F$, $SmA_F$, $SmC_F$) para triagem de respostas NLO aprimoradas.

Principais Resultados

• Seleção de Método: Nenhum método único reproduziu perfeitamente os valores absolutos, pois as magnitudes dependem fortemente dos parâmetros de volume assumidos (densidade, parâmetros de ordem). No entanto, as tendências relativas foram robustas.

  • Melhor Desempenho: A combinação wB97X-D/aug-cc-pVDZ foi selecionada como o protocolo de triagem ideal devido à menor dispersão nas razões previstas entre os materiais de referência.
  • Bases (Basis Sets): O aumento do tamanho da base de dupla-zeta (aug-cc-pVDZ) para tripla-zeta (aug-cc-pVTZ) não melhorou significamente a precisão, justificando o uso das bases de dupla-zeta menos dispendiosas para triagem de alto rendimento.
  • Média de Conformeros: A média explícita sobre conjuntos de conformeros não melhorou consistentemente a concordância com o experimento em comparação com o uso da geometria de mínimo global único. Para fins de triagem, cálculos de conformero único foram considerados suficientes.

• Regras de Design Identificadas:

  • Assimetria Doador-Aceitador: Os maiores valores de $d_{22}$ previstos ocorrem em sistemas com sistemas eletrônicos "push-pull" fortes. Doadores mais fortes (ex: N,N-dimetilamino, metiltio) e aceitadores fortes (ex: nitro) correlacionam-se com coeficientes mais altos.
  • Caminhos de Conjugação: Caminhos conjugados estendidos e planos ao longo do eixo longo molecular são críticos. Ligações acetilênicas superam significativamente as de ésteres.
  • Prejuízo Estrutural: Motivos que promovem a ordem polar, mas interrompem a conjugação, como a difluorometilenooxi (ponte Q) ou 1,3-dioxanos 2,5-dissubstituídos, resultam em baixas respostas NLO.
  • Grupos Laterais: Embora o flúor aumente a estabilidade térmica da ordem polar, sua presença como grupo lateral frequentemente reduz os coeficientes NLO calculados em comparação com análogos não fluorados.

• Resultados da Triagem: A aplicação do fluxo de trabalho a um conjunto de dados de materiais $N_F$ conhecidos identificou vários candidatos com valores de $d_{22}$ previstos excedendo o do RM734 (calculado como ~29,7 pm/V vs. experimental de 22,5 pm/V). Os valores mais altos previstos aproximaram-se de ~49,7 pm/V.

Significância e Alegações
O artigo afirma estabelecer um fluxo de trabalho computacionalmente eficiente e prático para triagem de cristais líquidos polares para aplicações de NLO de segunda ordem. Sua principal significância reside em demonstrar que:

1. Tendências relativas são mais confiáveis que valores absolutos: Embora os coeficientes $d$ absolutos sejam sensíveis a parâmetros de volume incertos (densidade, parâmetros de ordem), a classificação relativa de candidatos moleculares é robusta através de diferentes métodos DFT.
2. Protocolos simplificados são suficientes: A triagem de alto rendimento não requer bases de tripla-zeta caras ou a complexa média de conjuntos de conformeros; cálculos de geometria única com híbridos de intervalo separado (especificamente o wB97X-D) e bases de dupla-zeta fornecem métricas de triagem consistentes.
3. Design Racional: O estudo fornece regras de design concretas para a próxima geração de materiais $N_F$: maximizar a assimetria doador-aceitador e manter a conjugação estendida, evitando motivos estruturais que interrompam o sistema $\pi$.

Os autores mantêm uma postura modesta, observando que os coeficientes calculados servem como "métricas de triagem relativa" em vez de previsões precisas de desempenho macroscópico absoluto, uma vez que as propriedades de volume (densidade, parâmetros de ordem) variam entre os materiais e são difíceis de prever ab initio. O trabalho lança as bases para a descoberta orientada por dados de nemáticos ferroelétricos otimizados para NLO sem a necessidade de nova síntese experimental para cada candidato.