Population difference gratings created on vibrational transitions by nonoverlapping subcycle THz pulses

We study theoretically a possibility of creation and ultrafast control (erasing, spatial frequency multiplication) of population density gratings in a multi-level resonant medium having a resonance transition frequency in the THz range. These gratings are produced by subcycle THz pulses coherently interacting with a nonlinear medium, without any need for pulses to overlap, thereby utilizing an indirect pulse interaction via an induced coherent polarization grating. High values of dipole moments of the transitions in the THz range facilitate low field strength of the needed THz excitation. Our results clearly show this possibility in multi-level resonant media. Our theoretical approach is based on an approximate analytical solution of time-dependent Schrödinger equation (TDSE) using perturbation theory. Remarkably, as we show here, quasi-unipolar subcycle pulses allow more efficient excitation of higher quantum levels, leading to gratings with a stronger modulation depth. Numerical simulations, performed for THz resonances of the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$H_20$$\end{document}H20 molecule using Bloch equations for density matrix elements, are in agreement with analytical results in the perturbative regime. In the strong-field non-perturbative regime, the spatial shape of the gratings becomes non-harmonic. A possibility of THz radiation control using such gratings is discussed. The predicted phenomena open novel avenues in THz spectroscopy of molecules with unipolar and quasi-unipolar THz light bursts and allow for better control of ultra-short THz pulses.

www.nature.com/scientificreports/ The possibility of ultrafast creation and control of optical gratings by attosecond-long single-cycle and subcycle optical pulse trains was studied in [15][16][17] . However, generation of attosecond pulses with single-cycle and subcycle pulse duration requires a complicated setup. Besides, in the optical range, the field strengths needed to facilitate the effective interaction are extremely high ( 10 4 − 10 5 ESU ≈ 10 6 − 10 7 V/cm), which makes practical implementation even more complicated. In contrast, in the THz range, subcycle pulses can be generated much easier 1-4 . Next, the pump power needed to realize coherent interactions in THz range is several orders of magnitude lower, provided that high THz vibrational transitions possess huge dipole moments d 12 20,21 . The latter allows using THz pulses experimentally available to date, and field strengths, significantly lower than in the optical range. It worth to note that practically available subcycle THz pulses often contain a burst of single polarity and a big tail of opposite polarity and small amplitude [1][2][3][4][22][23][24][25][26][27][28][29][30][31][32][33][34][35][36][37][38] . Action of such quasi-unipolar pulses is nearly the same as the true unipolar ones 38 .
As we mentioned above, in previous works of the authors [15][16][17] grating dynamics was studied in the optical range using few-cycle pulses in a two-level medium. However, due to broad spectral content of few-cycle pulses taking into account multilevel character of the system is of high importance. Some remarks on the possibility of grating formation in multi-level systems in the optical range were shortly reported in a prior work 16,41 . However, a detailed analysis of grating dynamics under such circumstances was up to now still missing. In this paper, we study a possibility of population density grating formation on vibrational THz resonance transitions in three-and multi-level media, assuming high values of transition dipole moments (tens to hundreds Debyes), created by subcycle THz pulses. We consider a possibility to use short quasi-unipolar subcycle pulses for grating formation, because they allow more effective excitation and control of quantum state population 23,[36][37][38] . In previous studies, only bipolar few-cycle pulses were considered in this context 15,16 .
For the theoretical analysis of the multi-level system in question we use the standard perturbation approach for the time-dependent Schrödinger equation (TDSE) for "small" electric field amplitude. For numerical simulations, we use the Bloch equations for density matrix elements of a three-level medium. For analytical studies, we model the multi-level system as a quantum harmonic oscillator (HO) with equidistant levels assuming resonance frequency in THz range and high values of dipole moments of the resonant transitions. In spite of the fact, that our theoretical analysis is quite general, some particular examples will be mentioned.
To realize coherent light-matter interactions, the Rabi frequency � R = d 12 E 0 / ( d 12 is the transition dipole moment, E 0 is the electric field amplitude) should be larger than the inverse polarization relaxation time T 2 , � R > 1/T 2 and the pulse duration τ p < T 2 . First, most of molecules have vibration resonances in THz range [6][7][8][9][10] and low energy states can be modelled, to some level of approximation, using quantum HO 42 . In our analysis below we consider a water molecule H 2 O having resonances in the interval 0.5-1 THz 43,44 . Other examples are Rydberg atoms with large quantum numbers n >> 1-they also have large dipole moments 45 , can be arranged to possess long life times 46 and resonance transitions in THz range [32][33][34]45 . Next, it is shown 47 that the confinement potential of an electron in nanostructures and quantum dots (QD) for energy below the Fermi level is approximately parabolic, thus, the HO model can be used in this case. QD media can have extremely high values of transition dipole moments (tens-hundreds Debye) as well as T 2 at low temperatures as high as hundreds ns 48,49 ; they also can be used for THz generation 50 . Furthermore, quantum cascade lasers can generate THz radiation 51 . Semiconductor nanostructures used in such devices have very high dipole moments (tens of Debye) and can operate as three-level systems 52,53 . In some media like crystals, containing impurities of rare-earth ions, the values of T 2 at low temperatures can approach extremely high values -from several seconds to several hours 54 .
There is another advantage of creating grating on vibrations transitions in THz range: In quantum HO in the first order perturbation approach, only the transition from the ground state to the 1st excited state is possible, in contrast to, for instance, 1/r potential. As we see here, this significantly reduces the probability of excitation of higher order levels, even if we go into nonperturbative regime with HO.
Furthermore, creation of the gratings on vibrational transitions in different systems opens novel opportunities in THz ultra-fast spectroscopy. For instance, diffraction of a weak probe pulse on such gratings can be used for spectroscopic measurements of T 2 -a task which is otherwise proved to be quite difficult experimentally [12][13][14] . As we show below, such gratings can be also used for control of THz radiation, for instance ultrafast pulse reshaping using fast mirrors. This paper is organized as follows. First, in "Classical picture of grating formation in a system of classical HOs", for a deeper understanding of the idea of grating formation, we consider aclassical picture of the gratings arising in a medium consisting of classical harmonic oscillators (HO). In "Multi-level vibration system: weak THz field strength" we consider a grating formation in a multi-level system using perturbation theory for TDSE, when the THz field strength is low. Here we use a δ-pulse approximation, valid when the pulse duration is smaller than the inverse transition period of the medium. After this, pulses of finite duration are considered. In "Numerical simulations" we perform numerical simulations for three-level HO using parameters for H 2 O molecule. A possibility of THz radiation control by the gratings is considered in the "Control of long THz pulses". Finally, concluding remarks are drawn.

