Wavelet Spectra of a Chirped Gaussian Pulsed-Driven Harmonic Oscillator

: We investigate the nature of spectral lines of emitted radiation due to of a non-dissipative single-mode quantized harmonic oscillator (HO) with train of n-chirped Gaussian pulses. Speciﬁcally, we analyze the transient emitted spectrum through two window wavelet functions of the radiation detector, namely the Haar wavelet and Morlet wavelet. Computational display of the exact analytical results shows how the driving pulse parameters (strength and number of pulses, chirping, repetition time) act as control knobs to shape the detected emitted spectrum as desired.


Introduction
Driven quantized harmonic oscillator (HO) is a fundamental model of radiation-radiation interaction in the field of nonlinear and quantum optics, e.g., [1,2]. The case when the driving field is a pulsed laser has been investigated regarding the fluorescence Fourier transform (FT) [3][4][5][6][7] and wavelet transform (WT) [7,8] spectra of the emitted radiation for various shapes of the laser pulse. Briefly speaking, FT analysis is a time-averaging process over the entire time interval of the detected signal, and in turn some localized-time structure may be lost due to this time-averaging [9]. On the other hand, time-localized structure of different frequency components in the detected signal can be dealt with through adjustment of the time width of the detecting "mother window function". This is known as WT analysis, in which wavelets of high frequencies are narrower and vice versa [9]. Transient wavelet spectral analysis of the emitted radiation of a single 2-level atom driven by a rectangular pulse has been given in the cases of Morlet and Mexican hat [10] and Haar [11] wavelet window functions. Also, for the case of 2-level atom, Haar wavelet spectrum has been investigated where the driving pulse is of exponentially decaying shape [12] and of sin 2 -shape [13].
In this paper, we expand our results in [5] for the model of a single non-dissipative HO driven by a field that has the shape of a train of n-Chirped Gaussian pulses and calculate the wavelet transform spectra with two different wavelet window functions, namely Haar wavelet (HW) and Morlet wavelet (MW) windows. We investigate the WT spectra for cases of initial vacuum and number states of the HO.
The paper is presented as follows: in Section 2, we present the model equations for the driven HO with the expression for the average auto-correlation function of the HO amplitude operators. Analytical expressions for the transient WT spectra are derived in Section 3 for the two cases of HW and MW window functions. In Section 4, we present the computational results for both WT spectra, with initial vacuum and number states of the HO, followed by summary in Section 5.

