SciELO - Scientific Electronic Library Online

 
vol.43 número2Análisis y Simulación del Fenómeno de Dispersión por Modo de Polarización (PMD) en Fibras Ópticas ITU-T G.655 índice de autoresíndice de materiabúsqueda de artículos
Home Pagelista alfabética de revistas  

Servicios Personalizados

Revista

Articulo

Indicadores

Links relacionados

  • No hay articulos similaresSimilares en SciELO

Compartir


Revista Politécnica

versión On-line ISSN 2477-8990versión impresa ISSN 1390-0129

Rev Politéc. (Quito) vol.43 no.2 Quito may./jul. 2019

 

Articles

Propiedades Ópticas no lineales en Señales de Mezcla de Cuatro Ondas: Consideraciones Estocásticas del Solvente y Difusión Espectral

Nonlinear Optical Properties in Four-Wave Mixing Signals: Stochastic Considerations of Solvent and Spectral Diffusion

José Luis Paz 1   *  

Luis Lascano 1  

Cesar Costa-Vera 1  

1Escuela Politécnica Nacional, Departamento de Física, Quito, Ecuador


Resumen:

Usando las ecuaciones Ópticas de Bloch y considerando un sistema molecular que interactúa con un reservorio térmico, derivamos expresiones analíticas para las propiedades ópticas no lineales en la espectroscopia de mezcla de cuatro ondas. En este trabajo se introducen diferentes metodologías para el estudio de las polarizaciones macroscópicas y susceptibilidades ópticas no lineales. En un sistema de dos niveles, hemos considerado una frecuencia dependiente del tiempo para describir los efectos colectivos del reservorio térmico, cuya manifestación es el ensanchamiento del nivel superior según una función aleatoria prescrita. Aquí, las respuestas ópticas son calculadas mediante una metodología basada en el promedio de la coherencia obtenida por la resolución de las ecuaciones de Bloch ópticas estocásticas. Evaluamos las susceptibilidades ópticas a diferentes frecuencias, representadas como superficies numéricas para el comportamiento de las respuestas ópticas, en función del factor de desintonización de frecuencias bombeo-prueba, y de la intensidad del ruido. Se estudian los tiempos de relajación utilizando las propiedades de simetría de las respuestas ópticas y los efectos de la difusión espectral en la susceptibilidad no lineal de tercer orden para sistemas cuánticos de dos y tres estados. Los resultados presentados en este trabajo muestran aplicaciones potenciales para la caracterización de sistemas moleculares complejos a través de sus propiedades ópticas no lineales.

Palabras clave: susceptibilidad óptica; propiedades absortivas y dispersivas; mezcla de cuatro ondas

Abstract:

Using the Optical Bloch equations, and considering a molecular system interacting with a thermal reservoir in Four-wave mixing spectroscopy, we derive analytical expressions for the nonlinear optical properties. In this work, we introduce different methodologies for the study of the macroscopic polarizations and nonlinear optical susceptibilities. In a two-level system, we have considered the collective effects of the thermal reservoir as a time dependent frequency, whose manifestation is the broadening of the upper level according to a prescribed random function. Here, the optical responses calculated through a methodology based on the averaging of the coherence obtained by the resolution of the optical stochastic Bloch equations. We were able to calculate the optical susceptibilities at different frequencies, depicted as numerical surfaces for the behavior of the optical responses, as function of the pump- and probe-frequency detuning factor and the strength of the noise source. Relaxation times using symmetry properties of the optical responses, effects of spectral diffusion on the third order nonlinear susceptibility for two- and three-level quantum systems, are studied. The results presented in this work show potential applications for the optical characterization of complex molecular systems.

Keywords: Optical susceptibility; absorptive and dispersive properties; four-wave mixing

1. INTRODUCCIÓN

Las técnicas espectroscópicas de mezcla de cuatro ondas (FWM por sus siglas en inglés) degenerada y no-degenerada han sido empleadas de muy variadas maneras en las últimas décadas. Destacamos como esta técnica ha sido utilizada para estudiar polímeros (Singh y col. 1988), moléculas biológicas (Sugisaki y col. 2007), clusters de semiconductores de tamaño nanométrico (Wang, 1991), excitones (Abramavicius y col., 2007), siendo particularmente útil en la descripción de los datos no lineales cuando el sistema molecular tiene una fuerte interacción con su medio en presencia de campos electromagnéticos. Más recientemente, Al-Saidi y Abdulkareem (2016) han reportado cambios en las propiedades ópticas de una solución química, debido a la concentración molecular del colorante estudiado. En este caso, la respuesta del sistema al cambio de la concentración molecular se produce por la elevada absorción saturable no lineal inducida, de importancia para el desarrollo de dispositivos no lineales, como limitadores ópticos. Por otro lado, Sudheesh y col. (2013) han reportado un estudio sobre el tercer orden no lineal, propiedades ópticas de una nueva clase de fenilhidrazonas y la influencia de las nanopartículas de plata y oro en su respuesta no lineal. Experimentos de FWM han sido muy útiles para extraer información espectroscópica de los tiempos de relajación longitudinal y transversal (Yajima y Souma, 1978; Souma et al, 1980). En estos experimentos, las ondas incidentes que se propagan interactúan con el medio molecular, generando una señal no lineal. En mezcla de cuatro ondas (FWM), dos haces láser (campo de bombeo con frecuencia ω1 y campo de prueba con frecuencia ω2) se enfocan sobre el medio resonante, generando una señal con frecuencia ω3=2ω1-ω2 y vector de onda k32k1-k2, donde los subíndices 1,2 y 3 se refieren a los haces de bombeo, prueba y señal generada, respectivamente. Todas las frecuencias ópticas consideradas son cercanas a la frecuencia molecular de interés (García-Golding y Marcano, 1985; Paz y col., 1988). Diferentes estudios teóricos se han desarrollado en relación con esta técnica. Se han derivado formalismos generales para incluir el efecto Stark dinámico aplicado a muchos problemas, en particular, en la comprensión de las inestabilidades de los campos electromagnéticos a través de los medios de absorción (Boyd y Mukamel, 1984). Muchos fenómenos de interacción radiación-materia, se notan cuando los campos electromagnéticos pasan a través de una determinada trayectoria óptica (Paz y col., 2018). Se han desarrollado algunos métodos que utilizan procedimientos algebraicos Wei-Norman para el estudio de la propagación de campos a través de sistemas homogéneos, con una analogía entre el operador de evolución y las llamadas matrices de propagación óptica. Los estudios sobre el efecto de la propagación de los campos en una distribución espectral homogénea de la señal FWM, han considerado modelos en los que la intensidad del bombeo se considera constante a lo largo del trayecto óptico (Agarwal y Nayak, 1984; Boyd y col, 1981). En líquidos y cristales, muchas transiciones son ensanchadas de forma no homogénea mediante distribuciones aleatorias locales en los sitios de absorción. En otros casos en los que el mecanismo de relajación implica cambios en la frecuencia de transición del sistema molecular, las ecuaciones que gobiernan la evolución temporal del vector de Bloch se definen como ecuaciones de Bloch ópticas estocásticas (OSBE por sus siglas en inglés). La mayor parte del trabajo realizado en óptica no lineal con las ecuaciones de Bloch Ópticas Convencionales (OCBE por sus siglas en inglés) considera como parámetro la frecuencia de transición del sistema de dos niveles. Sin embargo, si uno está interesado en la generación de radiación a partir de un sistema de dos niveles impulsado por un campo monocromático e inmerso en un medio condensado (reservorio térmico), la suposición ya no es válida. En esas circunstancias, las colisiones moleculares inducen desplazamientos de frecuencia y debido a la naturaleza de las colisiones, conducen a la modulación estocástica de la frecuencia central. Este cambio en la frecuencia de transición viene dado por ξ(t)=ω_0+σ(t) , donde el término σ(t) comprende toda la estocasticidad del proceso, y ω_0es la frecuencia de transición del sistema molecular aislado de dos niveles. Estas fluctuaciones han sido tratadas por diferentes teorías: El modelo gaussiano Markoviano (Brewer y De Voe, 1984), los modelos no Markovianos (Apanasevich, 1984) y un "modelo telegráfico aleatorio" (Wodkiewics y Eberly, 1985).

