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

Resumen Técnico: La hiperpolarizabilidad molecular como descriptor de cribado para la óptica no lineal de segundo orden en cristales líquidos nemáticos ferroeléctricos

Planteamiento del Problema
Los cristales líquidos nemáticos ferroeléctricos ($N_F$) ofrecen una combinación única de fluidez y orden polar macroscópico, lo que los convierte en candidatos atractivos para aplicaciones ópticas no lineales (NLO) de segundo orden, como la generación de segundo armónico (SHG). Si bien el material $N_F$ arquetípico, RM734, exhibe grandes coeficientes de óptica no lineal, la mayoría de los nemáticos ferroeléctricos conocidos fueron diseñados para el ordenamiento polar más que para el rendimiento NLO. En consecuencia, las características estructurales que maximizan la hiperpolarizabilidad molecular ($\beta$) dentro de las restricciones de la organización de cristal líquido polar están insuficientemente establecidas. Este vacío limita el descubrimiento racional de nuevos materiales $N_F$ optimizados para NLO. El desafío central abordado es si los cálculos de estructura electrónica pueden cribar de manera fiable los materiales de cristal líquido polar conocidos para identificar candidatos prometedores de NLO y establecer reglas de diseño, a pesar de las complejidades de convertir las propiedades moleculares microscópicas en coeficientes macroscópicos.

Metodología
Los autores desarrollaron un flujo de trabajo computacional para conectar la estructura electrónica molecular con los coeficientes NLO macroscópicos ($d_{ij}$):

1. Cálculos de Estructura Electrónica:

   • Geometría: Las geometrías iniciales se optimizaron utilizando el método semiempírico GFN2-xTB. Se generaron confórmeros para capturar el 80% de la población Boltzmann acumulada.
   • Hiperpolarizabilidad ($\beta$): Se calcularon los tensores de hiperpolarizabilidad molecular dependientes de la frecuencia a 800 nm utilizando una gama de métodos de la Teoría del Funcional de la Densidad (DFT) (HF, híbrido GGA, híbrido meta-GGA y híbridos de rango separado) junto con varios conjuntos de bases (6-31+G(d), 6-311++G(d,p), aug-cc-pVDZ, aug-cc-pVTZ).
   • Promedio de Confórmeros: El estudio comparó los cálculos en geometrías de un solo mínimo global frente a conjuntos de confórmeros ponderados por Boltzmann para determinar si el promedio explícito mejoraba la precisión.

2. Conversión Macroscópica (Modelo de Gas Orientado):

   • Los tensores moleculares se convirtieron en coeficientes $d$ macroscópicos mediante un modelo de gas orientado.
   • Esto implicó el promedio térmico de los componentes del tensor utilizando parámetros de orden orientacional ($\langle P_1 \rangle$ y $\langle P_3 \rangle$) derivados de datos experimentales de polarización y birrefringencia para materiales de referencia (RM734, DIO, C1).
   • Se aplicaron correcciones de campo local utilizando la aproximación del índice de refracción efectivo de Vuks ($n_{eff}$), y la densidad numérica se estimó basándose en una densidad representativa de 1.3 g cm$^{-3}$.

3. Benchmarking y Cribado:

   • El flujo de trabajo se contrastó con los coeficientes experimentales $d_{22}$, $d_{15}$ y $d_{31}/d_{32}$ para RM734, DIO y C1.
   • Las métricas de desempeño incluyeron el sesgo (error porcentual medio con signo) y la dispersión (coeficiente de variación) de las relaciones calculadas-a-experimentales.
   • El protocolo con mejor desempeño se aplicó a un conjunto curado de cristales líquidos polares conocidos ($N_F$, $SmA_F$, $SmC_F$) para cribar la búsqueda de respuestas NLO mejoradas.

Resultados Clave