Classical picture of grating formation in a system of classical HOs
First, we consider the classical picture of a grating formation. Let a dielectric resonant medium consist of linear classical harmonic oscillators, non-interacting with each other and having eigen-frequency ω 0 , distributed along z-axis, see Fig. 1. The displacement of the oscillator at point z is governed by the following equation (neglecting damping of oscillations): Here q is the electron charge and m the oscillator mass. Let the medium be excited by a pair of extremely-short pulses with duration much smaller than the medium proper oscillations period T 0 = 2π/ω 0 . In this approximation, and for simplicity assuming the electric field in the form of two delta pulses E(t) = S 0 [δ(t) + δ(t − �)] , counter-propagating as shown in Fig. 1, with the delay �(z) depending on the oscillator position z. We assume low enough density of the oscillators on the string, i.e., changing of the pulse shape during the propagation is negligible. Under the action of such kick-like excitation, the medium starts to harmonically oscillate at the eigenfrequency ω 0 . It can be easily shown that the response of the medium can be represented as 55 From the last expression it is seen that after arrival of the second pulse, the amplitude of the displacement periodically depends on the pulse-to-pulse delay cos (ω 0 �/2) . This can be interpreted as periodic grating formation in the system of classical oscillators. In the next section, we consider the corresponding quantum picture.