En la mayoría de los casos, los efectos radiativos no son ciertamente la única fuente y, por lo general, no son la fuente dominante de fluctuaciones e incoherencias de relajación. Si el mecanismo de relajación implica cambios en la frecuencia de transición del sistema molecular que persisten tras una interacción entre la molécula-activa y la perturbación electromagnética, las ecuaciones que rigen la evolución temporal son de naturaleza integro-diferencial y las OCBE ya no son válidas. Normalmente la frecuencia de transición en el par de estados moleculares se considera como un parámetro. Esta suposición es válida sólo en casos de interacción débil entre solutos y disolventes, en cuyo caso puede generarse una línea de absorción homogéneamente ensanchada. Sin embargo, si uno está interesado en el desarrollo de la radiación de un sistema de dos niveles impulsado por un campo monocromático y sumergido en un medio condensado o reservorio térmico, la suposición ya no es válida. En el presente trabajo de revisión, y bajo estas consideraciones, y tras el cálculo de la coherencia y susceptibilidad del sistema, hemos analizado las respuestas ópticas para un modelo estándar de un colorante orgánico mediante espectroscopia de mezcla de cuatro ondas. Tanto en los casos de mezcla de cuatro ondas de frecuencias degeneradas como las no degeneradas, han sido ampliamente aplicadas durante varios años, y más recientemente, han sido objeto de renovado interés en el campo de la nanociencia, donde este tipo de espectroscopia se utiliza para determinar la estructura de los nanomateriales a través de sus respuestas ópticas no lineales. Estos estudios están aún en desarrollo, prometiendo abrir una intensa línea de investigación, donde los modelos que proponemos tienen una amplia aplicación. Nuestro grupo de investigación, ha empleado diferentes metodologías para resolver la OSBE y obtener expresiones que permitan la comprensión de los perfiles ópticos de absorción y dispersión con el fin de caracterizar los sistemas moleculares en estudio.

Para efectos del presente trabajo, nuestros estudios pueden ser presentados de la siguiente forma: a) métodos de perturbación como una forma de cuantificar los efectos de saturación por parte del campo intenso de bombeo, con el fin de obtener expresiones analíticas de los elementos de la matriz densidad reducida (Paz y Mendoza-García, 2012; Colmenares y col., 1995; Mastrodomenico y col., 2008; Paz y col., 2017 y Romero-Depablos y col., 2009); b) efectos de la distribución de probabilidad de la forma de línea de absorción (Mendoza-García y col., 2010); c) consideraciones sobre la estructura molecular de un modelo de molécula simple de dos estados. Distinguimos aquí, los efectos del cruce de curvas en acoplamientos vibrónicos en la aproximación diabática (Paz y col., 1994; García-Sucre y col., 1994; Squitieri y col., 1994; García-Sucre y col., 1993), y los efectos de la difusión espectral en estados electrónicos con alta densidad de estados vibro-rotacionales (García Golding, 1983; Mourou, 1975; Yajima y col., 1978). Precisamente entender la naturaleza de la estructura molecular y los procesos de relajación cuando un sistema molecular interactúa con una radiación láser, y caracterizar este último formalmente, permitió reconocer a Gerard Mourou en 2018 con el Premio Nobel en Física.

Dada la complejidad en la formulación de un modelo que considere simultáneamente estos estudios, presentaremos en este trabajo los aspectos más relevantes de cada una de estas propuestas. En todos nuestros resultados es evidente que las propiedades ópticas no lineales son muy sensibles a la presencia del solvente, al tipo de tratamiento perturbativo de los campos electromagnéticos, al mecanismo de relajación inducido, al cruce de las curvas de potencial de los estados electrónicos iniciales y, finalmente, a la alta densidad de estados en el mecanismo de difusión espectral. El efecto de cruce de curvas por Hamiltonianos residuales del tipo spin-orbita, el cual deriva en un acoplamiento intramolecular, no será considerado en este trabajo en vista de la complejidad que se presenta cuando se trata simultáneamente este efecto vibrónico con los derivados del baño térmico. Nuestros resultados permiten estudiar con mayor precisión la correlación entre los efectos de la estructura molecular con la respuesta óptica generada. Los efectos saturativos del campo electromagnético, simetrías en el coeficiente absortivo y procesos de difusión espectral intra e interbanda, son algunas de las consideraciones que permitirán una caracterización más completa de los sistemas moleculares complejos.

2. CONSIDERACIONES TEÓRICAS GENERALES

La dinámica de un sistema de dos niveles, definidas por los estados electrónicos a, b en presencia del reservorio térmico e interactuando con una señal FWM, es descrita fundamentalmente por las OSBE en la aproximación escalar, de la forma: tρ(t)=Wξ(t)ρ(t)+R donde ρ(t) y R son la matriz densidad reducida y el vector de relajación no radiativo, respectivamente, dadas por:

ρt=ρbatρabtρDt, R=00T1-1ρDeq, Wξ(t)=ξt0iHba0ξt*-iHab2iHab-2iHba-T1-1 (1)

Wξ(t) es la matriz radiativa que incluye la variable estocástica ξ(t); ξt=iξ(t)-T2-1; ρD(t)=ρaa(t)-ρbb(t); T1; y T2 con los tiempos de relajación longitudinal y transversal, respectivamente; el superíndice (eq) define el valor de equilibrio de la diferencia de poblaciones; el elemento de matriz del Hamiltoniano radiación-materia es definido como: Hba=μba.E(t) en la aproximación de dipolo-eléctrico; μ ⃗_ba es el momento dipolar molecular de transición del sistema molecular de dos estados; E(t) corresponde al campo clásico electromagnético total, considerado como: E(t)=m=13Em, donde:

Em=E0mcos(kmz-ωmt+ϕm)=E(ωm)exp(-iωmt)+c.c., m = 1, 2, 3. (2)

Hemos seleccionado esta forma de establecer la dinámica de los estados del sistema, ya que permite incluir la simultaneidad de los efectos del reservorio térmico de manera estocástica (con un desplazamiento de la frecuencia de Bohr) y la interacción radiativa por efectos del campo electromagnético. El desplazamiento seleccionado de la frecuencia de Bohr a una frecuencia dependiente del tiempo está sujeto a procesos multiplicativos y bajo consideraciones expresadas de acuerdo con el teorema de límite central.

3. METODOLOGIAS Y RESULTADOS

Con el fin de agrupar el conjunto de resultados o las distantes variantes en el cálculo de las propiedades ópticas no lineales, presentamos un esquema general el cual supone iniciar con la ecuación de Liouville y Von-Newmann para la matriz densidad, y resolver el problema de sus elementos de matriz tanto en el dominio de tiempo como en el de frecuencias, y a partir de allí seguir una ruta de cálculo relativamente equivalente.

Figura 1.  Bajo este esquema podemos incluir: formas perturbativas del campo, efectos de saturación, distintas anchos de líneas, efectos de difusión espectral y acoplamiento intramolecular en la descripción de la estructura 

3.1 Perturbación y Cálculos de Coherencias a Tercer Orden en las OSBE

En este estudio, las OSBE fueron desarrolladas hasta el tercer orden en la coherencia, obteniendo como resultado un conjunto de ecuaciones diferenciales acopladas, de la forma:

tρ~(t)=C(t)ρ~(t)+R1(t), (3)

Donde se han definido los siguientes elementos:

ρ~(t)=ρba3tρD2tρba1tρab1t;C(t)=A-tBt000-T1-12B*t2Bt00A-t0000A+t;R1(t)=00B(t)B*(t)ρDeq

Donde A±(t)=±iξ(t)T2-1 y B(t)=-iHba(t). Los superíndices 1, 2, 3 en la matriz ρ~(t) indican los órdenes totales de perturbación de los campos electromagnéticos participantes. Para la solución de la ec. (3) es necesario introducir la forma explícita del operador Hamiltoniano con la representación clásica de los campos incidentes. El nuevo sistema de ecuaciones integro-diferenciales obtenidas (solución homogénea) es:

tρ~(t)=C~(t)ρ~(t)+R2(t) (4)

Donde

C~t=0BtC2tφ00002B*tC1t2BtC1t00000000;R2t=00BtC2tφB*tC2tφ*,C2,1t=expT2,1-1t-t0  and  φ=exptotiξ(t)dt

Aplicamos la metodología propuesta por Colmenares y col. (1995) para la resolución de estas ecuaciones diferenciales estocásticas, consideramos la aproximación de la onda rotante (Paz y col, 2015) (RWA por sus siglas en inglés) y seleccionando únicamente términos con frecuencia cercana a ω_3, obtenemos las coherencias ρba(t) y ρab(t) a primer orden en la perturbación. Sustituyendo estos resultados de primer orden y la forma explícita del Hamiltoniano, obtenemos los elementos de matriz densidad para las diferencias poblacionales en el dominio del tiempo. Usando este elemento a segundo orden y resolviendo las ecuaciones homogéneas correspondientes, podemos deducir la expresión para los elementos de la matriz densidad reducida asociada a las coherencias entre estados moleculares. Con ello, calculamos el valor promedio sobre todas las realizaciones de la variable estocástica, y bajo la aproximación ergódica, podemos finalmente encontrar:

ρba(3)(t)σ=2iΩ12Ω*ρDeqf1+f2exp-iω3texpik3.r (5)

Donde las funciones dependientes de las desintonizaciones en frecuencia bombeo-prueba vienen dadas por