• Selección de Métodos: Ningún método único reprodujo perfectamente los valores absolutos experimentales, ya que las magnitudes dependían fuertemente de los parámetros de volumen asumidos (densidad, parámetros de orden); sin embargo, las tendencias relativas fueron robustas.

  • Mejor Desempeño: La combinación wB97X-D/aug-cc-pVDZ fue seleccionada como el protocolo de cribado óptimo debido a la menor dispersión en las relaciones predichas entre los materiales de referencia.
  • Conjuntos de Bases: Aumentar el tamaño del conjunto de bases de doble zeta (aug-cc-pVDZ) a triple zeta (aug-cc-pVTZ) no mejoró significamente la precisión, lo que justifica el uso de los conjuntos de doble zeta menos costosos para el cribado de alto rendimiento.
  • Promedio de Confórmeros: El promedio explícito sobre conjuntos de confórmeros no mejoró consistentemente el acuerdo con el experimento en comparación con el uso de la geometría de un solo mínimo global. Para fines de cribado, los cálculos de un solo confórmero se consideraron suficientes.

• Reglas de Diseño Identificadas:

  • Asimetría Donante-Aceptor: Los valores de $d_{22}$ predichos más altos ocurren en sistemas con sistemas electrónicos de tipo "push-pull" fuertes. Donantes más fuertes (p. ej., N,N-dimetilamino, metiltio) y aceptores fuertes (p. ej., nitro) se correlacionan con coeficientes más altos.
  • Vías de Conjugación: Las vías conjugadas, planas y extendidas a lo largo del eje largo molecular son críticas. Los conectores acetilénicos superan significativamente a los ésteres.
  • Detrimento Estructural: Los motivos que promueven el orden polar pero interrumpen la conjugación, como el difluorometilenooxi (puente Q) o los 1,3-dioxanos 2,5-disustituidos, resultan en bajas respuestas NLO.
  • Grupos Laterales: Aunque el flúor mejora la estabilidad térmica del orden polar, su presencia como grupo lateral a menudo reduce los coeficientes NLO calculados en comparación con los análogos no fluorados.

• Resultados del Cribado: La aplicación del flujo de trabajo a un conjunto de materiales $N_F$ conocidos identificó varios candidatos con valores de $d_{22}$ predichos que exceden al de RM734 (calculado como ~29.7 pm/V frente a los 22.5 pm/V experimentales). Los valores más altos predichos se acercaron a ~49.7 pm/V.

Significado y Reivindicaciones
El artículo afirma establecer un flujo de trabajo computacional práctico y eficiente para cribar cristales líquidos polares para aplicaciones de NLO de segundo orden. Su principal significado radica en demostrar que:

1. Las tendencias relativas son más fiables que los valores absolutos: Aunque los coeficientes $d$ absolutos son sensibles a los parámetros de volumen inciertos (densidad, parámetros de orden), el ranking relativo de los candidatos moleculares es robusto a través de diferentes métodos DFT.
2. Los protocolos simplificados son suficientes: El cribado de alto rendimiento no requiere conjuntos de bases triple zeta costosos ni complejos promedios de conjuntos de confórmeros; los cálculos de geometría única con híbridos de rango separado (específicamente wB97X-D) y conjuntos de bases de doble zeta proporcionan métricas de cribado consistentes.
3. Diseño Racional: El estudio proporciona reglas de diseño concretas para la próxima generación de materiales $N_F$: maximizar la asimetría donante-aceptor y mantener una conjugación extendida evitando simultáneamente los motivos estructurales que interrumpen el sistema $\pi$.

Los autores mantienen una postura modesta, señalando que los coeficientes calculados sirven como "métricas de cribado relativas" en lugar de predicciones precisas del rendimiento macroscópico absoluto, ya que las propiedades de volumen (densidad, parámetros de orden) varían entre materiales y son difíciles de predecir ab initio. El trabajo sienta las bases para el descubrimiento impulsado por datos de nemáticos ferroeléctricos optimizados para NLO sin requerir una nueva síntesis experimental para cada candidato.