Multi-level vibration system: weak THz field strength
Theoretical approach. In the previous studies we predicted population density gratings formation in a two-level medium [15][16][17] . However, subcycle pulses have broadband spectra and two-level approximation can be inapplicable. In this section we consider a multi-level vibration system, namely, a quantum harmonic oscillator (HO). The HO model is commonly used as a simplest model for theoretical description of molecular vibrations 42 .
Interaction of a subcycle pulse with a quantum system with arbitrary configuration of states is governed by the TDSE describing the evolution of the wave function ψ 42 : H 0 is the intrinsic Hamiltonian of the unperturbed system, and V (t) = −qxE(t) is the interaction potential with the excitating pulse in the dipole approximation. The wave functions of the eigenstates are the eigenfunctions of the intrinsic Hamiltonian Ĥ 0 and given by 42 : where H n x mω 0 -Hermite polynomials of order n. The general solution of the TDSE (3) can be written as a superposition of eigenstates ψ n with the amplitudes a n : H n x mω 0 , Figure 1. The simplest setup with four THz pulses counter-propagating in a medium to create a grating. Inset shows a 4-level HO molecular medium with the eigenfrequency ω 0 . This figure was created with Paint application for Windows and converted to eps using online convertor (https ://image .onlin e-conve rt.com/ru/ conve rt-to-eps). www.nature.com/scientificreports/ Let us assume that the system was initially in the ground state with n = 0 and the pulse amplitude is small. The amplitudes of the eigenstates can be easily calculated in the 1st order of the perturbation theory 42 : The transition probability of a HO from the ground state of the discrete spectrum to the n-th state w 0n is given by |a I n | 242 : V 0n is the matrix element of the perturbation operator, which is expressed through the matrix element of the dipole moment d n,n+1 as: For HO, the matrix d n,n+1 is given by 42 : Other matrix elements are zero, because in HO only transitions between neighbouring states are possible. That is, in the 1st order perturbation approach, only transition from the ground state to the 1st excited state has nonzero probability, i.e., This allows us to consider HO as an effective two-level medium, as least in perturbative approximation. This gives rise to high-contrast grating formation only on the first transition 0 − 1 of the vibrational spectrum.
Population density grating formation and their control using δ-pulses. Before providing detailed analysis of the grating dynamics let us consider for simplicity interaction of a HO with a δ-pulse. A validity of this approximation follows from the fact that experimentally obtained subcycle pulse shape contains a short burst of electric field (half-wave) of one polarity and a long damped tail of small amplitude of the opposite polarity [1][2][3][4][22][23][24][25][26][27][28][29][30] . In view of the small amplitude of the tail, the main contribution is given by the strong half-wave part 38 . Assuming, that the pulse duration is smaller than all resonant transition periods in the atom/molecule, we can describe our driving field as a set of delta-function-like pulses: having the time delays n and electric pulse area S E,n = E n (t)dt 56 , E n is the electric field of the n-th pulse.
In the case of a spatially extended medium, one can consider a string of atoms distributed along the z axis. If the concentration of atoms in such a string is relatively small, one may neglect nonlinear pulse shape modification due to propagation effects. A setup for a grating formation in the case with four pulses is shown in Fig. 1 (similar scheme was used earlier for a two-level medium 15,16 ). The first pulse propagates in the medium from left to right, whereas the second one travels from right to left, after the 1st one already left the medium, so that the pulses do not overlap in the medium. Under these conditions, the delay between two pulses is constantly changing with time. The 3rd pulse propagates after the 2nd one in the same direction. The 4th pulse travels in the same direction as the 1st one after 3rd pulse left the medium. To track this situation without directly solving the wave equation, we may consider every point separately. In such consideration, the difference from one point to another in z-direction is mapped to the delay between two pulses � ∼ z/c . Thus, we can consider a single-atom response and calculate the populations of the levels in dependence on the delay between two ultra-short pulses in order to show the existence of a grating as well as to estimate its properties.
First, consider the action of two pulses propagating in the opposite directions as shown in Fig. 1: where is the delay between two pulses. Using Eq. (7) and performing integration it is easy to obtain an expression for the transition probability w 01 (the only nonzero one, according to Eq. (8)): From Eq. (12), the periodic dependence of the transition probability on the delay between the pulses is seen. Since in the case of an extended medium, the delay � ∼ z/c determines the time of arrival of the second pulse to a point of the medium with coordinate z, Eq. (12) shows a possibility of a periodic grating formation in the (5) ψ(x, t) = n a n (t)ψ n (x). www.nature.com/scientificreports/ medium. Equation (12) is valid for any multi-level quantum system, when the concentration of particles is small and the change in the pulse shape during the propagation process can be neglected. However, formula Eq. (12) was obtained in the perturbation theory in the weak-field approximation. The case of an arbitrary strong driving field will be considered below. We now show a possibility of erasing and spatial frequency multiplication of the gratings. Suppose that the system interacts with three pulses, see Fig. 1: where 23 is the delay between the 2nd and the 3rd pulse. Using Eq. (7) for the transition probability it is easy to obtain From Eq. (14) it is seen that when e iω 01 � 23 = −1 , the 3rd pulse erases the grating since the probability of transition does not depend on the delay and thus on the spatial position.
Next, let the system interact with the 4th pulse (see Fig. 1): where 34 -delay between the 3rd and the 4th pulse. For the transition probability we obtain: We now show the possibility of the multiplication of the spatial frequency of the gratings analogously as it is done in a two-level system 15 . We choose the following values of the delays between pulses: � 23 = π/ω 01 + 2πk/ω 01 , that is e iω 01 � 23 = −1) and � 34 = � + 2πk/ω 01 , where k = 0, 2, 4, 6 Hence, from Eq. (16) it is easy to obtain: From Eq. (17) it is seen that the spatial frequency of the gratings is multiplied by a factor of two. Thus, we showed that by selecting the delay between the pulses it is possible to control period of the multi-level-based grating.

