A fully-quantized nonperturbative coherent four-wave (4W) mixing theory and a three-dimensional computer program to solve the model are presented here. A detailed algorithm for the solution, program extension for extra resonances, and plotting techniques are also given.
Nonlinear perturbative susceptibilities have been used successfully to analyze experiments in nonlinear optics for two decades. By its very nature, a perturbative approach is limited to a small value of the expansion parameter, and nonlinear susceptibilities are not useful for strong saturating fields. Presently, a theory is presented for the treatment of FWM which is applicable to arbitrary strength of the input fields, both on and off resonance. The field harmonic time dependences are removed by a transformation to the proper rotating frame, and the resulting eigenvalue problem is solved numerically. Several new effects are predicted. For strong fields resonance frequencies are shifted and lines are split (similar to the ac Stark effect). When such a split line is acted on by another strong field (a phenomenon not studied so far), a new effect of stirring is observed, which results in the broadening and merging of the two separate components. A new extra resonance appears which is induced by the strong field and is present even in the absence of any collisions or other dephasing processes or of laser phase fluctuation. This new extra resonance is a new member in the family of extra resonances which were discussed in various previous publications.
The nonperturbative approach to the calculation of nonlinear optical spectra of Seidner et al. [J. Chem. Phys. 103, 3998 (1995)] is extended to describe four-wave mixing experiments. The system-field interaction is treated nonperturbatively in the semiclassical dipole approximation, enabling a calculation of third order nonlinear spectroscopic signals directly from molecular dynamics and an efficient modeling of multilevel systems exhibiting relaxation and transfer phenomena. The method, coupled with the treatment of dynamics within the Bloch model, is illustrated by calculations of the two-dimensional three-pulse photon echo spectra of a simple model system-a two-electronic-level molecule. The nonperturbative calculations reproduce well-known results obtained by perturbative methods. Technical limitations of the nonperturbative approach in dealing with a dynamic inhomogeneity are discussed, and possible solutions are suggested. An application of the approach to an excitonically coupled dimer system with emphasis on the manifestation of complex exciton dynamics in two-dimensional optical spectra is presented in paper II Pisliakov et al. [J. Chem. Phys. 124, 234505 (2006), following paper].
A nonperturbative semiclassical theory is presented for the treatment of four-wave mixing (FWM) in a four-level system. This treatment is based on a transformation to a generalized rotating frame followed by a numerical diagonalization of the resulting time-independent Hamiltonian. For strong fields, all resonant frequencies are shifted from the atomic (molecular) frequencies. Lines are split up in a way similar to the ac Stark effect for a two-level system, and for even stronger fields, the effect of stirring is observed for the first time in four-wave mixing. In addition, a new extra resonance appears which is induced by the saturating field. This field-induced resonance is the new member in the family of extra resonances (pressure induced, or stochastic fluctuation induced), and it should be observed whenever saturating fields interact with a multilevel atomic or molecular system.