Model Equations and Correlation functions
The Hamiltonian model of a single-mode HO of frequency ω o driven by a laser pulse of complex envelope f (t) and of oscillating frequency ω L , within the rotating wave approximation (RWA ) and in the absence of damping processes (in units ofh = 1) has the form [14,15], The notations are: Ω 0 is the Rabi frequency associated with the driving pulse (assumed real), the annihilation and creation HO operators,â andâ † , respectively, obey the boson commutation relation [â,â † ] = 1, and h.a stands for Hermitian adjoint.
In the case, n-chirped Gaussian pulse, the envelope f (t) has the form [5], where τ = t τ o is the normalized time, τ 0 is the full width of one pulse at 1 2 -maximum, The solutions of (2) for arbitrary initial time t o are of the form, and have taken t o → −∞, since Gaussian pulses are of smooth switch-on shape, with ∆ = ω L − ω o is the atomic detuning. The quantity I(t) in (5) in this case is given by [5], where and erf(z) = 2 √ π z 0 e −z 2 dz is the error function of complex argument [16]. The average value of the auto-correlation function â † (t 1 )â(t 2 ) can be calculated from the solutions (4), As seen from the dynamical solutions for the HO amplitudes (4), (5) and the autocorrelation functions, Equation (8), that the Rabi frequency parameter Ω o does not contribute to any oscillatory behavior, i.e., no induced Rabi oscillations. This is unlike the case of driven atomic systems (of finite number of energy levels), e.g., [3,[10][11][12][13], where Rabi oscillations (particularly in the strong field limit) manifest itself in the bouncing of atoms between its energy levels. Oscillatory behavior in our driven HO model are soloy due to pulse shape and wavelet window functions parameters.

Transient Wavelet Spectra: Analytical Results
The transient spectrum of the scattered radiation is given by [17].
For the response (i.e., window) function of the photo-detector (of width Γ and frequency ω f ) is given by [8], where ψ(t − t 1 ) is the mother wavelet function. Here, we consider the following two cases of wavelet functions, ψ(t − t 1 ).

Haar Wavelet Spectrum
The mother wavelet function in this case has the bi-polar step shape [9,18] with α, β, γ are scaled time-dependent variables given by, and j o , k o are the dilation and shift parameters, respectively. The HW, although has a good time resolution, it has no frequency resolution, i.e., continuous signals are not approximated very well due to its time domain discontinuity [18,19].
Using Equation (12) into (10) and the use of Equation (9) into (8), we get the wavelet spectrum, in the case of initial number state |n o of the HO into the form, The expressions for the J 1,H (t), J 2,H (t) are given by, with,

Morlet Wavelet Spectrum
The mother wavelet function ψ M (t − t 1 ) is in this case [9,18], where j, k are the dilation and shift parameters, respectively, and ω M is the central frequency of the mother wavelet and the σ is full width at 1 2 -maximum of the Gaussian envelope. In what follows, for the Morlet spectrum, we put σ = 1 which means the parameters t, Ω o , ω o , ω f , ω M , Γ are normalized to σ. The Morlet (or Gaussian) wavelet is much suitable with signal analysis for two aspects: its oscillatory and decay nature allow for good frequency and time resolutions, respectively [18,19].
Using Equations (11), (9), (16) in Equation (10) we get the analytical expression for the Morlet wavelet spectrum with initial number state |n o of the HO into the form, The quantities for the J 1,M , J 2,M are given by, where The integral in (18c) is best treated numerically.

Computational Results
Here, we plot the analytical formulas for the wavelet spectra

Haar Spectrum
For initial vacuum state (n o = 0), the normalized spectrum S H (D H ) is governed by the quantity |J 2,H | 2 in the second term in Equation (14). At exact resonance (∆ o = 0) and for zero shift k o = 0 and no chirp (C = 0), the bi-polar nature of the HW wavelet function reflects itself in the spectrum as two symmetric peaks around D H = 0 with fading symmetric smaller peaks for one pulse (n = 1)case- Figure 1a. Increasing the pulse number n = 2, 10 tends to merge the two peaks of (n = 1) case to one central peak that is more narrowed for larger n ≥ 10 with very small two peaks around this central peak (Figure 1b). The interplay of the detuning parameter (∆ o = 0) at large shift parameter k o = 10 is shown in Figure 1b, where for one pulse (n = 1), the spectrum is asymmetric (one large peak near D H = 0 and very small flat peak near D H ∼ ∆ o + k o ). With increasing pulse number (n = 10), the spectrum has a central narrowed peak shifted to the right, D H > 0 with dense oscillations in its envelope and two flattening asymmetric peaks. For non-zero chirp parameter (C = 0), ∆ o , k o = 0, the asymmetry is more pronounced for c ≷ 0, with very dense oscillatory envelopes for the main peak shifted for D H > 0 with C > 0 (Figure 2a) or the main peak around D H ∼ 0 with C < 0 (Figure 2b).
For initial number state (n o = 1), the combination of the three parameters (∆ o , k o , C) is best seen with larger pulse number (n = 10) and large k o . In Figure 3, curve (a), for (∆ o = 3, k o = 10, C = 0) there is a single narrowed asymmetric peak to the right of D H = 0 with two asymmetric broader peaks. The non-zero chirp C > 0, causes the spectrum- Figure 3, curve (b), to have two main peaks asymmetrically displaced from D H = 0, with small spikes in between. The case of C < 0- Figure 3-curve (c), makes the spectrum less asymmetric, with central spike near D H = 0. Therefore, the oscillations in the Haar WT spectrum with initial vacuum state are due to non-zero values of atomic detuning (∆ o ) and the window shift parameter k o . These oscillations become denser with large pulse number (n ≥ 10) and non-zero chirp parameter (c = 0) Figure 2a,b and they are diffused in the case of initial number state |n o (the term promotional to n o in Equation (14), Figure 3).

Morlet Spectrum
For initial vacuum state (n o = 0), at exact resonance, the general feature of the Morlet (Gaussian) wavelet spectrum is a single symmetric peak, Figure 4a, that turns to dense oscillatory spikes due to the large shift parameter k = 10- Figure 4b. For large n = 5 pulse, and non-zero chirp C = 6, there is a main peak at D M = 0 with irregular and asymmetrical oscillatory envelopes, Figure 4c.
For initial number state (n o = 5), and weaker Rabi frequency (Ω o = 1) and large pulse n = 5, in the off-resonance case (∆ o = 10) the spectrum is a skewed single peak with little grooving on its envelope for D M > 0, Figure 5.
Therefore, the Morlet WT spectrum and with initial vacuum state and single pulse (n = 1) has dense oscillatory structure due to large shift parameter of the Morlet window function, Figure 4b or irregular oscillations due to large pulse number and non-zero chirped parameter.

Summary
The time-dependent wavelet spectra of the fluorescent radiation due to the coupled system of non-dissipative harmonic oscillator (HO) with a train of chirped Gaussian pulses are calculated analytically and investigated computationally using two wavelet window functions: the Haar and Morlet wavelet functions. Asymmetry, narrowing and dense oscillations in the spectra in the two cases of initial and number states of the HO are essentially associated with large shift parameter, non-zero chirp and large pulse numbers. In general, these features are absent in the Fourier transform spectrum [5], due to the averaging technique over the entire time interval of the detected signal (c.f. [9]). Study of pulsed-driven qubit system is very essential within the context of quantum information sciences and quantum computing. In particular, the spectrum function of the scattered radiation provides information about the splitting of the atomic energy levels, so the response of the atom for a certain frequency can be identified [20]. Finally, it is worth mentioning physically that the driven HO realized as radiation-radiation interaction may be also realized as radiation-matter interaction, with matter sought as N → ∞ number of atoms bathed in a thermal or squeezed vacuum radiation [21][22][23]. Therefore, our results provide deep insight into the giant quantum oscillator (formed as noted in this N → ∞ limit) and driven by the chirped Gaussian laser pulses.