Population density grating formation using subcycle THz pulses of finite duration. Using
Eq. (7), we show a possibility of existence of population density gratings in HO-type systems for the pulses of finite duration. As we mentioned above, experimentally obtained half-cycle THz pulses contain a strong halfwave of high amplitude and a long weak tail of the opposite polarity [1][2][3][4] . It was shown in 38 , that the impact of this tail can be neglected, provided that its amplitude is small, and duration is larger than the width of the main half-wave. Hence, we first neglect for simplicity this tail and consider the impact on the system of two Gaussian pulses given by: where is the time delay between two pulses. In the subsequent analysis we include the tail into consideration. Substituting Eq. (18) into Eq. (7) and taking account that e ibx−ax 2 dx = π a e −b 2 4a , we obtain where we have explicitly used Eq. (8) for d 01 as well as denoted the transition frequency 0 → 1 as ω 0 = ω 01 . Expression Eq. (19) is proportional to the square of the field strength E 0 and at the first glance it seems that probability w 01 can be very high at large field amplitudes. However, Eq. (19) was obtained using perturbation approach and is valid when the electric field amplitude is small, i.e. q 2 E 2 0 πτ 2 p /m ≪ ω 0 . From Eq. (19) one can also see that the periodic dependence of the transition probability (population of the first excited state) on the delay between pulses is similar to the case of δ-function-like pulses considered earlier. The modulation depth of the grating is proportional to the square of the electric field amplitude of the pulse. Formula Eq. (19) can be interpreted as a harmonic inversion grating created by a pair of counter-propagating subcycle pulses in a spatially-extended medium. Thus, in the case of the pulses with finite duration the existence of the inversion gratings is also possible.
We apply the theory developed above to a H 2 O molecule. It has isolated resonances in the THz range, in particular, around 0.97-0.99 THz according to experimental results of 43 . Below we take ω 0 /2π ≈ 0.97 THz . Using the value for the "characteristic radius" of water R = 4 Angstrom 44 , we estimate transition dipole moment as d 01 = qR = 19.2 Debye (q is the electron charge). The two-dimensional diagram in Fig. 2 illustrates the dependence of the 1st excited state population w 01 Eq. (19) for water molecules vs the pulse duration τ p and the pulse amplitude E 0 for = 3 fs. The maximal probability value w 01 at the fixed delay is determined by the term For small pulse durations (shorter than the transition period, ω 0 τ p ≪ 1 ) this term increases with increase of τ p , reaching its maximum value at certain point. With the particular parameters mentioned above, w 01 reaches its maximum at the pulse duration τ p of around 200 fs. When the pulse duration becomes larger than    www.nature.com/scientificreports/ the transition period, ω 0 τ p ≫ 1 , the population w 01 tends to zero. It means, that an unipolar pulse with a duration smaller than the inverse transition frequency ( ω 0 τ p ≪ 1 ) acts more effectively than a long unipolar one with a larger duration and the same amplitude. This statement is valid in general case of quasi unipolar ones, see below.
It can be also seen, that the modulation depth of the gratings significantly depends on the duration of the incident pulses and their amplitudes. It is also seen from Fig. 2 that for creation of gratings, electric fields of small amplitude of ∼ kV/cm are sufficient, what can be easily reached in practice [1][2][3][4] .
Since the HO model is used to describe molecular vibrations (more precisely, low excited vibration states), and the oscillation frequencies of the molecules could lie in the THz region, the analytical result Eq. (19) clearly shows the possibility of a grating creating in molecular systems using a pair of THz subcycle pulses with experimentally available strengths, several orders of magnitude lower than in the optical range. Now, we proceed further to take into account the tails of the Gaussian pulses. We assume, instead of Eq. (18), the following more general expression for the exciting field: with the carrier-wave frequency and the carrier-envelope phase (CEP) ϕ i . For clarity, we assume below CEP for all pulses in the sequence to be fixed ϕ i = ϕ = const , i.e., that the pulses are coming from a CEP-stabilized laser pulse source. The electric pulse area for such pulses is given as: so that these pulses can be safely treated as unipolar ones if �τ p ∼ 1 or even �τ p ≪ 1 and ϕ = ±π/2.
Substituting now Eq. (20) into Eq. (7) and performing integration we find: Equation (22) again yields a periodic dependence of the transition probability w 01 (population of the first excited state) on the delay similarly to Eq. (19) above. Moreover, Eq. (19) contains some interesting peculiarities. When CEP ϕ = ±π/2 , i.e. the exciting pulses are strictly bipolar according to Eq. (21) and S E = 0 , the population gratings still arise, even though they have smaller modulation depth. The dependence of the modulation amplitude on the duration and CEP of the pump pulses is determined from Eq. (22) by the following factor: where we have introduce the factor ω 2 0 to make the parameter κ dimensionless. Figure 3 illustrates the dependence of this parameter on CEP ϕ and the pulse duration τ p . One can see that κ achieves its maximal value at τ p ≈ 0.4 ps and cos ϕ = 1 . Also, it is worth noting that in the limit → 0 Eq. (22) coincides with Eq. (12).
It is also seen from Fig. 3 that the grating modulation depth is larger for subcycle pulses ( �τ p ≪ 1 ) than for multi-cycle ones, and it is smaller for bipolar pulses (when ϕ = π/2 and the pulse area S E = 0 ). This fact is in agreement with our previous results, which showed that quasi-unipolar subcycle pulses allow more efficient nonresonant excitation of quantum systems with respect to single-cycle and multi-cycle ones 37,38 .
Let us check now the results of the previous section for the sequence of three and four pump pulses. For three pulses we have the exciting electric field in the form: cosh ω 0 �τ 2 p + cos 2ϕ (1 + cos ω 0 �).

