En este trabajo se estudia de forma teórica las propiedades ópticas no lineales absortivas y dispersivas de sistemas moleculares complejos inmersos en un reservorio térmico que interactúan con una señal de mezcla de cuatro ondas, teniendo en consideración la estructura interna del sistema molecular mediante la inserción de Hamiltonianos residuales del tipo spin-órbita. Se hace uso de las ecuaciones de Bloch ópticas estocásticas y se calculan las componentes de Fourier de las coherencias y de poblaciones en matrices densidad reducidas. Como parámetros del sistema en los procesos de disipación, se consideran los tiempos de relajación transversal y longitudinal para los procesos estocásticos de interacción soluto-solvente, empleando procesos del tipo Ornstein-Uhlenbeck. Los efectos de acoplamiento intramolecular sobre las respuestas ópticas se estudiaron usando un modelo de ...