f1=T~2-1+iΔ1-1T1-1-iΔ-1T~2-1-iΔ3-1; f2=T2-1+iΔ2-1T1-1-iΔ-1T2-1-iΔ3-1; Δm=ωm-ω0; m=1,2; Δ=ω1-ω2; Ωj=μba.Ej(ωj)/, se define como la frecuencia de Rabi. Aquí, hemos definido un tiempo de relajación transversal efectivo de la forma T~2-1=T2-1+γ, siendo γ la intensidad del ruido. En la ec. (4) hemos empleado para la resolución del promedio un desarrollo de cumulantes (Paz y Mendoza-García, 2012), dada por:

expimx(t)=expimx(t1)+(i2m2/2!)x(t1)x(t2)+...+(inmn/n!)x(t1)...x(tn) (6)

Donde es posible definir x(t1)=x(t1), mientras x(t1)x(t2)=x(t1)x(t2). Aquí, los cumulantes de alto orden son cero, dado el proceso gaussiano considerado. La función de correlación seleccionada ξ(t1)ξ(t2)=2γδ(t1-t2) corresponde a ruido blanco (delta-correlacionada). Usando el estado estacionario y la aproximación escalar, podemos obtener una expresión cerrada para la polarización macroscópica no lineal en la condición tensorial, dada por:

P(ω3)=χSV(ω3)E(ω3)+χeffcoup(ω3)E12E2* (7)

En la ec (7) hemos considerado que la polarización no lineal macroscópica inducida está dada por P(ω3)=Nρba(3)(ω3)σμabθ, donde el corchete interno denota un promedio sobre todas las realizaciones de la variable aleatoria, mientras que el corchete externo denota un promedio sobre las orientaciones moleculares en el sistema de estudio. Aquí, N representa la concentración química de especies activas, las cuales interactúan con el campo electromagnético externo. En (7), χSV(ω3) y χeffcoup(ω3) representan la susceptibilidad eléctrica del solvente y la componente efectiva de la susceptibilidad asociada con el proceso de acoplamiento entre los campos incidentes, respectivamente. Considerando un solvente ópticamente transparente (frecuencia lejana no resonante) al campo incidente, y tomando χeffcoup(ω3) como un número complejo, obtenemos:

Reχeffcoupω3=2g1-g2Nμba4ρDeq3; Imχeffcoup(ω3)=2h1+h2Nμba4ρDeq3 (8)

g1=Δ1Δ3Δ-Δ1/T1T~2-Δ3/T1T~2-Δ/T~2T~2T~2-2+Δ12T1-2+Δ2T~2-2+Δ32; g2=Δ2Δ3Δ+Δ3/T1T2+Δ/T2T2-Δ2/T2T2T2-2+Δ22T1-2+Δ2T2-2+Δ32; h1=1/T1T~2T~2-Δ1Δ/T~2-Δ3Δ/T~2-Δ1Δ3/T1T~2-2+Δ12T1-2+Δ2T~2-2+Δ32; h2=1/T1T2T2+Δ3Δ2/T1-Δ3Δ/T2+ΔΔ2/T2T2-2+Δ22T1-2+Δ2T2-2+Δ32

Usando la ecuación de Maxwell 2E-(1/c2)2E/t2=4π/c22P/t2 con P como la polarización local inducida y considerando la dimensión z como la dirección de propagación en el trayecto óptico, tenemos las respuestas ópticas en términos del coeficiente de absorción no lineal y el índice de refracción no lineal, dados por:

α3ω3=2πω3η3ω3cImχeffcouplω3; η(3)(ω3)=1+4πReχSV+4πReχeffcoup(ω3)1/2 (9)

En ambos coeficientes hemos considerado la aproximación de variación suave para la envolvente del campo electromagnético. Para el estudio de las respuestas ópticas no lineales, seleccionamos el colorante orgánico Cloruro de Malaquita Verde, dado su amplio uso experimental, además del comportamiento a través de un esquema de dos niveles, y el conocimiento pleno de los parámetros que lo caracterizan relacionado con sus mecanismos de disipación. Empleamos para el modelaje y cálculo de las respuestas ópticas los siguientes parámetros: μab=2.81×10-18erg-1/2cm3/2 y una frecuencia de Bohr resonante de magnitud ω0=3.06×1015s-1. En las Figs. 1 y 2 hemos representado las susceptibilidades ópticas en términos de sus componentes imaginarias y reales, Imχeffcoupl(ω3) y Reχeffcoupl(ω3) las cuales definen los coeficientes absortivos y dispersivos, respectivamente. Ambos coeficientes han sido representados como función de la desintonización de la frecuencia de bombeo respecto a la frecuencia de resonancia del sistema de dos estados Δ1 y de los tiempos de relajación. El incremento de γ, implica una reducción en el tiempo de relajación efectivo, asociado este último con la coherencia inducida. Notamos la atenuación en la propiedad óptica como una consecuencia de la orientación de los momentos dipolares moleculares con el campo eléctrico.

Figura 1.  Coeficiente de absorción no lineal como función del desintonizamiento Δ_1 y de los tiempos de relajación transversal: a) γ=0 y b) γ=1.1×10^12 

Figura 2.  Índice de refracción no lineal como función del desintonizamiento de bombeo Δ_1 y de los tiempos de relajación: a) γ=0 and b) γ=1.1×10^12 

En la Figura 2 hemos considerado el agua como solvente, con índice de refracción 1,333. El perfil de dispersión varía asimétricamente alrededor de este valor, siendo coincidente únicamente en la región de resonancia. El tratamiento que hemos realizado refleja el comportamiento de los coeficientes dispersivos y absortivos no lineales en este espacio de desintonización y relajación, cuando la incidencia perturbativa es tan sólo al mínimo orden. Esto es, la señal generada de FWM requiere una oscilación de poblaciones al segundo orden. Hemos resuelto el problema en el dominio del tiempo para la matriz densidad reducida y finalmente, evaluamos las componentes de Fourier para el cálculo de polarización. Con esta forma, hemos evitado los incidentes de correlación que se presentan a dos tiempos distintos entre el campo electromagnético y la respuesta por parte de la susceptibilidad macroscópica, esto es, la integral de polarización dada por P(t)=χ(t-t´)E(t´)dt´. Este tipo de metodologías temporales son muy válidas a medida que la interacción radiación-materia es muy débil. Sin embargo, no siempre se dispone de generaciones de FWM para campos de bombeos muy poco intensos. Este tipo de situaciones experimentales son complicadas cuando se pretende el registro y detección de la señal para muy bajas intensidades del campo perturbador.

3.2 Efectos de Saturación por Parte del Haz de Bombeo

El formalismo que a continuación se presenta establece un tratamiento igualitario para los haces de débil intensidad (prueba y señal), pero define un tratamiento del haz de alta intensidad que perturba el sistema molecular, bajo un esquema perturbativo, con desarrollos en serie de muy alto orden. Esto define lo que en la literatura se conoce como efectos de saturación por parte del bombeo. No se pretende ir adicionando términos de perturbación en la serie e ir resolviendo el problema. Se ha diseñado una metodología en la cual no se establece a priori la condición para la amplitud del bombeo, sino que se restringe el orden total, dependiendo del proceso fotónico, sujeto a la no violación de la aproximación de onda rotante. Es esta última la que condiciona el término a considerar, lo cual incluye esquemas de desarrollos infinitos para las amplitudes en el haz de bombeo. Seguimos manteniendo la presencia del reservorio térmico a través del corrimiento de la frecuencia de resonancia determinista empleando el tratamiento de Liouville en la matriz densidad reducida. Los tratamientos de saturación inducidos por “Efectos de altas potencia de bombeo” High-pump-power effects” serán considerados en el dominio de Fourier. A diferencia del caso anterior y en vista de su complejidad de cálculo, la ecuación dinámica tρ(t)=Wξ(t)ρ(t)+R se resuelve empleando transformaciones de Fourier para el cálculo de las componentes diagonales y no diagonales de la matriz densidad en el dominio de frecuencias, lo cual resulta:

ρba(ω3,ξ)=iΩ3ϒ1-2iΩ31L1*+1L3Ω12-2iΩ12Ω2*1L2*+1L1ρDdcL3ϒ1 (10)

ρba(ω2,ξ)=iΩ2ϒ1-2iΩ21L1*+1L2Ω12-2iΩ12Ω3*1L3*+1L1ρDdcL2ϒ1 (11)

ρba(ω1,ξ)=iΩ1ρDdcL1 (12)

Para la componente de frecuencia cero de la diferencia poblacional ρD(t) se tiene la expresión de la forma ρDdc=ρDeqL12T22/L12T22+4S; Lj=Lj(ξ) dada por L2n+1=i(ξ-ω2n+1)+T2-1; S=Ω12T1T2 es el parámetro de saturación y la función ϒ_1 está dada por ϒ1=Γ1(Δ)+2Ω121/L2*+1/L3, donde Γ1(Δ)=T1-1-iΔ. Es importante destacar que la presencia del solvente se la describe en términos de la distribución Lorentziana escogida. Es allí donde reside el efecto colisional (Mastrodomenico y col. 2012) y el desplazamiento de la frecuencia de resonancia se describe en términos de una variable aleatoria σ(t).