Numerical simulations
To check our analytical theory and observe, what happens in the nonperturbative regime, we solve numerically a system of Bloch equations for density-matrix element, which models interaction of a three-level medium with THz pulses: Equations (27) Low-power THz fields. Equation (19) was obtained using perturbation approach valid for low field strength. In the 1st order perturbation approach, only probability transition w 01 is non-vanishing. This is because transition dipole moments in HO d 0,n = 0 for n larger than 1. Numerical simulations (see below) showed that populations of higher levels are 10-100 times lower than the populations of the 2nd level w 01 . This fact enables the possibility of creating the high contrast grating on the main transition 0 → 1 neglecting other transitions using weak THz field values.
To confirm this statement, we performed numerical simulations using system of density matrix equations (27)-(32) with a pair of Gaussian pulses (18). Figure 4 shows the behaviour of the populations ρ 11 (a), ρ 22 (b), and ρ 33 (c) versus the delay for the water molecule and E 0 = 15 kV/cm assuming infinite relaxation times. The influence of finite relaxation times will be studied in the next section. It can be seen from Fig. 4b that the grating obtained numerically has the same shape and roughly the same amplitude as predicted by Eq. (19). Besides, the population of the 3rd level ρ 33 is 10 times smaller than ρ 32 .

High-power THz fields.
To study the grating dynamics in the high-power field we performed a series of the numerical simulations of the Bloch equations described above at different pump field strengths and durations. Figure 5 shows the behaviour of the populations difference ρ 11 − ρ 22 (a), ρ 11 − ρ 33 (b), ρ 22 − ρ 33 (c) versus the delay between two Gaussian THz pulses and time t for the H 2 0 molecule ( ω 0 /2π = 0.97 THz, d 12 = 19.2 Debye) for the following parameters: E 0 = 150 kV/cm, τ p = 150 fs. Typical lifetimes of the vibrational levels T 1 and coherence times T 2 of vibration states are in the ps range 57,58 . We take relaxation population lifetimes T 22 = T 33 = 250 ps and coherence lifetimes T 21 = T 32 = T 31 = 5 ps 57,58 .
One can see noticeable differences to the case of low field strength described above. The grating shape is not harmonic anymore, and has complex peak structure in contrast to the weak-field case. Also in contrast to the case of the two-level system, where both analytical and numerical solutions predict harmonic shape only 15,16 . Next, due to finite values of the relaxation times, the grating modulation depth decreases dramatically with time.
To summarize the results of these sections, we see that using THz pulses and vibrational transitions introduces many differences and advantages with respect to the previously studied case of optical field excitation. It allows reducing field strength using experimentally available subcycle THz fields. One can use only the transition 0 → 1 for grating formation in the weak-field case. The shape of the grating can be controlled by increasing the field strength, in contrast to the case of a two-level medium, where grating shape was always harmonic.

