Tuning Terahertz Optomechanics of MoS2 Bilayers with Homogeneous In-plane Strain

Resumen Técnico: Sintonización de la Optomecánica de Terahertzos de Bicapas de MoS2 con Deformación In-plane Homogénea

Problema y Motivación
Si bien los materiales de van der Waals (vdW) son explorados activamente para la generación o detección de radiación terahertz (THz), sus vibraciones de fonones Raman-activos en el régimen THz ofrecen una plataforma paralela para el control optomecánico. Una pregunta fundamental sigue siendo cómo la deformación mecánica in-plane se propaga a través de monocapas en sistemas de bicapa (BL) para afectar el acoplamiento interlaminar. Estudios previos han examinado respuestas intralayer y deformación compresiva isotrópica, pero la optomecánica del acoplamiento vdW interlaminar bajo un acoplamiento de deformación anisotrópico in-plane y out-of-plane apenas ha comenzado a ser explorada recientemente. Comprender este acoplamiento cruzado es crítico para el diseño de heteroestructuras de vdW donde las interacciones interlayered dictan fenómenos que van desde excitones de moiré hasta la superconductividad.

Metodología
Los autores sintetizaron bicapas de MoS2 con apilamiento hexagonal (2H) y romboédrico (2R) mediante deposición física de vapor (PVD) de alta temperatura sobre sustratos de SiO2 y Si3N4. La alta temperatura de crecimiento (1200 °C) combinada con el desajuste de la expansión térmica entre el MoS2 y los sustratos induce una deformación de tracción biaxial in-plane homogénea. Los niveles de deformación se estimaron en aproximadamente +0.36% en Si3N4 y +0.72% en SiO2.

Para caracterizar las muestras, los autores emplearon:

• Generación de Segundo Armónico (SHG): Utilizada para distinguir entre los órdenes de apilamiento 2H y 2R basados en la simetría (el 2H carece de un centro de inversión en el bulto pero muestra una SHG débil debido a la ruptura de la interfaz; el 2R muestra una SHG fuerte) y para confirmar la naturaleza biaxial de la deformación mediante la dependencia angular.
• Espectroscopía Raman: Se registraron espectros tanto de alta frecuencia (modos $E^1_{2g}$ e $A_{1g}$ intralayer) como de ultra baja frecuencia (modos de cizalladura interlaminar $E^2_{2g}$ y de respiración $B_{2g}$ cerca de 1 THz).
• Cálculos de Primeros Principios: Se realizaron cálculos de la Teoría del Funcional de la Densidad (DFT) utilizando el funcional optPBE-vdW y correcciones DFT-D3 para modelar las dispersiones de fonones, las actividades Raman y la respuesta geométrica del espaciado interlaminar bajo deformación.

Resultos Clave

1. Endurecimiento de los Modos Interlayer inducido por Deformación: El hallazgo más significativo es el endurecimiento del modo de respiración interlaminar ($B_{2g}$) bajo deformación de tracción in-plane. A medida que la deformación biaxial aumentó de 0% (muestra relajada, transferida) a +0.72% (crecida sobre SiO2), la frecuencia de $B_{2g}$ aumentó de 1.13 THz a 1.21 THz. Esto arroja una tasa de sintonización por deformación experimental de +0.11 THz/% (+3.7 cm⁻/%).
2. Mecanismo vía Contracción de Poisson: Este endurecimiento es contraintuitivo para osciladores armónicos estándar, pero se explica por el efecto Poisson. La deformación de tracción in-plane provoca una contracción del espaciado interlaminar (brecha vdW) para conservar el volumen. Esta contracción geométrica vuelve más pronunciado el potencial de confinamiento entre las capas, aumentando así la frecuencia de la oscilación out-of-plane.
3. Relación de Poisson Efectiva: Al correlacionar el desplazamiento de frecuencia experimental con la contracción interlaminar teórica, los autores derivaron una relación de Poisson out-of-plane efectiva ($\nu_{eff}$) para las bicapas de aproximadamente 0.19–0.24. Este valor es inferior al reportado para el MoS2 en bulto (~0.30) y se alinea con las predicciones modernas de DFT con funcionales híbridos.
4. Parámetro de Grüneisen Gigante: El sistema exhibe un parámetro de Grüneisen out-of-plane ($\gamma_{out}$) excepcionalmente grande de aproximadamente 14–20. Este valor supera incluso el valor "gigante" reportado para la fosforeno ($\gamma \approx 8.6$), lo que indica una sensibilidad extrema de la curvatura del potencial vdW a la compresión geométrica.
5. Simetría y Apilamiento: El estudio confirmó que el modo de cizalladura $E^2_{2g}$ no se divide, lo cual es consistente con una deformación biaxial uniforme que preserva la simetría hexagonal. Se encontró que las intensidades relativas de los modos $E^2_{2g}$ y $B_{2g}$ se invierten entre los órdenes de apilamiento 2H y 2R, reflejando sus diferentes simetrías.

Significancia y Reivindicaciones
Este artículo establece que la deformación de tracción biaxial in-plane homogénea puede utilizarse para sintonizar finamente la respuesta optomecánica de terahertzos de las bicapas de MoS2 sin la necesidad de perillas de presión externas. Los autores afirman que el acoplamiento vertical suave en estas bicapas de vdW actúa como un amplificador sensible de la deformación mecánica in-plane.

Los hallazgos redefinen el espaciado vdW en sistemas de bicapa como un grado de libertad dinámico y altamente sintonizable en el régimen de THz. Los autores postulan que el endurecimiento del enlace interlaminar implica una modulación directa del parámetro de salto electrónico interlaminar ($t_\perp$). Este acoplamiento electromecánico proporciona una base física para:

• Manipular dinámicamente la hibridación de capas de los excitones (controlando los tiempos de vida y la coherencia de valle).
• Ofrecer un control determinista para sintonizar el ancho de banda de las bandas planas de moiré en heteroestructuras torcidas sin alterar el ángulo de torsión.

El trabajo posiciona a estas bicapas de van der Waals como plataformas mecánicas no lineales altamente sintonizables direccionables en el régimen de THz, con aplicaciones potenciales en filtros acústicos de microondas y moduladores acusto-ópticos.