Escogemos ξ(t)=0 conforme el teorema de límite central. Tomando el promedio de la coherencia sobre la distribución ξ(t) en la aproximación escalar, obtenemos ρba(ωk,ξσ, para valores de k=1,2,3. Estas cantidades promediadas son subsecuentemente usadas en el cálculo de las susceptibilidades ópticas inducidas.

Para llevar a cabo los promedios mencionados anteriormente, Van Kampen (1981) ha propuesto un método en el que resuelve formalmente la ecuación diferencial estocástica asumiendo que es determinista, y luego toma un promedio sobre las realizaciones de la variable estocástica. Un enfoque diferente consiste en tomar el mismo promedio antes de resolver las ecuaciones ópticas de Bloch (OBE por sus siglas en inglés). En este último caso, el conjunto de ecuaciones diferenciales obtenidas puede describirse como un proceso Ornstein-Uhlenbeck (OUP). En este trabajo, resolvemos las OBE como si fueran deterministas y luego, reconociendo el hecho de que depende de las realizaciones de ξ, tomamos el promedio de su conjunto sobre la distribución de estados con frecuencia molecular entre ξ y ξ+dξ. Así, resolviendo las ecuaciones estocásticas y tomando un promedio de conjunto en lugar de un promedio de tiempo, hemos pasado por alto el problema de resolver los promedios que involucran el término de ruido multiplicativo. Con esto, la componente de Fourier de la polarización no lineal compleja local será de la forma:

P(ω3)=Nρba(ω3,ξ)σμabθ (13)

En la aproximación tensorial, obtenemos:

P(ω1)=χSV(ω1)E(ω1)+χ~eff(1,incoh)(ω1)E(ω1) P(ω2)=χSV(ω2)E(ω2)+χ~eff(1,incoh)(ω2)+χ~eff(3,coh)(ω2)E(ω2)+χ~eff(3,coup)(ω3)E(ω1)E(ω1)E*(-ω3) P(ω3)=χSV(ω3)E(ω3)+χ~eff(1,incoh)(ω3)+χ~eff(3,coh)(ω3)E(ω3)+χ~eff(3,coup)(ω3)E(ω1)E(ω1)E*(-ω2) (14)

χ~eff(3,coup)(ωj) es la susceptibilidad escalar efectiva compleja a la frecuencia ω_j, la cual refiere a los procesos de acoplamiento señal-bombeo (j=2) y prueba-bombeo (j = 3) , respectivamente; χ~eff(3,coh)(ωj)(j=2,3) y χ~eff(1,incoh)(ωk)(k=1,2,3) representan las componentes de las susceptibilidades asociadas a las contribuciones coherentes e incoherentes para la absorción y dispersión de los haces débiles, respectivamente. El superíndice en χeff representa el mínimo orden requerido para esta contribución no despreciable del proceso fotónico que tiene lugar. Es importante señalar la asimetría en las ecs. (14) para la prueba y señal en relación a la del haz de bombeo, misma que se debe al tratamiento asimétrico dado en teoría de perturbaciones de primer orden a los haces débiles y un tratamiento a todos los órdenes en el haz de bombeo. El factor de acoplamiento entre los haces de prueba y señal que podría generar fotones a la frecuencia de bombeo es muy débil y por ello, no aparece el proceso de acoplamiento en el lado derecho de la ec. (14.a). En lo sucesivo y por simplicidad, consideramos la condición de degenerancia ω1ω2, como condición relativa a los efectos de máximos de pulsación de población (Marquez y col. 1993). Considerando Yξ=4ST22Lξ2 en la expresión para ρDdc como un parámetro de desarrollo en la serie (1+Yξ)-1, expresamos las susceptibilidades ópticas de la forma:

χ~eff(1,incoh)(ωj,z)=iμba2NρDeqn=0M(-1)n4ST22nΦξ(1,n)σ(j=1,2,3) χ~eff(3,coh)(ωk,z)=-4iμba4NT13T2ρDeqn=0M(-1)n(n+1)4ST22nΦξ3,nσEω12 (k=2,3) χ~eff(3,coup)(ωk,z)=Eω1-2χ~eff(3,coh)(ωk,z) (k=2,3) (15)

Donde M es un orden requerido para garantizar la convergencia de esas expresiones. Hemos definido además Φξ(1,n)=Lξ-1Lξ-2n y Φξ(3,n)=Lξ-1Lξ-(2n+2). Aquí, las funciones Φξ(1,n) y Φξ(3,n) contienen toda la estocasticidad por la presencia del reservorio térmico. Estos desarrollos de Fourier son válidos para las funciones estocásticas estacionarias, es decir, valores medios definidos y función de correlación estacionaria. Los valores promedio del conjunto en las ecuaciones (15) se pueden calcular a partir de la forma Wξ=-WξP(ξ)dξ, donde P(ξ) es la densidad de probabilidad de los eventos con frecuencia ξ cercana a la establecida por Bohr ω0. En otros trabajos relacionados con esta temática, hemos considerado el uso de funciones del tipo Voigt como distribución de probabilidades, de mucha importancia en problemas espectroscópicos. En esos casos empleamos aproximantes asintóticos de Padé, definidos a través de desarrollos cuasi-fraccionales (Puerta y Martín, 1981; Martin y Puerta, 1981; Romero-Depablos y col. 2009; Paz y col. 2016). Sin embargo, en este trabajo y por simplicidad, consideramos la función de probabilidad del tipo P(ξ)=(2πγτ)-1/2exp-(ξ-ω0)2/2γτ, donde ξ es ahora considerada a través de un OUP con intensidad γ y una función de correlación exponencial con velocidad de decaimiento τ. El término n de la función compleja Φξ puede ser calculado si expresamos de la forma Φξ(1,n)=ReΦξ(1,n)-iImΦξ(1.n), donde ahora sus partes reales e imaginarias vienen dadas por:

ReΦξ(1,n)σ=T22nT2-(T2-1)ReΦξ(1,n-1)σ ImΦξ(1,n)σ=T22nΔ1ReΦξ(1,n-1)σ (16)

Donde el término para n = 0 viene dado por:

Φξ(1,0)σ=1γτReΘ0-iImΘ0 siendo Θ0=1-erfT2-1-iΔ12γτexpT2-1-iΔ122γτ (17)

erf(u) es la función error evaluada en u, y en la que hemos considerado que Φξ(3,n)σ=Φξ(1,n+1)σ. Una vez calculados estos términos asociados a las susceptibilidades complejas efectivas tanto incoherentes como coherentes, se obtienen para los coeficientes de absorción no lineal y el índice de refracción no lineal las siguientes expresiones:

αkωk,z=2πωkImλn=0MΛ~nSΦξ1,nσcηkωk,z; ηk(ωk,z)=ηo2+4πReλn=0MΛ~n(S)Φξ(1,n)σ1/2 (18)

Λ~n(S)=(n+1)(-1)n4ST22, λ=iμba2NρDeq/. El formalismo aquí presentado tiene aplicabilidad sólo cuando la expresión 4S/T22Lξ2 es menor que uno. Los resultados obtenidos en este trabajo son útiles para situaciones en las que se emplean haces láser de intensidad moderada o alta para la determinación de los parámetros de relajación asociados al sistema activo sumergido en un reservorio térmico. La intensidad moderada o alta queda establecida de forma experimental cuando los procesos fotónicos que comienzan a surgir tienen dependencias con términos de alto orden con la amplitud del haz de bombeo (normalmente superior al cuadrático).

3.3 Simetrías en las Respuestas Ópticas de un Sistema de Dos Niveles

Hemos visto la importancia que tiene el análisis de los parámetros de relajación en las propiedades ópticas. En este contexto óptico se propone una metodología para la estimación de los tiempos de relajación, la cual se basa en algunas propiedades de simetría de las respuestas ópticas no lineales. Estas simetrías son evidentes en el análisis de estas propiedades en función de la desintonización de los haces incidentes definidos en las superficies y las correspondientes curvas de nivel en el espacio de frecuencias. La simetría a la que nos referimos, corresponde al comportamiento de la intensidad de la señal FWM o a las propiedades ópticas no lineales del índice de absorción o refracción, en relación con las frecuencias ópticas asociadas a los haces de bombeo y de prueba. Nuestro modelo se basa en el cálculo de los tiempos de relajación utilizando los ejes mayor y menor de la elipse que representa el coeficiente de absorción después de las transformaciones unitarias, en función de la desintonización de frecuencia de los haces de bombeo y prueba. Consideramos la presencia del baño térmico a través del desplazamiento de la frecuencia del Bohr. De la ec. (4) obtenemos la componente de Fourier de la coherencia en tercer orden en la teoría de la perturbación, expresado como:

ρba(3)(ω3)ξ=2iΩ12Ω2*ρDeqexpik3rΨΔ1,Δ2 (19)

Siendo Ψ(Δ1,Δ2)=T2-2+T̑2-2-2iT2-1(Δ1-Δ2)-T̑2-1(3Δ1-Δ2)-Δ3(Δ1-Δ2)Q~1Q~3Q2*Q3Γ1(Δ). Q~k=T̑2-1-iΔk, Qk=T2-1-iΔk (k=1,2,3) y donde ahora el tiempo de relajación transversal efectivo queda definido como T̑2-1=T2-1+2γ. En la aproximación tensorial, la expresión para la polarización viene dada por la ec. (5). Usando la ecuación de Maxwell para definir la dinámica del campo electromagnético en la aproximación de envolvente suave, tenemos:

dE(ω3)dz=-α(3)(ω3)E(ω3) (20)

En donde el coeficiente absortivo α(3)(ω3) en la señal de FWM a la frecuencia ω3 es:

α(3)(ω3)=2πω3η(3)(ω3)cImχ(3)(ω3) (21)

En la Fig. 3 hemos representado el coeficiente de absorción (a) y el índice de refracción (b) no lineal, como función de las desintonizaciones de bombeo y prueba en referencia a la frecuencia de Bohr. Este tipo de representaciones son denominadas propiedades ópticas en el espacio de frecuencias. Para efectos de lo que pretendemos en esta sección, haremos uso de tales representaciones para elucidar de allí algunas propiedades de simetría inherentes a su comportamiento que permita evaluar el tiempo de relajación transversal, que como vimos anteriormente, juega un rol fundamental en la física asociada a estas respuestas ópticas. Con el fin de estudiar tales propiedades, representamos las mismas en función de las desintonizaciones en frecuencia de la forma α(3)(Δ1,Δ2) vs. (Δ1,Δ2) y (η(3)(Δ1,Δ2) vs. (Δ1,Δ2). Observamos en la Fig. 3(a) un único máximo en las coordenadas (0,0) en el plano (Δ1,Δ2).

Figura 3.  (a) Coeficiente de absorción no lineal y (b) el índice de refracción no lineal, en el espacio de las frecuencias (Δ_1,Δ_2). 

En la Fig. 3(b) para el índice de refracción se puede observar que la intensidad global varía poco a lo largo del plano. Realizando las correspondientes sustituciones, se pueden obtener expresiones tanto para el coeficiente de absorción como para el índice de refracción no lineal a la frecuencia indicada, en función de la desintonización de frecuencia, como se indica a continuación:

α(3)(Δ1,Δ2)=C1η0f00+k+j=2fkjΔ1kΔ2j (22)

En la ec. (22) hemos considerado que η(3)(Δ1,Δ2)=η01+Q2/η021/2η0, donde η01.333 correspondiente al solvente puro. Podemos observar que la superficie de dispersión presenta una inflexión precisamente en la coordenada (0,0) en el plano (Δ1,Δ2). Aquí hemos definido las funciones f00=ϕ23/a+ϕ~13/b f02=-3T2ϕ23-1T̑2ϕ~13; y los coeficientes elípticos de la forma a=T1T22 y b=T1T̑22. De igual forma, definimos las funciones ϕ23=L2LΔL3-1 y ϕ~13=L̑1LΔL̑3-1. En este caso L̑k=T̑2-2+Δk2; Lk=T2-2+Δk2 para k=1, 2, 3 LΔ=T1-2+Δ2 y el coeficiente C1=4πω3Nμba4ρDeq3c.

En este apartado se propone la determinación de los tiempos de relajación transversal y longitudinal pero en una forma distinta a lo habitual. Para colorantes orgánicos los tiempos de relajación longitudinal y transversal suelen hallarse definiendo el perfil de absorción, pero considerando los anchos medios a la altura media. Típicamente es lo que en óptica o espectroscopia se corresponde con el semi-ancho a la semi-altura (HWHM). Es importante destacar que nuestra propuesta se basa en el comportamiento de las propiedades ópticas no lineales. Por ejemplo, en el caso de los colorantes orgánicos, las relaciones siguientes 1Δj2Tk2 son válidas para los casos j, k = 1,2. Con esta aproximación ϕ23a2 y la ec. (22) se transforma en α~(Δ1,Δ2)=k+j=2fkjΔ1kΔ2j (Fig. 4). Este nuevo coeficiente de absorción no lineal puede transformarse en un espacio rotado de la forma α~(Δ~1,Δ~2)=λ1Δ~12+λ2Δ~22, donde aparecen los coeficientes λj=λj(T1,T2). En un espacio rotado ε~(Δ~1,Δ~2) se puede observar que las curvas de nivel generadas se pueden describir a través de una figura elíptica con eje mayor c2=α~/λ1 y eje menor d2=α~/λ2 (Fig. 5).

Como se puede ver, de los semiejes es posible extraer información de los tiempos de relajación. Este ejemplo de cálculo demuestra cómo al establecer las propiedades de simetría asociadas a las propiedades ópticas no lineales, podemos determinar información útil de los procesos de relajación (Paz y col. 2017). Ahora bien, dado que ambos valores propios dependen exclusivamente de los tiempos de relajación T1 y T2, a través de las expresiones aproximadas λ1T1T2-T22/2 y λ219T1T2-3T22/2, es posible determinar ambos tiempos (en unidades de s-1) directamente a partir de la curva de nivel, considerando los dos semiejes c y d, a través de un sistema de dos ecuaciones para λ1 y λ2.

Figura 4.  Representación (Curva 2D) del coeficiente de absorción en el espacio de frecuencias 

Figura 5.  Coeficiente de absorción en el espacio transformado 

Este método introduce una importante posibilidad de calcular los tiempos de relajación fundamentales que no requieren una identificación previa del sistema molecular en estudio, en las condiciones del modelo general aquí presentado (Paz y col. 2016, 2017).

3.4 Difusión Espectral y Propiedades Ópticas No Lineales

Los efectos de la difusión espectral sobre la susceptibilidad no lineal de tercer orden para sistemas cuánticos de dos y tres niveles, inmerso en un reservorio térmico, se evalúan con una señal de mezcla de cuatro ondas a través del formalismo de la matriz de densidad. Las ecuaciones ópticas de Bloch se resuelven para un sistema cuántico de dos niveles que permite la relajación en un estado cuántico del reservorio y para un sistema cuántico de tres niveles sin considerar la relajación en el reservorio; en ambos casos, se incluye la difusión espectral. Es decir, se comparan las modificaciones en la respuesta óptica no lineal para el caso en el que el tercer estado |C⟩ del sistema pertenece al reservorio y el caso en el que pertenece al sistema molecular activo. Además, en los modelos cuánticos se utiliza una distribución de la frecuencia natural, que sigue una función lorentziana. De esta manera, se puede investigar el efecto potencial de la relajación cruzada sobre la absorción óptica no lineal y las respuestas de dispersión. Consideramos un modelo de un sistema cuántico de dos niveles no homogéneamente ensanchado (Fig. 6).

Figura 6.  Diagrama de energía para un sistema de dos niveles, el cual incluye la relajación de cruce en los estados bases y excitados. El sistema relaja sobre un estado del reservorio 

Para considerar la difusión espectral en este tipo de sistemas y evaluar su incidencia en las respuestas ópticas, resolvemos las ecuaciones de Bloch ópticas convencionales. A diferencia de los casos anteriores, nuestro sistema de estudio contiene un estado del reservorio en el mecanismo de relajación. No explicitamos, como en los casos anteriores, un corrimiento en la frecuencia de Bohr por presencia del baño térmico. En otras palabras, hacemos uso de las ecuaciones de Bloch ópticas convencionales pero a diferencia de los casos anteriores, explicitamos los mecanismos de relajación hacia y dentro de la banda de estados. Necesitamos obtener la coherencia al tercer orden. A través del desarrollo perturbativo, es necesario consolidar las coherencias al primer orden, y con ellas los términos de diferencias poblacionales al segundo orden (Moncada y col. 2018). En este tratamiento, resolvemos los elementos de matriz densidad tanto diagonales como no-diagonales en el espacio de frecuencias. Las coherencias al tercer orden están dadas por:

ρba(3)(ω3)=-iμba2μbaE12E2*ρD(0)3D3(Γ3+γa+γb-iΔ)1D1+1D2*f2+g(ω0)g2 (23)

Γ3 es la velocidad de relajación o inverso al tiempo de relajación de difusión espectral,

f2=2+γbΓ3+γc-iΔ; g2=-2πΓ3Δ2+iΔΓ3+2γc+γb-γc(γc+Γ3)-γb2Γ3+γa+γb+γcγc-iΔγa+γc-iΔΓ3+γc-iΔ

En la coherencia dada por la ec. (23) observamos la función de distribución g(ω0) seleccionada en este caso como una lorentziana g(ω0)=Γ4(ω0-ωc)2+Γ42-1/π de amplio uso en espectroscopia. Luego de resolver la integral:

P(ω3)=N0μabρba(3)(ω3)+c.c.dω0 (24)

Obtenemos la susceptibilidad a tercer orden a la frecuencia ω_3, dada por:

χ2E+R(3)(ω3)=-iA(If2+I1g2)Γ3+γa+γb-iΔ (25)

Distinguiendo que esta respuesta se corresponde para un sistema de dos estados en presencia del reservorio térmico (2E+R). En la ec. (25) tenemos que:

A=Nμba4ρD(0)/3;

J=ilnΓ22+ω22Γ22+ω12+tan-1ω2Γ2+tan-1ω1Γ2π

I=0D3-1D1-1+(D2*)-1dω0πΓ2-iΔ

I1=0g(ω0)D3dω0ig(ω3)12lnω32Γ42+ωc2-ωc-ω3-iΓ2Γ4π2tan-1ωcΓ4.

Se trata de una banda de forma lorentziana centrada en ω_c y con un semi-ancho a la semi-altura dado por el parámetro Γ_4. Es importante además considerar que las cantidades encontradas anteriormente están sujetas al conjunto de condiciones Γ4,ω1,ω2Γ2,γaΔ. Estas condiciones surgen, porque en el proceso de FWM es conveniente fijar frecuencias de bombeo y de prueba cercanas entre sí y con respecto a la frecuencia de resonancia del par de estados moleculares. En las Figs. (7) consideramos las partes reales e imaginarias de la susceptibilidad óptica como función de la magnitud adimensional ΔT2.

Figura 7.  Parte real (a) y parte imaginaria (b) de la susceptibilidad óptica como función del parámetro ΔT2 para distintas fortalezas del proceso de difusión espectral en las condiciones γa=γb=10γc=Γ2  

Observamos cómo la propiedad dispersiva o absortiva disminuye en su intensidad a medida que aumenta la intensidad de la difusión hacia el estado del reservorio. Notamos que el efecto de contraste dispersivo se hace cada vez más plano con la mayor disipación desde el estado molecular superior. Así mismo, notamos que la absorción neta del estado inferior al superior se hace menor, dado que la energía que se absorbe se disipa inmediatamente hacia el estado del reservorio. En la coherencia a tercer orden a la frecuencia de la señal de FWM, podemos observar que en el caso que no exista difusión espectral entre la banda de los estados moleculares, ni exista el canal de relajación hacia el estado del reservorio térmico, la coherencia se reduce a

ρba(3)(ω3)=-2iμba2μbaE12E2*ρD(0)3D3Γ1(Δ)1D1+1D2* (26)

Donde Γ1(Δ)=T1-1-iΔ, es decir, el único mecanismo de disipación que tiene lugar es el que corresponde entre las poblaciones moleculares mediado por el tiempo de relajación longitudinal. Este resultado fue empleado por Franco y col. (1990) con algunos conceptos propios de teoría de grupos, para el estudio de las propiedades de simetría de la señal de FWM en el espacio de Fourier.

Una representación distinta y más realista corresponde a la representación de los procesos multifotónicos que tienen lugar en un esquema de tres niveles cuánticos. Es decir, a diferencia del caso anterior, consideramos tres estados cuánticos moleculares sin la presencia del estado del reservorio. Conceptualmente, a diferencia del caso anterior, en este caso debemos resolver un sistema matricial de 3x3 donde todos los estados son de la molécula y la presencia del reservorio térmico sólo se explicita a través de los tiempos de relajación longitudinal y transversal. Para ello tenemos el siguiente diagrama (Fig. 8):

Figura 8.  Diagrama de energía para un sistema cuántico de tres estados 

En este sistema consideramos como frecuencia de resonancia aquella que se determina desde los puntos medios de cada una de las bandas de los estados moleculares. Ahora, el estado |C⟩ es un canal de relajación de la molécula (estado sin ensanchamiento). Para efectos de cálculo, es necesario que este tercer estado se encuentre lejos de la región de resonancia, referida a las frecuencias de los haces de bombeo y de prueba, pero colocado exactamente a la mitad de la diferencia neta de energía entre los estados moleculares fundamental y excitado. A diferencia del caso anterior, donde es necesario resolver una matriz densidad 2x2 y explicitar la difusión espectral intrabanda y definir la cinética de relajación hacia el estado del reservorio, ahora tenemos que resolver el sistema de la matriz densidad empleando el siguiente conjunto de ecuaciones diferenciales en el dominio del tiempo:

tρbb=-iHbaρab-ρbaHab+Hcaρac-ρcaHac-γaρbb-γbρbb-ρbbT3+g(ωab)T30ρbb(ωab',t)dωab'tρjk=-iHjk(ρjj-ρkk)+Hklρlj-ρklHlj-Γjk+iωjkρjk (27)

Para j,k,l={a,b,c} y j≠k≠l

tρmn=tρnm* con n,m=a,b,c y nm

Las disminuciones de la tasa de población entre los niveles del sistema se etiquetan como se muestra en la Fig. 8. El inverso de los tiempos de relajación transversal del sistema se denomina Γij. En este modelo, la frecuencia natural entre los estados moleculares ensanchados es ωab=ω0 y el inverso del tiempo de relajación transversal es Γab=Γ2. La condición impuesta de estar distante de la región de resonancia para la frecuencia de los haces incidentes, permite no considerar explícitamente en el cálculo los términos anti-resonantes dados por sumas de frecuencia de la forma ωac+ωj y ωbc+ωj para valores de j=1,2,3. En la literatura se corresponde con el uso de la aproximación de la onda rotante. Cuando estos elementos de coherencia de la matriz de densidad desaparecen, las ecuaciones diferenciales se resuelven de la misma manera que las del sistema cuántico de dos niveles. El mismo método de aproximación sucesivo utilizado anteriormente, produce las coherencias a primer orden de aproximación y diferencias de poblaciones en el segundo orden, ambas referidas a la amplitud de los campos incidentes de bombeo y prueba. Con estas magnitudes es posible evaluar la coherencia al tercer orden responsable del proceso de FWM, esta es

ρba(3)(ω3)=-iμba2μbaE12E2*ρD(0)3D3Γ3+γa+γb-iΔ1D1+1D2*f3+g(ω0)g3 (28)

Donde:

f3=2+γbΓ3+γc-iΔ+iΓ3γcΔ (29)

g3=-2πΓ3Δ2+iΔΓ3+γc+γb-γcΓ3-γb2Γ3+γa+γbγa+γb-iΔγc-iΔΓ3-iΔ (30)

En este caso observamos que la coherencia presenta una similitud de carácter analítico con relación a la calculada para el sistema de dos estados en presencia del baño térmico a través del estado |C⟩. La diferencia implícita se corresponde con las funciones f y g en ambos casos. En la función f_3 encontramos un término adicional que da cuenta de la cinética de relajación hacia el tercer estado molecular |C⟩, no existente en el caso anterior. De igual forma, es posible establecer correlaciones entre ambos tratamientos a través de la función g de la forma:

g3(Γ3)Γ3g2(Γ~3)Γ~3-γc donde: Γ~3=Γ3+γc (31)

Esta correlación de funciones g, que a su vez puede representar una correlación entre los distintos modelos, es válida en condiciones donde γ_b>γ_c. Resolviendo la integral para la ecuación de polarización y expresando el resultado bajo la forma tensorial, es posible calcular la susceptibilidad óptica de la forma:

χ3E(3)(ω3)=-iAIf3+I1g3Γ3+γa+γb-iΔ (32)

Distinguiendo que esta respuesta está referida a un sistema de tres estados (3E). Notamos de igual forma que la forma funcional de esta susceptibilidad es la misma al comparar con el sistema descrito anteriormente. Las propiedades ópticas dispersivas y absortivas calculadas a través de las partes reales e imaginarias de esta susceptibilidad, se representan en las Figs. 9(a) y (b), respectivamente:

Figura 9.  Parte real (a) y parte imaginaria (b) de la susceptibilidad macroscópica en un sistema de tres estados, dos de los cuales se ensanchan para producir el efecto de difusión espectral, bajo las condiciones γa=γb=10γc=Γ2

La respuesta dispersiva en un sistema de tres niveles muestra dos variaciones repentinas en la región ΔT3, la cual podemos apreciar en el recuadro señalado en la parte inferior. Por otra parte, en la propiedad absortiva se muestra un efecto en la misma región de resonancia ΔT20, donde la absorción para elevados valores de difusión espectral pareciera ser de carácter negativo y convertirse a uno parecido a un proceso de amplificación parametrizada por las condiciones de cálculo. Este efecto de manera similar, se presenta en determinados sistemas moleculares expuestos a perturbaciones con campos eléctricos de alta intensidad, lo que comúnmente se define por efecto saturativo como un quemado de hueco “Hole Burning” (Mourou, 1975). Aquí, el efecto no se produce por saturación en la región de resonancia, sino por intensidad del proceso de difusión en sistema de tres estados. Pareciera que la no disipación hacia el solvente químico (reservorio térmico) genera efectos de amplificación paramétrica. Estos efectos en la propiedad dispersiva y en la absortiva son consecuencia directa del término iΓ3γc/Δ, el cual en la región de resonancia para la dispersión genera efectos similares a un “swiche” óptico, mientras que en la región de resonancia para la absorción genera una amplificación paramétrica. Esta última fue estudiada extensivamente por Boyd y col. 1990.

De los esquemas de cálculo presentados, observamos que las diferencias entre los dos modelos se deben al término iΓ3γc/Δ. Este término no es cero si y sólo si hay la consideración simultánea de un tercer estado de la molécula y difusión espectral. Aparece en el sistema cuántico de tres niveles debido al término en la función f que proviene de la teoría de relajación cuántica. El efecto no lineal anómalo no aparece en el sistema de dos niveles, porque en este modelo la población del estado del reservorio térmico fue considerada a través de la ρaa+ρbb+ρcc=g(ω0). Si en el sistema cuántico de tres niveles no consideramos explícitamente la relajación a un tercer estado o no tomamos en cuenta la difusión espectral, el término iΓ3γc/Δ0, y los efectos mostrados en los insertos en las Figs. (9) desaparecen. Los efectos anómalos no lineales aparecen únicamente en el caso donde la relación de cinéticas entre los estados cumple la condición 0<Γ3<10Γ2 en el límite cuando en Γ3Γ2. Ambos modelos predicen la misma respuesta debido a que en esta condición límite se cumple f2=f3=2 y g2=g3. Bajo tales condiciones ambos modelos convergen y la respuesta óptica de susceptibilidad está dada por la siguiente expresión:

χ(3)(ω3)=-2iNμba4ρD(0)3γa+γb-iΔΓ2-iΔ1+γb2γc-iΔ (33)

Es decir la ec. (33) se refiere a la correspondencia entre los modelos de sistemas cuánticos de dos estados en presencia del reservorio térmico y el sistema de tres estados moleculares. Si en el sistema de dos estados en presencia del baño térmico eliminamos los efectos de difusión, se recuperaría un sistema de tres estados sin difusión espectral equivalente al presentado por Yajima y col. (1978). De la misma manera, al establecer en ambos modelos la condición γb=γc=0, es decir, al considerar un sistema cuántico de dos niveles con relajación cruzada, se reproducen nuevamente los resultados correspondientes presentados por Yajima y col. (1978), de acuerdo con la siguiente ecuación:

χYajima(3)(ω3)=-2iNμba4ρD(0)3γa+Γ3-iΔ1+πΓ3ΔI1γa-iΔ (34)

Los modelos analizados de interacción radiación-materia donde se pretende estudiar propiedades ópticas no lineales y caracterizar a través de las mismas los sistemas moleculares complejos, pueden relacionarse adecuadamente a través de la inserción de mecanismos de relajación que conecten el sistema molecular con la presencia del baño térmico. Normalmente este último se lo incluye a través de los tiempos de relajación longitudinales y transversales. Sin embargo, explicitar la presencia del baño térmico bien a través del corrimiento de la frecuencia de Bohr por efectos colisionales, o definirlo como canal de relajación como los presentados en esta sección, ofrece una mayor potencialidad en los estudios, cuando se desea correlacionar las respuestas ópticas con efectos de amplia importancia, como son los sensores ópticos en dispositivos opto-electrónicos.

Finalmente es importante destacar con la metodología presentada que, bajo ciertas condiciones, es muy útil el trabajo con los estados cuánticos a través de la matriz densidad empleando ecuaciones ópticas de Bloch, tanto en forma determinista como estocástica. De igual forma y bajos ciertas condiciones relacionadas a los aspectos cinéticos de los mecanismos de disipación, que tratan (o simular) un sistema de tres estados moleculares, como un sistema de dos estados moleculares y definir la presencia del reservorio térmico como un tercer estado sólo para la relajación, es de suma importancia y de gran utilidad. Correlacionar estas metodologías de estudios de 2E+R con 3E es de importancia dado que reduce considerablemente la forma de cálculo en términos del número de ecuaciones diferenciales a resolver para los elementos de matriz densidad diagonales y no diagonales.

4. COMENTARIOS FINALES

El presente trabajo de revisión pretende dar al lector una visión de distintas formas de caracterización de sistemas moleculares complejos a través de las propiedades ópticas no lineales. En ciertos casos donde la interacción radiación-materia es débil, el sistema puede describirse como un sistema de dos estados y emplear teoría de perturbaciones para la amplitud del campo para tratar las respuestas ópticas en señales de mezclas de cuatro ondas. Sin embargo, cuando el campo perturbador es relativamente intenso, muchos procesos multi-fotónicos suelen presentarse y es necesario estudiar las respuestas ópticas bajo un esquema saturativo. Resolver los problemas de matriz densidad tanto en el dominio del tiempo como en el espacio de frecuencia, tiene cierta significancia cuando lo que se pretende es determinar parámetros específicos del sistema, como son los tiempos de relajación, de altísimo interés en diseños opto-electrónicos. De igual forma, cuando se incluye una variable aleatoria como es el caso presentado de corrimiento de la frecuencia de Bohr a una función dependiente del tiempo, evaluar los promedios en el conjunto estadístico sobre las realizaciones de esta variable, puede presentar un obstáculo complejo. Metodologías propuestas por Van Kampen y col. (1981) –aún vigentes- dan respuesta sencilla a estas dificultades. Hemos presentado una metodología para la determinación de los tiempos de relajación tanto longitudinal como transversal empleando propiedades de simetría intrínsecas en las propiedades ópticas. Equivalencias en ciertas respuestas al invertir una de las desintonizaciones de frecuencia bombeo o prueba, o variación simultánea de ambas desintonizaciones, pueden verse reflejadas en un espacio de frecuencias rotado. Aquí podemos, a partir de las determinaciones de los ejes elípticos mayores y menores, tener información de estos parámetros temporales. Brindar una metodología gráfica sencilla para su determinación, evita esquemas de cálculos más complejos, donde los tiempos de relajación se extraen de los semi-anchos a las semi-alturas máximas en los perfiles ópticos. Finalmente, poder comparar cálculos en sistemas de dos y tres estados con aquellos donde se explicita el estado de reservorio térmico, son de gran utilidad no sólo práctica en el cálculo de las matrices densidad y resolución de ecuaciones ópticas de Bloch, sino de importancia en diseños opto-electrónicos. En lo sucesivo, deseamos implementar un estudio de sistemas moleculares complejos a través de resoluciones temporales de la matriz densidad en presencia de baños térmicos que induzcan una aleatoriedad al sistema, pero considerando bases acopladas o cruce de curvas. En estos casos, buscamos resolver el asunto de los promedios de conjunto, empleando teoremas de Novikov al reconocer la coherencia como un funcional y empleo de ruidos multiplicativos para su solución. Considerar acoplamientos intramoleculares con previa inclusión de Hamiltonianos residuales en el esquema de cálculo, puede dar razón de una mayor y mejor comprensión de los procesos de absorción de dos fotones sin necesidad de abandonar la aproximación de onda rotante.

AGRADECIMIENTOS

JLP agradece al Vice-Rectorado de Investigación y Proyección Social de la Escuela Politécnica Nacional, Proyectos PII-DFIS-03-2017 y PII-DFIS-2019-02.

REFERENCIAS

Abramavicius D., Valkunas L., Mukamel S. (2007) Transport and correlated fluctuations in the nonlinear optical response of excitons. Europhys. Lett. 80, 17005 [ Links ]

Agarwal G.S., Nayak N. (1984). Multiphoton processes in two-level atoms in two intense pump beams. J. Opt. Soc. Am. B 1, 164-166 [ Links ]

Al-Saidi I.A-D.H., Abdulkareem S.A.-D. (2016). Nonlinear optical properties and optical power limiting effect of Giemsa dye. Opt. Laser Technol. 82, 150-156 [ Links ]

Apanasevich P.A., Kilin S. Ya, Nizovtsev A.P., Onischchenko N.S. (1984). On “anomalous” free induction decay rate. Opt. Commun. 52, 279-282 [ Links ]

Boyd R.W., Mukamel S. (1984) Origin of spectral holes in pump-probe studies of homogeneously broadened lines. Phys. Rev. A 29, 1973-1983 [ Links ]

Boyd R.W., Raymer M.G., Narum P., Harter D.J. (1981) Four-wave parametric interactions in a strongly driven two-level system. Phys. Rev. A 24, 411-423 [ Links ]

Brewer R.G., De Voe R.G. (1984). Brewer and DeVoe Respond. Phys. Rev. Lett. 52, 1354-1355 [ Links ]

Colmenares P.J., Almeida R., Paz J.L., (1995). Stochastic effects of the solvent on the absorptive and dispersive processes in a two-level system up to first order in the external field. J. Phys. B: At. Mol. Opt. Phys. 28, 4377-4385 [ Links ]

García-Golding F. (1983). Spectral-diffusion effects in polarization spectroscopy. J. Opt. Soc. Am. B 73, 59-62 [ Links ]

García-Golding F, Marcano A. (1985). High-order effects in Rayleigh-type optical mixing. Phys. Rev. A 32, 1526-1530 [ Links ]

García-Sucre M. Paz J.L., Squitieri E., Mujica V. (1993). Intramolecular Coupling Effect in the Refractive Index for a Simple Three-Level Model of Molecules Diluted in Water. Int. J. Quant. Chem. 27, 699-711 [ Links ]

García-Sucre M., Squitieri E., Paz J.L., Mujica V. (1994). Absorptive and dispersive processes in a two-level molecule with intramolecular coupling and non-zero permanent dipole moment J. Phys B: At. Mol. Opt. Phys. 27, 4945-4972 [ Links ]

Marquez L., Reif I., Franco H.J. Marcano A, Paz J.L. (1993). High-pump-power effects on resonant nearly degenerate four-wave-mixing signal for homogeneously and inhomogeneously broadened two-level systems, Phys. Rev. A 47, 4185-4192 [ Links ]

Martín P., Puerta J., (1981). Generalized Lorentzian approximations for the Voigt line shape. Appl. Opt. 20, 259-263 [ Links ]

Mastrodomenico A., Gorayeb M., Paz J.L. (2008). Collisional effect of the solvent on the optical responses of a two-level system. J. Nonlinear Opt. Phys. & Matter. 17, 213-224 [ Links ]

Mastrodomenico A, Izquierdo M, Paz J.L. (2012). Solvent collisional effects in Four-wave mixing, J. of Nonlinear Opt. Phys. & Matter 21, 1250016 [ Links ]

Mendoza-García A., Paz J.L., Romero Depablos A., Castro E., Martín P. (2010). Determination of nonlinear optical properties using the Voigt function: Stochastic considerations. J. Quant. Spect & Rad. Transf. 111, 155-159 [ Links ]

Moncada F., Paz J.L., Lascano L., Costa-Vera C. (2018). Effects of spectral diffusion on the nonlinear optical properties in two- and three-state quantum systems in a four-wave mixing signal. J. Nonlinear Opt. Phys. & Matter. 27, 1850038 [ Links ]

Mourou G. (1975). Spectral hole burning in dye solutions. IEEE J. Quant. Electron. 11, 1-8 [ Links ]

Paz J.L., Alvarado Y., Lascano L, Costa-Vera C. (2018). Three levels of propagation of the Four-wave mixing signal. Results in Phys. 11, 414-421 [ Links ]

Paz J.L., Franco H.J., Reif I., Marcano A., García-Golding F. (1988). Pump-power dependence due to parametric amplification of the Rayleigh-type optical-mixing signal. Phys. Rev. A 37, 3381-3385 [ Links ]

Paz J.L., García-Sucre, Squitieri E., Mujica V. (1994). The effect of intramolecular coupling on the optical susceptibilities of a two-level molecule. Chem. Phys. Lett. 217, 333-341 [ Links ]

Paz J.L., León-Torres J.R., Lascano L., Alvarado Y., Costa-Vera C., (2017). Two-Level Dipolar System in a Heat Bath: High-Pump Power Effects in the Nonlinear Optical Responses. Braz. J. of Phys. 47, 310-320. [ Links ]

Paz J.L., León-Torres J.R., Lascano L. Costa-Vera C., (2017). Relaxation times and symmetries in the nonlinear optical properties of a two-level system. Opt. Commun. 405, 238-243 [ Links ]

Paz J.L., Mastrodomenico A., Costa-Vera C., Cárdenas-García J., Rodríguez L. (2015). Rotating wave approximation effects on the nonlinear optical responses of complex molecular systems using a Four-wave mixing signal. J. of Mod. Opt. 62, 403-411. [ Links ]

Paz J.L., Mendoza-García A., (2012). Solvent influence on the nonlinear optical properties of molecular systems in the presence of degenerate and non-degenerate Four-wave mixing. J. Mod. Opt. 59, 71-82 [ Links ]

Paz J.L., Rodríguez L., Cárdenas-García J., Costa-Vera C. (2016). Solvent effects in the nonlinear optical properties using the Voigt function. J. Nonlinear Opt. Phys. & Matter. 25, 1650016. [ Links ]

Paz J.L., Rodríguez L, León-Torres J.R., Cárdenas-García J., Lascano L, Costa-Vera L. (2016). Symmetry studies in the nonlinear optical properties of molecular systems. Opt. Pur. Apl. 49, 205-218 [ Links ]

Puerta J., Martín P. (1981). Three and four generalized Lorentzian approximations for the Voigt line shape. Appl. Opt. 20, 3923-3928 [ Links ]

Romero-Depablos A., Paz J.L., Gorayeb M., Castro E., Martín P. (2009). Absortive and dispersive optical properties in molecular systems. Opt. Commun. 282, 1807-1814 [ Links ]

Singh B.P., Prasad P.N., Karasz F.E. (1988). Third-order nonlinear optical properties of oriented films of poly (p-phenylene vinylene) investigated by femtosecond degenerate four wave mixing. Polymer 29, 1940-1942 [ Links ]

Souma H., Yajima T., Taira Y., (1980). Ultrafast Relaxation Study by Resonant Rayleigh-Type Mixing Spectroscopy Using Picosecond Light Pulses. J. Phys. Soc. Jpn, 48, 2040-2047 [ Links ]

Squitieri E., García-Sucre M., Paz J.L., Mujica V. (1994). Refractive index in a dilute solution of molecules with intramolecular coupling up to third order in the external field. Mol. Phys. 82, 227-234 [ Links ]

Sudheesh P., Siji-Narendran N.K., Chandrasekharan K., (2013) Third-order nonlinear optical responses in derivatives of phenylhydrazone by Z-scan and optical limiting studies-influence of noble metal nanoparticles. Opt. Mater. 36, 304-309 [ Links ]

Sugisaki M., Yanagi K., Cogdell R.G., Hashimoto H. (2007) Unified explanation for linear and nonlinear optical responses in ß-carotene: A sub-20-fs degenerate four-wave mixing spectroscopic study. Phys. Rev. B 75, 155110 [ Links ]

Van Kampen N.G. (1981). Stochastic Processes in Physics and Chemistry: North-Holland N.Y. [ Links ]

Y. Wang, Nonlinear optical properties of nanometer-sized semiconductor clusters. (1991). Acc. Chem. Res. 24, 133-139 [ Links ]

Wodkiewics K., Eberly J.H. (1985). Shot noise and general jump processes in strong laser-atom interactions. Phys. Rev. A 31, 2314-2317 [ Links ]

Yajima T., Hirotoshi S., Ishida Y., (1978). Study of ultra-fast relaxation processes by resonant Rayleigh-type optical mixing. II. Experiment on dye solutions. Phys. Rev. A. 17, 324-334 [ Links ]

Yajima T., Souma H., (1978). Study of ultra-fast relaxation processes by resonant Rayleigh-type optical mixing. I. Phys. Rev. A 17, 309-323 [ Links ]

Published: 31/07/2019

Recibido: 26 de Marzo de 2019; Aprobado: 18 de Junio de 2019

* jose.pazr@epn.edu.ec

Luis Lascano Lascano. Doctorado en Física de Materiales, Universidad Autónoma de Madrid e Instituto de Cerámica y Vidrio, España. Físico, Escuela Politécnica Nacional (EPN), Quito. Profesor Principal del Departamento de Física de la EPN. Miembro del Grupo de investigación en Física de la Materia Condensada. Línea de investigación: Materiales Electrocerámicos. Trabaja en proyectos relacionados con: Materiales cerámicos piezoeléctricos, Materiales ferroeléctricos relaxores, Materiales multiferroicos, Síntesis por métodos químicos y caracterización de materiales electrocerámicos.

José Luis Paz. Licenciado en Química (1983), M.Sc en Química (1985) y Doctorado en Química (Universidad Central de Venezuela, 1988). Postdoctorado en Optica Cuántica UNAM. Profesor Jubilado Emérito (2014)-de la Universidad Simón Bolívar, Caracas, Venezuela. Individuo de Número de la Academia de Ciencias Físicas, Matemáticas y Naturales, Venezuela (2009). Investigador Senior Prometeo (2014-2017) Ecuador. Línea de investigación en Óptica No lineal y Óptica Cuántica. Docente a tiempo completo EPN. https://www.researchgate.net/profile/Jose_Paz3/research.

César Costa-Vera. Físico(1992), graduado en la Escuela Politécnica Nacional, Doctorado en Física de Iones (1999) en programa cooperativo entre la Escuela Politécnica Nacional y la Universidad de Uppsala en Suecia. Posdoctorado (Beca Georg Forster), Fundación Alexander von Humboldt, en Universidad de Giessen, Alemania. Miembro de la Academa de Ciencias del Ecuador. Líder del Grupo Mass Spectrometry y Optical Spectroscopy, EPN. Trabaja en óptica, fotónica, láseres y nanoespectroscopía. Profesor Principal del Departamento de Física de la EPN. https://www.scopus.com/authid/detail.uri?authorId=15041572100

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License