Control of long THz pulses
In this section, we describe an application of the gratings described above, to control propagation of THz pulses. As it follows from the Bloch equations, the modification of polarization in the linear regime corresponding to a density grating of full contrast ( |ρ 22  www.nature.com/scientificreports/ where N is the atomic density. Taking into account that δP = ε 0 χE , we have |χ| = 2d 2 12 NT 2 /( ε 0 ) , and therefore the refractive index modification A refractive index grating is a powerful tool to make highly frequency-selective reflection of THz radiation, fully controllable on the single-cycle (picosecond) time scale. As an example, in Fig. 6 we show the dependence of reflection on a wavelength resulting from a grating containing 1000 oscillation periods, each having the length 75 µ m in a gas with atmospheric pressure consisting of active molecules, which thus creates a density grating with the contrast δn ≈ 10 −3 . The reflection is calculated using the standard multilevel approach 59 : The ABCD matrix of a single layer of size h for the normal incidence is  www.nature.com/scientificreports/ where k 0 is the wavevector in vacuum, n is the refractive index, p = √ ǫ . For our calculation we have taken /4 layers with interleaving layers of the refractive index n 1 = 1 − δn/2 and n 1 = 1 + δn/2 . Multiplying the matrices for all layer gives the resulting reflection. Although an analytic expression exists for the whole such product 59 , we find more useful direct numerical multiplication.
The refractive index grating can be used in various ways. The key advantage is a possibility of fast switching on the single-cycle level. An example shown in the inset to Fig. 6 is a fast amplitude modulation of a THz CW wave allowing to encode information into it.

Conclusions
To summarize, we studied theoretically the possibility of population density grating formation created by subcycle THz pulses in a molecular medium having resonant transitions in THz range. The pulses interact with the molecules coherently and do not overlap in the medium. The amplitude of the grating can be controlled either by tuning the pump amplitude or by an additional pulse. The grating modulation depth is proportional to the square of the THz field amplitude but is also influenced by the coherence times in the medium.
In the case of a multi-level system, using the first order of the perturbation theory in the weak-field approximation, we studied the grating dynamics created by two, three and four subcycle pulses. The possibility of a grating creation and their ultrafast control (erasing and spatial frequency multiplication) in multi-level system excited by subcycle pulses has been shown. A possible scheme of a setup for a grating formation and their control using four counter-propagating pulses on vibrations transitions in a H 2 O molecule without overlap in the medium was considered.
Our analysis revealed significant advantages of a grating creation in the THz range with respect to the optical range considered earlier. Because of high dipole moments in the THz range, one needs significantly smaller field amplitudes than in the optical range, of the order of kV/cm. This opens a new opportunity for an ultrafast THz spectroscopy with subcycle THz pulses to determine the resonance frequencies and the corresponding dephasing T 2 and relaxation T 1 times via diffraction of a weak probe pulse on the induced gratings [12][13][14] . The investigated phenomenon can be also used for creation of ultrafast beam deflectors and switches on molecular systems in the THz range.
As we have seen, somewhat counter-intuitively, the HO model in the perturbative regime provides lower excitation probabilities of the higher states than, for instance, the atomic potential 1/r, and thus the quality of the grating in the molecular case is better than in the case of an atom. This was supported also by non-perturbative simulations with the Bloch equations.
Numerical simulations performed for the parameters of H 2 O molecules and analytical results revealed that in the low-field case gratings have a harmonic shape. But in the strong-field case, when the perturbation approach is not valid anymore, the grating shape differs dramatically from the harmonic one. It has complex peak structure, since the level populations in the medium undergo multiple Rabi floppings over the pulse duration.
The obtained results also indicate a possibility of effective grating creation by quasi-unipolar subcycle THz pulses with respect to bipolar few-cycle ones. This efficient grating creation occurs in spite of the broad spectrum of the subcycle pulses covering many resonances at once, and non-resonant interaction with the molecule transitions. This result opens novel opportunities in ultrafast spectroscopy in the THz range with unipolar half-cycle pulses.